有限元受力分析结构梁力计算.docx

1、有限元受力分析结构梁力计算绪论.2第一章有限元课程设计 . 4一.工程问题 .4二.简化模型 .4三.解析法求解 .5四.ANSYS求解 .8五.结果分析 .19第二章机械优化设计说明 .20一.题目及解析 . 20二.黄金分割法计算框图 . 23三.C语言程序 . 24四.运行结果 . 27 五.结果分析 . 27第三章设计感言 . 28第四章参考文献 . 28前言有限元法在解决圣维南扭转问题近似解时首先提出的。有限元在弹性力学平面问题的第一个成功应用是由美国学者于1956年解决飞机结构强度时提出的、经过几十年得发展，有限元一惊成为现代结构分析得有效方法和主要手段。它的应用已经从弹性力学的平

2、面问题扩展到空间问题和板壳问题。对于有限元法，从选择基本未知量的角度来看，他可以分为三种方法：位移法，力法，混合法。从推导方法来看，它可以分为直线法，变分法，加权余数法。但随后随着计算机的发展，有限元法如虎添翼。国内外已有许多大型通用的有限元分析程序，并已经出现了将人工智能技术引入有限元分析软件，形成了比较完善得专家系统，逐步实现了有限元的智能化。 优化设计是现代设计方法的重要内容之一。它以数学规划为理论基础以电子计算机为工具，在充分考虑多种设计约束的前提下，寻求满足预订目标的最佳设计。优化设计理论于方法用于工程设计是在六十年代后期开始的，特别是今年来，随着有限元素法，可靠性设计，计算机辅助设

3、计的理论与发展及优化设计方法的综合应用使整个工程设计过程逐步向自动化集成化智能化发展，其前景使令人鼓舞的。因而工程设计工作者必须适应这种发展变化，学习，掌握和应用优化设计理论与方法。 今年来随着计算机技术的普及和计算速度的不断提高，有限元分析在工程设计和分析中得到了越来越广泛的重视，已经成为解决复杂的工程分析计算问题的有效途径，现在从汽车到航天飞机几乎所有的机械制造都已离不开有限元分析计算，其再机械制造，材料加工，航空航天，汽车，土木建筑，电子电器，国防军土，船舶，铁道，石化能源，科学研究等各个领域的广泛使用已使设计水平发生了质的飞跃，主要表现在以下几个方面： 增加产品和工程的可靠性 在产品的

4、设计阶段发现潜在的问题 经过分析计算，采用优化设计 方案，降低原材料成本 缩短产品投向市场的时间模拟试验方案，减少试验次数，从而减少试验经费 ANSYS软件致力于耦合场的分析计算，能够进行结构，流体，热，电磁四种场的计算，已博得了世界上数千家用户的钟爱。ADINA非线性有限元分析软件由著名的有限元专家，麻省理工学院的K.J.Bathe教授领导开发，其单一系统即可进行结构，流体，热的耦合计算。并同时具有隐式和显式两种时间积分算法。犹豫其在非线性求解，流固耦合分析等发面的强大功能，迅速成为有限元分析软件的后起之秀，现已成为非线性分析得首选软件。第一章有限元课程设计 一、工程问题 材料力学辅导及习题

5、精解（陈平 第四版） P180受力如图所示的梁，试根据 q,Fs,M之间的微分关系和积分关系画出Fs,M图。二、简化模型 1、梁的参数： 长度 l=10m； 宽度 b=1m 厚度 h=0.5m 2、材料参数 梁选择线性、弹性、各向同性的材料。它的弹性模量EI=207e5 Pa.。 3、梁的边界条件 在节点A处梁受Y方向的约束；节点B受X方向和Y方向的约束。 4、梁的载荷 FB之间作用着均布载荷q=2kN/m，在节点E处作用着集中力偶Me=10kNm。D 处作用沿Y负方向的集中力 F=1 KN,在C处作用沿Y负方向的集中力F=3KN。 经上述分析，此外伸梁可简化成一个二维单元，其中一端固支，一端

6、外伸。另一个支撑点变化为一个绞支点，在梁的上部分布着均布载荷。如图所示(a) 3、解析法求解（1）先求支座反力，得Fa=6KN Fb=6Kn(2)做Fs图如图所示 从梁左端开始。由C截面有向下的集中力F1的作用，所以Fs图有零开始向下突变，其突变值为F1=-3KN。由于AC段内无分布载荷作用，所以C,A两截面的剪力增量为零。CA段得剪力图为一条水平直线，并从C点一直延伸到A点稍偏左的截面处。由于A截面处有向上的集中力Fa的作用，Fs图有一向上的突变，值为Fa=+6KN，所以A右截面的剪力值为-3+6=+3KN.在AD段，又因为没有均布载荷所以Fs图也为一条直线，从A点一直延伸到D点稍左的截面处

