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