结构矩阵分析及程序设计习题.docx

1、结构矩阵分析及程序设计习题如图所示，已知弹模E=2.1105MPa，各杆截面面积A=10.0 cm2，截面惯性距I=21.0 cm4 ，用平面刚架分析程序计算其内力，要求给出单元划分图，输入数据及分析结果。（注意单位的统一。将分析结果发送到 civil2007 邮箱） 答案：10 9 41 0.0 0.0 1 1 12 0.0 4.2 0 0 03 0.0 8.4 0 0 04 3.6 8.4 0 0 05 3.6 4.2 0 0 06 3.6 0.0 1 1 17 7.5 8.4 0 0 08 7.5 4.2 0 0 09 7.5 0.0 1 1 11 1 2 3 2 2 3 23 3 4

2、1 4 2 5 4 5 5 6 36 4 5 27 4 7 18 5 8 49 8 9 310 7 8 21 2.1E8 0.001 2.1E-72 2.1E8 0.001 4.2E-73 2.1E8 0.001 6.3E-74 2.1E8 0.001 3.15E-72 32 1 63 1 93 2 1.8 -107 1 3.9 -38 1 3.9 -5运行结果： INPUT FILENAME:11.txt OUTPUT FILENAME:22.txt NUMBER OF ELEMENT= 10 NUMBER OF NODE= 9 NUMBER OF ELEMENT TYPE= 4 NODAL

3、 MESSAGE NODE X Y XX YY ZZ 1 0.0000 0.0000 1 1 1 2 0.0000 4.2000 0 0 0 3 0.0000 8.4000 0 0 0 4 3.6000 8.4000 0 0 0 5 3.6000 4.2000 0 0 0 6 3.6000 0.0000 1 1 1 7 7.5000 8.4000 0 0 0 8 7.5000 4.2000 0 0 0 9 7.5000 0.0000 1 1 1 ELEMENT MESSAGE ELEMENT NODE-I NODE-J TYPE 1 1 2 3 2 2 3 2 3 3 4 1 4 2 5 4

4、5 5 6 3 6 4 5 2 7 4 7 1 8 5 8 4 9 8 9 3 10 7 8 2 ELEMENT TYPE MESSAGE NUMBER E A I 1 0.2100E+09 0.1000E-02 0.2100E-06 2 0.2100E+09 0.1000E-02 0.4200E-06 3 0.2100E+09 0.1000E-02 0.6300E-06 4 0.2100E+09 0.1000E-02 0.3150E-06 NPJ= 2 NPF= 3 NODAL LOADS NODE DIRECTION VALUE 2. 1. 6.0000 3. 1. 9.0000 NON-

5、NODAL LOADS ELEMENT IND a q 3. 2. 1.8000 -10.0000 7. 1. 3.9000 -3.0000 8. 1. 3.9000 -5.0000 MAX.HALF-BAND WIDTH-NW= 12 NODE U V CETA 1 0.000000E+00 0.000000E+00 0.000000E+00 2 0.431838E+00 0.887286E-04 0.981388E-01 3 0.100562E+01 0.572090E-04 0.124314E+00 4 0.100549E+01 -0.711760E-03 0.438206E-01 5

6、0.431792E+00 -0.465157E-03 0.984924E-01 6 0.000000E+00 0.000000E+00 0.000000E+00 7 0.100543E+01 -0.603449E-03 0.678640E-01 8 0.431758E+00 -0.447572E-03 0.867737E-01 9 0.000000E+00 0.000000E+00 0.000000E+00 ELEMENT N Q M 1 N1= -4.4364 Q1= 4.8374 M1= -13.2500 N2= 4.4364 Q2= -4.8374 M2= -7.0672 2 N1= 1

7、.5760 Q1= 1.5233 M1= -3.7486 N2= -1.5760 Q2= -1.5233 M2= -2.6493 3 N1= 7.4739 Q1= 1.5760 M1= 2.6493 N2= -7.4739 Q2= 8.4240 M2= 9.6772 4 N1= 2.6858 Q1= -6.0124 M1= 10.8158 N2= -2.6858 Q2= 6.0124 M2= 10.8288 5 N1= 23.2578 Q1= 4.8205 M1= -7.0206 N2= -23.2578 Q2= -4.8205 M2= -13.2256 6 N1= 12.3301 Q1= 3

8、.9263 M1= -9.3934 N2= -12.3301 Q2= -3.9263 M2= -7.0972 7 N1= 3.5558 Q1= 3.9061 M1= -0.2838 N2= -3.5558 Q2= 7.7939 M2= 7.8649 8 N1= 1.7905 Q1= 4.9153 M1= 3.2889 N2= -1.7905 Q2= 14.5847 M2= 15.5664 9 N1= 22.3786 Q1= 5.3472 M1= -8.4956 N2= -22.3786 Q2= -5.3472 M2= -13.9624 10 N1= 7.7939 Q1= 3.5561 M1=

