kjqx.m

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

subplot(1,2,1)
t=0:0.01:15;
x=sin(t); y=cos(t); z=t;
plot3(x,y,z)
view(-60,18)

subplot(1,2,2)
x1=[0.5,0.5,0,0.5]; y1=[0.5,0.5,0,0.5]; z1=[1,0,0,1];
x2=[0.5,0.5,1,0.5];  y2=[0.5,0.5,0,0.5]; z2=[1,0,0,1];
x3=[0,1]; y3=[0,0]; z3=[0,0];
plot3(x1,y1,z1,x2,y2,z2,x3,y3,z3)
hold on
fill3(x1,y1,z1,'g')
fill3(x2,y2,z2,'r')
view(24,30)