1 前言
多軸聯(lián)動的線性插補(bǔ)及其加減速算法是高檔數(shù)控系統(tǒng)的核心技術(shù)。加減速處理分前加減速處理和后加減速處理。前加減速處理在插補(bǔ)之前，其優(yōu)點(diǎn)在于對合成速度進(jìn)行控制，不影響位置精度，但要進(jìn)行減速點(diǎn)的預(yù)測。前加減速通常采用直線加減速。后加減速是對各插補(bǔ)軸分別進(jìn)行加減速控制，由于各軸之間沒有協(xié)調(diào)關(guān)系，因此，合成位置可能不準(zhǔn)確。
 
S曲線加減速通常用于后加減速處理，我們在開發(fā)基于開放式數(shù)控系統(tǒng)的多坐標(biāo)聯(lián)動纖維纏繞機(jī)時，將其應(yīng)用于前加減速處理，取得了很好的效果。
2 多軸聯(lián)動的線性插補(bǔ)S加減速
S曲線加減速規(guī)劃是指在加減速時，使其加速度的導(dǎo)數(shù)(Jerk) da/dt為常數(shù)，通過對加Jerk值的控制來最大限度地減小對機(jī)械系統(tǒng)造成的沖擊。另外，通過對加速度和Jerk兩個物理量的參數(shù)設(shè)定或編程設(shè)定，可實(shí)現(xiàn)柔性加減速控制，以適應(yīng)不同種類機(jī)床的工況。
設(shè)在n維線性插補(bǔ)數(shù)控程序中，任意一段插補(bǔ)數(shù)據(jù)為
| P1, P2, …，Pn, F |
其中：F為合成速度，P1～Pn為各插補(bǔ)軸當(dāng)前段的位移。
根據(jù)線性插補(bǔ)原理，各插補(bǔ)軸的位移與速度比相等，則有應(yīng)于各插補(bǔ)軸的分速度。令 |Pi| = P =TSEP      i=1, 2, …, n
 
Fi F
 (1)

圖1 “S加減速”規(guī)劃原理圖

 
式中：P=(Pi2)½表示合成位移；TSEP表示該線性插補(bǔ)段各軸同時到達(dá)終點(diǎn)所需時間；F1～F2為合成速度對應(yīng)于各插補(bǔ)軸的分速度。
令 Ki= Pi       i=1, 2, …, n
 
P
 (2)

則 Fi=KiF      i=1, 2, …, n (3)

在前加減速處理時，對給定速度進(jìn)行規(guī)劃，如圖1所示，整個加減速過程分為三段，即加速段(1, 2, 3區(qū))，勻速段(4區(qū))和速段(4, 5, 6區(qū))。在加速段和減速段，又分別包括變加減速區(qū)(1, 3, 5, 7區(qū))和恒加減速區(qū)(2, 6區(qū))：變加減速區(qū)，|da/dt|=J, Jerk為恒值；恒加減速區(qū)，|a|=A，加速度為恒值，勻速段(4區(qū))的速度為恒值Vc。
各軸的運(yùn)動參量成比例
對已規(guī)劃的合成速度v(t)在te點(diǎn)冪級數(shù)展開，令∆t=t-tx，有 v(t)=v(tx)+a(tx)∆t+½J(te)∆t2 (4)
同理，各插補(bǔ)軸對應(yīng)．點(diǎn)速度 vi(t)=vi(te)+ai(te)∆t+½Ji(tx)∆t2F      i=1, 2, …, n (5)

根據(jù)線性插補(bǔ)原理，合成速度與各插補(bǔ)軸速度有下列比例關(guān)系： vi(t)=Kiv(t)F      i=1, 2, …, n (6)

對于上面恒等式，應(yīng)有 vi(tx)=Kiv(tx), ai(tx)=Kia(tx), Ji(tx)=KiJ(tx) (7)

