有限元实例分析大作业一

合集下载

ansys有限元分析工程实例大作业

ansys有限元分析工程实例大作业————————————————————————————————作者：————————————————————————————————日期：辽宁工程技术大学有限元软件工程实例分析题目基于ANSYS钢桁架桥的静力分析专业班级建工研16-1班（结构工程）学号 471620445姓名日期 2017年4月15日基于ANSYS钢桁架桥的静力分析摘要：本文采用ANSYS分析程序，对下承式钢桁架桥进行了有限元建模;对桁架桥进行了静力分析，作出了桁架桥在静载下的结构变形图、位移云图、以及各个节点处的结构内力图（轴力图、弯矩图、剪切力图），找出了结构的危险截面。

关键词：ANSYS；钢桁架桥；静力分析；结构分析。

引言：随着现代交通运输的快速发展，桥梁兴建的规模在不断的扩大，尤其是现代铁路行业的快速发展更加促进了铁路桥梁的建设，一些新建的高速铁路桥梁可以达到四线甚至是六线，由于桥面和桥身的材料不同导致其受力情况变得复杂，这就需要桥梁需要有足够的承载力，足够的竖向侧向和扭转刚度，同时还应具有良好的稳定性以及较高的减震降噪性，因此对其应用计算机和求解软件快速进行力学分析了解其受力特性具有重要的意义。

1、工程简介某一下承式简支钢桁架桥由型钢组成，顶梁及侧梁，桥身弦杆，底梁分别采用3种不同型号的型钢，结构参数见表1，材料属性见表2。

桥长32米，桥高5.5米，桥身由8段桁架组成，每个节段4米。

该桥梁可以通行卡车，若只考虑卡车位于桥梁中间位置，假设卡车的质量为4000kg，若取一半的模型，可以将卡车对桥梁的作用力简化为P1，P2，和P3，其中P1=P3=5000N，P2=10000N,见图2,钢桥的形式见图1，其结构简图见图3。

图1钢桥的形式图2桥梁的简化平面模型（取桥梁的一半）图3刚桁架桥简图所用的桁架杆件有三种规格，见表1表1 钢桁架杆件规格杆件截面号形状规格顶梁及侧梁 1 工字形400X400X12X12桥身弦杆 2 工字形400X300X12X12底梁 3 工字形400X400X16X16所用的材料属性，见表2表2 材料属性参数钢材弹性模量EX泊松比PRXY 0.3密度DENS 78002 模型构建将下承式钢桁梁桥的各部分杆件，包括顶梁及侧梁，桥身弦杆，底梁均采用BEAM188单元，此空间梁单元可以考虑所模拟杆件的轴向变形; 定义了一套材料属性，各类杆件为钢材，其对应的参数如表2所示；根据表1中的杆件规格定义了三种梁单元截面，根据表1分别定义在相应的梁上；建模时直接建立节点和单元，在后续按照先建节点再建杆的次序一次建模。

有限元的大作业报告示例

1.题目概况
矩形板尺寸如下图1，板厚为5mm。

材料弹性模量为
松比μ= 0.27 。

施加约束和载荷并讨论：
图
1 计算简图
1.1基本数据
E = 2⨯105N/mm2，泊
序号载荷约束备注42 向下集中载荷F=800N, 作用于cd 边3/4 处（近d） c d 点简支
1.2分析任务/分析工况
讨论板上开孔、切槽等对于应力分布的影响。

（载荷约束组合不变）。

提示：各种圆孔，椭圆孔随大小、形状、数量，分布位置变化引起的应力分布变化；各种形状，大小的切槽及不同位置引起应力分布的变化等，选择二至三种情况讨论，并思考其与机械零部件的构型的相对应关系。

2.模型建立
2.1单元选择及其分析
由于平板长宽分别为300x100，故可取网格单元大小为1。

如图：
2.2模型建立及网格划分
模型按单元为1 划分后的网格大小如图所示：
2.3载荷处理
向下集中载荷F=800N, 作用于cd 边3/4 处（近d） c d 点简支
3.计算分析
3.1位移分布及其分析
（1）位移分布如图：。

有限元静力学分析实例

图26定义载荷随时间变化的幅值
（2）定义集中力在主菜单中选择Load——Manager，在弹出的对话框中进行如图27所示的设置，然后点击Continue，在弹出的对话框中进行如图28所示的设置，然后点击OK。
图27图28
在Load Manager对话框中可以看到，如图29所示，名为Load-Point的载荷在分析步Step-Load-1中开始起作用，并延续（propagate）到分析步Step-Load-2中。
图30图31
图32
（4）定义受剪力的面返回Load功能模块，在主菜单中选择Tools——Surface——Create，在弹出的对话框中，在Name后面输入Surf-Load，然后点击Continue，选择如图33红色区域。
图33
（5）定义面载荷在Load Manager对话框中再次点击Create，进行如图34所示的设置，在弹出的Edit对话框中，设置向量起始坐标（0.0,0.0,0.0），回车确认，输入向量的终点坐标（0.0,-10.0,0.0），再次回车确认。设置Magnitude（大小）为36，点击OK，得到的部件如图35所示。
2.由于关心的是应力集中部位的应力状态，所以在模型中使用C3D20R单元（20节点六面体二次减缩积分实体单元）。
3.基于结构和载荷的对称性，可以只取模型的1/2进行分析（见图1）。
4.由于圆孔处螺纹的应力应变状态不是所关心的重点，可以简化杆件和圆孔之间的连接关系（如图1所示），不对杆件和螺纹精确建模，而是在杆件一端的受力点和圆孔内表面之间建立分布耦合约束。
图8
五、定义装配体
在Module列表中选择Assembly功能模块。点击左侧工具区中的（Instance Part），接受默认参数，如图9所示，即类型为非独立实体，点击OK。

哈工大有限元大作业