7、。由于D截面处有一向下的集中力F2的作用，所以Fs图有一向下的突变，值为F2=-1KN，所以D右截面的剪力值为3-1=+2KN。在DEF段，因为均布载荷也为零所以Fs图为一水平直线，从D点一直延伸到F截面处。由于在FB段有向下的均布载荷作用，q=-2KN/m，所以Fs图为一下降的斜直线。BF两截面上的剪力之差等于载荷q(x)图的面积，即-2乘4=-8KN，所以B左截面上的剪力值为2-8=-6KN.最后在B截面由于有向上的集中力Fb的作用，所以Fs图有一个向上的突变，值为Fb=+6KN于是，Fs图画出如图所示(b)。(3)做M图（如图所示（c） 因为C截面为自由端，又没有集中力偶作用，所以 Mc

8、=0；弯矩图从0开始。在CA段，因为没有均布载荷，又因为Fs=-30，所以M图为一上升的斜直线；而D A两截面上的弯矩之差，等于AD段的Fs图的面积，即3乘1=3 KNm所以D截面的弯矩值为-6+3=-3 KNm。在DE段，均布载荷为0，Fs=2KN0，所以M图也为一个上增的斜直线；而E D两截面上的弯矩之差等于DE段Fs 图的面积，即2乘1=2 KNm，所以E左截面的弯矩值为-3+2=-1 KNm。在E截面，因为有一个顺时针转向的集中力偶作用，所以M图有一向上的突变，值为Me=5 KNm，所以E右截面的弯矩值为-1+5=+4 KNm。在EF段，因为没有均布载荷，又因为Fs=2KN0,所以M图

9、为一上增的斜直线；F E两截面上的弯矩之差等于EF段的Fs图面积，即2乘2=4 KNm，所以F截面的弯矩值为4+4=8 KNm。在FB段，由于q(x)=-2 KN/m,所以M图为一上凸的曲线。又因为在F截面以右的1m的FG段内Fs0，所以FG段得M图为一上增而向上凸的二次曲线。在G截面，因为Fs由正变负，所以，M图在这里有极大值。G F两截面上的弯矩之差等于FG段的Fs图的面积，即1/2乘2乘1=1 KNm所以，G截面的弯矩值为8+1=9 KNm。B G两截面上的弯矩之差等于GB段的Fs图的面积 即1/2乘（-6）乘3=-9 KNm所以，B截面的弯矩值为9-9=0，M图回到0.全梁的M图绘于（

10、c）中。如图四、求解将梁划分为20个单元，21个节点，用BEAM3来建立单元进行静力学分析。交互式的求解过程1 创建节点1.1 创建梁的各个节点1 Main Menu：PreprocessorModelingCreateNodeIn Active CS。2 在创建节点窗口内，在NODE后的编辑框内输入节点号1，并在X，Y，Z后的编辑框内输入0，0，0作为节点1的坐标值。3 按下该窗口内的Apply按钮。4 输入节点号21，并在X，Y，Z后的编辑框内输入10，0，0作为节点21的坐标值。5 按下OK按钮。6 Main Menu：Preprocessor-Modeling-CreateNodeFi

11、ll between Nds。7 在图形窗口内，用鼠标选择节点1和21。8 按下Fill between Nds窗口内的Apply按钮。9按下OK按钮，完成在节点1到节点21之间节点的充1.2 显示各个节点1 Utility Menu：PlotctrlsNumberings2 将Node numbers项设置为On。3 Utility Menu：PlotNodes4 Utility Menu：ListNodes5 对出现的窗口不做任何操作，按下OK按钮。6 浏览节点信息后，关闭该信息窗口。2定义单元类型和材料特性2.1 定义单元类型1 Main Menu：PreprocessorElement

12、 TypeAdd/Edit/Delete2 按下Element Type窗口内的Add按钮。3 在单元类型库中，选择左侧列表中的BEAM单元家族，及右侧列表中2D elastic 3类型。4 按下OK按钮完成选择。5 按下Close按钮关闭Element Type窗口。2.2 定义材料特性1 Main Menu：PreprocessorMaterial PropsMaterial Models。1. 在材料定义窗口内选择：StructuralLinearElasticIsotropic。2 在EX后的文本框内输入数值207e5作为弹性模量。3 按下OK按钮完成定义。2.3 定义几何参数1 Ma

