弹塑性固体.rar_变形_变形 有限元_塑性变形_塑性有限元_弹塑性变形

function [ke]=quad4ke2(iopt,ecoords,E,v,t); if iopt ~=1 %平面应变问题 t=1; end ke=zeros(8,8); %2个积分点及权重 [pointr,weightr]=glqd1(2);%r方向积分点及权系数 [points,weights]=glqd1(2);%s方向积分点及权系数 %计算D矩阵 D=matiso(ipot,E,v); for i=1:2 %对积分点循环 for j=1:2 %计算形函数及其导数 [N,dN]=Sfisot4(pointr(i),points(j)); %计算J矩阵 [detJ,ivsJ,dNxyz]=jacobian(dN,ecoords); %计算B矩阵 dNdx=dNxyz(1,1); dNdy=dNxyz(2,1); B1=[dNdx 0;0 dNdy;dNdy dNdx]; dNdx=dNxyz(1,2); dNdy=dNxyz(2,2); B2=[dNdx 0;0 dNdy;dNdy dNdx]; dNdx=dNxyz(1,3); dNdy=dNxyz(2,3); B3=[dNdx 0;0 dNdy;dNdy,dNdx]; dNdx=dNxyz(1,4); dNdy=dNxyz(2,4); B4=[dNdx 0;0 dNdy;dNdy,dNdx]; B=[B1 B2 B3 B4]; ke=ke+B'*D*B*t*detJ*weightr(i)*weights(j); end end