清晰;
clc;
一个= xlsread (“分析DE ESTRUCTURAS _ INPUT2”,“彭”,“B1: B6”);
我= xlsread (“分析DE ESTRUCTURAS _ INPUT2”,“彭”,“C1: C6”);
E = xlsread (“分析DE ESTRUCTURAS _ INPUT2”,“彭”,“D1: D6”);
W = xlsread (“分析DE ESTRUCTURAS _ INPUT2”,“彭”,“u1: u6”);
L = xlsread (“分析DE ESTRUCTURAS _ INPUT2”,“彭”,“Q1: Q6”);
一个= xlsread (“分析DE ESTRUCTURAS _ INPUT2”,“彭”,“R1: R6”);
n = xlsread (“分析DE ESTRUCTURAS _ INPUT2”,“彭”,V3的);
nGDL = xlsread (“分析DE ESTRUCTURAS _ INPUT2”,“NUDOS”,“J2”);
GDLG = xlsread (“分析DE ESTRUCTURAS _ INPUT2”,“彭”,“K2: P4”);
为i = 1: n
=(我);
我=我(我);
E = E(我);
L = L(我);
=(我);
吉隆坡(i) = [E * / L 0 0 - E * / L 0 0; 0 12 * E *我/ L ^ 3 6 * E *我/ L ^ 2 0 -12 * E *我/ L ^ 3 6 * E * / L ^ 2; 0 6 * E *我/ L ^ 2 *我/ L 0 6 * 4 * E E *我/ L ^ 2 E 2 * *我/ L; - E * / L 0 0 E * / L 0 0; 0 -12 * E *我/ L ^ 3 6 * E * / L ^ 2 0 12 * E *我/ L ^ 3 6 * E * / L ^ 2; 0 6 * E *我/ L ^ 2 *我/ L 0 6 * 2 * E E *我/ L ^ 2 4 * E *我/ L);
T (i) = (cos (a)罪(a) 0 0 0 0; sin (a)因为(a) 0 0 0 0; 0 0 1 0 0 0, 0 0 0 cos (a) sin () 0; 0 0 0 sin()因为(a) 0, 0 0 0 0 0 1);
公斤(i) = T (i)的(我)* T *吉隆坡(我);
G = [GDLG(我,1)GDLG(我,2)GDLG(我,3)GDLG(我,4)GDLG(我,5)GDLG(我,6)];
KT = 0 (nGDL nGDL);
KT (G, G) =公斤(我);
如果我= = 1
K = 0 (nGDL nGDL);
结束
K = K + KT;
结束
disp (K)