水平运动物体的建模分析及控制涉及 点击:6281 | 回复:237



    
  • 精华:2帖
  • 求助:0帖
  • 帖子:10帖 | 2393回
  • 年度积分:1
  • 历史总积分:5804
  • 注册:2006年1月06日
发表于:2011-03-12 09:58:52
楼主

 

1 控制对象的建模:
假定运动物体质量为M, 对物体进行受力分析 
水平方向的力:
1) 运动过程中受到的与速度相关的阻力与运行速度成正比, 比例系数为Kv,单位为N*s/m  (牛秒/米)
2) 滑动摩擦系数为u, 即受到的滑动摩擦力为uG
3) 牵引力F
竖直方向的力:
1) 地面的支持力N
2) 物体的重力G


物体的受力分析如下图:


竖直方向的力N与G相等, 不会发生运动。
水平方向的运动方程: F = Ma + Kv V + uG
       式中a为加速度m/s2, V为物体的运行速度 m/s,u为滑动摩擦系数,G= Mg

即:  F= M dv/dt + Kv * V + uMg
建立以下的数学模型:

控制对象的传递函数为:

 

M/Kv的单位:Kg / (N*s/m)= Kg * m / (N * s) = Kg * m / (Kg * m * s(-2 )*s) = s  单位为秒


因此控制对象等效为一阶环节, 表示如下:



式中 Tm, 单位为秒, 即控制对象的时间常数, 通常以毫秒ms表示
单位为Kv , 即牛秒/米

这一步已经推导出了控制对象的数学模型, 从数学模型上看, 控制对象已经被抽象成一阶环节, 但模型的参数是根据控制对象得来的, 每一个参数都有其特定的物理含义!

2。执行环节

要想对控制对象进行控制, 必须要有执行机构。

对任何执行机构而言,都有功率及输出作用力的限制,即执行机构不可能是功率无限大、输出作用力也不可能是无限大的。

从功率的计算公式P=FV来看, 特定的执行单元, 功率一定的情况下, 在运行达到平衡及加速度a = 0时, 其速度达到最大值!同样, 执行单元也不可能在速度v=0时, 输出无限大的作用力F。

这个特性由实际的执行机构所决定。

 

在这个模型中, 我们假定执行机构是线性放大环节, 输入为控制电压, 输出为作用力F, 因此执行机构的数学模型可以简单的描述为比例放大环节, 其单位为牛/伏, 即N/V。

上述简化中忽略了以下因素:

1)执行环节的时滞

2)执行环节的功率限制及作用力限制。

3)如果对执行机构采用内环输出力的闭环控制, 忽略内环的带宽和延迟带来的影响。

 

考虑到系统的快速响应及更好的控制性能, 可以对执行单元实行控制力的内环控制。在我们的模型中, 简单的认为其为比例放大环节。

执行机构的最大输出作用力会在控制器设计时进行限定!

执行机构的功率大小在后面的设计中也会予以考虑分析!

执行机构的特性可以用下面的图表示:

 在功率一定的情况下, 小于V0的速度时执行机构可以输出最大作用力, V0 = P/Fmax。

当速度高于V0时, 输出作用力与速度成反比。

在我们的模型设计中, 只限制最大输出作用力Fmax, 并将执行机构理想化为线性比例放大环节

 3. 控制系统设计

假设物体质量M=10Kg,滑动摩擦系数u=0.2, Kv=0.1,执行机构线性放大倍数为100

控制目标:物体运行速度10m/s, 阶跃信号给定,超调量小于5%

则开环控制时需要的作用力F=0.2*10*9.8+10 *0.1=20.6N

对控制过程仿真如下:

 图中可以看出,物体速度达到9m/s的时间超过200秒,控制对象存在很大时间延迟

采用下图的闭环控制:

整定PI控制器参数, 选取阻尼系数为0.707
kp =14.1411, Ti= 14.14ms 
不考虑执行机构的最大输出作用力限制时,仿真如下:
超调量4.3%, 上升时间小于40ms但是最大作用力F达到了将近5000牛顿。
考虑执行机构的最大输出作用力限制时,限制F最大为200N时,仿真如下:
达到指令速度时间为550ms,超调量小于4.3%
(转载请注明作者、出处)



刘志斌

  • 精华:5帖
  • 求助:0帖
  • 帖子:346帖 | 35827回
  • 年度积分:0
  • 历史总积分:75384
  • 注册:2006年7月30日
发表于:2011-07-13 08:39:30
81楼

下面看楼主的“3. 控制系统设计 ”:

1、楼主说:

“3. 控制系统设计

假设物体质量M=10Kg,滑动摩擦系数u=0.2, Kv=0.1,执行机构线性放大倍数为100

