nh.m

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

%%program nh.m
t=0:0.5:5;
s=2*t+0.4.*t.*t;
s1=[0.05  -0.03  -0.01  0.03  0.01 ...
    0.01  -0.02  -0.01  0.02  0.03  0.01];
ss=s+s1;
P=polyfit(t,ss,2);
y=polyval(P,t);
plot(t,y,t,ss,'+')