复合材料计算机模拟的组成,关于复合材料层合板结构力学性能数值仿真架构的讨论...

随着纤维增强树脂基复合材料(Fiber reinforced plastic, FRP)结构设计水平和制造能力的提升，FRP在结构上的应用由少及多、由次要变主要、由非承力到主承力、由静力覆盖到考虑冲击损伤及疲劳问题。FRP已成为“先进结构”不可或缺的部分。

科学基本原理大多源于试验、观察于试验、论证于试验，但试验费时、费力、费财。数值仿真尽管不能完全代替试验，但可为FRP结构性能评估和优化设计提供依据，从而大幅减少试验次数，提高试验的经济和时间效率。

诸多学者在FRP数值仿真领域已经获得积极的科研结果，Rose[较为全面地综述了FRP结构渐进损伤数值分析的若干问题及其研究进展。目前仿真领域的研究热点将趋向集中于两个方面，一是研究复杂载荷工况下的数值模拟。例如文献[[研究粘聚区长度对分层模拟效果的影响。Turon[数值分析了疲劳载荷下复材分层损伤的起扩展情况。

本文在多年研究工作的基础上，提出一种适用于模拟FRP综合力学性能的数值仿真架构，并总结其关键要素，最后给出多种载荷工况下的数值分析算例。

1 材料损伤的定义

FRP力学性能仿真本质上与载荷型式无关，其核心在于构造一个数值仿真架构来定量描述FRP的力学行为响应。力学行为响应是能量传递的过程，材料在吸收外载输入的能量时，其内部便伴随着损伤的起始与演化。FRP力学性能分析亦即等价于确定材料内部损伤的形成机理及其累积信息。

物理机理上，损伤的定义形式十分多样。宏观上有刚度、强度、寿命，细观上有局部硬度、透波性，微观上有裂纹密度等。损伤状态下各种性能参数与完好状态下相应参数的比值，皆可定义为损伤D

$

D = 1 - \frac{ {X\left( t \right)}}{ { {X_0}}}

$

(1)

式中:X(t)为t时刻的材料性能参数; X0为对应的初始性能参数。通常地，我们多用材料的刚度性能作为X来定义损伤。

目前的数值仿真中，该性能参数通常选择为刚度(损伤前)或者能量释放率(损伤后)。如可以将损伤定义为

$

D = 1 - \frac{ { {G^{\rm{d}}}}}{ { {G_{\rm{c}}}}}

$

(2)

式中: Gd为损伤状态时的能量释放率; Gc为临界能量释放率。

2 FRP仿真的物理架构

FRP仿真的物理架构流程图如

图 1 FRP仿真的物理架构

Figure 1 Analysis architecture of FRP mechanical behaviors simulation

2.1 物理场

物理场分为两类，一类是载荷场(代表外部信息)，另一类是材料场(代表内部信息)。

载荷场是多种载荷的线性叠加，其作用于FRP结构后，材料内部产生相应顺序的损伤累积。损伤不可逆原则意味着先期载荷相当于给后期载荷引入了不同程度的初始损伤。本征场主要包含FRP初始几何参数，它代表结构材料最基本的分布信息。损伤场与属性场都是本征场的另一类表征，属性场基于力学性能在本征场的基础上进一步丰富材料信息，损伤场基于损伤定义在本征场的基础上建立了无量纲信息场。应力应变场是整个循环的运行依据，仿真架构的核心模块集中于此。

通过力学分析，基于载荷场和力学属性场，可导出应变场。进而基于材料本构，推得应力场。材料完好本构导出的应力场配合损伤起始准则判断损伤萌生，材料损伤本构导出的应力场配合损伤扩展准则判断损伤质变。若未满足起始准则，损伤状态为0；若首次满足起始准则，损伤状态由式(1)确定；若已有损伤(包含初始损伤)且未满足损伤扩展准则，损伤状态由式(2)确定；若满足损伤扩展准则，则损伤质变为完全失效，即损伤状态为1。损伤状态的确定更新了损伤场，若损伤场出现贯穿性的失效材料分布，则整个架构停止分析；反之，则进行材料信息的更新与传递。对于完全失效的材料，将其等效为空洞，并由本征场进行拓扑信息的更新；对于部分失效的材料，由属性场结合损伤定义，对力学参数进行更新，更新完成后直接执行当前载荷步下的新一轮内循环分析。若所有材料均未损伤萌生，则由载荷场增加载荷增量步，进行新一轮的外循环分析。

2.2 核心模块

各物理场之间存在相互交错的个体关系及内外两个循环的整体关系。众多关系都会最终汇集到若干核心模块上：材料点本构、损伤起始准则、损伤演化及其扩展准则、初始材料属性场的建立(尺度效应与就位效应)。

2.2.1 材料点本构

材料区分层内本构及层间本构，层内(层间)本构主要通过应力应变(应力位移)关系进行描述。试验测得的本构曲线往往包括3个阶段：线性上升段、非线性上升段及非线性下降段。仿真模型常用的本构曲线为双线性本构，即包含线性上升段及线性下降段。本文认为这两种本构都不是真实的材料点本构，综合材料物理机理及定量描述的可行性，材料本构应该包含非线性上升段及线性下降段[，如

图 2 材料点本构示意图

Figure 2 Mechanical responses curve of material point

材料点本构中的非线性与微损伤无关