《现代设计方法》课程设计任务书

1、现代设计方法课程训练任务书学 院机械工程学院专 业机械设计制造及自 动化学生姓名班级、学号训练报告题目有限元和优化设计技术参数、内谷及要求：一、有限元课程训练1 .学习cae软件ansys,主要上机练习后（1）连杆的静力学分析 （2）桁架的有限元分析 （3）梁与曲轴结构的内力计 算（4）压力容器的静力学分析 （5）机翼模型的模态分析 （6）压杆稳定临界载 荷计算（7）过盈配合与拔销耦合分析2 .由学生通过调研，在工厂、企业或科研单位进行工程实践的基础上，结 合实际需要自己拟定的题目。二、优化课程 上机调试优化计算程序，并结合工程实际自找算例进行计算。 主要上机练习启：（1）二次插值法 （2）

2、powell法 （3）惩罚函数进度安排：一、有限元上机（20学时）第1次上机（4学时）：学习ansys软件，作练习1、2;第2次上机（4学时）：学习ansys软件，作练习3、4;第3次上机（4学时）：学习ansys软件，作练习5、6;第4次上机（4学时）：学习ansys软件，作练习7;第5次上机（4学时）：自拟题目上机。二、优化设计上机（12学时）优化计算程序的调试及计算算例注：利用业余时间撰写课程设计说明书。指导教师（签字）：2011年3月11日教研室主任（签字）2011年3月11日、 绪论 3二简支梁静力分析2.1 工程问题.42.2 力学模型42.3 有限元模型62.4 结果分析2.4.

3、1 有限元结果.152.4.2 分析方法结果 .152.4.3 结果比较与结论.16三优化设计1.1 进退法与黄金分割法简介 .161.2 进退法黄金分割法计算框图 .171.3 问题与结果.191.4 运行程序.201.5 结果分析.214 总结.225 参考文献.22绪论设计是人类一种复杂的思维过程，是创造性的劳动，是人类改造自然的基本活动之一。 人类在认识世界和改造世界的历史长河中， 一直在从事着设计活动。 从某种意义上说， 是人类不断进行的创新活动，推动了人类文明的进步。现代设计方法是随着当代科学技术的飞速发展和计算机技术的广泛应用而在设计领域发展起来的一门新兴的多元交叉学科， 它所涉

4、及的学科范伟非常广泛。 主要有设计方法学、 优化设计、 可靠性设计、计算机辅助设计、动态设计、有限元法、工业产品造型设计、人机工程学、并行工程、价值工程、反求工程设计、模块化设计、相似性设计、 虚拟设计、 疲劳设计、 三次设计等。 在运用它们进行工程设计时，一般都以计算机作为分析、计算、综合、决策的工具。设计人员 是新产品的重要创造者，对产品的设计和发展有重大影响。 为了适应当代科学技术发展的要求和市场经济对设计人才的需要， 必须加强设计人员的创新能力和设计素质的培养， 现代设计方法就是为达此目的而开设的。所以，今天学习现代设计方法，其目的不是要完全抛弃传统方法和经验， 而是要在掌握传统方法和

5、实践的基础上， 再掌握一些新的设计理论和技术手段， 从而推动设计方法和技术的进步，使传统设计方法产生质的变化。有限元课程设计简支梁静力分析1工程问题在下图中，外伸梁上均布载荷的集度为 q=3kn/m,集中力偶矩me=3kn.m.,列 出剪力方程和弯矩方程，并绘制剪力图和弯矩图。 （材料力学i,第四版，刘鸿 文，p121,例4.42力学模型1、梁的参数：长度l=8m ,宽度b=1m,厚度h=0.5m2 、材料参数梁选择线性、弹性、各向同性的材料。它的弹性模量e=207e5 pa.3 、梁的边界条件在节点a处受x方向、y方向的约束；节点b受y方向的约束。4 、梁的载荷ab之间作用着均布载荷 q=3

6、kn/m, d点处沿顺时针方向的集中力偶矩me=3kn.m。5、解析法求解由梁的平衡方程，求出支反力为fra=14.5kn, frb=3.5kn分ca、ad、db、三段考虑，列出剪力方程和弯矩方程，在ca段内，fs(x)= -qx = -3x(0<x< 2m)m(x)= -1/2qx 2= -3/2x2(0< x< 2m)在ad段内，fs (x) = f - qx = 14.5-3x(2m<x<6m )m(x) = fra(x-2) - 1/2qx2 = 14.5(x-2) 3/2x2(2m<x<6m )m(x)是x的二次函数，根据极值条件 dm

