《现代机械设计方法课程设计.doc》由会员分享，可在线阅读，更多相关《现代机械设计方法课程设计.doc（22页珍藏版）》请在三一办公上搜索。
1、目录绪论 3一. 有限元方法和优化方法的选择4二. 外伸梁梁受力分析42.1 工程问题42.2 解析法求解42.3 有限元计算说明62.4 结果分析182.4.1 有限元结果182.4.2 分析方法结果182.4.3 结果比较与结论18三. 黄金分割法19 3.1 黄金分割法简介19 3.2 黄金分割法计算框图19 3.3 问题与结果19附录20参考文献22绪论ANSYS软件是一个功能强大而灵活的设计分析软件包。该软件可运行于从PC机，NT工作站，UNIX工作站直至巨型机的各类计算机及操作系统中，数据文件在其他所有的产品系列和工作平台上兼容。ANSYN在PC机上生成的模型同样可运行在巨型机上，
2、这样就保证了所有的ANSYS用户的多领域，多变工程的问题。有限元法是依赖于计算机软硬件技术的一种数字模拟技术。由于这种方法的有效性，其理论与应用发展非常迅速。有限元法的应用已从固体力学发展到流体力学、热力学、电磁学、声学、光学、生物学等多项耦合场问题。有限单元法在某种程度上可以取代实验，从而加快了设计或者研究问题的进度，大大提高了其安全性、可靠性和经济性，为工业技术的进步起到了巨大的推动作用，常用有限元分析软件有ANSYS等。有限元法原理是将连续的求解域离散为一组单元的组合体，用在每个单元内假设的近似函数来分片的表示求解域上待求的未知场函数，近似函数通常由未知场函数及其导数在单元各节点的数值插
3、值函数来表达。从而使一个连续的无限自由度问题变成离散的有限自由度问题。有限元分析计算的基本思想：物体离散化，单元特性分析，选择位移模式，分析单元的力学性质，计算等效节点力，单元组集，求解未知节点位移。用有限元分析计算所获得的结果只是近似的。如果划分单元数目非常多而又合理，则所获得的结果就与实际情况相符合。有限元分析法的应用使设计水平发生了质的飞跃，在机械工程领域主要表现在以下几个方面： 1. 增加设计功能，缩短设计和分析的循环周期，减少设计成本，增加产品和工程的可靠性； 2. 采用优化设计，降低材料的消耗或成本； 3.在产品制造或工程施工前预先发现潜在的问题，模拟各种试验方案，减少试验时间和经
4、费，； 4. 进行机械事故分析，查找事故原因。这次课程设计使用有限元法对外伸梁进行分析，与解析法做对比；用黄金分割法解决最优解问题。 一.有限元方法和优化方法的选择有限元问题：外伸梁利用ANSYS软件进行分析。优化问题：利用黄金分割法解决最优解问题。二. 外伸梁受力分析2.1 工程问题有一外伸梁结构如下图所示，外伸梁上均布载荷的集度为q=3KN/m，集中力偶M=3KN/m，对该梁进行分析，画出弯矩图和剪力图（高等教育出版社刘鸿文编材料力学I第121页例4.4） 已知参数： 材料特性：弹性模量 E=2.07Gpa边界条件：此问题的边界条件为位移边界条件。 A处的约束条件是固定铰支，限制了x方向和
5、y方向两个自由度，即 UX=0,UY=0。 B处的约束条件是滑动铰支，限制了y方向一个自由度，即UY=0。2.2解析法求解由静力平衡方程 M=0 6xFM+qx6x5 = 0 M=0 Fx6 Mqx4x2qx2x1 = 0求得支反力为 F= 14.5 KNF= 3.5 KN以梁的左端为坐标原点，选取坐标系如图所示。梁在CA,AD,DB三段内，剪力和弯矩不能用同一方程式来表示，应分段考虑。在CA段内，Fs（x）= qx = 3x (0x2m) (a) M（x）= 1/2qxx = 3/2xx (0x2m) (b)在AD段内， Fs（x）= Fa qx = 14.5 3x (2mx6m) (c)
6、M（x）= Fa （x2）1/2qxx = 14.5(x2) 3/2xx ( 2x 6m) (d)M(x)是x的二次函数，根据极值条件可得， 14.53x=0, 解得x=4.82m,亦在这一截面上，弯矩为极值。代入（d）式得AD段内的最大弯矩为 M=6.04KNm。在DB段内, Fs（x）= F=3.5KN （6mx 8m） M（x）= F（8x）=3.5(8x) (6m x 8m) 依照剪力方程和弯矩方程，分段作剪力图和弯矩图。从图中看出，沿梁的全部长度，最大剪力为Fsmax=8.5KN,最大弯矩为Mmax=7KNm。2.3有限元计算说明将梁划分为16个单元，17个节点，用BEAM3来建立单
7、元，进行静力学分析。交互式的求解过程1.创建节点（1）创建梁的各个节点l Main Menu：PreprocessorModelingCreateNodeIn Active CS。l 在创建节点窗口内，在NODE后的编辑框内输入节点号1，并在X，Y，Z后的编辑框内输入0，0，0作为节点1的坐标值。l 按下该窗口内的Apply按钮。l 输入节点号17，并在X，Y，Z后的编辑框内输入8，0，0作为节点17的坐标值。l 按下OK按钮。l Main Menu：Preprocessor-Modeling-CreateNodeFill between Nds。l 在图形窗口内，用鼠标选择节点1和17。l
8、按下Fill between Nds窗口内的Apply按钮。l 按下OK按钮，完成在节点1到节点17之间节点的填充。（2）显示各个节点l Utility Menu：PlotctrlsNumberingsl 将Node numbers项设置为On。l Utility Menu：PlotNodesl Utility Menu：ListNodesl 对出现的窗口不做任何操作，按下OK按钮。l 浏览节点信息后，关闭该信息窗口.2定义单元类型和材料特性（1）定义单元类型l Main Menu：PreprocessorElement TypeAdd/Edit/Deletel 按下Element Type窗
9、口内的Add按钮。l 在单元类型库中，选择左侧列表中的BEAM单元家族，及右侧列表中2D elastic 3类型。l 按下OK按钮完成选择。l 按下Close按钮关闭Element Type窗口。（2）定义材料特性l Main Menu：PreprocessorMaterial PropsMaterial Models。l 在材料定义窗口内选择：StructuralLinearElasticIsotropic。l 在EX后的文本框内输入数值207e5作为弹性模量。l 按下OK按钮完成定义。（3）定义几何参数l Main Menu：PreprocessorReal ConstantsAdd/Ed
10、it/Delete。l 按下Real Constants for Element Type窗口内的OK按钮。l 依次输入1，1，0.02088，0.5。l 按下OK按钮完成定义。l 按下Real Constants窗口内的Close按钮。3创建单元（1）创建单元l Main Menu：PreprocessorCreateElementsAuto-NumberedThru Nodes。l 在图形窗口内，用鼠标点选节点1和2。l 按下按下OK按钮完成单元1的定义。l Main Menu：PreprocessorModelCopyElementsAuto-Numbered。用光标选择单元1，然后点A
11、pply。l 在ITIME后的编辑框内输入16（包括被复制的单元1）作为要复制的单元总数。l 按下按下OK按钮完成单元2到单元16的定义。（2）显示单元资料l Utility Menu：PlotCtrlsNumberingsl 在第一个下拉列表中，选择Elements numbers选项。l Utility Menu：PlotElementsl Utility Menu：ListElementsNodes+Attributesl 浏览单元信息后，关闭该窗口。4施加约束和载荷（1）节点自由度约束l Main Menu：SolutionDefine Loads ApplyStructural Di
12、splacement On nodes。l 用鼠标在图形窗口内选择节点5。l 按下选择窗口内的Apply按钮。l 选择自由度UX和UY，并在VALUE后为其输入数值0。l 按下Apply按钮。l 用鼠标在图形窗口内选择节点17。l 按下选择窗口内的Apply按钮。l 选择自由度UY，并在VALUE后为其输入数值0。l 按下OK按钮。（2）施加载荷施加节点13处的弯矩ml Main Menu：SolutionDefine Loads ApplyStructuralForce/Moment On nodes。l 用鼠标在图形窗口内选择节点13。l 按下选择窗口内的Apply按钮。l 在第一个下拉列
13、表中选择MZ，并在下面的文本框内输入其值-3（逆时针为正方向）（对照上面第4步）。l 按下OK按钮。施加单元1到单元12上的的分布载荷q。l Main Menu：SolutionDefine LoadsApplyStructural Pressure On Beams。l 用鼠标在图形窗口内选择单元1到单元12。l 按下选择窗口内的Apply按钮。l 在LKEY后的文本框内输入数值1。l 在VALI和VALJ后的编辑框内分别输入-3，l 按下OK按钮。5求解（1）定义分析类型l Main Menu：Solution Anslysis Type New Analysis。l 选中Static选项
14、。l 按下OK按钮。（2）求解l Main Menu：Solution SolveCurrent Ls。l 按下OK按钮关闭Solve Current Load Step窗口。l 按下Close按钮关闭求解结束后出现的Information窗口。l 浏览/STATUS Command窗口内的信息后，将其关闭6后处理（1）显示梁变形结果l Main Menu：General PostprocPlot ResultsContour Plot Nodal Solu. 选择DOF Solution下的Displacement vector suml 不改变对话框内的任何项，按下OK按钮。（2）建立单元
15、结果表创建单元表，计算节点弯矩。l Main Menu：General PostprocElement TableDefine Table。l 按下Element Table Data窗口内的Add按钮。l 在Lab后的文本框内输入IMOMENT。l 在左侧列表中选择By sequence num项。l 右侧列表中选择SMICS，项。l 在右侧列表下的文本框内输入SMICS，6。l 按下Apply按钮。l 在Lab后的文本框内输入JMOMENT。l 重复上面的步骤4和5。l 右侧列表下的文本框内输入SMICS，12。l 按下OK按钮。创建单元表，计算节点剪力。l Main Menu：Gener
16、al PostprocElement TableDefine Table。l 按下Element Table Data窗口内的Add按钮。l 在Lab后的文本框内输入ISHEAR。l 在左侧列表中选择By sequence num项。l 右侧列表中选择SMICS，项。l 右侧列表下的文本框内输入SMICS，2。l 按下Apply按钮。l 在Lab后的文本框内输入JSHEAR。l 重复上面的步骤4和5。l 右侧列表下的文本框内输入SMICS，8。l 按下OK按钮。（3）列出所有表格资料列出资料l Main Menu：General PostprocList ResultsElement Tabl
17、e Data。l 在List Element Table Data窗口内选择IMOMENT，JMOMENT，ISHEAR和JSHEAR。l 按下OK按钮并在浏览资料窗口内的信息后，将其关闭。（4）画剪力图和弯矩图l 在第一个下拉列表中选择ISHEAR，在第二个下拉列表中选择JSHEAR。l 按下OK按钮。l Main Menu：General PostprocPlot ResultsLine Elem Resl 在第Main Menu：General PostprocPlot ResultsLine Elem Resl 一个下拉列表中选择IMOMENT，在第二个下拉列表中选择JMOMENT。l
18、 按下OK按钮。剪力图和弯矩图及其相应数据如下：7退出程序l Toolbar：Quit。l 选择Quit-No Save！l 按下OK按钮。2.4 结果分析2.4.1有限元结果l 剪力结果分析 在A点有向下最大的剪力 ,数值为8.5KN. (由A到C) 在A点有向上最大的剪力 ,数值为6KN。 （由D到B） 在A点最大剪力从向上的6KN突变到向下的8.5KN. 从D点到B点剪力由向下的3.5KN不变l 弯矩结果分析 在D点有顺时针最大的弯矩 数值为7KNm 在A点有逆时针最小的弯矩8KNm 在D点弯矩由突变 从A点到D点由逆时针的6KNm变到顺时针的4KNm2.4.2分析方法结果 l 剪力结果
19、分析 从C点到A点有向下最大的剪力 数值为6KN 从D点到A点有向上最大的剪力 数值为8.5KN 从D点到B点剪力由向下的3.5KN不变 在A点剪力由向下的6KN突变到向上的8.5KNl 弯矩结果分析 在D点有顺时针最大的弯矩 数值为7KNm 在A点由逆时针最小弯矩6KNm 在D点弯矩由突变 从A点到D点由逆时针的6KNm到顺时针的4KNm2.4. 3结果比较与结论 按材料力学原理计算的结果与用Ansys软件运行出的结果进行比较，两者计算结果一样。三. 黄金分割法3.1 黄金分割法简介在实际工程优化设计中，目前应用最多的一维搜索方法是黄金分割法，又称作0.618法，黄金分割法属于区间消去法的范
20、围。3.2 黄金分割法计算框图3.3 问题与结果问题：已知一维目标函数f（）=*7*+10，利用黄金分割法求解极小点，精度e=0.35,已知初始区间2,8。解用黄金分割法求解1=a+0.382(b-a)=4.292, f1= 1.6227362= a+0.618(b-a)=5.708, f2=2.625264比较函数值，有f1 f2, 则b=2=5.708，2=1=4.292, f2= f1= 1.622736, 1=a+0.382(b-a)=3.416456f1= f（1）= 2.24302判断迭代终止条件：b-ae不满足迭代终止条件，比较函数值f2 ,f1, 继续缩短区间。直到满足迭代条件
21、为止。其中迭代了7次，可满足给定精度，则最优解为*=0.5(b+a)=3.583685 f*=f(*)=2.242997附录（源程序）#include stdafx.h#include#includevoid main()float a=2,b=8,e=0.35,y1,y2,r1=0.382,r2=0.618,a1,a2,u,w;int n=0;a1=a+r1*(b-a);a2=a+r2*(b-a);y1=a1*a1-7*a1+10;y2=a2*a2-7*a2+10;printf(第%d次迭代：na=%f,b=%f,y1=%f,y2=%fn,n,a,b,y1,y2);doif(y1=y2)a=
22、a1;a1=a2;y1=y2;a2=a+r2*(b-a);y2=a2*a2-7*a2+10;n+;printf(第%d次迭代：na=%f,b=%f,y1=%f,y2=%fn,n,a,b,y1,y2);elseb=a2;a2=a1;y2=y1;a1=a+r1*(b-a);y1=a1*a1-7*a1+10;n+;printf(第%d次迭代：na=%f,b=%f,y1=%f,y2=%fn,n,a,b,y1,y2);while(b-a)e);u=(b+a)/2;w=u*u-7*u+10;printf(u=%f,w=%fn,u,w);运行结果:参考文献1 刘鸿文主编，材料力学(第四版)；高等教育出版社2 谭浩强主编，c程序设计；清华大学出版社3 倪洪启、谷耀新主编，现代机械设计方法；化学工业出版社
