ssb.m

力学计算 力学计算 力学计算 力学计算 力学计算

function [kinmtps]=ssB(nnel,dhdx,dhdy,shape)

%------------------------------------------------------------------------
%  目的：
%     计算剪切几何函数矩阵
%
%  变量：
%     nnel - 每个单元节点数
%     dhdx - 形函数对迪卡尔坐标x的求导
%     dhdy - 形函数对迪卡尔坐标y的求导
%     shape - 形函数
%------------------------------------------------------------------------

 for i=1:nnel
 i1=(i-1)*3+1;  
 i2=i1+1;
 i3=i2+1;
 kinmtps(1,i1)=-shape(i);
 kinmtps(1,i3)=dhdx(i);
 kinmtps(2,i2)=-shape(i);
 kinmtps(2,i3)=dhdy(i);
 end