精密技工.docVIP

“有限元法及其应用”大作业 一、平面桁架系统静力分析。 已知：桁架结构如图，节点 D 作用载荷 F，a=30cm，b=50cm，各杆横截面均为正方形，横截面积 A=1cm2，材料为 45#钢，外载荷 F=2000N。 要求：给出有限元软件分析的主要步骤、参数设置，加载前后变形图，轴力及轴向应力图；并利用材料力学知识求出解析解，与有限元结果对比。 1、求解过程 主要步骤： 清除当前数据库并开始新的分析 指定工作文件名：这里设置为hengjia 定义标题名：这里设置为hengjiajinglifenxi 定义单元类型：选择LINK180单元 定义实常数：定义单元截面积为0.0001平方米 定义材料属性：选择结构材料的线性且各向同性材料 填写材料相关参数：泊松比设置为0.3，材料弹性模量设置为2.09e11 创建实体模型 网格划分，其中包括材料和单元的属性分配 加载 求解 加载前后的位移图 结果列表 加载前后的轴向应力图 结果列表 理论校核 如右图建立X,Y轴，其中AD杆记为L1，BD杆记为L2，CD杆记为L3，并得节点A的静力平衡方程为 Fx=0, Fn1sina-Fn2sina=0 Fn1=Fn2 (a) Fy=0, Fn3+2Fn1cosa-F=0 根据变形协调方程：L1=L3cosa (b) 综上可得 △L1=(Fn1*b)/(E1*A1*cosa) △L3=(Fn3*b)/(E3*A3) 这两个表示变形与轴力关系的形变可称为物理方程，将其代入（b）得： (Fn1*b)/(E1*A1*cosa)=(Fn3*b/E3*A3) (d) 这是在静力平衡方程外得到的补充方程，从（a）（b）两式得： Fn1=Fn2=(Fcos2a)/(2cos3a+E3*A3/E1*A1) Fn3=F/{1+2*E1*A1*cos3a)/(E3*A3)} 又F=2000N,a=30cm,b=50cm,A=1cm2 代入上式得 Fn1=Fn2=650.389N Fn3=884.721N △L3=2.11e-5m 分析： 理论计算的轴应力Fn1=Fn2=650.389N Fn3=884.721N与有限元结果Fn1’=Fn2’=650.409N Fn3’=884.577N相差不大，可认为有限元结果正确！ 理论计算的杆件位移△L1=1.81e-5m △L3=2.11e-5m与有限元结果△L3’=2.12e-5m几乎一样，可以认为结果正确！ 二、零件模型建立及静力分析 已知： 零件如图所示， 底面矩形长 0.3m， 宽 0.2m， 高 0.04m； 两小孔直径 0.025m， 呈对称分布，两孔中心距为 0.2m；竖直方柱长 0.08m，宽 0.08m，高 0.2m；水 平方柱参数与竖直方柱相同。材料的弹性模量 E=206GPa，泊松比 μ=0.3，密度 ρ=7800kg/m3。底端为固定端。重力加速度 g=9.8m/s2，外力作用于水平方柱右端 面且垂直端面向外，大小为 5000N。 要求： 给出有限元建模和分析的主要步骤、 参数设置， 变形图， 应力及应变云图， 并指出危险点。 求解过程 主要步骤： 清除当前数据库并开始新的分析 指定工作文件名：这里设置为3Dmoxing 定义标题名：这里设置为3Dmoxingjinglifenxi 定义单元类型：选择Solid185单元 定义材料属性：选择结构材料的线性且各向同性材料 填写材料相关参数：泊松比设置为0.3，材料弹性模量设置为2.06e11并设置材料密度为ρ=7800kg/m3 创建实体模型 网格划分，其中包括材料和单元的属性分配 网格划分后的图形 加载 求解 加载前后的变形图 应力云图 应变图 由应力云图可看出竖直长方立柱根部四个角点红色区域应力最大，为危险点！ 三、本课程学习的心得体会、意见及建议 在未学习ANSYS之前我就已经了解到其在工程领域中的应用具有重要作用，所以我对ANSYS的课程学习非常重视！ 一开始在图书馆找了一些自学教材并结合老师上课所讲知识进行了系统学习。当然在学习过程也遇到了不少问题，特别是语言上的问题，还有一些操作并不能完全理解其用意，其实应该是说没有理解其深处的含义，这也让我意识到学习一门专业软件并非一件容易的事，也不是一下就能学会的。所谓冰冻三尺非一日之寒，滴水穿石非一日之功！