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

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

【技术实现步骤摘要】

本专利技术是关于复合材料损伤预测领域，特别涉及预测低速冲击下复合材料层合板渐进失效的有限元方法。
技术介绍
当前，复合材料正广泛应用于航空航天、风力发电、压力容器、汽车等高新
但是低速冲击损伤对复合材料的强度、刚度及使用寿命都有较大的影响，因此必须清楚的了解在低速冲击下复合材料的渐进损伤演化过程。一些复合材料如T300/914，AS4/PEEK等，在横向和剪切加载下具有显著的非线性和塑性，复合材料的不可逆变形是由各种不同的失效机制导致的，比如纤维基体的损伤，纤维/基体的界面分离，塑性变形的累积等，此时弹性损伤本构模型预测这类复合材料的低速冲击响应并不准确，因此提出一种可以准确预测有塑性行为的复合材料层合板的低速冲击响应的弹塑性损伤本构模型尤为必要。近年来针对复合材料层合板提出的弹塑性损伤本构模型未能合理考虑面外应力的影响和合适的失效准则及损伤演化，且基本集中在复合材料的静态损伤研究，而将弹塑性损伤本构模型用来准确预测复合材料层合板在低速冲击下的损伤研究很少。ABAQUS软件可以对基于二维HASHIN失效准则对弹性复合材料进行低速冲击下的渐进失效研究，但无法直接运用较有优势的PUCK失效准则，更无法直接对有塑性特征的复合材料层合板进行低速冲击研究。
技术实现思路
本专利技术的主要目的在于克服现有技术中的不足，提供一种能准确预测含塑性特征的复合材料在低速冲击下的渐进损伤失效的方法。为解决上述技术问题，本专利技术的解决方案是：提供预测低速冲击下复合材料层合板渐进失效的有限元方法，包括下述过程：一、建立含冲锤、复合材料层合板以及支撑板的低速冲击有限元模型；二、建立复合材料弹塑性损伤本构模型；三、基于(使用FORTRAN语言编写的)ABAQUS-VUMAT用户动态材料子程序模块，运用后向欧拉算法实现提出的弹塑性损伤本构模型，求解应力、应变和损伤；四、对低速冲击进行计算，进一步获得冲击力、位移、速度和加速度；所述过程一中，建立含冲锤、复合材料层合板以及支撑板的低速冲击有限元模型：基于ABAQUS建立冲锤、复合材料层合板、支撑板部件，分别设置材料属性和划分网格，再进行ASSEMBLY对其组装之后设置分析步和通用接触属性；其中，复合材料层合板每层根据铺层角度进行铺层，不同铺层之间设置ABAQUS自带的双线性内聚力单元；所述过程二具体包括下述步骤：步骤(1)：建立含损伤复合材料层合板层内本构关系；复合材料应力-应变本构方程：S＝Cd：Ee；其中，为未损伤材料的有效应力；S为损伤材料的名义应力；Ep为格林-拉格朗日应变张量的塑性部分，Ee＝E-Ep为格林-拉格朗日应变张量的弹性部分；Cd(C，d1，d2)为含损伤材料的四阶弹性张量，C为未损伤材料的四阶弹性张量，d1和d2是分别对应纤维断裂和基体开裂的损伤变量；步骤(2)：建立塑性模型，具体建立方式为：(a)考虑面外应力的复合材料塑性流动准则为：其中，F为屈服方程；a66为描述各向异性材料与塑性相关的常数，由偏轴拉伸测试确定(其数值约为1.25左右)；f为塑性势函数；k和分别为塑性硬化应力和等效塑性应变；所述是指未损伤材料的有效应力，其中i，j用来确定应力方向；所述指层合板面内垂直于纤维方向的有效应力；所述指面外垂直于纤维的有效应力；所述指面外剪切有效应力；所述指面内剪切有效应力；(b)塑性硬化应力k为：其中，β和n是为满足实验硬化曲线的常数(取值分别为600MPa和0.272左右)；塑性变形假定发生在损伤材料的未损伤区域，塑性流动准则和硬化准则表达在有效应力空间中；(c)等效塑性应变率和塑形应变率分别为：其中，为一致塑性因子；步骤(3)：建立基于应变描述的PUCK失效初始判据和损伤演化准则，具体建立方式为：(d)对于纤维拉伸和压缩，损伤初始判据为：其中，和分别为纤维拉伸和压缩的初始失效应变；所述T，C分别指拉伸和压缩；所述E11是指纤维方向应变；所述分别指纤维拉伸和压缩失效判断因子；纤维拉伸和压缩的损伤演化准则为：其中，所述是指纤维拉伸和压缩损伤变量；所述E11是指纤维方向应变；所述是指纤维损伤变量达到1的纤维临界拉伸和压缩失效应变；所述是指纤维损伤变量为零的纤维初始拉伸和压缩失效应变；(e)对于基体拉伸损伤失效初始判据为：其中，所述指基体拉伸失效判断因子；所述E22是指基体方向应变；所述是指基体损伤变量为零的基体初始拉伸失效应变；基体拉伸损伤演化准则为：其中，为基体损伤变量达到1时基体临界拉伸失效应变；所述是指基体拉伸损伤变量；(f)对于基体压缩损伤初始判据为：其中，N是关于失效断裂面的法向方向，T和L是关于失效断裂面的切向方向；YC是横向压缩强度，断裂平面上的应力Sij(i，j＝L，T，N)由笛卡氏坐标系下的Piola-Kirchhoff应力Sij(i，j＝1，2，3)通过旋转矩阵T(α)旋转获得，T(α)为笛卡尔坐标系到断裂面坐标系的旋转矩阵；所述SNN是指断裂面的法向应力；SNT，SNL是指断裂面的切向应力，μNL，μTN为断裂面面内两个切向方向摩擦系数，θf为断裂面的断裂角；所述指基体压缩失效判断因子；所述是指笛卡尔坐标系下面内剪切强度；所述是指在断裂平面内的横向剪切强度；所述S123是指在笛卡尔坐标系下的六个Piola-Kirchhoff应力Sij(i，j＝1，2，3)；所述SLTN是指在断裂面坐标系下的六个Piola-Kirchhoff应力Sij(i，j＝L，T，N)；所述T(α)T是指T(α)的转置矩阵；所述90°是指采用角度制计量的90度；基体压缩损伤演化准则为：其中，所述γγ是指断裂面联合剪切应变；和是联合剪切应变的初始和最大应变，γNT和γNL是断裂面的剪切应变；是指基体压缩损伤变量；所述过程三具体包括下述步骤：步骤(4)：通过用户子程序VUMAT的SDV定义第n+1增量步的开始时的初始状态变量值，同时也是第n增量步结束时的状态变量值在第n+1增量步开始时，VUMAT读入；其中，所述n是指第n增量步，所述En是指第n增量步结束时的格林-拉格朗日总应变张量，所述是指第n增量步结束时的格林-拉格朗日塑性应变张量，所述是指第n增量步结束时的等效塑性应变，所述是指第n增量步结束时的未损伤材料的有效应力，所述kn是指第n增量步结束时的塑性硬化应力，所述dij，n是指第n增量步结束时的损伤变量；步骤(5)：VUMAT由应变增量驱动，计算试应力，将试应力代入到步骤(2)屈服方程中；该公式是步骤(2)中的公式在第n+1增量步的特定计算；如果Fn+1≤0，则处于弹性阶段，将所有试应力和应变更新为n+1增量步状态变量如果Fn+1＞0，塑性加载出现，根据后向欧拉隐式算法(采用Newton-Raphson迭代)，实现试应力到屈服面的最近点返回，将试应力赋予步骤(4)中的迭代初始条件；变量是塑性一致因子的增量Δλn+1的函数，运用牛顿-拉夫森算法求解再更新应变和有效应力，直到Fn+1≤0，结束迭代，得到n+1增量步状态变量步骤(6)：应变和有效应力得到更新后，根据有效应力和应变代入到步骤(3)的PUCK失效准则判断是否出现损伤，如有损伤再通过损伤演化公式获得损伤变量，再根据步骤(1)通过有效应力和损伤变量计算名义应力Sn+1；所述过程四具体为：将过程一建立的模本文档来自技高网...