7、(x)/dx=0得14.5-3x=0由此解出x=4.83m,即在这一截面上，弯矩为极值。代入上式得ad段内的最大弯 矩为m = 6.04 kn - m当界面=截面取在db段内时，用截面右侧的外力计算剪力和弯矩比较方便，结果为fs(x) = - frb = - 3.5 kn(6m<x<8m)m(x) = frb (8 x) = 3.5 (8-x)(6m<x<8m)依照剪力方程和弯矩方程，分段作剪力图和弯矩图如下m曾m/kn 臼从图中看出，沿梁的全部长度，最大剪力为8.5 kn,最大弯矩为7 kn m设梁宽b=2m,高h=0.5m则截面面积a=1m2面积惯性矩 izz=b*

8、h3/12=0.020833m?3 有限元模型利用 ansys 10.0 计算说明将梁划分为 16 个单元， 17 个节点，用 beam3 来建立单元，进行静力学分析。1 创建节点1.1 创建梁的各个节点1. main menu： preprocessomodeling fcreatenodef in active cs。2 在创建节点窗口内， 在 node 后的编辑框内输入节点号1， 并在x， y， z 后的编辑框内输入 0， 0， 0 作为节点 1 的坐标值。3 按下该窗口内的apply 按钮。4 输入节点号17，并在x， y， z 后的编辑框内输入8， 0， 0 作为节点 17 的坐标值

9、。5 按下 ok 按钮。6 main menu : preprocessormodeling-create - node - fill between nds。7 在图形窗口内，用鼠标选择节点 1 和 17。8 按下 fill between nds 窗口内的 apply 按钮。9 按下 ok 按钮，完成在节点 1 到节点 17 之间节点的填充。23e7iu1 cretu kod” b，l峰 w mm”版del 置！)套 fill b*5*qnf 117h71ll fwbu frf1* fill15fflcrcat c nod-于 bttteen 2 hordenstrtmtine telw

10、t«en fill.j no史菖spacespftciinf r tt i #mme ha »e fill*i in.rlnd.ing 4r ijginal)1现 必加 bw>tr incr«»4inl -! £&r tuh 工ucf玉工i <r fill «p«r9 »d)okcwtl1.2显示各个节点1 . utility menu : plotctrls -numberings2 .将 node numbers®设置为 on。3 . utility menu： plotnodes

11、4 . utility menu： list-nodes5 .对出现的窗口不做任何操作，按下 ok按钮6 .浏览节点信息后，关闭该信息窗口。2,定义单元类型和材料特性2.1 定义单元类型1 . main menu： preprocessors element type add/edit/delete2 .按下 element type9口内的 add 按钮。3 .在单元类型库中，选择左侧列表中的 beam单元家族，及右侧列表中 2d elastic 3 类型。4 .按下ok按钮完成选择5 .按下 close按钮关闭 element type0口。2.2 定义材料特性1 . main menu：

12、 preprocessor>material prop material models。2 .在材料定义窗口内选择：structural-linear- elastic isotropic。3 .在ex后的文本框内输入数值207e5作为弹性模量。4 .按下ok按钮完成定义。2.3定义几何参数1 . main menu: preprocessors real constant”add/edit/delete。2 .按下real constants窗口内的add按钮。3 .按下 real constants forelement type9口内的 ok 按钮。4 .依次输入 1,1, 0.0

13、2088, 0.55 . 按下ok按钮完成定义。6 . 按下 real constants窗 口内的 close 按钮。3.创建单元资料3.1 创建单元1. main menu： preprocessors create- elements auto-numberedthrunodes。2 .在图形窗口内，用鼠标点选节点 1和2。3 .按下ok按钮完成单元1的定义。4 . main menu : preprocessorf model f copy f elements - auto-numbered。用光标选择单元1,然后点apply。5 .在itime后的编辑框内输入16 （包括被复制的单