13、in Menu：PreprocessorReal ConstantsAdd/Edit/Delete。2 按下Real Constants窗口内的Add按钮。按下Real Constants for Element Type窗口内的OK按钮。3 依次输入1，1，0.02088，0.5。4 按下OK按钮完成定义。5 按下Real Constants窗口内的Close按钮。3创建单元3.1 创建单元1 Main Menu：PreprocessorCreateElementsAuto-NumberedThru Nodes。2 在图形窗口内，用鼠标点选节点1和2。3 按下按下OK按钮完成单元1的定义。4

14、 Main Menu：PreprocessorModelCopyElementsAuto-Numbered。用光标选择单元1，然后点Apply。5 在ITIME后的编辑框内输入20（包括被复制的单元1）作为要复制的单元总数。6 按下按下OK按钮完成单元2到单元20的定义。3.2 显示单元资料1 Utility Menu：PlotCtrlsNumberings2 在第一个下拉列表中，选择Elements numbers选项。3 Utility Menu：PlotElements4 Utility Menu：ListElementsNodes+Attributes5 浏览单元信息后，关闭该窗口。4

15、施加约束和载荷4.1 节点自由度约束1 Main Menu：SolutionDefine Loads ApplyStructural Displacement On nodes。2 用鼠标在图形窗口内选择节点5。3 按下选择窗口内的Apply按钮。4 选择自由度UY，并在VALUE后为其输入数值0。5 按下Apply按钮。6 用鼠标在图形窗口内选择节点21。7 按下选择窗口内的Apply按钮。8 选择自由度UX和UY，并在VALUE后为其输入数值0。9 按下OK按钮。4.2 施加载荷4.2.1施加节点1处的集中载荷F1。1 Main Menu：SolutionDefine Loads Appl

16、yStructuralForce/Moment On nodes。2 用鼠标在图形窗口内选择节点1。3 按下选择窗口内的Apply按钮。4 在第一个下拉列表中选择FY，并在下面的文本框内输入其值-3（向上为Y轴正方向）。5 按下Apply按钮。4.2.2施加节点7处的集中载荷F2。6 Main Menu：SolutionDefine Loads ApplyStructuralForce/Moment On nodes。7 用鼠标在图形窗口内选择节点7。8 按下选择窗口内的Apply按钮。9 在第一个下拉列表中选择FY，并在下面的文本框内输入其值-1（向上为Y轴正方向）。10 按下Apply按钮

17、。4.2.3施加节点3处的弯矩m。1 Main Menu：SolutionDefine Loads ApplyStructuralForce/Moment On nodes。2 用鼠标在图形窗口内选择节点9。3 按下选择窗口内的Apply按钮。4 在第一个下拉列表中选择MZ，并在下面的文本框内输入其值-5（逆时针为正方向）（对照上面第4步）。5 按下OK按钮。4.2.4施加单元13到单元20上的的分布载荷q。1 Main Menu：SolutionDefine LoadsApplyStructural Pressure On Beams。2 用鼠标在图形窗口内选择单元13到单元20。3 按下选

18、择窗口内的Apply按钮。4 在LKEY后的文本框内输入数值1。5 在VALI和VALJ后的编辑框内分别输入10，6 按下OK按钮。5求解5.1 定义分析类型1 Main Menu：Solution Anslysis Type New Analysis。2 选中Static选项。3 按下OK按钮。5.2 求解1 Main Menu：Solution SolveCurrent Ls。2 按下OK按钮关闭Solve Current Load Step窗口。3 按下Close按钮关闭求解结束后出现的Information窗口。4 浏览/STATUS Command窗口内的信息后，将其关闭。6后处理6

19、.1 显示梁变形结果1 Main Menu：General PostprocPlot ResultsContour Plot Nodal Solu. 选择DOF Solution下的Displacement vector sum2 不改变对话框内的任何项，按下OK按钮。6.2 建立单元结果表6.2.1创建单元表，计算节点弯矩。1 Main Menu：General PostprocElement TableDefine Table。2 按下Element Table Data窗口内的Add按钮。3 在Lab后的文本框内输入IMOMENT。4 在左侧列表中选择By sequence num项。5

