MATLAB多目标优化计算.pptVIP

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