TLAB及其在理工课程中的应用指南 第三版 十二五陈怀琛 第5－9章 第7章.pptVIP

　　◆ 程序运行结果 　　运行程序所得的结果如图7-12所示。注意几根曲线之间的积分关系。本题之所以简单，是因为在x=0处，转角和挠度都为零，因此两次积分的积分常数恰好都为零。如果它们不为零，程序中就得有确定积分常数的语句，这可在例7-2-3中看到。 图 7-12 悬臂梁弯矩(M)、转角（A）和挠度（Y）曲线 　　【例 7-2-3】 简支梁左半部分受均匀分布载荷q作用，右边L/4处受集中力偶M0作用(如图7-13所示)， 求其弯矩、转角和挠度。设L=2m，q=1000 N/m ，M0＝900 N·m，E=200×109N/m2，I=2×10-6m4。 图 7-13 简支梁受力图 　　解： 　　◆ 建模 　　此题解法基本上与例7-2-2相同，主要差别是要处理积分常数问题。 　　支撑反力Na和Nb可由平衡方程求得，设　　　　，则 各段弯矩方程为: 　　对M/EI作积分，得转角A，再作一次积分，得到挠度Y，每次积分都要出现一个待定积分常数 此处设A0(x)=cumtrapz(M)*dx/EI。 此处设Y0(x)=cumtrapz(A0)*dx。 　　两个待定积分常数Ca和Cy可由边界条件Y(0)=0及Y(L)=0确定: Y(0)=Y0(0)+Cy=0 Y(L)=Y0(L)+Ca·L+Cy=0 于是可得 即 　　◆ MATLAB程序 　%输入已知参数L，q，M0，E，I后，先求两铰链的支撑反力Na和Nb 　L=2; q=1000; M0=900; E=200e9; I=2e-6; 　Na=(3*q*L^2/8-M0)/L; Nb=(q*L^2/8+M0)/L; %求支撑反力 　x=linspace(0，L，101); dx=L/100; %将x分为100小段 　M1=Na*x(1∶51)-q*x(1∶51).^2/2; %分三段用数组列出M的表达式 　M2=Nb*(L-x(52∶76))-M0; 　M3=Nb*(L-x(77∶101)); M=［M1，M2，M3］; %列写完整的M数组 　A0=cumtrapz(M)*dx/(E*J); %由M积分求转角（未计积分常数） 　　Y0=cumtrapz(A0)*dx; %由转角积分求挠度（未计积分常数） 　C=［0，1; L，1］＼［-Y0(1); -Y0(101)］; %由边界条件求积分常数Ca，Cy 　Ca=C(1)，Cy=C(2)， 　A=A0+Ca; Y=Y0+Ca*x+Cy; %求出转角与挠度的完整值 　subplot(3，1，1)，plot(x，M)，grid %绘图 　subplot(3，1，2)，plot(x，A)，grid 　subplot(3，1，3)，plot(x，Y)，grid 　　◆ 程序运行结果 　　执行本程序的结果如图7-14所示。 　　梯形积分累加函数cumtrapz与定积分函数trapz的不同在于cumtrapz类似于不定积分，逐点给出积分的值，因而得出一个数列，而trapz只给出积分到终点的一个值。这些函数都假定步长为1，因此累加的值必须乘以dx才与积分等价。 　　用A=cumtrapz(M)来求面积，长度M为101，只能形成100个A。而cumsum则是把101个点逐个相加，相当于多算了一个点。准确地说，可以推导出 cumtrapz(M) =cumsum(M) - M(1)/2 -M/2 实际上只要点取得足够多，直接用cumsum(M)代替cumtrapz(M)，在工程上也是可以接受的。 图 7-14 例7-2-3的弯矩、转角和挠度曲线 (a) 弯矩曲线; (b) 转角曲线; (c) 挠度曲线 　　【例 7-2-4】 拉弯合成部件的截面设计。 这一设计计算将归结为解一个三次代数方程,过去要用试凑法反复运算,本例显示了用MATLAB求解的简洁。钻床立柱如图7-15所示。设P＝15 kN， 许用拉应力[σ]=35 MPa，钻头轴与立柱轴距离为0.4 m，试求立柱直径。 　　 图7-15 钻床受力图 　　解： 　　◆ 建模 　　立柱受到拉力P和弯矩Pl作用，两者产生的拉应力之和为最大拉应力，令它小于［σ］，即 把　　　　　　 　　代入上式后，得到求直径d的方程 这个三次代数方程可用MATLAB多项式求根的roots函数求解。 　　◆ MATLAB程序 　P=input(′P=′)，l=input(′l=′)， %输入力和偏心距 　asigma=input(′［σ］=′)， 　 %输入许用拉应力