大坝混凝土动态性能细观力学及并行计算研究ppt课件

大坝混凝土动态性能细观力学及并行计算研究,马怀发2006年6月24日,一、研究的意义和试验成果分析二、混凝土细观力学研究方法三、混凝土试件数值模拟的基本方程四、混凝土试件的二维数值模拟五、混凝土三维细观有限元并行计算研究六、结束语,汇报大纲,1、研究的目的和意义,1）重大工程抗震安全评价的主要环节地震作用工程结构地震响应材料动态抗力大坝混凝土动态抗力是大坝抗震安全评价的薄弱环节：动态强度静态强度提高30%，依据国外早期少量湿筛试件的试验结果:大坝混凝土的湿筛试件和全级配试件动态性能存在差异。2）拱坝受力特点：(1)高拱坝是一种空间结构，处于偏心受压状态。在强震时，坝体薄弱部位主要因坝体混凝土的动弯拉应力引起开裂，因此，动弯拉强度是高拱坝抗震设计中的控制性指标之一;(2)在地震作用时，坝体各部位都已经存在不同程度的静态应力。因此,有必要研究:在不同预静载作用下，混凝土动态性能的变化规律。,1）小湾四级配混凝土和湿筛混凝土静动态弯拉试验成果：,2、大坝混凝土的动态试验成果,2)小湾拱坝三级配混凝土，以冲击加载和三角波加载两种方式，在不同预静载水平下得到的动弯拉强度：,3）试验结果发现：a)全级配大坝混凝土的弯拉强度低于湿筛混凝土的弯拉强度,纯动弯拉强度增强系数都低于30%；b)全级配混凝土的动弯拉强度增强系数高于湿筛混凝土的动弯拉强度增强系数;c)一定预静载水平对动弯拉强度产生强化作用。4）数值摸拟的出发点：a)反映混凝土材料的内部结构特征，研究全级配混凝土与湿筛混凝土的细微结构差异；b)考虑混凝土材料的动载强化作用；c)反映预静载对混凝土试件动损伤破坏过程的影响。,一、研究的意义和试验成果分析二、混凝土细观力学研究方法三、混凝土试件数值模拟的基本方程四、混凝土试件的二维数值模拟五、混凝土三维细观有限元并行计算研究六、结束语,汇报大纲,(1)微观层次，材料的结构单元尺度在原子、分子量级，即从小于10-7厘米至10-4厘米。(2)细观层次，尺度范围在10-4厘米至几厘米,或更大一些。在这个层次上，混凝土被认为是一种由粗骨料、硬化水泥砂浆和它们之间的过渡区（粘结带/界面）组成的三相材料。(3)宏观层次，特征尺寸大于几厘米。材料假定为均质。要研究混凝土级配对混凝土宏观力学特性的影响，必须考虑混凝土细观结构的差异。,1、混凝土材料的细观结构层次,2、细观数值模型,为了建立混凝土细微观结构特性与其在宏观力学特性的关系，自上个世纪70年代末，人们发展了混凝土细观力学研究方法。细观力学将混凝土看作由粗骨料、硬化水泥胶体以及两者之间的界面粘结带组成的三相非均质复合材料，在细观层次上划分单元。考虑骨料单元、固化水泥砂浆单元及界面单元材料力学特性的不同，以其细观各相材料的破坏准则或损伤模型反映单元刚度的退化，利用数值方法模拟混凝土试件的裂缝扩展过程及破坏形态。,(1)二维随机骨料模型瓦拉文用三维富勒级配曲线推导出了试件内截面任一点具有某种骨料直径的内截圆出现的概率，即瓦拉文公式。对于一个特定的混凝土拌和物，使用这个公式即可产生试件截面上骨料的颗粒数。利用瓦拉文将三维问题简化为二维问题,1)二维随机骨料模型,2）二维随机骨料模型细观有限元网格,(2)三维随机骨料模型1)随机球骨料模型分布将骨料简化为球体，按照实际级配的骨料粒径和各粒组骨料的比例，确定各种骨料的颗粒数。,2)随机凸多面体骨料模型分布在球形骨料模型的基础上形成多面体骨料基，以球骨料含量为控制参数，进行随机延凸，得到随机凸体骨料模型。,3)随机骨料细观有限元剖分,一、研究的意义和试验成果分析二、混凝土细观力学研究方法三、混凝土试件数值模拟的基本方程四、混凝土试件的二维数值模拟五、混凝土三维细观有限元并行计算研究六、结束语,汇报大纲,1、混凝土材料的应变率强化现象（试验成果）,应变率效应是固体材料的基本特性。非均质材料较均质材料受应变率的影响更为显著。混凝土在变形可逆的弹性阶段，几乎无应变率敏感现象。仅当加载至内部出现损伤时才呈现加载速率的影响。应变率对混凝土动弹模的影响有类似的强化规律，但其影响弱于对动强度的影响。混凝土动拉压强度随应变速率增长，二者规律相似。,混凝土抗拉强度和弹性模量与拉（或压）应变率的关系：(1)式中fts为混凝土材料的静拉强度，Es为静弹性模量。及表示混凝土动抗拉强度和动弹性模量是拉（或压）应变率的函数。Ht为强度强化系数；HE为弹性模量强化系数，二者与应变率关系为：(2)式中：At、Bt、Ct为强度强化参数，AE、BE、CE为弹性模量强化参数。,2、混凝土的应变率强化定律,3、混凝土损伤演化定律,混凝土拉伸变形起控制作用,采用了各向同性弹性损伤力学的本构关系来描述混凝土材料的力学性质。按照勒梅特应变等价原理，受损材料的本构关系：式中E0为初始弹性模量，为损伤后的弹性模量,D为损伤变量。,（0D1）,4、混凝土试件数值模拟的动力学方程,在动力学方程中，考虑到损伤弱化、应变率强化以及初始预静载对动态性能的影响。1）动力学微分形式：其中,动载增量,在动位移上因损伤产生的应力转移而产生的附加荷载,在静位移上因损伤产生的应力转移而产生的附加荷载,2）动力学方程的积分形式：,一、研究的意义和试验成果分析二、混凝土细观力学研究方法三、混凝土试件数值模拟的基本方程四、混凝土试件的二维数值模拟五、混凝土三维细观有限元并行计算研究六、结束语,汇报大纲,计算模型（湿筛试件）,1、预静载及加载速率对混凝土动弯强度的影响,1）加载速率对动弯强度的影响,动弯拉强度不但与加载速率有关，而且还与初始静荷载有关。,2）初始预静载对动弯强度的影响,计算表明：混凝土梁的动弯拉强度不仅与加载速率有关，而且与应变率、损伤积累密切相关。,初始静预载对混凝土静动综合弯拉强度的影响,1）二维随机骨料随机参数模型(RAPM)(1)骨料随机分布（同随机骨料模型。按骨料级配及含量）(2)混凝土及其细观各相单元的抗拉强度和弹性模量随机分布这里所说的不均匀性包括混凝土骨料颗粒级配的差异和材料参数分布的统计随机性；随机骨料随机参数模型是在随机骨料模型的基础上，在细观有限元剖分后，再生成各类单元弹模和拉强度参数分布。,2、混凝土细观结构差异对其动弯强度的影响,2)细观结构对弯拉强度的影响,混凝土细观各相材料力学特性的不均匀性，使其动弯拉强度增强系数对初始静预载的敏感性增强。,湿筛混凝土细观有限单元模型（单位:mm）,细观不均匀性对静动综合弯拉强度的影响,3）混凝土级配对动弯拉强度的影响,全级配混凝土的动弯拉强度增强系数高于湿筛混凝土的动弯拉强度增强系数。,全级配混凝土骨料颗粒分布及单元剖分,混凝土级配对静动综合弯拉强度的影响,4)骨料形态对混凝土对弯拉强度的影响,1)湿筛混凝土试件,湿筛混凝土两种细观骨料模型结算结果比较,2)全级配混凝土试件,3、梁裂缝扩展分析1)静载作用下弯拉失稳过程,2)静动弯拉失稳破坏机理对比分析(a)静力:作用时间相对较长，材料内部的损伤主要表现为微裂缝的扩展、联结和贯穿然后导致破坏；(b)动力:因此材料内部的损伤主要表现为大量的微空洞和微裂纹生成，但来不及扩展。,一、研究的意义和试验成果分析二、混凝土细观力学研究方法三、混凝土试件数值模拟的基本方程四、混凝土试件的二维数值模拟五、混凝土三维细观有限元并行计算研究六、结束语,汇报大纲,1、并行计算研究的意义1）混凝土黏结带区域薄，一般在40-50微米；2）全级配混凝土试件尺寸大；全级配试件三维数值模型的节点自由度达到几十万甚至上百万，要进行大尺度全级配试件三维数值模拟，必须想办法解决计算速度这个瓶颈问题，必须进行并行计算研究。,1、并行计算的两种方案1）对现有混凝土细观力学分析串行程序进行了并行算法优化设计（国防科大合作）：首先，从整体上提出了一个将有限单元分布与未知量分布有机结合的整体并行算法设计方案。在进行并行算法设计时，整体上按有限元个数平均分配来进行任务划分，假设一共有p个处理器，将所有有限元分成p部分，每一部分中的有限元分配到一个处理器上。对一个给定的有限单元，如果分配到第k个处理器，则与之有关的计算在第k个处理器上进行。在对各个单元进行计算时，各个单元的计算之间是互不相关的，可以完全并行。然后分别针对刚度矩阵装配、双门槛不完全Cholesky分解预条件、稀疏矩阵与向量相乘、稀疏向量相加等核心算法，提出了相应的高效并行算法。模拟标准湿筛试件静加载损伤破坏过程，其三维模型共有71013节点，78800单元，求解自由度约为21万。计算约需要8个小时。计算平台为水科究院网络中心Sun6800（4个主频1GCPU、8G内存、540G的磁盘整列存储空间）,2）基于并行有限元程序自动生成系统PFEPG（1）并行有限元程序自动生成系统PFEPG的基本思想：元件化程序设计程序自动生成技术区域分解法将有限元程序分成：前处理分区、start、bft、solv、E和U等六个模块,其中E和U元件程序根据用户给出的偏微分方程描述文件（vde或pde），算法文件（nfe）由系统自动生成，其余几个元件程序由系统给出，不随表达式的变动而变动。,（2）并行计算的程序结构,（3）混凝土细观数值模拟的FEPG文件gio命令流文件stdya位移场对应单元计算文件为stdya.vdestdyb应力场对应单元计算文件为stdyb.vde空一行#elemtypec8g2单元类型，8节点六面体2阶高斯积分3dxyz三维直角坐标系,gcn命令流文件defi关键字astdya&a场求解位移，对应a场算法文件为stdya.nfe&表示生成的文件不重新命名bstdybb场求解应力，对应算b场法文件为stdyb.nfe空一行，表示定一段结束，以下通过命令流组织计算流程startilua对应预处理共轭梯度法不存单刚方式位移场初始化ifexiststopdelstop:1bfta边值计算ifexistenddelend:2solvilua预处理共轭梯度法求解位移ifnotexistendgoto2控制非线性迭代stressb最小二乘法求解应力场posta按时间步输出结果ifnotexiststopgoto1控制时间循环,nfe算法文件（1）位移算法stdya.nfe文件按照动静力学方程求解步骤，用FEPG有限元语言编写的求解位移的nfe文件。（2）应力算法stdyb.nfe文件在求得位移后，用最小二乘法求应力。表示由已知变形算得的应力，为待求的光滑后的应力。,vde文件a)单元位移计算的stdya.vde文件b)单元应力计算的stdyb.vde文件,5、PFEPG程序的生成和运行1)填写Partition文件2)运行mpi命令：mpigcn文件的前缀名3)把全部并行fortran源程序（*.f文件），fegen.b、runmpi、time0、nbefile、partition.dat、*.io传送到并行机4)在并行机上运行命令:cshxfegen.b5)运行前处理6)修改hosts文件和runmpi文件7)在并行机上运行如下命令:cshxrunmpi8)运行后处理,(1)并行计算环境1）硬件环境：4台以主频3G、内存1G的PC机以千兆以太网连接的工作站。2）软件环境：Linux操作系统,PFEPG并行有限元程序自动生成系统平台采用MPICH通信。,6、并行计算分析,(2)细观模型选取及并行程序生成1）计算模型：整个剖分区域共有111988个单元，76810个节点。,图55三维随机骨料分布,a）一种单元类型对应的一个vde文件；b）跨中细观剖分区域与其它区域独立分开；c）执行gio命令流文件后修改pre文件和对应各物理场以io为扩展名的文件；d）再执行pre命令和其它vde文件（除gio文件中vde文件）e）生成单机串行程序，再按PFEPG程序步骤生成该问题的并行程序。f）修改材料参数（修改elem0和elemb0中的材料参数），进行并行计算。,模型并行计算分区图,partition.dat文件为424个分区，2表示按坐标分区221x方向剖分数为2，y方向剖分数为2，z方向剖分数为1,(3)预静载作用下三维细观模型并行计算分析,图中2dp1为第三章二维串行程序的,三维模型并行计算是在主频为3G、内存为1G的4个PC节点的并行集群上进行的。最大控制位移为5.010-4，非线性迭代位移误差范数控制为10-15，最大迭代次数为15次。静力计算耗时约15.6小时；动力计算耗时约21小时；其它有预静载的情况介于两者之间。,图59预静载对静动综合强度的影响,一、研究的意义和试验成果分析二、混凝土细观力学研究方法三、混凝土试件数值模拟的基本方程四、混凝土试件的二维数值模拟五、混凝土三维细观有限元并行计算研究六、结束语,汇报大纲,全级配混凝土的动态性能研究是当前我国高拱坝抗震分析中需要深入研究的高难度的前沿科研课题。已取得的主要研究成果有以下几个方面：1）目前我们已经初步建立了混凝土细观静、动力学分析体系的框架。这个体系由以下几部分组成：(1)混凝土材料力学特性的应变率强化定律；(2)混凝土损伤演化关系；(3)既考虑细观损伤弱化又计及应变率强化效应的混凝土非线性