14、元1）作为要复制 的单元总数。6 .按下ok按钮完成单元2到单元16的定义。3.2显示单元资料1. utility menu： plotctrls-numberings2. 在第一个下拉列表中，选择 elements numbers选项。3. utility menu： plot-elements4. utility menu： list一elementsnodes+attributes5. 浏览单元信息后，关闭该窗口。4,施加约束和载荷4.1 节点自由度约束1 . main menu ： solution - define loads - apply - structural 一 displ

15、acement- on nodes2 .用鼠标在图形窗口内选择节点 5。3 .按下选择窗口内的apply按钮。4.选择自由度ux和uy ,并在value后为其输入数值0。5.按下apply按钮。6. 用鼠标在图形窗口内选择节点17。7. 按下选择窗口内的apply按钮。8. 选择自由度uy,并在value后为其输入数值09. 按下ok按钮。4.2施加载荷4.2.1 施加节点5,17处的集中载荷f。1 . main menu ： solution - define loads - apply - structural 一force/ momentf on nodes2 .用鼠标在图形窗口内选择节

16、点53 .按下选择窗口内的apply按钮。4 .在第一个下拉列表中选择fy,并在下面的文本框内输入其值14.5,3.5 （向上为y轴正方向）。5 .按下apply按钮。4.2.2施加节点13处的弯矩m。1. main menu ： solution - define loads - apply - structural 一force/ momentf on nodes2 .用鼠标在图形窗口内选择节点 13。3 .按下选择窗口内的apply按钮。4 .在第一个下拉列表中选择mz ,并在下面的文本框内输入其值-3（逆 时针为正方向）（对照上面第4步）。4.2.3施加单元1到单元12上的的分布载荷q

17、。1 . main menu： solutionfdefine loadaapplyfstructural pressurefon beams2 .用鼠标在图形窗口内选择单元 1到单元12。3 .按下选择窗口内的apply按钮。4 .在lkey后的文本框内输入数值1。5 .在vali和valj后的编辑框内分别输入3,6.按下ok按钮。5 .求解5.1 定义分析类型1 . main menu： solution- anslysis typ门 new analysis02 .选中static选项。3 .按下ok按钮。5.2 求解1 . main menu： solution- solve-curr

18、ent ls。2 .按下 ok 按钮关闭 solve current load step®口。3 .按下close按钮关闭求解结束后出现的information窗口4 .浏览/status command®口内的信息后，将其关闭。6 .后处理6.1 显示梁变形结果1 . main menu: general postproc plot resultscontour plot nodal solu选择 dof solution 下的 displacement vector sumansysjjiu j1 1011 03502:572 .不改变对话框内的任何项，按下 ok按钮。

19、uptal sciunrmrrxi*-i sct -toe-1»sts-0工icc .41e-c4sbc -.48jt-04am-04. u 5e - 0-i3 *2e-m.414e-a4帖鸵m6.2 建立单元结果表6.2.1 创建单元表，计算节点弯矩。1 . main menu: general postprocelement tab-define table。2 .按下 element table data窗口内的 add 按钮。3 .在lab后的文本框内输入imoment 。4 .在左侧列表中选择 by sequence num项。5 .右侧列表中选择smics,项。6 .在右

20、侧列表下的文本框内输入 smics, 6。7 .按下apply按钮。8 .在lab后的文本框内输入jmoment 。9 .重复上面的步骤4和5。10.右侧列表下的文本框内输入 smics, 1211.按下ok按钮。6.2.2创建单元表，计算节点剪力。1 . main menu: general postproeelement tab-define table。2 .按下 element table data窗口内的 add 按钮。3 .在lab后的文本框内输入ishear。4 .在左侧列表中选择 by sequence num项。5 .右侧列表中选择smics,项。6 .右侧列表下的文本框内输

21、入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 results- element table data2 .在 list element table data 窗口 内选择 imoment , jmoment ,ishear 和 jshear。3 .按下ok按钮并在浏览资料窗口内的信息后，将其关闭6.3.2 画剪力图1 . main men

22、u: general postproeplot results line elem res2 .在第一个下拉列表中选择ishear ,在第二个下拉列表中选择jshear。ansysj ax 19初i工 02：oft;0?3 .按下ok按钮。lime stressstep-1sub -1tike-1i3heaf j3kem?miu 7. 5ele1u61ux $elm$-e,5-5.278-2.056l.1674.3a9一吼酹9-3.667-.4444442.77®6.3.3 画弯矩图1 . main menu: general postprocplot results- line e

