FEM-Course-Matlab/4.1 薄板厚板模态分析/Jacobian.m

10 lines
319 B
Matlab

function[jacobian]=Jacobian(shapeFunction,elemCoordinates)
% Compute Jacobian and calculate global cordinates
jacobian=struct();
jacobian.matrix=elemCoordinates'*shapeFunction.dfun;
jacobian.invmatrix=inv(jacobian.matrix);
jacobian.globalDerivatives=(jacobian.invmatrix*shapeFunction.dfun')';%对物理坐标求偏导
return