13矩阵位移法教学资料.docVIP

矩阵位移法教学资料 一 单元刚度矩阵形成程序 1单元刚度矩阵各元素的算式 将局部坐标系下的单元刚度矩阵、坐标变换矩阵代入变换式，得结构坐标系下单元刚度矩阵各元素的算式： 在计算机程序中，单元刚度矩阵各元素一般直接按其算式确定，然后用一个二维数组KE(6,6)存放，程序中： EL─单元杆件的长度 EF(N,3),EF(N,4) ─单元的抗拉刚度EA、抗弯刚度EI 2 形成单元刚度矩阵的子程序 SUBROUTINE SETKE(N,EF,EL,SI,CO,KE) REAL KE DIMENSION EF(100,4),KE(6,6) A=EF(N,3)*CO*CO/EL+12*EF(N,4)*SI*SI/EL**3 KE(1,1)=A KE(1,4)=-A KE(4,4)=A B=(EF(N,3)/EL-12*EF(N,4)/EL**3)*CO*SI KE(1,2)=B KE(1,5)=-B KE(2,4)=-B KE(4,5)=B C=6*EF(N,4)*SI/EL**2 KE(1,3)=-C KE(1,6)=-C KE(3,4)=C KE(4,6)=C D=EF(N,3)*SI*SI/EL+12*EF(N,4)*CO*CO/EL*