《有限元受力分析--结构梁-力-计算》由会员分享，可在线阅读，更多相关《有限元受力分析--结构梁-力-计算（29页珍藏版）》请在金锄头文库上搜索。
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点稍左的截面处。由于D截面处有一向
7、下的集中力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=0；弯矩图从0开始。
8、在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图为一上增的斜直线；F
9、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图绘于（c）中。如图四、求解将
10、梁划分为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-CreateNodeFill between N
11、ds。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 TypeAdd/Edit/
12、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 Main Menu：Preproce
13、ssorReal 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 Main Menu：PreprocessorModelCopyElementsAuto-Numbered。用光标选择单元1，然后点Ap
