15 lines
451 B
Matlab
15 lines
451 B
Matlab
function N=ShapeFun(s,t) % 4级子程序，矩形单元形函数矩阵
|
||
%% 四边形八结点等参单元形函数矩阵
|
||
% 角点
|
||
N1=1/4*(1-s)*(1+t)*(-s+t-1);
|
||
N2=1/4*(1-s)*(1-t)*(-s-t-1);
|
||
N3=1/4*(1+s)*(1-t)*(s-t-1);
|
||
N4=1/4*(1+s)*(1+t)*(s+t-1);
|
||
|
||
% 边中点
|
||
N5=1/2*(1-t^2)*(1-s);
|
||
N7=1/2*(1-t^2)*(1+s);
|
||
N6=1/2*(1-s^2)*(1-t);
|
||
N8=1/2*(1-s^2)*(1+t);
|
||
N=[N1 0 N2 0 N3 0 N4 0 N5 0 N6 0 N7 0 N8 0;
|
||
0 N1 0 N2 0 N3 0 N4 0 N5 0 N6 0 N7 0 N8]; |