17

31 17 1387 3 4

3 2 1

(1387/11/19 : 1386/07/04 :)

of

SI

. .

.
.
.
.
:

ive

Non-linear Vibrational Analysis of an Axially Moving Drillstring

System with Time Dependent Axial Load and Velocity in an
Inclined Well
S.M. Sahebkar, M.R. Ghazavi, S. Emailzadeh-Khadem, and M.H. Ghayesh
ABSTRACT

ch

Mech. Eng. Dept., School of Eng., Tarbiat Modarres Univ.

Ar

A non-linear model of a drill string system in an inclined well with axially moving motion and loading is
established and analyzed by perturbation technique. The drillstring is modeled as a simply supported axially
moving rotor. Based on Hamilton's principle, the dynamic model developed consists of kinetic and potential
energy. Non-linear coupling terms are kept in the formulation, which leads to full coupling between the axial
and the transverse vibrational modes. The method of multiple scales is employed to solve the non-linear
equations, to analyze the steady state response, and to find the stability region of the system. Analytical and
numerical results reveal the non-linear phenomena, such as primary and parametric resonance that have not
been previously reported. Our results are compared with some existing experimental date showing relatively
good agreements.

Key Words: Drillstring, Non-linear Vibration, Axially Moving Rotor, Method of Multiple Scales, Dynamic
Stability

[email protected] :( ) -
-
-
-

www.SID.ir

 4 3 1387

18

-1

] .[10-11

].[4-5, 12-14

] .[1-3

].[15

SI

2 ] .[4

] .[16

].[17-19

of

] [20 .

] [5

] .[16

] [6

] 17 [21

] [22

ive

. -

ch

] .[2

].[7

Ar

] .[8-9

)1-Rate of Penetration (ROP

)2-Bottom Hole Assembly (BHA

3-Jeffcott

www.SID.ir

19

...

" "

) (BHA

-2

SI

] .[11 BHA

.

of

. ) (BHA

1 .

BHA

x .

ch

ive

BHA

BHA .

Ar

1 xyz

G . C O

)( .

e0

P l

u w v
O yz .

)1-Weight on Bit (WOB

www.SID.ir

 4 3 1387

20

) (t ) (ROP
.

]:[23
t2

(T + Wc + Wnc ) = 0 ,

)(1

t1

T
Wc

Wnc . Wc
Wc = U : U
.
]:[16 11

2
u

u
+ (1 + ) +

t
x

l v
v
1

+
+ e (x ) cos(t ) + dx
T = A
0
x
2

t
2
w
w

+
) e (x ) sin(t

x
t

2
w v
2

1 l
+ 2J p
( +Jy
) + J z ( ) 2 dx
x t x
t
t
2 0

SI

)(a

)(2

of

+ 2J p

l A
.

A-A .

ive

Section A-A
)(b
) :(a) :(1 ) :(b

2
1 l 2J p w v dx

x t x

] .[11 ) (2 y

) ( = w x z

ch

) ( = v x.

] 15 :[24

Ar

) (3

Dh2 + Do2
),
Dh2 Do2