20、 右侧列表中选择SMICS，项。6 在右侧列表下的文本框内输入SMICS，6。7.按下Apply按钮。在Lab后的文本框内输入JMOMENT。 8. 在Lab后的文本框内输入JMOMENT。9 重复上面的步骤4和5。10右侧列表下的文本框内输入SMICS，12。11按下OK按钮。6.2.2创建单元表，计算节点剪力。1 Main Menu：General PostprocElement TableDefine Table。2 按下Element Table Data窗口内的Add按钮。3 在Lab后的文本框内输入ISHEAR。4 在左侧列表中选择By sequence num项。5 右侧列表中选

21、择SMICS，项。6 右侧列表下的文本框内输入SMICS，2。7 按下Apply按钮。8 在Lab后的文本框内输入JSHEAR。9 重复上面的步骤4和5。10右侧列表下的文本框内输入SMICS，8。11按下OK按钮。6.3 列出所有表格资料6.3.1列出资料1 Main Menu：General PostprocList ResultsElement Table Data。2 在List Element Table Data窗口内选择IMOMENT，JMOMENT，ISHEAR和JSHEAR。3 按下OK按钮并在浏览资料窗口内的信息后，将其关闭。6.3.2画剪力图1 Main Menu：Gen

22、eral PostprocPlot ResultsLine Elem Res2 在第一个下拉列表中选择ISHEAR，在第二个下拉列表中选择JSHEAR。3 按下OK按钮。6.3.3画弯矩图1 Main Menu：General PostprocPlot ResultsLine Elem Res2 在第一个下拉列表中选择IMOMENT，在第二个下拉列表中选择JMOMENT。3 按下OK按钮。弯矩图及其相应数据如下（图形已经过反色处理）：7退出程序1 Toolbar：Quit。2 选择Quit-No Save！3 按下OK按钮。五、结果分析 梁的弯矩在有集中力偶的地方会发生跳变，而剪力是在有集中力

23、的地方会有跳变的两种方法的求解结果一样，证明在运用正确的方法，选用正确的单元与节点进行有限元的分析，能得到与实际相符的结果，所以在。工程实际中将实际问题转化成数学与物理模型，能得到想要的结果第二章机械优化设计说明一题目及解析用黄金分割法求f(x)=2X3+3X2-6X+1的最优解.设初始点0=0,初始步长h=1,取迭代精度=0.01解:a1=a0=0 , f1=f(a1)=1 a2=a1+h=1 , f2=f(a2)=0 f1f2 , 作前进运算 h=2h=2 , f3=f(a1)=17 出现了“两大头,中间小”的情况 初始搜索区间为a,b=0,2 a1=a+0.382(b-a)=0.764

24、, f1=-2.941025a2=a+0.618(b-a)=1.236 , f2=1.9435525 f1 继续比较f1 继续比较 省略求解步骤 得结果 a=0.5(a+b)=0.617934f=f(a)=-1.090262二黄金分割法计算框图三C语言程序#include#include#include#define e 0.0001#define tt 1float function(float x)float y=2*pow(x,3)+3*pow(x,2)-6*x+1; return(y);void finding(float a3,float f3) float t=tt,a1,f1,i

25、a;int i; a0=0; f0=function(a0); for( i=0;i+) a1=a0+t;f1=function(a1); if(f1=e) t=-t;a0=a1;f0=f1; else t=t/2;for(i=0;i+)a2=a1+t;f2=function(a2); if(f2f1) break;t=2*t;a0=a1;f0=f1;a1=a2;f1=f2;if(a0a2)a1=a0;f1=f0;a0=a2;f0=f2;a2=a1;f2=f1;return;float gold(float *ff)int i;float a13,f13,a4,f4;float aa;find

26、ing (a1,f1);a0=a10;f0=f10;a3=a12;f3=f12;a1=a0+0.382*(a3-a0);a2=a0+0.618*(a3-a0);f1=function(a1);f2=function(a2);for( i=0;i+)if(f1=f2)a0=a1;f0=f1;a1=a2;f1=f2;a2=a0+0.618*(a3-a0);f2=function(a2);else a3=a2;f3=f2;a2=a1;f2=f1;a1=a0+0.382*(a3-a0);f1=function(a1);if(a3-a0)e)aa=(a1+a2)/2;*ff=function(aa);break;return(aa);void main()float xx,ff;xx=gold(&ff);printf(nThe Optimal Design Result Is:n);printf(ntx*=%fntf*=%f,xx,ff);getch();四运行结果五 结果分析程序运行结果与实际计算结