19 lines
745 B
Mathematica
19 lines
745 B
Mathematica
|
function Nst=DiffShapeFun(s,t)
|
|||
|
|
|
|||
|
|
N1s=(s/4 - 1/4)*(t + 1) + ((t + 1)*(s - t + 1))/4;
|
|||
|
|
N2s=-(s/4 - 1/4)*(t - 1) - ((t - 1)*(s + t + 1))/4;
|
|||
|
|
N3s=((t - 1)*(t - s + 1))/4 - (s/4 + 1/4)*(t - 1);
|
|||
|
|
N4s=(s/4 + 1/4)*(t + 1) + ((t + 1)*(s + t - 1))/4;
|
|||
|
|
N5s=t^2/2 - 1/2;
|
|||
|
|
N6s=s*(t-1);
|
|||
|
|
N7s=1/2 - t^2/2;
|
|||
|
|
N8s=-s*(t + 1);
|
|||
|
|
% <EFBFBD><EFBFBD><EFBFBD><EFBFBD> t <EFBFBD><EFBFBD><EFBFBD><EFBFBD>
|
|||
|
|
N1t=(s/4 - 1/4)*(s - t + 1) - (s/4 - 1/4)*(t + 1);
|
|||
|
|
N2t=-(s/4 - 1/4)*(t - 1) - (s/4 - 1/4)*(s + t + 1);
|
|||
|
|
N3t=(s/4 + 1/4)*(t - s + 1) + (s/4 + 1/4)*(t - 1);
|
|||
|
|
N4t=(s/4 + 1/4)*(t + 1) + (s/4 + 1/4)*(s + t - 1);
|
|||
|
|
N5t=t*(s - 1); N6t=s^2/2 - 1/2;
|
|||
|
|
N7t=-t*(s + 1); N8t=1/2 - s^2/2;
|
|||
|
|
Nst=[N1s N1t;N2s N2t;N3s N3t;N4s N4t;N5s N5t;N6s N6t;N7s N7t;N8s N8t;];
|
|||
|
|
end
|