含大量随机分布细小纤维的单向纤维增强复合材料横向弹性性能参数计算的内嵌区域模型法研究.pdfVIP

第27卷专辑 固体力学学报 V01．27S．Issue 2006年12月 ACTAMECHANICASOLIDASINICA December2006 含大量随机分布细小纤维的单向纤维增强复合材料 横向弹性性能参数计算的内嵌区域模型法 王晓春h2 冯西桥2 余寿文2 (1复且大学力学与工程科学系，上海，9,00433)(。请华大学航天航空学院工程力学系．北京．100084) 擒要 利用有限元方法求取单向纤维增强复合材料的横向弹性性能参数的计算模型包括三维模型、两维 平面应变模型、单胞模型等等．由于单胞模型仅仅适用于纤维规则排列情况．在纤维随机分布且纤维大小亦为随机 时，单向纤维增强复合材料横向弹性性能参数必须通过对于复合材料块体的计算才能获得．同时在随机分布纤维 的数量增大时，三维模型和二维平面应变模型的计算量急尉增加．模型的处理能力不强．该文提出一种利用内嵌区 域模型来计算古大量随机大小、随机分布细小纤维的单向纤维增强复合材料块体的横向弹性性能参数的方法，有 效降低了计算量．在较低的计算费用下，能够快速获得单向纤维增强复合材料的横向弹性性能参数． 关■词 单向纤维增强复合材料，有限单元法，内嵌区域模型，横向弹性性能参数，F．shelby张量。应力(变)集 中张量 0 引育 维模型时，即使计算仅仅只含一根纤维的单胞模型， 单向纤维增强复合材料是有着重要商业和技术 其计算量也是很大的．对于长纤维增强复合材料来 应用的复合材料．由于层板复合材料可以通过单向 说，二维模型为平面应变模型．与三维模型相比较， 纤维增强复合材料在不同方向上堆叠粘结而成，而 二维模型大大降低了计算量，或者说有效地增强了 层板复合材料又是目前应用最为广泛的复合材料之 模型处理相关问题的能力．不过即使如此，当单向纤 一．因此研究单向纤维增强复合材料的力学性能十 维增强复合材料基体中含有大量随机分布的细小纤 分重要．目前计算单向纤维增强复合材料弹性性能 维时，计算单向纤维增强复合材料块体弹性性能参 参数的方法很多，主要分为三类：解析法，经验法和 数时，计算量仍然很大．二维平面应变模型所能处理 数值法[1]．如果纤维在单向纤维增强复合材料中随 的块体中含有的纤维数依然是相当有限的．近来纳 机分布，纤维增强复合材料可被视作横向各向同性 米技术获得了飞速发展和广泛的应用，致使复合材 材料，具有五个独立的弹性常数：轴向弹性模量F。， 料的细观尺度逐渐向纳米量级过渡．这使得制备含 轴向剪切模量G，。，横向弹性模量Ec，横向剪切模量 有大量随机分布的微纳米尺度增强相复合材料成为 G，以及离面泊松比y7。．大量文献和实验证实单向 可能．因此，寻找一种简单而且快速的方法，在最经 长纤维增强复合材料的轴向弹性性能参数E，。，G，。 济的条件下通过有限元计算预测含大量随机大小、 和离面泊松比I，7。，与复合材料的弹性性能参数理论 随机分布细小单向纤维增强复合材料的弹性性能参 值Voigt上限相当接近[2]．但是横向弹性性能参数 数的方法，是至为重要且有相当实用意义的． 并不能从Reuss下限[3】得到，它们之间差异甚大r4j．按照Eshelby理论，嵌于无限大均匀介质中的 一些半经验公式常被用来计算纤维增强复合材料的 有限小区域，如果发生特征应变，那么其实际所发生 横向弹性性能参数，如Halpin-Tsai方程口]等．随着 的应变应为Eshelby张量与特征应变的双点 计算机技术的发展．利用有限单元法计算复合材料 积Ⅱ川o]．利用Eshelby张量，也可以求得嵌于无限大 性质的独特优势慢慢显示出来．如果纤维规则排列 基体中的小区域在承受远场均匀应力场时的应力集 在基体中，通常选择包含一根纤维的单元的