sjc.m

来自「力学计算 力学计算 力学计算 力学计算 力学计算」· M 代码 · 共 23 行

function [jacob2]=sjc(nnel,dhdr,dhds,xcoord,ycoord)

%------------------------------------------------------------------------
% 目的：
%     求二维雅可比行列式
%  变量：
%     jacob2 - 二维雅可比行列式
%     nnel - 每个单元节点数
%     dhdr - 自然坐标r对形函数的求导
%     dhds - 自然坐标s对形函数的求导
%     xcoord - 节点的x坐标值
%     ycoord - 节点的y坐标值
%------------------------------------------------------------------------

 jacob2=zeros(2,2);

 for i=1:nnel
 jacob2(1,1)=jacob2(1,1)+dhdr(i)*xcoord(i);
 jacob2(1,2)=jacob2(1,2)+dhdr(i)*ycoord(i);
 jacob2(2,1)=jacob2(2,1)+dhds(i)*xcoord(i);
 jacob2(2,2)=jacob2(2,2)+dhds(i)*ycoord(i);
 end