ssf.m

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

function [shapeq4,dhdrq4,dhdsq4]=ssf(rvalue,svalue)

%------------------------------------------------------------------------
%  目的：
%      在自然坐标下计算形函数和对其相应的求导
%
%  变量：
%     shapeq4 - 四节点单元形函数
%     dhdrq4 - 形函数对r求导
%     dhdsq4 - 形函数对s求导
%     rvalue - 对应点的r坐标值
%     svalue - 对应点的s坐标值　
%
%  说明：
%     第一个点自然坐标(-1,-1), 第二个点自然坐标 (1,-1)
%     第三个点自然坐标(1,1), 第四个点自然坐标(-1,1)
%------------------------------------------------------------------------

% 形函数

 shapeq4(1)=0.25*(1-rvalue)*(1-svalue);
 shapeq4(2)=0.25*(1+rvalue)*(1-svalue);
 shapeq4(3)=0.25*(1+rvalue)*(1+svalue);
 shapeq4(4)=0.25*(1-rvalue)*(1+svalue);

% 相应的求导

 dhdrq4(1)=-0.25*(1-svalue);
 dhdrq4(2)=0.25*(1-svalue);
 dhdrq4(3)=0.25*(1+svalue);
 dhdrq4(4)=-0.25*(1+svalue);

 dhdsq4(1)=-0.25*(1-rvalue);
 dhdsq4(2)=-0.25*(1+rvalue);
 dhdsq4(3)=0.25*(1+rvalue);
 dhdsq4(4)=0.25*(1-rvalue);