作业一一．计算程序和结果展示1.程序clearsyms a b c x E lD=E*pi*(b*x+c)^4/64;B(:,:,1)=[-6/l^2+12*x/l^3 4/l-6*x/l^2];B(:,:,2)=[6/l^2-12*x/l^3 ,2/l-6*x/l^2];n=0;for i=1:2for j=1:2n=n+1;f=B(:,:,i)*D*transpose(B(:,:,j));k(:,:,n)=int(f,x,0,l);endendk11=k(:,:,1);k12=k(:,:,2);k21=k(:,:,3);k22=k(:,:,4);K=[k11 k12;k21 k22];K=simple(K);2.结果（1）bx+cK =[ (3*pi*E*(11*b^4*l^4 + 49*b^3*c*l^3 + 84*b^2*c^2*l^2 + 70*b*c^3*l +35*c^4))/(560*l^3), -(pi*E*(19*b^4*l^4 + 84*b^3*c*l^3 + 147*b^2*c^2*l^2 + 140*b*c^3*l + 105*c^4))/(1120*l^2), -(3*pi*E*(11*b^4*l^4 + 49*b^3*c*l^3 + 84*b^2*c^2*l^2 + 70*b*c^3*l + 35*c^4))/(560*l^3), -(pi*E*(47*b^4*l^4 + 210*b^3*c*l^3 + 357*b^2*c^2*l^2 + 280*b*c^3*l + 105*c^4))/(1120*l^2)][ -(pi*E*(19*b^4*l^4 + 84*b^3*c*l^3 + 147*b^2*c^2*l^2 + 140*b*c^3*l +105*c^4))/(1120*l^2), (pi*E*(3*b^4*l^4 + 14*b^3*c*l^3 + 28*b^2*c^2*l^2 +35*b*c^3*l + 35*c^4))/(560*l), (pi*E*(19*b^4*l^4 + 84*b^3*c*l^3 + 147*b^2*c^2*l^2 +140*b*c^3*l + 105*c^4))/(1120*l^2), (pi*E*(13*b^4*l^4 + 56*b^3*c*l^3 +91*b^2*c^2*l^2 + 70*b*c^3*l + 35*c^4))/(1120*l)][ -(3*pi*E*(11*b^4*l^4 + 49*b^3*c*l^3 + 84*b^2*c^2*l^2 + 70*b*c^3*l +35*c^4))/(560*l^3), (pi*E*(19*b^4*l^4 + 84*b^3*c*l^3 + 147*b^2*c^2*l^2 + 140*b*c^3*l +105*c^4))/(1120*l^2), (3*pi*E*(11*b^4*l^4 + 49*b^3*c*l^3 + 84*b^2*c^2*l^2 +70*b*c^3*l + 35*c^4))/(560*l^3), (pi*E*(47*b^4*l^4 + 210*b^3*c*l^3 + 357*b^2*c^2*l^2 + 280*b*c^3*l + 105*c^4))/(1120*l^2)][ -(pi*E*(47*b^4*l^4 + 210*b^3*c*l^3 + 357*b^2*c^2*l^2 + 280*b*c^3*l +105*c^4))/(1120*l^2), (pi*E*(13*b^4*l^4 + 56*b^3*c*l^3 + 91*b^2*c^2*l^2 +70*b*c^3*l + 35*c^4))/(1120*l), (pi*E*(47*b^4*l^4 + 210*b^3*c*l^3 + 357*b^2*c^2*l^2 +280*b*c^3*l + 105*c^4))/(1120*l^2), (pi*E*(17*b^4*l^4 + 77*b^3*c*l^3 +133*b^2*c^2*l^2 + 105*b*c^3*l + 35*c^4))/(560*l)]（2）ax^2+bx+c（将程序中D的直径换成“ax^2+bx+c”）K =[ (pi*E*(518*a^4*l^8+2233*a^3*b*l^7+2420*a^3*c*l^6+3630*a^2*b^2*l^6+7920*a^2*b*c*l^5 + 4356*a^2*c^2*l^4 +2640*a*b^3*l^5+8712*a*b^2*c*l^4+9702*a*b*c^2*l^3+3696*a*c^3*l^2 + 726*b^4*l^4 + 3234*b^3*c*l^3 + 5544*b^2*c^2*l^2 + 4620*b*c^3*l +2310*c^4))/(12320*l^3), -(pi*E*(938*a^4*l^8+4004*a^3*b*l^7+4290*a^3*c*l^6+6435*a^2*b^2*l^6+13860*a^2*b*c*l^ 5+7524*a^2*c^2*l^4+4620*a*b^3*l^5+15048*a*b^2*c*l^4+16632*a*b*c^2*l^3+6468*a*c^3*l ^2 +1254*b^4*l^4+5544*b^3*c*l^3+9702*b^2*c^2*l^2+9240*b*c^3*l+6930*c^4))/(73920*l^2), -(pi*E*(518*a^4*l^8 + 2233*a^3*b*l^7 + 2420*a^3*c*l^6 + 3630*a^2*b^2*l^6 +7920*a^2*b*c*l^5 + 4356*a^2*c^2*l^4 + 2640*a*b^3*l^5 + 8712*a*b^2*c*l^4 +9702*a*b*c^2*l^3 + 3696*a*c^3*l^2 + 726*b^4*l^4 + 3234*b^3*c*l^3 + 5544*b^2*c^2*l^2 + 4620*b*c^3*l + 2310*c^4))/(12320*l^3), -(pi*E*(2170*a^4*l^8 + 9394*a^3*b*l^7 +10230*a^3*c*l^6 + 15345*a^2*b^2*l^6 + 33660*a^2*b*c*l^5 + 18612*a^2*c^2*l^4 +11220*a*b^3*l^5 + 37224*a*b^2*c*l^4 + 41580*a*b*c^2*l^3 + 15708*a*c^3*l^2 +3102*b^4*l^4 + 13860*b^3*c*l^3 + 23562*b^2*c^2*l^2 + 18480*b*c^3*l +6930*c^4))/(73920*l^2)][ -(pi*E*(938*a^4*l^8 + 4004*a^3*b*l^7 + 4290*a^3*c*l^6 + 6435*a^2*b^2*l^6 + 13860*a^2*b*c*l^5 + 7524*a^2*c^2*l^4 + 4620*a*b^3*l^5 + 15048*a*b^2*c*l^4 +16632*a*b*c^2*l^3 + 6468*a*c^3*l^2 + 1254*b^4*l^4 + 5544*b^3*c*l^3 + 9702*b^2*c^2*l^2 + 9240*b*c^3*l + 6930*c^4))/(73920*l^2), (pi*E*(434*a^4*l^8 + 1848*a^3*b*l^7 + 1980*a^3*c*l^6 + 2970*a^2*b^2*l^6 + 6435*a^2*b*c*l^5 + 3564*a^2*c^2*l^4 +2145*a*b^3*l^5 + 7128*a*b^2*c*l^4 + 8316*a*b*c^2*l^3 + 3696*a*c^3*l^2 + 594*b^4*l^4 + 2772*b^3*c*l^3 + 5544*b^2*c^2*l^2 + 6930*b*c^3*l + 6930*c^4))/(110880*l),(pi*E*(938*a^4*l^8 + 4004*a^3*b*l^7 + 4290*a^3*c*l^6 + 6435*a^2*b^2*l^6 +13860*a^2*b*c*l^5 + 7524*a^2*c^2*l^4 + 4620*a*b^3*l^5 + 15048*a*b^2*c*l^4 +16632*a*b*c^2*l^3 + 6468*a*c^3*l^2 + 1254*b^4*l^4 + 5544*b^3*c*l^3 + 9702*b^2*c^2*l^2 + 9240*b*c^3*l + 6930*c^4))/(73920*l^2), (pi*E*(1946*a^4*l^8 + 8316*a^3*b*l^7 +8910*a^3*c*l^6 + 13365*a^2*b^2*l^6 + 28710*a^2*b*c*l^5 + 15444*a^2*c^2*l^4 +9570*a*b^3*l^5 + 30888*a*b^2*c*l^4 + 33264*a*b*c^2*l^3 + 12012*a*c^3*l^2 +2574*b^4*l^4 + 11088*b^3*c*l^3 + 18018*b^2*c^2*l^2 + 13860*b*c^3*l +6930*c^4))/(221760*l)][ -(pi*E*(518*a^4*l^8 + 2233*a^3*b*l^7 + 2420*a^3*c*l^6 +3630*a^2*b^2*l^6 + 7920*a^2*b*c*l^5 + 4356*a^2*c^2*l^4 + 2640*a*b^3*l^5 +8712*a*b^2*c*l^4 + 9702*a*b*c^2*l^3 + 3696*a*c^3*l^2 + 726*b^4*l^4 + 3234*b^3*c*l^3 +5544*b^2*c^2*l^2 + 4620*b*c^3*l + 2310*c^4))/(12320*l^3), (pi*E*(938*a^4*l^8 + 4004*a^3*b*l^7 + 4290*a^3*c*l^6 + 6435*a^2*b^2*l^6 + 13860*a^2*b*c*l^5 +7524*a^2*c^2*l^4 + 4620*a*b^3*l^5 + 15048*a*b^2*c*l^4 + 16632*a*b*c^2*l^3 +6468*a*c^3*l^2 + 1254*b^4*l^4 + 5544*b^3*c*l^3 + 9702*b^2*c^2*l^2 + 9240*b*c^3*l + 6930*c^4))/(73920*l^2), (pi*E*(518*a^4*l^8 + 2233*a^3*b*l^7 +2420*a^3*c*l^6 + 3630*a^2*b^2*l^6 + 7920*a^2*b*c*l^5 + 4356*a^2*c^2*l^4 +2640*a*b^3*l^5 + 8712*a*b^2*c*l^4 + 9702*a*b*c^2*l^3 + 3696*a*c^3*l^2 + 726*b^4*l^4 + 3234*b^3*c*l^3 + 5544*b^2*c^2*l^2 + 4620*b*c^3*l + 2310*c^4))/(12320*l^3),(pi*E*(2170*a^4*l^8 + 9394*a^3*b*l^7 + 10230*a^3*c*l^6 + 15345*a^2*b^2*l^6 +33660*a^2*b*c*l^5 + 18612*a^2*c^2*l^4 + 11220*a*b^3*l^5 + 37224*a*b^2*c*l^4 +41580*a*b*c^2*l^3 + 15708*a*c^3*l^2 + 3102*b^4*l^4 + 13860*b^3*c*l^3 +23562*b^2*c^2*l^2 + 18480*b*c^3*l + 6930*c^4))/(73920*l^2)][ -(pi*E*(2170*a^4*l^8 + 9394*a^3*b*l^7 + 10230*a^3*c*l^6 + 15345*a^2*b^2*l^6 +33660*a^2*b*c*l^5 + 18612*a^2*c^2*l^4 + 11220*a*b^3*l^5 + 37224*a*b^2*c*l^4 +41580*a*b*c^2*l^3 + 15708*a*c^3*l^2 + 3102*b^4*l^4 + 13860*b^3*c*l^3 +23562*b^2*c^2*l^2 + 18480*b*c^3*l + 6930*c^4))/(73920*l^2), (pi*E*(1946*a^4*l^8 +8316*a^3*b*l^7 + 8910*a^3*c*l^6 + 13365*a^2*b^2*l^6 + 28710*a^2*b*c*l^5 +15444*a^2*c^2*l^4 + 9570*a*b^3*l^5 + 30888*a*b^2*c*l^4 + 33264*a*b*c^2*l^3 +12012*a*c^3*l^2 + 2574*b^4*l^4 + 11088*b^3*c*l^3 + 18018*b^2*c^2*l^2 + 13860*b*c^3*l + 6930*c^4))/(221760*l), (pi*E*(2170*a^4*l^8 + 9394*a^3*b*l^7 + 10230*a^3*c*l^6 +15345*a^2*b^2*l^6 + 33660*a^2*b*c*l^5 + 18612*a^2*c^2*l^4 + 11220*a*b^3*l^5 +37224*a*b^2*c*l^4 + 41580*a*b*c^2*l^3 + 15708*a*c^3*l^2 + 3102*b^4*l^4 +13860*b^3*c*l^3 + 23562*b^2*c^2*l^2 + 18480*b*c^3*l + 6930*c^4))/(73920*l^2),(pi*E*(2282*a^4*l^8 + 9933*a^3*b*l^7 + 10890*a^3*c*l^6 + 16335*a^2*b^2*l^6 +36135*a^2*b*c*l^5 + 20196*a^2*c^2*l^4 + 12045*a*b^3*l^5 + 40392*a*b^2*c*l^4 +45738*a*b*c^2*l^3 + 17556*a*c^3*l^2 + 3366*b^4*l^4 + 15246*b^3*c*l^3 +26334*b^2*c^2*l^2 + 20790*b*c^3*l + 6930*c^4))/(110880*l)]（3）b=0（等截面梁)K =[ (3*pi*E*c^4)/(16*l^3),-(3*pi*E*c^4)/(32*l^2),-(3*pi*E*c^4)/(16*l^3),-(3*pi*E*c^4)/(32*l^2)] [ -(3*pi*E*c^4)/(32*l^2), (pi*E*c^4)/(16*l), (3*pi*E*c^4)/(32*l^2), (pi*E*c^4)/(32*l)][ -(3*pi*E*c^4)/(16*l^3), (3*pi*E*c^4)/(32*l^2), (3*pi*E*c^4)/(16*l^3),(3*pi*E*c^4)/(32*l^2)] [ -(3*pi*E*c^4)/(32*l^2), (pi*E*c^4)/(32*l), (3*pi*E*c^4)/(32*l^2), (pi*E*c^4)/(16*l)]总结：将结果（3）与教材等截面梁刚度矩阵比较，发现表达式一样，侧面证明了程序的正确性。

