ts.m

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

[X,Y]=meshgrid(-8:.5:8);
R=sqrt(X.^2+Y.^2)+eps;
Z=sin(R)./R;
colormap(pink)

subplot(2,2,1)
surfl(X,Y,Z)

subplot(2,2,2)
surfl(X,Y,Z,[40,60],[0.1,0.1,0.7,1])

subplot(2,2,3)
surfl(X,Y,Z)
shading flat

subplot(2,2,4)
surfl(X,Y,Z,[40,60],[0.1,0.1,0.7,1])
shading interp