有圆形的等截面销轴.docVIP

系统建模与仿真实验报告 姓名 孙阳 班级 流体11-1 学号 1107240117 一、MATLAB在优化设计中的应用 有一圆形的等截面销轴,一端固定在机架上,另一端受集中载荷P=10kN和转矩M=100N.m的作用,其简化模型如图1所示。由于结构的需要,轴的长度l??8cm,已知销轴材料的许用弯曲应力[??w]=120MPa;许用切应力[??]=80MPa;允许挠度[f]=0.01cm;密度??=7.8t/m3;弹性模量E=2??105MPa。现要求设计这根销轴,在满足使用要求的前提下,使其质量最轻[5]。 1.1数学模型的建立 优化设计追求的指标为销轴质量Q最轻,Q的计算式为: Q(X)=Q(x1,x2)= 1 4 ??x21x2??=14??d2l??=0.00613d2l式中:设计变量X=[x1,x2]T=[d,l]T。 显然,d、l越小,Q值越小。但二者的取值要受强度、刚度 及结构尺寸等各种使用要求的限制,将以上各限制条件整理后可得: (1)抗弯强度条件: ??max= Pl 0.1d 3 ??[??w]????d3-8.33l??0(2)抗扭强度条件: ??max= M 0.2d3 ??[??]????d3-6.25??0(3)刚度条件: fmax= Pl33EJ=64Pl3 3E??d4 ??[f]????d4-0.34l3??0(4)结构尺寸限制: l??lmin????l-8??0 据此,可写出优化设计的数学模型为: minQ(X)=0.00613d2l????X=[x1,x2]=[d,l]Ts.t??d3-8.33l??0 d3-6.25??0d4-0.34l3??0l-8??0 这是一个单目标非线性二维约束优化问题。 1.1.2??求优结果 为了便于对采用不同手段得到的优化结果进行比较,以及确认MATLAB分析结果的可信度,特意选择了已经采用其它方法得到最优结果的例子(引自文献[5],并假定其结果正确),最优方案为:d*=4.309cm,l*=8cm,Q*=0.911kg。1.2??采用MATLAB软件对优化问题示例的分析1.2.1??约束函数的性态分析 应用MATLAB语言编程: [x,y]=meshgrid(linspace(0,17,19),linspace(0,11,19));??%根据函数的定义划分网格区域 Q=0.00613*y.^2*x; mesh(x,y,Q);??%通过三维网格模拟目标函数图形xlabel(??l??);ylabel(??d??); zlabel(??Q??);??%标注三坐标轴 title(??目标函数的图像??)??%标注图形名称 ??收稿日期:2003-09-25;修订日期:2004-01-30 作者简介:王春香(1962-),女,内蒙古包头人,副教授,硕士,主要研究方向:机械优化设计方法及CAD/CAM/CAI技术。 ????运行程序可画出目标函数的三维图形,如图2所示 同理,可作如下编程: ezplot(??x^3-8.33*y??,[0,10,0,20])??%弯曲强度条件的图形holdony=0??0.01??20;x=6.25^(1/3); plot(x,y,??k-??)??%扭转强度条件的图形holdon ezplot(??x^4-0.34*y^3??,[0,10,0,20])??%刚度条件的图形holdonx=0??0.01??10;y=8; plot(x,y,??k-??)??%长度边界条件的图形title(??各约束函数的图像??)text(6.5,11,??可行域??)holdoff 执行程序绘出各约束函数的图形,如图3所示。 从图2可以看出:目标函数的图像规则,即性态好,对于多数优化方法均适用;但由图3可行域的构成分析可知:d3 -6.25??0 (即抗扭强度条件)不是起作用约束, 设计时可不予考虑。?应用MATLAB软件的优化工具箱对优化问题示例求解 a.取设计变量的初值为:x0=[4??10]t 首先,编写目标函数的m文件:Objfun.m,返回x处的函数值f。 functionf=Objfun(x) f=0.00613*x(1)^2*x(2); b.因设计约束含3个非线性约束,故需编写一个描述非线性约束的m文件:NonLinConstr.m function[c,ceq]=NonLinConstr(x)c(1)=-x(1)^3+8.33*x(2); c(2) =-x(1)^3+6.