有限元实例分析

作业一：有限元分析实例实例：请对一个盘轴配合机构进行接触分析。

轴为一等直径空心轴，盘为等厚度圆盘，其结构及尺寸如图所示。

盘和轴为一种材料，材料参数为：弹性模量Ex=2.5E5，泊松比NUXY=0.3，摩擦系数MU=0.25，试采用有限元计算方法分析轴和盘在过盈配合时的应力应变分布以及将轴从盘心拔出时轴和盘的接触情况。

问题分析说明（1）本题主要分析装配过程中结构的静态响应，所以分析步选择通用静态分析步。

由于为过盈配合，属于大变形，故应考虑几何非线性的影响。

（2）模型具有轴对称性，所以可以采取轴对称模型来进行分析，先建立二维模型计算，再转换为三维模型计算，这样可以节省计算时间。

分析过程由两个载荷步组成, 第一个载荷步为过盈分析, 求解过盈安装时的情况。

第二个载荷步为将轴从盘心拔出时的接触分析, 分析在这个过程中盘心面和轴的外表面之间的接触应力。

它们都属于大变形问题, 属于非线性问题。

在分析时需要定义一些非线性选项来帮助问题的收敛。

（3）接触面之间有很大的相对滑动，所以模型要使用有限滑移。

模型建立的分析说明（1）进定义单元类型此项实例分析的问题中涉及到大变形, 故选用So li d185 单元类型来建立本实例入部件模块，的模型。

盘轴接触问题属于面面接触, 目标面和接触面都是柔性的,将使用接触单元T ARGET 170 和CO NTAT17 4来模拟接触面。

分别创建名为为part1、part2的部件。

（2）定义材料属性，在线性各向同性材料属性对话框中的EX (弹性模量) 文本框中输入 2 . 5E5,PRX Y (泊松比) 文本框中输入0 . 3,并将定义的材料属性赋予给part1和part2。

如下图所示。

（3）进入装配模块，创建两者间的装配关系。

（4）进入分析步模块定义名为step1和step2的两个分析步。

（5）进入相互作用模块，创建相互作用属性，设置摩擦系数；然后定义接触关系。

如下图所示。

（6）进入载荷模块，创建边界条件，依次定义名为BC -2（类型为：完全固定）、BC -3（类型为：位移/转角，约束U1、UR3），分析步均为Initial 。

有限元分析大作业

超静定梁的有限元分析本文分别通过材料力学解法和有限元解法，求出了超静定梁的支反力、最大位移及最大位移出现位置，并对两者进行了比较和误差分析。

一、超静定梁的材料力学解法梁的约束反力数目超过了有效平衡方程数，单纯使用静力平衡不能确定全部未知力的梁称为超静定梁。

超静定梁比静定梁有许多优点，如可用较少材料获得较大的刚度和强度，个别约束破坏后仍可工作等。

因而超静定梁在工程中得到较多的应用。

超静定梁的解法有很多种，本文采用力法的一种——变形比较法求解未知量。

图1图2选取C 点的支座为多余约束，Rc 为多余支座反力，则相应的基本静定梁为一外伸梁，如图2所示，其上受集中载荷P 、均布载荷q 和多余支座反力Rc 的作用。

相应的变形条件为：c cP cq cRc f f f f =++=其中316cP B Pl f l EI θ=⨯= 4724cq ql f EI =-323c cRc R l f EI =则316Pl EI 4724ql EI -+323c R l EI=0 将已知数据带入可求得 6.25c R =- 负号表示c R 的方向与假设的方向相反。

再列出平衡方程：0X =∑AX R =0A M =∑ 232022B C ql Pl R l R l ---=0C M =∑ 232022AY B ql PllR R l +--=带入已知条件求得：AX R = 393.75AY R = 812.5B R =用叠加法求最大位移：最大的向下位移在A 与B 两点中间：334410.7910481632C R l Pl ql f EI EI EI -=-++=-⨯最大的向上位移在B 与C 两点中间：3344213490.22525103248512C R l Pl ql f EI EI -=--=⨯二、超静定梁的有限元解法在ANSYS 平台上，求解超静定梁。

建模、单元划分、加载后结果如图3所示。

图3求解后可以通过图形和列表两种方式查看结果。

有限元分析报告大作业

《有限元分析》大作业基本要求:1. 以小组为单位完成有限元分析计算，并将计算结果上交；2. 以小组为单位撰写计算分析报告；3. 按下列模板格式完成分析报告；4. 计算结果要求提交电子版，一个算例对应一个文件夹，报告要求提交电子版和纸质版。

