预测低速冲击下复合材料层合板渐进失效的有限元方法

摘要

本发明涉及复合材料损伤预测，旨在提供预测低速冲击下复合材料层合板渐进失效的有限元方法。该预测低速冲击下复合材料层合板渐进失效的有限元方法包括过程：建立含冲锤、复合材料层合板以及支撑板的低速冲击有限元模型；建立复合材料弹塑性损伤本构模型；基于ABAQUS‑VUMAT用户动态材料子程序模块，运用后向欧拉算法实现提出的弹塑性损伤本构模型；对低速冲击进行计算，进一步获得冲击力、位移、速度和加速度。本发明利用ABAQUS‑VUMAT用户子程序来数值实现所建立的将塑性和损伤联合的弹塑性损伤本构模型，该模型同时考虑塑性和材料性能退化的影响，能准确预测含塑性特征的复合材料在低速冲击下的渐进损伤失效。

说明书

技术领域

本发明是关于复合材料损伤预测领域，特别涉及预测低速冲击下复合材料层合板渐进失效的有限元方法。

背景技术

当前，复合材料正广泛应用于航空航天、风力发电、压力容器、汽车等高新技术领域。但是低速冲击损伤对复合材料的强度、刚度及使用寿命都有较大的影响，因此必须清楚的了解在低速冲击下复合材料的渐进损伤演化过程。

一些复合材料如T300/914，AS4/PEEK等，在横向和剪切加载下具有显著的非线性和塑性，复合材料的不可逆变形是由各种不同的失效机制导致的，比如纤维基体的损伤，纤维/基体的界面分离，塑性变形的累积等，此时弹性损伤本构模型预测这类复合材料的低速冲击响应并不准确，因此提出一种可以准确预测有塑性行为的复合材料层合板的低速冲击响应的弹塑性损伤本构模型尤为必要。近年来针对复合材料层合板提出的弹塑性损伤本构模型未能合理考虑面外应力的影响和合适的失效准则及损伤演化，且基本集中在复合材料的静态损伤研究，而将弹塑性损伤本构模型用来准确预测复合材料层合板在低速冲击下的损伤研究很少。

ABAQUS软件可以对基于二维HASHIN失效准则对弹性复合材料进行低速冲击下的渐进失效研究，但无法直接运用较有优势的PUCK失效准则，更无法直接对有塑性特征的复合材料层合板进行低速冲击研究。

发明内容

本发明的主要目的在于克服现有技术中的不足，提供一种能准确预测含塑性特征的复合材料在低速冲击下的渐进损伤失效的方法。为解决上述技术问题，本发明的解决方案是：

提供预测低速冲击下复合材料层合板渐进失效的有限元方法，包括下述过程：

一、建立含冲锤、复合材料层合板以及支撑板的低速冲击有限元模型；

二、建立复合材料弹塑性损伤本构模型；

三、基于(使用FORTRAN语言编写的)ABAQUS-VUMAT用户动态材料子程序模块，运用后向欧拉算法实现提出的弹塑性损伤本构模型，求解应力、应变和损伤；

四、对低速冲击进行计算，进一步获得冲击力、位移、速度和加速度；

所述过程一中，建立含冲锤、复合材料层合板以及支撑板的低速冲击有限元模型：基于ABAQUS建立冲锤、复合材料层合板、支撑板部件，分别设置材料属性和划分网格，再进行ASSEMBLY对其组装之后设置分析步和通用接触属性；

其中，复合材料层合板每层根据铺层角度进行铺层，不同铺层之间设置ABAQUS自带的双线性内聚力单元；

所述过程二具体包括下述步骤：

步骤(1)：建立含损伤复合材料层合板层内本构关系；

复合材料应力-应变本构方程：S＝C_d：E^e；

其中，为未损伤材料的有效应力；S为损伤材料的名义应力；E^p为格林-拉格朗日应变张量的塑性部分，E^e＝E-E^p为格林-拉格朗日应变张量的弹性部分；C_d(C，d₁，d₂)为含损伤材料的四阶弹性张量，C为未损伤材料的四阶弹性张量，d₁和d₂是分别对应纤维断裂和基体开裂的损伤变量；

步骤(2)：建立塑性模型，具体建立方式为：

(a)考虑面外应力的复合材料塑性流动准则为：

