复合细观有限元专题讲述.pptVIP

专题 复合材料细观有限元分析 §A-1 短纤维复合材料力学行为分析 二、循环变形行为的模拟 采用如图所示的代表性体积单元来分析复合材料的整体行为。 在有限元分析中可采用2-D轴对称8节点单元进行分析（如图中所示），边界条件为对称位移边界条件。 各个尺寸之间满足如下关系： L/2=l/2+2.5d，ld2 /LD2=Vf Al2O3纤维为弹性材料，弹性模量Ef=300 GPa，泊松比?f=0.20。 Al合金基体为弹塑性材料，其循环本构模型描述如下。 §A-2 颗粒增强复合材料力学行为分析 * 一、有效模量分析 1、短纤维复合材料的形貌 面内分布 厚度方向分布 单向短纤维复合材料的理想化模型 2、有限元分析模型 层板比拟法 p为离散的单元号，n为单元总数。 只需求出了 ，即可得： 三维代表性体积单元 所有的计算都是基于上述代表性体积单元。对随机分布短纤维复合材料的处理方法与前一致。 斜向纤维等效模型 不同的方法得到的结果不同，见下表。 -- 78 83 88 -- 78 83 88 -- 78 88 夹杂理论 -- 76 80 84 -- 76 80 84 -- 76 84 H-T方程 -- 85 93 102 -- 85 93 102 -- 85 102 混合律 70 78.1~80.2 85.2~89.8 94.2~97.2 70 78.9 87.4~89.2 94.8~95.6 70 73.6~75.0 80.6 -- 81.4 87.7 93.9 -- 81.4 87.7 93.9 -- 81.4 93.9 0 10 15 20 0 10 15 20 0 10 20 ?-Al2O3f/Al-5.5Mg ?-Al2O3f/Al-5.5Zn ?-Al2O3f/Al-12Si 测量 FEM Vf 复合材料 （单位：GPa） 二、有效性能分析（拉伸应力-应变曲线） 基体材料的弹塑性力学性能参数 不同纤维位向角下的拉伸应力-应变曲线 复合材料代表性体积单元以及有限元分析网格（Vf=15%，l/d=20） 主控方程： 屈服函数： 背应力演化方程： (k=1, 2, …, M) 为材料常数 该模型能够合理描述未增强基体的循环变形行为。 基体材料的材料参数值 ξ(1)=5000, ξ(2)=1250, ξ(3)=800, ξ(4)=400, ξ(5)=200, ξ(6)=100, ξ(7)=50, ξ(8)=33.3, ξ(9)=25, ξ(10)=20; r(1)=4.6, r(2)=3.5, r(3)=4.5, r(4)=2.4, r(5)=5.2, r(6)=5.4, r(7)=2.6, r(8)=2.1, r(9)=2.0, r(10)=12.5 (MPa). M=10, E=70GPa, v=0.33, Q0=15MPa; 上述模型中的材料参数可以通过单拉实验确定，对纯Al的值如表1所示。 未增强基体的?r~N模拟结果 复合材料的?r~N模拟结果 1、颗粒增强复合材料的形貌 Vp=14% Vp=21% 2、颗粒增强复合材料的力学行为 单调拉伸 应变循环 棘轮行为 3、单颗粒有限元模型 轴对称边界条件 Kang GZ, et al: Mater. Sci. Eng. A 426(2006), 66; Mater. Sci. Eng. A 458(2007), 170. 单拉应力-应变曲线 单轴应变循环曲线 棘轮行为 (c) body-centered cube (BC) single particle (SP) (b) simple cube (SC) (d) face-centered cube (FC) (e) edge-centered cube (EC) (f) body-face-edge-centered cube (BFE) 4、多颗粒有限元模型 晶体点阵结构 单胞及有限元网格（规则排列）