关于材料力学中平截面假定的检验.docVIP

材料力学中平截面假定的检验 工程力学 林彦英 201120105181 定义 平截面假定：垂直于杆件轴线的各平截面（即杆的横截面）在杆件受拉伸、压缩或纯弯曲而变形后仍然为平面，并且同变形后的杆件轴线垂直。 目的 检验平截面假定是否合理。 对于图1所示的受力梁，用弹性力学理论知识，取以下五个截面 作检验。 图1梁的受力简图 方法 分别用材料力学和弹性力学知识，计算梁截面变形后的水平和竖直位移公式； 用matlab软件画出位移图像； 分析两种计算方法之间的偏差。 设定弹性模量，剪切模量，泊松比，长度，高度，宽度的值并计算出惯性矩的值分别如下： 分析 对于图1所示的受力梁，在材料力学中，根据平截面假定，我们可以得到挠曲线的近似微分方程 对于图1所示的梁，有 且有边界条件，当时，有，. 由此可以求得挠度方程为 挠度值也就是竖直方向上的位移，即有 （1） 图2 梁的变形示意图 如图2所示，又有几何关系有 因为是小变形，有 又 所以有， 因此，水平方向的位移 （2） 由弹性力学可以推导出 水平方向位移为 （3） 竖直方向位移为 （4） 将各参数的值分别代入公式（1）（2）（3）（4）中，用matlab软件画出位移图像，图像中，蓝色线（或）为弹性力学的结果，红色线（或）为材料力学的结果。横坐标为高度值（即y值）的变化，纵坐标为位移值。 截面 1.水平位移图像 程序 x=0; y=(-0.1:0.01:0.1) u1=-(10^4*x^2*y)/(2*208*10^9*1/(3*10^3))-(0.3*10^4*y.^3)/(6*208*10^9*1/(3*10^3))+(10^4*y.^3)/(6*80*10^9*1/(3*10^3))+((10^4*2^2)/(2*208*10^9*1/(3*10^3))-(10^4*0.1^2)/(2*80*10^9*1/(3*10^3)))*y; u2=-(10^4*x^2*y)/(2*208*10^9*1/(3*10^3))+(10^4*2^2*y)/(2*208*10^9*1/(3*10^3)); plot(y,u1,b,y,u2,r:); 其它截面，只需将第一个程序段改为相应的值即可。 图3 截面上的水平位移 图4 截面上的水平位移 图5 截面上的水平位移 图6 截面上的水平位移 图7 截面上的水平位移 截面 2.竖直位移图像 程序 x=0; y=(-0.1:0.001:0.1); v1=(0.3*10^4*x*y.^2)/(2*208*10^9*1/(3*10^3))+(10^4*x^3)/(6*208*10^9*1/(3*10^3))-(10^4*2^2*x)/(2*208*10^9*1/(3*10^3))+(10^4*2^3)/(3*208*10^9*1/(3*10^3)); v2=(10^4*x^3)/(6*208*10^9*1/(3*10^3))-(10^4*2^2*x)/(2*208*10^9*1/(3*10^3))+(10^4*2^3)/(3*208*10^9*1/(3*10^3)); plot(y,v1,b,y,v2,r:) 其它截面，只需将第一个程序段改为相应的值即可。 图8 截面上的竖直位移 图9 截面上的竖直位移 图10 截面上的竖直位移 图11 截面上的竖直位移 图12 截面上的竖直位移 从以上位移图像看出，弹性力学和材料力学所计算得出的结果相差很小，最大的偏差数量级也达到，在五个截面上，随着值的增大，水平位移偏差越来越大，竖直位移偏差越来越小。 3.绝对误差（偏差） 直接将公式（1）和（4）作差，（2）和（3）作差，即可得到位移偏差的公式 水平方向的位移偏差公式 竖直方向的位移偏差公式 4.相对误差 水平方向位移的相对误差，竖直方向位移的相对误差为. 用matlab软件绘制出相对误差图像 (由于误差非常小，为了能从图像中显示出来，所以将其扩大10^6倍) 4.1水平方向位移相对误差计算 x=0; y=(-0.1:0.0001:0.1); a=((-(0.3*10^4*y.^3)/(6*208*10^9*1/(3*10^3))+(10^4*y.^3)/(6*80*10^9*1/(3*10^3))+(-(10^4*0.1^2)/(2*80*10^9*1/(3*10^3)))*y)/