( Di2 + Do2

A = d A + f

( D D ,),
2
i

2
o

=A

4
D0 Di

Dh d f
.

) :(2

.

www.SID.ir

e( x) = e0
J = J y = J z J p = 2J J p

.[11] J p = d ( Do4 Di4 ) 32

21

...

.
]:[231621
= ) (Wc + Wnc
2w
2v
1 l
EAexx2 + EI ( 2 ) 2 + ( 2 ) 2 dx

x
x
2 0

)(4

1 l
[ P + q(l x) cos( )] exx dx

0
2
,
1 l

q sin( )v dx + Wh

2 0

)(9

) (2 )
(4 ) (1
.
) (8
:

SI

E I
y A z P q

) .(5

q = A( d f ) g

)(5

)(10

) q sin( u

* z * z .
) (10 :

)(8

]:[16

www.SID.ir

) P + ql (1 x* ) cos(
EI
*
*
,
= I = 2 ,P
P0l
P0

ive

x u w
,u = , w = ,
l
l
l
P0
v
A
, =
.
v = , t * = t
2
l
P0
Al
= *x

J
EA
= 2
= *, J
= ,
,
2
Al
2 P0
l AP0

)(11

f CD Dol
= ,C
,
= e
2 A
P0
) ql cos( ) * ql sin(
*
Ale0 2

*
D

ch

CD .

. :

J p

P0

= , q2

P0

= q1

* z
]:[21
)(12

Ar

)(7

v
w
f CD Do ( ) 2 + ( ) 2

1 l
t
t
dx,
Wh =
0
2 w
v

* ( ) w + ( ) v
t

t

* 4 z
* z * z * z
* z
* *+ CD* ( * )( * )( * ) q1
*4
x
t t t
x
* 2
*
*
1
z
z z
= e* exp(it ) + *2 ( * )( * )dx* iq2* ,
x 0 x x

z * = w* + iv*

v Wh . .

) ( exx
]:[16
u 1 w 2 1 v 2
)(6
= exx
+ ( ) + ( ).
2 x
x 2 x

]:[11

* 2 z
* 2 z * * z
* 2 z
* * * + *2 *2 + * * + 2
*2
t
x
t x
x t
* 3
* 4
* 2

z
+i 2 * *2 J * *2 *2 + P* *2
t x
t x
x
* +I

of

)(4

]) [ P + q(l x) cos(

* 1 1 w* 2 v* 2
) + ( * ) dx .
(
x
*2 0 x

= *ex* x

* 2 z
* 2 z
(0, t ) = *2 (1, t ) = 0.
*2
x
x

= ) z* (0, t ) = z * (1, t

) (P WOB
) (t
] 21 :[22
)(13

P (t ) = P0 + P0 sin( p t ),

)(14

(t ) = 0 + 0 sin( t ),

P0
WOB .
p ,

 4 3 1387

22

.
) (111314
:
)(15

P* = 1 + sin( p t ) + q1* (1 x* ),

= 0 + 0 sin( t ).
)(16

A
P0

)(17

0 = 0

) O (

] 17 . z * = z : [21

2
2 z
2 z

+
+

(sin

)
t

0
0
t *2
x*2
z
3 z
+ 0 cos( t ) * + i 2 * *2
x
t x
2
z
4 z
+2 0 + 0 sin( t ) * * J *2 *2
x t
t x
2
z
+ 1 + sin( p t ) + q1 (1 x* ) *2
x
4 z
z z z
z
* + I * *4 + CD ( * )( * )( * ) q1
x
t t t
x
= e exp(it ) i q2

* z
) dx .
*x

SI

) :(1
].[311
Drillstring
Drilling Fluid
E = 210 GPa
f = 1500 kg/m3

Bore hole

Do = 0.2286 m
)(9 inch
Di = 0.0762 m
)(3 inch
e0 = 0.00508 m
)(0.2 inch
)l = 10 m (65 ft
0 = 0.01 m/s

ch

n = 0,1, 2,..... ,

Tn = nt * ,

2 = 2 , e* = 2 e , CD* = 2 CD.
) (10
:

Ar

] 3 .[11

q2* = 2 q2 , q1* = q1 , J * = J

z ( x, t , ) = z0 ( x, T0 , T1 ,....) + z1 ( x, T0 , T1 ,....) + ...

www.SID.ir

) (17
:

.

.
] 25 [26 .
] 25 [26 )(17
:
)(18

P0 = 100 kN
= 5 rad/s
)(47.7 rpm

-3

ive

Dh = 0.4445 m
)(17.5 inch
= 0.1 rad

( x

z
x*2

of

d = 7850 kg/m3

()

T0 , T1 , T2 ,.... T0 .
T1 , T2 ,...
:
T T1

= 0
+
+ ... = D0 + D1 + ...
t * t * T0 t * T1

)(19

= , D1
.
T0
T1

= D0

2
= D02 + 2 D0 D1 + ...,
t *2

) (18 ) (19
:
)(20

23

...

O ( 0 ) :

2 z0
z
4 z
+ 2 0 D0 0* + I * *40 = 0
*2
x
x
x

)D02 z0 + ( 0 2 + 1

O( ) :
1

2 z1
z
4 z
= + 2 0 D0 1* + I * *41
x*2
x
x
2

z
z
*2 D0 D1 z0 2 0 2 sin( t ) *20 0 cos( t ) 0
x
x
z
z
2 z
2 0 sin( t ) D0 0* 2 0 D1 0* i 2 D0 *20
x
x
x
2
2
2 z0
2 z0
* z0
+ JD0 *2 sin( p t ) *2 q1 (1 x ) *2
x
x
x
z0
2
) + q1 * CD D0 z0 z0 D0 z0 + e exp(it
x
* 2 z0 1 z0 z0
+ *2 * * dx iq2 .
x 0 x x

)(26

) (12 )(25
:

)D02 z1 + ( 0 2 + 1

)(21

SI

)(22

)(27

. Cn 2 , Cn 3 , Cn 4 ) (27
) (25
:
= )Z n ( x

z0 ( x, t ) = [ An (T1 ) exp(inT0 ) ] Z n ( x),

n

n n
An (T1 ) .

)I * Z n'''' ( x) + ( 0 2 + 1) Z n'' ( x

) ( )(e e

1 +
in 3
in 2
( )(e e ) in 4 x

+
e .
) ( 2 2 )(ein 2 ein1

n4
n1
+ 2

i
i
2

) (n 4 n3 )(e n 2 e n 3

in1

in 3

2
n1
2
n2

2
n4
2
n4

ch

z1
. ) (22 ) (20
:

)(28

in1x (n24 n21 )(ein 3 ein1 ) in 2 x

e
e 2
) (n 4 n22 )(ein 3 ein 2

Cn1
in 2
in1
2
2
(n 4 n1 )(e e ) ein 3 x
) (n24 n23 )(ein 2 ein 3

ive

Z n (x ) ) An (T1 T1

1
1
1
1

2
2
2

1
2
3
4
n
n
n
n

) exp(in1
) exp(in 2
)exp(in3
) exp(in4
2

2
2
2
) n1 exp(in1) n2 exp(in 2 ) n 2 exp(in3) n2 exp(in 4
1
0
C
0
2
n
= * Cn1
Cn3
0 .

C
0
n4

Cn 2 , Cn 3 , Cn 4

of

) O( 2
. ) (20
]:[25

I *ni4 ( 0 2 + 1)ni2 20nni n2 = 0.

* x
]:[21

n . .
) (28 ) (22 ) (21

=(...)

) sin( t ), cos( t ) sin( p t

)(23

)(24

),(...
=
,
*x
x*2
3
4
(...)'''= *3 , (...)''''= *4 .
x
x

) (23 ]:[21
Z n ( x) = Cn1[ein1 x + Cn 2 ein 2 x

)(25
+ Cn 3ein 3 x + Cn 4 ein 4 x ].
) (25 )(23
:

www.SID.ir

Ar

+2i 0n Z n' ( x) n2 Z n ( x) = 0.

n ni ) (26
) (27

= D z + z + 2 0 D z + z + I z
'''' *
1

)(29

''
1

'
0 1

'' 2
0 1

2
0 1

) n1 exp(inT0 ) + n 2 exp(i (n + )T0

) + n 3 exp(i (n )T0 ) + n 4 exp(i(n + p )T0
+ n 5 exp(i (n p )T0 ) + e exp(it ) iq2 .

n 5 , n 4 , n 3 , n 2 , n1
.

 4 3 1387

24

.
) (29
) exp(inT0 ) exp(inT0 .


.[7] = n :
n .

.

.
n=3
n=1 ] .[7

:

An Z n An Z n.

-4

n1 , n 2 , n3 .

)(31

1
an (T1 ) exp [i n (T1 ) ] ,
2

= ) An (T1

an (T1 ) ) n (T1
T1 . ) (31
) (30
:
)(32

an
1
+ Re( n1 )an + Re( n 2 ) an3
4
T1

SI

1
+ Re( n 3 ) an2 = 0,
2
n
1
+ Im( n1 )an + Im( n 2 ) an3
an
4
T1

)(33
1 2
+ Im( n3 ) an = 0,
2
Re Im

of

-1 : n
) (= n 0 2n

) A(T1 ]:[2526

.
an n ) (32 )(33

) (= n 0 2n

1 ) (32 )(33

-3 ) (= n 2n

) (= n 0 2n .

] [26


) A(T1 :
)(30

) An (T1
) + n1 An (T1 ) + n 2 An2 (T1 ) An (T1
T1
+ n 3 An (T1 ) An (T1 ) An (T1 ) = 0,

www.SID.ir

1 2

Im( n1 ) + Im( n 2 ) 4 an 0
)(34
n =
T1 + n 0 .
+ Im( ) 1 a

n3
n0
2

Ar

-1-4 : n

a = 0 an = an 0
'
n

ch

2n .
] [26
) (29
.

ive

-2 n

. n1 , n 2 , n 3

T1 = T0 = t
) (34 ) (18 ) (22
:
1
)an exp [i ( n + nT0 ) ] Z n ( x
2

) Im( n1

i + Im( ) a 2 t
1

n
n2
n0

= an exp
4 Z n ( x),

+ Im( n 3 ) an 0

+i

n0

= ) z0 ( x, t

)(35

25

...

1 2

Im( n1 ) + Im( n 2 ) 4 an 0
= n
.
+ Im( ) 1 a

n3
n0
2

n , nonlinear

) (38 ) (39 .
n ) (38 ) (39

) (an ,
:
1
2

4 Re( n 4 ) 2 + 4 Im( n 4 ) 2

)(40
-2-4 n

) (= n 0 2n .

SI

n
n
:
)(36
= n + ,
) (36

1 2

Re( n1 ) + Re( n 2 ) 4 an .

+ Re( ) 1 a

n3
n
2

an

) (40 "-" .
"-"
.
) ( an 0 , n 0 ) (38 ) (39
an T1 = n T1 = 0 .

.

of

n.
) (36 ) (29
:

1
4

= Im( n1 ) Im( n 2 ) an2 Im( n 3 ) an

+ n 3 An (T1 ) An (T1 ) An (T1 ) + n 4 exp(i T0 ) = 0.

n1 , n 2 , n 3 n 4 .

an
1
+ Re( n1 )an + Re( n 2 ) an3
T1
4

n
1
+ Im( n1 )an + Im( n 2 ) an3
T1
4

)(39

an

) Im( n 4 ) cos( n
1
*+ Im( n 3 ) an2 + 2
an = 0,
2
) + Re( n 4 ) sin( n

) G (an , n

an
.

) G (an , n

n

an T1 = n T1 = 0
) (39 38
an = 0 .

www.SID.ir

n
= G (an , n ),
T1

) F (an , n

a
n
= J I
) F (an , n

) (31
) A n (T1 ]:[26
1
i T
)(42
An = (rn + isn ) exp( 1 ).
2

. n = T1 n

an
= F (an , n ),
T1

Ar

)(38

) Re( n 4 ) cos( n
1
*+ Re( n 3 ) an2 + 2
an = 0,
2
) Im( n 4 ) sin( n

)(41

ch

) A n (T1 ) (31
:

.
) (
:

ive

)(37

) An (T1
) + n1 An (T1 ) + n 2 An2 (T1 ) An (T1
T1

J I = 0

) (42 ) (37

.

 4 3 1387

26

rn = sn = 0
:

= 0.

)(43

)n = 0

