基于ANSYS 的穿透裂纹三维板数值分析

2012-08-01李成英

山西建筑 2012年29期

高健李成英

(辽宁工业大学，辽宁锦州 121000)

0 引言

裂纹广泛存在于工程构件中，其对工程结构的破坏作用越来越受到重视，各国学者在大量实验、分析和数据基础上建立了一些裂纹模型和判定准则，在工程中得到了比较广泛的应用。随着近代计算机技术的发展，有限元方法逐渐渗透到了构件裂纹分析当中。带裂纹的构件在外力作用下，材料内部会产生应力集中，裂尖附近材料发生塑性屈服[1]，并在构件最薄弱的裂尖处首先开裂，使构件在低于实际强度很多的情况下发生破坏，所以研究含裂纹材料的开裂条件是问题的关键，由此引入了断裂力学中应力强度因子的概念。

应力强度因子K反映了裂纹尖端应力场的强弱，通过应力强度因子与材料断裂韧性KIC的比较作为裂纹扩展发生的判据。各国对应力强度因子的研究不断深入，我国航空研究院编制的应力强度因子手册对各种规律裂纹形式给出了较为精确的解析解。本文以不同位置穿透裂纹为例，运用ANSYS中KCALC命令和J积分命令流程序来模拟推导应力强度因子计算。

1 穿透裂纹板模型

矩形板尺寸 2a=254，2b=304.8，c=12.7，裂纹 2t=50.8，受压荷载 σ=0.01 MPa，弹性模量 E=24 GPa，泊松比 γ=0.2。图1表示穿透裂纹板建立的位置、尺寸及荷载信息。

图1 中心穿透裂纹矩形板

2 裂纹板建模

在线弹性问题中，裂纹尖端附近的应力场可统一为:

且有:

其中，r为任一点到裂纹尖端距离;σ为该点应力张量，根据式(1)当r趋于0时，对应的σ趋于无穷大，所以裂纹尖端处的应力与应变存在奇异性。ANSYS在进行模拟时，为了选取应变奇异点，相应的裂纹面需与其一致，即围绕裂纹顶点的有限元单元应选取带中间节点的二次奇异单元，并将其中节点放到距尖端1/4 处[2]。

ANSYS建模主要思路:对三维八节点实体Solid65单元通过添加1/4位置的中间节点转化成Solid95单元，编辑宏程序*CREATE，FRACT，MAC[3]来实现。穿透裂纹矩形板模型为 X，Y 轴双向对称，故取1/4建模分析即可。划分网格时对裂纹周围需采用特殊划分处理。

建模具体过程如下[4]:

1)在裂纹前缘上，单元的大小取决于局部曲率的数值。本文采用沿圆环状弯曲裂纹前缘，在15°～30°内设置一个单元。裂纹尖端周围单元建立如图2～图5所示。

图2 裂尖一层单元前视图

图3 裂尖一层单元等轴侧视图

图4 裂尖沿圆环多层单元前视图

图5 裂尖沿圆环多层单元等轴侧视图

2)剩余节点中每8个相邻节点形成一个单元，将1/4板模型建立完整，如图6，图7所示。

图6 1/4板单元平视图

图7 1/4板单元等轴侧视图

图8 1/4板单元施加约束及荷载平视图

图9 1/4板单元施加约束及荷载等轴侧视图

3)1/4板单元建立后施加双向对称约束，并对上边和右边施加UX，UY，UZ简支约束。对裂纹处不施加任何约束来模拟穿透裂纹的力学特性。板表面均匀压力荷载及约束如图8，图9所示。

3 计算方法及理论依据

3.1 J积分

J积分[5]为围绕裂纹尖端的围线积分，这个积分值与积分路径无关，为一常数，并反映了裂纹尖端应力应变场强度的量。

J积分表达式:

式中:Γ——围绕裂纹尖端的一条任意反时针回路，起端始于裂纹下表面，末端终于裂纹的上表面;

W——弹塑性条件下，回路Γ上任意点(x，y)的应变能密

Ti——回路Γ上任意点(x，y)处的应力分量;

ui——回路Γ上任意点(x，y)处的位移分量;

ds——回路Γ上的弧元，在ANSYS中通过编写APDL命令流来完成任意路径下的积分求值过程。

3.2 KCALC 命令

ANSYS软件本身具有计算断裂力学参数的功能。可以采用后处理器POST1中的KCALC命令(Main Menu＞General Postproc＞Nodal Calcs＞ Stress Int Factr)或编写 KCALC，KPLAN，MAT，KCSYM，KLOCPR命令流来计算复合型断裂模式中的应力强度因子(KⅠ，KⅡ，KⅢ)。后者操作更为简便。