?? stiffnessmatrix.m
字號:
function k= StiffnessMatrix( ie ) % (定義單剛計算函數(shù))
global node element material K
k = zeros( 6, 6 ) ; % (零矩陣)
E = material( element(ie, 3), 1 ) ; %(彈性模量)
A = material( element(ie, 3), 2 ) ; %(截面積)
I=material(element(ie,3),3); %慣性據(jù)
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) ; %(單元長度)
k = [ 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 ];
t=trans(ie);
k = t*k*transpose(t) ; %(坐標轉(zhuǎn)換)
return
?? 快捷鍵說明
復(fù)制代碼
Ctrl + C
搜索代碼
Ctrl + F
全屏模式
F11
切換主題
Ctrl + Shift + D
顯示快捷鍵
?
增大字號
Ctrl + =
減小字號
Ctrl + -