其中，F为屈服方程；a₆₆为描述各向异性材料与塑性相关的常数，由偏轴拉伸测试确定(其数值约为1.25左右)；f为塑性势函数；k和分别为塑性硬化应力和等效塑性应变；所述是指未损伤材料的有效应力，其中i，j用来确定应力方向；所述指层合板面内垂直于纤维方向的有效应力；所述指面外垂直于纤维的有效应力；所述指面外剪切有效应力；所述指面内剪切有效应力；

(b)塑性硬化应力k为：

其中，β和n是为满足实验硬化曲线的常数(取值分别为600MPa和0.272左右)；塑性变形假定发生在损伤材料的未损伤区域，塑性流动准则和硬化准则表达在有效应力 空间中；

(c)等效塑性应变率和塑形应变率分别为：

其中，为一致塑性因子；

步骤(3)：建立基于应变描述的PUCK失效初始判据和损伤演化准则，具体建立方式为：

(d)对于纤维拉伸和压缩，损伤初始判据为：

其中，和分别为纤维拉伸和压缩的初始失效应变；所述T，C分别指拉伸和压缩；所述E₁₁是指纤维方向应变；所述分别指纤维拉伸和压缩失效判断因子；

纤维拉伸和压缩的损伤演化准则为：

其中，所述是指纤维拉伸和压缩损伤变量；所述E₁₁是指纤维方向应变；所述>

(e)对于基体拉伸损伤失效初始判据为：

其中，所述指基体拉伸失效判断因子；所述E₂₂是指基体方向应变；所述是指基体损伤变量为零的基体初始拉伸失效应变；

基体拉伸损伤演化准则为：

其中，为基体损伤变量达到1时基体临界拉伸失效应变；所述是指基体拉伸损伤变量；

(f)对于基体压缩损伤初始判据为：

其中，N是关于失效断裂面的法向方向，T和L是关于失效断裂面的切向方向；Y^C是横向压缩强度，断裂平面上的应力S_ij(i，j＝L，T，N)由笛卡氏坐标系下的Piola-Kirchhoff应力S_ij(i，j＝1，2，3)通过旋转矩阵T(α)旋转获得，T(α)为笛卡尔坐标系到断裂面坐标系的旋转矩阵；所述S_NN是指断裂面的法向应力；S_NT，S_NL是指断裂面的切向应力，μ_NL，μ_TN为断裂面面内两个切向方向摩擦系数，θ_f为断裂面的断裂角；所述>123是指在笛卡尔坐标系下的六个Piola-Kirchhoff应力S_ij(i，j＝1，2，3)；所述S_LTN是指在断裂面坐标系下的六个Piola-Kirchhoff应力S_ij(i，j＝L，T，N)；所述T(α)^T是指T(α)的转置矩阵；所述90°是指采用角度制计量的90度；

基体压缩损伤演化准则为：

其中，所述γ_γ是指断裂面联合剪切应变；和是联合剪切应变的初始和最大应变，γ_NT和γ_NL是断裂面的剪切应变；是指基体压缩损伤变量；

所述过程三具体包括下述步骤：

步骤(4)：通过用户子程序VUMAT的SDV定义第n+1增量步的开始时的初始状态变量值，同时也是第n增量步结束时的状态变量值 在第n+1增量步开始时，VUMAT读入；

其中，所述n是指第n增量步，所述E_n是指第n增量步结束时的格林-拉格朗日总应变张量，所述是指第n增量步结束时的格林-拉格朗日塑性应变张量，所述是指第n增量步结束时的等效塑性应变，所述是指第n增量步结束时的未损伤材料的有效应力，所述k_n是指第n增量步结束时的塑性硬化应力，所述d_ij，n是指第n增量步结束时的损伤变量；

步骤(5)：VUMAT由应变增量驱动，计算试应力，将试应力代入到步骤(2)屈服方程中；该公式是步骤(2)中的公式 在第n+1增量步的特定计算；

如果F_n+1≤0，则处于弹性阶段，将所有试应力和应变更新为n+1增量步状态变量>

如果F_n+1＞0，塑性加载出现，根据后向欧拉隐式算法(采用Newton-Raphson迭代)，实现试应力到屈服面的最近点返回，将试应力赋予步骤(4)中的迭代初始条件；变量>n+1的函数，运用牛顿-拉夫森算法求解>n+1≤0，结束迭代，得到n+1增量步状态变量

步骤(6)：应变和有效应力得到更新后，根据有效应力和应变代入到步骤(3)的PUCK失效准则判断是否出现损伤，如有损伤再通过损伤演化公式获得损伤变量，再根据步骤(1)通过有效应力和损伤变量计算名义应力S_n+1；

