计算结构力学实验报告.doc

文档介绍：《计算结构力学》
计算报告
姓名:

学号:

班级:
指导:
实验一: 形成单元定位向量的程序设计
求图示结构各单元的单元定位向量。
输入数据文件
9 , 8, 2
5,7,6,8,3,5,4,6,1,3,2,4,2,3,4,5,6,7
1,1,1,1,2,1,1,1
计算结果:

平面杆系分析
单元数NE= 9
结点数NJ= 8
特殊结点数NJT= 2
单元两端结点号:
单元起点终点单元起点终点单元起点终点单元起点终点
1 5 7 2 6 8 3 3 5 4 4 6
5 1 3 6 2 4 7 2 3 8 4 5
9 6 7
特殊结点信息(结点号,结点各自由度的约束特征数):
结点号特征数特征数特征数结点号特征数特征数特征数结点号特征数特征数特征数
1 1 1 1 2 1 1 1
未知量编号:
0 0 0 0 0 0 1 2 3 4 5 6 7 8 9 10
11 12 13 14 15 16 17 18
单元定位向量MW(6)
单元编号: 1 MW = 7 8 9 13 14 15
单元编号: 2 MW = 10 11 12 16 17 18
单元编号: 3 MW = 1 2 3 7 8 9
单元编号: 4 MW = 4 5 6 10 11 12
单元编号: 5 MW = 0 0 0 1 2 3
单元编号: 6 MW = 0 0 0 4 5 6
单元编号: 7 MW = 0 0 0 1 2 3
单元编号: 8 MW = 4 5 6 7 8 9
单元编号: 9 MW = 10 11 12 13 14 15

实验分析:
①是二力杆,是轴力杆单元;
、截面性质、材料都无关。
实验日期:二〇一二年六月二十三日星期六
实验二: 形成单元刚度矩阵和结构刚度矩阵的程序设计
图示结构竖杆的EA=468000KN,EI=14040KNm2,其余杆件EA=390000KN,EI=8125KNm2。求结构的总载荷向量。
输入数据文件
9,8,2,2
5,7,6,8,3,5,4,6,1,3,2,4,2,3,4,5,6,7
1,1,1,1,2,1,1,1
2,1,1,1,1,1,1,2,2
468000,390000
14040,8125
0,6,0,6,0,6,0,6
0,0,4,4,7,7,10,10
计算结果:

平面杆系分析
单元数NE= 9
结点数NJ= 8
特殊结点数NJT= 2
截面特性分组数NAI= 2
单元两端结点号:
单元起点终点单元起点终点单元起点终点单元起点终点
1 5 7 2 6 8 3 3 5 4 4 6
5 1 3 6 2 4 7 2 3 8 4 5
9 6 7
单元截面特性:
2 1 1 1 1 1 1 2 2
特殊结点信息(结点号,结点各自由度的约束特征数):
结点号特征数特征数特征数结点号特征数特征数特征数结点号特征数特征数特征数
1 1 1 1 2 1 1 1
截面特性EA:

截面特性EI:

结点坐标:
结点 X坐标 Y坐标结点 X坐标 Y坐标结点 X坐标 Y坐标
1 .000 .000 2 .000 3 .000
4 5 .000 6
7 .000 8
未知量编号:
0 0 0 0 0 0 1 2 3 4 5 6 7 8 9 10
11 12 13 14 15 16 17 18
单元定位向量MW(6)
单元编号: 1 MW = 7 8 9 13 14 15
单元编号: 2 MW = 10 11 12 16 17 18
单元编号: 3 MW = 1 2 3 7 8 9
单元编号: 4 MW = 4 5 6 10 11 12
单元编号: 5 MW = 0 0 0 1 2 3
单元编号: 6 MW = 0 0 0 4 5 6
单元编号: 7 MW = 0 0 0 1 2 3
单元编号: 8 MW = 4 5 6 7 8 9
单元编号: 9 MW = 10 11 12 13 14 15
单元刚度矩阵:单元号 1
.00 - - .00 -
.00 .00