MATLAB理论力学应用

应用实例为三维空间中的机构受力
代码无误可以复现，其中代码编写参考书籍：《MATLAB理论力学》——程燕平

力学模型的建立与求解

对于本例的力学分析，涉及带轮、轴、链条、桁架、电机、拨叉等重要构件。其中也含有轴承的校核，对于主传动端，即电机（万向连轴）—>轴—>带轮（履带)。

结构示意图（单侧）如下：

图1 主传动端的结构示意图

其中履带的拉力设定为

按具体要求可取定

即

轴的一端连接采用万向联轴器与电机（动力源）结合，整体受到M大小的扭矩，带轮的直径在设计时定为D=210mm，履带的上下侧与铅垂线之间的夹角分别为

和

经过测量得到

从而来计算轴承座的约束力与履带的拉力。

以整个轴为研究对象，轴承A处与万向联轴器B处为一对正交分力，根据静力学定理，可按空间任意力系来构建平衡方程。具体方程组如下

所选联轴器：十字轴式万向联轴器

可以允许两轴间有较大的夹角最大夹角为

且在机器运转时，夹角发生改变仍能正常转动。夹角过大，会使传动效率显著降低。

计算联轴器的计算转矩：由于机器启动时代动载荷和运转中可能出现过载现象，故应当按轴上的最大转矩作为联轴器的计算转矩T_ca

如下式

其中T为公称转矩，单位为N.m

K_A为工作情况系数。

由于原动机为电动机，工作情况为转矩变换和冲击载荷中等的场合。

故可取

即

其中公称转矩

单位为N.mm，且有

其中驱动力矩大概为

用MATLAB进行求解，可以求得：

其中负号表示与图示方向相反。

下述代码是MATLAB语言编写的，作用是求解力学模型 :

 clear;clc;%clean the workspace%理论力学——空间力系（静力学）tictry%specify the symbolic variablessyms  M D a b F1 F2 F_AX F_AZ F_BX F_BZ L_1 L_2%write out the three equations as iseq1 = D*F2/2-D*F1/2+M == 0;eq2 = -L_1*cos(a)*F1-L_1*cos(b)*F2+L_2*F_BZ == 0;eq3 = L_1*sin(a)*F1-L_1*sin(b)*F2+L_2*F_BX == 0;eq4 = F_AZ+F_BZ+cos(a)*F1+cos(b)*F2 ==0;eq5 = F_AX+F_BX+sin(a)*F1+sin(b)*F2 ==0 ;eq6 = F2 == 0.01*F1;%solve the equations symbolicallysol = solve(eq1,eq2,eq3,eq4,eq5,eq6,F1,F2,F_AX,F_AZ,F_BX,F_BZ);%substitute values into the resultsvar = [M D L_1 L_2 a b];val = [18145 210 134 93 85*pi/180 pi/2];F1=subs(sol.F1,var,val);F2=subs(sol.F2,var,val);F_AX = subs(sol.F_AX,var,val);F_AZ = subs(sol.F_AZ,var,val);F_BX = subs(sol.F_BX,var,val);F_BZ = subs(sol.F_BZ,var,val);%使得matlab的计算结果默认显示为小数F1 =  eval(F1);F2 =  eval(F2);F_AX =  eval(F_AX);F_AZ =  eval(F_AZ);F_BX =  eval(F_BX);F_BZ =  eval(F_BZ);disp('该工程问题利用Matlab求解输出的F1为：');disp(F1)disp('该工程问题利用Matlab求解输出的F2为：');disp(F2)disp('该工程问题利用Matlab求解输出的F_AX为：');disp(F_AX)disp('该工程问题利用Matlab求解输出的F_AZ为：');disp(F_AZ)disp('该工程问题利用Matlab求解输出的F_BX为：');disp(F_BX)disp('该工程问题利用Matlab求解输出的F_BZ为：');disp(F_BZ)toccatch %程序异常处理disp('代码运行过程中出现错误, 请检查！')end

MATLAB理论力学应用相关推荐

matlab 理论力学,MATLAB理论力学
绪论静力学静力学引言第一章静力学公理和物体受力分析 1-1 静力学公理 1-2 约束和约束力 1-3 物体的受力分析和受力图习题第二章汇交力系 2-1 汇交力系合成与平衡的几何法(图解法 ...
matlab 碗形状,付星020806-用matlab模拟小球撞碗的运行轨迹_理论力学
用MATLAB 模拟小球撞碗的运行轨迹仪21 付星020806 关键词:斜抛,小球撞碗,MATLAB模拟参考书目:彭芳麟. 理论力学计算机模拟. 清华大学出版社. 2002 <力学与实践&g ...
matlab 开普勒方程,理论力学中非线性问题的 MATLAB数值解
第20卷第10期大学物理Vol.20No.10 2001年10月COLL EGE PHYSICS Oct.2001 基础物理教学现代化问题理论力学中非线性问题的MAT LAB 数值解收稿日期 ...
理论力学类毕业论文文献包含哪些？
本文是为大家整理的理论力学主题相关的10篇毕业论文文献,包括5篇期刊论文和5篇专利论文,为理论力学选题相关人员撰写毕业论文提供参考. 1.[期刊论文]新工科背景下"理论力学"课程建 ...
matlab力学实验,Matlab在力学课程课堂教学和虚拟实验中的应用
龙源期刊网 http://www.doczj.com/doc/31524c5aeffdc8d376eeaeaad1f34693daef10b7.html Matlab在力学课程课堂教学和虚拟实验中的应 ...
[渝粤教育] 四川农业大学理论力学参考资料
教育 -理论力学-章节资料考试资料-四川农业大学[] 第2讲单元测试 1.[单选题]以下说法中错误的是( ) A.理论力学是研究物体机械运动一般规律的科学. B.理论力学与物理中力学部分的主要区别在 ...
[渝粤教育] 西南科技大学理论力学在线考试复习资料
理论力学--在线考试复习资料一.单选题 1.已知W=100kN,P=80kN,摩擦系数f=0.2,物块将( ). A.向上运动 B.向下运动 C.静止不动 2.已知...为作用刚体上的平面共点力系, ...
【渝粤题库】陕西师范大学500001 理论力学作业（专升本）
<理论力学>作业一. 填空在平面极坐标系中,速度的径向分量为______ ,横向分量为_____,加速度的径向分量为_____,横向分量为_____. 在平面自然坐标系中,的方向为__ ...
哈工大理论力学第八版电子版_校史上的这些天（37）| 和你一起在“岁月”中读懂哈工大...
1954年5月26日中央高等教育部委托哈工大组织举办理论力学.机械原理.机械零件3门课程教学座谈会. 1990年5月26至27日哈尔滨建筑工程学院召开第六次党员代表大会.大会选举王玉林.朱华.刘剑 ...
惯性力偶矩公式中j_理论力学(I)习题详解达朗贝尔原理.pdf
您所在位置:网站首页 > 海量文档 &nbsp>&nbsp高等教育&nbsp>&nbsp习题/试题理论力学(I)习题详解达朗贝尔原理.pdf15页 ...

MATLAB理论力学应用

力学模型的建立与求解

MATLAB理论力学应用相关推荐

最新文章

热门文章