专利名称：一种分析复合材料层合板层间损伤的方法
技术领域：
本发明属于计算力学领域，具体涉及一种用于计算模拟复合材料层合板层间损伤的方法。
背景技术：
复合材料层合板已在许多领域被广泛应用，尤其在航空结构中十分常见。层合板是由多层单层板粘合在一起组成整体的结构板，每个单层板称为一个子层。根据需要可将单层板的材料主方向按不同方向和不同顺序铺设，得到不同力学性能的层合板，因此层合板的力学性能具有可设计性。层间分层是层合板主要的损伤形式之一，层间强度相对面内强度较弱使得层合板容易产生层间损伤，尤其在冲击作用下更明显。层间损伤的出现使得结构稳定性和剩余强度大大降低，并可能导致灾难性的后果。如何准确预测和分析复合材料层合板的层间损伤，对合理设计使用复合材料结构十分重要。内聚力单元法(Cohesive Zone Model s, CZM)是目前分析层合板层间损伤十分有效的有限元方法。CZM法通过引入内聚力界面单元来模拟层合板的层间胶层，在层合板各层之间，或在结构可能发生分层的位置布置内聚力界面单元(Cohesive Element)，当内聚力界面单元在周围其他单元作用下的应力值达到界面强度时层间损伤起始；当内聚力界面单元进一步承受载荷，且能量释放率达到临界值时层间损伤扩展。CZM法易于有限元程序实现，不需要给出初始损伤位置就可以预测层间损伤的起始和扩展。然而，CZM法的缺点之一是为了获得合理准确的结果，需要有限元模型划分非常精细的网格才能使计算结果足够合理。目前，大多数研究者分析层合板层间性能时都使用实体单元和内聚力界面单元来建模子层面内用实体单元模拟，在上下相邻子层的实体单元之间插入零厚度或厚度很小的内聚力界面单元来模拟层合板的层间胶层。实体单元相比板壳单元本身计算量就大，当划分网格很细的时候使得计算规模很大，这就限制了 CZM法的工程应用前景。
发明内容
为克服现有技术中存在的计算规模很大限制工程应用前景的不足，本发明提出了一种分析复合材料层合板层间损伤的方法。本发明的具体过程是步骤1，建立有限元模型；步骤2，提取各组合界面单元的结点位移和坐标；提取的过程包括a.提取结点位移；从已建立好的有限元模型中，提取出各组合界面单元主结点的位移；该主结点位移用向量形式表示为Mu = ^Uu,--,Mm , ,Μ'μ , ^.γ, ^Mv )40χ1(夏)公式(I)中，uM是组合界面单元的主结点位移这,M+症分别表示第i个主结点对沿X，y，Z方向的平动位移(i = 1,2,3,4)分别表示第i个主结点对绕着X，y轴的转动位移，上标负号表示下主结点的位移分量，加号表示上主结点的位移分量，结点编号如图2所示；所述主结点对是指同一个组合界面单元中x，y坐标相同，z坐标不同的两个主结点；b.提取结点坐标；从已建立好的有限元模型中，提取出各组合界面单元主结点的坐标；该主结点坐标用向量形式表示为
权利要求
1.一种分析复合材料层合板层间损伤的方法，其特征在于，具体步骤是 步骤I，建立有限元模型； 步骤2，提取各组合界面单元的结点位移和坐标；提取的过程包括 a.提取结点位移；从已建立好的有限元模型中，提取出各组合界面单元主结点的位移；该主结点位移用向量形式表示为 公式
全文摘要
一种分析复合材料层合板层间损伤的方法，通过推导组合界面单元的刚度方程，缩聚掉组合界面单元副结点的自由度，使组合界面单元在有限元模型中是一个整体，只有八个主结点。采用本发明的组合界面单元和有限元模型计算模拟结果较好，计算得到曲线在线性段与试验曲线基本重合，在峰值之后的下降段与试验曲线的趋势一致，层间损伤起始和损伤演化过程的模拟都符合实际。计算得到的载荷峰值62.8N，试验值63.1N，误差-0.4%。与现有技术方法相比，本发明具有计算规模小，计算效率高的特点，同时，由于组合界面单元具有厚度，因此前处理建模时相比零厚度的内聚力单元来说容易实现。
文档编号G06F17/50GK102930079SQ20121037673
公开日2013年2月13日 申请日期2012年10月8日 优先权日2012年10月8日
发明者李彪, 李亚智, 胡博海, 樊振兴, 苏杰 申请人:西北工业大学