功能梯度材料层合板内跨界面I型裂纹数值分析.doc

功能梯度材料层合板内跨界面I型裂纹数值分析 于红军，吴林志，果立成，贺启林 哈尔滨工业大学复合材料与结构研究所，哈尔滨 150001 E-mail：yuhongjun@hit.edu.cn 传真：0451电话：0451要 基于非均匀单元的有限元方法是在形成单元刚度矩阵时，采用高斯积分点处的实际材料属性，使得在同一单元内部也能够反映出实际材料的非均匀特征，因此能够更加精确、高效地模拟材料的非均匀性。利用该方法研究了功能梯度材料层板内跨层间界面的I型裂纹问题，针对各层弹性模量沿厚度方向为不同指数或线性形式的层合板，分析了裂纹位置和材料参数对应力强度因子的影响。发现裂纹位置对应力强度因子有很大影响，且裂纹尖端在层间界面处应力强度因子出现极值现象，这与弹性模量在层间界面处的导数有关。 关键词 功能梯度材料 应力强度因子 非均匀单元 材料属性随着空间坐标连续变化的功能梯度材料具有很好的耐热、耐磨、强韧等性能，它能克服传统单一或组合材料属性不连续的缺点，已经在航空航天等领域引起了人们的兴趣。材料在使用过程中失效主要来源于断裂破坏，因此功能梯度材料的断裂性能倍受关注。 上世纪80年代以来，许多学者对功能梯度材料断裂性能进行了理论研究。1988年，Delale和Erdogan[1]首先研究了无限大非均匀板中的内部裂纹问题。发现功能梯度材料中裂纹尖端的应力场仍具有逆平方根奇异性，并且泊松比对应力强度因子的影响很小。1994年，Jin和Noda[2]对比了均匀材料和功能梯度材料裂纹前沿附近的物理场，并指出：对于一般的非均匀材料，只要材料属性连续且分段可微，则非均匀材料裂纹尖端场的奇异性与均匀材料的一致。Erdogan和Wu[3]研究了不同载荷下含有垂直于边界的裂纹的各向同性功能梯度条，给出了应力强度因子随着裂纹长度等参数变化的规律。Wang[4]等考虑了含有垂直于边界裂纹的功能梯度条，他们通过把材料沿着梯度方向分成了许多均匀材料层进行分析。Kim和Paulino领导的科研组对功能梯度材料数值计算进行了系统的研究，Jin和Paulino[5]提出了将单元节点处的材料属性使用与位移相同的形函数进行插值，得到单元内部任意点处的材料信息。 本文采用基于非均匀单元的有限元方法[6]研究了含有与梯度方向平行的I型裂纹的各向同性功能梯度层合板的断裂问题，利用位移法求解了功能梯度材料层合板内跨界面I型裂纹前沿的应力强度因子，给出了裂纹尖端到达层间界面处时，应力强度因子的变化规律。 基于非均匀单元的有限元方法 有限元方法中形成单元刚度矩阵时，需要如下积分： (1) 这里为应变转换矩阵，为材料矩阵，为单元体积。本文采用20节点等参单元，如图1(a)所示。有限元计算中需要将积分离散，本文采用高斯积分，公式(1)中的单元刚度矩阵经过坐标变换后可以写成： (2) 其中，、、为等参坐标系O-中的坐标，为坐标系O-xyz变换到坐标系O-的Jacobian矩阵，为Gauss积分点，、和分别是三个坐标方向的权重值，m为单轴方向上的积分点个数。 传统有限元整个单元取同一材料属性。而这里，对于梯度材料属性，令公式(2)中材料矩阵直接取Gauss积分点处对应直角坐标系中处的实际材料属性[6]。每个单元内部含有m3个高斯积分点，同一单元内部可以模拟材料的梯度变化，如图1(b)所示。文[6]取一个单层无限长的功能梯度板，杨氏模量为指数形式，材料参数，将本方法计算出的应力强度因子值与Kadioglu[7]和Guo[8]的理论结果进行了对比，结果吻和很好，如图3所示。 梯度层合板模型及数值求解 含有垂直于纸面方向（z方向）贯穿裂纹的各向同性功能梯度板，裂纹平行于梯度方向。如图2所示（仅显示x-y平面），功能梯度材料板长为2L，厚度分别为H1和H2。厚度方向界面处的材料杨氏模量分别为、和，由于泊松比对功能梯度材料裂纹尖端应力强度因子的影响很小[1]，因此泊松比取为常数。裂纹尖端的x坐标为分别为a和b。定义材料梯度参数仅为x的函数，如公式(3)所示。 (3) 本文所研究的层板厚度的无限长板，取平面应变情况。长度与厚度的比值，它足以模拟无限长板。定义无量纲裂纹中心和裂纹长度分别为和，及无量纲应力强度因子，其中，为裂纹面上均匀压载荷。 由于本程序为三维程序，因此在垂直于纸面方向（z方向）取一个单元，限制所有节点的z方向自由度。裂纹前沿布置应力具有奇异性的楔形单元，如图4所示。使用位移法计算裂纹前沿的应力强度因子。 取，