控制目标:物体运行速度10m/s, 阶跃信号给定,超调量小于5%

则开环控制时需要的作用力F=0.2*10*9.8+10 *0.1=20.6N ”

刘志斌

  • 精华:5帖
  • 求助:0帖
  • 帖子:346帖 | 35827回
  • 年度积分:0
  • 历史总积分:75384
  • 注册:2006年7月30日
发表于:2011-07-13 08:46:24
82楼

2、楼主说:“开环控制时需要的作用力F=0.2*10*9.8+10 *0.1=20.6N ”,也就是说这个作用力F20.6N 只是维持系统匀速10m/s运动时的力;

3、运动物体的速度由0到10m/s的加速度的力没有;

刘志斌

  • 精华:5帖
  • 求助:0帖
  • 帖子:346帖 | 35827回
  • 年度积分:0
  • 历史总积分:75384
  • 注册:2006年7月30日
发表于:2011-07-13 09:16:08
83楼

4、如果运动物体从0到10m/s期间的作用力始终是20.6N ,那么启动时的加速度的力是:

      20.6N - 0.2*10*9.8 =  20.6N - 19.6N = 1N

当时明月

  • 精华:0帖
  • 求助:0帖
  • 帖子:0帖 | 138回
  • 年度积分:0
  • 历史总积分:533
  • 注册:2006年6月21日
发表于:2011-07-13 09:43:29
84楼
我嘞个去~~~~

刘志斌

  • 精华:5帖
  • 求助:0帖
  • 帖子:346帖 | 35827回
  • 年度积分:0
  • 历史总积分:75384
  • 注册:2006年7月30日
发表于:2011-07-13 11:23:07
85楼

5、如果运动物体从0到10m/s期间的作用力始终是20.6N ,那么启动时的加速度的力是: 

       20.6N - 0.2*10*9.8 = 20.6N - 19.6N = 1N

       最大加速度是:dv/dt=F/M=1N/10Kg=0.1m/ss

刘志斌

  • 精华:5帖
  • 求助:0帖
  • 帖子:346帖 | 35827回
  • 年度积分:0
  • 历史总积分:75384
  • 注册:2006年7月30日
发表于:2011-07-13 11:31:29
86楼

6、如果运动物体从0到10m/s期间的作用力始终是20.6N ,那么启动到速度为10m/s时的加速度是0:

7、平均加速度是    ( 最大加速度+0)/2=0.05m/ss

8、如果运动物体从0到10m/s期间的作用力始终是20.6N ,那么启动到速度为10m/s时,共用时间是:

          时间=V/加速度=10/0.05=200(s)

luedong

  • 精华:0帖
  • 求助:0帖
  • 帖子:2帖 | 607回
  • 年度积分:0
  • 历史总积分:1203
  • 注册:2011年5月26日
发表于:2011-07-13 11:33:19
87楼
51楼问题呢?s/M是什么意思呢? 1/Ms可以理解,1/s是个积分环节,表达的意思是对加速度的积分,就是速度v。 推导出s/M的公式,有什么物理含意?秒/千克?恐怕要改写传递函数、微分方程、拉氏变换的历史发展过程了!

刘志斌

  • 精华:5帖
  • 求助:0帖
  • 帖子:346帖 | 35827回
  • 年度积分:0
  • 历史总积分:75384
  • 注册:2006年7月30日
发表于:2011-07-13 11:41:55
88楼
9、如果解微分方程   dv/dt=(20.6 -19.6 - 0.1V) /10    也可以得出时间约为200s 

刘志斌

  • 精华:5帖
  • 求助:0帖
  • 帖子:346帖 | 35827回
  • 年度积分:0
  • 历史总积分:75384
  • 注册:2006年7月30日
发表于:2011-07-13 14:07:08
89楼

对控制过程仿真如下:

上图中传输函数10/(1+100s)是错误的,应该为10/(1+100/s);

 

mengx9806

  • 精华:0帖
  • 求助:0帖
  • 帖子:1帖 | 57回
  • 年度积分:0
  • 历史总积分:258
  • 注册:2011年5月21日
发表于:2011-07-13 14:41:47
90楼

分析中,理解中,自愧自控原理没学好中~~~~~我想哭。

天高云淡

  • 精华:0帖
  • 求助:0帖
  • 帖子:0帖 | 46回
  • 年度积分:0
  • 历史总积分:435
  • 注册:2010年12月21日
发表于:2011-07-13 21:04:42
91楼

刘老师丢人了吧!

这么简单的问题你都搞错了!

还在这自娱自乐12345的!

 

刘志斌

  • 精华:5帖
  • 求助:0帖
  • 帖子:346帖 | 35827回
  • 年度积分:0
  • 历史总积分:75384
  • 注册:2006年7月30日
