MATLAB多标优化计算.ppt

function y=barass(kk,k,i,j) %计算节点i，j整体刚度矩阵 kk(i,i)=kk(i,i)+k(1,1); kk(i,j)=kk(i,j)+k(1,2); kk(j,i)=kk(j,i)+k(2,1); kk(j,j)=kk(j,j)+k(2,2); y=kk; 杆单元 function y=barelefor(k,u) %计算单元节点力矢量 y=k*u%u位移 function y=barelestr(k,u,a) %计算单元应力矢量 y=k*u/a 杆单元 E=210GPa,A=0.003m^2,P=10kN,右端点位移0.002m,L12=1.5m,L23=1m 1 2 P 3 杆单元 e=210e6；a=0.003；l1=1.5；l2=1； k1=bareal(e,a,l1) %单元刚度矩阵 k2=bareal(e,a,l2) 杆单元 k1 = 420000 -420000 -420000 420000 k2 = 630000 -630000 -630000 630000 kk=zeros(3,3); %整体 kk=barass(kk,k1,1,2) kk=barass(kk,k2,2,3) 杆单元 kk = 420000 -420000 0 -420000 420000 0 0 0 0 kk = 420000 -420000 0 -420000 1050000 -630000 0 -630000 630000 P=10kN,节点3位移0.002m 0 u2 0.002 f1 -10 f3 k=kk(2,2);k0=kk(2,3); u0=0.0020;f=[-10];%边界条件 f0=f-k0*u0; u=k\f0%节点2位移 uu=[0;u;u0]; ff=kk*uu%支反力 杆单元 u = 0.0012 ff = -500.0000 -10.0000 510.0000 u1=[0;uu(2)] f1=barelefor(k1,u1) s1=barelestr(k1,u1,a) u2=[uu(2);uu(3)] f2=barelefor(k2,u2) s2=barelestr(k2,u2,a) f1 = -500.0000 500.0000 s1 = 1.0e+005 * -1.6667 1.6667 f2 = -510.0000 510.0000 s2 = 1.0e+005 * -1.7000 1.7000 杆单元 半径为r的轮子沿直线轨道无滚动的滑动．设轮子转角ψ=ωt，ω为常量。求轮缘上—点M的运动方程，并求该点的速度和加速度。 w=5;t=0:0.00005*pi:4*pi; r=0.5; x=r*w*t-r*sin(t); y=r-r*cos(t); x1=diff(x); y1=diff(y); x2=diff(x1); %x方向加速度 y2=diff(y1); % y方向加速度 c=sqrt(x1.^2+y1.^2);%速度 d=sqrt(x2.^2+y2.^2); figure(1); subplot(4,1,1),plot(x,y);ylabel(m);xlabel(m); title(运动轨迹,fontsize,16); subplot(4,1,2),plot(t(1:(length(t)-1)),c);ylabel(rad/s);xlabel(时间(t));title(速度图,fontsize,16); subplot(4,1,3),plot(t(1:(length(t)-2)),x2);ylabel(rad/s2);xlabel(时间(t));title(x方向加速度图,fontsize,16); subplot(4,1,4),plot(t(1:(length(t)-2)),y2);ylabel(rad／s2);xlabel(时间(t));title(y方向加速度图,fontsize,16) * 作业要求： （1）对该问题进行分析，写出该问题的物理模型； （2）将物理模型转化为优化模型（包括设计变量、目标函数、约束条件）； （3）将优化模型转化为matlab程序（m文件）； （4）利用matlab软件求解该优化问题，写出最优解。 （