My_coding – CellularBeam_Mesh ©
function matC2 = LocalMatrixC(Noded,GPoint)
% NUMERICAL INTEGRATION...
S35 = sqrt(3.0/5.0);
SN = [-S35 -S35; S35 -S35; S35
0 -S35; S35
0;
0
S35; -S35
S35; -S35
S35; ...
0; 0 0 ];
% LOCAL C MATRIX...
N2 = 2*Noded;
for i = 1:GPoint
s = SN(i,1); n = SN(i,2);
% SHAPE FUNCTION...
matC1 = [
1/4*(n-1)*s*n+(1/4*s-1/4)*(n-1)*n ...
-(1/4*s-1/4)*s*n-(1/4*s-1/4)*(n-1)*s ...
1/4*(n-1)*s*n+(1/4*s+1/4)*(n-1)*n ...
-(1/4*s+1/4)*s*n-(1/4*s+1/4)*(n-1)*s ...
1/4*(n+1)*s*n+(1/4*s+1/4)*(n+1)*n ...
-(1/4*s+1/4)*s*n-(1/4*s+1/4)*(n+1)*s ...
1/4*(n+1)*s*n+(1/4*s-1/4)*(n+1)*n ...
-(1/4*s-1/4)*s*n-(1/4*s-1/4)*(n+1)*s ...
-1/2*(s+1)*(n-1)*n+(-1/2*s+1/2)*(n-1)*n ...
-(-1/2*s+1/2)*(s+1)*n-(-1/2*s+1/2)*(s+1)*(n-1) ...
-1/2*s*(n-1)*(n+1)+(-1/2*s-1/2)*(n-1)*(n+1) ...
-(-1/2*s-1/2)*s*(n+1)-(-1/2*s-1/2)*s*(n-1) ...
-1/2*(s+1)*(n+1)*n+(-1/2*s+1/2)*(n+1)*n ...
-(-1/2*s+1/2)*(s+1)*n-(-1/2*s+1/2)*(s+1)*(n+1) ...
-1/2*s*(n-1)*(n+1)+(-1/2*s+1/2)*(n-1)*(n+1) ...
-(-1/2*s+1/2)*s*(n+1)-(-1/2*s+1/2)*s*(n-1) ...
(s+1)*(n-1)*(n+1)+(s-1)*(n-1)*(n+1) ...