由于tx為任意一點(diǎn)，此式表明段內(nèi)加減速過程中各插補(bǔ)軸的速度、加速度和Jerk分別與合成的速度、加速度和Jerk對應(yīng)成比例。當(dāng)對合成速度按S曲線規(guī)劃時，各插補(bǔ)軸在保證空間軌跡的同時，也按S曲線進(jìn)行加減速，即S曲線加減速可用于前加減速控制。同時，上述關(guān)系可用于各插補(bǔ)軸的速度、加速度和Jerk的極限值檢查。
S加減速的插補(bǔ)遞推公式
設(shè)插補(bǔ)周期為T，則在第k個插補(bǔ)周期結(jié)束時的合成位移Sk為 Sk = ∫ tk v(t)dt= ∫ tk-1 v(t)dt+ ∫ tk-1+T v(t)dt=Sk-1 ∫ t (vk-1+ak-1t+½Jt2)dt =Sk-1+vk-1T+½ak-1T2+(1/6)JT3
0 0 tk-1 0
 
 (8)

第k個插補(bǔ)周期內(nèi)的合成位移增量為 ∆Sk=vk-1T+(1/2)ak-1T2+(1/6)JT3=vk-1T+(1/2)(ak-1+(1/3)JT)T2=vk-1+(1/2)akT2=(vk-1+(1/2)akT)=vkT (9)
ak=ak-1+(1/3)JT (10)
vk=vk-1+(1/2)akT (11)
注意，上述遞推公式是分區(qū)適應(yīng)的，即 J= { J, T∈[t0，t1]∪t6，t7
0, t∈(t1，t2)∪(t3，t4)∪(t5，t6)
-J, t∈[t2，t3]∪[t4，t5]
 (12)

只要初始條件ak-1和vk-1給定，則可推導(dǎo)出各插補(bǔ)周期的合成位移增量。進(jìn)而得到各插補(bǔ)軸在插補(bǔ)周期內(nèi)的位移增量，其公式為 ∆Pik= Pi ∆Sk=Ki∆Sk
 
P
 (13)
區(qū)間的判別
段內(nèi)加減速時，每程序段伺服電動機(jī)速度總要減到零后再執(zhí)行下一程序段。因此其加速段和減速段的位移相等，見圖1。
1區(qū)(t0-t1)的初始速度和初始加速度為0，則在t1時刻的位移Pti=(1/6)Jt13，其加速度a1=A=Jts，速度Vt1=(1/2)At12=(1/2)Ats，則 ts=t1=A/J (14)

由圖1中的加速度圖線可以看出 V=(1/2)Ats+Atl+(1/2)Ats=A(ts+tl) (15)

則 tl=(V/A)-(A/J) (16)
ta=2ts+tl=(V/A)+(A/J) (17)
由ts、tl、tm可計(jì)算出t0～t7的時間，并以此進(jìn)行區(qū)間判別。
終點(diǎn)判別
終點(diǎn)距離∆S及各軸距終點(diǎn)距離∆Si

圖2 插補(bǔ)計(jì)算流程框圖

 
∆S=P- k ∆Sk
∑
1
 (18)
∆Si=Pi- k ∆Pik
∑
1
 (19)

3 算法描述
加減速的規(guī)劃在插補(bǔ)預(yù)處理中進(jìn)行。
讀入A、J、V、Pi。
計(jì)算ki，計(jì)算出的各軸的速度和加速度將與各軸的極限值比較，如果超出極限值，修正V。
掃描速度倍率K，令V=KV。
初始化，t0=O，∆Si=Pi。
速度規(guī)劃，計(jì)算tl、ts、tm。
置速度規(guī)劃完成標(biāo)志。
插補(bǔ)。
插補(bǔ)的計(jì)算流程框圖如圖2所示，插補(bǔ)預(yù)處理的加減速規(guī)劃中，已得出各區(qū)的時間范圍，插補(bǔ)時只需判斷當(dāng)前插補(bǔ)周期所在區(qū)間，即可按插補(bǔ)迭代公式計(jì)算出各插補(bǔ)軸與速度規(guī)劃適應(yīng)的位移增量，從而實(shí)現(xiàn)其加減速。
4 結(jié)語
本文給出了基于速度S曲線加減速規(guī)劃的插補(bǔ)迭代算法，并應(yīng)用于前加減速處理。雖然只考慮了一般情況下的段內(nèi)加減速，但其思想和方法對更復(fù)雜的段內(nèi)加減速和段間加減依然適用。本文給出的算法可采用加速度和加速度導(dǎo)數(shù)編程，能最大限度地適應(yīng)各類數(shù)控機(jī)床的加工工況，同時，本算法也非常簡潔和易于實(shí)現(xiàn)。