《现代设计方法课程设计有限元分析课程设计》

1、现代设计方法课程设计课题名称有限元分析课程设计系 别机械工程系专 业 机械制造及自动化班 级 08 机制四（3）班姓名学号指导老师完成日期指导教师评语建议成绩：优 良 中 及格 不及格指导教师签字年 月 日最终评定成绩：优 良 中 及格 不及格系主任签字年 月 日目录1 . 前言 42 .ANSYSa件介绍442.2.ANSYS软件主要功能563 . 有限元结构分析 77794 .机械优化设计910111214145 . 课程设计感言 14一前言有限元方法是力学、计算数学和现代计算技术相结合的产物，是一种求解微分方程边界值问题和初值问题的强有力的数值方法， 是求解各种复杂数学物理问题的重要方法

2、， 同时也是处理各种复杂工程问题的重要手段也是进行科学研究的重要工具有限单元法自 20 世纪 50 年代诞生以来， 其理论体系不断完善， 新的应用软件程序不断涌现， 现已成为工程分析的必备工具。 我国广大科技工作者经过数十年的不懈努力， 在理论研究和工程应用两方面都取得了丰硕的成果和宝贵的实践经验，解决了许多重大工程的问题，为国民经济建设做出了贡献。随着计算机的广泛应用， 结合矩阵分析方法， 有限元法不到半个世纪的发展，几乎渗透到了科学、 工程的方方面面从最初的结构力学位移法到固体的弹性力学分析、 固体的弹塑性、 粘弹性、 粘塑性分析、 流体力学和流固耦合作用分析；从静力分析到静动力分析、 稳

3、定性分析、 优化设计； 从单纯的力学分析到物理分析、材料设计、机械制造等多学科。总之，有限元法的应用范畴从力学领域扩展到了物理、 材料、 机械等多种工程领域， 并且在工程及科学研究中具有重要意义。本课程设计选取了受均布载荷作用和集中力偶力矩的外伸梁问题，用ANSY驮件进行应力和结构安全分析。 所得结果于弹性力学相比较， 二者的误差在允许误差 范围之内，结论令人满意。二.ANSY歆件介绍2.1 ANSYS 简介：ANSY歆件是融结构、流体、电场、磁场、声场分析于一体的大型通用有限 元分析软件。由世界上最大的有限元分析软件公司之一的美国ANSYSF发。它能与多数CAD软件接口，实现数据的共享和交换

4、，如 Pro/Engineer,NASTRAN, Alogor, I -DEAS, AutoCAD等， 是现代产品设计中的高级CAD工具之一。软件主要包括三个部分：前处理模块，分析计算模块和后处理模块。前处理模块提供了一个强大的实体建模及网格划分工具， 用户可以方便地构造有限元模型；分析计算模块包括结构分析 （可进行线性分析、 非线性分析和高度非线性分析） 、流体动力学分析、电磁场分析、声场分析、压电分析以及多物理场的耦合分析，可模拟多种物理介质的相互作用，具有灵敏度分析及优化分析能力；后处理模块可将计算结果以彩色等值线显示、 梯度显示、 矢量显示、 粒子流迹显示、立体切片显示、透明及半透明显

5、示（可看到结构内部）等图形方式显示出来，也可将计算结果以图表、曲线形式显示或输出。2.2 ANSYS软件主要功能软件提供了 100 种以上的单元类型， 用来模拟工程中的各种结构和材料。 该软件有多种不同版本，可以运行在从个人机到大型机的多种计算机设备上，如PC，SGI, HP； SUN DEC IBM, CRAW。ANSYSg布10.0新版本美国宾夕法尼亚州 6 月 2 日延续了 ANSYS-贯强大的耦合场技术，10.0版本为复杂的流周耦合（FSI） 问题提供了更完善的解决方案。 该版本整合了世界一流的应力分析和流体分析技 术， 形成了一套完整的 FSI 解决方案。 通过适合于特定场要求的网格

6、划分， 一个 单一的几何体可以应用于两种场。 该版本提供了有效地解决FSI 动力学分析的信息交换功能。 目前市场上没有任何其他的 FSI 软件可以提供如此强大的稳健性和 高度的精确性分析。另外，该版本可以在多个机群进行并行处理解决超大模型。"ANSYS 10.0代表了最先进的CAE整合技术，较9.0有了显著的提高”， ANSY宓司总裁兼首席执行官Jim Cashman说，"我们一直致力于拓展 ANSY助 真技术的广度和深度， 同时建立各种类型的仿真分析软件的空前大连盟。 得益于 ANSYSWorkbench整合CAEg术的架构，我们创建了建模、仿真、分析、前后处 理的一系列

