15 lines
451 B
Mathematica
15 lines
451 B
Mathematica
|
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];
|