结构矩阵程序设计完成书Word下载.docx
《结构矩阵程序设计完成书Word下载.docx》由会员分享，可在线阅读，更多相关《结构矩阵程序设计完成书Word下载.docx（65页珍藏版）》请在冰豆网上搜索。
1.2E6
（4）直接结点荷载：
编号
i
结点号
mj（i）
XD
qj（i,1）
YD
qj（i,2）
MD
qj（i,3）
10.0
（5）非结点荷载：
单元号
mf（i）
荷载类型
ind（i）
a
aq（i）
B
bq（i）
q1
q1（i）
q2
q2（i）
2.0
-15.0
-10.0
（6）可动支座中的约束条件：
编号i
位移分量号ibd（i）
位移值bd（i）
10
11
5
12
0.01745
根据上面的准，在数据文件fr.txt中将其输入，其填写格式如下：
4,3,0,2,5,2,2
1,0.0,0.0
2,0.0,4.0
3,6.0,4.0
4,6.0,0.0
1,1,2,7.2E6,3.0E5
2,2,3,1.05E7,1.2E6
3,3,4,7.2E6,3.0E5
1,1,10.0,0.0,0.0
2,3,20.0,0.0,0.0
1,1,1,2.0,0.0,-15.0,0.0
2,2,2,0.0,6.0,-10.0,-10.0
1,2,0.0
2,3,0.0
3,10,0.0
4,11,0.0
5,12,0.01745
2.结果输出
程序运行后得到结果文件fw.txt,如下：
planeframestructuralanalysis
****************************
=====
StructuralContro1Data
-----------------------------
nnnenfndndfnpjnpen
430252212
NodalCoordinates
Nodexy
100
204
364
460
ElementInformination
Ele.Nojljreaeinpjal
11272000003000004
2231050000012000006
33472000003000004
NodeLoad
imjxdydmd
111000
232000
ElementLoad
imfindapbqq1q2
11120-150
22206-10-10
BoundaryConditions
iibdbd
120
230
3100
4110
512.01745
outputdata
===============
nodaldisplacement
---------------
nodeno.uvfai
14.18411862990811E-021.81822125813453E-18-1.01387430071926E-18
23.97597437193853E-021.01012292118585E-04-8.85165734292348E-04
33.97454580050996E-02-1.34345625451919E-042.02272900254976E-03
44.5E-19-2.41822125813454E-18.01745
elementno.&
member-enfforce:
eleno.n（l）q（l）m（l）n（r）q（r）m（r）
----------------------
1-181.822125813449-10.000000000000293.8874300719265181.82212581344925.0000000000002-23.8874300719258
225.0000000000052-181.8221258134523.887430071926-25.0000000000052241.822125813451247.04532480877
3241.8221258134545.0000000000001-1247.04532480877-241.82212581345-45.00000000000011067.04532480877
M图（KN*m）
Q图（KN）
N图（KN）
习题2.9
1.数据准备与输入
对结构的结点、单元编号、并取如图所示坐标系。
根据程序中数据输入的顺序进行数据准备和输入（单位kN；
m）。
nn，nf，nd，ndf，ne，npj，npe，其中
结点总数nn=10，固定支座总数nf=4，
可动支座总数nd=0可动支座中的约束总数ndf=0
单元总数ne=12有直接结点荷载的结点数npj=2
有非结点荷载作用的单元数npe=6
（2）结点坐标：
6
7
8
9
14.0
20.0
8.0
（3）各单元式末端结点号及EA、EI的值：
末端号jl（i）
4.2E6
5.5E4
（4）直接接点荷载
编号（i）
结点号mj（i）
XDqj（i,1）
YDqj（i,2）
MDqj（i,3）
25.0
-30.0
40.0
（5）非结点荷载
单元号mf（i）
荷载类型ind（i）
b
7.2
-8.3
-5.2
-50.0
3.6
-16.7
11.1
（6）
-------------------------------------------------------------------------------------------------------
由上面的准备，在数据文件fr.txt中将其输入，其填写格式如下：
10,12,4,0,0,2,6
1,0.0,4.0
2,6.0,8.0
4,14.0,8.0
5,14.0,4.0
6,20.0,4.0
7,0.0,0.0
8,6.0,0.0
9,14.0,0.0
10,20.0,0.0
1,1,2,4.2E6,5.5E4
2,1,3,4.2E6,5.5E4
3,1,7,4.2E6,5.5E4
4,2,3,4.2E6,5.5E4
5,3,8,4.2E6,5.5E4
6,2,4,4.2E6,5.5E4
7,3,5,4.2E6,5.5E4
8,4,5,4.2E6,5.5E4
9,5,9,4.2E6,5.5E4
10,4,6,4.2E6,5.5E4
11,5,6,4.2E6,5.5E4
12,6,10,4.2E6,5.5E4
1,1,25.0,0.0,-30.0
2,6,0.0,0.0,40.0
1,1,2,0.0,7.2,-8.3,-8.3
2,1,7,0.0,7.2,-5.6,-5.6
3,6,2,0.0,8.0,-10.0,-10.0
4,7,1,4.0,0.0,-50.0,0.0
5,10,1,3.6,0.0,-16.7,0.0
6,10,6,3.6,0.0,11.1,0.0
结果输出
程序运行结果文件fw.txt，如下：
10124002618
104
268
4148
5144
6204
700
860
9140
10200
1124200000550007.21110255092798
2134200000550006
3174200000550004
4234200000550004
5384200000550004
6244200000550008
7354200000550008
8454200000550004
9594200000550004
10464200000550007.21110255092798
11564200000550006
126104200000550004
11250-30
260040
11207.2-8.3-8.3
21707.2-5.6-5.6
36208-10-10
47140-500
51013.60-16.70
61063.6011.10
19.89777466232346E-04-3.71795616661629E-051.55515552182729E-04
21.03821933489288E-03-1.59823385641561E-049.58715061721365E-05
39.84626485381254E-04-9.36194978810051E-053.95657894187572E-04
49.93304144363291E-04-9.86950097866573E-05-3.71760105709024E-04
59.99813052217335E-04-6.06138803765017E-05-1.918215853485E-04
61.01507461024866E-03-2.01501549286852E-054.43881541222527E-04
7000
8000
9000
10000
member-endforce:
136.338753751327428.5576412537383-30.57779402051083.9812462486725631.202358746261740.4452027706483
23.60568659576393-4.8799678977379612.4385988915028-3.605686595763934.8799678977379616.8412084949249
339.0385397494716.99957185675229-11.8608048709921-39.038539749471-6.99957185675229-16.1374825560171
469.5140821485833-9.5851176218303915.0481724084485-69.51408214858339.5851176218303923.292298078873
598.30047277505541.993516562875521.45326291932808-98.3004727750554-1.99351656287552-9.42732917083015
623.580475028032341.3437522943616-55.4933751790968-23.580475028032338.656247705638444.7433568242039
7-7.9729475889420223.9064227287341-41.5867694931267.9729475889420226.093577271265950.335387663253
839.985185880663411.5567492658166-25.5876531865904-39.9851858806634-11.5567492658166-20.639343876676
963.644574395326814.2668922988041-31.17133139615-63.6445743953268-14.2668922988041-25.8962377990663
109.267182293031467.77530643934579-19.1557036376135-20.36718229303158.9246935606542123.3925841390978
11-10.6830906219296-2.434188756602471.47528760957310.68309062192962.4341887566024713.1298449300418
1221.157********951.31290013047473.47757093086035-21.157********95-1.3129001304747-8.72917145275916
例题3.1
1.数据准备与输入
对结构的结点、单元编号，并取结构坐标系，根据程序中数据输入的顺序进行数据准备和输入。
nn,nf,nd,ndf,ne,npj,npe,其中
节点总数nn=8固定支座总数nf=0
可动支座总数nd=4可动支座中的约束总数ndf=5
单元总数ne=15有直接结点荷载的结点数npj=5
有非结点荷载作用的单元数npe=0
（2）节点坐标：
坐标x（i）
坐标y（i）
3.0
（3）各端点始末端结点号及EA/EI值：
ea（i）
ei（i）
1.8E8
0.0062
0.003
13
14
15
（4）直接结点荷载：
Xd
qj（i，1）
Yd
qj（i，2）
30.0
-40.0
（5）约束条件
16
19
20
22
根据上面的准备，在数据文件tr3.2.txt中将其输入，其填写格式如下：
8,15,0,4,5,5,0
1,0.0,6.0
2,0.0,3.0
3,4.0,6.0
4,4.0,3.0
5,8.0,6.0
6,8.0,3.0
8,4.0,0.0
1,1,2,1.116E6,1.0E-20
2,2,7,1.116E6,1.0E-20
3,1,3,1.116E6,1.0E-20
4,1,4,5.4E5,1.0E-20
5,2,3,5.4E5,1.0E-20
6,2,4,1.116E6,1.0E-20
7,2,8,5.4E5,1.0E-20
8,4,7,5.4E5,1.0E-20
9,3,4,1.116E6,1.0E-20
10,4,8,1.116E6,1.0E-20
11,3,5,1.116E6,1.0E-20
12,3,6,5.4E5,1.0E-20
13,4,5,5.4E5,1.0E-20
14,4,6,1.116E6,1.0E-20
15,5,6,0.558E6,1.0E-20
1,1,0.0,-30.0,0.0
2,2,30.0,0.0.0,0.0
3,3,0.0,-50.0,0.0
4,5,0.0,-30.0,0.0
5,6,0.0,-40.0,0.0
1,13,0.0
2,16,0.0
3,19,0.0
4,20,0.0
5,22,0.0结果输出
程序运行后得到结果文件tw3.2.txt，打印如下
8150455024
106
203
346
443
586
683
700
840
11211160001E-203
22711160001E-203
31311160001E-204
4145400001E-205
5235400001E-205
62411160001E-204
7285400001E-205
8475400001E-205
93411160001E-203
10481116000
升级会员 