《有限元分析》大作业小组成员：储成峰李凡张晓东朱臻极高彬月Job name :banshou完成日期：2016-11-22一、问题描述（要求：应结合图对问题进行详细描述，同时应清楚阐述所研究问题的受力状况和约束情况。

图应清楚、明晰，且有必要的尺寸数据。

）如图所示，为一内六角螺栓扳手，其轴线形状和尺寸如图，横截面为一外接圆半径为0.01m的正六边形，拧紧力F为600N,计算扳手拧紧时的应力分布＜图1扳手的几何结构（要求：针对问题描述给出相应的数学模型，应包含示意图，示意图中应有必要的尺寸数据;数学模型图2数学模型如图二所示，扳手结构简单，直接按其结构进行有限元分析。

三、有限元建模3.1单元选择（要求：给出单元类型，并结合图对单元类型进行必要阐述，包括节点、自由度、实常数等。

）图3单元类型扳手截面为六边形，采用4节点182单元,182单元可用来对固体结构进行如进行了简化等处理，此处还应给出文字说明。

）Figure 1B2.1: PLANE1S2 Geometry二维建模。

182单元可以当作一个平面单元，或者一个轴对称单元。

它由4个结点组成，每个结点有2个自由度，分别在x,y方向。

扳手为规则三维实体，选择8节点185单元，它由8个节点组成，每个节点有3个自由度，分别在x，y，z方向。

3.2实常数（要求：给出实常数的具体数值，如无需定义实常数，需明确指出对于本问题选择的单元类型，无需定义实常数。

）因为该单元类型无实常数，所以无需定义实常数3.3材料模型（要求：指出选择的材料模型，包括必要的参数数据。

）对于三维结构静力学，应力主要满足广义虎克定律，因此对应ANSYS^的线性，弹性，各项同性，弹性模量EX 2e11 Pa,泊松比PRXY=0.33.4几何建模由于扳手结构比较简单，所以可以直接在ANSYS软件上直接建模，在ANSYS建立正六边形，再创立直线，面沿线挤出体，得到扳手几何模型VULUHES 1TYPE NUMEZ图4几何建模3.5网格划分方案（要求：指出网格划分方法，网格控制参数，最终生成的单元总数和节点总数，此外还应附上最终划分好的网格截图。

有限元分析大作业

有限元大作业一题目要求：图1所示为一悬臂梁，在端部承受载荷，材料弹性模量为E，泊松比为1/3，悬臂梁的厚度（板厚）为t，若该粱被划分为两个单元，单元和节点编号如图所示，试按平面应力问题计算各个节点位移计支反力。

