计算结构力学.doc

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

学号: 08201010318
指导: 苏少卿晏燕
安徽建筑工业学院土木工程学院
实验一: 形成单元定位向量的程序设计
求图示结构各单元的单元定位向量。
输入数据文件
8,10,10
9,2,9,4,10,6,10,8,1,3,3,5,5,7,9,10
1,1,1,0,2,1,1,10001,3,10001,0,0,4,3,3,10001,5,10001,0,0,6,
5,5,10001,7,10001,1,0,8,7,1,10001,9,0,0,10001,10,0,0,10001
计算结果:
平面杆系分析

单元数NE= 8
结点数NJ= 10
特殊结点数NJT= 10
单元两端结点号:
单元起点终点单元起点终点单元起点终点单元起点终点
1 9 2 2 9 4 3 10 6 4 10 8
5 1 3 6 3 5 7 5 7 8 9 10
特殊结点信息(结点号,结点各自由度的约束特征数):
结点号特征数特征数特征数结点号特征数特征数特征数结点号特征数特征数特征数
1 1 1 0 2 1 1 10001 3 10001 0 0
4 3 3 10001 5 10001 0 0 6 5 5 10001
7 10001 1 0 8 7 1 10001 9 0 0 10001
10 0 0 10001
未知量编号:
0 0 1 0 0 -1 -1 2 3 -1 2 -1 -1 4 5 -1
4 -1 -1 0 6 -1 0 -1 7 8 -1 9 10 -1
单元定位向量MW(6)
单元编号: 1 MW = 7 8 -1 0 0 -1
单元编号: 2 MW = 7 8 -1 -1 2 -1
单元编号: 3 MW = 9 10 -1 -1 4 -1
单元编号: 4 MW = 9 10 -1 -1 0 -1
单元编号: 5 MW = 0 0 1 -1 2 3
单元编号: 6 MW = -1 2 3 -1 4 5
单元编号: 7 MW = -1 4 5 -1 0 6
单元编号: 8 MW = 7 8 -1 9 10 -1
实验分析:
,其余皆为二力杆;
、截面性质、材料都无关。
实验日期:二〇一一年六月二十一日星期二
实验二: 形成单元刚度矩阵和结构刚度矩阵的程序设计
图示为梁桁,桁式杆的EA=390000,横梁的EI=14040,且不考虑轴向变形。求单元的单元刚度矩阵和结构刚度矩阵。
输入数据文件
8,10,10,2
9,2,9,4,10,6,10,8,1,3,3,5,5,7,9,10
1,1,1,0,2,1,1,10001,3,10001,0,0,4,3,3,10001,5,10001,0,0
6,5,5,10001,7,10001,1,0,8,7,1,10001,9,0,0,10001,10,0,0,10001
1,1,1,1,2,2,2,1
390000,0
0,14040
0,0,2,2,6,6,8,8,2,6
2,2,2,2,2,2,2,2,0,0计算结果:
平面杆系分析
单元数NE= 8
结点数NJ= 10
特殊结点数NJT= 10
截面特性分组数NAI= 2
单元两端结点号:
单元起点终点单元起点终点单元起点终点单元起点终点
1 9 2 2 9 4 3 10 6 4 10 8
5 1 3 6 3 5 7 5 7 8 9 10
单元截面特性:
1 1 1 1 2 2 2 1
特殊结点信息(结点号,结点各自由度的约束特征数):
结点号特征数特征数特征数结点号特征数特征数特征数结点号特征数特征数特征数
1 1 1 0 2 1 1 10001 3 10001 0 0
4 3 3 10001 5 10001 0 0 6 5 5 10001
7 10001 1 0 8 7 1 10001 9 0 0 10001
10 0 0 10001
截面特性EA:
.000
截面特性EI:
.000
结点坐标:
结点 X坐标 Y坐标结点 X坐标 Y坐标结点 X坐标 Y坐标
1 .000 2 .000 3
4 5 6
7 8 9 .000
10