基于积分点的材料分区方法研究.docVIP

基于积分点的材料分区方法研究 邱晓峰 河海大学结构工程，南京 （210098） E-mail： uiuxf0203@ 163.com 摘 要：在有限元计算屮，通常将儿何模型按照材料区域的不同来分块，以此来实现模型屮 不同材料的分区。而基于高斯积分点的材料分区可以将材料与有限元网格分离，避免父杂分 区带来的模型网格划分上的困难。本文运用有限元软件ABAQUS对相同的模型分别使用两 种材料分区方法进行计算。分析结果发现基于高斯积分点的材料分区方法计算精度合理，nJ 以用于具有复杂材料分区模型的有限元计算屮。关键词：材料分区；高斯积分点；自定义 材料 1引言 我国水利水电发展迅速，在西南地区水资源十分的丰富，建设了大量的拱坝、重力坝。 但西南地区地属地憲带，地质条件非常的复杂，岩体常常带有节理、断裂层，断层破碎带、 大裂隙、软弱夹层等。由于基岩的成因和构造复杂，岩性多样，往往存在着多条的软弱结构 血和多纟H.的节理裂隙，其物理力学参数具有明显的分散性和不确定性，特别是其上作用荷载 时其结构面的形态和性质对坝体的稳定性起着重要的作用⑴。若岩体被多组断层、裂隙切割, 失稳的岩体大小、形状、滑动形式以及结构面的物理特征等具有复杂性，要确定单元的类型, 准确的模拟,基岩的受力及分析其位移和稳定性，在材料的分区建模上就存在相当大的困难。 基于高斯点的材料分区是根据积分点的位置的判断赋予材料屈性，省去了在建模时的按照不 同的材料分区，从而大大降低了建模和划分网格的困难。但是此方法是需要细化单元的网格 划分，来模拟真实的应力应变、位移，要求比较多的计算时间。 2计算方法 2.1平面应力问题 在平面应力状态下，运用ABAQUS分析模型的应力和位移，进行分区与不分区的对比。 2.1.1采用对不同材料进行分区域建立平面模型，材料一：E, = 2e9 , =0.3；材料二： E2 = k8 , v2 =0.3o荷载在顶部施加均布压应力，其值为10000c在底部施加X、Y方向约 朿，其余备边白由。 计算mises应力如图,位移如图?2所示: 3, Mxses(Ave? Crxt?:75%)2211111111876236e*04 102e+04 967e+04 832e+04 697e+04 562e+04 427e+04 293e+04 158e+04 023e+04 881e+03 533e+03 184e*03ODB: pl-1?odbStep: 3tep-l Increment Primary Var:Deformed Var:ABAQUS/STANDARD Version 6?5-1 Mon Mar 31 21:30:10 GMT+08:00 2008 3, Mxses (Ave? Crxt?:75%) 2211111111876 236e*04 102e+04 967e+04 832e+04 697e+04 562e+04 427e+04 293e+04 158e+04 023e+04 881e+03 533e+03 184e*03 ODB: pl-1?odb Step: 3tep-l Increment Primary Var:Deformed Var: ABAQUS/STANDARD Version 6?5-1 Mon Mar 31 21:30:10 GMT+08:00 2008 1: Step Time = 丄?000 3, MisesU Def orma-tion Scale Factor: +3 ? 775e+03 图对材料进行分区的平Ifu* miscs应力 ABAQUS/STANDARD Version 6?5-1-+1130e03-+1036色03414e04-+8 ABAQUS/STANDARD Version 6?5-1 -+1 130e 03 -+1 036色 03 414e 04 -+8 473e 04 ? 47 531e 04 ■ +6 590e 04 -+5 64 9e 04 707e 04 +3 766e 04 +2 824e 04 883e 04 *9 414o 05 *0 000e+00 Ux Magnitude LODB: pl-1?odb kep: Step-l Mon Mar 31 21:30:10 GMT+08:00 2008 Increment 1: Step Time =Primary Var: U, Magnitude 1.000 Deformed Var: U Deformation Scale Factor:? 775e-*-03 图?2对材料进行分区下的位移分布 2?1?2采用对材料的高斯点进行材料的划分，建立整体模型。在ABAQUS屮采用H定义材 MUser Material:E, =2e