一种基于细观力学模型的横观各向同性纤维材料性能的预测方法
【技术领域】
[0001] 本发明设及横观各向同性纤维材料性能的预测，适用于复合材料细观力学方面的 研究，具体设及一种基于细观力学模型的横观各向同性纤维材料性能的预测方法。
【背景技术】
[0002] 复合材料由于具有比强度、比刚度高及性能可设计等优点而越来越广泛地应用于 航空航天飞行器结构中，其应用范围已经从最初的非承力结构、次承力结构扩展到主承力 结构。研究表明，在设计与分析复合材料力学结构时，对于复合材料自身的力学性能的研究 是十分重要的。
[0003] 在研究复合材料力学性能的方法中，除了通过试验的手段来测得外，还有学者提 出了基于细观力学的方法来预测复合材料的材料性能，该其中又包含了理论公式法W及基 于代表性体积单元的有限元方法。然而，通过该两种基于细观力学的方法预测复合材料的 材料性能，均需要已知复合材料的细观材料性能，对于纤维增强复合材料而言，即纤维与基 体的材料性能。然而，由于纤维的尺寸过小，难W通过试验的方法获得纤维的材料性能，特 别当纤维不是各向同性材料时。
[0004] 随着复合材料越来越多地应用于飞机结构，基于细观力学的损伤及失效分析受到 了越来越多的注视。对于纤维增强复合材料而言，在进行基于细观力学的损伤及失效分析 时，需要考虑组分各自的退化，因此纤维与基体的材料性能是必须的。另外，复合材料的细 观材料性能对于材料的热弹性应力分析也是必要的。
【发明内容】
[0005] 本发明要解决的技术问题为：克服现有技术的不足，W复合材料基本材料性能、基 体的材料性能W及纤维长度方向的弹性模量为输入参数，通过细观力学理论模型反演，得 到横观各向同性纤维的材料性能的初值，再利用基于代表性体积单元的有限元优化方法计 算得到横观各向同性纤维的材料性能，给出了一种求解横观各向同性纤维材料性能的优化 方法。
[0006] 本发明解决上述技术问题采用的技术方案为；一种基于细观力学模型的横观各向 同性纤维材料性能的预测方法，实现步骤如下：
[0007] 步骤A，根据复合材料的基本材料性能参数、基体的材料性能W及纤维长度方向的 弹性模量，利用细观力学理论模型，反演得到横观各向同性纤维材料性能参数的初值；
[000引步骤B，建立基于细观力学的代表性体积单元有限元模型，带入纤维材料性能参数 的初值W及基体的材料性能，施加周期性边界条件及四种载荷；
[0009] 步骤C，基于步骤B中的有限元模型，预测复合材料的材料属性，获得复合材料基 本材料性能参数的模拟值；
[0010] 步骤D，将模拟值带入优化模型的目标函数，判断是否收敛；
[0011] 步骤E，如果优化没有收敛，则通过ANSYS自带优化算法改变横观各向同性纤维材 料的性能参数值，转入步骤C;
[0012] 步骤F，如果优化收敛，计算结束，则认为此时横观各向同性纤维的材料性能参数 值为结果值。
[0013] 进一步的，所述步骤A中细观力学理论模型为化amis公式，它的反演公式为：
【主权项】
1. 一种基于细观力学模型的横观各向同性纤维材料性能的预测方法，其特征在于包括 以下步骤： 步骤A，根据复合材料的基本材料性能参数、基体的材料性能以及纤维长度方向的弹性 模量，利用细观力学理论模型，反演得到横观各向同性纤维材料性能参数的初值； 步骤B，建立基于细观力学的代表性体积单元有限元模型，带入纤维材料性能参数的初 值以及基体的材料性能，施加周期性边界条件及四种载荷； 步骤C，基于步骤B中的有限元模型，预测复合材料的材料属性，获得复合材料基本材 料性能参数的模拟值； 步骤D，将模拟值带入优化模型的目标函数，判断是否收敛； 步骤E，如果优化没有收敛，则通过ANSYS自带优化算法改变横观各向同性纤维材料的 性能参数值，转入步骤C ; 步骤F，如果优化收敛，计算结束，则认为此时横观各向同性纤维的材料性能参数值为 结果值。
2. 根据权利要求1所述的一种基于细观力学模型的横观各向同性纤维材料性能的预 测方法，其特征在于所述步骤D中的优化模型为： min f (p) s. t. g(p) =O ρ1^： ρ p u (1) P是设计变量，它可以表示为如下的列向量的形式： P 一 {Ef2, Gfl2, vfl2, vf23} 而P1和P 是根据设计变量的初值所确定的设计变量的取值范围； (2) g(p)是约束方程，它可以表示成如下的形式：
其中，E为材料的弹性模量，G为材料的剪切模量，V为材料的泊松比，上标"C"表示相 应的计算值，下标"f"代表纤维，"m"代表基体。而下标" 1，2, 3"代表了材料坐标系的主方 向，其中" 1"代表了纤维长度方向，平面" 0-2-3 "是各向同性面； 另外，约束方程中的4个h函数则表示代表性体积单元在4种边界条件及载荷状况下 的作用效果； 输入参数除了设计变量P以及值为已知量的Efl A111和V111外，Gf23由以下的关系式确定：
(3) f(p)是目标函数，它的表达式如下：
其中，E为材料的弹性模量，G为材料的剪切模量，v为材料的泊松比，上标"c"表示相 应的计算值，"e"表示相应的试验值。
【专利摘要】本发明涉及一种基于细观力学模型的横观各向同性纤维材料性能的预测方法，以复合材料基本材料性能、基体的材料性能以及纤维长度方向的弹性模量为输入参数，通过细观力学理论模型反演，得到横观各向同性纤维的材料性能的初值，再利用基于细观力学模型的有限元优化方法计算得到横观各向同性纤维的材料性能。本发明给出了一种求解横观各向同性纤维材料性能的优化方法，利用基于细观力学的有限元模型得到了难以试验测得的横观各向同性纤维的材料性能，为复合材料的细观力学研究提供了必要的材料性能参数。
【IPC分类】G06F19-00
【公开号】CN104537259
【申请号】CN201510015046
【发明人】赵丽滨, 李杨, 张建宇, 刘丰睿, 陈玉丽
【申请人】北京航空航天大学
【公开日】2015年4月22日
【申请日】2015年1月12日