昆明理工大学机械优化设计作业.docxVIP

1.二力杆的受力分析钢管所受的压力： 钢管所受的压力：Fl= 钢管所受的压力：Fl=FV(B2+h 钢管所受的压力：Fl= FV(B2+h2) 压杆失稳的临界力：F「= 兀2 EI 13 A其中：I为钢管截面惯性矩：1=1（72+。2） 其中：A为钢管截面面积：A=7iTD钢管所受的压应力：。 钢管所受的压应力：。 钢管所受的压应力：。_Ft _ 钢管所受的压应力：。 _Ft _ fV(S2+/i2) nTDh 钢管的临界应力：（Jc= 钢管的临界应力：（Jc 钢管的临界应力：（Jc= Fc n2E(T2+D2) A 8(F2+/i2) 有题目中给出数据：2F=300000N 2B=1520mm T=2.5mm 有题目中给出数据：2F=300000N 2B=1520mm T=2.5mm即 F=150000N 有题目中给出数据：2F=300000N 2B=1520mm T=2.5mm 即 F=150000N 即 B=760mm p - 8300 kg!nr crs = 700 Mpa crs crs = 700 Mpa E=2Ax\Q]] N/m2 解：1.建立数学模型设计变量： x1 1X == X2 J目标函数： f（x） = 2K兀dXb? +H?= 1.225x 1001^577600+xf约束条件： 1）圆管杆件中的压应力。应小于或等于内，即F^B F^B2 +H2 a = F^B2 +H F^B2 +H2 a = 7fTDH by 于是得? ?(X ) = 19098.59 ?(X ?(X ) = 19098.59 T577600+7J 22 2）圆管杆件中的压应力。应小于或等于压杆稳定的临界应力/，由欧拉公式得钢管的压杆 温度应力』 °C =7T2EIl2a72石(。 °C = 7T2EI l2a 72石(。2+72)8(B2+H2) = 2.6xl05 Xj2 + 6.25577600+xf a/577600+x^s 99g2(x) = b -ae= 19098.59^ — 2.6x 1()5(% ： + 6.25)/(57-00 + x^)0 xxx23）设计变量的值不得小于或等于0 于是得 g3(X)= -X] 0g2(X)= —%2 40 2.从以上分析可知，该优化设计问题具有2个设计变量，4个约束条件，按优化方法程序 的规定编写数学模型的程序如下： /*用来解n维具有不等式约束优化问题，文件名funct.c*//*以下子程序用来构造惩罚函数*/ #include stdio.h #include “stdlib.h/* /*常量，定义约束条件个数*/ /*全局变量，惩罚因子*/ /*目标函数子程序*/ /*常量，定义约束条件个数*//* /*常量，定义约束条件个数*/ /*全局变量，惩罚因子*/ /*目标函数子程序*/ ff^l.225e4*x[0]*sqrt(577600.0+x[l]*x[l]); retum(ff);} void strain(double x口,double g[])/*约束条件子程序*/{g[0]= 19098.59*sqrt(577600.0+x[l]*x[l])/(x[0]*x[l])-700.0; g[l]= 19098.59*sqrt(577600.0+x[l]*x[l)]/(x[0]*x[l]))- 2.6e5*(x[0]*x[0]+6.25)/(577600.0+x[l]*x[l]);g[2-x[0]; gx(3)=-x[l];double objf(double p[])/*惩罚函数子程序*/ {int i; double ff,sg,*g; g=(doublc *)malloc(kkg *sizcof(doublc));sg=O; strain(p9g);fbr(i=O;ikkg;i++) {if(*(g+l)0)sg=sg+rO/(*(g+i)); elsesg=sg+r0*(le+10); }free(g); ff^f(p)+sg;return(fl); /*内点惩罚函数主程序，文件名mldfhs.c*/#include powell.c* /*调用 powell 法函数 */ void main() {int i;double p[]={72,700}; double fdm?fxo,c?x[2];c=0.4; ro=l;fbm=100; do {fxo=powell(p,0.4,0000001,2,x); if(fabs(fom-fxo)0.000001) {fbm=fxo;r0=c*r0; fbr(i=0;i2;i++)*(p+i)=x[i]; }else {printf(“输出最优点及其目标函数值：\n)；return;