弹性梁的大挠度形变问题问题有一根长度为b(m)的均匀的薄弹性梁.PDFVIP

弹性梁的大挠度形变问题 问题: 有一根长度为b(m) 的均匀的薄弹性梁，左端垂直地固接于墙体内，梁因重力而向 下弯曲，设固接处为直角坐标原点，x 轴向右为正，y 轴向下为正. 设梁右端的纵坐标为h1 ， 设地面为y=h0 ， 当h 小于h1 时, 梁的右端将与地面接触. 假定梁在受弯力矩为M (Nm) 时曲率为k=  M，  (N-1m-2)=1/EI, EI 为抗弯刚度. 设该梁单位长度的质量为q (kgm-1) ， 当地重力加速度为g （ms-2 ）, 假定梁的形变都不超出弹性限度. 对于任意给定的b ，，q， g 的具体值后，讨论给定不同的h 值时梁弯曲的形状，当hh1 时, 还要求确定梁对地面的压 力F(N)及确定梁与地面接触的部位的坐标. 理论分析: 以梁的曲线的弧长s 为自变量，x(s)，y(s)是弧长的函数，弧长从坐标原点起 算，故x(0)=0, y(0)=0, 设梁右端的位置为(x(b),y(b)), 记梁的曲线的切线的倾角为z(s), 则k = z(s) =  M, 当梁的右端不与地面接触时, 即hh1 时, 在曲线上点(x(s),y(s)处的弯矩M(s)等 于梁在点(x(s),y(s)) 的右边直到点(x(b), y(b)) 的部分由重力对点(x(s),y(s))产生的力矩., 则 M(s)=qg  int((x(t)-x(s)),t, s,b), 0stb. 记a= qg (m-3) ，则由问题的假设z(s)= a  int((x(t)-x(s)),t, s,b) ，这是一个微分-积分方程，为 了化为微分方程，两边对s 求导得 z(s) = - a (b-s) cos(z(s)), (因x(s)=cos(z(s)) ，y(s)=sin(z(s))), 用MATLAB 的dsolve 求解得知这个方程没有初等积分的解，故只能用数值计算求解，我们 采用MATLAB 的解常微分方程的初值问题的函数ode45 来求解. 因此问题的数学模型为一阶常微分方程组的边值问题 z(s)=w(s), w(s) =a (s-b) cos(z), z(0)=0, w(b)=0, (1) 边值问题的求解可以通过选取适当的初值 w(0)=w . 使得w(b)=0 而惟一得到，这方法称为 0 求边值问题的打靶法. 对于(1)的解z=z(s)，由初值问题y(s)=sin(z(s)), y(0)=0, 积分得y=y(s), 记y(b)=: h . 1 当 hh1 时，假设梁只有右端点与地面接触， 设 F0 是地面对梁的作用力，则对点 (x(s),y(s))产生的总力矩为 M(s)=qg  int((x(t)-x(s)),t, s,b) –F(x(b)-x(s)), 0stb. 这个力F 应不大于梁的重力，故可设 F=(1-p)qgb， 0p 1. 类似地得到微分方程 z(s) = (F - qg (b-s)) cos(z), 记a=qg, 则对应于常微分方程组的边值问题 z(s)=w(s), w(s) =a (s-pb) cos(z), z(0)=0, w(b)=0, y(b)=h (2) 其中p 由条件y(b)=int(sin(z(s)), s, 0, b)=h 确定. 当h 小于另一临界值h h 时，梁上的右边将有一段弧与地面重合，设梁的弧长的[c， 2 1 b]段与地面重合，这时梁的曲线在弧长的c 点是和地面相切的, 总力矩为 M(s)=qg  int((x(t)-x(s)),t, s,c) –int(f(t)(x(t)-x(s)), t, c, b), 0