一种桥梁拉索检测机器人吸力控制方法与流程

[0001]
本发明涉及一种磁吸式拉索爬升机器人的吸力控制方法，属于机器人控制技术领域。


背景技术：

[0002]
机器人在壁面上的吸附方式有多种，例如负压吸附，真空吸附，磁吸附等等，由于斜拉索一般都有导磁性物质设计，因此使用磁吸附设计吸附结构成本最低，且效果最好。磁吸附方式按照其吸附单元又分为电磁吸附和永磁吸附两种。
[0003]
永磁铁在任何情况下都有磁性，其与磁吸性物质表面接触在一起就可以产生磁力，然而磁力的大小是固定不变的，是一个定值。若遇到前文分析的恶劣天气时，机器人的机身发生抖动，却无法及时调整吸附力的大小，磁吸力不足以将机器人稳定在斜拉索上，则可能会出现坠机的情况，不仅会造成经济上的损失，还会出现很危险的情况。


技术实现要素：

[0004]
发明目的：为了克服现有技术中存在的不足，本发明提供一种桥梁拉索检测机器人吸力控制方法，可以及时的控制磁吸附力的大小，使用电磁铁来设计吸附结构，电磁铁磁力便于控制，控制通入电磁铁线圈电流的大小即可控制磁吸力的值。
[0005]
技术方案：为实现上述目的，本发明采用的技术方案为：
[0006]
一种桥梁拉索检测机器人吸力控制方法，包括以下步骤：
[0007]
步骤1，根据电流大小、线圈匝数、磁极与斜拉索的相对面积建立电磁力模型；
[0008]
步骤2，步骤2中根据电磁铁曲率半径得到磁极处曲线的数学方程为如下式：
[0009][0010]
根据绳索半径得到绳索的绳索的数学方程如下：
[0011][0012]
将磁极处曲线的数学方程和绳索的数学方程相减后积分就能得到气隙的总长度：
[0013]
隙的总长度为：
[0014][0015]
其中，表示气隙总长度，x表示直角坐标系横坐标的值，y
1
表示直角坐标系中磁极曲线纵坐标的值，r
1
表示磁极曲率半径，y
2
表示直角坐标系中绳索曲线纵坐标的值，r
2
表示绳索半径大小；
[0016]
步骤3，根据线圈匝数、铁心横截面积、饱和磁通量密度以及线圈电感值计算磁饱和电流和电感值：
[0017]
步骤4，机器人静止时电磁力的最小值：
[0018][0019]
其中，f
磁min
表示静止时电磁力的最小值，mg为整个机器人的重力，α表示斜拉索的倾斜角，θ表示机器人两个轮臂之间的夹角；
[0020]
步骤5，将机器人的振动作为干扰，使用震动传感器检测机器人的抖动状态；磁吸附机器人在斜拉索上运动，若出现振动，则振动传感器会检测到振动信号，并且用arduino单片机读取振动信号值，再根据信号大小，调整输出端口电压值的大小；最后利用升压模块将单片机输出端口的电压值放大，通入电磁铁两端；
[0021]
步骤6，采用位置式pid控制来控制单片机输出端口的电压从而控制通入电磁铁线圈的电流：
[0022][0023]
其中，u(t)表示pid控制输出值，k
p
表示比例控制系数，e(t)表示实时误差值，k
i
表示积分控制系数，k
d
表示微分控制系数，e表示误差值，setpoint表示输出值，input表示输入值，t表示当前时间；
[0024]
步骤7，直接增大电磁力的pid控制，直接增大电磁力就是将电压由某一值瞬间跳转到下一个值，此时电压的变化类似一个阶跃变化模型，单位阶跃函数模型：
[0025][0026]
其中，y(t)表示阶跃函数值，t表示当前时间，当t<0时，函数值为0；当t>0时，函数值为1；在t＝0这段时间内，函数值由0变为1；
[0027]
步骤8，得到最基本的pid控制的基本传递函数：
[0028][0029]
其中，c(s)表示传递函数值，s表示复频域系数；设置pid控制参数k
p
、k
i
、k
d
，将电压变化的阶跃函数作为输入，pid控制函数作为传递函数；
[0030]
步骤9，针对正弦变化抖动的pid控制，当机器人在斜拉索上发生正弦抖动时，振动传感器将会检测到一条呈正弦变化规律的振动等级变化曲线，因此跟随抖动变化，单片机端口的输出电压也会随之产生正弦变化：
[0031][0032]
其中，u(t)表示电压值，a表示偏移，b表示幅度值，t表示振动发生时开始计时的时间，t
0
表示开始振动的时间
[0033]
对函数进行拉氏变换得：
[0034][0035]
其中，u(s)表示复频域函数值，s表示复频域参数；
[0036]
步骤10，针对不规则变化抖动的pid控制，设置一个随机函数的信号源，将生成的电压变化随机函数作为信号源，使用pid控制。
[0037]
优选的，步骤1中的电磁力模型：
[0038][0039]
其中，f:电磁力大小，i表示电流大小，n表示线圈匝数，表示保险系数，其值在0.05～0.15之间，μ
0
表示空气的导磁系数，其值为4π
×
10-7
h/m，s表示磁极与斜拉索的相对面积，δ表示气隙平均长度；。
[0040]
优选的：步骤2中根据电磁铁曲率半径得到磁极处曲线的数学方程为如下式：
[0041][0042]
根据绳索半径得到绳索的绳索的数学方程如下：
[0043][0044]
其中，x表示直角坐标系横坐标，y
1
表示直角坐标系磁极曲线纵坐标，r
1
表示磁极曲率半径，y
2
表示直角坐标系绳索曲线纵坐标，r
2
表示绳索半径。
[0045]
优选的：步骤3中磁饱和电流和电感值公式如下：
[0046][0047][0048]
式中，l表示螺旋线圈的电感值，l表示螺旋线圈的长度约，n表示线圈匝数，s表示线圈横截面积，μ表示线圈内部磁芯导磁率，k表示长冈系数值；，i
sat
表示磁饱和电流值，ae表示铁心横截面积，b表示饱和磁通量密度大。
[0049]
本发明相比现有技术，具有以下有益效果：
[0050]
本发明机器人能根据检测到的振动信号，自动选择效果最好的pid参数进行控制。当pid控制参数合适时，电磁力的控制效果可以满足使用。
附图说明
[0051]
图1为本发明机器人的结构示意图；
[0052]
图2为曲率半径直角坐标系；
[0053]
图3为绳索半径r
2
与气隙间隙关系；
[0054]
图4为阶跃变化电压仿真曲线图；
[0055]
图5为调整pid参数后电压变化曲线图；
[0056]
图6为电压变化曲线；
[0057]
图7为pid控制后电压正弦变化仿真曲线；
[0058]
图8为调整pid参数后电压仿真曲线；
[0059]
图9为随机生成无规律电压变化曲线图；
[0060]
图10为pid控制电压变化曲线。
具体实施方式
[0061]
下面结合附图和具体实施例，进一步阐明本发明，应理解这些实例仅用于说明本发明而不用于限制本发明的范围，在阅读了本发明之后，本领域技术人员对本发明的各种等价形式的修改均落于本申请所附权利要求所限定的范围。
[0062]
一种桥梁拉索检测机器人吸力控制方法，机器人为桥梁拉索检测机器人，如图1所示，本机器人包括伸缩结构1、扭转结构2、电磁吸附结构3、驱动电机4、驱动轮5以及支架6，所述伸缩结构1一端安装在支架6上，另一端与电磁吸附结构3固定连接，所述驱动电机4、驱动轮5安装在扭转结构2上，所述扭转结构2安装在支架6上，所述驱动电机4与驱动轮5驱动连接。
[0063]
所述伸缩结构1包括伸缩臂11、调节螺母一12、调节螺母一13，所述伸缩臂11顶部设置有外螺纹，底部与电磁吸附结构3固定连接，所述伸缩臂11与支架6滑动连接，所述调节螺母一12、调节螺母一13通过外螺纹与伸缩臂11固定连接，且所述支架6位于调节螺母一12与调节螺母一13之间。
[0064]
机器人需要吸附在斜拉索上，吸附面为曲面，机器人整体结构可设计为v型结构，即机器人的两个轮臂之间夹成一个v型，且驱动轮与轮臂垂直。所述扭转结构2包括轮臂21、旋转轴和定位螺丝22，所述轮臂21通过旋转轴与支架6转动连接，所述定位螺丝22用于将支架6与旋转轴固定连接。
[0065]
所述电磁吸附结构3背离伸缩结构1的另一端面为弧形面。所述电磁吸附结构3为电磁铁。为了确保电磁结构提供的吸附力能均匀分布到机器人的每个部位，将电磁结构设计在机器人的机身正中央，且机身底部靠近头部和尾部部分分别安放一个电磁结构。
[0066]
为了保证足够的驱动力，机器人设计成四轮驱动，即每个轮胎处都安装一个驱动电机，每个驱动轮都单独控制，因此即便机器人的某一个驱动轮失效，其他的三个驱动轮也能提供足够的动力，确保机器人的运动。
[0067]
当需要调整电磁吸附结构3和缆索位置时，先把调节螺母一12与调节螺母一13拧松，调整伸缩臂11，使电磁吸附结构3与缆索贴近，再把调节螺母一12与调节螺母一13拧紧即可。当需要调整两个轮臂之间的角度时，拧松定位螺丝22，把轮臂21之间的角度调整到需要的角度时，再将定位螺丝22拧紧。为了让吸附轮与拉索面垂直，可调整轮臂与机器人机身的角度。电磁吸附结构3提供电磁力，让将机器人能稳定吸附在斜拉索面上。驱动电机和驱动轮为机器人的运动提供动力。
[0068]
本发明运用的就是作为起重电磁铁部分的功能，将电磁铁放在机器人底端，和机器人连接在一起。把斜拉索看成铁磁性材料，当电磁铁导线通电时，即可产生电磁力将机器人和斜拉索吸附在一起，以实现吸附目的。
[0069]
由于在运动和静止时，机器人在垂直于斜面方向受力分析不变，因此我们分析在静止时的受力分析即可。
[0070]
由于在圆柱面上各切面支持力处处相等，因此有：
[0071]
n
0
＝n
1
+n
2
＝n
3
+n
4
ꢀꢀꢀꢀꢀ
(2-1)
[0072]
n
1
＝n
2
＝n
3
＝n
4
ꢀꢀꢀꢀꢀꢀꢀ
(2-2)
[0073]
此时对垂直于斜面进行动力学分析得到：
[0074][0075]
式(2-3)为电磁力和机器人收到的重力，与斜面对机器人支持力间的关系。
[0076]
机器人在平行于斜面方向上受到重力的分力，轮胎提供的静摩擦力。由于电磁结构与斜面接触是类似一条直线，接触面可忽略不计，因此摩擦力也可忽略。因此有如下关系式：
[0077]
mgcosα＝f
1
+f
2
+f
3
+f
4
ꢀꢀꢀꢀꢀꢀꢀ
(2-4)
[0078]
因为轮胎的材质一样，设其摩擦系数为μ则：
[0079]
f
1
＝f
2
＝f
3
＝f
4
＝μn
1
ꢀꢀꢀꢀꢀꢀꢀ
(2-5)
[0080]
结合式(2-1),(2-4),(2-5)得：
[0081]
mgcosα＝2μn
0
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(2-6)
[0082]
结合(2-1),(2-3)得：
[0083][0084]
因此根据式(2-6)式(2-7)我们可以得出静止时电磁力的最小值
[0085][0086]
由式(2-8)可以看出，当轮胎夹角和绳索倾斜角固定时，电磁力最小值和机器人重力成正相关。
[0087]
下面进行线圈匝数和电阻的推导：
[0088]
线圈纵向匝数为：
[0089][0090]
线圈横向匝数为：
[0091][0092]
电磁铁宽度为：
[0093]
d＝d
0
+n
横
×
l
0
×
2
ꢀꢀꢀꢀꢀꢀꢀꢀ
(3-3)
[0094]
线圈平均匝长为：
[0095][0096]
线圈电阻为：
[0097]
r＝ρln
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(3-5)
[0098]
其中l为铁心长度，l
0
为线圈直径，d为电磁铁宽度，d
0
为铁心直径，l为线圈平均匝长，r为线圈电阻，n为线圈总匝数。
[0099]
磁吸附机器人通过电磁力将机器人吸附在斜拉索上，控制机器人在缆索上的运动最重要的就是对电磁力的控制。因此首先需要针对设计的电磁结构进行电磁力的数学建模。
[0100]
麦克斯韦方程组是电磁学中的最基本方程。它由四个方程组成，分别是：麦克斯韦-安培定律、高斯磁定律、高斯定律和法拉第感应定律、。利用麦克斯韦吸力公式建立吸力
的数学模型。
[0101]
从吸力公式可知：
[0102][0103]
式中：
[0104]
f：电磁吸力；
[0105]
b：磁感应强度；
[0106]
s：磁极的接触表面积；
[0107]
μ
0
：空气的导磁系数，其值为4π
×
10-7
h/m；
[0108]
由于电磁铁与斜拉索间的相对面积和空气导磁系数是定值，因此磁吸力大小和磁感应强度成正相关。下面分析电流和磁感应强度间的关系，根据磁势公式可得：
[0109]
in＝(in)
δ
+(in)
cm
+(in)
h
ꢀꢀꢀꢀꢀꢀꢀꢀ
(4-2)
[0110]
式中：
[0111]
(in)
δ
：消耗在气隙中的磁势；
[0112]
(in)
cm
：消耗在铁心磁阻上的磁势；
[0113]
(in)
h
：消耗在非工作气隙上的磁势；
[0114]
其中(in)
cm
和(in)
h
的大小相对(in)
δ
可忽略不计，但是保险起见，取其总和为(in)
δ
的倍。
[0115]
(in)
δ
和磁感应强度大小之间关系如下：
[0116][0117]
代入保险系数得到磁感应强度与电流关系得：
[0118][0119]
其中δ为平均气隙长度。
[0120]
再结合式(4-4)和式(4-1)得：
[0121][0122]
式(4-5)为电磁力模型，式中：
[0123]
f:电磁力大小；
[0124]
i：电流大小；
[0125]
n：线圈匝数；
[0126]
保险系数，其值在0.05～0.15之间；
[0127]
μ
0
：空气的导磁系数，其值为4π
×
10-7
h/m；
[0128]
s：磁极与斜拉索的相对面积；
[0129]
δ：气隙平均长度；
[0130]
由式(4-5)可知，电磁力大小和电流，绕铁心线圈匝数，电磁铁与斜拉索接触面积，气隙长度有关。其中线圈匝数，电磁铁与斜拉索之间接触面积是固定的，而气隙长度与斜拉
索的直径和电磁铁磁极形状有关，下面分析斜拉索直径和平均气隙的关系。
[0131]
电磁铁长度为50mm，且其曲率半径为45mm。而斜拉索的直径在60mm到180mm间，下面分析气隙和斜拉索直径的关系。建立如下图2所示的直角坐标系。
[0132]
磁极曲率半径为45mm，因此当缆索的半径大于45mm时，磁极的两端和缆索接触。当缆索半径小于45,mm时，磁极的中央与缆索接触。
[0133]
磁极的形状为弧形，设其曲率半径为r
1
，绳索半径为r
2
。磁极处曲线的数学方程为如下式：
[0134][0135]
绳索的数学方程如下：
[0136][0137]
因为气隙的长度处处不相同，所以我们需要先将气隙的总长度求出来。将磁极处曲线的数学方程和绳索的数学方程相减后积分就能得到气隙的总长度。因此气隙的总长度为：
[0138][0139]
将y
1
,y
2
代入其中得：
[0140][0141]
将得到的气隙总长度除以磁极的长度则可以得到气隙的平均长度。根据条件我们知道磁极曲率半径r
1
＝45mm，而绳索的半径：30mm<r
2
<90mm。得到拉索与机器人之间的平均间隙，随着斜拉索半径变化的曲线如下图3所示。
[0142]
由变化曲线可知，当拉索曲率半径小于45mm即直径小于90mm时，平均气隙随直径增大而减小，当拉索曲率半径大于45mm即直径大于90mm时，平均气隙随直径增大而增大。当拉索曲率半径和磁极曲率半径相同时，能完全贴合在一起，气隙为0。因此当机器人在半径不同的斜拉索上运动时，其电磁力的模型会根据斜拉索的直径变化而发生变化。
[0143]
由前面的分析我们得到了电磁力的模型如式(4-5)所示：
[0144][0145]
我们取保险系数0.1，根据实物得知线圈匝数为3050匝，电磁铁与斜拉索接触面积为1500mm
2
。则电磁力模型如下式(4-10)所示。
[0146][0147]
取气隙δ＝0.2mm则电磁力大小如下式(4-11)所示
[0148][0149]
磁饱和现象即亚铁磁性物质和铁磁性物质处于磁极化强度或磁化强度不随磁场强度增强而显著增大的状态。由于电磁铁存在磁饱和现象，即通入电磁铁线圈的电流一直增大到某一值时，电磁铁的磁感应强度不会继续增大，此时电磁力的大小也不会发生改变。
因此需要判断在对机器人的控制过程中，是否会出现电磁铁的磁饱和现象。
[0150]
线圈电阻的大小为37.75ω，利用arduino单片机作为核心控制单元口，其端口的输出电压大小为0到5v。
[0151]
用端口最大输出电压除以电阻可以得到电流的最大值为：
[0152]
磁饱和电流计算公式如下式(4-12)所示：
[0153][0154]
式中：i
sat
：磁饱和电流值
[0155]
n：线圈匝数3050匝
[0156]
ae：铁心横截面积大小为707mm
2
[0157]
b：饱和磁通量密度大小为841mt
[0158]
l：线圈电感值
[0159]
因此要计算出饱和电流的大小，还需要计算出线圈的电感值，本次毕业设计制作的电磁铁是螺旋线圈，其电感为：
[0160][0161]
式中：l：螺旋线圈的电感值
[0162]
l：螺旋线圈的长度约为400m
[0163]
n：线圈匝数为3050匝
[0164]
s：线圈横截面积大小为0.196mm
2
[0165]
μ：线圈内部磁芯导磁率大小为300h/m
[0166]
k：长冈系数值为0.85
[0167]
把上面的值代入式(4-13)得到线圈电感值为：
[0168][0169]
将电感的值代入式(4-12)求出饱和电流的大小为：
[0170][0171]
可以看出计算得到的让电磁铁出现磁饱和现象时的电流的值，远远大于控制电磁力时所需要通入电磁铁的最大电流值，因此在整个电磁力的变化范围内，电磁铁都不会出现磁饱和现象。
[0172]
根据求得的电磁力的变化公式，以电流为横坐标，电磁力大小为纵坐标，电磁力的最大值为60n左右。
[0173]
假设斜拉索的倾斜角α＝60
°
，机器人两个轮臂之间的夹角θ＝120
°
，而由式(2-8)得到将机器人稳定在斜拉索上最小值：
[0174][0175]
当机器人的质量较大时，只靠单片机输出端口的电压可能不能提供足够的吸附力让机器人稳定在斜面上。因此，为了解决电磁力不足的问题，可以在单片机的输出端口串联上一个dc-dc升压电路，让电压的最大值增大，从而让电磁力增大。
[0176]
将使用arduino-uno单片机作为核心控制单元实现电流的自动控制。arduino是一款方便使用、操作便捷的开源电子平台。主要包含两个主要的部分:硬件部分是arduino电路板，可以用它来实现电路连接，与其他元器件相连；另外一个是arduino ide，它为你的计算机中的程序提供了开发环境。只需要在ide中撰写相应程序的代码，将程序输入到arduino电路板后，程序就会告诉arduino开发板需要完成什么任务。该单片机输出端口可以提供0-5v的电压。将机器人的振动作为干扰，使用震动传感器检测机器人的抖动状态。
[0177]
该传感器使用压电陶瓷产生振动，当压电片振动时，产生电信号。当振动程度不同时，输出的振动等级数值不同，振动等级共分为10级，而单片机输出端口电压值最大为5v，因此可以得出振动等级可电压值得关系为式(5-1)：
[0178]
v＝振动等级
÷
10
×
5＝0.5
×
振动等级
ꢀꢀꢀꢀꢀꢀꢀꢀ
(5-1)
[0179]
利用振动和电压的关系即可实现振动越大电流越大，吸附力越大，保证机器人的稳定性。
[0180]
磁吸附机器人在斜拉索上运动，若出现振动，则振动传感器会检测到振动信号，并且用arduino单片机读取振动信号值，再根据信号大小，调整输出端口电压值的大小，振动信号值越大，输出端电压值越大。最后利用升压模块将单片机输出端口的电压值放大，通入电磁铁两端。即可实现自动控制吸附力大小。
[0181]
位置式pid控制输出的即为当前的控制量，无需换算，可以直接用来控制电流，因此本次任务使用位置式pid控制，其基本公式为式(5-2)：
[0182][0183]
由式(5-2)我们很清楚地知道相应的控制器参数包括比例系数k
p
，积分控制系数k
i
，微分控制系数k
d
。这三个参数的取值直接影响到了pid控制系统整体的好坏。
[0184]
利用pid控制来控制单片机输出端口的电压从而控制通入电磁铁线圈的电流。因此式(5-2)中的input是当前端口的电压值，setpoint为目标电压值，即用振动等级计算后的电压值。偏差值e(t)即为目标电压值与当前电压的差值。输出u(t)为使用pid控制后的电压值。
[0185]
当机器人在斜拉索上面运动时，若机器人发生大幅度振动，即可根据振动传感器检测到的振动等级，利用公式(5-1)转换为相应等级的电压值。再使用pid控制，将输出端口电压由当前电压调整至振动等级相应的电压值以供使用。
[0186]
在pid控制的过程中，输出端口的电压值不断变化，因此电磁铁线圈内通入的电流值也会不断变化，电磁吸力也随之变化，慢慢使机器人稳定下来。当电磁吸力到达了某一值机器人停止抖动时，输出端的端口电压停止发生变化，此时电磁力保持稳定。振动传感器依旧在检测是否产生振动，若机器人机身又发生了振动，则重复上述控制流程。
[0187]
当机器人在斜拉索上运动时，若机身产生抖动，此时有下面两种控制方法使抖动消除。
[0188]
(1)直接增大电磁力，不管机器人抖动的规律如何，直接将电磁力增大到某一固定值，从而使机器人稳定下来。
[0189]
(2)跟随机器人抖动规律来改变电磁力，当机器人发生抖动时，将机器人的抖动规律模拟成一条曲线，此时调整电磁力跟随抖动的曲线变化而变化，使机器人稳定。
[0190]
两种方法都可以使机器人稳定，但是第一种方法对电能的消耗大，而且由于电磁力过大，会使机器人的运动编的不灵活。针对第二种方法，我将机器人抖动分为正弦变化和无规则变化两种变化函数，并且对其进行仿真。
[0191]
1.直接增大电磁力的pid控制仿真
[0192]
直接增大电磁力就是将电压由某一值瞬间跳转到下一个值，此时电压的变化类似一个阶跃变化模型。单位阶跃函数模型如下式(5-3)所示：
[0193][0194]
当t<0时，函数值为0；当t>0时，函数值为1；在t＝0这段时间内，函数值由0变为1。下面对机器人在斜拉索上面的抖动做一个模拟仿真。
[0195]
假设机器人一开始在斜拉索上运动时单片机输出端口的输出电压值为2.5v，此时因为受到外界影响机器人机身发生抖动，此时直接调整单片机输出端口电压为3.5v。那么电压的变化函数是由2.5v到3.5v的阶跃变化函数，因此电压值变化函数如式(5-4)所示：
[0196][0197]
对电压变化函数进行拉氏变换得：
[0198][0199]
由式(5-2)我们知道，pid控制公式为：
[0200][0201]
因此得到最基本的pid控制的基本传递函数为式(5-6)所示：
[0202][0203]
我们随机设置pid控制参数k
p
＝k
i
＝k
d
＝0.3，将电压变化的阶跃函数作为输入，pid控制函数作为传递函数。利用matlab/simulink进行数学仿真，得到pid控制后的电压变化的仿真曲线如图4所示。
[0204]
由图4可以看出，电压变化曲线的超调量，调节时间和稳态误差都非常大，因此当设置pid参数为k
p
＝k
i
＝k
d
＝0.3时，控制系统的动态性能非常差，不能满足设计的要求。因此我们根据前面分析的结论，适当增大比例控制系数k
p
和积分控制系数k
i
,减小微分控制系数k
d
的值。因此再重新设置一组新的pid参数：k
p
＝1、k
i
＝1.2、k
d
＝0.005。此时pid仿真曲线如图5所示。
[0205]
可以看出，当调整了pid参数后，系统的调节时间减小了，且超调量和稳态误差几乎为0。整个控制系统的动态性能可以达到要求。因此当机器人机身发生抖动时，若直接增大电磁力到某一固定值，可以将机器人稳定住。
[0206]
2.跟随抖动规律的电磁力pid控制仿真
[0207]
当机器人在斜拉索上发生正弦抖动时，振动传感器将会检测到一条呈正弦变化规律的振动等级变化曲线。因此跟随抖动变化，单片机端口的输出电压也会随之产生正弦变化如式(5-7)所示：
[0208][0209]
该式表示的是，在t
0
时刻前机器人保持电压为a对应的电磁力匀速运动，在t
0
时刻机器人发生抖动，且抖动的幅度为b。
[0210]
假设机器人一开始在斜拉索上运动时单片机输出端口的输出电压值为2.5v，此时因为受到正弦变化风力影响机器人机身发生正弦抖动，且振动传感器检测到的振动等级变化范围是6级到7级，则由式(5-7)可得，此时的目标电压值在2.5v到3.5v间变化。电压变化如图6所示：
[0211]
因此电压值变化函数如式(5-8)所示：
[0212][0213]
对函数进行拉氏变换得：
[0214][0215]
利用之前的pid参数k
p
＝1、k
i
＝1.2、k
d
＝0.005，利用matlab/simulink进行数学仿真，得到pid控制后的电压变化的仿真曲线如图7示。
[0216]
由图6可以看出利用之前的一组pid参数进行仿真，得到的结果误差非常大，我们设计一组新的pid参数此时k
p
＝0.05、k
i
＝0.2、k
d
＝0.01。用这一组pid参数得到的电压仿真曲线如图8所示。
[0217]
3.针对不规则变化抖动的pid控制仿真
[0218]
当机器人机身在绳索上发生不规则抖动时，振动传感器接收到的数据将是一组无规则的曲线，因此在matlab-simulink设置一个随机函数的信号源，将采样时间设置为0.3s，采样范围为1-5v，生成的电压变化曲线如图9所示。
[0219]
将生成的电压变化随机函数作为信号源，使用pid控制，其参数随机设置为k
p
＝0.01、k
i
＝0.09、k
d
＝0.1，利用matlab/simulink进行数学仿真，得到pid控制后的电压变化的仿真曲线如图10所示。
[0220]
由图10可以看出，加入pid控制后，电压的变化曲线基本能与机器人抖动曲线变化趋势相同，能在一定程度上控制机身抖动，但是控制的准确性有待提高，因此面对不同的变化函数，需要使用不同的pid控制参数。
[0221]
由前面的仿真实验结果来看，当机器人遇到不同形式的抖动，电压变化为不同的函数时，使用同一组pid参数进行pid控制的效果差别很大。因此电磁力pid控制的参数选择对磁吸力的控制尤为重要，这时就需要引入pid的参数整定了。因此在后续的研究中可以设计出pid参数整定的方法，让机器人能根据检测到的振动信号，自动选择效果最好的pid参数进行控制。当pid控制参数合适时，电磁力的控制效果可以满足使用。
[0222]
以上所述仅是本发明的优选实施方式，应当指出：对于本技术领域的普通技术人员来说，在不脱离本发明原理的前提下，还可以做出若干改进和润饰，这些改进和润饰也应视为本发明的保护范围。