qtzt.m

四面体有限元程序,单元刚度的计算及总刚度矩征集成

subplot(1,2,1)
sphere(40)
axis equal
subplot(1,2,2)
t=0:0.1*pi:2*pi;
[X,Y,Z]=cylinder(4*cos(t));
mesh(X,Y,Z)
axis square