一、单元划分1.计算简图及单元划分如下所示：2.进行节点及单元编号节点i j m单元① 2 3 4② 3 2 13.节点坐标值节点号1 2 3 4坐标值X 2 2 0 0Y 1 0 1 0二、计算单元刚度矩阵1、计算每个单元面积△以及i b ，i c （m j i i ,,=） ①②单元的面积相等，即12121=⨯⨯=∆ 单元①的i b ，i c⎩⎨⎧=--==-=0)(1m j i m j i y x c y y b ⎩⎨⎧=--==-=2)(0i m ji m j x x c y y b ⎩⎨⎧-=--=-=-=2)(1j i mj i m y x c y y b 对平面应力问题，其表达式为[]⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡-+-+-+-+∆-=s r s r sr s r s r s r s r s r b b uc c cb u b uc b c u c ub c c u b b u Et Krs 21212121)1(42 然后对单元①求解单元刚度子矩阵2==i r 2==i s []⎥⎦⎤⎢⎣⎡=3/1001329)1(22Et K 2==i r 3==j s []⎥⎦⎤⎢⎣⎡=03/23/20329)1(23Et K2==i r 4==m s []⎥⎦⎤⎢⎣⎡----=3/13/23/21329)1(24Et K 3==j r 3==j s []⎥⎦⎤⎢⎣⎡=4003/4329)1(33Et K 3==j r 2==i s []⎥⎦⎤⎢⎣⎡=03/23/20329)1(32Et K 3==j r 4==m s []⎥⎦⎤⎢⎣⎡----=43/23/23/4329)1(34Et K 4==m r 4==m s []⎥⎦⎤⎢⎣⎡=3/133/43/43/7329)1(44Et K 4==m r 2==i s []⎥⎦⎤⎢⎣⎡----=3/13/23/21329)1(42Et K 4==m r 3==j s []⎥⎦⎤⎢⎣⎡----=43/23/23/4329)1(43Et K由子矩阵[]e rs K 合成单元刚度矩阵[]⎥⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎢⎢⎣⎡----------------=3/133/443/23/13/23/43/73/23/43/2143/24003/23/23/403/43/203/13/203/23/103/213/2001329)1(Et K将单元①的单元刚度矩阵补零升阶变为单元刚度矩阵，其在总体刚度矩阵中的位置为：节点号→单元②的i b ，i c⎩⎨⎧=--=-=-=0)(1m j im j i y x c y y b ⎩⎨⎧-=--==-=2)(0i m ji m j x x c y y b ⎩⎨⎧=--==-=2)(1j i mj i m y x c y y b 然后对单元求解单元刚度子矩阵：3==i r 3==i s []⎥⎦⎤⎢⎣⎡=3/1001329)2(33Et K 3==i r 2==j s []⎥⎦⎤⎢⎣⎡=03/23/20329)2(32Et K 3==i r 1==m s []⎥⎦⎤⎢⎣⎡----=3/13/23/21329)2(31Et K 1 2 3 412[])1(22K[])1(23K[])1(24K3[])1(32K[])1(33K[])1(34K4[])1(42K[])1(43K[])1(44K2==j r 2==j s []⎥⎦⎤⎢⎣⎡=4003/4329)2(22Et K 2==j r 3==i s []⎥⎦⎤⎢⎣⎡=03/23/20329)2(23Et K 2==j r 1==m s []⎥⎦⎤⎢⎣⎡----=43/23/23/4329)2(21Et K 1==m r 1==m s []⎥⎦⎤⎢⎣⎡=3/133/43/43/7329)2(11Et K 1==m r 3==i s []⎥⎦⎤⎢⎣⎡----=3/13/23/21329)2(13Et K 1==m r 2==j s []⎥⎦⎤⎢⎣⎡----=43/23/23/4329)2(12Et K 由子矩阵[]e rs K 合成单元刚度矩阵[]⎥⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎢⎢⎣⎡----------------=3/133/443/23/13/23/43/73/23/43/2143/24003/23/23/403/43/203/13/203/23/103/213/2001329)2(Et K将单元②的单元刚度矩阵补零升阶变为单元贡献矩阵，其在总体刚度矩阵中的位置为：节点号→1 2 3 41 [])2(11K[])2(12K[])2(13K2 [])2(21K[])2(22K[])2(23K3 [])2(31K [])2(32K [])2(33K 4三、计算总体刚度矩阵总体刚度矩阵是由各单元的贡献矩阵迭加而成)2()1(][][][][K K K K e +==∑四、进行节点约束处理根据节点约束情况，在总刚矩阵中可采用划行划列处理约束的方法，由题目易知，节点3和4的已知水平位移和垂直位移都为零，划去其相对应的行和列，则总刚矩阵由8阶变为4阶，矩阵如下：⎪⎪⎭⎪⎪⎬⎫⎪⎪⎩⎪⎪⎨⎧--=⎪⎪⎭⎪⎪⎬⎫⎪⎪⎩⎪⎪⎨⎧⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡--------2/02/03/13043/203/73/23/443/23/133/43/23/43/43/73292211p p v u v u Et329][Et K =1 2 3 413/133/43/43/743/23/23/4----3/13/23/21----000243/23/23/4----3/13003/73/43/403/13/23/21----33/13/23/21----3/43/403/13003/743/23/23/4----40003/13/23/21----43/23/23/4----3/133/43/43/7化简⎪⎪⎭⎪⎪⎬⎫⎪⎪⎩⎪⎪⎨⎧--=⎪⎪⎭⎪⎪⎬⎫⎪⎪⎩⎪⎪⎨⎧⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡--------Et p Et p v u v u 3/1603/160130122072412213424472211 五、求解线性方程组方法：采用LU 分解法 1.求解矩阵[]U 各元素⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡-----⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡------⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡--------75/10775/640075/6475/353007/767/27/7502447~7/877/87/7607/87/337/207/767/27/7502447~13012207241221342447⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡----353/44900075/6475/353007/767/27/7502447~ 得到的[]U 矩阵如下:[]⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡----=353/44900075/6475/353007/767/27/7502447U 2.求解矩阵[]L 各元素⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡--------=⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡----⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡----13012207241221342447353/44900075/6475/353007/767/27/75024471353/6475/767/20175/27/40017/40001 得到的[]L 矩阵如下：[]⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡--------=13012207241221342447L3.进行求解⎪⎪⎭⎪⎪⎬⎫⎪⎪⎩⎪⎪⎨⎧---=⇒⎪⎪⎭⎪⎪⎬⎫⎪⎪⎩⎪⎪⎨⎧--=Et p Et p Et p y Et p Et p Ly 79425/850800225/323/1603/1603/160⎪⎪⎭⎪⎪⎬⎫⎪⎪⎩⎪⎪⎨⎧---=⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡----⇒=⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡Et p Et p Et p v u v u y v u v u U 79425/850800225/323/160353/44900075/6475/353007/7675/27/750244722112211 解得Et p v /422.82-= Et p u /497.12-= Et p v /028.91-= Et p u /897.11=于是求得各节点的位移为：⎩⎨⎧-==Etp v Etp u /028.9/897.111 ⎩⎨⎧-=-=Etp v Etp u /422.8/497.122 ⎩⎨⎧==033v u ⎩⎨⎧==044v u 六、求解相应的支反力（运用静力学的平衡方程进行求解）3号节点和4号节点的支反力如下图所示：。

有限元分析大作业

《有限元分析及应用》大作业——齿根弯曲应力计算报告班级：无可奉告姓名：无可奉告学号：无可奉告指导老师：无可奉告目录目录 (2)1.概述 (3)1.1工程问题描述 (3)1.2问题分析 (3)2.建模过程 (4)2.1几何建模 (4)2.2CAE网格划分与计算 (5)2.3后处理 (8)3.多方案比较与结果分析 (9)3.1多方案比较 (9)3.2结果分析 (11)1.概述1.1工程问题描述我在本次作业中的选题为齿根弯曲应力的计算与校核。

通过对机械设计的学习，我们可以知道，齿轮的失效形式主要是齿面接触疲劳和齿根弯曲断裂，而闭式传动硬齿面齿轮的失效形式以齿根弯曲断裂，这个时候进行齿根弯曲应力的校核才比较有意义，在设计问题的时候应当选取这种类型的算例。

设计计算的另一个主要思路是将有限元计算的结果与传统机械设计的结算结果进行对比，以从多方面验证计算结果的准确性。

综上，我们最终选取了《机械原理》（第三版）P50例3-1中的问题进行校核计算。

已知起重机械用的一对闭式直齿圆柱齿轮，传动，输入转速n1=730r/min,输入功率P1=35kW，每天工作16小时，使用寿命5年，齿轮为非对称布置，轴的刚性较大，原动机为电动机，工作机载荷为中等冲击。

z1=29，z2=129，m=2.5mm，b1=48mm，b2=42mm，大、小齿轮均为20CrMnTi，渗碳淬火，齿面硬度为58~62HRC，齿轮精度为7级，试验算齿轮强度。

齿面为硬齿面，传动方式为闭式传动。

根据设计手册查出的许用接触应力为1363.6Mpa，计算结果为1260Mpa，强度合格。

根据设计手册查出的许用弯曲应力为613.3MPa，计算结果为619Mpa，强度略显不够。

1.2问题分析大小齿轮啮合，小齿轮受载荷情况较为严峻，故分析对象应当为小齿轮。

可以看出，由于齿轮单侧受载荷，传动过程中每个齿上载荷的变化过程是相同的，故问题可被简化为反对称问题，仅需研究单个齿。

有限元热分析第1次大作业(16-17)

课程名称：有限元与热分析数值仿真
2016—2017学年第二学期
第一次大作业（任课教师：钱作勤）
1、简述工程热力学的三大基本定律，并深刻阐述其重要意义和应用领域。

2、一台小型化工装置采用水蒸汽再热循环。

透平蒸汽进口参数为8.5Mpa和480
度。

再热参数为1.2Mpa和440度。

凝汽压力为7kPa.。

透平和泵的效率分别为0.92和0.80。

画出该循环的T-S图。

并确定：1）每千克工质的净功；2）再热占总吸热量的百分比；3）循环热效率。

3、一稳定运行的理想蒸汽压缩制冷系统采用R134a作为工作流体，压缩机进口
是压力为0.16Mpa的饱和蒸汽。

冷凝器出口参数为0.9Mpa和32度。

质量流量为5kg/min。

压缩机等熵效率为80%。

试确定：1）压缩机功率；2）冷吨；
3）性能系数。

有限元分析及应用大作业

有限元分析及应用大作业作业要求：1）个人按上机指南步骤至少选择习题中3个习题独立完成，并将计算结果上交；也可根据自己科研工作给出计算实例。

2）以小组为单位完成有限元分析计算；3）以小组为单位编写计算分析报告；4）计算分析报告应包括以下部分：A、问题描述及数学建模；B、有限元建模（单元选择、结点布置及规模、网格划分方案、载荷及边界条件处理、求解控制）C、计算结果及结果分析（位移分析、应力分析、正确性分析评判）D、多方案计算比较（结点规模增减对精度的影响分析、单元改变对精度的影响分析、不同网格划分方案对结果的影响分析等）题一：图示无限长刚性地基上的三角形大坝，受齐顶的水压力作用，试用三节点常应变单元和六节点三角形单元对坝体进行有限元分析，并对以下几种计算方案进行比较：1）分别采用相同单元数目的三节点常应变单元和六节点三角形单元计算；（注意ANSYS中用四边形单元退化为三节点三角形单元）2）分别采用不同数量的三节点常应变单元计算；3）当选常应变三角单元时，分别采用不同划分方案计算。

解：1.建模：由于大坝长度>>横截面尺寸，且横截面沿长度方向保持不变，因此可将大坝看作无限长的实体模型，满足平面应变问题的几何条件；对截面进行受力分析，作用于大坝上的载荷平行于横截面且沿纵向方向均匀分布，两端面不受力，满足平面应变问题的载荷条件。

因此该问题属于平面应变问题，大坝所受的载荷为面载荷，分布情况P=98000-9800*Y；建立几何模型，进行求解；假设大坝的材料为钢，则其材料参数：弹性模量E=，泊松比σ=；2：有限元建模过程：进入ANSYS :程序→ANSYS APDL设置计算类型 :ANSYS Main Menu: Preferences →select Structural →OK选择单元类型 :ANSYS Main Menu: Preprocessor →Element Type→Add/Edit/Delete →Add →select Solid Quad 4node 182(三节点常应变单元选择Solid Quad 4node 182，六节点三角形单元选择Solid Quad 8node 183) →OK (back to Element Types window) →Option →select K3: Plane Strain →OK→Close (the Element Type window)定义材料参数：ANSYS Main Menu: Preprocessor →Material Props →Material Models →Structural →Linear →Elastic →Isotropic →input EX:, PRXY: →OK生成几何模型：生成特征点：ANSYS Main Menu: Preprocessor →Modeling →Create →Keypoints →In Active CS →依次输入四个点的坐标：input:1(0,0),2(10,0),3(1,5),4,5) →OK生成坝体截面：ANSYS Main Menu: Preprocessor →Modeling →Create →Areas →Arbitrary →Through KPS →依次连接四个特征点，1(0,0),2(6,0),3(0,10) →OK网格划分：ANSYS Main Menu: Preprocessor →Meshing →Mesh Tool→(Size Controls) lines: Set →依次拾取两条直角边:OK→input NDIV: 15 →Apply→依次拾取斜边:OK →input NDIV: 20 →OK →(back to the mesh tool window)Mesh:Areas, Shape: tri, Mapped →Mesh →Pick All (in Picking Menu) →Close( the Mesh Tool window)模型施加约束：给底边施加x和y方向的约束：ANSYS Main Menu: Solution →Define Loads →Apply →Structural →Displacement →On lines →pick the lines →OK →select Lab2:UX, UY →OK给竖直边施加y方向的分布载荷：ANSYS 命令菜单栏: Parameters →Functions →Define/Edit →1) 在下方的下拉列表框内选择x ,作为设置的变量；2) 在Result窗口中出现{X},写入所施加的载荷函数：98000-9800*{Y}； 3) File>Save(文件扩展名：func) →返回：Parameters →Functions →Read from file：将需要的.func文件打开，参数名取meng，它表示随之将施加的载荷→OK →ANSYS Main Menu: Solution →Define Loads →Apply →Structural →Pressure →On Lines →拾取竖直边；OK →在下拉列表框中，选择：Existing table →OK →选择需要的载荷为meng参数名→OK分析计算：ANSYS Main Menu: Solution →Solve →Current LS →OK(to close the solve Current LoadStep window) →OK结果显示：ANSYS Main Menu: General Postproc →Plot Results →Deformed Shape… → select Def + Undeformed →OK (back to Plot Results window)→Contour Plot →Nodal Solution→select:DOF solution, UX,UY,Def + Undeformed , Stress ,SX,SY,SZ, Def + Undeformed →OK：结果分析：（第一小题）图为采用三节点常应变单元，200个单元，左下图为位移变形图；右下图为应力变形图图为采用六节点三角形单元，200个单元，左下图为位移变形图；右下图为应力变形图根据位移和应力图可得：单元类型Min(位移)Max(位移)Min(应力)Max(应力)常应变三节点05461392364六节点三角形0607043对比分析可得：最大位移都发生在A点，即大坝顶端，最大应力发生在B点附近，即坝底和水的交界处，且整体应力和位移变化分布趋势相似，符合实际情况；结果显示三节点和六节点单元分析出来的最大应力值相差较大根据结果显示，最小三节点和六节点单元分析出来的最小应力值相差极为悬殊，结合理论分析，实际上A点不承受载荷，最小应力接近于零，显然六节点三角形单元分析在这一点上更准确。

