4维非均匀脆性材料破坏过程数值模拟.pdfVIP

第34卷第3期 力 学 学 报 Vol34．Nn：j 2002年5i ACTAMECHANICASIN／CA May，2002 三维非均匀脆性材料破坏过程的数值模拟 陈永强郑小平姚振汉 ‘清华大学工程力学最北京100084’ T酤 日 摘要采用有限元方j去模拟了三维非均匀固体材料的宏观等效力学性质和破坏过程．首先采用 格形(1attice)方法把试件离散成三维均匀网格，在每个单元格中将材料按照均匀处理，根据给 定的统计规则来确定不同单元格中的材料常数以反映材料的非均匀性．然后对非均匀脆性材料 选用简单的本构关系与断裂准则，采用自适应选取载荷步长对试件进行加载，通过非平衡迭代 技术对刚度矩阵进行不断修正，实现了非均匀脆性材料的弹性行为及破坏过程的数值模拟在 此基础上，通过数值计算研究了材料的非均匀分布对宏观等效力学性质和破坏过程的影响，给 出了破坏全过程的非线性载茼一位移曲线以及不同载荷阶段的三维损伤破坏的演化图． 关键词 三维非均匀，破坏过程，格形模型，统计方法，数值模拟 严格来讲自然界的一切材料都是非均匀的，只不过彼此非均匀程度不同而已．显然采用传 统的均匀模型对宏观和细观非均匀材料进行分析计算不再合适，而对于微观非均匀材料在破坏 过程中表现出的一些特殊力学行为和现象，均匀模型也不能给出满意的解释．所以建立非均匀 材料的力学模型和相应的分析手段是当前计算力学研究的重要课题之一． 对于非均匀材料问题，国内外学术界进行了长期探索，特别是近■十年来取得了很大的进 展，并逐渐形成了一个新型学科分支．从研究方法来看，比较流行的主要有两种类型：一种是 建立在各种各样自洽模型基础上的有限元与边界元方法，另一种是建立在统计技术上的格形化 有限元方法．这两种方法已被成功地应用于各种非均匀材料问题，特别在复合材料力学性重研 究和脆性材料破坏过程的演化方面取得了令人瞩目的成果． 采用简单格形模型模拟了混凝土等脆性材料的典型破坏机理； 1997年G．Frantziskonis等口’3J 假定材料的弹性性质为一个平稳随机场，分析了材料非均匀性对结构边界附近位移场的影响， 并且得到了一维情况下的解析解；1998年s．C．Blair和N．Cook【4j』采用统计方法研究了岩石 材料的非线性行为及由于颗粒形状和尺寸的变化而引起的应力场微观非均匀件；1998年P．K． Kaiser和C．A．Tang[6,7J用非均匀有限元模型模拟了岩石不稳定破坏现象． 但是以上所有工作均是按照二维情况进行模拟的，据作者所知，到目前为止关于三维毒均 匀材料的数值模拟很少有文献报道．最近，本文作者【8’q在这一领域进行了探索研究．在此基 础上，本文采用格形化方法(1attice)和统计技术建立了三维非均匀固体材料的有限元模型，然 后对非均匀脆性材料选用简单的本构关系与断裂准则，通过非平衡迭代技术对刚度矩阵进行不 断修正，实现了非均匀脆性材料的弹性行为及破坏过程的数值模拟．通过数值计算研究r材料 2000-02—26收到第一稿，200106_18收到修改稿 万方数据 352 力 学 学 撤 2002年第34卷 的非均匀分布对材料宏观等效弹性’陀质、破坏过程的影响，并且给出了材料破坏全过程的载茼 一位移曲线、以及在不同载荷阶段的二维损伤破坏的演化图． 1格形有限元模型 下边以三维空间简单试件受轴向载荷作用为例来建立问题的格形有限元模型．如图l所 示，试件侧表面自由，底部沿。轴方向固定，顶端施加垂直位移载荷，即 (o)=0，Ⅲ(f)=面 (1) J L面 采用格形模型对