用户中心

资讯 > 封面特写

模型与仿真基础

作者:Vance VanDoren,CONTROL ENGINEERING 编辑顾问2007.03.30阅读 4238

  依靠数学模型仿真,工程师可以对连续性的过程控制系统进行模拟测试,并根据测试推算系统以后的运行情况。
  可塑模型和陶土模型展现的是对象的外形特征,和这种模型不同,过程模型则反映受控系统的操作过程以及对外界干扰的响应情况。通常,过程模型都可以用数学方程来确定系统输入量和输出量之间的关系。

图1通过这个简单的模型可以推算出当水以每分钟F立方米的速率填充一个底半径为R米的圆柱体水箱t分钟后水面高度L(单位米)是多少
图1:通过这个简单的模型,可以推算出,当水以每分钟F立方米
的速率填充一个底半径为R米的圆柱体水箱,t 分钟后水面高度L
(单位:米) 是多少。

  模型可以有效地帮助我们分析PID回路,也可以为过程控制系统的设计提供参照。举个例子来说,如图1所示CONTROL ENGINEERING China版权所有,是给一个圆柱形水箱充水的过程。假设水流入的速率是F,在t分钟后,水箱里的水面高度为:
  [1] 
  其中,R是水箱底部圆的半径。这样一个模型可以用来推算,若要装满整个水箱,控制器需要运行多长时间。
  模型的要素
  过程模型越复杂,包括的变量就越多,他们之间的数学关系也更加复杂,但不管系统多复杂,凡是连续过程的模型都包括以下四个基本要素:
  ■输入变量;
  ■输出变量;
  ■常数项
  ■算法
  输出是在开发模型时,通过输入量的值分析得出的

预期值。在填充水箱的例子中,输出量L就可以由输入量t和F的值推算出来。
  常数R的值取决于水箱的大小。一般情况下,常数是反映各基础科学的基本原理的常量,如物理学、化学、经济学、几何学等CONTROL ENGINEERING China版权所有,它决定着过程的基本性能,而且这些常数是恒定不变的,不会因输入输出值随时间的变化而变化。
  算法包括多种数学运算方法,主要是用来根据输入和常量值推算输出值。简单的像乘法和乘方,如方程[1]所示;复杂的如拉普拉斯变换和统计分配定律。
  控制模型
  过程模型可以为设计、执行、反馈控制测试提供有效的帮助。大多数解析法需要一个含增益和时间常数在内的模型,通过监测系统的反馈和响应速度,来判断该控制系统的控制性能。当对模型参数有一定了解后,工程师可以给过慢的进程增加一个响应较快的控制器,或是给过快的进程增加一个较慢的控制器,以达到预期的效果。
  同样,模拟预测控制器可以用数学模型来测定控制系统的性能控制工程网版权所有,以判断受控变量是否精确的达到了控制要求。在运行过程中通过比较得出的反馈,可以进行自我调整。
  即使某控制方案是基于某特殊应用而开发的,在进行实际应用之前控制工程网版权所有,也可以利用过程模型在虚拟环境下对该控制方案进行测试。模型的控制方程可以直接编程写入,也可以几种常用仿真语言中的任意一种或者自定义码写入专用的调试机中。在计算机上进行仿真,可以更快的发现该控制方案中的不足,也不用担心会对真实的工程造成任何损失。
  模型开发实例
  不管用途是什么控制工程网版权所有,建立模型的诀窍就是把系统的运行状态用一系列的控制方程表示出来,如方程[1]所示。以"倒置的单摆模型"为例,如图所示,在地面上固定一个弹簧,弹簧顶端固定安放一个重物。水平方向的力使该物体在一个弧形范围内来回摆动。
  这个模型的运行状态经过变形可以被应用于多个仿真研究。可以把它看成一个儿童玩具,也可以把它看成某高层建筑物在风中摇晃的缩影。加上一个接头,也可以近似看成是人在迈腿大步行走时腿部的动作。
  无论是哪种情况,基本原理都是相同的。弹簧对物体的弹力和物体自身的重力是相对的。如果将物体和弹簧合起来,看成是一个质心离地面H米,质量为m千克的联合装置,整个装置的运动状态就可以用含有弧度q的方程式表示,如图所示。
  实际应用检测
  正如这个装置的简化过程一样,这个模型是逼近现实的模拟。假定不存在任何外力影响物体的运动,如摩擦力等,而弹簧对物体的弹力也与角坐标θ完全成比例。另外还要假定,该运动过程是从一个绝对垂直位置开始的(即θ的初始值为零)。
  由于当θ无穷小时,q和sin(θ)的值近似相等。用q替换sin(θ)方程式[2]可以简化成方程式[3],根据式[3]画出关于θ的趋势图,得到一条斜率为-(k- mg)/mh的直线,因此式[3]为关于θ的线性方程。
  倒置单摆模型
  [2]   简化后的倒置单摆模型  [3]  对简化模型方程求解  [4]  由
  简化后的倒置单摆模型
  [3]
  对简化模型方程求解
  [4]
  由上式解得
  [5]   [6]  角速度变量的θ”是角速度θ的二次导数t为该物体从初始位置到现在的运动时间θ'0是物体在外力作用下的初始角速度A是…
  [6]
  角速度变量的θ”是角速度θ的二次导数,t为该物体从初始位置到现在的运动时间,θ'0是物体在外力作用下的初始角速度,A是物体摆动的振幅,ω为振动的频率。A和ω是两个值取决于k,m,g,h和θ'的常量,如等式[5]和等式[6]所示。

  方程式[3]真正的价值在于,通过它可以得到一个θ关于时间的趋势图,根据趋势图,我们可以预测该过程以后的运行状态。这个碰巧是一个只用方程式[3]这种线性方程就可以解决的问题。实际上,对方程式[3]进行求解,我们也可以得到关于θ(t)的表达式,如方程式[4]所显示。

图2该弹簧装置的总质量为m千克质心离地高度为h米且正好在弹簧支BR
图2:该弹簧装置的总质量为m千克,质心离地高度为h 米,且正好在弹簧支





































版权声明:版权归控制工程网所有,转载请注明出处!

频道推荐

关于我们

控制工程网 & CONTROL ENGINEERING China 全球工业控制、自动化和仪器仪表领域的先锋媒体

CE全球

联系我们

商务及广告合作
任小姐(北京)                 夏小姐(上海)
电话:010-82053688      电话:18616877918
rendongxue@cechina.cn      xiashuxian@cechina.cn
新闻投稿:王小姐

关注我们的微信

关于我们 | 网站地图 | 联系我们
© 2003-2020    经营许可编号:京ICP证120335号
公安机关备案号:110102002318  服务热线:010-82053688