大坝混凝土三维细观力学数值模型研究pdf.docVIP

第25卷第2期 计算力学学报 V01．25，No．2 2008年4月 Chinese Journal of Computational Mechanics April 2008 文章编号：1007—4708(2008)02一0241一07 大坝混凝土三维细观力学数值模型研究 马怀发“， 陈厚群1， 吴建平2， 黎保琨3 (1．中国水利水电科学研究院工程抗震中心，北京100044}2．并行与分布处理国家重点实验室，长沙410073； 3．北京工业大学，北京100044) 摘 要：在细观结构层次上将大坝混凝土作为骨料、固化水泥砂浆及其粘结界面组成的复合材料．建立了大坝混 凝土三维细观力学数值模型。该模型既能够反映混凝土及其细观各相材料在荷载作用下的损伤演化过程．又考 虑了动载作用的应变率强化效应．给出了该数值模型求解方法，并编制出能够在普通PC机上运行的串行程序。 加载过程既可按荷载控制又可按位移控制。同时，为了减少求解自由度应用了分尺度方法以使最小骨料和固化 水泥砂浆混合后其力学性能与一种复合介质等效。通过混凝土湿筛和三级配试件的静、动(冲击)弯拉数值计算 验证了本文计算方法和程序正确有效。另外．在串行程序的基础上，优化了刚度矩阵的存储方式，采用双门槛不 完全Cholesky分解(ICT)预条件的共轭梯度法(CG)，完成了能够在Sun Fire 6800服务器实现并行计算的并行程 序改造。从而大大提高了计算效率． 关键词：大坝混凝土；三维随机骨料模型，分尺度方法；预条件；并行计算  中图分类号：TV313 文献标识码：A ． j． 一 J 引 置 混凝土是由水、水泥和粗细骨料组成的复合材 料。骨料尺寸及其级配、水泥用量、水灰比、固化强 度、养护条件、环境湿度和混凝土的发热量等因素 直接影响混凝土力学特性。不均匀性是混凝土材 料的本质特点，微裂缝是决定其性能的主导因素。 为了反映混凝土细观结构不均匀性特征，对其损伤 断裂过程进行细观力学数值模拟，人们提出了许多 细观力学模型，如格构模型，随机力学特性模型和 随机骨料模型等。但到目前为止，对混凝土细观力 学的数值模拟研究工作主要限于平面问题[1-]。尽 管有三维骨料模型的研究[5]，但基本上停留在对 “凸型”多面体骨料随机投放算法的研究，并且骨料 含量低，没有考虑实际骨料级配。文献[6]提出了 按照混凝土实际级配的骨料含量和比例生成球形 骨料模型的随机投放方法，并在球形骨料模型的基 础上生成了接近于实际骨料形态的多面体骨料模 型。本文将在其基础上，针对混凝土三维细观力学 收稿日期：Z006—03—15f修改稿收到日期：2006—10—16， 基金项目：国家自然科学基金重点资助项目． 作者简介：马怀发’(1962·)，男，高级工程师 (E-mail；mahf(园iwhr．tom)l 陈厚群(1932一)。男．中国工程院院士．  数值分析方法及其编程技巧进行分析研究。 2基本方程 2．1动力学方程 混凝土及其细观各相材料的抗拉强度的强化 系数H。和弹性模量强化系数H￡与应变率的关系 表示为 H，=exp{[A，(19㈦+B，)]c，) 魄一exp{[A￡(19l三I+B￡]％)) (1) 式中A，，且，Cf为强度强化参数；AE，BE，CE为弹 性模量强化参数。动拉强度、动弹性模量可表示为 五(三)=H，Z，，E(三)=H￡E， (2) 式中正，为混凝土材料的静态抗拉强度，E为静态 弹性模量。^(；)及EG)表示混凝土动抗拉强度和 动弹性模量，是拉(或压)应变率；的函数，并认为 泊松比不受率效应的影响。 动力学增量方程‘73为 [明(△机(￡)}+[c](△抚(￡)}十 [KJ(￡)](AUa(￡))= (△Pd(￡)}+{△P：(D))+{△P，(D)} (3) 式中[明为质量阵，[C]为阻尼阵，[圮(￡)]表示  万方数据 计 算 力 学 学 报 第25卷  动刚度阵；{△吼(￡))，(△讥(￡))，{AUa(t)}分别为 结点加速度增量、速度增量和动位移增量， (△R(￡)}为结点动荷载列阵增量。{APd(D)}为 在动位移上因损伤产生的应力转移而产生的附加 荷载，{△P，(D))为在静位移上因损伤产生的应力 转移而产生的附加荷载。 采用文献[3]给出的双折线弹性损伤演化模 型，残余强度系数A=^，／^，^为混凝土及其细观 各相材料的抗拉强度，厶破坏单元的抗拉残余强 度；与残余强度对应的应变系数叩一Cr／e。，e。为单