工程结构优化设计编程作业.docVIP

工程结构优化设计—编程作业 1.用黄金分割法求方程 x3-3x+1=0 在MATLAB中没有专门的函数实现黄金分割法求解线性方程，可通过编写glodf.m函数实现黄金分割法求解，其代码如下: function x=glodf(f,a,b,eps) if(nargin==3) eps=1.0e-4; end f1=subs(sym(f),findsym(sym(f)),a); f2=subs(sym(f),findsym(sym(f)),b); if(f1==0) x=a; end if(f2==0) x=b; end if(f1*f20) disp(’两端点函数值乘积大于0!’); return; else t1=a+(b-a)*0.382; t2=a+(b-a)*0.618; f_1=subs(sym(f),findsym(sym(f)),t1); f_2=subs(sym(f),findsym(sym(f)),t2); tol=abs(t1-t2); while(toleps) if(f_1*f_20) a=t1; b=t2; else fa=subs(sym(f),findsym(sym(f)),a); if(f_1*fa0) a=t2; else b=t1; end end t1=a+(b-a)*0.382; t2=a+(b-a)*0.618; f_1=subs(sym(f),findsym(sym(f)),t1); f_2=subs(sym(f),findsym(sym(f)),t2); tol=abs(t2-t1); end x=(t1+t2)/2; end 其实现的MATLAB代码如下: clear all; tic x=glodf(x^3-3*x+1,-1,1,1e-6) x = 0.3473 2.用抛物线法求解y=x 在MATLAB中没有专门的函数实现抛物线法求解线性方程，可通过编写paowuxian.m函数实现抛物线法求解，其代码如下: function x=paowuxian(f,x0,x1,x2,e) if nargin4,e=1e-4; end x=x2;y=x1;z=x0; while abs(x-y)e h1=y-z; h2=x-y; c1=(feval(f,y)-feval(f,z))/h1; c2=(feval(f,x)-feval(f,y))/h2; d=(c1-c2)/(h2+h1); w=c2+h2*d; xi=x-(2*feval(f,x))/(w+(w/abs(w))*sqrt(w^2-4*feval(f,x)*d)); z=y; y=x; x=xi; end 其实现的MATLAB代码如下: fun=inline(x^3+2*x-0.2); paowuxian(fun,0,0.5,1,1e-6),format short ans = 0.0995 3.用齿形法对下图的结构进行杆件截面优化设计： 节点处作用两种工况：（1）P1=2000kN，P2=0kN；（2）P1=0kN，P2=2000kN。各杆采用同一材料制成，弹性模量E为常数，容许应力分别为：[σ+]=2╳107kPa，[σ-]=1.5╳107kPa。 此结构为对称结构，工况也对称，所以A1=A3，计算时只取两个变量A1和A2，一种工况P1=2000kN，P2=0kN。 目标函数取为结构总重量： 在MATLAB中没有专门的函数实现齿形法求解线性方程，可通过编写Zigzag.m函数实现齿形法求解，其代码如下: function [newA1,newA2,newA3,W,exitflag] =Zigzag( A1,A2,maxstep) if nargin3 maxstep=96; end A=[A1;A2]; U=[(1e-4)*(sqrt(2)*A1+A2)/(sqrt(2)*A1^2+2*A1*A2);(1e-4)*(sqrt(2)*A1)/(sqrt(2)*A1^2+2*A1*A2)]; w=sqrt(2)*2*A1+A2; w1=w; w2=w; newA1=A1; newA2