矩阵位移法大作业
结构力学(9.14.1)--矩阵位移法习题2
5kN m
8m 8m
8m
三 . 整体分析
12. 试求图示结构 ( 不计轴变 ) 的荷载列阵 ( 先处理法 ).
1(1,0,2) 2(1,0ቤተ መጻሕፍቲ ባይዱ3) 3(1,0,3)
X1
X2
4(0,0,0)
P
X
1
0
X
2
0
四 . 求杆端力
1. 连续梁在一般荷载作用下 , 单元杆端力由下式计算 . 是否正确 ?
6
48
4
2
1(0,0,0)
12
1 6
k
6
48
4(1,0,3)
3
2(0,0,0)
3
1
2
3
例 . 不计轴变 , 作弯矩图
已知 : 各杆长均为 12m, 线刚度均为 12
P 10kN, q 5kN / m
P 10kN, q 5kN / m
解 : 1 6 1 6
k
1
6
1
48 6
6 1
24
6
6
24
6
48
3(1,0,2)
2
1
1 6 1 6 1 0
k
1
6 1
48 6
6 1
24
2
0
63 1
6 24
EI
EI
EA 2l
2 2
l
l
三 . 整体分析
4(1,0,0)
5(1,0,0)
矩阵位移法题目及答案
1．作图示刚架的F、S F、M图，已知各杆截面均为矩形，柱截面宽N0.4m,高0.4m, 大跨梁截面宽0.35m,高0.85m，小跨梁截面宽0.35m,高0.6m，各杆E=3.0×104 MPa。
10分2、计算图示桁架各杆的轴力。
已知A=2400mm2,E=2.0×105MPa。
5分3．作图示连续梁的F、M图，已知各梁截面面积A=6.52m,惯性矩SI=5.504m，各杆E=3.45×104MPa。
5分答案1******************************************************************** ** 1 composite beam 2012.10.17 ** ********************************************************************3E10 16 13 9 11 2 0.2975 17.912E-32 3 0.2975 17.912E-33 4 0.21 6.3E-31 5 0.16 2.133E-33 6 0.16 2.133E-34 7 0.16 2.133E-35 6 0.2975 17.912E-36 7 0.21 6.3E-35 8 0.16 2.133E-36 9 0.16 2.133E-37 10 0.16 2.133E-38 9 0.2975 17.912E-39 10 0.21 6.3E-38 11 0.16 2.133E-39 12 0.16 2.133E-310 13 0.16 2.133E-3 0 10.93.8 10.97.6 10.911.4 10.90 7.77.6 7.711.4 7.70 4.57.6 4.511.4 4.50 07.6 011.4 0111 0112 0113 0121 0122 0123 0131 0132 0133 041 100E3 0 02 0 0 -15E33 0 0 -15E35 100E3 0 0712 2 -26E3 3.813 2 -26E3 2.77 4 -36E3 7.68 4 -36E3 3.812 4 -36E3 7.613 4 -36E3 3.814 3 20E3 4.5第一题结果******************************************************************* * * * 1 composite beam 2012.10.17 * * * *******************************************************************The Input DataThe General InformationE NM NJ NS NLC3.000E+10 16 13 9 1The Information of Membersmember start end A I1 12 2.975000E-01 1.791200E-022 23 2.975000E-01 1.791200E-023 34 2.100000E-01 6.300000E-034 15 1.600000E-01 2.133000E-035 36 1.600000E-01 2.133000E-036 47 1.600000E-01 2.133000E-037 5 6 2.975000E-01 1.791200E-028 6 7 2.100000E-01 6.300000E-039 5 8 1.600000E-01 2.133000E-0310 6 9 1.600000E-01 2.133000E-0311 7 10 1.600000E-01 2.133000E-0312 8 9 2.975000E-01 1.791200E-0213 9 10 2.100000E-01 6.300000E-0314 8 11 1.600000E-01 2.133000E-0315 9 12 1.600000E-01 2.133000E-0316 10 13 1.600000E-01 2.133000E-03The Joint Coordinatesjoint X Y1 .000000 10.9000002 3.800000 10.9000003 7.600000 10.9000004 11.400000 10.9000005 .000000 7.7000006 7.600000 7.7000007 11.400000 7.7000008 .000000 4.5000009 7.600000 4.50000010 11.400000 4.50000011 .000000 .00000012 7.600000 .00000013 11.400000 .000000The Information of SupportsIS VS111 .000000112 .000000113 .000000121 .000000122 .000000123 .000000131 .000000132 .000000133 .000000( NA= 357 )( NW= 1167 )Loading Case 1The Loadings at JointsNLJ= 4ILJ PX PY PM1 100000.0000 .0000 .000002 .0000 .0000 -15000.000003 .0000 .0000 -15000.00000 5 100000.0000 .0000 .00000The Loadings at MembersNLM= 7ILM ITL PV DST12 2 -26000.0000 3.80000013 2 -26000.0000 2.7000007 4 -36000.0000 7.6000008 4 -36000.0000 3.80000012 4 -36000.0000 7.60000013 4 -36000.0000 3.80000014 3 20000.0000 4.500000The Results of CalculationThe Joint Displacementsjoint u v phi1 1.845349E-02 -1.982711E-04 -1.263100E-042 1.841771E-02 -3.424398E-04 -8.180773E-063 1.838193E-02 -5.043591E-04 -1.356524E-044 1.836317E-02 -3.892198E-04 -1.683554E-045 1.608566E-02 -2.069957E-04 -9.278065E-046 1.601139E-02 -5.147233E-04 6.593305E-057 1.599555E-02 -3.701310E-04 -4.819689E-048 1.132049E-02 -1.535800E-04 -1.283845E-039 1.131820E-02 -3.849935E-04 3.869225E-0510 1.130585E-02 -2.796765E-04 -9.193725E-0411 7.105234E-18 -1.638186E-17 -1.781240E-1712 9.610834E-18 -4.106598E-17 -2.156936E-1713 7.783932E-18 -2.983216E-17 -1.882119E-17The Terminal Forcesmember N(st) Q(st) M(st) N(en) Q(en) M(en)1 84035.890 -13086.980 -41569.990 -84035.890 13086.980 -8160.5532 84035.890 -13086.980 -6839.447 -84035.890 13086.980 -42891.1003 31099.080 -28633.300 -52776.720 -31099.080 28633.300 -56029.8104 -13086.980 15964.110 41569.990 13086.980 -15964.110 9515.1395 -15546.320 52936.810 80667.820 15546.320 -52936.810 88729.9806 28633.300 31099.080 56029.810 -28633.300 -31099.080 43487.2307 87221.910 93210.620 -62622.290 -87221.910 180389.400 -268657.0008 26256.560 29751.550 -2861.126 -26256.560 107048.500 -144003.0009 80123.640 28742.190 53107.150 -80123.640 -28742.190 38867.85010 194594.600 113902.200 182788.200-194594.600-113902.200 181698.70011 135681.800 57355.640 100515.700-135681.800 -57355.640 83022.30012 2689.851 83695.000 -146729.300 -2689.851 215905.000 -355668.70013 20483.680 160.171 -42824.000 -20483.680 162639.800 -245087.30014 163818.600 26052.340 107861.500-163818.600-116052.300 211874.00015 410659.800 96108.340 216794.000-410659.800 -96108.340 215693.50016 298321.600 77839.320 162065.000-298321.600 -77839.320 188211.900( NA= 357 )单位（N m）( NW= 1195 )第二题答案******************************************************************* * * * 2 composite beam 2012.10.17 * * * ******************************************************************* 2E11 14 9 4 11 2 2.4E-3 1E-102 3 2.4E-3 1E-103 4 2.4E-3 1E-104 5 2.4E-3 1E-101 8 2.4E-3 1E-101 6 2.4E-3 1E-102 6 2.4E-3 1E-103 6 2.4E-3 1E-103 7 2.4E-3 1E-104 7 2.4E-3 1E-105 7 2.4E-3 1E-105 9 2.4E-3 1E-106 8 2.4E-3 1E-107 9 2.4E-3 1E-100 62 64 66 68 62 36 30 08 081 082 091 092 051 0 -50E3 02 0 -50E3 03 0 -50E3 04 0 -50E3 05 -10E3 -50E3 0第二题结果******************************************************************* * * * 2 composite beam 2012.10.17 * * * *******************************************************************The Input DataThe General InformationE NM NJ NS NLC2.000E+11 14 9 4 1The Information of Membersmember start end A I1 12 2.400000E-03 1.000000E-102 23 2.400000E-03 1.000000E-103 34 2.400000E-03 1.000000E-104 45 2.400000E-03 1.000000E-105 1 8 2.400000E-03 1.000000E-106 1 6 2.400000E-03 1.000000E-107 2 6 2.400000E-03 1.000000E-108 3 6 2.400000E-03 1.000000E-109 3 7 2.400000E-03 1.000000E-1010 4 7 2.400000E-03 1.000000E-1011 5 7 2.400000E-03 1.000000E-1012 5 9 2.400000E-03 1.000000E-1013 6 8 2.400000E-03 1.000000E-1014 7 9 2.400000E-03 1.000000E-10The Joint Coordinatesjoint X Y1 .000000 6.0000002 2.000000 6.0000003 4.000000 6.0000004 6.000000 6.0000005 8.000000 6.0000006 2.000000 3.0000007 6.000000 3.0000008 .000000 .0000009 8.000000 .000000The Information of SupportsIS VS81 .00000082 .00000091 .00000092 .000000( NA= 270 )( NW= 907 )Loading Case 1The Loadings at JointsNLJ= 5ILJ PX PY PM1 .0000 -50000.0000 .000002 .0000 -50000.0000 .000003 .0000 -50000.0000 .000004 .0000 -50000.0000 .000005 -10000.0000 -50000.0000 .00000The Loadings at MembersNLM= 0The Results of CalculationThe Joint Displacementsjoint u v phi1 -1.052370E-04 -9.375000E-04 -7.026900E-052 -1.746814E-04 -1.193735E-03 1.087907E-043 -2.441259E-04 -8.137530E-04 -3.230397E-054 -3.552370E-04 -1.302888E-03 -1.022944E-045 -4.663480E-04 -9.375000E-04 1.412285E-046 3.860367E-04 -8.812350E-04 1.226822E-047 -7.938921E-04 -9.903881E-04 -6.211760E-058 -3.833334E-18 -1.325000E-17 -2.257494E-049 2.833334E-18 -1.175000E-17 3.510750E-04The Terminal Forcesmember N(st) Q(st) M(st) N(en) Q(en) M(en)1 16666.660 .009 .007 -16666.660 -.009 .0112 16666.660 -.009 -.008 -16666.660 .009 -.0113 26666.660 .011 .011 -26666.660 -.011 .0104 26666.660 -.010 -.012 -26666.660 .010 -.0075 75000.000 -.001 -.003 -75000.000 .001 -.0046 -30046.240 -.002 -.004 30046.240 .002 -.0027 49999.980 -.002 -.003 -49999.980 .002 -.0038 39060.150 -.002 -.005 -39060.150 .002 -.0039 21032.390 .002 .004 -21032.390 -.002 .00310 49999.980 .002 .002 -49999.980 -.002 .00311 -30046.240 .002 .005 30046.240 -.002 .00212 75000.000 .001 .003 -75000.000 -.001 .00413 69106.410 .003 .008 -69106.410 -.003 .00414 51078.650 -.004 -.009 -51078.650 .004 -.004( NA= 270 )( NW= 907 )第三题答案******************************************************************** ** 3 composite beam 2012.10.17 ** ******************************************************************** 3.45E10 4 5 6 11 2 6.5 5.52 3 6.5 5.53 4 6.5 5.54 5 6.5 5.50 040 060 080 0120 011 012 013 022 042 052 013 0 -320E3 -100E341 4 -10.5E3 402 4 -10.5E3 203 4 -10.5E3 204 4 -10.5E3 40第三题结果******************************************************************* * * * 3 composite beam 2012.10.17 * * * *******************************************************************The Input DataThe General InformationE NM NJ NS NLC3.450E+10 4 5 6 1The Information of Membersmember start end A I1 12 6.500000E+00 5.500000E+002 23 6.500000E+00 5.500000E+003 34 6.500000E+00 5.500000E+004 45 6.500000E+00 5.500000E+00The Joint Coordinatesjoint X Y1 .000000 .0000002 40.000000 .0000003 60.000000 .0000004 80.000000 .0000005 120.000000 .000000The Information of SupportsIS VS11 .00000012 .00000013 .00000022 .00000042 .00000052 .000000( NA= 66 )( NW= 299 )Loading Case 1The Loadings at JointsNLJ= 1ILJ PX PY PM3 .0000 -320000.0000 -100000.00000The Loadings at MembersNLM= 4ILM ITL PV DST1 4 -10500.0000 40.0000002 4 -10500.0000 20.0000003 4 -10500.0000 20.0000004 4 -10500.0000 40.000000The Results of CalculationThe Joint Displacementsjoint u v phi1 0.000000E+00 3.713943E-18 4.951923E-172 0.000000E+00 -2.916827E-17 -5.219418E-053 0.000000E+00 -1.405865E-03 1.038816E-064 0.000000E+00 -3.431731E-17 4.276883E-055 0.000000E+00 6.771635E-18 5.239688E-05The Terminal Forcesmember N(st) Q(st) M(st) N(en) Q(en) M(en)1 .000 172860.600 904807.700 .000 247139.400-2390385.0002 .000 359543.300 2390385.000 .000-149543.300 2700481.0003 .000-170456.700-2800481.000 .000 380456.700-2708654.0004 .000 277716.300 2708654.000 .000 142283.700 .000( NA= 66 )( NW= 315 )。
结构力学 矩阵位移法作业
1．作图示刚架的N F 、S F 、M 图，已知各杆截面均为矩形，柱截面宽0.4m,高0.4m, 大跨梁截面宽0.35m,高0.85m ，小跨梁截面宽0.35m,高0.6m ，各杆E=3.0×104 MPa标号及分单元划分计算输出结果************************************************************************* * * wang cheng hao 200901 * * * ************************************************************************The Input DataThe General InformationE NM NJ NS NLC 3.000E+07 10 9 9 1The Information of Membersmember start end A I1 12 1.600000E-01 2.133000E-03 2 2 3 1.600000E-01 2.133000E-03 3 3 6 2.975000E-01 1.790000E-024 25 2.975000E-01 1.790000E-02 5 4 5 1.600000E-01 2.133000E-036 5 6 1.600000E-01 2.133000E-037 58 2.100000E-01 6.300000E-03 8 69 2.100000E-01 6.300000E-03 9 7 8 1.600000E-01 2.133000E-03 10 8 9 1.600000E-01 2.133000E-03The Joint Coordinatesjoint X Y1 .000000 .0000002 .000000 4.5000003 .000000 7.7000004 7.200000 .0000005 7.200000 4.5000006 7.200000 7.7000007 10.000000 .0000008 10.000000 4.5000009 10.000000 7.700000The Information of SupportsIS VS11 .00000012 .00000013 .00000041 .00000042 .00000043 .00000071 .00000072 .00000073 .000000( NA= 180 )( NW= 673 )Loading Case 1The Loadings at JointsNLJ= 1ILJ PX PY PM9 .0000 .0000 -15.00000The Loadings at MembersNLM= 7ILM ITL PV DST1 3 20.0000 4.5000002 3 20.0000 3.2000003 4 -196.0000 7.2000004 4 -36.0000 7.2000007 2 -26.0000 2.7000007 4 -36.0000 3.8000008 4 -196.0000 3.800000The Results of CalculationThe Joint Displacementsjoint u v phi1 3.034552E-21 -7.568894E-20 -7.494658E-212 4.656375E-03 -7.095838E-04 -4.689976E-043 6.414487E-03 -1.138627E-03 -3.252309E-034 4.091261E-21 -1.284215E-19 -9.106066E-215 4.698061E-03 -1.203951E-03 6.981090E-056 6.309349E-03 -1.976041E-03 2.021588E-037 3.774187E-21 -5.368958E-20 -8.623830E-218 4.687578E-03 -5.033399E-04 -9.276283E-059 6.303058E-03 -7.395401E-04 -5.394127E-04The Terminal Forcesmember N(st) Q(st) M(st) N(en) Q(en) M(en)1 756.889 75.346 108.697 -756.889 14.654 27.8582 643.565 -66.328 -84.600 -643.565 130.328 -230.0493 130.328 643.565 230.049 -130.328 767.635 -676.7014 -51.673 113.325 56.742 51.673 145.875 -173.9265 1284.215 40.913 91.061 -1284.215 -40.913 93.0466 1158.135 116.174 146.849 -1158.135 -116.174 224.9087 23.588 -19.795 -65.969 -23.588 182.595 -115.1788 14.154 390.500 451.794 -14.154 354.300 -28.7159 536.896 37.742 86.238 -536.896 -37.742 83.60010 354.300 14.154 31.578 -354.300 -14.154 13.715( NA= 180 )( NW= 701 )2、计算图示桁架各杆的轴力。
矩阵位移法大作业
①
1
2
ql
②
2
3
③
4
q
1① 2 ② 3 ③
4
y M, x
(a)
(b)
第 3 题图
各 杆 EI 、l 相 同，杆长也相同，具体数值可自己给定。
四．采用程序计算图示结构
i
跨长 l（m）
层高 h（m）
集中力（KN）
1
6
7
30
2
10
4
100
3
8
3
50
其他：
柱刚度：EA=105KN，EI=1.5×104KN.m2 梁刚度：EA=106KN，EI=1.0×104KN.m2 支座沉降 C=0.01m
四．采用程序计算图示结构，并作出弯矩图。 已知各杆 E=3.0×106KN/m2，A1=0.16m2，I1=0.012m4，I2=2I1， A2=2A1，I3=3I1，A3=3A1
第 3 题图
20KN
40KN．m
I1，A1 50KN
40KN
15KN/m I3，A3
25KN
I2，A2 40KN．m
4m
4m
3m
3m
五．编写一段程序，实现“将已知支座位移转化为等效节点荷载”。 六．采用程序计算图示结构，并作出最后内力图。已知各杆 E=3.2×106KN/m2，A=0.16m2，I=0.012m4。
36KN
8KN/m
12KN/m
36KN 54KN.m 3m
3m
q=10KN/m
6KN/m
36KN
3m
3m
六．不修改源程序，计算图示结构。
10KN
35KN
6m
15KN
矩阵位移法例题1
50 3 10 15 57 . 5
3 . 891 50 6 . 228 15 79 . 625 57 . 5
2 . 2387 10 6 m 7 2 . 6993 10 m 4 . 2905 10 6 rad
矩 阵 位 移 法（例题）
结构刚度方程为
F K
即
50 202 . 667 3 8 10 15 10 57 . 794 57 . 5 14 . 425 57 . 794 129 . 422 12 . 948 14 . 425 1 12 . 948 2 127 . 306 3
1 (0,0,0)
5m
y
(2)
(1 )
( 2 )
o
x
5m
(0,0,0) 3
2.5m
矩 阵 位 移 法（例题）
单元（1）
0
168 0 0 8 10 168 0 0
0
0 8 . 064 20 . 16 0 8 . 064 20 . 16
(2)
k
(2)
矩 阵 位 移 法（例题）
结构刚度矩阵
168 34 . 667 8 K 10 57 . 794 14 . 425 202 . 667 8 10 57 . 794 14 . 425 57 . 794 8 . 064 121 . 358 20 . 16 7 . 212 20 . 16 7 . 212 67 . 2 60 . 106 14 . 425
矩阵位移法例题
0
2 1 2
0
0
4 1 3
00 2 00 3
0
0
K③
41
3
0
0
0
00 3 000
5 集成总刚度矩阵
第8章矩阵位移法
4 2 2 2
0 1 8 4 0
K 2 2 4 2 4 1
21
2
4
12
2
0
2 1 4 1 4 1 3 0 2 8
1
2
3
6 形成荷载向量
P 60 190 62.5T
2 结点位移编号矩阵 3 各单元旳定位向量
0 0 0 C 0 0 1
0 0 2 0 0 0
2 3T
U1 0 0 0 0 0 1 U2 0 0 1 0 0 2 U3 0 0 2 0 0 0
-90 250
-250 187.5 -112.5
1
2
3
4
第8章矩阵位移法
4 各单元旳刚度矩阵
单元旳刚度矩阵与解法一相同
2 12i 2 BCx l2 Cy
12i (B l2 )CxC y
2 12i 2
BC Y
2 l
Cx
6i l Cy 6i l Cx
2 12i 2 BCx 2 C y
l 12i (B 2 )CxC y l
12i (B 2 )CxC y
l 2 12i 2 BCy 2 Cx
l
6i l Cy 6i l Cx
(e)
K
6i
4i
l Cy
6i l Cx
2i
2 12i 2 BCx 2 C y
l
12i (B 2 )CxC y
l
6i
结构力学习题集(下)矩阵位移法习题及答案
第八章 矩阵位移法一、(O) 二、(X) 3、(O) 4、(X) 五、(X) 六、(O) 7、(O) 八、(X) 九、(O) 10、(O) 1一、(A)一、判定题：1、单元刚度矩阵反映了该单元杆端位移与杆端力之间的关系。
二、单元刚度矩阵均具有对称性和奇异性。
3、局部坐标系与整体坐标系之间的坐标变换矩阵T 是正交矩阵。
4、结构刚度矩阵反映了结构结点位移与荷载之间的关系。
五、结构刚度方程矩阵形式为：[]{}{}K P ∆=，它是整个结构所应知足的变形条件。
6、图示结构用矩阵位移法计算时（计轴向变形）未知量数量为8个。
7、在直接刚度法的先处置法中，定位向量的物理意义是变形持续条件和位移边界条件。
8、等效结点荷载数值等于汇交于该结点所有固端力的代数和。
九、矩阵位移法中，等效结点荷载的“等效原那么”是指与非结点荷载的结点位移相等。
10、矩阵位移法既能计算超静定结构，也能计算静定结构。
11、已知图示刚架各杆EI = 常数，当只考虑弯曲变形，且各杆单元类型相同时，采纳先处置法进行结点位移编号，其正确编号是：(0,1,2)(0,0,0)(0,0,0)(0,1,3)(0,0,0)(1,2,0)(0,0,0)(0,0,3)(1,0,2)(0,0,0)(0,0,0)(1,0,3)(0,0,0)(0,1,2)(0,0,0)(0,3,4)A.B.C.D.2134123412341234( )二、计算题：12、用先处置法计算图示结构刚度矩阵的元素133322,,K K K 。
123ll4l5EI2EIEA(0,0,0)(0,0,1)(0,2,3)(0,0,0)(0,2,4)(0,0,0)EI13、用先处置法计算图示刚架结构刚度矩阵的元素153422,,K K K 。
EI ，EA 均为常数。
l14、计算图示结构整体刚度矩阵的元素665544,,K K K 。
E 为常数。
l l1342A , I AA /222A I , 2A1五、写出图示结构以子矩阵形式表达的结构原始刚度矩阵的子矩阵[][]K K 2224,。
习题课1 矩阵位移法(含答案作业)_518706462
4
5
6
7
8
k
i = 2,3 (1) 54
+ k
i = 2,3 (1) 55
(2) (3) (3) (3) k16 k15 k16 k14 0 (2) (3) (3) (3) k26 k25 k26 k24 0 (2) (3) (3) (3) k36 k34 k35 k36 0
+ k
+
(i ) 33
k
3EIa 2 a 3 + b3
A
3EIab a 3 + b3
B A
3EIab a 3 + b3
3EIb 2 a 3 + b3
B
3EIa a 3 + b3
e θA =1
−3EIa a 3 + b3
3EIb a 3 + b3
e θB =1
−3EIb a 3 + b3
[k ]
e
=
a2 ab
ab b2
e
3EI a 3 + b3
{F }
u2
v2 θ 2 θ 3 ]
−M 0 ]
[0 M 0
0 0 2M 0
T
4
3
3
4
5
0
0
6
2 2 2 2 2 2 k12 k13 k14 k15 k16 k11
2 2 2 2 2 2 k22 k24 k25 k21 k23 k26 2 2 2 2 2 2 k32 k34 k35 k31 k33 k36 2 2 2 2 2 2 k42 k45 k44 k41 k46 k43
y
x
解： T 用位移法求解，未知量为 {∆} = [θ 2 v3 ] 。 1） 杆端弯矩表达式
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
矩阵位移法大作业
学号：151210122
姓名：谭逸天
班级：土木一班
编制原理：
使用Math Work公司开发的科学与工程计算机软件——MATLAB，
利用其矩阵运算的便利性，将题目要求结构的基本信息编入脚本命令文件中，并编入求解步骤。
加上刚度信息的输入指令，以及提取解答要求信息并输出的指令。
令使用者只需输入结构材料相关信息便可计算题目对应悬索—拱组合体系的信息，并直接在命令窗口输出。
利用计算套路的重复性，程序开发时进行模块化设计。
再由重复单元完成多次、重复的运算。
从整体性考虑，数据储存采用“算后集装，装后回收”对变量及数组重复使用，由配音进行简单命名，提高可辨识度。
由于计算套路及程序本身高度模块化，并且题目所需个体信息相对于整体极少，提取个体化的信息只需简单改造命令模块，从整体信息中提取处理得出。
编程所需的“数据化”“编码”等预处理由人工在编程开始前完成，由左下斜索基座作原点，正右向为X轴正向，正上为Y轴正向，建立右手系。
编码顺序从左倒右由上及下，并用先处理法处理基座。
（如下图所示）
6
7
共45个单元，32个结点编号，71个位移编号。
本人学号对应节间数m=14;f1=7L/4;f2=7L/10;h=7L/2;以上数据
为编程中人工设定值，结构的其余信息根据用户的输入进行计算得出。
程序说明：
初始计算结构在坐标系中的坐标信息，手动编入悬索与拱的曲线关键点信息，代入方程求解。
随后由循环语句模块计算并存储结构中各类杆件的角度、长度信息，采用以直代曲的方法处理曲线。
由于先处理法，两端各四个单元不与其余单元通用编码递进规律，采用单独的语句进行计算并集装入总体信息储存矩阵中，其余规律性单元信息由循环的语句模块进行集装，便于之后的计算。
定位向量统一装至71行6列的矩阵“dingwei”中，单元的长度与夹角信息统一装至71行2列的矩阵“danyuan”中，第一列为长度，第二列为角度。
使两个信息矩阵的行序号对应单元序号，便于之后使用。
之后进入单元分析部分。
先是对上部悬索进行单元分析，此部分为桁架单元，从“danyuan”矩阵中提取长度信息与角度信息，结合
开始时输入的刚度信息组装单刚矩阵与坐标变换矩阵，进行坐标变换后直接提取定位向量进行集装部分总刚矩阵的步骤。
集装命令通过循环嵌套配合判断语句，对单刚矩阵进行二维遍历，并提取合格的元素填充至对应位置。
随后，通过少量改动实现对斜索、吊杆、拱、主塔的处理。
之后保留基本结构，进行单元结点荷载的分析，并集装出结构结点荷载矩阵。
之后通过简单矩阵运算即得结构结点位移列阵。
进入单元后处理。
将集装循环语句进行改造，达成逆向提取单元结点位移的功能。
提取之前存储的单元信息进行坐标变换。
最后算出
单元杆端力。
其中认为塔底截面同时包含拱底。
输出位移直接从结点位移矩阵中提取对应位移，杆端力则从后处理内容中提取。
其中加入对塔位移方向描述的判断语句。
源代码：
算例与结果：1
2。