yi00si
yicosyi0000cosyisi
yi00si
yicosyie
d悬索的转换矩阵s2zeros101fori110s2isqrtf2i1f2i2L2e
dKbO2zeros440
ffori110KbO2144i34iEAbs2i1010000010100000e
d悬索的单元刚度矩阵f110f129L20f134L5f1421L20f156L5f165L4f176L5f1821L20f194L5f1109L20f1110zzeros101fori110ziata
f1i1f1iLe
dT4zeros660fori110T46i56i6i56icoszisi
zi0000si
zicoszi0000001000000coszisi
zi0000si
zicoszi0000001e
d拱的转换矩阵s3zeros101fori110s3isqrtf1i1f1i2L2e
dKaOzeros660fori110KaO166i56iEAas3i00EAas3i00012EIas3is3is3i6EIas3is3i012EIas3is3is3i6EIas3is3i06EIas3is3i4EIas3i06EIas3is3i2EIas3iEAas3i00EAas3i00012EIas3is3is3i6EIas3is3i012EIas3is3is3i6EIas3is3i06EIas3is3i2EIas3i06EIas3is3i4EIas3ie
d拱的单元刚度矩阵T50100100000010010吊杆的转换矩阵s4zeros91
ffori19s4if2i1f1i1e
dKtOzeros436fori19KtO144i34iEAts4i1010000010100000e
d吊杆的单元刚度矩阵KcT1KcOT1总体坐标下主塔的单元刚度矩阵Kb1T2KbO1T2Kb11T21KbO1T21总体坐标下斜索的单元刚度矩阵Kb2zeros440fori110T3OT3144i34iKb2144i34iT3OKbO2144i34iT3Oe
d总体坐标下悬索的单元刚度矩阵Kazeros660fori110T4OT46i56i6i56iKa166i56iT4OKaO166i56iT4Oe
d总体坐标下拱的单元刚度矩阵Ktzeros436fori19KtOOKtO144i34iKt144i34iT5KtOOT5e
d总体坐标下吊杆的单元刚度矩阵定义51阶0矩阵K1zeros5151K2zeros5151K3zeros5151K4zeros5151K5zeros5151Xzeros5151Yzeros5151Zzeros5151把主塔整合到整体刚度矩阵中：K11313KcO4646K122242224KcO4646把斜索整合到整体刚度矩阵中：K21212Kb13434K222232223Kb111212把悬索整合到整体刚度矩阵中K31212KbO21212K31245KbO21234fori210
fX2i2i32i2i3KbO2144i34iK3K3Xe
d把拱整合到整体刚度矩阵中：K425272527KaO4646K449514951KaO135557fori29Y3i193i243i193i24KaO166i56iK4K4Ye
d把吊杆整合到整体刚度矩阵中：fori19Z2i22i32i22i3KtO1212Z2i22i33i223i23KtO1234Z3i223i232i22i3KtO3412Z3i223i233i223i23Kr