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

《计算结构力学》 计算报告 姓名: 吴亮亮 学号: 09290010402 班级: 09城建土木四班 指导: 苏少卿老师 安徽建筑工业学院土木工程学院 实验一: 形成单元定位向量的程序设计 求图示结构各单元的单元定位向量。 输入数据文件 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 实验分析： 1．单元①是二力杆，是轴力杆单元； 2．单元定位向量与结点坐标和杆件长度、截面性质、材料都无关。 实验日期：二〇一二年六月二十三日星期六  实验二: 形成单元刚度矩阵和结构刚度矩阵的程序设计 图示结构竖杆的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