1
1
旋轉圓盤受週期性外緣力時之參數共振研究
Par ametr ic Resonance of a Spinning Disk Under Pulsalting Edge Loads 計畫編號: NSC 87-2212-E-002-042 執行期限: 86年8月1日至87年7月31日 計畫主持人: 陳振山 jschen@w3.me.ntu.edu.tw 執行機關: 國立台灣大學機械系 摘要 本計畫以解析的方法研究旋轉圓盤受外緣 週期性外力時之參數共振現象。我們假設 徑向外緣力的分佈可以傅立葉級數的方式 展開。利用圓盤特徵函數的正交特性,我 們將偏微分方程式離散化成一組一階的希 爾方程式,再利用多重尺度法決定組合式 參數共振發生的條件。 關鍵詞:參數共振,旋轉圓盤,外緣力 Abstr act
The parametric resonance of a spinning disk under a space-fixed pulsating edge load is investigated analytically. We assume that the radial edge load can be expanded in a Fourier series. With use of the orthogonality properties among the eigenfunctions of a gyroscopic system, the partial differential equation of motion is discretized into a system of generalized Hill’s equations in the first order form. The method of multiple scale is employed to determine the conditions for combination resonances to occur.
Key Wor ds: Parametric Resonance, spinning disk, edge load
Intr oduction
The vibration analysis of a spinning disk under a space-fixed edge load attracts attention because of its possible application in such fields as circular saw cutting and grind wheel operation. Recently Chen (1994) formulated this problem with emphasis on the effects of relative motion between the disk and the edge load on the stability and natural frequencies of the loaded disk. The edge loads in this work are assumed to be independent of time. In the case when the
edge loads are periodically varying, parametric resonance may be induced. Tani and Nakamura (1978; 1980) studied the dynamic instability of an annular disk under periodically varying in-plane edge traction. The disk considered in their papers is clamped on both the inner and outer radii, and the edge traction is uniform in the circumferential direction. Zajaczkowski (1983) investigated the parametric resonance of a clamped-free disk under both uniform and concentrated pulsating torques. In these papers both the annular disk and the periodic loading are fixed in space.
A natural extension of these previous analyses is to study the dynamic instability of a spinning annular disk under periodically varying edge load which is fixed in space. This investigation may find application in the wood cutting industry as it represents a more general model for the cutting process by a circular saw.
Equation of Motion
We assume the disk is clamped at the inner radius and subjected to a periodic radial traction at the outer radius. The edge in-plane traction can be expanded in a Fourier series
ε
cos
γ
t
f
kcos
k
θ
k=
∞
∑
0
. The disk is rotating with constant speed Ω, while the edge load is fixed in space. The dimensionless equation of motion of the system, in terms of the transverse displacement w and with respect to the stationary coordinate system (r,θ), is
M
w
t
G
w
t
K
K w
∂
∂
∂
∂
2 2+
+
(
+
$
)
=
0
(1) where2
2
M
≡
1
,G
≡
2Ω
∂
∂θ
,K
≡
Ω
+ ∇ +
L
2 2 2 4∂
∂θ
,∃
cos
∃
(
) cos
K
t
f
L
r
k
r
k k k≡
+
−
= ∞∑
ε
γ
δ
θ ∂
∂
1
0Equation (1) can also be cast in the first-order operator form
(
A
+
A
∃
) (
∂
x
−
B
+
B x
∃
)
=
∂t
0
(2)by defining the state vector
x
≡
∂
∂
w
t
w
and the matrix differential operators
A≡ A B B ≡ ≡− − ≡ − M K K G K K K K 0 0 0 0 0 0 0 0 , ∃ ∃, , ∃ ∃ ∃
Figure 1 shows the dimensionless natural frequencies
ω
mn of a freely spinning disk as functions of dimensionless rotation speedΩ
.Discr etization
Since a closed form solution for Eq.(1) does not exist in general, we use an expansion in terms of finite number of eigenfunctions of the freely spinning disk to approximate the true solution of Eq.(1),
x
( , ; )
r
t
c
pq( )
t
x
pq( , )
r
q N N p Nθ
≈
θ
=− =∑
∑
2 2 1 0 (3)N
1 andN
2 are the maximum numbers of nodal circles and nodal diameters, respectively, of the modes used in the expansion. Substituting Eq.(3) into (11) and taking the inner product betweenx
mn and both sides of Eq. (1), we obtain a system of generalized Hill equations in the first order form, dc dt i c t A A dc dt B A c mn mn mn pq k mn mn mn pq pq kmn mn mn pq q N N p N k − + − = =− = = ∞ ∑ ∑ ∑ ω εcosγ ∃ ( ) ∃ ( ) 2 2 1 0 0 0 (4) whereA
mnmn=
4
mn mn+
n
∫
R
mn2r rdr
1πω
ω
η(
Ω
)
( )
Whenn q
− ≠ ±
k
,A
$
pq kmn( )=
B
$
pq kmn( )=
0
. On the other hand, whenn q
− = ±k
, we obtain$
$
( ) ( )A
pq kmn=
A
mn kpq=
=α σ σθ σθ σθ η k rk mn pq k r k mn pq r k mn pq r dR dr dR dr q n r d dr R R n q R dR dr dr [~ ( ~ ~ ) ( )~ ] + − − + ∫ 1$
(
) $
( ) ( )B
pq kmni
A
mn pq pq k mn=
ω
+
ω
where
α
0=2π
f
0, andα
k =π
f
k whenk
≠
0
. .Per tur bation Technique
The method of multiple scale assumes an expansion of the form
c
mn( )
t
=
c
mn( )0( , )
t T
1+
c
mn( )( , )
t T
+
O
(
)
1 1
2
ε
ε
(5)where
T
1≡ε
t
. Substituting (5) into (4) andequating coefficients of like powers of
ε
yieldε
0 :D c
0 mni
mn mnc
0 00
( )−
ω
( )=
(6)ε
1 : D c i c D c t A A D c B A c mn mn mn mn pq k mn mnmn pq pq k mn mnmn pq q N N p N k 0 1 1 1 0 0 0 0 0 0 2 2 1 ( ) ( ) ( ) ( ) ( ) ( ) ( ) cos $ $ − = − − − =− = = ∞ ∑ ∑ ∑ ω γ (7) whereD
t
0≡
∂
∂
, andD
1T
1≡
∂
∂
. The generalsolution of Eq.(6) can be written in the form
c
mn( )0H
mn( )
T e
i mnT1
0
=
ω (8)Substituting (8) into (7) yields
D c
mni
mn mnc
D H e
mn i mnT 0 1 1 1 0 ( )−
ω
( )= −
ω + + + − =− = = ∞∑
∑
∑
ωmn pq k ω γ ω γ mn mn mn pq i T i T q N N p N k A A H e e pq pq 2 0 0 2 2 1 0 0 $ ( ) ( ) ( ) (9) In the case whenγ
is nearω
rs+ω
uv, we assume thatγ ω
= rs+ω
uv+ες
, whereς
is a detuning parameter. The secular terms in Eq.(9) are eliminated if3
3
D H
i
A
A
H e
rs rs u v k rs k rs rs uv i T 1 02
0
1−
= −=
∞∑
ω
ς$
, ( ) (10)D H
i
A
A
H e
uv uv r s k uv k uv uv rs i T 1 02
0
1−
= −=
∞∑
ω
ς$
, ( ) (11)It follows that the solution is bounded if and only if
ς
2≥ Λ
uv rs whereΛ
uv rs rs uv u v k rs k rs rs uv uvA
A A
=
− = ∞∑
ω ω
∃
, ( ) 0 2 (12) On the other hand, combination resonance may occur whenς
2< Λ
rsuv. Frequencyω
rs+ω
uv is called the center frequency andΛ
uv rsis called the width parameter of the parametric resonance.
Concentr ated Edge Loading
We consider the case when the spinning disk is subjected to a space-fixed concentrated edge load. In this case the Fourier decomposition of the edge load would include all the harmonics of
cos
kθ
, and the summationA
$
urs, v k( )k − = ∞
∑
0 in Eq.(12) reduces to a single non-zero termA
$
u, v k( )rs
− ,
where
s v
+ = ±
k
. In other words, for any two modes (r,s) and (u,v), there always exists a specific Fourier component which renders the summationA
$
urs, v k( )k= −
∞
∑
0
non-zero. Bearing the relation
ω
mn= −
ω
m n,− in mind and speaking of only positive naturalfrequencies, we can conclude that
combination resonance of the sum type occurs when both modes are non-reflected or both modes are reflected. On the other hand, combination resonance of the difference type can occur only when one mode is reflected and the other is non-reflected. Single mode parametric resonance can occur for any mode, reflected or non-reflected.
The center frequency and the width parameter of a spinning disk under a pulsating concentrated edge load are shown in Fig.2. Solid lines represent the single mode parametric resonance, while dashed lines and dotted lines represent the cases of combination resonance of the sum type and difference type respectively. In the frequency range of Fig.3 only the modes with zero nodal circle contribute to the parametric resonance. The mode labels in Fig.2 are simplified by neglecting the number of nodal circle. For instance,
0 2
+
frepresents the combination resonance of the sum type involving modes
( , )
0 0
and( , )
0 2
f, and 2b−3r represents the combination resonance of the difference type involving modes ( , )0 2 b and ( , )0 3 r. It can be seen from Fig.2 that the rotation speed tends to squeeze the width of parametric resonance region.In the special case when the excitation frequency
γ
approaches zero, which corresponds to the case of constant concentrated edge load, we can see from Fig.2 that there exist two rotation speeds at which parametric resonance can occur. The rotation speed at which the solid line labeled3b+3b intersects the horizontal line
ω
rs+ω
uv=0 corresponds to the firstcritical speed. Therefore, the divergence instability reported by Chen (1994) can be considered as a limit case of single mode parametric resonance. On the other hand, the rotation speed at which the dotted line labeled 2b−3r intersects the line
ω
rs+ω
uv=0 is a speed at which thebackward wave ( , )0 2 b and reflected wave
( , )0 3 r become degenerate. Therefore, the flutter instability reported by Chen (1994) can be considered as a limit case of combination resonance of the difference type.
Conclusions
Dynamic stability of a spinning annular disk under periodically varying in-plane loading on the outer rim is studied analytically. The results show that combination resonance is possible only when there exists a specific Fourier
4
4
component
cos
kθ
in the edge load, where k equals the sum of the number of nodal diameters of these two modes. Sum type resonance occurs when both modes are non-reflected or both modes are non-reflected. On the other hand, difference type resonance occurs when one mode is reflected and the other is non-reflected.Refer ences
Chen, J.-S., 1994, "Stability Analysis of a Spinning Elastic Disk Under a Stationary Concentrated Edge Load", ASME Journal of Applied Mechanics, Vol.61, pp.788-792.
Tani, J., and Nakamura, T., 1978, ynamic Stability of Annular Plates Under Periodic Radial Loads”, Journal of Acoustical Society of America, Vol.64, pp.827-831.
Tani, J., and Nakamura, T., 1980, ynamic Stability of Annular Plates Under Pulsating Torsion”, ASME Journal of Applied Mechanics, Vol.47, pp.595-600.
Zajaczkowski, J., 1983, tability of Transverse Vibration of a Circular Plate Subjected to a Periodically Varying Torque”, Journal of Sound and Vibration, Vol.89, pp.273-286.
Fig.1 Dimensionless natural frequency versus dimensionless rotation speed for a freely spinning disk.
Fig.2 Center frequencies and width parameters of a spinning disk under a space-fixed pulsating concentrated edge load.