水平运动物体的建模分析及控制涉及 点击:6276 | 回复: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-12 18:29:16
41楼
引用luedong 的回复内容:

……

刘老师 , 你知道s是什么吗?

F=Ma = M dv/dt

dv/dt = F/M

用拉氏函数表示就是:

sV(s) = F(s)/M

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

传递函数就是:

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


F=Ma = M dv/dt

dv/dt = F/M

用拉氏函数表示就是:

sV(s) = F(s)/M

dv/dt =V(s)/s

dv/dt ≠sV(s)   你比你的朋友强一点,还能辩解,你的朋友干脆一句辩解的话也说不出!

所以应该是:V(s)/s =F(s)/M    即  V(s) / F(s) = s /M

  • 精华:2帖
  • 求助:0帖
  • 帖子:10帖 | 2393回
  • 年度积分:1
  • 历史总积分:5804
  • 注册:2006年1月06日
发表于:2011-07-12 18:35:21
42楼
引用luedong 的回复内容:

刘老师还是不要点评了!

越点评错的越多!

楼主这个帖子还是写的不错的, 基于自控原理的根本知识!

算是自控原理入门级内容!


 


 

解释什么呢?

老刘说我这帖子抄的, 笑话!

就这个帖子中的小儿科问题, 我还用去抄?

  • 精华:2帖
  • 求助:0帖
  • 帖子:10帖 | 2393回
  • 年度积分:1
  • 历史总积分:5804
  • 注册:2006年1月06日
发表于:2011-07-12 18:36:35
43楼

老刘, 懒得搭理他!

 

刘志斌

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

其实从传输函数的定义出发,直接可以写成:

V(s) / F(s) = (m/s)/(kg×m/ss)=(m/s)/(M×m/ss)=s/M

刘志斌

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

引用luedong 的回复内容: 

这个式子有什么错误?

分子分母同时Ms/Kv, 不就得到后面的式子了?

刘老师看出哪里错了?



1、你乘一下,看看错了没有!

2、分母应该是 1+Ms/Kv。

刘志斌

  • 精华:5帖
  • 求助:0帖
  • 帖子:346帖 | 35827回
  • 年度积分:0
  • 历史总积分:75384
  • 注册:2006年7月30日
发表于:2011-07-12 18:49:37
46楼
引用luedong 的回复内容:

这帖子又白瞎了!

呵呵!


1、是啊!只要说他错误就“帖子又白瞎了”;

2、看看这一对一唱一和的骂人高手的不讲理的狼狈相!!!!!!!!

刘志斌

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

引用luedong 的回复内容: braverider在27楼已经指出你的错误了! 一阶环节都看不出来! 自控原理怎么学的?

--------------------------------------------------------------------------------

1、自控原理,关于系统传输函数,因该是各级的传输函数的乘积;

2、做为控制对象的一级传输函数,应该是图中的




 

 

 


luedong

  • 精华:0帖
  • 求助:0帖
  • 帖子:2帖 | 607回
  • 年度积分:0
  • 历史总积分:1203
  • 注册:2011年5月26日
发表于:2011-07-12 19:07:22
48楼
引用刘志斌 的回复内容:

引用luedong 的回复内容: 

这个式子有什么错误?

分子分母同时Ms/Kv, 不就得到后面的式子了?

刘老师看出哪里错了?

 


 


1、你乘一下,看看错了没有!

2、分母应该是 1+Ms/Kv。




刘老师。 你没看到人家下面还写了一个公式?

G(s) = (1/Kv) / (1+ Tm s)

对比两个公式, 你看不出来什么意思吗?


luedong

  • 精华:0帖
  • 求助:0帖
  • 帖子:2帖 | 607回
  • 年度积分:0
  • 历史总积分:1203
  • 注册:2011年5月26日
发表于:2011-07-12 19:09:20
49楼
引用刘志斌 的回复内容:引用luedong 的回复内容:

这帖子又白瞎了!

呵呵!

 


 

1、是啊!只要说他错误就“帖子又白瞎了”;

2、看看这一对一唱一和的骂人高手的不讲理的狼狈相!!!!!!!!




人家骂你是垃圾, 我可没骂你刘老师, 呵呵

咱只分析问题, 就问题说事!

你说说上面的那些式子错了, 咱学习学习!


luedong

  • 精华:0帖
  • 求助:0帖
  • 帖子:2帖 | 607回
  • 年度积分:0
  • 历史总积分:1203
  • 注册:2011年5月26日
发表于:2011-07-12 19:17:46
50楼
引用刘志斌 的回复内容:

其实从传输函数的定义出发,直接可以写成:

V(s) / F(s) = (m/s)/(kg×m/ss)=(m/s)/(M×m/ss)=s/M




刘老师,真心请教:

上面的公式代表什么意思?

