平纹机织碳纤维复合材料随机力学性能的预测.docxVIP

平纹机织碳纤维复合材料随机力学性能的预测 面向平等的不确定性预测 碳纤维增强聚合物基（cfrp）是一种碳纤维材料。密度低、硬度高、强度高、设计能力强，密度低、比刚强、设计能力强等优点。结构的量化和高性能设计可以得到很好的应用。 为了实现对考虑随机性的复合材料力学性能的准确预测,国内外学者开展了广泛的研究工作.与有限元技术相结合的计算细观力学方法可以在各个尺度上建立考虑随机性的精细仿真模型来进行力学性能预测 基于计算细观力学的方法本质上是建立多个能够满足统计学表征的有限元模型并统计仿真分析的结果,因此计算代价过高.而解析细观力学方法由于具有较高的计算效率,在进行考虑不确定性的预测研究时具有很大的优势.Chakraborty等 现有的面向平纹机织CFRP的研究主要存在以下几个问题:一是所考虑的随机参数不够全面,大多研究仅考虑部分尺度的随机参数.且采用的不确定性传递方法大多具有局限性:随机有限元法,作为一种嵌入式的不确定性算法,对问题不具备很好的适应性;摄动法则对概率分布形式具有要求;二是在力学性能预测中对随机参数之间的相关性考虑不够充分.针对以上这些问题,本文提出一种基于PCE和Vine Copula的平纹机织CFRP随机力学性能预测方法,该方法考虑了微观和介观尺度中不确定性材料、结构参数,基于自下而上层级传递的策略逐层研究平纹机织CFRP的各尺度力学性能参数的不确定性.所采用的非嵌入式PCE方法能够高效准确地实现不确定性地传递.此外,Vine Copula函数能准确构造相关随机参数之间的高维联合概率分布,从而有效掌握相关性对性能响应的影响,提高平纹机织CFRP力学性能参数预测的精度和可靠性. 1 尺度确定性随机弹性力学性能预测 采用解析细观力学方法实现平纹机织CFRP的多尺度确定性随机弹性力学性能预测,该方法能充分考虑影响平纹机织CFRP力学性能的因素,包括组分材料的力学性能、各尺度的几何参数、体积分数-等 1.1 纤维束刚度矩阵 平纹机织CFRP中,每一根浸润过的纤维束都包含了嵌入在基体中的许多单向纤维,浸润纤维束被认为是单向复合材料.若纤维束的纤维丝体积分数为V 式中,下标f和m分别表示纤维丝和基体,V表示体积分数,E表示弹性模量,G表示剪切模量,u表示泊松比. 纤维丝的体积分数可以通过下式计算 其中,S 由于纤维束通常被假设为横向各向同性的,故其刚度矩阵可以定义为 其中,刚度矩阵各元素可通过下式获得 由于平纹机织CFRP的纤维束存在卷曲变形,对于卷曲的纤维束,如图3所示. 本文采用正弦函数来描述其构型 式中,T表示周期,A表示幅度.综合图2和图3可以得到T=2(a+l),A=(h-b)/2.将卷曲纤维束划分为无数的微段,各微段沿纤维主方向均可视作单向纤维复合材料,如图4所示,x 定义第N微段在其局部坐标系下的刚度矩阵为C 式中,上标L表示微段纤维局部坐标系(x 式中 式中 式中θ 1.2 刚度矩阵转换 对于CFRP,在宏观尺度上选取单胞来表征其宏观性能.由于平面机织CFRP是由经纱和纬纱机织而成,单胞整体的刚度矩阵可根据各纤维束的平均性能计算得到,由于经纱和纬纱的方向不同,需通过坐标转换,将纤维束刚度矩阵转化为单胞坐标系下的纤维束刚度矩阵 转换矩阵H 式中 式中,θ 式中,V 通过对宏观刚度矩阵进行求解,可获得单胞的弹性模量、剪切模量和泊松比. 1.3 cfrp的力学性能 本文所制备的平纹机织CFRP碳纤维丝和基体材料力学性能参数如表1所示,碳纤维丝选用台丽碳纤维TC33,碳纤维以3K的形式集合成碳纤维束.基体选择Huntsman公司的环氧树脂,牌号为LY1564SP,对应的固化剂牌号为Aradur 3486,两者的体积比为5:2.采用真空导入工艺进行力学性能试验所用的平纹机织CFRP的成型.通过上述工艺制备得到的CFRP的密度为1.47 g/cm 为了从宏观上探究CFRP的力学性能,分别根据试验标准ASTM D638 纤维丝和纤维束的几何参数通过对纤维束界面显微图像和Micro-CT图像(如图6和图7所示)进行数据统计可获得,纤维丝的直径为d=6.23μm,纤维束的几何参数如表2所示.表中平均值为统计得到的均值,上下限值为统计得到的±3σ置信区域数值.纤维丝占纤维束体积分数的均值为65.32%,纤维束占单胞体积分数的均值为67.90%. 以各变量的均值为输入参数,经过前述的解析模型预测得到的单向纤维和单胞的确定性力学性能如表3所示,通过与有限元仿真结果的对比可看出解析法的相对误差较小.特别地,解析法得到的结果与试验结果的相对误差分别为0.82%和6.33%,预测结果与试验结果有较好的一致性. 2 cfrp的层级结构 本文主要考虑组分材料力学性能和几何参数不确定性对单向纤维力学性能预测结果的影响,组分材料的力学性能参