当前位置: 首页 > news >正文

倒立摆建模

前言

系统由一辆具有动力的小车和安装在小车上的倒立摆组成,系统是不稳定,我们需要通过控制移动小车使得倒立摆保持平衡。

具体地,考虑二维情形如下图,控制力为水平力FFF,输出为角度θ\thetaθ以及小车的位置xxx

在这里插入图片描述

力分析和系统方程

设两个物体在水平和竖直方向上的相互作用力分别为NNNPPP
将小车水平方向上的力相加,得到牛顿力学方程:
Mx¨+bx˙+N=FM\ddot{x}+b\dot{x}+N=FMx¨+bx˙+N=F
将倒立摆水平方向上的力相加,得到:
N=mx¨+mlθ¨cosθ−mlθ˙2sinθN=m\ddot{x}+ml\ddot{\theta}cos\theta-ml\dot{\theta}^2sin\thetaN=mx¨+mlθ¨cosθmlθ˙2sinθ
NNN消去,得
Mx¨+bx˙+mx¨+mlθ¨cosθ−mlθ˙2sinθ=FM\ddot{x}+b\dot{x}+m\ddot{x}+ml\ddot{\theta}cos\theta-ml\dot{\theta}^2sin\theta=FMx¨+bx˙+mx¨+mlθ¨cosθmlθ˙2sinθ=F
将垂直于摆的力相加,得
Psinθ+Ncosθ−mgsinθ=mlθ¨+mx¨cosθ(a)Psin\theta+Ncos\theta-mgsin\theta=ml\ddot{\theta}+m\ddot{x}cos\theta(a)Psinθ+Ncosθmgsinθ=mlθ¨+mx¨cosθ(a)
为了消去P,NP,NP,N两项,将两个物体质心的力矩相加,得
Plsinθ−Nlcosθ=Iθ¨(b)Plsin\theta-Nlcos\theta=I\ddot{\theta}(b)PlsinθNlcosθ=Iθ¨(b)
(a)×l+(b)(a)\times l+(b)(a)×l+(b)
(I+ml2)θ¨+mglsinθ=−mlx¨cosθ(I+ml^2)\ddot{\theta}+mglsin\theta=-ml\ddot{x}cos\theta(I+ml2)θ¨+mglsinθ=mlx¨cosθ

线性化

将上述方程线性化,θ=π\theta=\piθ=π,并假设系统保持在这个平衡的小邻域内。设ϕ\phiϕ为摆对平衡位置的偏差,满足θ=π+ϕ\theta=\pi+\phiθ=π+ϕ,若偏差十分小,使用以下近似:
cosθ=cos(π+ϕ)≈−1sinθ=sin(π+ϕ)≈−ϕθ¨2=ϕ˙2≈0cos\theta=cos(\pi+\phi)\approx-1 \\ sin\theta=sin(\pi+\phi)\approx-\phi \\ \ddot{\theta}^2=\dot{\phi}^2\approx0cosθ=cos(π+ϕ)1sinθ=sin(π+ϕ)ϕθ¨2=ϕ˙20

将上述近似公式应用至前面的非线性控制方程,得到两个线性化的方程,并使用uuu替代FFF
(I+ml2)ϕ¨−mglϕ=mlx¨(M+m)x¨+bx˙−mlϕ¨=u(I+ml^2)\ddot{\phi}-mgl\phi=ml\ddot{x} \\ (M+m)\ddot{x}+b\dot{x}-ml\ddot{\phi}=u(I+ml2)ϕ¨mg=mlx¨(M+m)x¨+bx˙mlϕ¨=u

传递函数

假设初始条件为0,对系统方程应用拉普拉斯变换:

todo

http://www.lryc.cn/news/33592.html

相关文章:

  • SpringSecurity支持WebAuthn认证
  • 深度学习技巧应用3-神经网络中的超参数搜索
  • 【信号量机制及应用】
  • 围棋高手郭广昌的“假眼”棋局
  • 学成教育-统一异常处理实现
  • JNI内通过参数形式从C/C++中传递string类型数据至Java层
  • 自动化测试——执行javaScript脚本
  • 常用十种算法滤波
  • IO多路复用
  • Python中的错误是什么,Python中有哪些错误
  • 记录自己开发一款小程序中所遇到的问题(uniapp+uview)(持续更新)
  • 华为机试 HJ43 迷宫问题
  • 数据结构|链表
  • 计算机写论文时,怎么引用文献? - 易智编译EaseEditing
  • 实验三:贪心
  • MySQL日志文件
  • Intel8086处理器使用NASM汇编语言实现操作系统08-关于负数的相关处理idiv/cbw/cwde/cdqu/cwd/cdq/cdo/
  • JavaScript 混淆技术
  • 安装库报错:No CUDA runtime is found, using CUDA_HOME=‘/usr/local/cuda-11.3‘
  • CVTE前端面经(2023)
  • 基于EB工具的TC3xx_MCAL配置开发02_ICU模块配置
  • jmeter高阶系列--beanshell返回值中提取参数
  • 面向对象
  • mpi4py 运行过程中出现Read -1, expected xxx, errno = 1 解决方案
  • PMP考前冲刺3.07 | 2023新征程,一举拿证
  • 60条Python日常工作中的高频写法,收藏
  • (小甲鱼python)函数笔记合集七 函数(XI)总结 python函数的函数文档、类型注释、内省详解
  • Leetcode是什么
  • 2023-03-07 MySQL—基于规则优化-子查询优化
  • Rocketmq技术详解