有限元分析大作业

一、有限元方法的手工计算结果与ansys分析结果的对比1分析的问题描述如图1所示，桁架的杆截面面积为8，由钢制成（E=200GPa）。

用有限元法计算出每个节点的位移以及反作用力。

(1)(2)(3)图1对于上述问题，本文将用手工计算和ansys软件分别计算出结果，对计算出来的结果进行对比。

2手工计算2.1桁架结构的有限元计算方法对于桁架结构，每个单元的刚度矩阵为，（2-1）YX图2其中，为桁架单元在整体坐标系中与X轴的夹角；，A为桁架的截面积，E 为弹性模量，L为桁架长度。

在固体力学问题中，有限元公式通常由如下的一般形式，Ku=F（2-2）其中，K为刚度矩阵，u为位移矩阵，F为载荷矩阵。

运用公式（2-3），就能求出反作用力，R=Ku-F（2-3）其中，R为反作用力矩阵。

2.2计算过程计算每个桁架单元的刚度，用公式（2-1）计算每个每个桁架单元的刚度矩阵，将每个单元放入总刚度矩阵，他们的位置分别为：10-100000 00000000 -10100000 00000000 00000000 00000000 00000000 0000000000000000 00000000 0010-1000 0000000000-101000 00000000 00000000 000000003.9-4.90000-3.9 4.9-4.9 6.10000 4.9-6.1 00000000 00000000 00000000 00000000-3.9 4.90000 3.9-4.9 4.9-6.10000-4.9 6.100000000 00000000 00000000 000 1.28000-1.28 00000000 00000000 00000000 000-1.28000 1.2800000000 00000000 00000000 00000000 0000 3.9 4.9-3.9-4.9 0000 4.9 6.1-4.9-6.1 0000-3.9-4.9 3.9 4.9 0000-4.9-6.1 4.9 6.1将个刚度矩阵相加得到总刚度矩阵为，19.9-4.90-16000-3.9 4.9 -4.9 6.100000 4.9-6.1 -160320-16000 00012.8000-12.8 00-16019.9 4.9-3.9-4.9 0000 4.9 6.1-4.9-6.1 -3.9 4.900-3.9-4.97.80 4.9-6.10-12.8-4.9-6.1025应用边界条件施加载荷，将总刚度矩阵带入式（2-2）得：19.9-4.90-16000-3.9 4.9Ux1 -4.9 6.100000 4.9-6.1Uy1 -160320-16000Ux2 00012.8000-12.8Uy200-16019.9 4.9-3.9-4.9Ux3 0000 4.9 6.1-4.9-6.1Uy3 -3.9 4.900-3.9-4.97.80Ux4 4.9-6.10-12.8-4.9-6.1025Uy4带入边界条件解得：将结果带入（2-3）得：=Fx1Fy1Fx2Fy2Fx3Fy3Fx4Fy43用ansys软件求解（单位统一N,mm，Mpa）（1）选择单元（图3）图3（2）附材料属性（图4）图4（3）创建模型（图5）图5（4）施加载荷（图6）图6（5）求解每个节点的位移（图7）图7节点的反力（图8）图8（6）模型变形图（7）位移等值线分布图4结果对比及分析手算结果ansys 计算结果位移(mm)Ux100Uy100Ux2-0.0016-0.0016Uy2-0.0468-0.0468Ux300Uy300Ux4-0.0066-0.0066Uy4-0.0317-0.0317表1手算结果ansys计算结果节点反力(N)Fx1-1027.8-1027.8 Fy11608.31608.3 Fx2 5.60 Fy2-100 Fx32066.72063.1 Fy32257.12255.4 Fx48.80 Fy4-2.70表2由表1和表2可以看出，手工计算的结果与ansys计算的结果基本一致。

有限元受力分析大作业

１. 定义分析类型
Main Menu >Solution > Analysis Type > New Analysis, 选择Static 2. 定义位移约束
Main Menu >Solution> Define Loads >Apply >Structural>Displacement > On Areas, 在弹出对话框中选Circle,点击OK后，在弹出的对话框中填入：0
后处理
1. 查看总体变形
Main Menu >General Postproc >Plot Results > Deformed shape
（2）
（2）
四. 加载和求解
１. 定义分析类型
Main Menu >Solution > Analysis Type > New Analysis, 选择Static 2. 定义位移约束
Main Menu >Solution> Define Loads >Apply >Structural>Displacement > On Areas, 在弹出对话框中选Circle,点击OK后，在弹出的对话框中填入：0
七. 后处理
1. 查看总体变形
Main Menu >General Postproc >Plot Results > Deformed shape
S
作业三：轴承座盈利分析
1、题目
2、建模
3、ansys前处理
导入
ansys前处理
1、定义单元类型 Main Menu>Proprocessor>Element Type >Add/Edit/Delete 弹出对话框中后，点“Add” 。双弹出对话框，选“Solid”和“10node 92”,点 “OK”，退回到前一个对话框。

有限元分析技术课程大作业

有限元分析技术课程大作业1 工程介绍现需要对某露天大型玻璃平面舞台的钢结构进行分析，该钢结构布置在xy 平面内。

学生序号为079，分格的列数（x向分格）=0×10+7+5=12，分格的行数（y向分格）=9+4=13，共有156个分格。

每个分格x方向尺寸为1m，y方向尺寸为1m。

钢结构的主梁为高160宽100厚14的方钢管；次梁为直径60厚10的圆钢管（单位为毫米），材料均为碳素结构钢Q235；该结构固定支撑点位于左右两端主梁和最中间的次梁的两端。

