FEM-Course-Matlab/12+13框架结构动力分析(模态分析+时程分析)/FrameElementKeLocal.m

12 lines
577 B
Matlab

%¼ÆËãµ¥Ôª¸Õ¶È¾ØÕóº¯Êý
function [Ke_local] = FrameElementKeLocal(A,E,I,R,BarLength)
L=BarLength;
Ke_local=[E*A/L 0 0 -E*A/L 0 0
0 12*E*I/L^3 6*E*I/L^2 0 -12*E*I/L^3 6*E*I/L^2;
0 6*E*I/L^2 4*E*I/L 0 -6*E*I/L^2 2*E*I/L ;
-E*A/L 0 0 E*A/L 0 0
0 -12*E*I/L^3 -6*E*I/L^2 0 12*E*I/L^3 -6*E*I/L^2;
0 6*E*I/L^2 2*E*I/L 0 -6*E*I/L^2 4*E*I/L ];
end