9、-7.8649 N2= -7.7939 Q2= -3.5561 M2= -7.0707答案：28 20 21 0.0 0.0 1 1 12 0.0 3.3 0 0 03 0.0 6.6 0 0 04 0.0 9.9 0 0 05 0.0 13.2 0 0 06 6.6 13.2 0 0 07 6.6 9.9 0 0 08 6.6 6.6 0 0 09 6.6 3.3 0 0 010 6.6 0.0 1 1 111 9.0 13.2 0 0 012 9.0 9.9 0 0 013 9.0 6.6 0 0 014 9.0 3.3 0 0 015 9.0 0.0 1 1 1 16 15.0 13.2

10、 0 0 017 15.0 9.9 0 0 018 15.0 6.6 0 0 019 15.0 3.3 0 0 020 15.0 0.0 1 1 11 1 2 12 2 3 13 3 4 14 4 5 15 5 6 26 4 7 27 3 8 28 2 9 29 9 10 110 8 9 111 7 8 112 6 7 113 6 11 214 7 12 215 8 13 216 9 14 217 14 15 118 13 14 119 12 13 120 11 12 121 11 16 222 12 17 223 13 18 224 14 19 225 19 20 126 18 19 127

11、 17 18 128 16 17 11 3.0E7 0.25 0.00542 2.6E7 0.18 0.00524 02 1 30.223 1 58.734 1 89.095 1 109.59 运行后结果： INPUT FILENAME:22.txt OUTPUT FILENAME:33.txt NUMBER OF ELEMENT= 28 NUMBER OF NODE= 20 NUMBER OF ELEMENT TYPE= 2 NODAL MESSAGE NODE X Y XX YY ZZ 1 0.0000 0.0000 1 1 1 2 0.0000 3.3000 0 0 0 3 0.0000

12、 6.6000 0 0 0 4 0.0000 9.9000 0 0 0 5 0.0000 13.2000 0 0 0 6 6.6000 13.2000 0 0 0 7 6.6000 9.9000 0 0 0 8 6.6000 6.6000 0 0 0 9 6.6000 3.3000 0 0 0 10 6.6000 0.0000 1 1 1 11 9.0000 13.2000 0 0 0 12 9.0000 9.9000 0 0 0 13 9.0000 6.6000 0 0 0 14 9.0000 3.3000 0 0 0 15 9.0000 0.0000 1 1 1 16 15.0000 13

13、.2000 0 0 0 17 15.0000 9.9000 0 0 0 18 15.0000 6.6000 0 0 0 19 15.0000 3.3000 0 0 0 20 15.0000 0.0000 1 1 1 ELEMENT MESSAGE ELEMENT NODE-I NODE-J TYPE 1 1 2 1 2 2 3 1 3 3 4 1 4 4 5 1 5 5 6 2 6 4 7 2 7 3 8 2 8 2 9 2 9 9 10 1 10 8 9 1 11 7 8 1 12 6 7 1 13 6 11 2 14 7 12 2 15 8 13 2 16 9 14 2 17 14 15

14、1 18 13 14 1 19 12 13 1 20 11 12 1 21 11 16 2 22 12 17 2 23 13 18 2 24 14 19 2 25 19 20 1 26 18 19 1 27 17 18 1 28 16 17 1 ELEMENT TYPE MESSAGE NUMBER E A I 1 0.3000E+08 0.2500E+00 0.5400E-02 2 0.2600E+08 0.1800E+00 0.5200E-02 NPJ= 4 NPF= 0 NODAL LOADS NODE DIRECTION VALUE 2. 1. 30.2200 3. 1. 58.730

15、0 4. 1. 89.0900 5. 1. 109.5900 MAX.HALF-BAND WIDTH-NW= 24 NODE U V CETA 1 0.000000E+00 0.000000E+00 0.000000E+00 2 0.262574E-02 0.412457E-04 0.958458E-03 3 0.638787E-02 0.694240E-04 0.902553E-03 4 0.951571E-02 0.846701E-04 0.651327E-03 5 0.114831E-01 0.902267E-04 0.349485E-03 6 0.113527E-01 0.239451

16、E-03 0.283434E-03 7 0.940898E-02 0.234213E-03 0.485875E-03 8 0.631477E-02 0.202410E-03 0.635406E-03 9 0.261004E-02 0.125802E-03 0.610692E-03 10 0.000000E+00 0.000000E+00 0.000000E+00 11 0.113239E-01 -0.224007E-03 0.276990E-03 12 0.938598E-02 -0.219650E-03 0.477397E-03 13 0.629942E-02 -0.190421E-03 0