玻璃采用四点支撑与钢结构连接（采用四点支撑表明垂直作用于玻璃平面的面载荷将传递作用于玻璃所在钢结构分格四周的节点处，表现为点载荷；试对在垂直于玻璃平面方向的22/KN m的面载荷（包括玻璃自重、钢结构自重、活载荷（人员与演出器械载荷）、风载荷等）作用下的舞台进行有限元分析（每分格面载荷对于每一支撑点的载荷可等效于0.5KN的点载荷）。

作业提交的内容：（1）屏幕截图显示该结构的平面布置结构；（2）该结构每个支座的支座反力；（3）该结构节点的最大位移及其所在位置；（4）对该结构中最危险单元（杆件）进行强度校核。

2有限元模型的建立该钢结构中每一分格x方向尺寸为1m，y方向尺寸为1m，x方向分格数量为12，y方向分格数量为13。

该钢结构由主梁和次梁构成，其中主梁为高160mm、宽100mm、厚14mm的方钢管，次梁为直径60mm、厚10mm的圆钢管。

由于在该结构中所有构件均为梁单元，而Ansys程序中提供了多种梁单元，以模拟不同场合的应用，且对于每种梁单元类型都有特定的算法。

在本次建模过程中，考虑到需要对该结构中的危险单元进行强度校核，因此，选择了BEAM188单元类型来建立本钢架结构，进而对其进行有限元分析。

BEAM188为三维线性有限应变梁单元，该单元基于铁木辛哥的梁结构理论，考虑了剪切变形的影响，能够满足本次分析的需求。

以下为基于ANSYS图形界面(Graphic User Interface , GUI)的菜单操作流程(1) 进入ANSYS（设定工作目录和工作文件）程序→ ANSYS → ANSYS Interactive → Working directory（设置工作目录）→Initial Jobname（设置工作文件名）:Analysis → Run → OK(2) 设置计算类型ANSYS Main Menu：Preferences → Structural → OK(3) 定义单元类型ANSYS Main Menu：Preprocessor → Element Type → Add/Edit/Delete... → Add → Beam: 3D 2node 188 → OK（返回到Element Types窗口）→ Close(4) 定义材料参数ANSYS Main Menu: Preprocessor → Material Props → Material Models → Structural → Linear → Elastic → Isotropic → input EX: 2.0E5, PRXY: 0.3(定义泊松比及弹性模量) → OK → Close（关闭材料定义窗口）(5)定义梁单元截面ANSYS Main Menu：Preprocessor →Sections→Beam→Common Sections→Beam Tool(6) 构造梁模型生成舞台几何模型ANSYS Main Menu：Preprocessor → Modeling → Create → Keypoints → In Active CS → NPT Keypoint number：1，X，Y，Z Location in active CS：0，0，0 → Apply 通过复制关键点操作，形成14行13列的关键点。

滑轮有限元分析-ANSYS FEM 大作业

滑轮有限元分析-ANSYS FEM 大作业1.问题描述某滑轮结果如下图所示,试分析结构在实际工作中的受力情况,并利用FEM类软件校核材料的强度是否满足要求。

其中天车最大钓钩载荷为3150KN，游动系统以及钢丝绳总重为150KN。

材料为Q345。

2.问题分析天车最大钓钩载荷为3150KN，游动系统以及钢丝绳总重为150KN，游车与天车选用6x7轮系，钢丝绳实际最大拉力F=（3150+150）/12=275KN。

滑轮受力图如下图所示，当钢丝绳两端拉力平行，滑轮受力最大为2F=550KN。

图1 滑轮受力分析滑轮上端面与绳索接触，所有滑轮外表面的上半面受力，且载荷不是均匀分析，而是按照正弦函数分析。

同时滑轮内表面的上半面受力，下半面为自由状态。

在有限元分析中，需要注意选择合适的边界条件和载荷加载。

有限元分析（FEA，Finite Element Analysis）利用数学近似的方法对真实物理系统（几何和载荷工况）进行模拟。

利用简单而又相互作用的元素（即单元），就可以用有限数量的未知量去逼近无限未知量的真实系统。

有限元分析是用较简单的问题代替复杂问题后再求解。

它将求解域看成是由许多称为有限元的小的互连子域组成，对每一单元假定一个合适的（较简单的）近似解，然后推导求解这个域总的满足条件（如结构的平衡条件），从而得到问题的解。

因为实际问题被较简单的问题所代替，所以这个解不是准确解，而是近似解。

由于大多数实际问题难以得到准确解，而有限元不仅计算精度高，且能适应各种复杂形状，因而成为行之有效的工程分析手段3.求解步骤滑轮材料为Q345，根据API Spec 8C-2012第5版4.7规定滑轮的设计安全系数不小于3，所以滑轮的许用应力为115Mpa，其中弹性模量为2.1e11Pa，泊松比为0.3。

对滑轮结构进行有限元网格划分，滑轮存在较多倒角过度细节，所以采用四面体网格进行划分，对倒角圆孔区域进行局部加密，有限元网格模型如下图所示。

有限元实例分析01

有限元实例练习分析学号：0905010888姓名：刘广宇专业：材料成型及控制工程2012年5月2日引言有限元方法发展到今天。

已经成为一门相当复杂的实用工程技术。

有限元分析的最终目的是还原一个实际工程系统的数学行为特征。

即分析必须针对一个物理原型准确的数学模型。

模型包括所有节点、单元、材料属性、实常数、边界条件以及其他用来表现这个物理系统的特征。

Marc是一种融结构、热、流体、电磁和声学于一体的非线型有限元分析软件，可广泛应用于航空航天、汽车、造船、石油化工、铁道、能源、电子元件、机械制造、材料工程、土木工程、医疗器材、冶金工艺和家用电器等。

该软件功能特色具体包括：多种物理场的分析能力、复合场的耦合分析能力、强大的非线性分析能力、最先进的接触分析功能、并行计算功能、丰富的单元库、开放的用户环境、强大的网格自适应功能和全自动三维网格重划分；Marc的学习、应用是一个系统、复杂的工程。

由于它涉及到多方面的知识，所以在学Marc的过程中一定要对Marc所涉及到的一些理论知识有一个大概的了解，以加深对Marc的理解。

目录引言一、目的 (4)二、软件应用介绍 (4)三、实例内容 (6)四、求解步骤 (6)1. 建立有限元模型 (6)2. 加载求解 (11)3、后处理 (12)五、总结 (16)参考文献有限元实例练习分析一、目的1、熟悉有限元建模、求解及结果分析步骤和方法。

2、能利Marc软件对实例结构进行静力有限元分析。

3、加深有限元理论关于网格划分概念、划分原则等的理解。

二、软件应用介绍有限元分析是对于结构力学分析迅速发展起来的一种现在计算方法。

它是50年代首先在连续体力学领域—飞机结构静、动态特性分析中应用的一种有效的数值分析方法，随后很快广泛的应于求解热传导、电磁场、流体力学等连连续行问题。

想要解答，必须先简化结构，采用数值模拟方法分析。

（一）有限元软件发展特点1. 单一场计算向多物理耦合场问题的求解发展2. 由求解线性问题发展到求解非线性问题3. 与CAD/CAM等软件的集成4. 提高自动化的网格处理能力5．软件面向专业用户的开放性6. 软件开发强强联合(二)、分析研究过程1、前处理（1）建模有限元分析的最终目的是还原一个实际工程系统的数学行为特征，即分析必须针对一个物理原型准确的数学模型。

有限元分析题及大作业题答案

有限元分析及应用作业报告试题10一、问题描述确定图示扳手中的应力, E=210Gpa，μ=0.3, 假设厚度为10mm；并讨论采用何种处理可降低最大应力或改善应力分布。

图1为扳手的基本形状和基本尺寸图二、数学建模与分析由图1及问题描述可知，板手的长宽尺寸远远大于厚度，研究结构为一很薄的等厚度薄板，满足平面应力的几何条件；作用于薄板上的载荷平行于板平面且作用在沿厚度方向均匀分布在办手柄的左边缘线，而在两板面上无外力作用，满足平面应力的载荷条件。

