将机器人整体开源，同时总结一下机器人搭建过程中遇到的坑和未来的改进方向。在分享的文件里包含了结构设计、程序控制、电路设计以及其他模块相关资料供大家参考。

    文末扫码可以获取相关资料。

机器人原理分析

    首先来看成品图：

, `* m" q4 M  c; T. E0 ]- t

    如图所示，该机器人根据陀螺仪的位姿数据，通过三个全向轮驱动底部球体调整自己在球上的位置，保持动态平衡的同时实现全向移动。

    保持动态平衡过程需要对机器人进行运动学分析，这里参考了平衡小车之家的运动学方程：

    自平衡控制问题转化为三步：输入X、Y角度—控制器计算—输出A、B、C电机转速的控制模型。

控制器设计

    首先考虑参考平衡车控制，球上自平衡机器人本质上依然是一个一阶倒立摆问题。

    这里参考了飞思卡尔直立车的控制方法，采用串级PID控制器，外环PD角度环，内环速度PI环。相关文章：PID算法原理介绍。

    由于我的驱动方案选择的是42步进电机，在速度闭环的时候有些问题。正常的直流电机+编码器的控制方案可以通过编码器将轮子的真实速度计算出来，从而和控制器的理想转速作差，实现速度控制。

    而我这里的速度闭环是通过计算上一个时钟周期时给步进电机的控制量，通过运动学方程分解，得到机器人的虚拟速度，与理想转速作差控制。我认为这种速度闭环方式还是存在一定缺陷的，但是在网上查看论文的时候我发现有很多自平衡机器人都是用42步进电机来实现速度闭环的，不知道是什么方法。

    这里还可以好好思考一下为什么角度环要用PD控制，速度环要PI控制，角度环的P部分和D部分对机器人控制有什么影响？在很多CSDN调试平衡车的博客中都有解释，这里就留给大家思考了。

硬件及结构设计

    自平衡机器人的硬件清单有：

56mm全向轮 45元/个

42步进电机 25/个

42步进闭环模块 59.8元/个

LM2596S降压模块 20元

STM32F103C8T6-4飞控板 59.8元

GY-521六轴陀螺仪 25元

    用到的模块大致如上所示，C8T6的价格随着最近芯片涨价直线上升，我白嫖了实验室的两块板子，现在买一块实在太贵，可以等芯片价格稳定一些再买。其余开关排针等常见元件不再赘述。
) I0 s6 W2 x8 j0 `

    电路原理图如下所示：
2 j9 e: k) u; h7 g% f- U  V

    机器人使用solidworks设计整体结构，底板可在某宝定制6050太空铝切割，蓝色件为正常3D打印件。

程序部分

    在keil 5中开发STM32。

    控制程序采用定时器0.5ms定时中断的方式进行计算，每触发两次中断计算对电机控制一次，这里还是推荐大家采用外部中断读取GY-521上的INT引脚的方式，控制计算周期。GY-521上的INT引脚每5ms触发一次跳变，采用外部中断的方式可以严格保证读取位姿数据与计算处理同步。

# W$ {. c$ U* R  j1 o* B  S6 [2 v2 p

1. int TIM1_UP_IRQHandler (void)               
   
2. {
   
3.     u8 key_cal;
   ; i& \/ A" k1 t/ L
4.     if(TIM_GetITStatus(TIM1,TIM_IT_Update) != RESET)
   
5.     {
   
6.         TIM_ClearITPendingBit(TIM1,TIM_IT_Update);
   % P# S" {9 f# o
7.         flag_target=!flag_target;
   + B0 O; C+ p9 }6 c3 A4 j" @! H
8.         key_cal=KEY_Scan(0);
   ) I" z9 k6 Z! r. S
9.         if(state_flag==1)//矫正结束
   
10.         {
    
11.             if(flag_target==1)//每读取两次陀螺仪控制一次
    
12.             {
    
13.                         Read_DMP();                      //===读取倾角
    
14.                         scope();
    
15.                         return 0;
    
16.             }
    1 V# u6 a. J4 c, W
17.         }
      V' T% s# r7 h( _9 L
18.         if(key_cal==1)//矫正按键
    8 |: ?7 x; D, s7 h) ^, q/ x