J I (r =s
n

. 1

:

1 = 16.154

-3-4 ) (= n 2n
2n .

.


] 21 :[27
)(44
= 2n + .

:
) An (T1
) + n1 An (T1 ) + n 2 An2 (T1 ) An (T1
T1
)(45
) + n 3 An (T1 ) An (T1 ) An (T1

2 = 65.627

. 3 = 147.81 rad/s
1 , 2
n / n .


. n=1
= n

SI

= 2n . n=3
:

of

= n 3 . = 2n 3



.
) (ROP

.
WOB .

+ n 5 exp(i T0 ) An (T1 ) = 0.

) An (T1 ) (39
:

n
1
+ Im( n1 )an + Im( n 2 ) an3
T1
4

an

) Im( n 5 ) cos( T1
1
+ Im( n 3 ) an2 +
an = 0,
2
) Re( n 4 ) sin( T1

an n ) (47 46
.
-4

.
) (26
) (27

www.SID.ir

Ar

)(47

) Re( n 5 ) cos( T1
1
+ Re( n 3 ) an2 +
an = 0,
2
) + Im( n 5 ) sin( T1

ch

)(46

an
1
+ Re( n1 )an + Re( n 2 ) an3
T1
4

ive

n1 , n 2 , n 3 n 5 .

) :(3
) .( 0
3
.

27

...

l = 10m .
5
.


.

D
SI

) :(4

of

) * .( I

-
.

ive

) :(6

) .( P0 / P0

ch

) (WOB
. WOB

Ar

.
ROP WOB . 6
.
6
1 P0 = 100kN
) :(5
) .( l l

.
.

BHA . 5

www.SID.ir

 4 3 1387

28

) :(7
) :

