计算结构力学实验报告

《计算结构力学》

实验报告

姓名:彭鹏

学号:08201010534

班级:08土木5班

指导老师: 苏少卿 晏燕

安徽建筑工业学院土木工程学院

实验一: 形成单元定位向量的程序设计

一、 求图示结构各单元的单元定位向量。

8 9

5 7 10

6

4

1 2 3

二、 输入数据文件

8 10 5

8 9 5 7 6 10 5 4 4 6 1 5 2 4 3 6

1 1 1 1 2 1 1 1 3 1 1 1 8 7 7 10001 9 10 10 10001

三、 计算结果：

平面杆系分析

单元数NE= 8

结点数NJ= 10

特殊结点数NJT= 5

单元两端结点号:

单元 起点 终点 单元 起点 终点 单元 起点 终点 单元 起点 终点

1 8 9 2 5 7 3 6 10 4 5 4

5 4 6 6 1 5 7 2 4 8 3 6

特殊结点信息(结点号,结点各自由度的约束特征数):

结点号 特征数 特征数 特征数 结点号 特征数 特征数 特征数 结点号 特征数 特征数

1 1 1 1 2 1 1 1 3 1 1 8 7 7 10001 9 10 10 10001

未知量编号:

0 0 0 0 0 0 0 0 0 1 2 3 4 5 6 7

特征数 1

8 9 10 11 12 10 11 -1 13 14 -1 13 14 15

单元定位向量MW(6)

单元编号: 1 MW = 10 11 -1 13 14 -1

单元编号: 2 MW = 4 5 6 10 11 12

单元编号: 3 MW = 7 8 9 13 14 15

单元编号: 4 MW = 4 5 6 1 2 3

单元编号: 5 MW = 1 2 3 7 8 9

单元编号: 6 MW = 0 0 0 4 5 6

单元编号: 7 MW = 0 0 0 1 2 3

单元编号: 8 MW = 0 0 0 7 8 9

四、 实验分析：

1．单元①是二力杆，是轴力杆单元；

2．单元定位向量与结点坐标和杆件长度、截面性质、材料都无关。

五、 实验日期：二〇〇八年三月十四日星期五

实验二: 形成单元刚度矩阵和结构刚度矩阵的程序设计

一、 图示结构竖杆的EA=220KN,EI=550KNm2,其余杆件

EA=150KN,EI=350KNm2。求各单元的单元刚度矩阵和结构刚度矩阵。

8 9

5 7 10

6

4

1 2 3

二、 输入数据文件

8 10 5 2

8 9 5 7 6 10 5 4 4 6 1 5 2 4 3 6

1 1 1 1 2 1 1 1 3 1 1 1 8 7 7 10001 9 10 10 10001

1 1 1 1 1 2 2 2

150 220

350 550

0 2 5 2 0 5 1.5 1.5 4 4

0 0 0 2 4 4 5 5 5 5

三、 计算结果：

平面杆系分析

单元数NE= 8

结点数NJ= 10

特殊结点数NJT= 5

截面特性分组数NAI= 2

单元两端结点号:

单元 起点 终点 单元 起点 终点 单元 起点 终点

单元 起点 终点

1 8 9 2 5 7 3 6 10 4 5 4

5 4 6 6 1 5 7 2 4 8 3 6

单元截面特性:

1 1 1 1 1 2 2 2

特殊结点信息(结点号,结点各自由度的约束特征数):

结点号 特征数 特征数 特征数 结点号 特征数 特征数 特征数 结点号 特征数 特征数 特征数 1 1 1 1 2 1 1 1 3 1 1 1 8 7 7 10001 9 10 10 10001

截面特性EA:

150.000 220.000

截面特性EI:

350.000 550.000

结点坐标:

结点 X坐标 Y坐标 结点 X坐标 Y坐标 结点 X坐标 Y坐标

1 .000 .000 2 2.000 .000 3 5.000 .000

4 2.000 2.000 5 .000 4.000 6 5.000 4.000

7 1.500 5.000 8 1.500 5.000 9 4.000 5.000

10 4.000 5.000

未知量编号:

0 0 0 0 0 0 0 0 0 1 2 3 4 5 6 7

8 9 10 11 12 10 11 -1 13 14 -1 13 14 15

单元定位向量MW(6)

单元编号: 1 MW = 10 11 -1 13 14 -1

单元编号: 2 MW = 4 5 6 10 11 12

单元编号: 3 MW = 7 8 9 13 14 15

单元编号: 4 MW = 4 5 6 1 2 3

单元编号: 5 MW = 1 2 3 7 8 9

单元编号: 6 MW = 0 0 0 4 5 6

单元编号: 7 MW = 0 0 0 1 2 3

单元编号: 8 MW = 0 0 0 7 8 9

单元刚度矩阵：单元号 1

60.00 .00 .00 -60.00 .00 .00

.00 268.80 336.00 .00 -268.80 336.00

.00 336.00 560.00 .00 -336.00 280.00

-60.00 .00 .00 60.00 .00 .00

.00 -268.80 -336.00 .00 268.80 -336.00

.00 336.00 280.00 .00 -336.00 560.00

单元刚度矩阵：单元号 2

278.17 -292.45 -358.42 -278.17 292.45 -358.42

-292.45 521.88 537.63 292.45 -521.88 537.63

-358.42 537.63 776.58 358.42 -537.63 388.29

-278.17 292.45 358.42 278.17 -292.45 358.42

292.45 -521.88 -537.63 -292.45 521.88 -537.63

-358.42 537.63 388.29 358.42 -537.63 776.58

单元刚度矩阵：单元号 3

795.50 689.43 -742.46 -795.50 -689.43 -742.46

689.43 795.50 -742.46 -689.43 -795.50 -742.46

-742.46 -742.46 989.95 742.46 742.46 494.97

-795.50 -689.43 742.46 79 …… 此处隐藏：1988字，全部文档内容请下载后查看。喜欢就下载吧 ……