练习45材料力学-Read.DOC

练习45 材料力学 知识背景 在力学问题中，还有一部分非常重要的内容：材料力学。顾名思义，材料力学的研究对象是各种材料的力学性质。在工程建设中，我们第一个考虑的问题就是 建筑材料的性能。比如要建造大桥，我们就会对指定跨度小备选钢材能否承受负重进行调查。如果不合格，就不能用于建设。材料力学很大程度建立在实验的基础上，要想确切饿了解一个材料的性质，必须反复实验，归纳出他的性质，进而确定它在工程上的应用。 主要内容 【】×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=－Smin=3.8396 Txy=11.1490 改变参数后得到正应力都为正值。 有了这个程序，输入任意参数，我们都可以求出正应力和剪应力，以后遇到这类问题，就可以MATLAB轻松求解。 【】Contrapz主要用来做较为精确的积分运算。在积分时，要注意灵活应用。将数据可视化始终是解决工程实际问题的目标之一，比如绘制莫尔图。就是工程中常见任务。希望读者对照例题，仔细体会程序的意义。材料力学中的问题常常要涉及绘图，希望读者建立起以图形方法解决问题的思想，灵活运用学过的知识，以最优方法解决疑难问题， 【】Axis equsl语句的作用是什么？在什么情况下用到它？ 请自己编制一个用以绘制莫尔图的函数M文件，并进行调用。重新求解本练习中的问题，对照一下结果。 绘图时，hold命令有什么作用？还有没有实现此功能的其他方法？