[发明专利]一种基于细观力学模型的横观各向同性纤维材料性能的预测方法

说明书

技术领域

本发明涉及横观各向同性纤维材料性能的预测，适用于复合材料细观力学方面的研究，具体涉及一种基于细观力学模型的横观各向同性纤维材料性能的预测方法。

背景技术

复合材料由于具有比强度、比刚度高及性能可设计等优点而越来越广泛地应用于航空航天飞行器结构中，其应用范围已经从最初的非承力结构、次承力结构扩展到主承力结构。研究表明，在设计与分析复合材料力学结构时，对于复合材料自身的力学性能的研究是十分重要的。

在研究复合材料力学性能的方法中，除了通过试验的手段来测得外，还有学者提出了基于细观力学的方法来预测复合材料的材料性能，这其中又包含了理论公式法以及基于代表性体积单元的有限元方法。然而，通过这两种基于细观力学的方法预测复合材料的材料性能，均需要已知复合材料的细观材料性能，对于纤维增强复合材料而言，即纤维与基体的材料性能。然而，由于纤维的尺寸过小，难以通过试验的方法获得纤维的材料性能，特别当纤维不是各向同性材料时。

随着复合材料越来越多地应用于飞机结构，基于细观力学的损伤及失效分析受到了越来越多的注视。对于纤维增强复合材料而言，在进行基于细观力学的损伤及失效分析时，需要考虑组分各自的退化，因此纤维与基体的材料性能是必须的。另外，复合材料的细观材料性能对于材料的热弹性应力分析也是必要的。

发明内容

本发明要解决的技术问题为：克服现有技术的不足，以复合材料基本材料性能、基体的材料性能以及纤维长度方向的弹性模量为输入参数，通过细观力学理论模型反演，得到横观各向同性纤维的材料性能的初值，再利用基于代表性体积单元的有限元优化方法计算得到横观各向同性纤维的材料性能，给出了一种求解横观各向同性纤维材料性能的优化方法。

本发明解决上述技术问题采用的技术方案为：一种基于细观力学模型的横观各向同性纤维材料性能的预测方法，实现步骤如下：

步骤A，根据复合材料的基本材料性能参数、基体的材料性能以及纤维长度方向的弹性模量，利用细观力学理论模型，反演得到横观各向同性纤维材料性能参数的初值；

步骤B，建立基于细观力学的代表性体积单元有限元模型，带入纤维材料性能参数的初值以及基体的材料性能，施加周期性边界条件及四种载荷；

步骤C，基于步骤B中的有限元模型，预测复合材料的材料属性，获得复合材料基本材料性能参数的模拟值；

步骤D，将模拟值带入优化模型的目标函数，判断是否收敛；

步骤E，如果优化没有收敛，则通过ANSYS自带优化算法改变横观各向同性纤维材料的性能参数值，转入步骤C；

步骤F，如果优化收敛，计算结束，则认为此时横观各向同性纤维的材料性能参数值为结果值。

进一步的，所述步骤A中细观力学理论模型为Chamis公式，它的反演公式为：

对于其中缺少的v_f23则根据Christensen公式得到：

其中，E为材料的弹性模量，G为材料的剪切模量，v为材料的泊松比，V_f代表纤维体积分数。下标“f”代表纤维，“m”代表基体。而下标“1，2，3”代表了材料坐标系的主方向，其中“1”代表了纤维长度方向，平面“O-2-3”是各向同性面，如图1所示。

进一步的，所述步骤D中，根据本发明所提出的基于细观力学的有限元优化模型为：

min f(p)

s.t. g(p)＝0

p^l≤p≤p^u

(1)p是设计变量，它可以表示为如下的列向量的形式：

p＝{E_f2,G_f12,v_f12,v_f23}^T

而p^l和p^u则是根据设计变量的初值所确定的设计变量的取值范围。

(2)g(p)是约束方程，它可以表示成如下的形式：

其中，E为材料的弹性模量，G为材料的剪切模量，v为材料的泊松比。上标“c”表示相应的计算值，下标“f”代表纤维，“m”代表基体。而下标“1，2，3”代表了材料坐标系的主方向，其中“1”代表了纤维长度方向，平面“O-2-3”是各向同性面，如图1所示。