计算结构力学程序设计报告材料Word下载.docx
《计算结构力学程序设计报告材料Word下载.docx》由会员分享,可在线阅读,更多相关《计算结构力学程序设计报告材料Word下载.docx(12页珍藏版)》请在冰豆网上搜索。
![计算结构力学程序设计报告材料Word下载.docx](https://file1.bdocx.com/fileroot1/2023-1/28/f2f368e1-f397-4efd-8e7e-f4ddd2d165bd/f2f368e1-f397-4efd-8e7e-f4ddd2d165bd1.gif)
.00.00.00.00.00.00
-37500.00.00100000.0037500.00.0050000.00
-18750.00.0037500.0018750.00.0037500.00
-37500.00.0050000.0037500.00.00100000.00
单元号2
.0018750.0037500.00.00-18750.0037500.00
.0037500.00100000.00.00-37500.0050000.00
.00-18750.00-37500.00.0018750.00-37500.00
.0037500.0050000.00.00-37500.00100000.00
单元号3
18750.00.0037500.00-18750.00.0037500.00
37500.00.00100000.00-37500.00.0050000.00
-18750.00.00-37500.0018750.00.00-37500.00
37500.00.0050000.00-37500.00.00100000.00
结构刚度矩阵:
37500.037500.0.037500.037500.0100000.0
.0.0.0.0100000.050000.0
37500.0.050000.0200000.0
实验二形成单元刚度矩阵和结构刚度矩阵的程序设计
一求图示结构各单元的单元定位向量。
结构竖杆的EA=0EI=14040KN.m其余杆件EA=0EI=8125KN.m
9822
314253647586123456
71118111
111111222
4680039000
140408125
06060606
1010774400
三程序编写过程
!
输入数据说明(数据之间建议用空格作为分隔符)
!
同一项内容的数据可以在同一行输入,也可分多行输入,建议分多行输入,便于阅读
不同的数据项必须分行输入。
1.单元数NE、结点数NJ、特殊结点数NJT、材料特性组数NAI、结点荷载数NPJ、单元荷载数NPM
2.杆端结点号码JH
3.特殊结点信息JTX
4.杆件材料分组信息JMH
5.材料EA
6.材料EI
7.结点X坐标
8.结点Y坐标
9.如果NPJ>
0则读入结点荷载PJZ(NPJ,2)
结点号.方向号(1=X,2=Y,3=M)荷载值
力以沿坐标轴正向为正,弯矩以逆时针方向为正:
整体坐标系下
10.如果NPM>
0则读入单元荷载PMZ(NPM,3)
单元号.荷载类型号(教材P63表中的L值)荷载距单元始结点的距离荷载值
局部坐标系,与单元两端结点编号的顺序有关.
MODULEJGLX
INTEGERNE,NJ,NJT,NJZ,NAI,N,MW(6),NPJ,NPM
定义可分配大小的数组
INTEGER,ALLOCATABLE:
:
JH(:
:
),JW(:
),JZH(:
),JMH(:
),JTX(:
)
DOUBLEPRECISION,ALLOCATABLE:
EA(:
),EI(:
),X(:
),Y(:
),SL(:
),CX(:
),SY(:
),FLZ(:
),P(:
PJZ(:
),PMZ(:
),ZK(:
),XSA(:
ENDMODULEJGLX
形成结构总荷载向量
USEJGLX
CHARACTER*20INAME,ONAME
INTEGER*2ID
CALLGETARG(1,INAME,ID)
IF(ID<
0)THEN
WRITE(*,"
(1X,A\)"
)"
请指定输入数据文件:
"
READ(*,'
(A)'
)INAME
ENDIF
OPEN(1,FILE=INAME,STATUS='
OLD'
ERR=2)!
打开输入数据文件
CALLGETARG(2,ONAME,ID)
请指定输出数据文件:
)ONAME
OPEN(2,FILE=ONAME,STATUS='
REPLACE'
打开输出数据文件
读入单元数、结点数、特殊结点数、材料特性组数、结点荷载数、单元荷载数
READ(1,*)NE,NJ,NJT,NAI,NPJ,NPM
分配数组内存
ALLOCATE(JH(2,NE),JW(3,NJ),JTX(4,NJT),JMH(NE))
ALLOCATE(EA(NAI),EI(NAI),X(NJ),Y(NJ),SL(NE),CX(NE),SY(NE),XSA(NE,7))
READ(1,*)((JH(I,J),I=1,2),J=1,NE)!
读入杆端结点号码
READ(1,*)((JTX(l,J),l=1,4),J=1,NJT)!
读入特殊结点信息
READ(1,*)(JMH(I),I=1,NE)!
读入杆件材料分组信息
READ(1,*)(EA(I),I=1,NAI)!
读入材料EA
READ(1,*)(EI(I),I=1,NAI)!
读入材料EI
READ(1,*)(X(I),I=1,NJ)!
读入结点X坐标
READ(1,*)(Y(I),I=1,NJ)!
读入结点Y坐标
CALLQJW!
各结点求知量在结构位移列阵中的编号JW数组
CALLQJZH
ALLOCATE(FLZ(NJZ,3))
WRITE(2,1)NE,NJ,NJT,NJZ,NAI
WRITE(2,"
(/1X,'
单元两端结点号:
'
/4(4X,'
单元'
2X,'
起点'
2X,'
终点'
))"
(4(4X,I4,I6,I6))"
)(J,(JH(I,J),I=1,2),J=1,NE)
单元截面特性:
/20(1X,I3))"
)JMH
特殊结点信息(结点号,结点各自由度的约束特征数):
&
/3(1X,'
结点号'
特征数'
(3(1X,I6,2X,I6,2X,I6,2X,I6))"
)((JTX(I,J),I=1,4),J=1,NJT)
截面特性EA:
/6(1X,F14.3))"
)EA
截面特性EI:
)EI
结点坐标:
/3(2X,'
结点'
5X,'
X坐标'
5X,'
Y坐标'
(3(1X,I5,2F10.3))"
)(I,X(I),Y(I),I=1,NJ)
未知量编号:
/(16I5))"
)JW
ALLOCATE(ZK(N,N),P(N))!
结构总刚度矩阵与总荷载列阵,N是结构位移未知量总数
(/1X,A)"
单元定位向量MW(6)"
输出标题
DOM=1,NE!
求各单元的单元定位向量
CALLQMW(M)!
求单元M的单元定位向量
(1X,'
单元编号:
I5,4X,'
MW='
6I5)"
)M,(MW(I),I=1,6)
ENDDO
CALLDCH!
计算各单元长度及cosα,sinα
CALLQXS!
确定各单元的七个常数
CALLKJX1!
结构刚度矩阵(满阵)
)'
结构刚度矩阵:
(6F12.1)"
)((ZK(I,J),I=1,N),J=1,N)
(/2X,'
结点荷载数'
单元荷载数'
)"
(6X,I6,6X,I6)"
)NPJ,NPM
ALLOCATE(PJZ(NPJ,2))
ALLOCATE(PMZ(NPM,3))
CALLYDX!
计算结构荷载向量
CLOSE
(1)
CLOSE
(2)
DEALLOCATE(JH,JW,JTX,JMH,EA,EI,X,Y,SL,CX,SY,XSA,ZK,FLZ,P,PJZ,PMZ)!
释放数组内存
1FORMAT(//30X,'
平面杆系分析'
//1X,'
单元数NE='
I11,/1X,'
结点数NJ='
特殊结点数NJT='
I6,&
/1X,'
支座结点数NJZ='
I6,/1X,'
截面特性分组数NAI='
I2)
2END
四计算结果
单元数NE=9
结点数NJ=8
特殊结点数NJT=2
支座结点数NJZ=2
截面特性分组数NAI=2
荷载组数NL=0
131242353464
575686712834
956
001111112
71118111
12345.00013245.000
468000.000390000.000
10.00010.00026.00010.00030.0007.000
46.0007.00050.0004.00066.0004.000
70.0000.00086.0000.000
12345678910111213141516
1718000000
1MW=789123
2MW=101112456
3MW=131415789
4MW=161718101112
5MW=000131415
6MW=000161718
7MW=123456
8MW=789101112
9MW=131415161718
结构刚度矩阵各行带宽
9879879879876543
21
结构刚度矩阵对角元位置
1101825344249586673829097103108112
115117
0.000.000.000.000.000.00
208000.000.00-312000.00-208000.000.00-312000.00
0.004115.000.000.00-4115.000.00
-312000.000.00624000.00312000.000.00312000.00
-208000.000.00312000.00208000.000.00312000.00
0.00-4115.000.000.004115.000.00
-312000.000.00312000.00312000.000.00624000.00
单元号4
-312000.000.00312000.00312000.000.00624000.00
单元号5
87750.000.00-175500.00-87750.000.00-175500.00
0.003086.250.000.00-3086.250.00
-175500.000.00468000.00175500.000.00234000.00
-87750.000.00175500.0087750.000.00175500.00
0.00-3086.250.000.003086.250.00
-175500.000.00234000.00175500.000.00468000.00
单元号6
单元号7
2057.500.000.00-2057.500.000.00
0.0026000.0078000.000.00-26000.0078000.00
0.0078000.00312000.000.00-78000.00156000.00
-2057.500.000.002057.500.000.00
0.00-26000.00-78000.000.0026000.00-78000.00
0.0078000.00156000.000.00-78000.00312000.00
单元号8
单元号9
2207.500.000.00-2207.500.000.00
0.0021666.6765000.000.00-21666.6765000.00
0.0065000.00260000.000.00-65000.00130000.00
-2207.500.000.002207.500.000.00
0.00-21666.67-65000.000.0021666.67-65000.00
0.0065000.00130000.000.00-65000.00260000.00
结构刚度矩阵
2057.5000.0000.000-2057.5000.0000.000
0.0000.0000.00026000.00078000.0000.000
-26000.00078000.0000.0000.0000.000312000.000
0.000-78000.000156000.0000.0000.0000.000
2057.5000.0000.0000.0000.0000.000
0.0000.0000.00026000.000-78000.0000.000
0.0000.0000.0000.0000.000312000.000
0.0000.0000.0000.0000.0000.000
210057.5000.000312000.000-2057.5000.0000.000
-208000.0000.000312000.00030115.00078000.0000.000
-26000.00078000.0000.000-4115.0000.000936000.000
0.000-78000.000156000.000-312000.0000.000312000.000
210057.5000.000312000.0000.0000.0000.000
-208000.0000.000312000.00030115.000-78000.0000.000
0.0000.0000.000-4115.0000.000936000.000
0.0000.0000.000-312000.0000.000312000.000
297957.5000.000-136500.000-2207.5000.0000.000
28867.91765000.0000.000-21666.66765000.0001352000.000
0.000-65000.000130000.000297957.5000.000-136500.000
28867.917-65000.0001352000.000
质量矩阵:
1500.00001500.00004500.00001500.00001500.00004500.0000
2250.00002250.00005062.50002250.00002250.00005062.5000
2650.00002650.00004595.83332650.00002650.00004595.8333
迭代次数=1误差=NaN振型系数R
(1)=NaN
固有频率W
(1)=NaNHz周期T
(1)=NaN秒
振型向量:
NaNNaNNaNNaNNaNNaN
迭代次数=1误差=NaN振型系数R
(2)=NaN
固有频率W
(2)=NaNHz周期T
(2)=NaN秒