补充材料一-断裂力学ansys方法.ppt

Click to edit Master title styleClick to edit Master subtitle styleSeptember 30, 19981*断裂力学的ANSYS求解技术断裂参数1、应力强度因子K2、J积分3、能量释放率G September 30, 2003Structural Nonlinearities 断裂力学参数的求解 先进行裂纹模型的线弹性或弹塑性静力分析 然后用特殊的后处理命令或宏命令计算断裂参数September 30, 2003Structural Nonlinearities 裂纹区域的模拟 断裂模型中最重要的区域是围绕裂纹边缘的部位裂纹的边缘，在二维模型中为裂纹顶端，在三维模型中为裂纹的前缘 September 30, 2003Structural Nonlinearities 弹性问题中，在裂纹顶端附近（或裂纹前缘）某点的位移随 而变化，r是裂纹尖端到该点的距离，裂纹顶端处的应力和应变是奇异的，随 变化为选取应变的奇异点，相应的裂纹面须与它一致围绕裂纹顶点的有限元单元是二次奇异单元，它是把单元边上的中节点放到1/4处 September 30, 2003Structural Nonlinearities 2D-断裂模型 适用于二维断裂模型的单元是PLANE2，为六节点三角形单元，围绕裂纹尖端的第一行单元，必须具有奇异性，PREP7中的KSCON命令（Main MenuPreprocessor-Meshing-Shape&Size -Concentrat KPs-Create)是指定围绕关键点的单元大小和分割排列。

本命令自动围绕指定的关键点产生奇异单元命令的其他选项是控制第一行单元的半径，以及控制周围单元的数目等为获得理想的计算结果，围绕裂纹顶端的第一行单元，其半径应该是八分之一裂纹长或更小在裂纹四周每一单元最好有3040度角度，最好是等腰三角形September 30, 2003Structural Nonlinearities 三维模型推荐使用的单元类型为SOLID95，20节点块单元围绕裂纹前缘的第一行单元应该是奇异单元，单元的KLPO面退化成KO线产生三维断裂模型要比二维模型复杂，命令KSCON不能用于三维模型，必须要保证裂纹前缘沿着单元的KO边 推荐的单元尺寸与二维模型一样，此外在所有的方向上，单元相临边之比不超过1：4 对曲线的裂纹前缘，单元的大小决定于局部曲率的数值，大概应该使沿曲裂纹前缘中每个单元只有1530度的角度所有单元的边（包括裂纹前缘）都应该是直线3D-断裂模型September 30, 2003Structural Nonlinearities 计算应力强度因子 静态分析完成后，用POST1中的KCALC命令计算复合型断裂中的应力强度因子 、 、 ，该命令仅适用于在裂纹区域附近具有均匀的各向同性材料的线弹性问题。

步骤如下： 1、定义裂纹尖端或裂纹前缘的局部坐标系，以X轴平行于裂纹面，Y 轴垂直于裂纹面，如下图所示使用KCALC命令时，坐标必须是激 活的局部坐标系(模型坐标系和结果坐标系) GUI: Utility MenuWorkplaneLocal Coordinate SystemsCreate Local CSAT Specified LocSeptember 30, 2003Structural Nonlinearities 2、定义沿裂纹面的路径，应以裂纹顶端作为路径的第一点，对于半个裂纹模型而言，沿裂纹面需有两个附加点，对于整个模型，则应包括两个裂纹面共需四个附加点，两个点沿一个裂纹面，其他两个点沿另一个裂纹面，下面图形给出了二维模型的情况 GUI: Main MenuGeneral PostprocPath OperationDefine PathSeptember 30, 2003Structural Nonlinearities 3、计算 、 和 ,KCALC命令中的KPLAN区域指定模型是平面应力或平面应变除了薄板的分析，在裂纹顶端附近和它的渐近位置，其应力一般是考虑为平面应变，KSCYM域用来指定半裂纹模型是否具有对称边界条件，还是反对称边界条件，或是整体裂纹模型。

命令：KCALC GUI: Main MenuGeneral PostprocNodal CalcsStress Int FactrSeptember 30, 2003Structural Nonlinearities 计算J积分 J积分的最简单形式可以定义为与路径无关的曲线积分，它能标志裂纹顶端附近的奇异应力和应变强度，下式是二维情况下的J积分表达式： 其中， 围绕裂纹顶端的积分路径 应变能密度 设X轴的拉力向量= 设Y轴的拉力向量= 应力分量 路径 的外法向线分量 位移向量 路径 的距离 September 30, 2003Structural Nonlinearities J积分的计算在ANSYS中的General Postproc中进行，主要是围绕上面的公式计算，其计算流程图如下图所示：September 30, 2003Structural Nonlinearities 二维模型计算J积分的步骤1.读入所要的结果命令SETGUI:Main MenuGeneral PostprocLast Set2.存储每个单元的应变能和体积命令：ETABLEGUI: Main MenuGeneral PostprocElement Table Define Table3.计算每个单元的应变能密度命令：SEXPGUI: Main MenuGeneral PostprocElement Table Exponentiate September 30, 2003Structural Nonlinearities 4.定义积分线路径命令：PATHGUI:Main MenuGeneral PostprocPath Operations Define Path5.将步骤3存储在单元表中的应变能密度映射到积分路径上命令：PDEFGUI:Main MenuGeneral PostprocPath Operations Map onto PathSeptember 30, 2003Structural Nonlinearities 6.对总体Y轴积分命令：PCALC GUI:Main MenuGeneral PostprocPath Operations Integrate7.将积分的最后值赋值给一个参数，它就是J积分的第一项。

