倒立摆模型分析_倒立摆系统建模方法

倒立摆模型分析_倒立摆系统建模方法模型对小车的水平受力分析Mx¨=F−bx˙−NM\ddot{x}=F-b\dot{x}-NMx¨=F−bx˙−Nx¨\ddot{x}x¨代表对运动距离的二阶微分,即小车在外力作用下的加速度。FFF是外部施加给系统的外力。x˙\dot{x}x˙代表小车当前的运动速度,小车所受到的摩擦力为摩擦系数与小车运动速度之积,即f=bx˙f=b\dot{x}f=bx˙。NNN为倒单摆作用给小车水平方向的力。对倒单摆的受力分析摆杆做平面运动,可以分解为质心的平动和绕质心的转动。N=mx¨−(mlψ

大家好,又见面了,我是你们的朋友全栈君。如果您正在找激活码,请点击查看最新教程,关注关注公众号 “全栈程序员社区” 获取激活教程,可能之前旧版本教程已经失效.最新Idea2022.1教程亲测有效,一键激活。

Jetbrains全系列IDE使用 1年只要46元 售后保障 童叟无欺

参考及致谢

一阶倒立摆的PID控制和LQR控制
由拉普拉斯变换到传递函数
函数f(t)二阶导数的拉普拉斯变换是什么?
[1]翟龙余.一级倒立摆仿真模型的建立[J].大众科技,2011(8):268-270.

模型建立

在这里插入图片描述
在这里插入图片描述

对小车的水平受力分析

M x ¨ = F − b x ˙ − N M\ddot{x}=F-b\dot{x}-N Mx¨=Fbx˙N

  1. x ¨ \ddot{x} x¨代表对运动距离的二阶微分,即小车在外力作用下的加速度。
  2. F F F是外部施加给系统的外力。
  3. x ˙ \dot{x} x˙代表小车当前的运动速度,小车所受到的摩擦力为摩擦系数与小车运动速度之积,即 f = b x ˙ f=b\dot{x} f=bx˙
  4. N N N为倒单摆作用给小车水平方向的力。

对倒单摆的水平受力分析

摆杆做平面运动,其质心在外力作用下,在一段时间内的水平位移为 s = x − l × s i n ψ s= x-l×sin\psi s=xl×sinψ (因为倒立摆的倒向与外力 F F F的方向相反,所以中间用负号),其加速度可以表示成
s ¨ = d 2 s d t 2 = d 2 ( x − l × s i n ψ ) d t 2 = x ¨ − l ( d 2 s i n ψ d t 2 ) \ddot{s}=\frac {d^2s}{dt^2}=\frac{d^2( x-l×sin\psi)}{dt^2}=\ddot{x}-l(\frac{d^2sin\psi}{dt^2}) s¨=dt2d2s=dt2d2(xl×sinψ)=x¨l(dt2d2sinψ)
= x ¨ − l c o s ψ d ψ d t d t = x ¨ − l [ − s i n ψ ( d ψ d t ) 2 + c o s ψ d 2 ψ d t 2 ] =\ddot{x}-l\frac{cos\psi\frac{d\psi}{dt}}{dt}=\ddot{x}-l[-sin\psi(\frac{d\psi}{dt})^2+cos\psi \frac{d^2\psi}{dt^2}] =x¨ldtcosψdtdψ=x¨l[sinψ(dtdψ)2+cosψdt2d2ψ]
更换符号后即可得到:
s ¨ = x ¨ + l s i n ψ ( ψ ˙ ) 2 − l c o s ψ ( ψ ¨ ) \ddot{s}=\ddot{x}+lsin\psi(\dot{\psi})^2-lcos\psi(\ddot{\psi}) s¨=x¨+lsinψ(ψ˙)2lcosψ(ψ¨)
根据牛顿第二定律,此时摆质心的受力与加速度的关系为:
N = m x ¨ − ( m l ψ ¨ ) c o s ψ + m l ψ 2 ˙ s i n ψ N=m\ddot{x}-(ml\ddot{\psi})cos\psi+ml\dot{\psi ^2}sin\psi N=mx¨(mlψ¨)cosψ+mlψ2˙sinψ
联立关于倒立摆与小车的受力分析,替换掉相互作用力 N N N,得到:
( M + m ) x ¨ + b x ˙ − m l ψ ¨ c o s ψ + m l ψ ˙ 2 s i n ψ = F (M+m)\ddot{x}+b\dot{x}-ml\ddot{\psi}cos\psi+ml\dot{\psi}^2sin\psi=F (M+m)x¨+bx˙mlψ¨cosψ+mlψ˙2sinψ=F

对倒单摆的垂直受力分析

倒单摆的质心在一段时间内垂直方向上移动的距离可以表示成: h = l c o s ψ h=lcos\psi h=lcosψ
式中 ψ \psi ψ为单摆绕轴心转动的角度。
摆的质心在垂直方向的加速度可以表示为(注意,此时加速度方向与重力方向一致):
h ¨ = d 2 ( l c o s ψ ) d t 2 = l d ( − s i n ψ d ψ d t ) d t = − l c o s ψ ( d ψ d t ) 2 − l s i n ψ ( d 2 ψ d t 2 ) \ddot{h}=\frac{d^2(lcos\psi)}{dt^2}=l\frac{d(-sin\psi\frac{d\psi}{dt})}{dt}=-lcos\psi(\frac{d\psi}{dt})^2-lsin\psi(\frac{d^2\psi}{dt^2}) h¨=dt2d2(lcosψ)=ldtd(sinψdtdψ)=lcosψ(dtdψ)2lsinψ(dt2d2ψ)
垂直方向有重力 m g mg mg和小车对摆的支持力 P P P,另外单摆会有一个与重力方向一致的加速度。
垂 直 向 上 的 分 量 = 垂 直 向 下 的 分 量 垂直向上的分量=垂直向下的分量 =
P = m g + m h ¨ P=mg+m\ddot{h} P=mg+mh¨
P = m g − m l c o s ψ ( d ψ d t ) 2 − m l s i n ψ ( d 2 ψ d t 2 ) P=mg-mlcos\psi(\frac{d\psi}{dt})^2-mlsin\psi(\frac{d^2\psi}{dt^2}) P=mgmlcosψ(dtdψ)2mlsinψ(dt2d2ψ)
替换符号之后可以得到:
P = m g − m l ( c o s ψ ) ψ ˙ 2 − m l ( s i n ψ ) ψ ¨ P=mg-ml(cos\psi)\dot{\psi}^2-ml(sin\psi)\ddot{\psi} P=mgml(cosψ)ψ˙2ml(sinψ)ψ¨
假设摆受力不平衡,会有以铰链为圆心的角加速度,将 P P P N N N分别在转动方向上投影,根据倒单摆平衡时的力矩方程方程得到:
I ψ ¨ = P l s i n ψ + N l c o s ψ I\ddot{\psi}=Plsin\psi+Nlcos\psi Iψ¨=Plsinψ+Nlcosψ
观察上面的式子,你可能会发现里面少了一个分量,这个分量就是重力在垂直于摆方向的分力 m g s i n ψ mgsin\psi mgsinψ,很多博客和论文上也是直接这么写,没有解释原因。只有以质心为参考点时,重力不产生力矩,上式成立,但这显然是背离事实的,个人理解,这里在小角度时为了方便分析做了近似。
其中 I I I为摆的转动惯量。将 P P P N N N的表达式与力矩平衡方程联立,消去中间变量 P P P N N N,得到:
( I + m l 2 ) ψ ¨ − m g l s i n ψ = m l x ¨ c o s ψ (I+ml^2)\ddot{\psi}-mglsin\psi=ml\ddot{x}cos\psi (I+ml2)ψ¨mglsinψ=mlx¨cosψ

线性化

至此,我们通过受力分析得到了两个非常重要的式子:
( M + m ) x ¨ + b x ˙ − m l ψ ¨ c o s ψ + m l ψ ˙ 2 s i n ψ = F (M+m)\ddot{x}+b\dot{x}-ml\ddot{\psi}cos\psi+ml\dot{\psi}^2sin\psi=F (M+m)x¨+bx˙mlψ¨cosψ+mlψ˙2sinψ=F
( I + m l 2 ) ψ ¨ − m g l s i n ψ = m l x ¨ c o s ψ (I+ml^2)\ddot{\psi}-mglsin\psi=ml\ddot{x}cos\psi (I+ml2)ψ¨mglsinψ=mlx¨cosψ
考虑到倒单摆在实际工作时,偏转角 ψ \psi ψ通常比较小,于是有:
{ c o s ψ = 1 s i n ψ = ψ ψ ˙ = 0 \left\{ \begin{aligned} cos\psi & = & 1\\ sin \psi & = &\psi \\ \dot{\psi} & = & 0 \end{aligned} \right. cosψsinψψ˙===1ψ0
u u u来代表作用于受控对象的外力 F F F,结合上述近似结果,有:
{ ( M + m ) x ¨ + b x ˙ − m l ψ ¨ = u ( I + m l 2 ) ψ ¨ − m g l ψ = m l x ¨ \left\{ \begin{aligned} (M+m)\ddot{x}+b\dot{x}-ml\ddot{\psi}&=& u \\ (I+ml^2)\ddot{\psi}-mgl\psi &=& ml\ddot{x} \end{aligned} \right. {
(M+m)x¨+bx˙mlψ¨(I+ml2)ψ¨mglψ==umlx¨

求系统传递函数

由上一节,我们最终得到了一个关于系统状态的微分方程组。而拉普拉斯变换可以将微分方程转化为代数方程进行运算,使求解大为简化。

( M + m ) x ¨ + b x ˙ − m l ψ ¨ = u → ( M + m ) X ( s ) s 2 + b X ( s ) s − m l Ψ ( s ) s 2 = U ( s ) (M+m)\ddot{x}+b\dot{x}-ml\ddot{\psi}= u \rightarrow(M+m)X(s)s^2+bX(s)s-ml\Psi(s)s^2=U(s) (M+m)x¨+bx˙mlψ¨=u(M+m)X(s)s2+bX(s)smlΨ(s)s2=U(s)
( I + m l 2 ) ψ ¨ − m g l ψ = m l x ¨ → ( I + m l 2 ) Ψ ( s ) s 2 − m g l Ψ ( s ) = m l X ( s ) s 2 (I+ml^2)\ddot{\psi}-mgl\psi = ml\ddot{x}\rightarrow(I+ml^2)\Psi(s)s^2-mgl\Psi(s)=mlX(s)s^2 (I+ml2)ψ¨mglψ=mlx¨(I+ml2)Ψ(s)s2mglΨ(s)=mlX(s)s2
现在我们系统的输入变量是 U ( s ) U(s) U(s),而我们关心的是小车当前的位置 X ( s ) X(s) X(s)以及倒单摆的角度 Ψ ( s ) \Psi(s) Ψ(s)
经过整理,可以得到下面的系统传递函数。
摆角度的传递函数:
P p e n d ( s ) = ψ ( s ) U ( s ) = m l s q s 3 + b ( I + m l 2 ) s 2 − ( m + M ) m g l s − b m g l P_{pend}(s)=\frac{\psi(s)}{U(s)}=\frac{mls}{qs^3+b(I+ml^2)s^2-(m+M)mgls-bmgl} Ppend(s)=U(s)ψ(s)=qs3+b(I+ml2)s2(m+M)mglsbmglmls
小车位置的传递函数:
P c a r t ( s ) = X ( s ) U ( s ) = ( I + m l 2 ) s 2 − m g l q s 4 + b ( I + m l 2 ) s 3 − ( m + M ) m g l s 2 − b m g l s P_{cart}(s)=\frac{X(s)}{U(s)}=\frac{(I+ml^2)s^2-mgl}{qs^4+b(I+ml^2)s^3-(m+M)mgls^2-bmgls} Pcart(s)=U(s)X(s)=qs4+b(I+ml2)s3(m+M)mgls2bmgls(I+ml2)s2mgl
上面两式中 q q q为公因数项:
q = ( m + M ) ( I + m l 2 ) − m 2 l 2 q=(m+M)(I+ml^2)-m^2l^2 q=(m+M)(I+ml2)m2l2

版权声明:本文内容由互联网用户自发贡献,该文观点仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 举报,一经查实,本站将立刻删除。

发布者:全栈程序员-用户IM,转载请注明出处:https://javaforall.cn/171030.html原文链接:https://javaforall.cn

【正版授权,激活自己账号】: Jetbrains全家桶Ide使用,1年售后保障,每天仅需1毛

【官方授权 正版激活】: 官方授权 正版激活 支持Jetbrains家族下所有IDE 使用个人JB账号...

(0)
blank

相关推荐

  • Java之父:C语言是撑起一切的基石[通俗易懂]

    Java之父:C语言是撑起一切的基石[通俗易懂]为了忘却的纪念:无论任何领域,伟大的人物,永远不死!2011年10月9日,患有前列腺癌和心脏病的C语言之父丹尼斯·里奇离世。不知不觉,这位伟大的程序员已经离开我们已经有6年了。有人曾说,没有C语言之父,就没有乔布斯。C语言之父影响力其实不止这些。毫不夸张的说,没有C语言之父,就没有微软的Windows10和SurfaceBook,也没有安卓智能手机

  • assertEquals 方法「建议收藏」

    assertEquals 方法「建议收藏」assertEquals  函数原型1:assertEquals([Stringmessage],expected,actual)参数说明:message是个可选的消息,假如提供,将会在发生错误时报告这个消息。  expected是期望值,通常都是用户指定的内容。actual是被测试的代码返回的实际值。  函数原型2:assertEquals([Stringmessa…

  • (建议收藏)TCP协议灵魂之问,巩固你的网路底层基础

    (建议收藏)TCP协议灵魂之问,巩固你的网路底层基础

  • Linux 安装 rabbitmq 3.7.16 详细教程[通俗易懂]

    前提:下载依赖,安装rabbitmq需要下载以下依赖,这里可以提前下载上。依次执行以下命令安装依赖类库:#yum-yinstallmakegccgcc-c++kernel-develm4ncurses-developenssl-devel#yuminstallxmlto-y执行成功后如下(complete)一、安装erlang(最新版本22….

  • java 除法取整_java 除法运算只保留整数位的4种方式

    java 除法取整_java 除法运算只保留整数位的4种方式1.情景展示根据提供的毫秒数进行除法运算,如果将毫秒数转换成小时,小时数不为0,则只取整数位,依此类推…2.情况分析可以使用3个函数实现Math.floor(num)  只保留整数位Math.rint(num)  余数四舍五入Math.ceil(num)  取整位,再+1举例:doublenum=3.1415926;System.out.println(Math.floor…

  • EJB学习笔记

    EJB学习笔记Chapter01企业级开发背景知识一、什么是企业级程序(EnterpriseApplication)? 具有以下特点的程序: (1)围绕商业目的; (2)分布分层的程序架构。二、企业级应用的架构发展历史Host/Terminal(主机/终端,终端不具备处理能力,数据由主机处理)终端为哑终端,可接受命令,不处理命令优点:安全(不能直接访问)

发表回复

您的电子邮箱地址不会被公开。

关注全栈程序员社区公众号