FEM-Course-Matlab/18. 材料非线性问题matlab有限元编程/elastic/ElementArea.m

19 lines
519 B
Matlab

function area = ElementArea( ie )
% 셕炬데禱충생
% 渴흙꽝鑒:
% ie ---- 데禱뵀
% 럿쀼令:
% area ---- 데禱충생
global gNode gElement gMaterial
xi = gNode( gElement( ie, 1 ), 1 ) ;
yi = gNode( gElement( ie, 1 ), 2 ) ;
xj = gNode( gElement( ie, 2 ), 1 ) ;
yj = gNode( gElement( ie, 2 ), 2 ) ;
xm = gNode( gElement( ie, 3 ), 1 ) ;
ym = gNode( gElement( ie, 3 ), 2 ) ;
ai = xj*ym - xm*yj ;
aj = xm*yi - xi*ym ;
am = xi*yj - xj*yi ;
area = abs((ai+aj+am)/2) ;
return