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

图3 求选取点的mises应力程序：

clc

A=xlsread('D:\data','ansys');

c=[24,29,47,58,78,79,137,149,160,166,186]';

d=[166,167,168,169,170,171,172,173,174,175,176,177,178,179,180,181,182,183,184,185]';

e=size(c,1) f=size(d,1); for i=1:e

ansysr(i)=(0.5*((A(c(i),1)-A(c(i),2))^2+A(c(i),1)^2+A(c(i),2)^2+6*(A(c(i),3))^2))^0.5;

end for i=1:f

ansyst(i)=(0.5*((A(d(i),1)-A(d(i),2))^2+A(d(i),1)^2+A(d(i),2)^2+6*(A(d(i),3))^2))^0.5;

end ansysr ansyst