trans.m

是个好东西

function t=trans(ie)
global node element
xi = node( element( ie, 1 ), 1 ) ;  %(两端座标)
yi = node( element( ie, 1 ), 2 ) ;  %(两端座标)
xj = node( element( ie, 2 ), 1 ) ;  %(两端座标)
yj = node( element( ie, 2 ), 2 ) ;  %(两端座标)
L = ( (xj-xi)^2 + (yj-yi)^2)^(1/2) ;  %(单元长度)
l1 = (xj-xi)/L ;    %(cosα)
m1 = (yj-yi)/L ;    %(cosβ)
      t=[  l1   -m1   0     0     0     0
           m1   l1    0     0     0     0
           0     0    1     0     0     0
           0    0     0     l1   -m1    0
           0    0     0     m1    l1    0
           0    0     0     0     0     1] ;    %(转换矩阵)
return;