所述过程四具体为：将过程一建立的模型主文件和过程三建立的ABAQUS-VUMAT用户子程序联合，使用ABAQUS/EXPLICT方法对低速冲击进行计算，进一步获得冲击力、位移、速度和加速度；即完成低速冲击载荷下弹塑性复合材料层合板渐进失效特性的预测。

与现有技术相比，本发明的有益效果是：

本发明利用ABAQUS-VUMAT用户子程序来数值实现所建立的将塑性和损伤联合的弹塑性损伤本构模型，该模型同时考虑塑性和材料性能退化的影响，能准确预测含塑性特征的复合材料在低速冲击下的渐进损伤失效。

对于ABAQUS内嵌的基于HASHIN失效准则的弹性损伤本构，本发明基于更准确的PUCK失效准则的弹塑性损伤本构对含塑性的复合材料层合板的低速冲击载荷下的渐进失效过程的预测更加准确。

附图说明

图1为本发明实施例复合材料层合板在低速冲击下有限元模型图。

图2为本发明对所提出的弹塑性损伤本构模型的VUMAT数值实现流程图。

图3为实施例中复合材料层合板在5J能量低速冲击下冲击力-时间数值模拟结果与实验结果对比示意图。

图4为实施例中复合材料层合板在5J能量低速冲击载荷下中心位移-时间数值模拟结果与实验结果对比示意图。

图5为实施例中复合材料层合板在10J能量低速冲击载荷下冲击力-时间数值模拟结果与实验结果对比示意图。

图6为实施例中复合材料层合板在10J能量低速冲击载荷下中心位移-时间数值模拟结果与实验结果对比示意图。

图7为实施例中复合材料层合板在5J能量低速冲击结束后基体拉伸，基体压缩，层间分层损伤云图。

图8为实施例中复合材料层合板在10J能量低速冲击结束后基体拉伸，基体压缩，层间分层损伤云图。

具体实施方式

首先需要说明的是，本发明是计算机技术在复合材料损伤预测领域的一种应用。在本发明的实现过程中，会涉及到多个软件功能模块的应用。申请人认为，如在仔细阅读申请文件、准确理解本发明的实现原理和发明目的以后，在结合现有公知技术的情况下，本领域技术人员完全可以运用其掌握的软件编程技能实现本发明。凡本发明申请文件提及的均属此范畴，申请人不再一一列举。

下面结合附图与具体实施方式对本发明作进一步详细描述：

在ABAQUS/CAF中建立包括石墨/环氧树脂复合材料层合板，冲锤以及支撑板的低速冲击有限元模型，如图1所示，平板的铺层顺序为平板大小为100×100×2.1mm，密度为ρ＝1600kg/m³，总共分为20层，每层均用减缩积分三维实体单>

层合板上下各有一个正方形支撑框架模型，框架的每边宽度为10mm，厚度为2mm，冲锤和支撑框架的刚度和泊松比分别为207.6GPa和0.26，密度为ρ＝7830kg/m³，下框架的下表面固定，上支撑框架的上表面施加350kPa的压力。

利用ABAQUS/EXPLICT计算模拟冲锤冲击层合板的过程，如图2利用用户子程序VUMAT首先判别材料点是否进入塑性，得到有效应力，根据有效应力以及应变判别单元是否进入损伤，进入损伤阶段则按损伤演化准则计算损伤变量，从而得到名义应力。

图3和图4分别为在5J冲击能量下冲击力-时间和中心位移-时间的曲线图，图5和图6分别为在10J冲击能量下的冲击力-时间和中心位移-时间的曲线图，均与试验值十分准确的吻合。在冲击力-时间曲线中冲击力为0时对应的中心位移为永久中心位移，在中心位移-时间的最低点为冲击过程中的最大位移，可见，层合板在冲击过程中的最大位移和永久位移均与试验值吻合的很好。所以本发明提出的弹塑性损伤本构模型可以准确的捕捉有塑性特征的复合材料层合板在低速冲击下的渐进损伤失效特征。

本发明在ABAQUS软件的基础上进行用户子程序的开发，提出的弹塑性损伤本构模型同时考虑塑性和材料性能退化的影响，能准确预测具有塑性特征的复合材料在低速冲击下的渐进损伤失效，为深入阐明复合材料结构的损伤失效特性、提升轻量化强度设计水平提供了技术支撑。

最后，需要注意的是，以上列举的仅是本发明的具体实施例。显然，本发明不限于以上实施例，还可以有很多变形。本领域的普通技术人员能从本发明公开的内容中直接导出或联想到的所有变形，均应认为是本发明的保护范围。