19.         {
    ' w8 d" x/ B/ P; x
20.                         Angle_Zero_X=Angle_Balance_X;
    ! C5 V  Q# @7 W! o; Z7 ~: ^
21.                         Angle_Zero_Y=Angle_Balance_Y;
      `1 Y7 j, Y. N- \* y
22.                         key_cal=0;
    ( B3 C* H# @/ \& t6 n
23.                   Flag_Stop=0;
    
24.         }
    
25.         if(key_cal==2||key_cal==3)//矫正按键
    
26.         {
    
27.         Flag_Stop=1;//关闭速度环I积分
    " R, c3 @' L: q, g
28.         key_cal=0;
    
29.         }
    % B) s& f* C1 c: o/ `" Q" Y
30.         Angle_Bias_X =Angle_Balance_X-Angle_Zero_X;     //获取Y方向的偏差
    
31.         Angle_Bias_Y =Angle_Balance_Y-Angle_Zero_Y;     //获取Y方向的偏差
    
32.         if(control_mode==0)//PID控制模式
    
33.         {
    
34.         Encoder_Analysis(Motor_A,Motor_B,Motor_C);  //正运动学分析，得到X Y方向的速度
    * e+ D& E) T) \, Y) D4 s
35.         Balance_Pwm_X= balance_X(Angle_Bias_X,Gyro_Balance_X);//X方向的倾角控制
    
36.         Balance_Pwm_Y=-balance_Y(Angle_Bias_Y,Gyro_Balance_Y);  //Y方向的倾角控制
      s8 P7 [9 a/ U4 a4 `
37. //      if(++flag_target_2==4)//速度环频率慢于加速度环 但是还没加速度环
    ; p! V3 t4 s- }( x9 t* R
38. //          {
    
39.             Velocity_Pwm_X=velocity_X(compute_X);      //X方向的速度控制
    
40.             Velocity_Pwm_Y=velocity_Y(compute_Y);     //Y方向的速度控制  
    
41. //          flag_target_2=0;
    
42. //          }
    
43.             Move_X =Balance_Pwm_X+Velocity_Pwm_X;   //===X方向控制量累加   
    - {- p4 q/ o" G% C  @% T
44.             Move_Y =Balance_Pwm_Y+Velocity_Pwm_Y;   //===Y方向控制量累加   
    
45.             Move_Z=0;               
    7 S, D4 g4 {0 B
46.             Kinematic_Analysis(Move_X,Move_Y,Move_Z);//逆运动学分析得到ABC电机控制量
    
47.         }
    6 ~3 a( A- f1 d
48.             Motor_A=Target_A;//直接调节PWM频率
    
49.             Motor_B=Target_B;//直接调节PWM
    
50.             Motor_C=Target_C;//直接调节PWM
    
51. //以下都是为了速度连续化处理防止突变
    
52.             if(Motor_A==0)  Motor_A=motor_a_last;
    
53.             if(Motor_B==0)  Motor_B=motor_b_last;
    
54.             if(Motor_C==0)  Motor_C=motor_c_last;
    9 N! Z9 W* q$ \$ b
55.             Xianfu_Pwm(2000);
    9 s8 w0 g% F, n$ i
56.             Set_Pwm(Motor_A,Motor_B,Motor_C);
    / W5 _4 e& A4 e) x# |
57.             Gyro_Balance_X_last=Gyro_Balance_X;
    
58.             Gyro_Balance_Y_last=Gyro_Balance_Y;
    
59.             Gyro_Balance_Z_last=Gyro_Balance_Z;
    6 `0 e7 d7 z1 p2 c$ S) e, A
60.             Angle_Balance_X_last=Angle_Balance_X;
    
61.             Angle_Balance_Y_last=Angle_Balance_Y;
    
62.             Angle_Balance_Z_last=Angle_Balance_Z;
    % g. b: U; Z7 i8 Q3 [4 e
63.             motor_a_last=Motor_A;
    # W3 W. b; i/ u- n( o: `# s
64.             motor_b_last=Motor_B;
    6 r2 }! n; N! o* e8 V1 e" s
65.             motor_c_last=Motor_C;
    / B6 o* R& F3 H0 R' L
66.     }
    
67.       return 0;
    
68. }

复制代码

5 Z# {8 s9 v' Z* r7 G/ t对于电机控制，由于采用的驱动方案是步进电机，调速的方式是改变驱动步进电机的脉冲频率。我这里选择了三个定时器，动态调节定时器的频率，具体方式是在初始化时设定好定时器的预分频系数psc的值，然后在程序里动态更改ARR寄存器的值，从而改变定时器的定时频率。
7 P' e9 c$ s9 A1 `  q

1. 
   ! M" m& _! ?, `+ H, K
2. //这里以A电机的速度控制为例 输入为 电机方向和电机速度
   7 f9 T! \3 K9 G! R! Z! Z
3. void set_motorA_speed(u8 dir,u16 speed)
   
