创新实践实习报告石家庄铁道大学.docVIP

石家庄铁道大学 创新实践实习报告 学院: 班级: 姓名: 学号: 指导老师: 完成时期 2012年 7月 2日 题目:两跨连续钢梁受 10KN/M的均布力, 已知截面尺寸 0.1m*0.1m,每跨长度 10m, 钢的密度 :每立方米 7800kg, 左跨弹性模量 E=210GPa, 右跨弹性模量 E=420GPa. 求跨中挠度。 一、理论求解 根据结构力学求解超静定结构内力的力法和求解位移的图乘法,可以求得跨中挠 度,解题如下: (1求原结构弯矩图 基本结构如图为两跨简支梁。 基本结构弯矩图如下: 基本结构受均布力的弯矩图如下: p1= 0.5?10EI EI 51233 21=???? p2=EI EI 625132211251032=?????- 带入公式 p1?M+p2=0得 M=KNm p p 1251 2=- 叠加得出原结构弯矩图如下 (2图乘法求跨中挠度 在所求位置和方向加上单位一的力,弯矩如图: 上图和原结构弯矩图相乘即可得出跨中挠度,但是原结构弯矩图比较复杂,因 此可以将原结构弯矩图分解为两个简单的弯矩图,再分别和上图相乘,最后叠加即可 得出跨中挠度。 2?(32?125000?5?85? EI 5. 2=EI 12125000125? -21?10?2.5? 21?EI 125000= -EI 412500025? 左跨跨中挠度 EI 12125000125? - EI 412500025?=0.2976m 同理可得右跨的跨中挠度为 0.1488m. 二、利用有限元法和 matlab 求解跨中挠度 (1计算总刚 由于钢梁截面尺寸较小,跨度较大,所以可以利用“平面刚架的有限元法”求解跨中 挠度。 首先将原结构划分为四个杆单元如图: 本题不用考虑坐标变换,所以杆件 1、 2的单元刚度矩阵相同;杆件 3、 4的单元刚度 矩阵相同。可以首先建立求解单元刚度矩阵的方程: function k=dygd(E,A,I,L k=E*I/(L^3*[A*L*L/I 0 0 -A*L*L/I 0 0 0 12 6*L 0 -12 6*L 0 6*L 4*L*L 0 -6*L 2*L*L -A*L*L/I 0 0 A*L*L/I 0 0 0 -12 -6*L 0 12 -6*L 0 6*L 2*L*L 0 -6*L 4*L*L] 1、由程序语言 k=dygd(210000000000,0.01,0.0001/12,5得出杆件 1、 2的单元刚度矩 阵 k = 420000000 0 0 -420000000 0 0 0 168000 420000 0 -168000 420000 0 420000 1400000 0 -420000 700000 -420000000 0 0 420000000 0 0 0 -168000 -420000 0 168000 -420000 0 420000 700000 0 -420000 1400000 2、由程序语言 K=dygd(420000000000,0.01,0.0001/12,5得出杆件 3、 4的单元刚度矩 阵 k = 420000000 0 0 -420000000 0 0 0 168000 420000 0 -168000 420000 0 420000 1400000 0 -420000 700000 -420000000 0 0 420000000 0 0 0 -168000 -420000 0 168000 -420000 0 420000 700000 0 -420000 1400000 再将四个单元刚度矩阵扩阶成 15 15 k1=zeros(15,15 k1(1:6,1:6=k k2=zeros(15,15 k2(4:9,4:9=k k3=zeros(15,15 k3(7:12,7:12=K k4=zeros(15,15 k4(10:15,10:15=K 最后形成总刚:kk=k1+k2+k3+k4 总刚划行划列:kkk=kk([3 5 6 11 12 15],[3 5 6 11 12 15] 结果为 KKK= 1400000 -420000 700000 0 0 0 -420000 336000 0 0 0 0 700000 0 2800000 0 0 0 0 0 0 672000 0 840000 0 0 0 0 5600000 1400000 0 0 0 840000 1400000 2800000 (2节点力等效 共有两种基本单元类型:两端固定梁和一端固定一端铰支梁 单元 1 单元 2 单元 3、 4同理可得 最后得出综合结点向量为 f=?????? ?????????????????????-??-???-02