壓染機輥輪系統之動態建模與分析

(1)

,It Wd~fl~ ~Ntzi/] ~ ~;f~~>t;flT Dynamic modeling and analysis of a

padder roller system

$~ '" ~* Ft. ~#

Abstract

This paper is intended to deal with the dynamic modeling and analysis of a padder roller system. First, the method of dynamics is utilized for the development of the kinetic and potential energy ofthe system and the governing equations of motion are derived from Hamilton's principle. Afterwards, the open loop transfer function is obtained from the actuator and sensor's locations, and the IMSL . subroutine is used to find all the poles and zeros of the corresponding open loop system. The dynamic characteristics such as natural frequencies and mode shapes are analyzed by the frequency equation acquired by assume-mode method. Through the transfer function and dynamic characteristics developed by this research, it is expected to be used as the key for the design of controllers.

Keywords: Dynamic modeling, Padder roller, Assume-made-method, Natural frequencies, critical speeds

'"

xj~*

*~~.#~~

••

^-*~~

••

^~.~~

~o~~~ffl.h.~~*.~

••

^~~~.~~

{lltit' .it]:)'i.*~JtJ$.l!(Hamiiton'sprinciplc)1.f.t'J ~

~iJ.,ti fIJ • ~;f.l i\. ' .<.&' tl:dft t-;f!L $if

**

^4MLapiace

transform) ;j(.lli ~fo/t.1lt9 ~JliI.(Open-loop)~;f~ JUt' ;f'l ffl IMSL ~!!Ukft ;t<.~{U.s $4 ~,t5 0 .it J,:).{FHt:m • ii;-(Assume-mode-method).;R1.fflb_

*'

^iJt^{t.J!i~^f^~

.i\.'~#~.~• •~~~~~.f'

••

^#

~·*~~*M~~.~• •*~{t.J.8~.A•

• ~tl°~~• • •~~~B{t.Ja~A.b~#~

4>'t 1# IIItf{t.J

*"

^~1t<At •

~-~: ••^~.·M*••• ·.a••^~·~

t~f

00- *jl~1tJll~~~~-rtoo

•••

A-.~.~L.t'ffl~.M~.~

.~~fr*~.~~L.m'~~.A~~~~~

~*~.'~~A*~*~~~.m·~~~.~

~••t.~• •{t.J&~.,.~.a~*•••

~~.~A.~~~.{t.J••~~t,~~*}.

~•• ,~*~~~~*~.~~,.#~.~W

*k••{t.J_It~.'~~~~*~~m~~~

··~A~~*t~~.~~fr*~£

••

^l!'M

~.~ffl~

.•

~~fr~*~.~.l!·

~• •~~~i\..~~L.m~.t~~a

*~{t.J~~'*~~•••~a~~.~~.*~

a*{t.J~h'.~*••••^~~~~~m~• •

.~t~,~~~~~~t'~~• • •±~m~

~~'~~~m~.±~~~~.&I~~~~' AI~R%~4~._~~~••M·~••*~

M.~t-;r;if'-{iT.MBt' ~1t51~ttAItJ";\l!lt' ~~

.~*.MA._~

.•

^'~~~~

•••

^·~~

~~~~•••~~%.k~~4&&I~*~

~*·~M~-.~~*,~{t.J~~'*.~~.M

;fITM;~(analytical solution)[1,2] • ;fH!l7tt;k(finite element method,FEM)[3,4J • ~.;fHEF"';~(transfer

matrix method)[5,6] , {'!:f~M.'Z. M@'~1tlttiaJl ⁰ A

T • •~~• •*~~a~~~B~~~••~

••~%.'*-~.*~*,~~.3a••*~

lbjp;k(Root locus method)~·fT• •~;f1T ' JUJ,1KH9:.

;f}t

.)t;;

(Assume-mode-method) .;Ri~M;,. ~ .~/t.1tJ ~ f~.i\.'~#~.~• •^*,~~~t••^4.~

.It'••

^4t,ti⁰_{.=. •}_;f~_j·Uit_~

2.1 f,'t• •~*~;f}t1'!~~

*~~M.#~~~.A-~I~.~~.

•••~*,~{t.J~~~••• ' • •^A~.~.~

if'-~'M~{t.J*h6~~W~if'-~~'A• •^*~

~ ~h 't~%-Jco 00 -=-I'IT -r :

00.::. r.t._+{:h-r;'too

AT~~£

•••••

^+~~~•••1'!' • • • _~~M~{t.J*h~W.if'-~~'.~••'~A -fi~~._~~~.*'~'~~h-Jco~OO~M

-r:

z III -=. *~.

**

⁺

*'

^{ft.Jt. ,}^:;;t^~^j]⁰⁰

2.2 f~_.+

*

^.~

••

;Q';f1.i\.4it l,f

-Jco1ll~M-rA-._.+*~~~.~~~

~1I5.z.;F t mJ :

x

155

(2)

·-*~~-~.tl~'~~~~AOO~~.·OO t~xnA-oo~~.*~'~4• •~~~~

fh!4~~' ~Jt• •i£1£-;t-z ^1.ll.'lL-f~~1ii:l:.+

-f,cg , ~00 li.1'IT:iT- :

1Illi. - f~1Jti + *-g • •!4 $,.1: il n-dJt

t • •~Q~~~~.M'ftAi£X-Y~~_~

U-V{ji]~*}Bi-'$,.C&hti~7i-Jl'JA~ --e 'xyzA

;f!l~-t~7 ~i~*,Yt • ,,~ G A~J,I; , P J.Hl#.!-fi

,j)'iH, *-g J:..~1£-;t - J,I; , x' y' z ,dJ7i-Jjll +1t~ X ' Y - Z ~ , i • j , k A x • y • z ~ ~ ~.f1.ll ~Y , M

P.L~~g~~m.h.~M;f!l.~• •~m .t/ti,ltfj:

r. = (U + ~z-Q8z-Qy)i + ci; + Qz-Qcpz - 8z)j +(8y+Qcpy-~x 1-Q8x)k (1)

'l!t Jil:. - t"li +^*-^{g • •}^iliJ^~^,^$,..~E.(kinetic energy)

A:

.:lis:

1 if f' ,

Tk = - rprppdAdZ 20A

=~^S(tjl⁺V')PAdZ+~ JI,[(a2~), ^+(a²^;^)2}jZ

20 20 az az

+~tfI Q2[(a2u)2+(a'v)2Jdz+~lf ^J^Q²^dZ

20 ' az' az' 2 0 '

11 evau auav

+-fJ O ( - - - ) d Z (3)

2 0 ' az az az az

",,-h""'''108-\'..b Ilfl avau auav)lZ","~<i*

'l-'O"'J~bJ!b'l[JO\'T' - Q ( - - - ( -<.-'><"""

20' azaz DZaZ

# ~licr .it It(Coriolis acceleration)!'1T _ Ii. {ji] # ~

~' • •#'A-~~tl~'{f.L~t~A~~~

~'m'ti.9£(diametral inertia) , 1, A.m.~:tMHi.9£

(polar inertia) , ~Jjll A :

I, = fix' +z')pdA= fey' +z')pdA pn(r; _r;4)

I

^A ^A ⁴ ⁽⁴⁾

lJ,

= lex' +y')pdA pncr:2^-r;')

~t~·~7i-~.~~*~~+~·iliJ._.~

(bending) ~ ~1'1T '!II; ~5;p.,ti.1ft liE. (elasticity potential energy)A :

V. =~'fEI[tU)l"-(a'v)~Z ⁽⁵⁾

rL 20 az' az'

r

il.:fllHamilton'sprinciplc[8], $,.~a~;t~A :

t2 t2

f8(Tk - VPE )dt + foWnedt = 0 (6)

tl tl

$,. t oWnc A 4H,f-';j' h 1t ~,Yt!'1T1-t {ji]:<!J ' Wf J\(3)&

~(5)1\A~(6):i!$-, or1~r _7'J JiB :;7a~1$,f/t.. k

,~Jt.-**~t/) Ii! ~:U.~ :

EIU'''' +o'I,D' + pAD -I;D" -J,QY" = 0 (7) {

EIV' + Q'I,V" + pAY - I, Y" + J ,QU" = Q(t)

1Ht-m,E}j :

• a .. 8² "a' '. a'

( ) = Ot ' ( ) = ~ , ( ) = az' ' C) = 8Z' U=U(z,t) : X ~ ~J ~1H,P'J.t(m)

V=V(z,t) : Y ~~~1ft~~J.t(m)

A : m~~lIlIidii~( m' ) E : ~~~)tit(

%3 )

I : ,tJ'rHE( m' )

£ : .mt.ii.ptlt(m)

p : fIt.~~\fg(

k%l )

Q : fl.~~it(rad/s)

Is : t't.~t1ii1§}'f.J'IJHi.9£(Kg-m) J, : fH.,~iit1i.t.·rHE(Kg-m) 2,3 f/t... -f- * ,Yt#ilt 1A!4#1lt J!1it

fIt..M-f-*~~CAOO~ffl~.'$,.*~~.

.rr-{lKt1H^ur:

U(O,I)=U(£,I)=O 'EIU"(O,t)=EIU"(R,t)=O (9) V(O,t)=V(£,t) 0 'EIV"(O,t)cEIV"(f,t)=O (0)

~4-.{ji]• • •-f-~A*lt~*~IO~R {.llt~( Euler's beam)' JL-Jj~i£&~u&tlJ&J.t t ' 1':

~~~.#~~{ji]~~'m.A~Mit~-**,~~

.i£mt,i;t X ~~(7j(+fuJ)~1ft~1~+ ⁰ "4~JJa~.1t

(Coriolis acceleration)Ptf _A {ji] # ~~{ji] ~ ~-t!!. ~

~~~~~~,~~,a~MJ\mMkA:

EIV"" + Q2I, V" + pAY - I,y' = Q(t) (11 ) 4'~aK~m1t1ft~ V ~ :

V(z, t)~ <j> (z)v(t) (12)

A T -t-~Ml1!i"I±M(homogeneous solution) , 4'#J#

~lX, Q(t)=Oo ~f~(12)1\A~(1l)&' jU{f-~it8;t

M$t :!l Fa' ~.. • "Til

r

⁷¹^J^~*1. ~ :

Elf (z) + Q 2I,cp"(z) = vet) = a ' 00' (13) pAcjl(z)- I,cp"(z) - V(I)

Wi((13)<if.jj4l :

Elf' (z) + I,(Q' + (2)f(z) -pA<j>Cz) =0 (14)

v(t)+co'v(t)=O (15)

~~.lJ\(l4)~i!MA :

<P (z)=c l cos~!z~c2sini3lz+c)coshj3,z+C4sinh~2z (16)

$,.4'

-

^{I,(Q2 +(}²^)+JC1,(Q2^{+(0 2}

^»'

^+4ElpA

I l - '

^{13, -} JCI,(Q2 + 00'»2 1- 4EIpA- I.(Q2 ^2EI + CD 2) ⁽¹⁷⁾

13, - 2EI

Jt t C" C2,C),C ⁴A~~ .:rr..{:£-;t

'*"

^it^~it~~1! .rr-fl,tt 14- ⁰ tinlJi-1',tt14-J\(l 0) , Of il :

cp(O) =<l>U) =Elcp"(O) = Elcp(R) = 0 (18)

~fcp(0)=<P"ce)=01\AJ\(l6)or4l :

(19) iHf cp(e) = cp"(R) = 0 1\AJ\(l6)~1l :

c,Sin i3 IP +C4 sinhi32e= 0 (20)

{-C2i3~ sini3,l +c.i3; sinhi32e = 0

~A~,1l-1f *M(non-trivial solution) , C 2 - C 4 ~1~A

* '

~VHHi.~fIJ\ :

I ^{sin i3 le} ^{sinh p,e}I₌₀ ₍₂₁₎

-i3~sini3,P i3;sinhi3,e

(Pi +i3~)sinJ3lfsinhi32f =0 (22)

156

(3)

~mAa~~.#'~M*#• • w.·~~A

~~A-.~~~.'~*.I~.M'~• •~

~~~#~.wn~~I~~'.~~~~~~1 J.&~•• ^·~~.~#^{• •}^~~~~.A#

.fa Ji tt$\(. #tr~.i1 ~:it :

cjI(z) = sinj31nz +sinhj32nz (23) 2.4 flt~#-f

*'

^~~*h~~:it

1iL~ • -"[1ltH5.1U5.:t....~t1! A : Tk = .!.{fCt)Y[AW(t»)+.!.{fct)y[BHf(t)}

2 2

n' Q ' l

+T{f(t»)T[n]{fCt)}+T

p,dZ

⁽³¹⁾

VpE ~ {f(t)Y[C]{f(t)} (32)

~i.~PITMt~~flJ.~;7;f1.A ' A(l4)· &JtY,H}

1t-A(10) , jJ¥.;#Lt;#LM#J~(Laplace trans from) , M

M• •~;7.A'~~*~~~M~._$~.·

ttA(14)ljit.fli.An;f~~1f ;

ElY"' (z,s) + I,Q²V"(z,s) + pAS'V(z,s) - I,S'V' = Q(s) (24) 4'-S2 =-')..'1-\'A!tIl!1l ;

ElY"" (z,s) + I, (Q' + ')..'lV'(z,s) -p~}V(z,s) = Q(s) (25) tlJ ffl ~.iU*"J:..:t..;7)}[9] , ~ M~~ V(z,S) : Vez,S) = AI COSQIZ + A, sin alz + A3 cosh a,z + A, sinh a,z

Q(s) - pA~,' (26) J+o/

I,(Q' + ')..')+ J(I,(Q' +')..'))' +4ElpA')..' al =

J

2EI (27)

Ia. _ J0.(Q' +/...')' +4EIpA')..' -I,(Q' +')..')

l ,- ^2EI

~Jt-rr-1,*1hVI0)1\AA(26) , JJ¥.i&i5l.-ROtlJ- ' ~

>"jdl AI ' A, ' A3 ' A4 ' 1i1\AA(26) , !!puT-*

~l*,.~MJf'.#;fhfuit G(s) :

0(5) =^Y(z,s)= ZeAl (28)

Q(s) P(')..) J+o/

rC')..)=2pA')..'(Q~ +a;)sinalfsinha,R (29) zCAl = Ca~ +a~ -2a~ cosa,() sinh o-,Csin alz

+ 2o:~ sina,esinhcl,icosalz

+(a~ +a~ - 2at cosh a,()sina"sinh a,z +2at sina,(sinha,(cosha,z

-2Ca; +a~lsina,esinha,f (30)

tt~P(A) , Z(')..) ,jklA-1Hi<~~~&t:it, J<tft~

~~~~~tt~M' ~'Hm Fortran;f£AaiJiUi~

A$ IMSL [101M1l.iE.4.~-t~,\t5.14~J.I;fkb{• • 2.5 -titt;;,#-f

*'

~.:t...1iH1t~.~ :B-tiT

{)'UU~~*(Assume-modc-method) [11]k-ft

X~:i!~ ~ .~(continuous system).tft~t1t. 'fit m~;7

* '

^Jt-^~^i>-.l¥-^{l!. '}^11-

*'

^{tA. A 2p}ffl B;t , 4'-

*'

^iA.^~^11L

~ ~:itA-1JlZ--#r-~tt(admissible function)*"J.;.Ct!!

it~zatM.f1l ~~ 1& -f($\(.~}l n±A~ tifutt)a0~:tH!L 1;- , ~it liE, it )~ Ji.. ~ fA. ~1ar i!1f-11} 14-(geometri c boundary conditions)Ji1i ~ f.A.~e.Il1(domain)P'J t"&~

P*·~~M~Mz*A&.&ttt~~M~*~.

• • • •tt~~~8~• •t,~x~.:t...A.~0

-f(participation factor) ⁰ ~HI:.-1~H~.:t...1Uj;~ii1't A fit t ^Atf . ^,jl}~Jll Lagrange' equation ' !lp or ^1~

~a

••

;7.A·~*M.fflaiJ.~'~*_.~

~~.:t..• • 'MI.:t...4:itM~~,-.~.m~4 i( , -=-JlJ &~.~Jt t(BesseIMit* . li;!r,otr1OJ 4fJii

~+~~;7;f1.A~I-*-'11-2_3n~~m*,~• •

;7;f1.i(&i!rr-.#M$~.:t..• •~.,~~~.~

.iiiA.~_-f*,~~.~'~~fflM~~.~~•

• •§'~~~~Ji~.&aiJ~.M*·.a• •

~'d,~.Z1U~ ~ tt~ ;

V(Z,t)=

I

^IjIn(z)fn (t)= [1j1(z)J'[f(tl] (3])

n=-!

J't-0/ IjInCz) A-tlt~_flJ.:t...~ n 1~.~ • .;J:.ttka1r, PIT:?;jU'1l a~fH"t~{~ ~ , fn(t) M Att ~ n 11iJf~ .~BtM

*.0-} '1\Afi5f-A(4)' (5) . 3t,t·~7J:.if:;5"(ol

J't-o/

fA" l:AM.dZ

B= fI.IjI~IjI~dZ

e ⁼

lElljI~~<dZ

⁽³³⁾

I

^D⁼^fI,1jJ~ ljI~dZ o

~ffflE.i:.i(31)' (32)1\ALagrange'sequation' ~1l -tit~*,.~A.~.~;7.A :

C[A]+ [Bb{t'ct) J+ dej-n'[D]){f(t)} = Q(t) (34) iJt A ~ -oT it IJ/. ' (c1^A^fibt;;, .:t... 5Jf-'Ii. ^'ij;J1 (elastic stiffness) , . .(E1f\_~at ' .% 3=.-t-i1±:.1iL ' 1]2. 'f#~-Ei

'liA.Htat' ,t;tt',->fJPfj-:i1~.:t...Q2[D14r1!'l:!lin~~HtI!'iJJ

J1 ⁰ I]l,..tf.

a

11"l4'- :

{f(t)} = {f lsin(wt +0) (35)

,J:t;t~~J'Ij-* fn(t) (n=1,2._ .. ,N):tf'1't~-!V.ilflJ· J.;l.th ff.~~.:t...'~~~~• •~'J'Ij-*MI*.~tt :tf'~~~~~M• • '~l.~~·E~Otl• .:t...~

H;'3E:f:~.t~· J'~1f;j!fA(35) 1-\.AA(34)3E.t,~1!f<.

Q(t) 0/ ~~l :

[- w' c[A]+ [BD + dcl-Q'[DDXf(t)} = 0 (36)

~/t-*.1f ~~j!IJ ;

j-W'c[Aj+[BD+ C[Cl-Ql[DDj =0 (37)

J:. A llr A • • ~ M, .:t... lPJi $. ;7 f~ A (frequency equation)' ;fJi$.;7;f1.A~Mw ' !lpA~.~!t~I!J~lPJj

$.

.E-. ' .*?i *~ ttt~

i>-.~~.~£• • •~~M,~*~.ff.·3t

.~• • *,M,.:t...M• • •~M.~~fi*,~~#

~~#&• • •mJ:..:t...lo~m• •*~~.~

~~~,Mt.*,~• •;7;f1.A'3E-*-~*,~~M~

.n8~tt·~~.~*,~.:t..._8~:it~~~~.

*~#J+*,~#tl'A~.a~~~B~~.-l

**.'~*,~~~~m~.14~Mtl·~~'3E~

M• •~*-*-*.t;;,*,~~• •;7ftA'~#• • .~~~~I!J~• •*~.#~·i>-.~~A.ft£

.~~#;7*~• •~.~~#.tt~.t;;,.:t..~

~• • • •~#~.*,~.:t...~.~~~«~.

~'&.*1i• •~Ml~~t*.'1iL;f£.m J:.~~~~.~M.~.~M.~~M.a~&

&6A££I'*~~~~~~Al.'~~~~

A. J'a °T:l J1&:Ilg!Io A. ~t.:t... § tJiJ •

157

(4)

~~~

1. Chives, D. R. and Nelson, H. D., "The Natural Frequencies and Critical Speed of a Rotating, Flexible Shaft-Disk System" , ASME Journal of Engineering for Industry, pp. 881-886,(1975).

2. Ozguven, H Nevzat, "On the Critical Speed of Continuous Shaft-Disk System" , ASME Journal of Vibration, Acoustics, Stress, and Reliability in Design, 106, pp. 59-61(1984).

3. Nelson, N. D. and McVaugh, J. M., "The Dynamics of Rotor-Bearing System Using Finite Element" , ASME Journal of Engineering for Industry, pp. 593-600,(1976).

4, Child, D. Wand Graviss, K" "A Note on Critical Speed Solutions for Finite-Element-Based Rotor models" , ASME Journal of Mechanical Design, 104, ppA12-416(1982).

5. Lund, J. W, "Stability and Damped Critcal Speeds of a Flexible Rotor in Fluid-Film Bearings" , ASME Journal of Engineering for Industry, 96, pp. 593-600,(1976).

6. Chiau, S. Wand Huang, S, c., "On the Flexural Vibrations of Rotor System Using a Modified Transfer Matrix Method" , Proceedings of the 6^th Chinese Society of Mechanical Engineering, pp.

1607-1618(1989),

7, Meriam, J. L. and Kraige, L. G., "Engineering Mechanics Dynamics" , Fourth Edition ,1998, pp.90-91.

8. L. Meirovitch, "Dynamics and Conyrol of structures",John Wiley and Sons,NY, 1992.

9. E. Kreyszig, "Advanced Engineering Mathematics", Seven Edition ,1993,pp.141-146.

10. IMSL Fortran Library Version 5,0," IMSL Fortran Library User's Guide Stat/Library Volume 2 of2", http://www.vni.com.tw/tw Ibooks! docs!

11. L. Meirovitch, "Elements of Vibration Analysis", Fourth Edition ,McGraw-Hill, 1 986,pp.282-285 , 12. M, Lalanne and G.. Ferraris, "Rotordynamics

Prediction in Engineering", John Wiley and Sons,NY, 1990

13. Andrew D. D. and Stephen A. P., "Analytical Methods in Rotor Dynamics" ,Applied Science Publishers ,1983

14, Fredric F. E., "Handbook of rotordynamics" , McGraw-HiIl,l992

158