专利名称：一种含冲击损伤复合材料层压板剩余压缩强度的测算方法
技术领域：
本发明属于复合材料技术，涉及一种含冲击损伤复合材料层压板剩余压缩强度的测算方法。
背景技术：
文献[Hiroshi Suemasu, Wataru Sasaki. A Numerical Study on Compressive Behavior of Composite Plates with Multiple Circular Delaminations Considering Delamination Propagation[J]. Composite Science and Technology，2008，68(12) 2562-2567.]比较典型的代表了目前估算含冲击损伤复合材料层压板剩余压缩强度的方法。该方法是通过建立有限元模型方法来计算的。具体做法是I)根据设想的分层损伤，将层压板划分为若干子层。每层子层用体元(solid element)模拟；2)在子层与子层之间使用粘接单元(cohesive element)或虚拟裂纹闭合技术 (Virtual Crack Closure Technique, VCCT)或其他方法,模拟冲击造成的分层损伤；3)通过模拟分层损伤的扩展来计算层压板的剩余压缩强度。然而使用这种方法存在如下缺点I)冲击造成的分层损伤几乎存在于任意相邻两铺层之间，且形状不规则。要在有限元模型中模拟这些分层损伤，是非常困难的；2)在层压板模型内部设想的位置模拟冲击损伤，这与实际的分层损伤情况差别很大，影响了计算结果的可信性；3)按照这种方法建立的有限元模型，需要同时考虑几何非线性、结构稳定性、接触问题和破坏扩展，导致计算结果收敛困难，计算不易成功。
发明内容
本发明的目的提供一种能真实反映分层损伤分布情况、建模方便、计算速度快的含冲击损伤复合材料层压板剩余压缩强度的测算方法。本发明的技术方案是一种含冲击损伤复合材料层压板剩余压缩强度的测算方法，其包括如下步骤步骤I:无损检测对含冲击损伤的复合材料层压板受到冲击的表面，进行无损检测，获取该处损伤在层压板面内的位置信息和在层压板厚度方向上的位置信息；步骤2 :建立层压板有限元模型建立与被检测复合材料层压板相同构型且与压缩破坏试验相同边界条件的有限元模型；步骤3 :将无损检测信息输入有限元模型根据损伤的位置信息在有限元模型中找到对应的位置，将其所处的单元定义为含损伤单元，将其所处的铺层定义为含损伤铺层，从而将无损检测信息输入了有限元模型；步骤4:软化材料所述有限元模型只含单一铺层，将含损伤铺层的材料弹性常数乘以一个折减系数 Rw，使其达到材料软化，降低其承载能力的作用，定义折减系数为
权利要求
1.一种含冲击损伤复合材料层压板剩余压缩强度的测算方法，其特征在于，包括如下步骤步骤I:无损检测对含冲击损伤的复合材料层压板受到冲击的表面，进行无损检测，获取该处损伤在层压板面内的位置信息和在层压板厚度方向上的位置信息；步骤2 :建立层压板有限元模型建立与被检测复合材料层压板相同构型且与压缩破坏试验相同边界条件的有限元模型；步骤3 :将无损检测信息输入有限元模型根据损伤的位置信息在有限元模型中找到对应的位置，将其所处的单元定义为含损伤单元，将其所处的铺层定义为含损伤铺层，从而将无损检测信息输入了有限元模型；步骤4 :软化材料所述有限元模型只含单一铺层，将含损伤铺层的材料弹性常数乘以一个折减系数#' 使其达到材料软化，降低其承载能力的作用，定义折减系数为其中N(J)为含损伤铺层j的剩余压缩承载能力，tj为该铺层厚度，由层压板制造者提供，σ。 为该层压板的无损伤压缩强度；步骤5 :计算该有限元模型，得到剩余压缩强度值。
2.根据权利要求I所述的含冲击损伤复合材料层压板剩余压缩强度的测算方法，其特征在于步骤2中所建的有限元模型采用壳单元构建，壳单元的截面属性为复合材料属性， 即其材料体系参数、铺层数、铺层角度、铺层厚度、铺层顺序与被检测层压板相同。
3.根据权利要求2所述的含冲击损伤复合材料层压板剩余压缩强度的测算方法，其特征在于材料体系的破坏判据选用Hashin判据，所述Hashin判据包括四种破坏模式纤维拉伸破坏、纤维压缩破坏、基体拉伸破坏和机体压缩破坏。
4.根据权利要求3所述的含冲击损伤复合材料层压板剩余压缩强度的测算方法，其特征在于步骤5中使用牛顿迭代法求解有限元模型，提取该模型承受的最大压缩载荷，除以施加压缩载荷的面积，便得到该层压板的剩余压缩强度计算值。
全文摘要
本发明属于复合材料技术，涉及一种含冲击损伤复合材料层压板剩余压缩强度的测算方法。本发明先对含冲击损伤的复合材料层压板受到冲击的表面，进行超声C扫描无损检测，然后建立与被检测复合材料层压板相同构型且与压缩破坏试验相同边界条件的有限元模型，将无损检测信息输入有限元模型，通过将含损伤铺层的材料弹性常数乘以一个折减系数R(j)，使其达到材料软化，降低其承载能力的作用，最后计算该有限元模型，得到剩余压缩强度值。本发明测算方法能真实反映分层损伤分布情况、建模方便、计算速度快，具有较大的实际应用价值。
文档编号G01N3/00GK102607935SQ20111044269
公开日2012年7月25日 申请日期2011年12月27日 优先权日2011年12月27日
发明者杨宇, 杨胜春, 柴亚南, 沈真, 肖迎春 申请人:中国飞机强度研究所