FEM-Course-Matlab/10.四节点八节点四边形单元悬臂梁的Matlab有限元编程/4_nodes/Q4_G.m

12 lines
371 B
Mathematica
Raw Permalink Normal View History

2024-01-28 16:46:36 +00:00
function [G] = Q4_G(kesi,yita)
%Q4_G <EFBFBD><EFBFBD><EFBFBD><EFBFBD>G<EFBFBD><EFBFBD><EFBFBD><EFBFBD>
% <EFBFBD><EFBFBD><EFBFBD><EFBFBD><EFBFBD>κ<EFBFBD><EFBFBD><EFBFBD><EFBFBD><EFBFBD><EFBFBD><EFBFBD><EFBFBD><EFBFBD>ƫ<EFBFBD><EFBFBD>ֵ
% <EFBFBD><EFBFBD><EFBFBD><EFBFBD><EFBFBD><EFBFBD>ԪG<EFBFBD><EFBFBD><EFBFBD><EFBFBD>*(1/4)
G=[ yita-1, 0, 1-yita, 0,yita+1, 0,-yita-1, 0;
kesi-1, 0,-kesi-1, 0,kesi+1, 0, 1-kesi, 0;
0,yita-1, 0, 1-yita, 0,yita+1, 0,-yita-1;
0,kesi-1, 0,-kesi-1, 0,kesi+1, 0, 1-kesi];
G=G/4;
end