平面四节点等参单元matlab实现(7)

if elem(j,1)-c(i)==0 g=g+1;

sigmar(:,i)=sigmar(:,i)+sigma1(:,j); end

if elem(j,2)-c(i)==0 g=g+1;

sigmar(:,i)=sigmar(:,i)+sigma2(:,j); end

if elem(j,3)-c(i)==0 g=g+1;

sigmar(:,i)=sigmar(:,i)+sigma3(:,j); end

if elem(j,4)-c(i)==0 g=g+1;

sigmar(:,i)=sigmar(:,i)+sigma4(:,j); end end

sigmar(:,i)=sigmar(:,i)/g; %求应力绕节点平均

msigmar(:,i)=(0.5*((sigmar(1,i)-sigmar(2,i))^2+sigmar(1,i)^2+sigmar(2,i)^2+6*(sigmar(3,i))^2))^0.5; %求节点处的mises应力 end

msigmar %mises应力

for i=1:f %同上 g=0; for j=1:m

if elem(j,1)-d(i)==0 g=g+1;

sigmat(:,i)=sigmat(:,i)+sigma1(:,j); end

if elem(j,2)-d(i)==0 g=g+1;

sigmat(:,i)=sigmat(:,i)+sigma2(:,j); end

if elem(j,3)-d(i)==0 g=g+1;

sigmat(:,i)=sigmat(:,i)+sigma3(:,j); end

if elem(j,4)-d(i)==0 g=g+1;

sigmat(:,i)=sigmat(:,i)+sigma4(:,j); end end