1
M U
“Trong bn tóm tt lun án này, các hình v và công thc đã đc đánh s theo đúng
nh trong lun án”.
H truyn đng là mt tp hp các c cu ghép ni c khí phc v bin đi tc
đ, moment. Lun án quan tâm ti lp h truyn đng qua bánh rng đ truyn
moment quay và thay đi vn tc góc quay.
Nhim v ca bài toán điu khin h truyn đng qua bánh rng là phi xác đnh
đc quy lut thay đi moment dn đng to ra t đng c
dn đng đ h có đc
tc đ góc quay ca ti đu ra luôn bám n đnh đc theo mt qu đo đt và cht
lng đó không b nh hng bi khe h gia các bánh rng, ma sát, moment ti cng
nh đ không cng vng ca vt liu làm bánh rng.
Các phng pháp điu khin hin có
ã có nhiu công trình nghiên cu v k thut đ gii quyt các vn đ nêu trên.
Chúng đc chia thành hai loi là bin pháp c khí và đin.
Bin pháp c khí hin dùng là lp thêm bánh đà, nâng cao đ chính xác khi ch
to các chi tit, điu chnh và lp ráp theo các quy trình nghiêm ngt, chp hành các
ch đ bo qun bo dng và bôi trn Các gii pháp c khí ch thích hp vi ch
đ làm vic xác lp c
a h thng cng nh h thng có tính đng hc bin đi chm.
Bin pháp v đin, điu khin có th tóm tt nh sau:
iu khin vi mô hình xp x tuyn tính bng b điu khin PI: ây là phng pháp
ph thông nht và trc đây cng đc s dng nhiu nht. Nó đc s dng khi h
truyn đng là mô t xp x tuyn tính đc di dng tuyn tính tham s hng. B
điu khin đc s dng là b điu khin PI có hàm truyn:
1
() 1
p
I
Rs k
Ts
⎛⎞
=+
⎜⎟
⎝⎠
(0.1)
Ngi ta cng có th b sung thêm tính chnh đnh thích nghi cho tham s b
điu khin nhm nâng cao hn na tính bn vng ca h thng điu khin. Song theo
nhiu tài liu tham kho thì vic ch s dng PI t chnh thích nghi không đ đ làm
gim dao đng xon trên trc mt cách hiu qu mà phi s dng thêm các b điu
khin phn h
i trng thái.
iu khin h truyn đng có khe h: Xác đnh khe h và điu khin loi b s nh
hng ca khe h ti cht lng truyn đng là bài toán thng gp nht trong các bài
toán điu khin h truyn đng.
2
Khe h có mô hình toán nh sau:
khi 0 và ( )
( , , ) khi 0 và ( )
0 ngoài ra
mu u m u a
Buu u mua
τ
ττ τ
>=−
⎧
⎪
== <=+
⎨
⎪
⎩
(0.4)
iu khin bù khe h đc hiu là xác đnh mô hình ngc ca (0.4).
1
(, , )uB uuu
−
=
(0.5)
iu khin thích nghi bù khe h bng mng neural và h m: Vn đ tn ti là ta
không có mô hình (0.4) tuyt đi chính xác cho khe h. Nh vy ch cn mt sai lch
nh trong mô hình (0.4) s dn ti mt sai s rt ln trong phép nghch đo (0.5). Bi
vy có th nói k thut điu khin bng hàm ngc là không kh thi trong thc t.
Hình 0.6: iu khin bù khe h bng mng neural
Trên c s suy lun nh vy, ch trong vài nm gn đây nhiu công trình đã đc
công b cho vic thay hàm ngc (0.5) bng vic xp x nó nh h m hay mng
neural nh mô t hình 0.6.
Mc dù vy nhng phng pháp điu khin bù xp x này cng có mt hn ch
ca nó, nh:
− Vic xp x hàm phi tuyn nh mng neural hay h m ch có th có đc kt
qu xp x vi sai lch nh tùy ý trong min gii hn cho phép, nu nh hàm phi
tuyn cn đc xp x đó là liên tc.
− Vic bù bng mng neural hình 0.6 ch thc s có ý ngha khi tín hiu
τ
bên
trong h truyn đng có khe h là đo đc.
iu khin h truyn đng có khe h, ma sát và đ đàn hi: Theo mt vài nghiên cu
gn đây thì phn ln h truyn đng có khe h luôn tách đc thành hai khâu phi
H truyn đng có khe h
u
τ
y
H truyn
đng lý
tng
B điu
khin PI
,
pI
kT
w
e
1
s
3
tuyn mc ni tip gm khâu mô t khe h đng trc và mt khâu phi tuyn dng
affine truyn ngc cht đng phía sau:
1
khi 1 1
() (,) ()
kk
n
xx kn
xf dtg
τ
+
=≤≤−
⎧
⎨
=+ +
⎩
xx x
(0.7)
vi (,)dtx là hàm bt đnh. Ngoài ra, nhng tài liu này còn khng đnh vic nâng cao
cht lng bù khe h nh c cu chnh đnh thích nghi PI có th thay đc bng b
điu khin phn hi trng thái gán đim cc. Bi vy khi s dng mô hình trng thái
(0.7) ta có đc cu trúc điu khin bù khe h cho h truyn đng bng phn hi trng
thái đc mô t
hình 0.9.
Hình 0.9: Bù khe h moment ma sát và moment xon bng phn hi trng thái
Chính t cu trúc điu khin bù khe h bng b điu khin phn hi trng thái
thay vì b điu khin phn hi đu ra PI thích nghi đó mà ngi ta đã hoàn toàn d
dàng b sung vào cu trúc điu khin bù khe h thêm mt khâu phn hi trng thái th
hai có nhim v nhn dng đ bù các thành phn hàm bt đnh (,)dtx này, đc xem
nh hàm mô t moment ma sát ()
ms
Mt và đàn hi, đ điu khin h truyn đng va
có khe h ma sát và đ đàn hi ca vt liu (hình 0.9).
iu khin h truyn đng không theo nguyên lý bù: Mt ý tng rt khác v điu
khin h truyn đng khp ni mm có khe h, so vi nhng phng pháp nêu trên, đã
đc gii thiu là phng pháp điu khin d báo (MPC).
phng pháp MPC này, các thành phn bt đnh nh khe h, ma sát, moment
xon không đc nhn dng. Thay vào đó đu ra ca h luôn đc so sánh vi đu ra
ca mt mô hình xp x tuy
n tính có các thành phn bt đnh đc lý tng hóa bng
các giá tr c th. Tín hiu điu khin s đc xác đnh trên c s cc tiu hóa hàm sai
lch gia hai tín hiu đu ra này.
H truyn đng có khe h
u
τ
H truyn
đng lý
tng
B điu khin
phn hi
trng thái
d
x
e
1
s
x
Nhn dng
moment ma sát
và moment xon
4
Mc dù phng pháp này vn cn đn mô hình xp x tuyn tính, song cng nên
đc tham kho và phát trin tip sau này vì tính đc đáo ca nó.
Tính cp thit ca đ tài
Hình 0.13 mô t cu trúc ca mt h truyn đng qua bánh rng. Lun án s
nghiên cu s dng phng pháp điu khin mà đó tt c các thành phn không th
xác đnh đc chính xác s đc xem nh nhng thành phn bt đnh ca mô hình
toán, thay vì nhn dng và quan sát chúng. Vi nhim v đt ra đó, lun án hng ti
phng pháp điu khin t
ng quát, có th khng ch, loi tr các thành phn bt đnh
cng nh sai lch mô hình toán bng phng pháp điu khin phn hi ph mà không
cn nhn dng cng nh xác đnh xp x khe h, ma sát cng nh moment xon trên
trc, đ không cng vng tuyt đi ca vt liu làm bánh rng đ điu khin bù các
nh hng ca thành ph
n bt đnh đó.
Hình 0.13: Bài toán điu khin h truyn đng qua bánh rng
Mc tiêu nghiên cu
Trên c s nhim v nh vy, lun án đã đt ra mc tiêu:
− Xây dng mô hình tính toán đng lc hc đi vi mt h truyn đng c khí ca
máy t hp nói chung, trong đó có tính đn yu t đàn hi (moment xon), ma
sát tnh, ma sát đng và khe h gia các bánh rng di dng các bt đnh hng
s và bt đnh hàm s.
− Xây d
ng phng pháp điu khin thích hp trên nguyên tc kt hp các
phng pháp điu khin hin có nh thích nghi, bn vng, logic m và mng
neural, cho h truyn đng có c ba yu t bt đnh nêu trên.
2ms
M
1ms
M
4ms
M
3ms
M
M
d
ϕ
1
M
3∼
Bin tn
B điu khin
ϕ
2
M
c
Ti
5
CHNG 1: LÝ THUYT IU KHIN THÍCH NGHI
BN VNG H PHI TUYN
1.1 Các khái nim mô t
1.1.1 Khái nim n đnh Lyapunov
Xét h phi tuyn có mô hình thay đi theo thi gian, còn gi là h không dng:
(,)t=
xfx
(1.1)
Khi đó h s đc gi là:
a) n đnh, ∀ 0
ε
> ,
0
0t > ,
0
(, ) 0
δ
ε
∃>t sao cho:
00
(, )t
δε
<x ⇒
00
(, , )tt
ε
<
xx vi mi
0
≥tt (1.4)
trong đó
00
() (, , )ttt
=
xxx là nghim ca h phng trình vi phân (1.1) tha mãn
điu kin đu
00000
() (, ,)ttt==xxx x.
b) n đnh tim cn, nu nó va n đnh, va tha mãn
00
lim ( , , )
t
tt
→∞
= 0xx .
1.1.2 Tiêu chun xét tính n đnh Lyapunov
nh lý LaSalle: Xét h phi tuyn (1.1) cân bng ti gc. Ký hiu:
()
12
() , ()Vt
γγ
≤≤xx x vi
12
,
γ
γ
∞
∈
K (1.8)
(,) (,) ()
VV
Vt t W
t
∂∂
=+ ≤−
∂∂
xfxx
x
(1.9)
Khi đó:
− H s là n đnh, nu
()W x
là hàm bán xác đnh dng, tc là
() 0, W ≥∀xx
.
− H s là n đnh tim cn, nu ()W x là hàm xác đnh dng trong mt lân
cn O ca gc, tc là () 0, W ≥∀∈Oxx và () 0W
=
x khi và ch khi
=
0x .
1.2 Các phng pháp điu khin phi tuyn
1.2.1 iu khin n đnh thích nghi và nguyên tc certainty equivalence
áp dng đc phng pháp Lyapunov gián tip cho vic tng hp b điu
khin h bt đnh
(, ,,)t=
xfxu
θ
(1.22)
ngi ta đã làm nh sau:
6
1. Trc tiên ta gi s là đã có các tham s hng bt đnh
θ
. Khi đó vi phng pháp
Lyapunov gián tip ta tìm hàm CLF (,,)Vtx
θ
và b điu khin GAS tng ng
(, ,,)t=urxw
θ
.
2. Tip theo, ta thay thành phn bt đnh
θ
trong (,,)Vtx
θ
và (, ,,)t=urxw
θ
bi p ,
mà trong nhiu tài liu thng ký hiu là
θ
, đ có
(, ,)Vtxp
,
(, , ,)trxwp
, ri tìm
cách hiu chnh
p
, tc là tìm quy lut thay đi cho
(, , , ,)t
=
pcxwup
, đ vn có
đc tính xác đnh âm ca (, ,)Vt
xp .
1.2.2 iu khin thích nghi theo mô hình mu
có th ch đng to ra đc thêm các cht lng n đnh khác tt hn, ngi
ta thng áp dng phng pháp điu khin thích nghi theo mô hình mu, nu nh h
phi tuyn bt đnh không nhng có dng truyn ngc mà còn là truyn ngc cht:
2
1
1
1
(, , ) ()
n
n
T
n
n
x
x
x
x
x
fx x u
−
⎛⎞
⎛⎞
⎜⎟
⎜⎟
⎜⎟
⎜⎟
=
⎜⎟
⎜⎟
⎜⎟
⎜⎟
⎜⎟
+
+
⎝⎠
⎝⎠
… hx
θ
(1.25)
Hình 1.3 biu din cu trúc điu khin thích nghi theo mô hình mu này, trong đó:
− B điu khin là:
()
()
01 1
( ) , , , , , trong đó
TT
T
nn
ufv aa xx
−
=− − + = =……ax x a x (1.26)
vi các h s
01
, ,
n
aa
−
… đc chn sao cho mô hình mu:
mm
w=+A
xxb,
012 1
010 0
0
001 0
,
0
000 1
1
n
aaa a
−
⎛⎞
⎛⎞
⎜⎟
⎜⎟
⎜⎟
⎜⎟
⎜⎟
==
⎜⎟
⎜⎟
⎜⎟
⎜⎟
⎝⎠
⎜⎟
−−− −
⎝⎠
A
b
(1.27)
là n đnh vi cht lng đt trc.
− B bù tín hiu
vwz=− có:
1
()
()
T
T
z
−
⎧
=
⎪
⎨
=
⎪
⎩
EP
phxbe
phx
(1.28)
trong đó P là nghim đi xng xác đnh dng ca phng trình Lyapunov:
T
+=−AP PA Q
và E, Q là hai ma trn đi xng xác đnh dng tùy chn, s làm h (1.25) bám
theo đc mô hình mu (1.27), theo ngha sai lch mô hình
m
=−exx gia
chúng luôn b chn và tin tim cn v 0.
7
1.2.3 iu khin trt
Nguyên lý điu khin trt liên quan đn h có cu trúc truyn ngc:
1
khi 1 1
() (,)
kk
n
xx kn
xf dtu
+
=≤≤−
⎧
⎨
=+ +
⎩
xx
vi
(
)
1
, ,
T
n
xx= …x (1.35)
trong đó (,)dtx là hàm bt đnh tha mãn điu kin b chn:
,
sup ( , )
t
dt
δ
≤
x
x (1.36)
và ()f x là hàm h thng đc gi thit là đã bit trc. Nhim v điu khin là phi
to ra đc tín hiu ra
1
yx= ca h bám n đnh tim cn theo tín hiu mu ()wt
mong mun cho trc, tc là phi làm cho h có:
lim ( )
t
t
→∞
=
0e và ()t
<
∞e vi
(
)
(1
,, ,
T
n
ee e
−
=
…e , ewy
=
− (1.37)
Bn cht ca b điu khin trt là không trc tip to ra đc cht lng mong
mun (1.37) cho h thng, mà gián tip qua mt trt, đc hiu là mt:
(
)
(2) (1)
12 1
( ) , 1
T
nn
n
saeae ae e
−−
−
=++ + + =
eae (1.38)
có các h s
()
12 1
,, ,
T
n
aa a
−
= …a đc chn sao cho đa thc đc tính ca nó:
21
12 1
()
nn
n
ps a as a s s
−
−
−
=+ + + +
là Hurwitz. Vi mt mt trt (1.38) nh vy, thì đ có đc cht lng mong mun
(1.37), ta ch cn thit k b điu khin sao cho có đc:
lim ( ) 0
t
s
→∞
=e và
()st
<
∞
Mi b điu khin tha mãn:
()
()
( ) khi ( ) 0
(, )
( ) khi ( ) 0
nT T
nT T
wfks
uw
wfks
⎧
>+−−+ >
⎪
=
⎨
<− − − − − <
⎪
⎩
aw ax x e
x
aw ax x e
(1.39)
vi
()
2
, ,
T
n
xx= …x
, tc là
(
)
1
,
T
T
x=xx và k
δ
> , s làm h tha mãn (1.37).
Hình 1.3: Minh ha nguyên tc thích
nghi theo mô hình mu
u x
v
i tng
(1.25)
B điu khin
(1.26)
a)
m
x
z
w
u x
v
e
Mô hình mu
(1.27)
C cu bù
(1.28)
i tng
(1.25)
B điu khin
(1.26)
b)
8
1.3 iu khin thích nghi vi h m
Xét h kín có s đ mô t hình 1.7, trong đó đi tng điu khin đc gi
thit là tuyn tính và mô t đc bng hàm truyn
()Ss
. B điu khin m có quan h
truyn đt
()ufe=
cng đc gi thit là trong ch đ xác lp xp x đc bi hng s
khuch đi
FC
k . B điu khin m này có thêm h s khuch đi
p
k
đu ra là thay
đi đc. Nu ch xét riêng ch đ xác lp, sai lch
()et
s là:
1
1(0)
FC p
ew
kkS
=
+
⇒
()
2
(0) (0)
1(0)
1(0)
=− =−
+
+
FC FC
pFCp
FC p
kS kS
de
we
dk k k S
kkS
Nu nh h còn có
(0) 1
FC p
kkS >> thì ta có th xp x tip đc:
pp
de e
dk k
=−
S dng hàm CLF
2
1
()
2
Ve e= s đc:
2
pp
pp
de e
Veee k k
dk k
== =−
Bi vy, đ có 0V
<
, ta chn c cu chnh đnh:
2
p
p
e
k
k
γ
=
vi
0
γ
>
tùy chn. (1.43)
1.4 Kt lun
Chng 2 đã trình bày mt s phát trin chính ca lý thuyt điu khin các h
phi tuyn trong các nm gn đây. Lý thuyt điu khin các h phi tuyn phát trin t
khái nim mô t đn các phng pháp thit k có h thng. Các khái nim thit k điu
khin các h phi tuyn da trên c s các hàm Lyapunov đn các khái nim n đnh
đu vào-đu ra, tr
ng thái-đu vào đã đc trình bày. Phng pháp thit k các b điu
khin bn vng, thích nghi và ti u gián tip đc trình bày ngn gn.
w u
y
e
i tng
điu khin
B điu
khin m
Hình 1.7: iu khin h truyn
đng bng b điu
khin m thích nghi
p
k
Chnh đnh
thích nghi
9
CHNG 2: XÂY DNG MÔ HÌNH TOÁN CHO H
TRUYN NG QUA BÁNH RNG
2.1 Mô hình tng quát
Sau đây ta thc hin xây dng mô hình thc nghim h truyn bánh rng có tính
đn yu t đàn hi và hiu ng khe h đ tin hành nghiên cu cht lng ca b điu
khin. Hn na, s là không mt tính tng quát nu nh đây ta ch xây dng mô hình
toán cho h có mt cp bánh rng.
Hình 2.1: H nhiu cp bánh rng là h truyn ngc ca nhiu h mt cp bánh rng
2.1.1 Cu trúc vt lý và các đnh lut cân bng
Hình 2.2a) mô t hình thc ghép ni ca cp bánh rng, đc đánh s bánh rng
1, bánh rng 2 và hình 2.2b) biu din li cu trúc vt lý ca nó. Vi gi thit vt liu
làm các trc bánh rng là có đ cng tuyt đi, còn vt liu làm các bánh rng có b
bin dng trong quá trình làm vic, các rng ca hai bánh rng đang n khp vi nhau
ti đim n khp
P
, vt liu có đ cng tuyt đi thì t s truyn ca chúng đc vit:
02
11 2 2
12
22 1 01 1
L
L
r
rz
i
rrz
ϕ
ω
ϕω
= = =− =− =− (2.1)
trong đó
112 2
,
ω
ϕω ϕ
==
là vn tc góc tng ng ca hai bánh rng,
12
,
LL
rr là bán
kính ln tng ng ca hai bánh rng (bán kính ngoài),
01 02
, rr là bán kính và
12
, zz là
s rng ca hai bánh rng.
a) b)
y
x
M
d
M
c
M
ms
2
M
ms
1
1
2
J
1
J
2
DC
J
d
Hình 2.2: Cu trúc vt lý ca h truyn đng qua
mt cp bánh rng
1
2
M
d
M
d
ϕ
2
ϕ
1
M
2
ϕ
3
M
3
ϕ
4
M
c
Ti
10
2.1.2 Mô hình tng quát có tính đn hiu ng khe h, đ đàn hi ca vt liu
và moment ma sát
Trên c s h thng truyn đng hình 2.2, ta đã có đc mô hình đng lc hc
có tính ti yu t đàn hi ca cp bánh rng và ma sát trong các trc nh mô t trên
hình 2.3.
Dùng mt ct n-n, trên đó chu mt moment đàn hi ca hai bánh rng nh trên
hình 2.3. Gi J
1
là moment quán tính ca phn bên trái bao gm moment quán tính ca
rotor đng c dn đng, moment quán tính ca trc và bánh rng 1 và bên đó chu tác
đng ca moment dn đng ca đng c đin là M
d
, lc ma sát trong là M
ms
1
. Do có
bôi trn nên lc ma sát t l vi vn tc góc ca trc dn. Còn phn bên phi chu tác
đng ca mt moment đàn hi có chiu ngc li cng nh moment ma sát. Gi J
2
là
moment quán tính ca phn bên phi ca bánh rng b dn 2.
Moment đàn hi trên hai bánh rng là:
101011022
()Mcrrd rd
ϕ
ϕ
=+,
202022011
()
ϕ
ϕ
=
+Mcrr r
Suy ra:
22
11 1 1 12 2 1
22
22 2 2 211 2
cos ( )
cos ( )
LL dms
LL cms
Jcr i MM
Jcr i MM
ϕαϕϕ
ϕαϕϕ
⎧
++=−
⎪
⎨
−+=−−
⎪
⎩
(2.8)
trong đó:
−
11d
JJJ=+
và
12
,
LL
rr bán kính vòng tròn ln ca hai bánh rng.
−
L
α
góc n khp ca hai bánh rng đánh giá khe h gia các bánh rng. Trong
trng hp hai bánh rng tiêu chun thì
0
20
L
αα
== .
− khi có
0c =
h s đang ch đ khe h. Khi cc
=
thì h đang ch đ n khp
bánh rng và
c
là đi lng đánh giá đ cng ca bánh rng. Giá tr
c
càng nh,
đ mm do ca bánh rng càng ln.
−
12
,
ms ms
MM là thành phn moment ma sát trong các đ trc.
c
J
1
J
2
M
d
M
c
M
ms
2
M
ms
1
M
dh
M
dh
c
Hình 2.3: S đ đng lc hc
1
x
2
y
J
d
0
n
n
n
n
11
2.2 Mô t h ch đ xác lp
Xét riêng cho trng hp h có ch đ xác lp (chy đu), tc là:
111ms
Mb
ϕ
=
và
212ms
Mb
ϕ
=
(2.10)
Lúc đó, cùng vi gi thit này, mô hình (2.8) tr thành:
22
11 1 L 1 12 2 11
22
22 2 L 2 211 22
cos ( )
cos ( )
ϕαϕϕ ϕ
ϕ
αϕ ϕ ϕ
⎧
++=−
⎪
⎨
−+=−−
⎪
⎩
Ld
Lc
Jcr i Mb
Jcr i Mb
(2.11)
Ngoài ra, vì có thêm (2.1) nên, khi ký hiu tip:
22 2 2
1L12L2
cos , cos
αα
==
LzL z
cr c cr c
ri thay vào phng trình (2.11) ta có:
11 11 1 1 12 2
22 22 2 2 211
()
()
zd
zc
Jbc i M
Jbc i M
ϕϕ ϕ ϕ
ϕϕ ϕ ϕ
⎧
+++=
⎪
⎨
+− + =−
⎪
⎩
(2.12)
Mô hình cui cùng (2.12) trên chính là dng tng đng ca (2.8) nu có đc
thêm gi thit (2.10) v moment ma sát.
2.3 Kt lun
Trong chng này ta đã đa ra đc mô hình (2.8) cho h truyn đng qua mt
cp bánh rng và t mô hình đó, di các gi thit b sung thêm, ta còn có đc mô
hình (2.12) đn gin hn đ mô t ch đ chy đu ca h.
12
CHNG 3: IU KHIN THÍCH NGHI VÀ BN VNG
H TRUYN NG QUA BÁNH RNG
Chng này trình bày kt qu ca lun án xung quanh bài toán điu khin h
truyn đng có mô hình toán tng quát (2.8) đc xây dng t chng 2 nh sau:
22
11 1 1 12 2 1
22
22 2 2 211 2
cos ( )
cos ( )
LL dms
LL cms
Jcr i MM
Jcr i MM
ϕαϕϕ
ϕαϕϕ
⎧
++=−
⎪
⎨
−+=−−
⎪
⎩
(3.1)
3.1 iu khin m bng phng pháp bù sai lch mô hình
3.1.1 Xây dng c cu bù theo nguyên tc cân bng vi mô hình mu
Gi thit đ áp dng phng pháp bù này, hay còn gi là thích nghi theo mô hình
mu, là h truyn đng gn nh đã đi vào ch đ n đnh vi mt tc đ bng hng s.
Vi gi thit này, hai moment ma sát
12
,
ms ms
MM khi đó s đc xp x bi:
111 2 22
,
ms ms
MbM b
ϕ
ϕ
==
(3.2)
vi
12
, bb là hai hng s không cn phi bit trc. Vi gi thit này, mô hình (3.1)
ban đu tr thành:
1234
()ay ay ay ay u dt+++=+
(3.4)
vi các tham s hng và các ký hiu tín hiu vào, ra, tp nhiu (disturbance) nh sau:
12
1
22
221
(cos)
LL
JJ
a
cr i
α
=
,
12 1 2
2
22
221
(cos)
LL
Jb bJ
a
cr i
α
+
=
,
22 2
12 1 2 1 2
3
22 2
221221
(cos)
(cos)
LLL
LL L
bb J cr r J
a
cr i r i
α
α
−
=+
222
12 1 2
4
222
221 2 21
(cos)
(cos)
LLL
LLL
rb bcr
a
ri cr i
α
α
=−
,
d
uM
=
,
2
y
ϕ
=
và
2
1
2
221
()
Lc
L
rM
dt
ri
=−
T mô hình tng đng (3.4) có đu vào u , đu ra y , cng nh khi b qua
đc các moment ma sát gia tc ca h truyn đng qua bánh rng, ta nhn thy ngay
rng h có hàm truyn dng quán tính bc 3:
32
1234
1
()Gs
as a s as a
=
+++
(3.5)
là mt h tuyn tính tham s hng bt đnh vi thành phn tp nhiu
()dt
đu vào
ph thuc moment ti.
Nu ta ký hiu h kín cn đt đc phi là h có cht lng mong mun, đc
th hin qua hàm truyn cùng cu trúc vi mô hình đi tng:
13
32
123
1
()
1
m
Gs
as a s as
=
+++
(3.6)
vi
123
, , aa a
đc chn trc phù hp vi cht lng đt ra, thì do phn tp nhiu
()dt có ln trong đi tng (3.4) ch tác đng đu vào nên ta có th bù s nh hng
ca nó mt cách rt đn gin theo nguyên tc điu khin theo mô hình mu (3.6) nh
biu din hình 1.3. Nh vy, b điu khin phn hi đu ra lúc này ch còn phi đa
đc sai lch đu ra
m
eyy=− ca h tim cn v 0.
3.1.2 Thit k b điu khin m thích nghi cho h đã đc bù bt đnh
Hình 3.5: M hóa tín hiu vào ra và lut hp thành ca b PI m
Vi lp đi tng quán tính bc 3 (3.5) và mô hình mu (3.6) ta có phng trình
mô t tín hiu sai lch
e
gia hai tín hiu ra hình 3.1 và sai lch mô hình
()dt
:
1234
()ae ae ae ae dt+++=
vi
2
ey
ϕ
=
−
(3.7)
Bi vy, nu b điu khin trong hình 3.1 là b điu khin PI:
/
/
0
1
t
pp
I
ukke edt
T
⎛⎞
=+
⎜⎟
⎝⎠
∫
(3.8)
ta s có lut chnh đnh thích nghi (1.43) cho
p
k . Khi m ca b điu khin trong
hình 3.1 đ to hàm phi tuyn xác đnh
//
,
pI
kT có tín hiu vào
1
E ,
2
E và ra
U
đc m
hóa bng các hàm thuc
1
i
mf
E
,
2
i
mf
E
,
i
mf
U
, 1, , 7i
=
… dng hình chuông hình 3.5.
1
i
mf
E
2
i
mf
E
i
mf
U
ji
1 2 3 4 5 6 7
1 1 1 1 1 2 3 4
2 1 1 1 2 3 4 5
3 1 1 2 3 4 5 6
4 1 2 3 4 5 6 7
5 2 3 4 5 6 7 7
6 3 4 5 6 7 7 7
7 4 5 6 7 7 7 7
14
B điu khin m có lut hp thành:
R
ij
: Nu
1
1 i
mf=
E
E và
2
2 j
mf=
E
E thì
k
mf=
U
U
, ,1, ,7ij
=
… (3.11)
vi
4 khi 8 3
4 3sgn( 8) khi 8 3
ij ij
k
ij ij
⎧
+− +− ≤
⎪
=
⎨
++− +−>
⎪
⎩
3.2 iu khin thích nghi bn vng trong không gian trng thái
3.2.1 Xây dng b điu khin bám thích nghi bn vng trên nn điu khin
trt và gi đnh rõ cho h phi tuyn truyn ngc cht
H truyn ngc cht, bt đnh có mô hình:
1
khi 1 1
(,) (,) (,)
kk
TT
nf g
xx kn
xtdttu
+
=≤≤−
⎧
⎪
⎨
=++
⎪
⎩
fx x gx
(3.12)
trong đó:
− (,), (,)tt
f
xgx là hai vector hàm rõ (đã bit),
− ,
fg
θ
θ
là hai vector tham s hng s bt đnh và
−
(,)dtx
là thành phn tp nhiu (hàm) bt đnh tha mãn điu kin b chn:
,
(,) sup (,)
t
dt dt
δ
∞
==<∞
x
xx (3.13)
nh lý 3.1: Xét h truyn ngc cht, bt đnh, mt đu vào (3.12). Khi đó, vi mi
tín hiu mu
()wt
kh vi
n
ln cho trc, b điu khin phn hi trng thái:
1
() ()
1
( , ) sgn( )
,
(,)
n
knT
kf
k
T
g
ae w t s
u
t
λ
λ
δ
−
=
+− +
=>
∑
fx
gx
θ
θ
(3.14)
trong đó
()k
w là ký hiu ca đo hàm bc k ca ()wt , cùng c cu chnh đnh:
1
1
() ( ,)
() ( ,)
f
g
se t
se tu
−
−
⎧
=−
⎪
⎨
=−
⎪
⎩
F
G
fx
gx
θ
θ
(3.15)
s làm cho qu đo trng thái
()tx ca h bám n đnh đc theo qu đo mu:
()
(1)
,, ,
n
d
ww w
−
=
…x (3.16)
trong đó , FG là hai ma trn đi xng xác đnh dng tùy chn,
(2) (1)
12 1
()
nnT
n
se ae ae a e e
−−
−
=++ + + =
… ae
(3.18)
là mt trt có
()
(1)
, , ,
T
n
ee e
−
=
…
e ,
1
ewx
=
− và các h s
()
11
, , ,1
T
n
aa
−
= …a
là nhng s dng tùy chn đ đa thc:
15
21
12 1
()
nn
n
paa a
α
ααα
−
−
−
=+ + + + (3.19)
là Hurwitz (có tt c các nghim nm bên trái trc o).
3.2.2 ng dng cho h truyn đng qua bánh rng
Khi bin đi mô hình tng quát (3.1) v dng:
1
4
khi 1 3
(,)
kk
T
fg
xx k
xdtu
θ
+
=≤≤
⎧
⎪
⎨
=+ +
⎪
⎩
xx
θ
(3.37)
Ta có b điu khin (3.14) ca nó:
3
1()(4)
1
sgn( )
θλ
−
=
⎛⎞
=+−+
⎜⎟
⎝⎠
∑
kT
gk f
k
uaew sx
θ
(3.38)
vi
123
, , aa a đc chn bt k, min là vi chúng đa thc:
23
12 3
()
α
ααα
=+ + +paaa
(3.39)
là Hurwitz và
λ
phi đc chn tha mãn
λ
δ
> , trong đó
δ
là giá tr chn trên ca
hàm bt đnh
(,)dtx
tính theo (3.13).
Mt trt (3.18) ca b điu khin s là:
123
()=+++
se ae ae ae e
trong đó
1
=−ewx và ()wt là tín hiu đt kh vi ít nht bn ln mà đu ra
12
ϕ
=
x ca
h phi bám theo.
C cu chnh đnh thích nghi (3.15) s là:
1
()
() , 0
θξξ
−
⎧
=−
⎪
⎨
=− >
⎪
⎩
f
g
se
se u
F x
θ
(3.40)
vi ma trn
33
×
∈RF đi xng xác đnh dng tùy chn và 0
ξ
> .
3.2.3 Kt qu mô phng trên MatLab
Tin hành mô phng h truyn đng bánh rng (3.37) vi thành phn bt đnh
hàm s
(,)dtx
là nhiu n trng có chun vô cùng
0.5d
∞
= và tín hiu mu
( ) sin(0.1 )wt t= , thì khi chn các tham s sau cho b điu khin:
3
30=FI, 0.1
ξ
= ,
123
125, 75, 15aaa
=
==, 0.2
λ
=
b điu khin bn vng s làm cho h kín bám theo đc tín hiu mu ()wt . Hình 3.17
là đ th biu din tín hiu đu ra thc
12
()xt
ϕ
=
so sánh vi tín hiu mu ()wt .
16
0 20 40 60 80 100
-150
-100
-50
0
50
100
150
Time (s )
theta mu f
Hình 3.15: Kt qu chnh đnh các tham s
f
θ
0 10 20 30 40 50 60 70 80 90 100
-1
-0.5
0
0.5
1
1.5
Time (s )
tin hieu dau ra
tin hieu dat
Hình 3.17: Kt qu tín hiu đu ra thc
1 2
x
ϕ
=
ca h và tín hiu mu w
Hình 3.15 biu din giá tr các tham s ,
fg
θ
θ
ca b điu khin đc chnh đnh
thích nghi, tc là đu ra ca b chnh đnh (3.40). Ta có th thy khi vào ch đ xác
lp, các giá tr tham s này cng s tin ti mt hng s c đnh. Tuy nhiên các hng
s đó không bt buc phi là hng s
,
fg
θ
θ
thc ca đi tng điu khin.
3.3 iu khin thích nghi bn vng vi phn hi tc đ
3.3.1 Mô hình phn hi tc đ
Do ch quan tâm ti tc đ
2
ϕ
nên t mô hình đi tng, cng nh khi s dng
li các ký hiu v hng s bt đnh ,
fg
θ
θ
và hàm s bt đnh (,)dtx , ta s có vi ký
hiu trng thái mi thay th:
12
22
32
x
x
x
ϕ
ϕ
ϕ
⎛⎞
⎛⎞
⎜⎟
⎜⎟
==
⎜⎟
⎜⎟
⎝⎠
⎝⎠
x
,
d
uM=
[2]
f
θ
[3]
f
θ
[1]
f
θ
17
mô hình trng thái bt đnh tng đng vi (3.37), nhng là bc 3:
1
3
, 1,2
(,)
θ
+
==
⎧
⎪
⎨
=+ +
⎪
⎩
kk
T
fg
xx k
xdtu
θ
xx
(3.42)
Mô hình trng thái bc 3 này là hoàn toàn tng đng vi mô hình toán tng
quát (3.1) bc 4 ca h truyn đng bánh rng, nu nh bài toán điu khin thích
nghi bn vng cho nó ta ch quan tâm ti tín hiu mu ()wt cho trc là tc đ ca
góc quay
2
()t
ϕ
.
3.3.2 Thit k b điu khin thích nghi bn vng
Tng ng vi mô hình bc 3 (3.42) ca đi tng điu khin, ta cng có ngay
đc b điu khin phn hi trng thái thích nghi bn vng tng đng, đc xây
dng theo đnh lý 3.1 nh sau:
1
()
() , 0
θξξ
−
⎧
=−
⎪
⎨
=− >
⎪
⎩
F
f
g
se
se u
θ
x
và
(
)
1
12
sgn( )
θλ
−
=++−+
T
gf
uaeaew s
θ
x (3.43)
3.3.3 Kt qu mô phng trên MatLab
0 10 20 30 40 50 60 70 80 90 100
-25
-20
-15
-10
-5
0
5
10
15
20
25
Time (s )
tin hieu mau
tin hieu dau ra
Hình 3.28: Kt qu mô phng tín hiu mu ()wt và tín hiu đu ra
12
x
ϕ
=
0 10 20 30 40 50 60 70 80 90 100
-0.4
-0.2
0
0.2
0.4
0.6
0.8
1
1.2
1.4
Time (s)
Hình 3.31: So sánh tham s bt đnh [1]( )
f
t
θ
vi tham s chnh đnh
[1]( )
f
t
θ
B điu khin phn hi trng thái đng (3.43) có các tham s
1
1.8a = ,
2
1.08a
=
.
và:
18
100
020
003
⎛⎞
⎜⎟
=
⎜⎟
⎜⎟
⎝⎠
F
,
0.1
ξ
=
, 0.9
λ
=
Kt qu mô phng v cht lng bám n đnh gia tín hiu đt ()wt tùy ý kh vi
3 ln và tín hiu đu ra
12
x
ϕ
=
đc biu din trong hình 3.58.
Ngoài ra, trên nguyên tc ca đnh lý 3.1 thì b điu khin thích nghi bn vng
(3.43) cho mô hình trng thái bc 3 ca h truyn đng bánh rng (3.42) ch có th
đm bo đc cht lng bám n đnh bn vng cho tín hiu ra
12
x
ϕ
=
theo tín hiu
mu ()wt nu nh các tham s bt đnh ,
fg
θ
θ
ca mô hình là hng s.
Tuy nhiên ta vn có th s dng b điu khin này ngay c khi h truyn đng có
các thành phn bt đnh bin đi chm theo thi gian
(), ()
fg
tt
θ
θ
. Hình 3.28 và 3.31
biu din kt qu mô phng xác nhn khng đnh trên.
3.4 Kt lun
Chng này trình bày hai phng pháp gii quyt cho cùng mt bài toán đt ra
ban đu ca lun án là điu khin bám n đnh h truyn đng qua bánh rng, có đ ý
đn hiu ng khe h, moment ma sát, moment xon và moment ti mà không cn phi
nhn dng các thành phn bt đnh đó.
C hai phng pháp trên cng đã đc thc hin kim chng bng mô phng
trên MatLab cho mô hình h truyn đng qua bánh r
ng (3.1) vi các kt qu mô
phng trình bày các hình 3.28 - 3.31.
19
CHNG 4: XÂY DNG MÔ HÌNH VT LÝ H
TRUYN NG QUA BÁNH RNG VÀ
CÁC KT QU THC NGHIM
4.1 Xây dng mô hình thc nghim
Cu trúc mô hình bàn thí nghim h truyn đng đc biu din hình 4.1 vi:
− Máy tính Pentum IV, phn mm Matlab và phn mm ControlDesk Version 5.0.
− Card điu khin DS1104, Driver Servo motor Midi-Maestro 140x14/28 và đng
c, khp ni hai bánh rng và ti.
Hình 4.1 Cu trúc h thng thc nghim
4.2 Kt qu thí nghim vi b điu khin PI
4.2.1 Kt qu thí nghim vi b điu khin PID kinh đin
Không ti Thay đi ti
Hình 4.15: Tc đ
2
ϕ
khi có tín hiu tc đ đt
()
() 200sinwt t
π
=
Hình 4.15 biu din kt qu thí nghim vi hai b PID. B PID th nht có tác
dng n đnh dòng cho c cu chp hành là đng c dn đng vi
100, 7
p
pI
I
k
kk
T
===
Tc đ đt
Tc đ đt
20
và b điu khin PID th hai có tham s
1200
p
k
=
và 3
I
k
=
. Tín hiu đt là hàm điu
hòa:
()
() sin 2wt A ft
π
= (4.1)
4.2.3 Kt qu vi b điu khin PI m
Không ti Thay đi ti
Hình 4.18: Tc đ
2
ϕ
khi có tín hiu tc đ đt
()
() 50sin 2wt t
π
=
Hình 4.18 biu din kt qu thí nghim ng vi tín hiu tc đ đt cng có dng
hình sin cho bi (4.1), tc là khi tc đ đt thay đi giá tr liên tc. Thí nghim đc
tin hành vi b điu khin m có hai tín hiu vào và mt tín hiu ra. C ba bin ngôn
ng ca b điu khin m đc m hóa bng 7 giá tr m (tp m) cho hình 3.5.
4.3 Kt qu thí nghim khi có thêm khâu chnh đnh thích nghi theo
mô hình mu
Nhm nâng cao cht lng cho h thng, ta s đa thêm vào h hình 4.4 mt
khâu khuch đi
p
k theo s đ cu trúc đã trình bày hình 1.7 cùng vi lut chnh
đnh thích nghi (4.1) cho nó đ đu ra ca h bám theo đc tín hiu ra ca mô hình
mu. S đ khi mô t h thng thí nghim m thích nghi này trên mô hình vt lý
đc biu din hình 4.8.
Hình 4.21 là kt qu thc nghim thu đc vi khâu chnh đnh thích nghi (1.43)
và mô hình mu có hàm truyn
3
1
()
(1 20 )
m
Gs
s
=
+
. Kt qu thí nghim cho thy vic
đa thêm khâu chnh đnh thích nghi theo mô hình mu cng phn nào đã ci thin
đc cht lng điu khin, song không nhiu. Nói cách khác dao đng trong h vn
tn ti và không th loi b đc mt cách trit đ, mc dù nghiên cu sinh đã tin
hành th nghim vi rt nhiu các b tham s khác nhau.
Tc đ đt
Tc đ đt
21
Không ti Thay đi ti
Hình 4.21: Tc đ
2
ϕ
khi có tín hiu tc đ đt
()
() 50sin 2wt t
π
=
4.4 Kt lun
T thc nghim ta thy rng khi cha có b điu khin, h truyn đng qua bánh
rng dao đng, đ n rt ln. S dng b điu khin PI hay PI m cho h thng truyn
đng có s tham gia ca bánh rng đã cho phép ta gim đáng k nhng dao đng gây
nên bi khe h, đàn hi và ma sát ca bánh rng. Khi có tác đng điu khin h thng
chy êm, ting n c khí gim đi rt rõ rt. Kt qu trên đã khng đnh tính đúng đn
ca thut toán và cho phép áp dng vào điu khin các h thng truyn đng trong
thc t.
Tc đ đt
Tc đ đt
22
Kt lun và nhng hng nghiên cu tip theo
Lun án đã thc hin đc các công vic sau:
− Xây dng mô hình toán cho h truyn đng mt cp bánh rng.
− Xây dng đc b điu khin m thích nghi theo mô hình mu phn hi đu ra
cho h khi làm vic ch đ chy đu.
− xut
đnh lý 3.1 và trên nn đnh lý 3.1, xây dng đc b điu khin thích
nghi bn vng phn hi trng thái.
Lun án cng còn đã xác nhn bng mô phng rng b điu khin thích nghi bn
vng phn hi trng thái đc thit k theo ni dung đnh lý 3.1 vn áp dng đc khi
các vector tham s ,
fg
θ
θ
thay đi chm theo thi gian, mc dù khi thit k ta cn phi
có gi thit chúng là hng s.
Tuy nhiên, đ hoàn thin hn na cht lng điu khin cho h truyn đng, mt
s hng m rng sau nên đc nghiên cu tip tc:
1. Th nht, do b điu khin thích nghi bn vng ca đnh lý 3.1 đc xây dng trên
nn phng pháp điu khin trt nên không th tránh khi hin tng rung trong
h.
0 10 20 30 40 50 60 70 80 90 100
-150
-100
-50
0
50
100
150
Time (s )
Tin hieu dau ra u
Hình 5.1: Hin tng rung trong h bám thích nghi bn vng
nâng cao cht lng cho h thng, cn thit phi làm gim hin tng rung này.
Trc đây, khi đa ra đnh lý 3.1 làm nn tng cho vic thit k b điu khin thích
nghi bn vng, ta có đ cp ti kh nng làm gim hin tng rung nh b sung
thêm khâu xp x hàm phi tuyn bt đnh
(,)dtx
bng
(,)dt
x
.
2. Th hai, trong trng hp s dng mô hình bc 3 (3.42) ca h truyn đng và gi
s rng ta có th xp x đc:
23
12
23
3
(,) (,,)
T
fg
xx
xx
xdtudutuxx x
θ
⎧
=
⎪
⎪
=
⎨
⎪
=+ += +
⎪
⎩
θ
(5.1)
Khi đó, ch vi mt khâu bù đn gin:
(,,)uvd utx=−
ta s có ngay đc:
1
3
, 1, 2
kk
xx k
xv
+
==
⎧
⎨
=
⎩
(5.2)
Bi vy vn đ nghiên cu đây là có th hay không biu din đc mô hình bt
đnh (3.42) ca h truyn đng v dng (5.1). Nu có th thì cn phi b sung thêm
gi thit gì?
3. Th ba là kh nng s dng b điu khin thích nghi m. Hin ti lun án mi ch
s dng mt b điu khin m Mamdani c đnh, đc b sung thêm khâu chnh
đnh thích nghi bên ngoài (gi là b điu khin m thích nghi).
Tuy nhiên phng pháp m thích nghi này ch yu vn ch là thay đi khâu thích
nghi bên ngoài b điu khin m, ch cha cho thy đc kh nng thích nghi ca
bn thân b điu khin m, tc là cha áp dng đc kh nng t chnh đnh thích
nghi các giá tr ngôn ng, lut hp thành hay gii m trong b điu khin m.
iu này đã không cho thy ht đc tính u vit ca h thích nghi m. Bi vy
mt hng m tip theo là nghiên cu s dng thích nghi m, thay cho m thích
nghi. Cu trúc b thích nghi m phù hp trong trng hp này s là b điu khin
m Takagi-Sugeno, thay cho m Mamdani.
Vi thích nghi m Takagi-Sugeno ta s
có c hi không cn s dng thêm bt c
mt khâu chnh đnh thích nghi h s khuch đi nào bên ngoài b điu khin m.
24
Các công trình đã công b
[1] Hà,L.T.T. và Lãi,L.K: Hai gii pháp nâng cao cht lng h truyn đng có khe h. Tp
chí Khoa hc & Công ngh HTN, S 4-2009, trang 34-37, 2009.
[2] Hà,L.T.T.; Lãi,L.K. và Nguyt,L.T.M: Kho sát cht lng ca h truyn đng có khe h.
Tp chí Khoa hc & Công ngh HTN, S 3-2009, trang 124-130, 2009.
[3] Lãi,L.K. và Hà,L.T.T: Mt phng pháp nâng cao cht lng h truyn đng qua bánh
rng. Tuyn tp hi ngh toàn quc ln th 5 v c đ
in t, trang 134-137, 10.2010.
[4] Lãi,L.K. và Hà,L.T.T: Nghiên cu thc nghim điu khin m áp dng cho h truyn đng
qua bánh rng. Tuyn tp báo cáo Hi ngh toàn quc v iu khin và T đng hóa,
VCCA-2011, trang 759-763, 11. 2011.
[5] Li Khc Lãi, Lê Th Thu Hà, Lê Th Minh Nguyt và Nông Lê Huy: Mt phng pháp
điu khin thích nghi h truyn đng qua bánh rng. Tp chí Khoa hc & Công ngh i
hc Thái nguyên, T
p 88 - S 12-2011, trang 163-167, 2011.
[6] Ha,L.T.T. and Phuoc,N.D.: A Design of an Adaptive SM Tracking Controller for Two Wheel
Gearing Transmission Systems. Submitted and accepted for ISTS-2012, ThaiLand, 2012.
[7] Lê Th Thu Hà và Nguyn Doãn Phc: Thit k b điu khin bám thích nghi bn vng
cho h phi tuyn bt đnh và ng dng vào điu khin h truyn đng qua bánh rng. Tp
chí Khoa hc & Công ngh à nng, S 10(59) -2012, trang 1-6, 2012.
[8] Hà,L.T.T. và Phc,N.D: iu khin bám thích nghi h phi tuyn bt đnh có đ ý ti t
p
nhiu và ng dng vào điu khin h truyn đng qua bánh rng. Tuyn tp báo cáo Hi
ngh C hc toàn quc ln th 9, 2012.
[9] Ha,L.T.T. and Phuoc,N.D.: Robust and Adaptive Tracking Control of Two Wheel Gearing
Transmission Systems. Proceeding of 6
th
Vietnam National Conference on Mechatronic
VCM-2012.
[10] Phuoc,N.D. and Ha,L.T.T.: Robust and Adaptive Tracking Controller Design for Gearing
Transmission Systems by Using its Reduced Order Model. Journal of Science and
Technology. Technical Universities, No. 90, 2013.
[11] Phuoc,N.D. and Ha,L.T.T.: Model Reference Adaptive Controller Design for Gearing
Transmission System. Journal of Science and Technology. Technical Universities, No. 91,
2013.