sxy.m

有限元法计算平板受集中载荷的程序

function [dhdx,dhdy]=sxy2(nnel,dhdr,dhds,invjacob)

%------------------------------------------------------------------------
%  目的：
%     求dhdr,dhds在迪卡尔坐标下的值
%
%  变量：
%     dhdx - 形函数对迪卡尔坐标x的求导
%     dhdy - 形函数对迪卡尔坐标y的求导
%     nnel - 每个单元节点数
%     dhdr - 形函数对自然坐标r的求导
%     dhds - 形函数对自然坐标s的求导
%     invjacob - 求雅可比行列式的逆
%------------------------------------------------------------------------

 for i=1:nnel
 dhdx(i)=invjacob(1,1)*dhdr(i)+invjacob(1,2)*dhds(i);
 dhdy(i)=invjacob(2,1)*dhdr(i)+invjacob(2,2)*dhds(i);
 end