发表于:2011-07-13 23:18:13
92楼
引用天高云淡 的回复内容:

刘老师丢人了吧!

这么简单的问题你都搞错了!



1、欢迎把错误贴出来!

2、说出错误是什么?正确的是什么?

刘志斌

  • 精华:5帖
  • 求助:0帖
  • 帖子:346帖 | 35827回
  • 年度积分:0
  • 历史总积分:75384
  • 注册:2006年7月30日
发表于:2011-07-14 00:55:03
93楼

引用luedong 的回复内容:51楼问题呢?s/M是什么意思呢? 1/Ms可以理解,1/s是个积分环节,表达的意思是对加速度的积分,就是速度v。推导出s/M的公式,有什么物理含意?秒/千克?恐怕要改写传递函数、微分方程、拉氏变换的历史发展过程了!

 


 


1、开环,即Kv = 0 ;

2、F=定值;

3、动力学方程:    dv/dt =F/M

                                    dv=(F/M)dt

                                   ∫dv=(F/M)∫dt

                                      V=(F/M)s

                                   V/F=s/M

4、所以开环传输函数为 s/M  ;

5、所以“征”的开环传输函数 1/Ms 是错误的;

刘志斌

  • 精华:5帖
  • 求助:0帖
  • 帖子:346帖 | 35827回
  • 年度积分:0
  • 历史总积分:75384
  • 注册:2006年7月30日
发表于:2011-07-14 01:06:42
94楼
引用征 的回复内容:

你懂 s 是什么吗?


1、s是时间,是运动体加速的时间;

2、你认为s是什么?

刘志斌

  • 精华:5帖
  • 求助:0帖
  • 帖子:346帖 | 35827回
  • 年度积分:0
  • 历史总积分:75384
  • 注册:2006年7月30日
发表于:2011-07-14 01:09:45
95楼
引用征 的回复内容:

输出/输入=V/F=(m/s)/(M×m/ss)=s/M

自控原理理没有这样的东西!


1、开环传输系数的简单计算,你应该认识;

2、如果不认识,请看93楼的推导过程;

刘志斌

  • 精华:5帖
  • 求助:0帖
  • 帖子:346帖 | 35827回
  • 年度积分:0
  • 历史总积分:75384
  • 注册:2006年7月30日
发表于:2011-07-14 01:50:14
96楼

引用征 的回复内容:

……借用35楼的描述

根据上图的微分性质:

L[dv/dt] = sV(s) - v(0)

初值为v(0) = 0时, dv/dt = F/M 拉氏变换后可以表示为:

sV(s) = F(s)/M

即 V(s) / F(s) = 1/Ms


1、上述过程错误:

2、直接分离变量  

        dv/dt = F/M

             dv =  F/M  dt      (F、M 均为常量)

            ∫dv = F/M ∫dt     (积分上限分别为 V 、 s;积分下限分别为 0 、0)

               V = (F/M) s

            V/F = s /M

刘志斌

  • 精华:5帖
  • 求助:0帖
  • 帖子:346帖 | 35827回
  • 年度积分:0
  • 历史总积分:75384
  • 注册:2006年7月30日
发表于:2011-07-14 02:00:38
97楼

3、根据你的举例:

     开环传输函数 V/F = s /M 验证如下:

     左边=  V/F=10/(20.6-19.6)= 10

     右边=  s /M = 100/10 =10

 

刘志斌

  • 精华:5帖
  • 求助:0帖
  • 帖子:346帖 | 35827回
  • 年度积分:0
  • 历史总积分:75384
  • 注册:2006年7月30日
发表于:2011-07-14 02:12:59
98楼

4、如果用你“征”的结果是:

开环传输函数 V/F = 1/M s 验证如下:

左边= V/F=10/(20.6-19.6)= 10

右边= 1/M s =1/(10×100)=1/1000=0.001

刘志斌

  • 精华:5帖
  • 求助:0帖
  • 帖子:346帖 | 35827回
  • 年度积分:0
  • 历史总积分:75384
  • 注册:2006年7月30日
发表于:2011-07-14 02:16:45
99楼
5、“征”你错了,把你骂人的话自己删掉吧!

luedong

  • 精华:0帖
  • 求助:0帖
  • 帖子:2帖 | 607回
  • 年度积分:0
  • 历史总积分:1203
  • 注册:2011年5月26日
发表于:2011-07-14 06:35:30
100楼
对刘老师的解析彻底无语!传递函数中的s是加速时间?一代数学家、天文学家拉普拉斯的拉普拉斯变换如下: http://baike.baidu.com/view/132034.htm

热门招聘
相关主题

官方公众号

智造工程师