命令：*GET,name,Path,LastGUI:Main MenuParametersGet Scalar DataSeptember 30, 2003Structural Nonlinearities 8.将应力分量SX,SY,SXY映射到积分路径上命令：PDEFGUI:Main MenuGeneral PostprocPath Operations Map onto Path9.定义路径法向量命令：PVECTGUI:Main MenuGeneral PostprocPath Operations Unit VectorSeptember 30, 2003Structural Nonlinearities 10.计算TX,TY命令：PCALCGUI:Main MenuGeneral PostprocPath Operations Operation11.沿X轴正方向和负方向将路径移动一小段距离，计算位移向量的导数 这涉及到以下步骤：计算路径移动的距离DX，一般情况下取路径的总长度的1%，可以通过以下命令得到路径总长度：*GET,Name,Path,Last,S沿X轴的负方向移动DX/2距离PCALC,ADD,XG,XG,-DX/2(即从原点处移动DX/2距离） ，将UX和UY映射到路径上，取名UX1,UY1。

September 30, 2003Structural Nonlinearities 沿X轴的正方向移动DX距离PCALC,ADD,XG,XG,DX(即从原点处移动DX/2距离） ，将UX和UY映射到路径上，取名UX2, UY2.把路径移回原点（距离-DX/2)，然后采用PCALC计算（UX2-UX1)/DX和(UY2-UY1)/DX, 它就分别代表12. 采用第10步和第11步计算得到的数据，计算J积分的第二项，并对路径的距离S积分，得到J积分第二项13.采用7和12步所获得的数据，求和得J积分的值September 30, 2003Structural Nonlinearities 能量释放率能量释放率用于计算裂纹张开或闭合时所作的功计算能量释放率的方法是虚拟裂纹扩展方法在虚拟裂纹扩展方法中，必须作两次分析，一次是裂纹长度为a，另一次是裂纹长度为a+a，如果两种情况下的位能被存储，则：其中B是断裂模型的厚度在第二次分析中，把裂纹扩展a，选择裂纹附近的所有节点，并在X方向给予缩放，以得到扩展了的裂纹模型Main MenuPeoprocessorModeling-OperateScale裂纹附近通常指裂纹尖端在半径a/2 内的全部节点，同样对节点的缩放因子a， 一般应取裂纹长度的0.5%2%。

September 30, 2003Structural Nonlinearities 如图所示薄平板，100mmX100mmX10mm，长条板有单边裂纹，裂纹长度为a=20 mm，泊松比为0.3，杨氏模量为200e9，在1MPa拉伸应力作用下，求平板的应力强度因子K和J积分实例 step by stepSeptember 30, 2003Structural Nonlinearities 一、建立模型1、在前处理中建立五个关键点1(0,0,0)、2(0.02,0,0)、3(0.1,0,0) 、4(0.1,0.05,0)、5(0,0.05,0) GUI:MainMenuPreprocessorModelingCreateKeypoints In Active CS2、建立直线,如图所示 GUI:Main MenuPreprocessorModelingCreateLines Linesstraight line September 30, 2003Structural Nonlinearities 3、通过线产生面 GUI: Main MenuPreprocessorModelingCreateAreas ArbitraryBy Lines二、定义单元类型和材料属性1、定义单元类型为PLANE2,并且在Options选项中定义Element behavior 为Plane strs w/thk。

GUI:Main MenuPreprocessorElement TypeAdd/Edit/Delete AddSeptember 30, 2003Structural Nonlinearities 2、定义实常数，板厚为0.01m. GUI: Main MenuPreprocessorRealConstants Add/Edit/DeleteAdd3、定义材料属性，E=200e9,NUXY=0.3 GUI: Main MenuPreprocessorMaterial Props Material ModelsSeptember 30, 2003Structural Nonlinearities 三、定义应力集中点1、定义关键点2为应力集中点，在出现的对话框中依次输入0.02/10、1、15、skewed 1/4pt GUI: Main MenuPreprocessorMeshingSize Cntrls Concentrat KPsCreateSeptember 30, 2003Structural Nonlinearities 四、划分网格1、用三角形单元划分网格，在裂纹处可以适当的。