101基于MSCMARC的表面裂纹扩展特性的仿真_2010MSC年会论文002-刘青峰

基于MSCMARC的表面裂纹扩展特性的仿真刘青峰1，谢基龙2（1.海军后勤技术装备研究所，北京100072；2.北京交通大学机械与电子控制工程学院，北京100044）摘要：基于MSCMarc软件，本文通过编程实现了参数化、模块化建模与自动化分析；采用模块化有限元建模技术和接触技术仿真半椭圆表面裂纹(前缘)，采用在裂纹前缘形成辐射状奇异单元网格的建模方法和接触技术来施加边界条件的方法，实现裂纹前缘奇异应力场和裂纹模型远场应力的模拟。在此基础上，运用裂纹虚拟扩展技术，计算裂纹前缘的J积分,通过线弹性有限元理论，换算出裂纹前缘的K因子分布，对裂纹前缘离散节点实行正交扩展，通过拟合裂纹前缘，实现了裂纹连续扩展仿真。仿真结果与试验结果基本吻合，说明本文所采用的方法具有很好的实用性。关键词：接触技术；奇异单元；裂纹前缘；半椭圆表面裂纹；MSCMarcSimulationResearchofContinualExtensionRuleof3DSurfaceCrackUsingMSCMARCSoftwareLIUQingfeng1，XIEJilong2（1．NavyLogisticTech．&EquipmentInstituteofPLA，Beijing100072，China；2．SchoolofMechanicalandElectronicControlEngineering，BeijingJiaotongUniversity，Beijing100044，China)Abstract：TheautomatizationofparameterizationandmodularizedmodelingandanalysisarerealizedthroughprogramusingMSCMarcsoftware.Themodularizedfiniteelementmodelingandcontacttechniqueswereadoptedtosimulatethesemi-ellipsesurfacecrack.Radialsingularunitmeshesinthecrackfrontwereusetosimulatethesingularstressfieldofthecrackfront.Atthesametime,theboundaryconditionsofcrackmodelbyfiniteelementcontacttechniquewereappliedtosimulatethefarfieldstressofthecrackmodel.Inaddition,thecrackvirtualextensiontechniquewasusedtocalculatetheJintegraloncrackfront.Thenthestressintensityfactordistributioncanbetransformedbasedonthelinearelasticfiniteelementtheories.Ontheotherside,thecrackfrontwasdispersedtorealizeorthogonalextensionandthenewcrackfrontisfittedtogetcontinualextensionsimulation.theextensionofsemi-ellipsesurfacecrackwassimulated.Theresultwasinaccordancewiththefracturedataofthesamplewithbeachmarkings,themodelofsemi-ellipsesurfacecrackextensionsimulationbuiltinthepaperwasverifiedtobereasonableandfeasible.Keywords：contacttechnique；singularunit；crackfront；semi-ellipsesurfacecrack；MSCMarc国内外对于裂纹的研究，多数集中在理论推导和试验上[1]，只能通过理论或经验公式指导工程结构的可靠性设计。许多研究者[2-6]介绍了表面裂纹的应力强度因子的各种理论计算方法和分布特征来表征裂纹的扩展。随着计算机技术和有限元技术的快速发展，使裂纹及其扩展仿真成为可能，为含裂纹结构可靠性研究创造了条件[7-11]。裂纹虚拟扩展技术和J积分能量定义是有限元法计算裂纹前缘应力强度因子分布的理论基础；而有限元法计算裂纹前缘K因子分布规律是实现裂纹连续扩展的前提条件，从计算裂纹前缘的J积分（能量释放率）分布着手，在线弹性情况下可以换算裂纹前缘的应力强度因子分布。本文通过软件实现裂纹前缘的不同步扩展。事先不给出裂纹前缘曲线的形状假设，而将裂纹前缘离散成系列点，计算裂纹前缘应力强度因子的分布，即局部（各个离散点）应力强度因子，并对裂纹前缘的各个离散点单独实行正交扩展，拟合得到的代表新裂纹前缘的系列离散点，来获得裂纹前缘曲线的形状表达，实现裂纹前缘的不同步扩展(裂纹前缘各部分扩展速度不同)。这样可获得更加精确的典型裂纹的扩展特性和寿命。1裂纹连续扩展模型的建立为了真实地仿真裂纹尖端应力场，本文采用MSCMARC仿真裂纹、施加载荷与边界条件，采用参数化、模块化建模技术建立裂纹连续扩展模型。图1是含有半椭圆表面裂纹体。图1含有半椭圆表面裂纹块体1.1裂纹的参数化、模块化建模在MSCMARC软件中，含有半椭圆表面裂纹的四分之一模型被划分为7个模块，各模块分别表示为：裂纹载体tets，裂纹体包括三个子模块：crack1模块、crack2模块、fill模块，z方向对称面cracksym，x方向对称面symmx，载荷施加面moving。它们之间的位置关系如图2所示。把裂纹体镶嵌入裂纹载体中，从而可单独对裂纹体划分特定的映射网格，裂纹载体划分为四面体网格。图2各个模块之间的位置关系1.2裂纹前缘及其节点的定义通过接触技术定义，z方向对称面cracksym只对裂纹韧带部分而不对裂纹自由表面运用对称边界条件，即cracksym仅与裂纹的tets和crack1模块粘连在一起，而与crack2和fill模块没有连接关系。在模块crack1和crack2交界的曲线上，两模块的节点是重合的，但是crack1的节点属于裂纹韧带，与对称面cracksym粘连（Glue）在一起；crack2的节点属于裂纹自由表面与cracksym没有连接关系；这样在模块crack1和模块crack2交界曲线处形成了三维的裂纹前缘。相应的裂纹前缘曲线的离散节点用来求解局部J积分。2半椭圆裂纹模型的前缘离散正交扩展与拟合图3-1、图3-2、图3-3、图3-4给出了半椭圆裂纹前缘扩展拟合过程的图形表示。图3-1表示按弧长均匀离散裂纹前缘，使用离散节点代替裂纹前缘；图3-2表示前缘离散节点正交扩展得到新前缘的离散点，新获得的点不能直接用来表示裂纹前缘，要获得裂纹前缘的数学描述，需要进行数据拟合；图3-3表示椭圆拟合新的离散节点，获得椭圆前缘参数；图3-4表示按弧长均匀离散新拟合的椭圆裂纹前缘；重复图3－2、图3－3、图3－4所示的操作就实现了裂纹的连续扩展。图3-1按弧长均匀离散裂纹前缘图3－2前缘节点正交扩展图3-3椭圆拟合新的裂纹前缘图3-4按弧长离散新裂纹前缘2.1正交扩展计算裂纹前缘节点坐标裂纹前缘的离散节点数在20～30个节点之间。为了清楚说明问题，图形中裂纹前缘的节点取为5个。设椭圆裂纹长半轴为A，与x轴平行，短半轴为B，与y轴平行；坐标原点、椭圆的中心为O，如图4所示。图4裂纹前缘节点正交扩展椭圆方程为：12222ByAx（1）裂纹前缘系列离散点的坐标序列记为nniiyxyxyx,,,,,,,11。在点iiyx,作椭圆的切线iL1，则点iiyx,的导数就是切线iL1的斜率，记为ik1，对椭圆方程12222ByAx两边求x的导数，求得：ikyAxBdxdy122（2）过iiyx,点与切线iL1垂直的直线iL2的斜率记为iikk121；在点iiyx,，裂纹的扩展量为ia，需要求解在直线iL2上，距点iiyx,长ia的点的坐标，方法如下：①直线iL2的方程为：iiiixxxByAxxkyy222（3）②过iiyx,点半径长为ia的圆的方程：222iiiayyxx（4）将直线iL2方程平方代入圆方程中，设iikk231，解出：2221iiikaxx（5）2321iiikayy（6）在程序中有效点的判断方法为：将11,iiyx，22,iiyx代入椭圆方程，如果大于1就是点在椭圆外，即点11,iiyx；如果小于1就是点在椭圆内，即点22,iiyx。裂纹是向外扩展，点11,iiyx为有效点。2.2拟合裂纹前缘及其步长控制对裂纹前缘系列离散点进行椭圆拟合，变量代换可将椭圆拟合转化为线性拟合：设2xX、2yY，椭圆函数模型就可转化为：222BXABY（7）设22ABb，2Ba，直线方程变为：abXY，由最小二乘法直线拟合可以算出a和b，这样可以解出椭圆的长短轴A和B。图5椭圆拟合裂纹前缘的误差分析如图5所示，坐标系中心为O，椭圆中心也为O0,0，裂纹前缘正交扩展后节点的坐标序列记为nniiyxyxyx,,,,,,,11，过裂纹正交扩展获得新节点与椭圆原点作直线，新节点到椭圆原点的距离为22iiyxOQ，拟合椭圆的交点到椭圆原点的距离OP。P点的坐标是直线与椭圆的交点。定义裂纹前缘节点的相对偏差：相对残差=OQOQOP（8）相对残差平方和均值开方=niiOQOQOPn121（9）通过控制（8）式表示的误差不超过1％来确定裂纹扩展的步长，同时保证了拟合裂纹前缘的精度。3半椭圆表面裂纹连续扩展仿真试验3.1仿真试验参数仿真试验采用ZG25材料作裂纹块体，其扩展速率为：da/dN=1.44×1010（△K）2.79；门槛值△Kth=6.3246Mpa·m1/2。取含有半椭圆表面裂纹块体的四分之一为有限元模型，厚度t为10mm，宽度B为15mm，长度L为15mm，参见图1；a为裂纹深度（短半轴），c为裂纹长度（长半轴），a0和c0为初始值。裂纹前缘划分为22份，即离散为23个节点；每个节点设置4个J积分半径，对应可以计算4个J积分值。按照脉动载荷计算K因子范围。仿真试验的载荷参数如表1所示。载荷条件为拉伸、弯曲、拉弯三种，确定载荷的大小是使裂纹前缘有4个离散点局部K因子大于门槛值。采用相同初始半椭圆表面裂纹模型：02mma，00/0.8ac；进行等步长裂纹扩展仿真。表1仿真试验载荷参数载荷大小0a00/ac步长扩展步数前缘个数拉伸16000N20.8不变2930弯矩60000Nm20.8不变2829拉弯8000N&15000Nm20.8不变30313.2仿真试验结果仿真试验结果如图6所示。由图6-1、图6-2、图6-3中裂纹前缘变化规律可见，在裂纹扩展相同长度（长轴c方向）时，拉伸载荷作用下，裂纹沿深度方向（短轴a方向）扩展较快，而弯曲载荷作用下则较慢，拉弯载荷作用下介于两者之间；当裂纹深度超过一定尺寸后，三种载荷条件下裂纹沿深度方向的扩展均逐渐放慢，而沿长度方向逐渐加快。图6-4的a/t－a/c关系也反映出在拉伸载荷作用下，裂纹沿深度方向扩展较快，在弯曲载荷作用下裂纹沿长度方向扩展较快，拉弯载荷作用下介于两者之间。图6-1拉伸载荷下裂纹前缘变化规律图6-2弯曲载荷下裂纹前缘变化规律图6-3拉弯载荷下裂纹前缘变化规律图6-4三种载荷下半椭圆表面裂纹的a/t－a/c曲线3.3仿真结果与试样实验结果比较本文进行了拉伸载荷下的裂纹扩展实验，将仿真结果与实验结果进行比较，得出图7所示结果。文献[12]是对含