【技术保护点】
预测低速冲击下复合材料层合板渐进失效的有限元方法，其特征在于，包括下述过程：一、建立含冲锤、复合材料层合板以及支撑板的低速冲击有限元模型；二、建立复合材料弹塑性损伤本构模型；三、基于ABAQUS‑VUMAT用户动态材料子程序模块，运用后向欧拉算法实现提出的弹塑性损伤本构模型，求解应力、应变和损伤；四、对低速冲击进行计算，进一步获得冲击力、位移、速度和加速度；所述过程一中，建立含冲锤、复合材料层合板以及支撑板的低速冲击有限元模型：基于ABAQUS建立冲锤、复合材料层合板、支撑板部件，分别设置材料属性和划分网格，再进行ASSEMBLY对其组装之后设置分析步和通用接触属性；其中，复合材料层合板每层根据铺层角度进行铺层，不同铺层之间设置ABAQUS自带的双线性内聚力单元；所述过程二具体包括下述步骤：步骤(1)：建立含损伤复合材料层合板层内本构关系；复合材料应力‑应变本构方程：S＝Cd：Ee：其中，为未损伤材料的有效应力；S为损伤材料的名义应力；Ep为格林‑拉格朗日应变张量的塑性部分，Ee＝E‑Ep为格林‑拉格朗日应变张量的弹性部分；Cd(C，d1，d2)为含损伤材料的四阶弹性张量，C为未损伤材料的四阶弹性张量，d1和d2是分别对应纤维断裂和基体开裂的损伤变量；步骤(2)：建立塑性模型，具体建立方式为：(a)考虑面外应力的复合材料塑性流动准则为：其中，F为屈服方程；a66为描述各向异性材料与塑性相关的常数，由偏轴拉伸测试确定；f为塑性势函数；k和分别为塑性硬化应力和等效塑性应变；所述是指未损伤材料的有效应力，其中i，j用来确定应力方向；所述指层合板面内垂直于纤维方向的有效应力；所述指面外垂直于纤维的有效应力；所述指面外剪切有效应力；所述指面内剪切有效应力；(b)塑性硬化应力k为：其中，β和n是为满足实验硬化曲线的常数；塑性变形假定发生在损伤材料的未损伤区域，塑性流动准则和硬化准则表达在有效应力空间中；(c)等效塑性应变率和塑形应变率分别为：其中，为一致塑性因子；步骤(3)：建立基于应变描述的PUCK失效初始判据和损伤演化准则，具体建立方式为：(d)对于纤维拉伸和压缩，损伤初始判据为：其中，和分别为纤维拉伸和压缩的初始失效应变；所述T，C分别指拉伸和压缩；所述E11是指纤维方向应变；所述分别指纤维拉伸和压缩失效判断因子；纤维拉伸和压缩的损伤演化准则为：其中，所述是指纤维拉伸和压缩损伤变量；所述E11是指纤维方向应变；所述是指纤维损伤变量达到1的纤维临界拉伸和压缩失效应变；所述是指纤维损伤变量为零的纤维初始拉伸和压缩失效应变；(e)对于基体拉伸损伤失效初始判据为：其中，所述指基体拉伸失效判断因子；所述E22是指基体方向应变；所述是指基体损伤变量为零的基体初始拉伸失效应变；基体拉伸损伤演化准则为：其中，为基体损伤变量达到1时基体临界拉伸失效应变；所述是指基体拉伸损伤变量；(f)对于基体压缩损伤初始判据为：其中，N是关于失效断裂面的法向方向，T和L是关于失效断裂面的切向方向；YC是横向压缩强度，断裂平面上的应力Sij(i，j＝L，T，N)由笛卡氏坐标系下的Piola‑Kirchhoff应力Sij(i，j＝1，2，3)通过旋转矩阵T(α)旋转获得，T(α)为笛卡尔坐标系到断裂面坐标系的旋转矩阵；所述SNN是指断裂面的法向应力；SNT，SNL是指断裂面的切向应力，μNL，μTN为断裂面面内两个切向方向摩擦系数，θf为断裂面的断裂角；所述指基体压缩失效判断因子；所述是指笛卡尔坐标系下面内剪切强度；所述是指在断裂平面内的横向剪切强度；所述S123是指在笛卡尔坐标系下的六个Piola‑Kirchhoff应力Sij(i，j＝1，2，3)；所述SLTN是指在断裂面坐标系下的六个Piola‑Kirchhoff应力Sij(i，j＝L，T，N)；所述T(α)T是指T(α)的转置矩阵；所述90°是指采用角度制计量的90度；基体压缩损伤演化准则为：其中，所述γγ是指断裂面联合剪切应变；和是联合剪切应变的初始和最大应变，γNT和γNL是断裂面的剪切应变；是指基体压缩损伤变量；所述过程三具体包括下述步骤：步骤(4)：通过用户子程序VUMAT的SDV定义第n+1增量步的开始时的初始状态变量值，同时也是第n增量步结束时的状态变量值在第n+1增量步开始时，VUMAT读入；其中，所述n是指第n增量步，所述En是指第n增量步结束时的格林‑拉格朗日总应变张量，所述是指第n增量步结束时的格林‑拉格朗日塑性应变张量，所述是指第n增量步结束时的等效塑性应变，所述是指第n增量步结束时的未损伤材料的有效应力，所述kn是指第n增量步结束时的塑性硬化应力，所述dij，n是指第n增量步结束时的损伤变量；步骤(5)：VUMAT由应变增量驱动，计算试应力，将试应力代入到步骤(2)屈服方程中；该公式是步骤...