17、.623704E-03 14 0.260305E-02 -0.118673E-03 0.597791E-03 15 0.000000E+00 0.000000E+00 0.000000E+00 16 0.113002E-01 -0.105671E-03 0.332895E-03 17 0.936844E-02 -0.992332E-04 0.630449E-03 18 0.628954E-02 -0.814135E-04 0.874927E-03 19 0.258111E-02 -0.483750E-04 0.925637E-03 20 0.000000E+00 0.000000E+00 0.

18、000000E+00 ELEMENT N Q M 1 N1= -93.7403 Q1= 56.4902 M1= -140.2605 N2= 93.7403 Q2= -56.4902 M2= -46.1573 2 N1= -64.0415 Q1= 37.4043 M1= -58.9727 N2= 64.0415 Q2= -37.4043 M2= -64.4616 3 N1= -34.6502 Q1= 30.5060 M1= -38.0020 N2= 34.6502 Q2= -30.5060 M2= -62.6679 4 N1= -12.6287 Q1= 17.0970 M1= -13.3924

19、N2= 12.6287 Q2= -17.0970 M2= -43.0278 5 N1= 92.4920 Q1= -12.6287 M1= 43.0277 N2= -92.4920 Q2= 12.6287 M2= 40.3217 6 N1= 75.6812 Q1= -22.0215 M1= 76.0603 N2= -75.6812 Q2= 22.0215 M2= 69.2819 7 N1= 51.8317 Q1= -29.3912 M1= 102.4636 N2= -51.8317 Q2= 29.3912 M2= 91.5186 8 N1= 11.1342 Q1= -29.6988 M1= 10

20、5.1300 N2= -11.1342 Q2= 29.6988 M2= 90.8821 9 N1= -285.9135 Q1= 86.6812 M1= -113.0445 N2= 285.9135 Q2= -86.6812 M2= -173.0034 10 N1= -174.1095 Q1= 89.1845 M1= -145.9413 N2= 174.1095 Q2= -89.1845 M2= -148.3677 11 N1= -72.2792 Q1= 67.2987 M1= -118.3835 N2= 72.2792 Q2= -67.2987 M2= -103.7024 12 N1= -11

21、.9056 Q1= 36.4775 M1= -70.1259 N2= 11.9056 Q2= -36.4775 M2= -50.2500 13 N1= 56.0143 Q1= -24.5343 M1= 29.8043 N2= -56.0143 Q2= 24.5343 M2= 29.0782 14 N1= 44.8577 Q1= -82.3950 M1= 99.3516 N2= -44.8577 Q2= 82.3950 M2= 98.3965 15 N1= 29.9472 Q1= -131.2215 M1= 158.1250 N2= -29.9472 Q2= 131.2215 M2= 156.8

22、066 16 N1= 13.6370 Q1= -141.5028 M1= 170.5301 N2= -13.6370 Q2= 141.5028 M2= 169.0766 17 N1= 269.7105 Q1= 87.4543 M1= -114.9536 N2= -269.7105 Q2= -87.4543 M2= -173.6458 18 N1= 163.0636 Q1= 90.9280 M1= -148.7591 N2= -163.0636 Q2= -90.9280 M2= -151.3033 19 N1= 66.4301 Q1= 68.6863 M1= -120.5147 N2= -66.

23、4301 Q2= -68.6863 M2= -106.1501 20 N1= 9.9028 Q1= 37.5000 M1= -71.7132 N2= -9.9028 Q2= -37.5000 M2= -52.0368 21 N1= 18.5128 Q1= -14.6316 M1= 42.6350 N2= -18.5128 Q2= 14.6316 M2= 45.1545 22 N1= 13.6746 Q1= -25.8679 M1= 74.1551 N2= -13.6746 Q2= 25.8679 M2= 81.0526 23 N1= 7.7065 Q1= -34.5879 M1= 98.102

24、8 N2= -7.7065 Q2= 34.5879 M2= 109.4246 24 N1= 17.1105 Q1= -34.8559 M1= 97.1803 N2= -17.1105 Q2= 34.8559 M2= 111.9552 25 N1= 109.9432 Q1= 57.0054 M1= -48.6185 N2= -109.9432 Q2= -57.0054 M2= -139.4993 26 N1= 75.0874 Q1= 39.8946 M1= -68.3156 N2= -75.0874 Q2= -39.8946 M2= -63.3368 27 N1= 40.4995 Q1= 32.1883 M1= -65.1124 N2= -40.4995 Q2= -32.1883 M2= -41.1091 28 N1= 14.6316 Q1= 18.5136 M1= -45.1547 N2= -14.6316 Q2= -18.5136 M2= -15.9403