4. {
   & B( x! t$ k  u! `/ I
5.     u32 arr;
   * D) `/ h$ y& q! s' ]8 ]( v% I% A; Z
6.     arr=speed;
   - y% \  V6 X* q/ v, J9 U: x( d; X
7.         TIM_ARRPreloadConfig(TIM3,DISABLE);
   8 v" {' W1 ^* \4 X
8.         TIM3->ARR=arr;//计数到10000在归零重新计数
   ( o% t- A/ R; @# Q, t* k+ L
9.         TIM3->CCR4=arr/2;//保持占空比为50%
   
10.         TIM_ARRPreloadConfig(TIM3,ENABLE);
    & s, J1 P6 }' J3 Y9 B
11.         TIM_Cmd(TIM3,ENABLE);
    
12.     if(dir==0)
    
13.     {
    
14.     GPIO_SetBits(GPIOA,GPIO_Pin_1);
    ! J$ V4 U3 ~2 F& L7 w# t
15.     }
    & b/ P  r, R2 q
16.   else
    7 X3 g, T* X% s2 A* f: N: o4 T; z( g
17.   {
    
18.     GPIO_ResetBits(GPIOA,GPIO_Pin_1);
    
19.   }
    ( `/ d  R& N6 K4 n
20. }

复制代码

    小车的运动学分解代码实现如下，参考了平衡小车之家的代码：

1 F0 ?# |5 k2 @

1. 
   ' z, y/ k# Q7 B, e* M
2. /**********************************************************
   + ^# v+ ~% \' f! `) y& G
3. 函数功能：小车运动数学模型
   ; H+ m5 h' b  _- K% O* B
4. 入口参数：X Y Z 三轴速度或者位置
   
5. 返回  值：无
   5 m% _  B$ S" C4 u% H  z
6. ***********************************************************/
   
7. void Kinematic_Analysis(float Vx,float Vy,float Vz)
   8 n" v5 o# j  j% `$ u7 k
8. {
   
9.         Target_A   = Vx + L_PARAMETER*Vz;
   3 j5 X4 f& M% m- f) x
10.         Target_B   = -X_PARAMETER*Vx + Y_PARAMETER*Vy + L_PARAMETER*Vz;
    ( t& ~4 S1 S6 S+ t3 U
11.         Target_C   = -X_PARAMETER*Vx - Y_PARAMETER*Vy + L_PARAMETER*Vz;
      N' o# W* v  v  V( j  |% q7 t6 y: G
12. }
    % v+ l" A- |( W$ e/ {5 O! V+ S% `
13. /*****************************************************************
    
14. 函数功能：小车运动 正运动学分析
    ; `( d, _  m, _# |; \, e
15. 入口参数：A B C三个电机的速度
    * H0 v% h9 r% ~6 N  l1 h
16. 返回  值：无
    
17. ******************************************************************/
    
18. void Encoder_Analysis(float Va,float Vb,float Vc)
    
19. {
    
20.     compute_X=(Va*2-Vb-Vc);
    ! D' Z( c7 U# `, R( ]+ T
21.     compute_Y=((Vb-Vc)*sqrt(3));
    
22.     compute_Z=(Va+Vb+Vc);  
    
23. }

复制代码

) p5 C% P1 V: a6 \, Z

总结与展望

    球上自平衡机器人可以作为算法试验平台， 输入输出固定，更换不同控制器，将数据导入MATLAB进行分析即可比较控制器性能。

    个人认为结构有两个改进方向，一方面参考以下论文：余义. 单球驱动自平衡机器人位姿解算与控制系统研究[D].武汉科技大学,2019。论文中采用的四足式驱动结构更有利于机器人自平衡控制。
' Y: k# E. L" d3 B2 R# {7 g

    另一方面可以增加球体和机器人固定装置，利用机械结构将机器人与底部驱动球结合成一个整体防止机器人跳轮等问题。同时驱动球对于机器人平衡的影响较大，最好还是定制空心钢球，然后喷漆增大摩擦力，最有利于机器人自平衡控制。
7 W% X, s$ O% r/ _/ R7 h  J" Q

    控制部分的改进，首先是控制原理，本文是针对建立好的运动学方程进行分析，通过串级PID算法来实现自平衡运动。该机器人的控制问题本质上是一阶倒立摆问题，可以采用动力学建模的方式，通过动力学分析算出平衡需要的虚拟力矩，再对电机进行力矩控制。相关推荐：演示PID三个参数的控制作用。

    其次是控制器，PID控制算法应用广泛但也有一定的缺点，可以考虑采用模糊PID，ADRC自抗扰控制器，强化学习等智能控制算法对机器人自平衡进行控制。

) L  ^! \$ i, {8 a! ^+ E
: Z* d3 Q! ^# `; A' _5 V4 N