故该问题属于平面应力问题，薄板所受的载荷为面载荷，分布情况及方向如图1所示，建立几何模型，并进行求解。

薄板的材料为钢，则其材料参数：弹性模量E=2.1e11，泊松比σ=0.3三、有限元建模1、单元选择：选取三节点常应变单元来计算分析薄板扳手的位移和应力。

由于此问题为平面应力问题，：三节点常应变单元选择的类型是PLANE42（Quad 4node42），该单元属于是四节点单元类型，在网格划分时可以对节点数目控制使其蜕化为三节点单元。

2、定义材料参数：ANSYS Main Menu: Preprocessor →Material Props →Material Models →Structural →Linear →Elastic →Isotropic →input EX:2.1e11, PRXY:0.3 →OK3、生成几何模型：a.创建关键点点：ANSYS Main Menu: Preprocessor →Modeling →Create →Keypoints →In Active CS →依次输入16个点的坐标→OKb、将这16个关键点有直线依次连起来，成为线性模型4、生成实体模型：ANSYS Main Menu: Preprocessor →Modeling →Create →Areas →Arbitrary →Through KPS →连接特征点→生成两个area→Operate→Subtract→拾取整个扳手区域→OK→生成扳手模型5、结点布置及规模6、网格划分方案ANSYS Main Menu: Preprocessor →Meshing →Mesh Tool →Mesh: Areas, Shape: Tri,Free →Mesh →Pick All (in Picking Menu) →Close( the Mesh Tool window)7、载荷及边界条件处理8、求解控制A、模型施加约束给模型施加x方向约束ANSYS Main Menu: Solution →Define Loads →Apply →Structural →Displacement →On Lines →拾取模型左部的竖直边→OKB、给模型施加载荷ANSYS Main Menu: Solution →Define Loads →Apply →Structural →force→on keypoints→拾取上面左端关键点→700N/mm→okC、分析计算：ANSYS Main Menu: Solution →Solve →Current LS →OK(to close the solve Current Load Step window) →OK6)结果显示：ANSYS Main Menu: General Postproc →Plot Results →Deformed Shape… →select Def + Undeformed →OK (back to Plot Results window) →Contour Plot →Nodal Solu →select: DOF solution →displacement vector sum，von mises stress→OK四、计算结果及结果分析1、三节点常应变单元1）三节点单元的网格划分图2 常应变三节点单元的网格划分平面图图3 常应变三节点单元的网格划分立体图2）三节点单元的约束受载情况图4 常应变三节点单元的约束受载图3)三节点单元的位移分析图5 常应变三节点单元的位移分布图4)三节点单元的应力分析图6 常应变三节点单元的应力分布图2、六节点三角形单元1)六节点三角形单元网格划分图7 六节点三角形单元网格划分图2)六节点三角形单元约束和受载情况分析图8 六节点三角形单元约束受载图3）六节点三角形单元位移分析图9 六节点三角形单元的变形分布图4) 六节点三角形单元的应力分析图9 六节点三角形单元的应力分布图图10 六节点三角形单元的局部应力分布图根据以上位移和应力图，可以得出常应变三节点单元和六节点三角形单元的最小最大位移应力如表1-1所示。

有限元大作业

有限元大作业第一篇：有限元大作业有限元应力分析报告大作业机械与运载工程学院车辆四班龙恒 20110402415 2014年8月30日一、问题描述桦木板凳材料参数如图形状参数：长40mm，宽30mm，高45mm（其他详细参数见零件图）通过施加垂直于板凳上表面的均匀载荷600N，分析板凳的应变和应力？二、使用inventor进行建模及应力分析1、通过inventor建立板凳3D模型利用草图拉伸等方法建立与零件图中尺寸一致的三维立体板凳模型2、点选环境下的应力分析开始对板凳进行应力分析3、根据所给条件设置材料等参数、将安全系数设为屈服强度，因为板凳主要受压变形点开“木材（桦木）”根据前面所给参数对其进行参数设置4、固定约束如图板凳的4个脚底面设置为固定约束，使得板凳受载后，脚底面不会沿垂直方向位移，模拟真实情况5、施加载荷在板凳上表面施加大小为600N的垂直均布载荷（这里是模拟一个成人坐上去的重力）6、划分网格通过设置网格的尺寸参数来划分出5种不同网格数量，从而得出5种不同网格数划分得出的应力应变分布图，最后分析划分不同网格数对结果的影响。

（1）网格最大（2）网格较大（3）网格一般大小（4）网格较小（5）网格最小7、求解得出结果得出5组不同网格数所得数据（应力云图，应变云图，所有结果数据）（1）网格数1437根据应力云图可知，红色地方所受的应力最大，最大应力为：15.48Mpa 根据应变云图可知，红色地方的应变最大，最大应变为：0.001434μl（2）网格数8651根据应力云图可知，红色地方所受的应力最大，最大应力为：18.88Mpa 根据应变云图可知，红色地方的应变最大，最大应变为：0.001755μl（3）网格数20484根据应力云图可知，红色地方所受的应力最大，最大应力为：22.62Mpa 根据应变云图可知，红色地方的应变最大，最大应变为：0.002103μl（4）网格数41578根据应力云图可知，红色地方所受的应力最大，最大应力为：23.76Mpa 根据应变云图可知，红色地方的应变最大，最大应变为：0.002206μl（5）网格数68788根据应力云图可知，红色地方所受的应力最大，最大应力为：25.97Mpa 根据应变云图可知，红色地方的应变最大，最大应变为：0.002454μl综合上述5种请况可知随着网格的细分，所得的应变以及应力的结果是收敛的。

有限元分析大作业

基于ＡNSYS软件的有限元分析报告机制1２05班杜星宇U201２10671一、概述本次大作业主要利用ANSYS软件对桌子的应力和应变进行分析,计算出桌子的最大应力和应变.然后与实际情况进行比较，证明分析的正确性,从而为桌子的优化分析提供了充分的理论依据,并且通过对ANSYS软件的实际操作深刻体会有限元分析方法的基本思想,对有限元分析方法的实际应用有一个大致的认识。

二、问题分析已知：桌子几何尺寸如图所示,单位为mm。

假设桌子的四只脚同地面完全固定,桌子上存放物品,物品产生的均匀分布压力作用在桌面,压力大小等于３00Ｐa，其中弹性模量E＝9。

3GPa,泊松比μ=０。

３5,密度ρ=560kg/ｍ3，分析桌子的变形和应力．将桌脚固定在地面，然后在桌面施加均匀分布的压力，可以看作对进行平面应力分析，桌脚类似于梁单元。

由于所分析的结构比较规整且为实体，所以可以将单元类型设为八节点六面体单元。

操作步骤如下：1、定义工作文件名和工作标题（1)定义工作文件名:执行UｔｉｌitｙＭenu/ File/ChangｅＪoｂnａme,在弹出Chaｎｇe Ｊobnａme 对话框修改文件名为Tａblｅ。

选择Ｎew log aｎｄerrｏｒｆiｌeｓ复选框。

(2)定义工作标题：Ｕtility Ｍenu／File／Ｃhaｎｇe Tｉtle，将弹出ChangｅTit ｌｅ对话框修改工作标题名为Tｈｅ aｎａlysiｓ of taｂｌe。

（3)点击:Pｌｏt/Ｒｅplot。

2、设置计算类型（１）点击:Main Menu/Preferencｅs，选择Structurａl，点击OK。

3、定义单元类型和材料属性(1)点击:Ｍain Meｎu／Prｅprｏcｅssｏr／Ｅlement Type/Add/Edit/Deleｔｅ,点击Aｄd，选择Solid〉Brｉｃk ８ｎode 185,点击OK,点击Clｏse。

（2)点击Main ｍenu/preｐｒｏceｓsｏr/Matｅｒiaｌ Prｏps/Maｔerial Models / Struｃtural／ Liｎｅar/ Elａstiｃ／Ｉsotropic，设置EＸ为9.３e9,PRXY为0。