为了适应公司新战略的发展，保障停车场安保新项目的正常、顺利开展，特制定安保从业人员的业务技能及个人素质的培训计划复合材料杨氏模量计算 三维取向短纤维增强复合材料弹性模量的计算 赵延军1刘春太2李克华1史冬丽1 摘要：基于单向取向纤维增强复合材料的力学性能计算模型，借助于纤维取向分布函数及坐标转换，建立了三 维取向短纤维增强复合材料弹性模量的数值计算模型按该模型对短纤维增强树脂基复合材料的弹性模量进行计算，并将其结果与同类材料的实验结果比较验证结果表明，该模型的预测具有较好的准确性 关键词：纤维取向，弹性模量，短纤维增强复合材料短纤维增强树脂基复合材料由于具有较高的强度、弹性模量、刚度，以及抗蠕变性能好等优点，近年来得到了广泛的应用该类制品的成型过程中，由于模腔几何形状的不同，成型条件特别是填充速度的差异，纤维及基体材料的物理性能不同，使成型后的制品不同部位呈现不同的纤维取向一般来说，短纤维在制品中的取向是三维的而纤维增强复合材料的弹性模量等力学性能及热物理性能依赖于基体和纤维的本身性能、纤维含量、纤维的取向状态和纤维与基体间界面结合因此，研究短纤维增强复合材料的基体和纤维的性能、纤维含量、纤维的取向状态与弹性模量等力学性能之间的关系，对制品及模具设计、提高纤维增强复合材料制品的质量等具有重要的指导意义。
v12?G23?Gm? 界条件，Tendon-Weng模型推出单向取向纤维增强复合材料的弹性模量、剪切模量及泊松比分别为：E11?E22? Em 1?c(A1?2?mA2)/A Em 1?c[?2?mA3?(1?vm)A4?(1?vm)A5A]/2A G12?Gm? G0c Gm/(Gf?Gm)?2(1?c)E1212 Gmc Gm/(Gf?Gm)?2(1?c)E2323 vmA?c(A3?vmA4)A?c(A1?2vmA2) 其中，E、G、v分别表示弹性模量、剪切模量和泊松比，下标m和f分别代表聚合物基体和增强纤维，c为复合材料中纤维所占的体积百分含量，E1212、E2323为Eshelby张量的分量，A、A1、A2、A3、A4、A5分别为纤维体积百分含量、长径比及Eshelby张量的显式函数[1] 2．纤维取向分布函数 将短纤维增强的塑料熔体近似为浓悬浮液，假设纤维为圆柱形刚体，长度和直径均匀，则可用单位矢量P等效描述纤维取向矢量P的方向与纤维轴线重合而空间一点的纤维取向状态可用取向分布函数Ψ(P)描述，Ψ(P)定义为一根纤维取向于向量P方向的概率。
1.单向取向纤维增强复合材料力学性能的计算 许多研究者对单向取向纤维增强复合材料的力学性能预测方面已经做了大量的研究工作，并得到了很多在理论和工程上都比较有用的理论或半理论半经验模型其中包括自恰模型、Mori-Tanaka模型、Tendon-Weng模型及Halpin-tsai模型等在上述各种模型中，通过对比研究发现，Tendon-Weng充分考虑了纤维长径比、纤维体积百分含量、及复合材料各组分的物性参数对复合材料力学性能的影响及纤维之间的相互作用，无论在纤维体积百分含量较低或较高时都有较好的准确性因此，本文将以此模型为单向取向纤维增强复合材料力学性能的计算模型 在纤维在复合材料中呈单向取向、基体和纤维呈现横观各向同性的前提下，并结合一定的边 图1纤维空间取向角的定义 1 基金项目：科研院所技术开发研究专项资金项目 其中， ?sin?sin?(6)?sin?cos?,PP?cos?,P321 这样，就把三维取向的复合材料简化为具有不同θ角的材料层组成的层片组该转换过程可用图2直观的表示出来图2为具有空间随机取向的实际短纤维复合材料，从图中可以看到，实体各面都有纤维末端穿过，说明实际复合材料的纤维取向是随机的。
如要计算1方向的弹性模量，则以1方向为主方向，如图中坐标系所示θ表示纤维取向方向与轴1方向的夹角由于只需 考虑θ角，所以可用Ψ代替Ψ，实际随机纤维取向的复合材料简化为具有不同取向角φ的不同材料层 在上述假设下，单根纤维在空间某一点的取向行为可以用图1中的角度定义，空间中一点的取向状态可以采用概率分布函数Ψ(θ,φ)描述概率分布函数Ψ(θ,φ)的定义为：在角θ1和θ1+dθ，φ1和φ1+dφ区域内发现一个纤维的可能性，即， P(?1????1?d?,?1????1?d?) ??(?,?1)sin?1d?d? 1 由函数Ψ的定义可知，它必须满足一定的物理条件首先，纤维在角度θ，φ和(π-θ,φ-π)上取向无区别，即满足周期性条件： ?用Tendon-Weng模型求得各单向取向层材料的弹性模量； 2）由1 ）中所得的各层弹性模量，结合弹性 模量与刚度系数的关系求出各层的刚度系数Cij即，C11?E11(1?v12v21)(11)C12?v21C11(12)C12?v21C11(13) C22?E22(1?v12v21)(14)C26?0(15) C66?G12(16) 对模型进行描述，但在实际纤维取向的数值计算中，则采用了一个二阶取向张量[3]来代替取向分布函数。
即 aij?pipj?(P)dP(23) 由于该纤维取向张量是一个表征纤维取向状态的二阶实对称张量，故可用Yacabi方法计算出对应的特征向量及特征值，即： ??10aij???0?2 ??000? 0??:?e1e2?3?? e3?(24) 3)由转角公式将偏轴刚度系数转化为正轴刚 度系数即'?C11??m4n42m2n24m2n2??'??? n4m42m2n24m2n2??C11??C22?? ?'?22????m2n2m4?n4?4m2n2??C12??mn?C22? ?'??????22 m2n2?2m2n2(m2?n2)2??C12??C66??mn '???C16m3n?mn3mn3?m3n2(mn3?m3n)???C66???'??? ?mnmn?mn2(mn?mn)???mn??C26??? (17) 其中，m=cosθ，n=sinθ 4）将各层的刚度系数进行积分得到实际复合材料的刚度系数：Cij? 该特征向量的在物理上表示纤维可能的取向，与之对应的特征值表示纤维在特征向量方向取向的可能性分别由各特征向量求出其取向角θ，代入偏轴刚度系数与正轴刚度系数之间的转换矩阵并用离散的?值代替连续的纤维分布函数Ψ(θ)即可算出复合材料的弹性模量值。
于是，式18可变为： Cij? ??? ? ?max ' Cij??(25) min ? ?max ?min ' Cijg(?)d?(18) 其中，λθ表示取向角为θ的短纤维所对应特征向 量对应的特征值 利用上述三维纤维取向复合材料弹性模量计算模型，对具有完全随机三维取向的短纤维增强聚合物基复合材料的弹性模量进行了计算算例中所用纤维和基体的力学及其它性能参数见表1该算例分别计算了纤维质量百分比为30％和40％两种含量的E-玻璃纤维增强PPS的力学性能由于选定的复合材料的代表体积单元的纤维呈三维完全随机取向，所以，可认为该材料是各向同性的表2分别将这两种不同纤维含量的复合材料弹性模量的理论计算结果与Cheng[4]的实验结果进行了对比，并计算了理论值与实验值的误差由表中的数据对比可以看出，该模型预测值与实验值的最大误差率为％，平均误差为％，说明该模型的计算结果能满足计算要求 表1算例中所用纤维和基体材料的性能参数材料PPSE-glassfiberEG 476 v 平均长径比a /25 其中，0o≤θmin≤θ≤θmax≤90o 5)最后再由刚度系数与弹性模量的关系求出实际具有三维取向的短纤维增强复合材料的弹性模量： E11?(C11C22?C12)/C22(19) E22?(C11C22?C12)/C11(20)G12?C66(21) 22 v12?C12/C22(22)5.计算结果与分析 纤维取向分布函数具有很明显的物理意义，并能完全、清晰地描述纤维取向分布状态，但其不利之处是非常麻烦，利用数值分析几乎不可能完成这些计算任务。
故本文用纤维取向分布函数 表2弹性模量实验值与理论值的对比 质量百分比3040 体积百分比 E(实验)E(计算)E 聚合物基复合材料弹性模量的简单数值计算模型利用该模型对具有完全随机三维取向的短纤维增强PPS的弹性模量进行了计算，并将其结果与同类材料的实验结果进行了比较结果表明，该模型的预测具有较好的准确性 参考文献 1.G.P.Tandon,G.J.Weng.Propertiesof UnidirectionallyAlignedComposites[J].PloymerComposites,1984,5(4):327-333.2.Shao-YunFu&BerndLauke.Theelastic modulusofmisalignedshort-fiber-reinforcedpolymers[J].CompSciTechnol,1998,58:389-400. 3.AdvaniS.G,TuckerC.L..TheUseofTensors toDescribeandPredictFiberOrientationinShortFiberComposites[J],JournalofRheology,1987,31:98-101. 4.Chao-hsunChen,Chiang-HoCheng．Effective elasticmoduliofmisorientedShort-Fiber Composites[J],SolidsStructures,1996,33(17):2529-2539. 图3给出了利用该模型计算所得三维完全随机纤维取向玻璃纤维增强PPS的弹性模量E随纤维体积百分数的变化曲线。
由图可见，在该复合材料中，弹性模量E随纤维体积百分含量的增加而增大其变化趋势与实际复合材料的弹性模量随纤维体积含量增加时的变化规律一致 图3弹性模量E随纤维体积百分数的变化曲线6.结论 基于单向取向纤维增强复合材料力学性能计算的Tendon-Weng模型，借助于纤维取向分布函数及坐标转换，建立了具有三维取向短纤维增强 ElasticModulusPredictionofThreeDimensionaloriented Short-Fiber-ReinforcedComposites ZHAOYan-jun1LIUChun-tai2LIKe-hua1SHIDong-li1 (1.ZhengzhouResearchInstituteforAbrasives&Grinding,Zhengzhou,) (2.NERCofMold?Die,ZhengzhouUniversity,Zhengzhou,) Abstract:Anelasticmoduluspredictionmodelof3Dorientedshort-fiber-reinforcedcompositesisconcludedbasisonthemechanicalpropertypredictionmodelofunidirectionalorientedfiberreinforcedcomposites,fiberorie。