练习45材料力学.docVIP

练习45 材料力学 知识背景 在力学问题中，还有一部分非常重要的内容：材料力学。顾名思义，材料力学的研究对象是各种材料的力学性质。在工程建设中，我们第一个考虑的问题就是 建筑材料的性能。比如要建造大桥，我们就会对指定跨度小备选钢材能否承受负重进行调查。如果不合格，就不能用于建设。材料力学很大程度建立在实验的基础上，要想确切饿了解一个材料的性质，必须反复实验，归纳出他的性质，进而确定它在工程上的应用。 主要内容 【本练习讲述知识点】 本练习我们将着手解决一些常见材料力学问题。我们将用到材料力学中的正应力、剪力、弯矩，挠度等概念。希望读者首先查阅相关书籍，做好知识准备，本联系将介绍MATLAB中一个积分函数contrapz的用法，请注意他的用途。 练习过程 （！）悬臂梁 在材料力学中，悬臂梁是一个非常基本的概念。有关悬臂梁的计算也横多。我们来看这样的问题。有一个悬臂梁，入土45-1所示。 悬臂梁长为L。左端插入墙内，离固定端L1处受集中力P，求转角和挠度。 先来建模型，设弯矩为M,转角为A，挠度为Y。取E=150×109 N/m2,I=1.5×10-5m4,l列出力学方程如下： （） =0 显然需要做积分。我们使用MATLAB中的cumtrapz函数做积分。首先来简单介绍一下cumtrapz函数。 我们在命令区里输入： Y=[0 1 2;3 4 5] a=cumtrapz(Y,1) b=cumtrapz(Y,2) 得到的结果如图45-2所示。 从运行结果我们可以看出，cumtrapz(Y,1)与cumtrapz(Y,2)的计算结果不同，请注意各自的用法。 下面我们就用cumtrapz函数来解上面的问题。 取L=2m, P=1500N, L1 =1m, 在命令区里输入： L=2;P=1500;L1=1; E=150e9;I=1.5e-5; x=linspace(0,L,101);dx=L/100 n1=L1/dx+1; M1=-P*(L1-x(1:n1)); M2=zeros(1,101-n1); M=[M1,M2]; A=cumtrapz(M)*dx/(E*I); Y=cumtrapz(A)*dx; subplot(3,1,1),plot(x,M),grid subplot(3,1,2),plot(x,A),grid subplot(3,1,3),plot(x,Y),grid gtext(M) gtext(A) gtext(Y) 得到图45-3。 从图中可以明显地看出悬臂梁弯矩、转角、挠度沿梁分布关系。在工程上，经常要绘制这样的关系图。 （2）绘制莫尔图 已知两正交截面上的正应力、和剪应力，求其它斜截面的应力值。设斜截面与上的正应力和剪应力的计算公式是： 我们在命令区里输入程序： Sx=20; Sy=0; Txy=10; a=linspace(0,pi,36); Sa=(Sx+Sy)/2;Sd=(Sx-Sy)/2; s=Sa+Sd*cos(2*a)-Txy*sin(2*a); t=Sd*sin(2*a)+Txy*cos(2*a); plot(s,t,Sx,Txy); axis equal v=axis; line([v(1),v(2)],[0,0]) line([0,0],[(3),v(4)]) hold,plot(Sa,0,x) Smax=max(s),Smin=min(s),Tmax=max(t) 用gtext命令标注后，得到图45-4所示的莫尔图。 运算结果为： Smax=24.1101 Smin=－4.1386 Tmax=14.1101 若取Sx=20,Sy=20,Txy=10,运行上面的程序，得到图45-5。 得到的结果是： Smax=26.1753 Smin=3.8396 Txy=11.1490 改变参数后得到正应力都为正值。 有了这个程序，输入任意参数，我们都可以求出正应力和剪应力，以后遇到这类问题，就可以MATLAB轻松求解。 【练习小结】 本练习中我们遇到了材料力学中常见的问题：求悬臂梁转角和挠度；绘制莫尔图。通过联系我们知道了contrapz函数的用法。Contrapz主要用来做较为精确的积分运算。在积分时，要注意灵活应用。将数据可视化始终是解决工程实际问题的目标之一，比如绘制莫尔图。就是工程中常见任务。希望读者对照例题，仔细体会程序的意义。材料力学中的问题常常要涉及绘图，希望读者建立起以图形方法解决问题的思想，灵活运用学过的知识，以最优方法解决疑难问题， 【思考题】 在求解悬臂梁的程序中，如果不用cont