复合材料冲击压缩研究难点.docx

不对称锥型先进树脂基层合板的冲击和冲击后压缩实验的数值模拟摘要本篇文章阐述了一种借助离散模型来获得锥型复合材料层合板冲击和冲击后压缩实验的数值模拟的方法。三种损伤模式：基体开裂，分层，纤维断裂在该模型中被考虑到。由于应力场中铺层的不连续性现象，导致了模拟的困难。在两种实验中，对模拟结果的分析用来理解损伤和失效机制。模拟得到的力-位移曲线，分层形貌，冲击后压缩位移场，剩余压缩强度与实验结果对比发现，冲击实验模拟结果与实验结果吻合较好，剩余压缩强度预测值偏低，且依赖于过渡区网格划分的质量。本篇文章重点突出了冲击模拟时内层基体开裂的建模方法。1 绪论复合材料的冲击损伤容限在结构设计中非常重要，冲击损伤的主要特征是基体开裂，分层和纤维断裂，这三种特征通常在冲击点下方沿四周传播。此类损伤主要出现在层合板内部，外部很难检测到，[1,2]即使是低速低能冲击，剩余压缩强度也会受到严重影响。对于薄层合板，目测可视损伤对应的能量阈值对于剩余压缩强度是非常重要的损失[3,4]由于一些航空结构更倾向于关注在服役生涯中的一些工具下落和碎片造成的撞击情况，所以航空当局设定的要求和规则中包含了冲击损伤容限的概念。厚度的减小通常被用来降低复合材料结构的重量，由此导致了在层合板底层铺设树脂的方法的出现。由于材料的不连续性和连续板的局部曲率，最终的结构可能会在静载荷和疲劳载荷的作用下过早失效。[5,9]基体开裂和分层形式的损伤主要有由层合板所受的高压应力引起。在文献中，许多作者建立了复合材料冲击损伤数值模型。根据其复杂程度，这些模型可以再现部分或所有的损伤类型：基体开裂，分层和纤维断裂。最初由Ladevèze和Allix 引入的中尺度水平能够真实再现冲击损伤类型。[15]根据材料定律，失效损伤力学通常用于模拟层内损伤。基于损伤模型的弹性能量单元在许多工作中被用来模拟离散失效如分层，[17,18]一些作者也用它来模拟基体开裂并定义了基体开裂和分层的间接耦合。基于对剩余压缩强度的预测，分析模型[20,22]和有限元[23,25]模型都被提出，在一些文章中，初始冲击损伤在模拟压缩之前在模型中被引入，其他作者[24,25]则更信赖两步模型：冲击损伤先被模拟，然后再进行压缩测试的模拟。就我们所知，对于锥型复合材料，不管是冲击损伤模拟还是剩余压缩强度模拟都没有在文章中被报道过。本篇研究基于Bouvet等的研究工作，旨在提出一种离散层模型来预测冲击容限，模拟冲击和冲击后压缩。也是一种对业已提出了的锥形层压板的低速冲击响应的分析的延伸和拓展。[26]这篇文章揭示了锥型层合板与平板损伤模式的强烈一致性。在这篇文章中，离散层模型用来模拟锥型层合板的冲击和冲击后压缩，模拟的结果与实验结果对比，依次来验证模型的有效性和揭示损伤机制。2数值模拟模拟方法基于Bouvet等的研究工作，[27]他们建立了一种在微观尺度层面上对复合材料冲击损伤进行模拟的离散的三维有限元模型，后来又被推广至冲击后压缩实验的模拟。[24]在这篇文章中，模型的功能被拓展至同时模拟多层层合板的冲击损伤和冲击后压缩试验。为了达到这个目的，一种对于层合板内部耦合的网格划分方法被提出来，同时相比之下，一些材料定律也会发生变化。2.1网格划分层合板的网格划分如图1a所示，层合板按C3D8模型划分，每个单元在厚度方向上表示具有相同取向的两个层。层板单元与零厚度弹性单元连接以模拟分层和横向基体开裂。为减轻层间连接，结点网格需要在整个铺层中是均匀的，这样会导致45°和-45°层单元的扭曲（图1b）。层合板中0°和90°方向的单元尺寸是1.25×1.25。图1b建立冲击模型的难点之一在于揭示各损伤模型间复杂的交互作用，在离散层模型中，横向基体开裂和分层间的耦合关系通过层合板网格划分被建立，而它是模型中重要的建筑模块之一。基体开裂的内聚单元方向与层合板纤维方向平行，同时定义在每个体积元之间。结果，体积元在基体开裂的情况下沿层合板横向方向不连续，像这类几何不连续性可用来间接模拟基体开裂和分层间的耦合作用。下降层区域的模型如图2a所示，网格表示层板落下区域的子层压板如图2b所示。图2这个模型中，编号为2的铺层在中间中断，但编号为1和3的铺层是连续的，在较厚的截面上，层间弹性单元终止于下降层的端部。然后，一组新的弹性单元在较薄的截面上用来连接连续层。注意到在树脂楠中没有刚性单元和弹性单元，这个假设类似于假设这个区域在冲击之前已经损坏。实际上，由于固化残余应力，预应力会更高。因此，为了避免模型过于复杂和过长的计算时间，没有必要为区域建立详细的网格划分策略。在文献中，相同的假设也被用在锥型层合板受面内负载的情况中。在冲击作用的影响下，假设更为有效，因为重点在于分层的传播而不在于引发。2.2材料定律冲击模拟采用Abaqus?6.11显式/动态解算器