【技术特征摘要】
1.预测低速冲击下复合材料层合板渐进失效的有限元方法，其特征在于，包括下述过程：一、建立含冲锤、复合材料层合板以及支撑板的低速冲击有限元模型；二、建立复合材料弹塑性损伤本构模型；三、基于ABAQUS-VUMAT用户动态材料子程序模块，运用后向欧拉算法实现提出的弹塑性损伤本构模型，求解应力、应变和损伤；四、对低速冲击进行计算，进一步获得冲击力、位移、速度和加速度；所述过程一中，建立含冲锤、复合材料层合板以及支撑板的低速冲击有限元模型：基于ABAQUS建立冲锤、复合材料层合板、支撑板部件，分别设置材料属性和划分网格，再进行ASSEMBLY对其组装之后设置分析步和通用接触属性；其中，复合材料层合板每层根据铺层角度进行铺层，不同铺层之间设置ABAQUS自带的双线性内聚力单元；所述过程二具体包括下述步骤：步骤(1)：建立含损伤复合材料层合板层内本构关系；复合材料应力-应变本构方程：S＝Cd：Ee：其中，为未损伤材料的有效应力；S为损伤材料的名义应力；Ep为格林-拉格朗日应变张量的塑性部分，Ee＝E-Ep为格林-拉格朗日应变张量的弹性部分；Cd(C，d1，d2)为含损伤材料的四阶弹性张量，C为未损伤材料的四阶弹性张量，d1和d2是分别对应纤维断裂和基体开裂的损伤变量；步骤(2)：建立塑性模型，具体建立方式为：(a)考虑面外应力的复合材料塑性流动准则为：其中，F为屈服方程；a66为描述各向异性材料与塑性相关的常数，由偏轴拉伸测试确定；f为塑性势函数；k和分别为塑性硬化应力和等效塑性应变；所述是指未损伤材料的有效应力，其中i，j用来确定应力方向；所述指层合板面内垂直于纤维方向的有效应力；所述指面外垂直于纤维的有效应力；所述指面外剪切有效应力；所述指面内剪切有效应力；(b)塑性硬化应力k为：其中，β和n是为满足实验硬化曲线的常数；塑性变形假定发生在损伤材料的未损伤区域，塑性流动准则和硬化准则表达在有效应力空间中；(c)等效塑性应变率和塑形应变率分别为：其中，为一致塑性因子；步骤(3)：建立基于应变描述的PUCK失效初始判据和损伤演化准则，具体建立方式为：(d)对于纤维拉伸和压缩，损伤初始判据为：其中，和分别为纤维拉伸和压缩的初始失效应变；所述T，C分别指拉伸和压缩；所述E11是指纤维方向应变；所述分别指纤维拉伸和压缩失效判断因子；纤维拉伸和压缩的损伤演化准则为：其中，所述是指纤维拉伸和压缩损伤变量；所述E11是指纤维方向应变；所述是指纤维损伤变量达到1的纤维临界拉伸和压缩失效应变；所述是指纤维损伤变量为零的纤维初始拉伸和压缩失效应变；(e)对于基体拉伸损伤失效初始判据为：其中，所述指基体拉伸失效判断因子；所述E22是指基体方向应变；所述是指基体损伤变量为零的基体初始拉伸失效应变；基体拉伸损伤演化准则为：其中，为基体损伤变量达到1时基体临界拉伸失效应变；所述是指基体拉伸损伤变量；(f)对于基体压缩损伤...

【专利技术属性】
技术研发人员：廖斌斌，刘鹏飞，郑津洋，
申请(专利权)人：浙江大学，
类型：发明
国别省市：浙江;33