7、无缝链接。 10.0 新版本整合了世界上最优秀的结构、热、流体等分 析功能。”ANSYS10.6口入了旋转机械和叶片设计工具，丰富了 Workbench环境下的行 业化功能。即 ANSYS BladeModeler, 一款针对旋转机械叶片构件的高效的三维 设计工具； 以及 ANSYSTurboGrid ， 一款高质量的叶片设计六面体网格划分工具。”结合了 ANSYSCFX和涡轮专用白前后处理 CFM能，10.0版本提供了涡轮 机械设计和分析完整的解决方案，"ANSY%司副总裁兼总经理Chris Reid说，" 应力分析、 计算流体动力学分析或流固耦合分析的模型可以直接建立，

8、 通过 CAD 系统连通性，可以把模型扩展到上下游部件，最终完成整个模型的分析。 ANSYS Workbench 是提供此功能上独一无二的环境，借此空气动力学工程师可以进行 CFDS计，同时确认结构特征。这将大幅度缩短设计流程。”在机械应用领域，ANSYS 10.0包括了 ANSYS Workbench全部的热瞬态分 析功能。这不仅帮助用户进行非常复杂的时域仿真，同时ANSYSVorkbench也可自动完成很多建模和求解工作。 这样可以轻松快速地求解设备在一定运行时间内 的热性能。为了满足日益增加的对大型复杂问题及时有效的分析需求， ANSYS 10.0 的 并行求解器如今可增加了对 CPUf

9、fi通信技术的选择余地。除了支持 Ethernet和 Gigabit Ethernet ， ANSYS 10.0还支持Myrinet 和 InfiniBand 。相对于以前的架构， ANSYS 10.0能以最少的成本满足高性能的机群计算。本着以低成本硬件设备提供高性能解决方案的目标，ANSYSWorkbench现可支持 Windows XP 64位机的AM邯口 EMT64K片集。此项改革解决了许多用户在 Windows操作系统下运行大型模型所面临的 2GB内存限制。另外，它也使得ANSYS 用户不再需要写硬盘就能完成整个求解，从而节约求解时间。对于用户， 这将帮助他们更加经济有效地解决大型模型

10、问题， 如那些低频稳态和全瞬态电磁分析问题。 ANSYS 10.0 并行求解器可以解决高于一亿自由度的大型电磁问题，在CAE行业独树一帜。在高频电磁领域， 10.0 版本提供了一个新的模式端口。此端口大大简化了集成电路（IC）、射频识别（RFID）和射频微机电系统（MEMS殍多种设备分析传输 线端口的建模。标准算例显示，利用此端口建模，可以显著缩小模型尺寸，在保证精确的频域计算结果前提下，节约30到50的求解时间和内存需求。新版本增加了旋转机械的陀螺效应，它提高了 ANSYS寸涡轮机械和其他旋转结构的转子动力学分析的能力。 在耦合场领域，结构热电磁三场耦合分析中增加热弹阻尼（TED）, 一个在

11、金属、制陶及 MEM颔域非常重要的内耗装置。ANSYS昧续Workbench主旋律，提供我们的用户可供选择的全自动或个人控制的强大分析软件。我们在核心的网格处理技术上有十足的增强，在ANSYSWorkbench各个应用程序间共享网格。另外，双向参数互动的CADg口的稳健性也得到了提高。ANSYS® ICEM CFD? 10.0通过混合网格剖分新功能和 CAD奠 型细节处理功能， 提供了完整的一系列网格划分工具以模拟真实世界， 如汽车引 擎罩下的散热分析和汽车碰撞分析。2.3 ANSYS 软件主要特点主要技术特点：? 唯一能实现多场及多场耦合分析的软件?唯一实现前后处理、求解及多

12、场分析统一数据库的一体化大型FEA软件?唯一具有多物理场优化功能的FEA软件? 唯一具有中文界面的大型通用有限元软件? 强大的非线性分析功能? 多种求解器分别适用于不同的问题及不同的硬件配置? 支持异种、异构平台的网络浮动，在异种、异构平台上用户界面统一、数 据文件全部兼容? 强大的并行计算功能支持分布式并行及共享内存式并行多种自动网格划分技术，良好的用户开发环境三有限元结构分析m 80kN mq 20kN /m ，F 40kN 。对该梁进行分析，画出弯矩图和剪力图解：由梁的平衡方程，求出 A点和B点的支座反力：Ma 0 F*5 FA*4 M+20*4*2=0 FA=70KNMB 0F*1 M

