结构矩阵分析及程序设计习题
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
如图所示，已知弹模E=2.1×105MPa，各杆截面面积A=10.0 cm2，截面惯性距I=21.0 cm4 ，用平面刚架分析程序计算其内力，要求给出单元划分图，输入数据及分析结果。（注意单位的统一。将分析结果发送到civil2007@邮箱）
答案：
10 9 4
1 0.0 0.0 1 1 1
2 0.0 4.2 0 0 0
3 0.0 8.
4 0 0 0
4 3.6 8.4 0 0 0
5 3.
6 4.2 0 0 0
6 3.6 0.0 1 1 1
7 7.5 8.4 0 0 0
8 7.5 4.2 0 0 0
9 7.5 0.0 1 1 1
1 1
2 3
2 2
3 2
3 3
4 1
4 2
5 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
1 2.1E8 0.001 2.1E-7
2 2.1E8 0.001 4.2E-7
3 2.1E8 0.001 6.3E-7
4 2.1E8 0.001 3.15E-7
2 3
2 1 6
3 1 9
3 2 1.8 -10
7 1 3.9 -3
8 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 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
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-0
2 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 V ALUE
2. 1. 6.0000
3. 1. 9.0000
NON-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-0
4 0.124314E+00
4 0.100549E+01 -0.711760E-03 0.438206E-01
5 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.5760 Q1= 1.523
3 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