基于Ansys的颗粒增强复合材料力学行为分析.doc

基于Ansys的颗粒增强复合材料力学行为分析 陈 晨 力学与工程学院，结构2010-01班【摘要】本文中，将以T6热处理的SiCp/6061A1合金复合材料为研究对象，在细观分析模型的基础上，建立一定的体积代表性单元，将颗粒和基体作为两种材料建立模型，通过施加适当的边界条件，得出复合材料的应力场和应变场，然后进行体积平均，得到等效应力和等效应变，从而对该复合材料的单轴循环应变行为进行有限元分析。 【关键词】固体力学，复合材料力学，Ansys有限元分析 1. 问题描述 nsys软件将颗粒和基体作为两种材料建立模型，通过施加适当的边界条件得出复合材料的应力场和应变场，进行体积平均得到等效应力和等效应变，从而计算出复合材料的力学性能。 参数名 C1 C2 C3 C4 C5 C6 C7 C8 C9 EX NUXY Chaboche 206 8e4 1000 3e5 2e5 1600 400 0 0 7e4 0.3 在有限元模型中，各尺寸满足且。边界条件为对称边界条件。分析中， 图2 三维代表性体积单元 图3 加载方式（应力加载和循环应变加载） 2. 建立模型 2.1选择分析模式2.2 定义单元类型 图选择单元类型 图选择2.3定义 图 2.4定义模型参数 为了后续工作的方便，在建模之前，首先定义了建模所需参数。利用Utility Men|Parameters|Scalar Parameters进行定义。其各项参数如下表：  表2 模型参数 参数名称 代码 值 颗粒体积分数 F 0.14 基体长度 L 20 球颗粒半径 R 应变幅值 STR 0.6 图  2.4建立几何模型 图6 建立实体模型 2.5 网格划分 选择MeshTool网格划分工具，首先对线进行划分。其中，在接触处的L5进行加密，以便使结果更精确，如图7。 图7对线划网格 然后对基体和颗粒分别选材进行划分，划分方式为Sweep，划分后的网格如图8所示。 图8 划分好的网格 2.6 定义边界条件 建立有限元模型后，就需要定义边界条件并施加荷载，才能进行求解。本实例的位移边界条件是将模型左侧底面边界上节点的全部位移都约束住。具体步骤为：选择Main Menu|Preprocessor|Loads|Define Loads|Apply|Structural，定义点选Displacement|On Lines，弹出面选择对话框，要求选择欲对其施加位移约束的线。完成后的约束如图9。 图9 施加边界约束 3.求解 3.1求解 本节中将利用施加的循环压力载荷进行求解。单击 Main Menu |Solution |Analysis Type |Sol’n Control进行求解设置。在上边界施加位移，共五次，每施加一次，求解一次，写一个载荷步文件，进行循环求解。 图10 在顶端施加位移 3.2查看结果 下面将读入结果数据并在中查看模型的von Mises应力分布。通过点选Read Result，分别查看各个载荷状态下的结果。以下列出了应变控制下，体积分数为14%的复合材料在第1个循环结束的等效von Mises应力分布和等效塑性应变分布；第10个循环结束的等效von Mises应力分布和等效塑性应变分布。 图11 第1个循环的等效应力分布图 图13 第1个循环的等效塑性应变分布图 图14 第10个循环的等效应力分布图 图15 第10个循环的等效塑性应变分布图 3.3结果分析 点击Timehist Postproc，进入时间历程处理器，点选最大塑性应变节点155，分别添加其y向弹性、塑性应变，相加得到总应变，绘制其应力应变曲线如图17。 图16 时间历程处理器 图17 体积分数为14%时最大塑性应变处的应力应变曲线 3.4计算平均应力应变 为了绘制与复合材料应力-应变（有效性能）相关的曲线，还需求出复合材料的平均应力和平均应变。在这里，采用体积平均的方法求复合材料在各个时刻的平均应力和应变。相应计算公式为： 再次进入时间历程处理器，添加应力最小处的节点1，按公式计算出其平均应力与应变，绘制曲线如图20。 图19 添加节点 图20 体积分数为14%时单轴整体应力-应变循环曲线 4.扩展 4.1求解 将体积分数F改为21，即体积分数为21%重新建立模型进行计算，按照前述方法，再次获得应力应变曲线如下图： 图21 体积分数为21% 时最大塑性应变处的应力应变曲线 图20 体积分数为21% 时单轴整体应力-应变循环曲线 4.结语 由以上分析可知，由于颗粒增强体的引入，使得基体抵抗塑性变形的能力增强，体积分数越大，产生的塑性应变越小；最大塑性产