哑节点断裂单元重点分析.pptVIP

* 哑节点断裂单元 * 下面介绍基于等效区域积分法断裂单元的使用方法。如图所示，断裂单元具有n个节点，因此需要给出这n个节点的节点号（N1，N2，N3，…Nn）。由于是二维问题，每个节点有2个自由度（1,2），且只需给出两个坐标。有一个输出变量，即J积分。值得注意的是，节点N1和N2 可以具有相同的初始坐标，即它们分别是重合的。需要输入的参数有（n+2）个，其中包括两个方向上的弹簧刚度系数（Kx，Ky）以及材料的弹性模量和泊松比共4个，以及定义q函数所需的（n-2）个数据。ABAQUS的断裂单元调用格式如下： q=0 q=1 哑节点断裂单元 * q=0 q=1 哑节点断裂单元 * 积分路径1 积分路径2 积分路径5 积分路径4 积分路径3 哑节点断裂单元 * 哑节点断裂单元 * Fracture Mechanics 华中科技大学船海学院 袁锐 Fracture Mechanics 华中科技大学船海学院 袁锐 * HUST-XD-B-0111-2014 哑节点断裂单元 * 哑节点断裂单元 前面讨论了三个基本断裂参数相关的数值计算方法： （1）应力或应变外推法计算应力强度因子； （2）等效积分区域计算J积分 （3）全域/局部虚拟裂纹扩展法和虚拟裂纹闭合法计算能量释放率 基本做法是对有限元分析的结果进行后处理，以得到所需要的断裂参数。这些方法中有限元分析和计算过程是独立的，削弱了断裂参数计算和有限元分析之间的有机联系。 * 2.1 ABAQUS自定义单元子程序UEL 哑节点断裂单元 * 对于图2.1所示的平面中两个节点间的弹簧单元而言，其ABAQUS的位移矢量和节点力矢量分别为： 节点力和节点位移存在下列关系： 2.1 ABAQUS自定义单元子程序UEL 哑节点断裂单元 * 式中，Kx和Ky分别为沿总体坐标系（X，Y）的弹簧刚度系数。 公式2.2还可以用矩阵表示为： 式中，K为弹簧的系数刚度矩阵 因此，对于线性弹簧： 由于ABAQUS是一个通用的非线性有限元程序，因此它的求解过程是迭代法。在每一个分析步，ABAQUS求解平衡方程： 式中，P为外力矢量，被ABAQUS存放在数组RHS中。 在迭代法求解中，令 式中：R为残余力矢量，问题转化为求解R=0. 根据Newon-Raphson法，其切线刚度矩阵为： 2.1 ABAQUS自定义单元子程序UEL 哑节点断裂单元 * UEL: User subroutine to define an elment Product:Abaqus/Standard ABAQUS规定的程序开头 RHS An array containing the contributions of this element to the right-hand-side vectors of the overall system of equations. For most nonlinear analysis procedures, NRHS=1 and RHS should contain the residual vector. The exception is the *STATIC, RISK procedure (“static stress analysis, section 6.2.2”), for which NRHS=2 and the first column in RHS should contain the residual vector and the second column should contain the increments of external load on the element.RHS(K1,K2) is the entry for the K1th degree of freedom of the element in the K2th right-hand-side vector. 2.1 ABAQUS自定义单元子程序UEL 哑节点断裂单元 * AMATRX An array containing the contribution of this element to the Jacobian (stiffness) or other matrix of the overall system of equations. The particular matrix required at any time depends on the entries in the LFLAGS array (see below). All nonzero entries in AMATRX should be defined, e