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

sigmat(:,i)=sigmat(:,i)/g;

msigmat(:,i)=(0.5*((sigmat(1,i)-sigmat(2,i))^2+sigmat(1,i)^2+sigmat(2,i)^2+6*(sigmat(3,i))^2))^0.5; end msigmat

四、计算结果：

1.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]';

其在图中的位置为图1中红线标注：