n1 = n 2 = n 3 = 0 :
.( = 0.01 - 1

SI

) (35

7

) :(8 -

of

) : :
(.

= 0.01 . 7

www.SID.ir

ive

) ( an ,
.
) (42
.

ch

. 8

Ar

.



) ( n1 = n 2 = n3 = 0
.

.
) (40

) ( an ,
.
] [12

.


.
)
(
].[12
] [12
) (3-4

.



.

29

...

-6

.

.



.


) (

SI

of

. ) ( P


.


.
.



.

.

.

) (
.


.

.


.


.



.

.

.

ive

ch

Ar

www.SID.ir

1. Mason, J. and Sprawls, B., Addressing BHA

Whirl: The Culprit in Mobile Bay, SPE Drilling
and Completion, Vol. 13, No. 4, pp. 231236,
1998.
2. Vandiver, J.K., Nicholson, J., and Shyu, R.J.,
Case Studies of the Bending Vibration and
Whirling Motion of Drill Collars, SPE Drilling
Eng., Vol. 5, No.4, pp. 282290, 1990.
3. Jansen, J., Non-linear Rotor Dynamics as
Applied to Oil well Drillstring Vibrations, J.
Sound and Vibration", Vol. 147, No. 1, pp.115
135, 1991.
4. Mitchell, R. and Allen, M., Lateral Vibration:
The Key to BHA Failure Analysis, World Oil,
Vol. 200, No. 4, pp. 101104, 1985.
5. Chen, S. and Gradin, M., An Improved Transfer
Matrix Technique as Applied to BHA Lateral
Vibration Analysis, J. Sound and Vibration, Vol.
185, No.1, pp. 93106, 1995.
6. Vaz, M., and Patel, M. Analysis of Drill Strings
in Vertical and Deviated Holes, Using the
Galerkin Method, Eng. Struc., Vol. 17, No.6, pp.
437442, 1995.
7. Dunayevsky, V., Abbassian, F., and Judzis, A.,
Dynamic Stability of Drillstrings under
Fluctuating Weight on Bit, SPE Drilling and
Completion, Vol. 8, No. 2, pp. 8492, 1993.
8. Leine, R.I., Bifurcations in Discontinuous
Mechanical Systems of Filippov-Type, PhD

1387 3 4

30

10.

11.

12.

9.

Dissertation, Mech. Eng. Dept., Univ. of

Eindhoven, Eindhoven, The Netherlands, 2000.
Zamanian, M., Khadem, S.E., and Ghazavi, M.R.,
Stick-slip Oscillations of Drag Bits by
Considering Damping of Drilling Mud and Active
Damping System, J. Petroleum Science Eng.,
Vol. 59, No. 1, pp. 289-299, 2007.
Shyu, R.J. Bending Vibration of Rotating Drill
Strings, PhD Dissertation, Massachusetts
Institute of Technology, Dept. Ocean Eng.,
Cambridge Univ. 1989.
Christoforou, A. and Yigit, A. Dynamic
Modeling of Rotating Drillstrings with Borehole
Interactions, J. Sound and Vibration, Vol. 206,
No. 2, pp. 243260, 1997.
Berlioz, A., Hagopian, J.D., Dufour, R., and
Draoui, E., Dynamic Behavior of a Drill String:
Experimental
Investigation
of
Lateral
Instabilities, J. Vibration and Acoustics, Trans.
ASME, Vol. 118, No. 3, pp. 292298, 1996.
Zifeng, L., Xingrui, M., Wenhu, H., and Xisheng,
L., A 3-D Analysis of a Bottomhole Assembly
under Large Deflection, SPE Drilling and
Completion, Vol. 11, No. 2, pp.104110, 1996.
Khulief, Y.A., Al-Sulaiman, F.A., and Bashmal,
S., Vibration Analysis of Drillstrings with Selfexcited Stick Slip Oscillations, J. Sound and
Vibration, Vol. 299, No. 3, pp. 540558, 2007.
Heisig, G. and Neubert, M. Lateral Drillstring
Vibrations in Extended-Reach Wells, Paper
IADC/SPE 59235, IADC/SPE Drilling Conf.,
New Orleans, LA, 2000.
Wickert, J.A. and Mote, C.D. Classical Vibration
Analysis of Axially Moving Continua, J. Appl.
Mech., Vol. 57, No. 1, pp. 738744, 1990.
Pakdemirli, M. and Ulsoy, A.G., Stability
Analysis of an Axially Accelerating String, J.
Sound and Vibration, Vol. 203, No. 2, pp. 815832, 1997.
Pellicano, F., Fregolent, A., Bertuzzi, A., and
Vestroni, F., Primary and Parametric Non-linear
Resonances of a Power Transmission Belt:
Experimental and Theoretical Analysis, J. Sound
and Vibration, Vol. 244, No. 4, pp. 669684,
2001.
Chakraborty, G., Mallik, A.K., and Hatal, H.
Non-linear Vibration of Traveling Beam, Int. J.
Non-linear Mech., Vol. 34, No. 2, pp. 655670,
1999.
Ghayesh, M.H. and Khadem S.E. Non-linear
Vibration and Stability Analysis of a Partially
Supported Conveyor Belt by a Distributed
Viscoelastic Foundation, J. Struc. Eng. Mech.,
Vol. 27, No. 1, pp. 17-32. 2007.
Oz, H.R. On the Vibration of an Axially
Traveling Beam on Fixed Supports with Variable
Velocity, J. Sound and Vibration, Vol. 239, No.
3, pp. 556564, 2001.

SI

22. Marynowski, K. and Kapitaniak, T. Zener

Internal Damping in Modeling of Axially Moving
Viscoelastic
Beam
with
Time-dependent
Tension, Int. J. Non-linear Mech., Vol. 42, No. 1,
pp. 118-131, 2007.
23. Meirovitch. L. Analytical Methods in
Vibrations, New York, Macmillan, 1967.
24. Harris, C.M. and Piersol, A.G. Harris Shock and
Vibration Handbook, Fifth Ed., McGraw-Hill,
New York, 2002 .
25. Nayfeh, A.H. and Mook, D.T. Non-linear
Oscillations, John Wiley, New York, 1979.
26. Nayfeh, A.H. and Balachandran, B., Applied
Non-linear Dynamics: Analytical, Computational,
and Experimental Methods, John Wiley, New
York, 1995.
" " .. .27

.1387

13.

of

14.

Ar

ch

ive

15.

16.
17.

18.

19.

20.

21.

www.SID.ir

31

...

n2

'
n

'
n

'
n

{ Z C i Z ( x)Z ( x)Z ( x)dx}

=
{ Z 2i Z ( x) 2 Z ( x) dx}
1

n3

2
n

'
n

{ Z [e ] dx}
1

n4 =

{ Z 2i Z ( x) 2 Z ( x) dx}
1

'
n

2 ''
'
1 i 0 Z n ( x) 0 Z n ( x)
2
dx
0 Z n

'
''

+n 0 Z n ( x) iZ n ( x)

n5 = 1
'
Z n 2in Z n ( x) 20 Z n ( x) dx

Ar

ch

ive

n 5 = iA n (T 1 ) Z n'' (x )

''
n

of

A n (T 1 )
A n (T 1 ) '

2i n T Z n (x ) 2 0 T Z n (x )

1
1

q1 (1 x * )A n (T1 ) Z n'' (x ) + q1A n (T 1 ) Z n' (x )

2
''
2
''
n 1 = +n A n (T1 ) Z n (x ) J n A n (T1 ) Z n (x )

2
C D i n n A n (T1 ) Z n (x )A n (T1 ) Z n (x )A n (T 1 ) Z n (x )

1
2
''
'
'
*
+ A n (T 1 )A n (T1 ) Z n (x ) 0 Z n (x ) Z n (x ) dx

1
2

''
'
i A (T ) Z (x ) 0 A n (T1 ) Z n (x )
n2 = 0 n 1 n
2

A (T ) Z ' (x )

n 0 n 1 n

2
''
'
i 0 A n (T1 ) Z n (x ) 0 A n (T1 ) Z n (x )
n3 =
2

+ A (T ) Z ' (x )

n 0 n 1 n

''
n 4 = iA n (T1 ) Z n (x )

{ Z Z ( x) Z ( x )Z ( x ) dx dx}
=
{ Z 2i Z ( x) 2 Z ( x) dx}

SI

1 q1 (1 x* ) Z n'' ( x) + q1Z n' ( x)

dx
0 Z n
2 ''
2 ''
+n Z n ( x) J n Z n ( x)

n1 =
1
'
Z n 2in Z n ( x) 20 Z n ( x) dx

www.SID.ir