13、- FB*4 20*4*2=0FB=50KNCA 段内：Ms(x)40*xAB 段内：Ms(x): 40x+801/2q(x-1)基本数据：该梁选择45号优质碳素结构钢，进行调质处理，其机械性能查机械设计表 8-1 得：B 637MPa , S 353MPa ,1 268MPa ,i 155MPa,0.2,0.1。取此钢的横截面积为矩形，h=0.5m其横截面积为1,惯性矩为0.020833。将梁划分为10个单元，11个节点，用BEAM来建立单元，进行静力学分析。 在节点3到11处施加分布载荷，q=20KN/M在节点1处施加一竖直方向力ANUL 7 2OLD21 QIF= 40KN在节点3处施加

14、一集中力偶80KN*M依照上面所施加的各个约束及载荷，用ANSYS乍出如下所示受力图：ETEHEtnS exbm mmFH/STATUS CouandFilsSOLUTION OPT 1 OHSFROBLEH DIMENSIONALITY2-DDEGREES OF FREEDOMUX UV ROTZRNRLVSIS I?Pb.SIRItC (STfcflUV-SGLOBALLY ASCENDLED M由TRIMGVMKETR1CLOAD STEP 0 T I I 0 N SLOAD STEP NUMEER 1TIME AT EHD OF THE LOAD STEP. . .一 一一 一一一 1

15、.0Q0QNUMBER OF SUBSTEFS * .1STEP CHANGE B0UNIX1BV CONPITIOHG HOPRINT OUTPWr CONTROLSNO PRINTOUTDATADfiGE OUTPUT 8HTRQLU11LL DfiTfl UFtITTCN求解观察梁的变形图变形之后如下图所示HrKKTAU CwetiHFil*PHIMT ELEHENT TABLE ITEMS FER ELEMENT.一. FOSTA ELEMENT TABLE LISTING *STATCURRENTCURRENTCURRENTCURRENTkL&RJNUNHNiJbHEjAR1

16、0 .000B-2B.0B04B _ B0I340.0QB220=30040 a 00040 * 00040,00034WHMM&2.GBHTH . HHH-HW.HMM4.巧况阿_/词。20,000id. 5。阿560b00B&2 .&ER0TL回口网0b74607E-136O.SWRriR.RR«.4Kfi33E-131 6 H HR1R760h00052&50f»10,00020h000H5N.5。修立修 HMDZH.H（州946一函m22.500中。.0酬046一忖0。1022bS000.189eSE-1340.00050b000MI

17、NIHUMUALUESEL正M22J3ANSYS 软 件LINK STU33STBP-lSUB -1.TIME-1工 SHE第R JSHEACeHita =-3Q解工0-3wra 7 口 fJJLElT”所求剪力图与弯矩图如下3.3 .结果分析：根据以上分析得出如下结果：1根据结果，变形图符合工程规律。2误差小于3%符合要求。3结构强度符合设计要求按材料力学原理计算的结果与用 Ansys软件运行出的结果进行比较，可以看出两 者有一定的偏差，但是偏差不大。再将元格分开对比 ，两者在各个受力点处值几 乎一样。四.机械优化设计用黄金分割法求f(x)=2X 3+3X-6X+1的最优解.设初始点Q

18、76;=0，初始步长h=1,取 迭代精度e =0.01解:a1=a0=0 , f1=f(a1)=10a2=a1+h=1 , f2=f(a2)=0,.fl>f2 ,作前进运算h=2h=2 , f3=f(a1)=15 出现了 “两大头，中间小”的情况 初始搜索区间为a,b=0,2a1=a+0.382(b-a)=0.764 , f1=-2.941025a2=a+0.618(b-a)=1.236 , f2=1.9435525vf1<f2 b=a2=1.236a2=a1=0.764 ,f2=f1=-2.941025a1=a+0.382(b-a)=0.472 , f1=-0.950570判断迭

19、代终止条件：.b-a=1.23> £ .继续比较vf1<f2 a=a2=0.472 , a1=a2=0.746 , f1=f2=-2.941025 a2=a+0.618(b-a)=0.944 f2=f(a2)=-0.307366 判断迭代终止条件：v b-a=1.236-0.472=0.764> & 继续比较省略求解步骤 得结果a=0.5(a+b)=0.617934 f、f(a 尸-1.090262语言程序#include<stdio.h>#include<conio.h>#include<math.h>#define e

20、 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,ia;int i;a0=0;f0=function(a0);for( i=0;i+)a1=a0+t;f1=function(a1);if(f1<f0) break;if(fabs(f1-f0)>=e) t=-t;a0=a1;f0=f1; else t=t/2;for(i=0;i+)a2=a1+t;f2=func

21、tion(a2);if(f2>f1) break;t=2*t;a0=a1;f0=f1;a1=a2;f1=f2;if(a0>a2)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;finding (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(f