23、lem res2 .在第一个下拉列表中选择imoment ,在第二个下拉列表中选择jmoment 。3 .按下ok按钮。uwe stresssub >1time*!lokment jmomemthin -6elem-6kax -7ansysjak 19 2q1102:06:48elek=14-6-3,111-p 22z2t22用e575.5567/56-1.6471.2224. lil7 .退出程序1 .toolbar： quit。2 .选择 quit-no save !3 .按下ok按钮。4结果分析8 . 1有限元结果用ansys10.0勺求解结果:最大弯矩 m ma=7kn m最小弯

24、矩 m min=-6kn - m最大剪力f ma=8.5kn最小剪力 f min=-6kn8.2 分析方法结果最大弯矩出现在剪力等于0处,最大弯矩 m ma=7kn m最小弯矩 m min=-6kn - m最大剪力 f max=8.5kn最小剪力 f min=-6kn8.3 结果比较与结论比较两种方法可以看出， 在剪力图上有所不同， 理论与实际的剪力图是关于x 轴对称， 而其大小都是相同的， 证明实用的正确的方法，理论数据与实际数据基本一致。优化课程设计进退法简介进退法的今本思想是：由单峰函数的性质可知，在极小点 a*左边函数值应严格下降，而在极小点右边函数值应严格上升。因此，可以从某一给定的

25、初始点a0出发，以初始步长h0沿着目标函数值的下降方向，逐步前进(或后退) ，直至找到相继的 3 个试点的函数值保证“两头大，中间小”为止；把“两头大，中间小”的区间作为初始区间。黄金分割法简介黄金分割法又称 0.618 法，它通过不断缩短搜索区间的长度来寻求一维函数f(a) 的极小点。这种方法的基本原理是：在搜索区间【a,bl内按每次区间等比例缩短原则和对称性原 则取两点al和a2,符合这两个原则的计算公式为:a1=a+0.382(b-a), a2=a+0.382(b-a)计算它们的函数值f1=f(a1),f2=f(a2),比较f1,f2的大小，根据单峰函数特点，极小点在“两头大，中间小”的

26、区间内。进退法计算框图黄金分割法计算框图3.3 问题与结果题目：用黄金分割法求f(x)=x2-5x+5的最优解。设置初始点为：x0=0,初始 步长h=1,取迭代精度e=0.35.首先用进退法确定搜索区间a1=a0=0, f1=f(a1)=10a2=a1+h=1, f2=f(a2)=4因f2>f3,作前进运算：h=2h=4a1=a2=1,f1=f2=4d2=d3=2,f2=f3=0a3=a2+h=4f3=f(a3)=-2因f2>f3,再作前进运算；h=2h=4,a1=a2=2,f1=f2=0a2=a3=4,f2=f3=-2a3=a2+h=6f3=f(a3)=18此时，a21,a2,a

27、3三点的函数值出现了 “两头大，中间小”的情况，所以初始搜索 区间a,b = 2,6。然后利用黄金分割法求最优解，在初始区间a,b = 2,6中去两个计算点并计算起函数值a1=a+0.382(b-a)=4.292, f1=f(a1)=-1.622736a2=a+0.618(b-a)=5.078, f1=f(a1)=2.625264比较函数值，缩短区间。因有f1<f2,则b-a2=5.708,a2=a1=4.292,f2-f1=-1.622736a1=a+0.382(b-a)=3.416456,f1=f(a1)=-2.243020判断迭代终止条件：b-a=4-2.15=1.85>0.

28、30不满足终止条件，比较函数值fl、f2,继续缩短区间。区间缩短 4次之后，满足了给定精度，迭代终止，近似最优解为a*=0.5(b+a)=2.448379, f*=f(a*)=-1.247335 。3.4 运行程序#include "stdio.h"#include "math.h"#include "conio.h"#define e 0.35#define tt 1float function(float x)float y;y=pow(x,2)-5*x+5;return(y);void searching(float a3,float f3)float h=tt,a1,f1,ia,i;a0=0;f0=function(a0);for(i=0;i+)a1=a0