一种考虑复合材料参数三维各向异性非线性的有限元数值模拟方法技术

本发明专利技术提出一种考虑复合材料参数三维各向异性非线性的有限元数值模拟方法，利用ABAQUS子程序UMAT，在复合材料数值模拟过程中，对每一迭代步复合材料的各向异性参数进行计算更新，由此对复合材料在数值计算过程中各项异性参数的非线性做出模拟。本发明专利技术充分利用了有限元计算分析过程，将复合材料的有限元计算过程二次分解，对每一次迭代步计算的单元积分点的应力应变加以提取计算，由子程序不断计算更新每一次迭代步后材料属性的变化，包括弹性模量、泊松比等，弥补了目前数值模拟方法不能模拟复合材料的弹性参数在受力变形过程中出现各向异性非线性的缺陷。

A Finite Element Numerical Simulation Method Considering Three-Dimensional Anisotropic Nonlinearity of Composite Material Parameters

The present invention proposes a finite element numerical simulation method considering the three-dimensional anisotropic nonlinearity of composite material parameters. In the process of composite material numerical simulation, the anisotropic parameters of each iteration step are updated by using the ABAQUS subroutine UMAT, thus the nonlinearity of the anisotropic parameters of composite materials in the numerical calculation process is simulated. The invention makes full use of the finite element calculation and analysis process, decomposes the finite element calculation process of composite material twice, extracts and calculates the stress and strain of the element integral points calculated in each iteration step, calculates and updates the change of material properties after each iteration step by subroutine continuously, including elastic modulus, Poisson's ratio, etc., which makes up for the fact that the current numerical simulation method can not simulate the complex. The elastic parameters of composite materials are anisotropic and non-linear in the process of deformation.

【技术实现步骤摘要】
一种考虑复合材料参数三维各向异性非线性的有限元数值模拟方法
本专利技术属于复合材料的材料参数设置
，具体涉及对复合材料进行的有限元数值模拟方法。
技术介绍
随着科学技术的进步与发展，复合材料凭借其较高的强度-重量比、刚度-重量比、较好的抗疲劳性能与抗腐蚀性能，已经在很多领域开始崭露头角。目前在结构工程中使用复合材料的比例越来越高，各种各样的复合材料结构逐渐取代原先的钢材结构。但是研究表明，复合材料的材料参数以及力学性能远远比钢材更为复杂，材料参数的不稳定性造成了复合材料力学性能的不稳定。因此工程中在对复合材料结构进行有限元强度计算校核时，经常会高估复合材料的力学性能，给实际的工程应用带来隐患。再考虑到复合材料价格昂贵，因此对复合材料力学性能的高估，不仅会导致巨大工程隐患，更会造成不可估计的人员、财产损失。目前对复合材料结构进行结构强度分析时，主要采用数值模拟法、实验法、理论法。实验方法仅仅是针对简单的复合材料，对其进行拉伸、弯曲、冲击等实验操作，并通过大量实验，总结实验数据，拟合出相对可靠的经验公式。理论方法目前大多针对结构简单、形状常规的物体，进行理论推导计算结构响应。数值模本文档来自技高网...

【技术保护点】
1.一种考虑复合材料参数三维各向异性非线性的有限元数值模拟方法，其特征在于，包括以下步骤：步骤S1：利用ABAQUS子程序UMAT，调用主程序ABAQUS中当前迭代步下的积分点应力、应变与应变增量等参数；步骤S2：根据上述读取的材料应力应变，利用基体材料应力应变关系式，计算出此时基体的弹性模量，同时带入复合材料的各向异性参数计算公式，计算出当前迭代步下的复合材料的各项异性参数；步骤S3：由上述求得的当前迭代步下的各向异性参数，组装更新当前迭代步下复合材料的应力应变关系式；步骤S4：根据步骤S1读取的应变增量，以及步骤S3得到的应力应变关系式，计算求解当前迭代步下的应力增量；步骤S5：根据步骤S...

【技术特征摘要】
1.一种考虑复合材料参数三维各向异性非线性的有限元数值模拟方法，其特征在于，包括以下步骤：步骤S1：利用ABAQUS子程序UMAT，调用主程序ABAQUS中当前迭代步下的积分点应力、应变与应变增量等参数；步骤S2：根据上述读取的材料应力应变，利用基体材料应力应变关系式，计算出此时基体的弹性模量，同时带入复合材料的各向异性参数计算公式，计算出当前迭代步下的复合材料的各项异性参数；步骤S3：由上述求得的当前迭代步下的各向异性参数，组装更新当前迭代步下复合材料的应力应变关系式；步骤S4：根据步骤S1读取的应变增量，以及步骤S3得到的应力应变关系式，计算求解当前迭代步下的应力增量；步骤S5：根据步骤S4中的应力增量，结合步骤S1中读取的数据，更新应力、应变，最终将更新的应力、应变、应变增量带入到主程序ABAQUS进行收敛测试；在迭代步中完成以上5个步骤后，将当前迭代步的计算结果在主程序ABAQUS中进行收敛测试，如果当前迭代步的结果收敛，则该增量步结束，进入到下一个增量步；否则继续到下一个迭代步，重复上述步骤S1-步骤S5，如果迭代次数超过主程序最大迭代次数，则终止计算。2.根据权利要求1所述的有限元数值模拟方法，其特征在于：所述步骤S2中，基体材料应...

【专利技术属性】
技术研发人员：王树青，姚潞，何文涛，孟祥剑，张常友，魏慧宾，
申请(专利权)人：中国海洋大学，
类型：发明
国别省市：山东,37