《简支梁有限元分析.doc》由会员分享，可在线阅读，更多相关《简支梁有限元分析.doc（15页珍藏版）》请在金锄头文库上搜索。
1、现代设计方法课程训练任务书学 院机械工程学院专 业机械设计制造及自动化学生姓名董峰班级、学号0801012206训练报告题目平面外伸梁梁结构的内力计算以及黄金分割优化设计技术参数、内容及要求：一、 有限元课程训练1. 学习CAE软件ANSYS,主要上机练习有（1）连杆的静力学分析 （2）桁架的有限元分析 （3）梁与曲轴结构的内力计算 （4）压力容器的静力学分析 （5）机翼模型的模态分析 （6）压杆稳定临界载荷计算 （7）过盈配合与拔销耦合分析2. 由学生通过调研，在工厂、企业或科研单位进行工程实践的基础上，结合实际需要自己拟定的题目。 二、 优化课程 上机调试优化计算程序，并结合工程实际自找算
2、例进行计算。主要上机练习有：（1）二次插值法 （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二有限元
3、设计 42.1外伸梁受问题 42.2力学模型 42.3有限元模型 52.4结果分析 52.4.1 有限元结果 52.4.2 材料力学分析结果 72.4.3 结果比较与结论 8三优化设计 83.1黄金分割法简介 83.2黄金分割方框图 93.3 问题与结果 103.3.1 问题的解析过程 103.3.2 结果分析 11四附录 12五参考文献 15一、绪论有限元法在解决圣维南扭转问题近似解时首先提出的。有限元在弹性力学平面问题的第一个成功应用是由美国学者于1956年解决飞机结构强度时提出的。经过几十年得发展，有限元法已经成为现代结构分析得有效方法和主要手段。它的应用已经从弹性力学的平面问题扩展到空
4、间问题和板壳问题。对于有限元法，从选择基本未知量的角度来看，它可分为三种方法：位移法，力法，混合法。从推导方法来看，它可分为直接法，变分发，加权余数法。但随后随着计算机的发展，有限元法如虎添翼。国内外已有许多大型通用的有限元分析程序系统可供使用，在许多大型有限元分析软件都已配备了功能很强得前后处理程序，并已出现了将人工智能技术引入有限元分析软件，形成了比较完善得专家系统，逐步实现了有限元的智能化。 优化设计是现代设计方法的重要内容之一。它以数学规划为理论基础以电子计算机为工具，在充分考虑多种设计约束的前提下，寻求满足预定目标的最佳设计。优化设计理论于方法用于工程设计是在六十年代后期开始的，特别
5、是近年来，随着有限元素法，可靠性设计，计算机辅助设计的理论与发展及与优化设计方法的综合应用使整个工程设计过程逐步向自动化集成化智能化发展，其前景使令人鼓舞的。因而工程设计工作者必须适应这种发展变化，学习 掌握和应用优化设计理论与方法。 近年来随着计算机技术的普及和计算速度的不断提高，有限元分析在工程设计和分析中得到了越来越广泛的重视，已经成为解决复杂的工程分析计算问题的有效途 径，现在从汽车到航天飞机几乎所有的设计制造都已离不开有限元分析计算，其在机械制造、材料加工、航空航天、汽车、土木建筑、电子电器，国防军工，船舶， 铁道，石化，能源，科学研究等各个领域的广泛使用已使设计水平发生了质的飞跃。
6、本次课设研究简支梁的受力问题，最大剪力4.2KN,最大弯矩3.8KNm.有限元分析与材料力学分析结果相一致，误差少于3%。二 、有限元设计2.1外伸梁受力问题外伸梁的受力情况如下图所示，C点受集中力3KN，B点受支反力FB=8KN, A点受支反力FA=7.2KN。AD间受均布载荷q=2KN/m，D点受转矩6KN/m。梁总长6m，长度分布如下图。2.2力学模型2.3有限元模型（1） 根据物理模型载荷，支撑单元划分了12个单元13个节点，一个单元表示0.5m网格图如下；2.4结果分析2.4.1有限元分析1. 剪力图2. 弯矩图3. 位移变形图2.4.2材料力学分析依照题意，列出建立方程和弯矩方程，
7、方程中X以m 为单位，F(X)以KN为单位，M（X）以KNm为单位。在CA段和DB段内q=0.CA段F(X)=-3KN ，DB段F(X)=-3.8KN ，CA段M(X)=-3X M图为斜直线.DB段M(X)=22.8-3.8X M图为斜直线在AD段：F(X)=-3KN-1/2q（X-1）+7.2KN=-X+5.2KN M（X）=7.2(X-1)-3X-(X-1)(X-1) 依照剪力方程和弯矩方程，分段作剪力图和弯矩图，沿梁的全部长度，最大剪力为Fsmax=4.2kN, Mmax=3.8KNm1 剪力图和弯矩图2.4.3 结果比较与结论有限元分析与材料力学分析结果相一致。根据以上分析得出如下结果
8、：根据结果，变形图符合工程规律,误差小于3%符合要求,结构强度符合设计要求.三优化设计3.1.黄金分割法简介把一条线段分割为两部分，使其中一部分与全长之比等于另一部分与这部分之比。其比值是一个无理数，取其前三位数字的近似值0.618。由于按此比例设计的造型十分美丽，因此称为黄金分割，也称为中外比。这是一个十分有趣的数字，我们以0.618来近似，通过简单的计算就可以发现：1/0.618=1.618(1-0.618)/0.618=0.618在搜索区间【a,b】内按每次区间等比例缩短原则和对称性原则取两点1和2，符合这两个原则的计算公式为：1=a+0.382(b-a) , 2=a+0618(b-a)
9、3.2黄金分割法流程图给定:a,b,a+0.382(b-a)1,f(1)f1a+0.618(b-a)2,f(2)f2f1f2?2b,12,f1f2a+0.382(b-a) 1f(1)f11a,21,f2f1a+0.618(b-a) 2f(2)f2b-a输出:0.5(a+b)停3.3问题与结果用黄金分割法求f(x)=2X3+3X2-6X+1的最优解.设初始点0=0,初始步长h=1,取迭代精度=0.013.3.1问题的解析过程解:a1=a0=0 , f1=f(a1)=10a2=a1+h=1 , f2=f(a2)=0 f1f2 , 作前进运算 h=2h=2 , f3=f(a1)=15 出现了“两大头
10、,中间小”的情况 初始搜索区间为a,b=0,2a1=a+0.382(b-a)=0.764 , 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.0902623.3.2 结果分析程序运行结果与实际计算结果有些许误差，这是因为初选迭代步长不同而引起的,不可避免。四附录#include#include#include#define e 0.0001#define tt 1float function(float x)float y=2*pow(
11、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=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
12、=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(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=a
