[法律资料]A20004051 中南大学 辛婧华 高乐 甘磊.doc
《[法律资料]A20004051 中南大学 辛婧华 高乐 甘磊.doc》由会员分享,可在线阅读,更多相关《[法律资料]A20004051 中南大学 辛婧华 高乐 甘磊.doc(38页珍藏版)》请在三一文库上搜索。
1、2009高教社杯全国大学生数学建模竞赛承 诺 书我们仔细阅读了中国大学生数学建模竞赛的竞赛规则.我们完全明白,在竞赛开始后参赛队员不能以任何方式(包括电话、电子邮件、网上咨询等)与队外的任何人(包括指导教师)研究、讨论与赛题有关的问题。我们知道,抄袭别人的成果是违反竞赛规则的, 如果引用别人的成果或其他公开的资料(包括网上查到的资料),必须按照规定的参考文献的表述方式在正文引用处和参考文献中明确列出。我们郑重承诺,严格遵守竞赛规则,以保证竞赛的公正、公平性。如有违反竞赛规则的行为,我们将受到严肃处理。我们参赛选择的题号是(从A/B/C/D中选择一项填写): A 我们的参赛报名号为(如果赛区设置
2、报名号的话): 20004051 所属学校(请填写完整的全名): 中南大学 参赛队员 (打印并签名) :1. 高乐 2. 辛婧华 3. 甘磊 指导教师或指导教师组负责人 (打印并签名): 刘诚 日期: 2009 年 09 月 14 日赛区评阅编号(由赛区组委会评阅前进行编号):2009高教社杯全国大学生数学建模竞赛编 号 专 用 页赛区评阅编号(由赛区组委会评阅前进行编号):赛区评阅记录(可供赛区评阅时使用):评阅人评分备注全国统一编号(由赛区组委会送交全国前编号):全国评阅编号(由全国组委会评阅前进行编号):制动器试验台的控制方法分析摘要本文通过研究驱动电流与可观测扭矩、转速的关系,设计了制
3、动器试验台控制方法,属于机理型建模。针对问题一,我们依据能量守恒定律和刚体定轴转动的动能定理,建立了载荷存在时的平动能量与等效转动惯量的等价式,求解即得该车轮制动时等效转动惯量为52 kgm2;针对问题二,由转动惯量的积分定义式可计算出三个飞轮各自的转动惯量,任意选择组合可得到的八种机械惯量(见表1),按求得的机械惯量取值区间挑选得电动机补偿的惯量为12 kgm2或者-18 kgm2两种。针对问题三,由扭矩合成原理可得到主轴扭矩、制动扭矩和补偿扭矩三者的关系,再根据刚体定轴转动定理和动能定理,推导出补偿扭矩依赖于主轴扭矩(可观测)的函数,最后引用题目给出电流正比于补偿扭矩,即得电流依赖模型。同
4、时简单计算得已知初速度、且减速度为定值的制动过程中驱动电流大小为174.83(A)或者-262.24(A)。两外对正负电流的意义给出了解释:正值表示按电动机工作,负值表示按发电机工作。针对问题四,路试耗能可由刚体定轴转动定理得到;台试耗能的计算采用“全局连续问题离散化、局部离散问题连续化”方法,通过积分求得耗能。得到耗能相对误差为,进而据此评价该控制方法基本可控,但仍有待改进。针对问题五,考虑实际控制过程中尽量保证瞬时补偿惯量靠近电动机补偿惯量的原则,建立由其差值决定驱动扭矩的波动差模型,再基于第二步提出的电流依赖模型即得驱动电流值。将该法在每个小时间区段上重复使用即为计算机控制方法(见图3)
5、。最后计算该方法的残差,得知在可接受域内波动,且利用它可很好的解释第四问中主轴扭矩随时间变化趋势,故认为模型具有较好的实用性。针对问题六,考虑到问题五中所设计的方法未涉及实际惯性试验台的滞后、时变、非线性性等,为消除这些不足,重新建立了模糊控制系统,而模糊控制系统的鲁棒性强,干扰和参数变化对控制效果的影响可被大大减弱。针对带有模糊控制的电惯量系统模型通过Matlab/Simulink 工具对该系统进行仿真模拟。根据仿真结果(见图15),可看出所设计的计算机控制方法能够获得较好的效果。关键字:扭矩相对能力量误差波动差 模糊控制系统MTLAB仿真一、问题重述行车制动器的设计是车辆设计中最重要的环节
6、之一,而检验其性能的路试方法在设计阶段是无法实现的,故只能利用试验台研究。因此,准确计算设计制动器试验台的计算机控制方法,使其制动过程与路试的制动过程尽量一致显得至关重要。题目首先要求能够利用物理学中转动惯量定义计算公式、刚体定轴转动定理、转动动能定理和能量守恒定律,根据已知半径、载荷、初末速度等量,通过简单推理和计算得到等效转动惯量、飞轮的转动惯量、以及驱动电流对主轴扭矩、转速的依赖关系。其次要求能够计算出路试和台试两种情况下制动器耗能差或相对差,并以此评价计算机控制模型。最后在驱动电流对扭矩依赖关系的基础上,将连续制动时间离散化,设计出由上一时间段测量数据得到本阶段驱动电流的计算机控制方法
7、;同时指出此方法的不足并且改进、评价。二、问题分析针对问题一到问题三,我们可以简单的根据已有物理规律,如刚体定轴转动定理、动能定理、能量守恒以及基本概念计算和推理得到。问题四同样可以由刚体转动动能定理建模,利用先将连续制动时间离散化为小时间段,再将每个时间段内考虑为连续问题的方法计算台试时制动器耗能,路试简单用动能定理即可得到,二者求相对能量差,作为评价指标。问题五基本上建立在问题三的模型上面,只是在代入驱动扭矩时先采用波动差方法进行处理即可。问题六针对问题五中未考虑到计算机控制的时滞性等缺点,采用模糊控制方法即可优化方案。三、模型假设1. 制动过程中车轮承受载荷可以全部等效为车自重2. 电动
8、机的驱动电流正比于其产生的补偿扭矩3. 模型评价时不考虑观测误差、随机误差和连续问题离散化产生的误差四、 符号说明36-转速-角速度-能量相对误差-路试时制动器吸收能量-试验时制动器吸收能量-角加速度-制动扭矩-主轴扭矩(可观测)-机械惯量-等效转动惯量-瞬时补偿惯量-电动机补偿惯量五、模型建立与求解1.问题一:1.1问题的分析由假设1可知,制动时单个车轮承受的载荷可以等效的看作车座自重给该轮的分压力,则由此可以计算出等效质量。而车轮在平动过程中所具有的分能量,在忽略车轮自身转动能量后即为动能。再根据相关的物理公式即可将此部分动能转化为等效转动惯量。1.2模型的建立与求解物理公式查阅力学(第二
9、版)【6】得到公式并推导如下:由假设1可知,载荷可以全部转化为重力,即: 又汽车在载荷作用下平动的能量,忽略车轮转动产生的自身能量,剩下的仅为动能,有:另外,由刚体定轴转动的动能定理可知:两种模式下由能量守恒定律有: 以上四式联立可解得将代入上述公式得:,即等效转动惯量约为52kgm2 。2.问题二:2.1 问题分析每个飞轮在已知其半径、厚度和密度的前提下,其围绕中心轴旋转的转动惯量可以用物理积分方法得到。三个飞轮中可以任意选取其中的03个构成飞轮组,最后再与基础惯量共同得到可组合的8种机械惯量。经分析,上一问已经知道了等效转动惯量,又已知电动机补偿惯量区间,则易得机械惯量的取值范围,再从中挑
10、选。2.2 模型建立与求解各个飞轮转动惯量的计算公式为:同时我们有: 对于相同内、外径,相同密度的三个飞轮分别代入其厚度即可得惯量分别为:由上面分析可知,组合数为:种,分别表示如表1(用表示机械惯量各个飞轮):表1:各种飞轮组合的机械惯量组合指标NULL飞轮组惯量0306012090150180210基础惯量1010101010101010机械惯量104070130100160190220制动器试验台上,等效转动惯量是由机械转动惯量和电动机补偿的转动惯量之和,即:.由于电动机补偿转动惯量的范围是,而第一问已知等效转动惯量为,故机械惯量的取值区间为,从表中易知可以选择或飞轮单独使用,则电动机补偿
11、惯量分别为和。3.问题三:3.1 问题分析通过能量守恒原理推导可得到电动机驱动扭矩与制动扭矩的关系,而这两者之间可通过主轴扭矩联系,这样可将驱动扭矩与主轴扭矩的关系得到。根据假设2可知,驱动电流与驱动扭矩存在正比关系,这样即得模型。在已知驱动电流与驱动扭矩的比例常数、制动减速度亦为常数且可算的情况下,利用物理推导依次可得角加速度、电动机驱动扭矩,进而由正比系数可知两种飞轮组合情况下的驱动电流。3.2 模型建立与分析试验台制动过程中,制动器吸收的能量可由下表示: (1)式中为制动角速度, 为制动开始时间为制动结束时间而制动力矩为: (2)由第二问分析可知: (3)将(3)式带入(2)式中得: (
12、4)将(4)式带入(1)式中得: (5)由(5)可知制动器吸收的能量由两部分组成,其中制动器的负载中的部分由飞轮模拟,部分需要电机模拟。令: (6)因此,将(6)和(2)式中相除并整理可得: (7)同时对于主轴和飞轮组成的系统来说,共受到两个扭矩制动扭矩(负功)和电动机驱动扭矩(正功),则主轴扭矩(即可观测扭矩)为: (8)将(8)代入(7)式中可得到: (9)又已知驱动电流驱动扭矩成正比:因为题目中给出比例系数,故确定驱动电流为主轴扭矩的单变量函数: (10)有关分析:由(9)式可进行以下讨论:1) 如果,那么,电驱动系统按照电动机工作,在制动试验中提供扭矩以补充惯量;2) 如果,那么,电驱
13、动系统按照发电机工作,将制动试验中多余惯量对应的能量转化为电能。3.4 模型应用与求解计算角加速度: 计算电动机驱动扭矩:故由题目中假设驱动电流与驱动扭矩成正比可知:将等数据代入上述公式中,当时,计算可知 ;当时,计算可知。3.5 结果分析,表明此时为电动机产生驱动扭矩,以补充不足的惯量;而,表明此时为发电机,消耗多余惯量,并将其中储存的能量转化为电能。4. 问题四:4.1 模型的分析经分析可知评价控制方法优劣的重要数量指标为能量误差的大小,故现欲对题目中给出的控制方法进行评价,最基础的方法就是计算路试和试验两种制动情形下制动器消耗能量之差。而路试耗能可通过始末转速和等效转动惯量计算得到;试验
14、耗能可根据给出的数据每为周期分为机械惯量对应能量和电驱动系统能量两部分计算。然后计算二者的能量误差和能量相对误差,理论分析可知其值越小,说明模拟环境越接近实际路试环境,即该控制方法越好。4.2 模型的建立由第三问模型的建立过程,可知制动器在飞轮、电补偿的双分流载荷控制下消耗的能量: (11) 而实际路试的能量消耗为: (12) 这样,就可以知道模拟消耗的能量和路试消耗的能量误差为:4.3 模型的求解试验中的能量消耗的计算,主要是求出其中的,令:由上面的式(6)和式(9),可以得到: (13)其中,和都是关于时间变化的,可以看成都是关于时间的函数,分别为和。下面试图确定和:通过实际数据,做出扭矩
15、(如图1)和转速对时间的散点图(如图2)。图1:扭矩的观测点图2:转速的观测点图形分析:1) 由图可以知道,图一中扭矩与时间的关系无法用函数求出来。2) 通过图形可以知道,由于是以10ms为一个观测单位点,图形一和图形二都可以看成两条曲线,这样就可以把两个图形反过来看成曲线的离散化。3) 相邻两个点的距离很小,可以看成这段是线性关系。所以图形是有相邻线性关系。所以,从分析中可以看出和都不能整体由函数确定,故将整个制动时间分为n段,则有: ,而则为每一个分时间段对应的。由于相邻两个点间隔很短,可以简化成和与时间成线性关系的。取: (14) 将式(14)代回式(13)中有: 对上述结果通过编程即可
16、计算出各自的能量以及能量差。应用于具体题目中给出的数据,并且按照上述模型求解的思路,利用MATLAB7.0软件编程(见附录一)求解,得到数值如下表2:表2:能量指标评价法路试能量试验能量能量差相对能量差52.21754.9492.7325.23%4.4 结果分析利用模型计算得到表中数据可知,制动器在控制台试验时所消耗能量大于路试时所消耗能量,为了验证其科学性,我们做了如下分析:1) 制动器消耗能量多少取决于与之对应的惯量大小,而几个惯量中对吸收能量起决定性作用的是补偿惯量。所以我们只要比较控制台试验过程中瞬时补偿惯量的均值与补偿惯量的大小,即可从理论上判断上述模型的正确性。所以我们利用MATL
17、AB7.0编程(见附录二)求解得到瞬时补偿惯量的均值为,而实际系统电动机的补偿惯量应该为。显然有,所以制动器在控制台试验消耗的能量必然大于路试所消耗的能量,即上述模型准确性得到验证。2)路试和台试能量相对差为5.23%,可以说这种控制方法勉强可行,但是仍然存在不足,由待改进。5.问题五:5.1 问题分析虽然第三问已经建立了电流依赖于客观测量变化的模型,给出了电流推导公式,那么很容易根据上一阶段扭矩值得到本阶段驱动电流值,但鉴于该控制方法过分简化、涉及变量少、转速相互影响的实时效应被忽略了,所以考虑应用上述模型前做部分数据处理。电动机瞬时补偿惯量与被控制系统的需补偿惯量的差值越小,越贴近路试环境
18、,该控制台试验越成功。故先计算差值,然后确定需增加的驱动扭矩,代入第三问模型中得到相应驱动电流值。将该控制方法应用于可测扭矩和转速的每个时间段,即可得每个时间段驱动电流值,进而可利用计算机控制。5.2 模型建立与求解实际过程中的补偿惯量在模拟系统中已经定了,为:即:制动时等效转动惯量-机械惯量。由于这两个惯量都是由实验本身决定的,不受人为控制等产生影响。而控制电流的过程中,就是为了补偿机械惯量的不足而损失的能量,电流的主要作用就是补偿少的那部分惯量。在控制电流过程中,就得到需要的目标函数:在反映了实际需要补偿惯量与一个时间段产生的补偿惯量之差,换言之:现在电流产生的补偿惯量比实际中的多或者少的
19、量。它有三个结果:(1):当时,表示现阶段中产生的补偿惯量少于实际需要的,所以下一过程就需要通过改变电流来增加补偿惯量。(2):当时,表示现阶段中产生的补偿惯量刚好等于实验中需要的补偿惯量,这是模拟实验中最完美的结果。所以,下一阶段可以不用调节其中的电流,而使现有的状态持续到下一阶段。(3):当时,表示现阶段中产生的补偿惯量多于实验需要的,所以下一过程就需要通过改变电流来减少补偿惯量。控制过程中,从实际出发,的变动越小越好。因为如果的变动范围越大,下一个时间段需要改变的范围就越大。而电流具有滞后性,不可能马上就达到实验需要的值。这样,实验消耗的能量就很大,而且产生的误差也很大。根据第三问导出的
20、模型可以知道: (其中为电动机驱动扭矩;为观测到得主轴扭矩;为瞬时补偿惯量;制动时等效转动惯量;为机械惯量)从结果可以看出,由于和都是常数。而和都是实验中观测的值,故每一段时间的这个变化率是可以计算出来的。所以,这段时间中产生的就可以确定。根据公式可以最后的表达形式:从结果可以看出:和是成正关系的。那么就有:电流增加,补偿惯量增加;电流减小,得到的补偿惯量就减小。通过上面的分析后,得到控制过程:将连续的时间段离散化,把时间段分为n个区间,其中每一个区间的时间长度都趋向于0。对于每一个时间段,可以看成转速是均匀变化的,则有(当为负数时,只是方向问题)所以通过判断上一阶段的正负情况,来判断下一段时
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 法律资料 法律资料A20004051 中南大学 辛婧华 高乐 甘磊 法律 资料 A20004051 中南 大学
链接地址:https://www.31doc.com/p-1984678.html