luedong

  • 精华:0帖
  • 求助:0帖
  • 帖子:2帖 | 607回
  • 年度积分:0
  • 历史总积分:1203
  • 注册:2011年5月26日
发表于:2011-07-12 19:19:34
51楼

V(s) / F(s) = (m/s)/(kg×m/ss)=(m/s)/(M×m/ss)=s/M

请刘老师详细解释一下上面公式怎么得出的?

非常感谢!

luedong

  • 精华:0帖
  • 求助:0帖
  • 帖子:2帖 | 607回
  • 年度积分:0
  • 历史总积分:1203
  • 注册:2011年5月26日
发表于:2011-07-12 19:22:27
52楼

我姑且猜测:

V(s)    m/s       速度的表示单位:米/秒

F(s)  kg×m/ss    力的表示单位 Mg: 千克 * 米/(秒 * 秒)

然后一除:V(s) / F(s) = (m/s)/(kg×m/ss)=(m/s)/(M×m/ss)=s/M

是这样吗 ?

刘老师?

 

  • 精华:2帖
  • 求助:0帖
  • 帖子:10帖 | 2393回
  • 年度积分:1
  • 历史总积分:5804
  • 注册:2006年1月06日
发表于:2011-07-12 19:43:45
53楼
引用刘志斌 的回复内容:

其实从传输函数的定义出发,直接可以写成:

V(s) / F(s) = (m/s)/(kg×m/ss)=(m/s)/(M×m/ss)=s/M


 


 

其实从传输函数的定义出发,直接可以写成:

V(s) / F(s) = (m/s)/(kg×m/ss)=(m/s)/(M×m/ss)=s/M

这个公式写的就跟白痴一样!

  • 精华:2帖
  • 求助:0帖
  • 帖子:10帖 | 2393回
  • 年度积分:1
  • 历史总积分:5804
  • 注册:2006年1月06日
发表于:2011-07-12 19:50:01
54楼

拉氏变换的微分性质:

  • 精华:2帖
  • 求助:0帖
  • 帖子:10帖 | 2393回
  • 年度积分:1
  • 历史总积分:5804
  • 注册:2006年1月06日
发表于:2011-07-12 19:55:58
55楼
引用luedong 的回复内容:引用刘志斌 的回复内容:

 

上图中1/Ms是错误的,应该是s/M。


 


 


s/M ?

刘老师 , 你知道s是什么吗?

F=Ma = M dv/dt

dv/dt = F/M

用拉氏函数表示就是:

sV(s) = F(s)/M

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

传递函数就是:

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






借用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


  • 精华:2帖
  • 求助:0帖
  • 帖子:10帖 | 2393回
  • 年度积分:1
  • 历史总积分:5804
  • 注册:2006年1月06日
发表于:2011-07-12 19:58:05
56楼
引用刘志斌 的回复内容:

其实从传输函数的定义出发,直接可以写成:

V(s) / F(s) = (m/s)/(kg×m/ss)=(m/s)/(M×m/ss)=s/M





写这样的白痴公式, 咋说呢?

水平不是一般的高!

 

  • 精华:2帖
  • 求助:0帖
  • 帖子:10帖 | 2393回
  • 年度积分:1
  • 历史总积分:5804
  • 注册:2006年1月06日
发表于:2011-07-12 19:59:31
57楼
引用luedong 的回复内容:

我姑且猜测:

V(s)    m/s       速度的表示单位:米/秒

F(s)  kg×m/ss    力的表示单位 Mg: 千克 * 米/(秒 * 秒)

然后一除:V(s) / F(s) = (m/s)/(kg×m/ss)=(m/s)/(M×m/ss)=s/M

是这样吗 ?

刘老师?

 


 


 


歇歇吧!


  • 精华:2帖
  • 求助:0帖
  • 帖子:10帖 | 2393回
  • 年度积分:1
  • 历史总积分:5804
  • 注册:2006年1月06日
发表于:2011-07-12 20:09:35
58楼
欢迎讨论, 但是刘志斌这样的人, 不欢迎!

luedong

  • 精华:0帖
  • 求助:0帖
  • 帖子:2帖 | 607回
  • 年度积分:0
  • 历史总积分:1203
  • 注册:2011年5月26日
发表于:2011-07-12 20:13:52
59楼

不好意思,我还等刘老师给解释一下51楼的问题呢!

这么高深的传递函数公式, 咱得学习一下!

刘志斌

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

我来说说:

1、这个以控制对象为一级环节的传输函数是:

1)开环放大倍数可以视为:输出/输入=V/F=(m/s)/(M×m/ss)=s/M

2)闭环反馈系数可以视为:Kv=F/ V

3)则一级传输函数可视为:G(s)= (s/M) / 1+Kv(s/M) = (1/Kv) / (1+M/Kvs) 


热门招聘
相关主题

官方公众号

智造工程师