l}v=mt p=Ut QO Rt=}r}w |QU }=Q QkDsjce.journals.sharif.edu/article_4560_317ad6e910bcd9ef2aa97cd65b3... · 2 / |xQ =t u =DU @ = h} Q u= Qt a | U Ovy t " ^x; ; ( ; ); ; ( ; ) # "; ;

�1396

u=DU@=D�

'h}QWu

=Qta|U

Ovyt97

�87"X

'2/2|

xQ=tW'

33�2

|xQwO

Ori

gina

lArt

iclel}v=mt p�=Ut QO Rt=}r} w |QU ?}=Q[ ?} QkD

xO=iDU= =@ |]N |a=HDQ= CUmW=RHt CqO=at VwQ R=

�|QDmO |wHWv=O� |v=OR} |Oyt�O=DU=� �|H=N QY=v

TQOt C}@ QD x=oWv=O 'CU} R \}Lt w u=Qta |UOvyt |xOmWv=O

"CU=yxR=U QO Xkv w lQD OwHw '=yxR=U |L=Q] w p}rLD QO p�=Ut u} QDsyt R= |m}p�=Ut pL 'wQu}= R= &OvDU}v pL p@=k |r}rLD CQwYx@ 'OvQ=O lQD xm =yxR=U R= |Q=}U@

Q=DWwv "CU= xOW p}O@D syt p�=Ut R= |m} x@ |OOa |=yVwQ =@ CUmW l}v=mt'CU= xDN=OQB =RHt CqO=at VwQ s=v x@ O}OH |r}rLDxt}v VwQ l} |xaUwD x@ Q[=L|a=HDQ= u=O}t ?}=Q[ '|]N |a=HDQ= CUmW l}v=mt |x} Q_v R= xO=iDU= =@ u; QO xm

4 |x@U=Lt =@ w x@U=Lt 'OvQ=O CQyW Rt=}r} w |QU ?}=Q[ x@ xm 'lQD lwv |@v=Ht=@ 'xt=O= QO "CU= xOW xOR ?} QkD VvD w |}=Hx@=H u=O}t 'Rt=}r} w |QU pw= ?} Q[|=ywQ}v Q=OQ@ R= O}OH sQi l} h} QaD w lQD lwv QO `HQt C=YDNt x=oDUO h} QaD'|OOa p=Ft wO pL =@ w xOW |R=UxO=}B =RHt CqO=at VwQ QO lQD |xr�Ut '|}xQo

"CU= xDiQo Q=Qk |HvUCLY OQwt =RHt CqO=at VwQ

[email protected]@modares.ac.ir

|QU '|]N |a=HDQ= CUmW l}v=mt '=RHt CqO=at VwQ %|O}rm u=oS=w"|Oa)@ wO p�=Ut 'Rt=}r} w

xtOkt "1<=RH= VwQ R= xO=iDU= =@ VvD COW ?} Q[ uOQw; CUOx@ |=Q@ '|rmQw]x@QO [3�1]"OwW|t xO=iDU= �lQD lwv u=tr= =}� sQ=yJl} u}mD |=yu=tr= R= OwOLt|OOa |=yVwQ 'OwOLt <=RH= VwQ QO OwHwt CqmWt R= |[a@ Crax@ 'xt=O=|=ysD} Qwor= Ov}=Qi QO OwHwt |=y|oO}J}B p}rO x@ ,qFt "Ov=xDi=} xaUwD |Qo}OOwHwx@ xDi=}xaUwD OwOLt <=RH= VwQ 'OwOLt <=RH= VwQ QO |r=wDt |Ov@u=tr=|=ylQD lwv h=Q]= QO u=w=Qi |=yu=tr= x@ R=}v p}rOx@ 'u}vJty [6�4]"CU= xOt;xOW xO=O xaUwD |RQt u=tr= VwQ 'q=@ |v=tR w |D=@U=Lt |=yxv} Ry w R} Q Q=}U@u}rw= 'CUmW l}v=mt p�=Ut pL |=Q@ |RQt u=tr= VwQ |xaUwD [8�7]"CU=|tm CkO =@ VvD COW ?} Q[ u; QO xm [9]'CU= xOW x�=Q= 1972 p=U QO Q=@OwHw uWmQD j}kO sQi u; QO xm 'lQD u} Qo `@=D uOQw; CUOx@ =@ "Ow@ xOt; CUOx@Ow@y@ VvD ?}=Q[ pL CkO w CiQ u}@ R= lQD |=yx@r |R=UpOt x@ R=}v 'CW=Ou=tr= VwQ QO xm |}=Hv; R= "CU= xDi=} xaUwD Q=}U@ VwQ u}= RwQt= x@ =D w Ci=}x@ C@Uv |D=@U=Lt |=yxv} Ry ^=Lr R= 'OwW|t |Ov@u=tr= xRwL RQt \ki |RQt\=kv QO =yVvD '|RQt u=tr= VwQ QO "OQ=O |Q}osWJ Vy=m 'OwOLt <=RH= VwQpN=O ?=wH QO |@} QkD '|RQt u=tr= VwQ QO =Q} R 'OvQ=O |}q=@ CkO 'u=O}t pN=OVwQ u}=Q@=v@ "CU= xDUw}B w j}kO u=O}t pN=O QO ?=wH w OwW|tv p=ta= u=O}t'�lQD xr�Ut Ovv=ty� CU= O=} R VvD C=Q}}eD u; QO xm |r�=Ut QO |RQt u=tr=

"CU= ?U=vt Q=}U@

Crax@ =yxR=U QO xOWO=H}= |=y|oDN} QwQi QDW}@ xm OvyO|t u=Wv =yVywSB|oDUw}B=v "CU= VvD RmQtD uOt; OwHwx@ w =yxR=U |xUOvy QO |oDUw}B=v O=H}=w lQD 'M=QwU '|oOWR=@ '|UOvy pmW QO C=Q}}eD CQwYx@ Ov=wD|t xUOvy QOuDiQo Q_vQO '=yxR=U |L=Q] w p}rLD QO syt p�=Ut R= |m} u}=Q@=v@ "OW=@ h=mWp@=k |r}rLD CQwYx@ 'OvQ=O lQD xm |r�=Ut R= |Q=}U@ "CU= xR=U QO lQD OwHw QF=|xr�Ut pL "OvwW pL |OOa |=yVwQ =@ p�=Ut u}= CU= sRq w &OvDU}v pLxO}J}B |OL =D '|Ov@pwtQi ^=Lr R= |OOa |=yVwQ R= |NQ@ R= xO=iDU= =@ lQDOw@y@ w xaUwD 'wQu}= R= "OvQ=O |}q=@ ,=D@Uv |xv} Ry '|D=@U=Lt ^=Lr R= w CU=u} QDsyt "CU= Q}PB=v?=vDH= |Qt= 'p�=Ut p}@k u}= pL |=Q@ |OOa |=yVwQ'Ov=xOW xO=O xaUwD CUmW l}v=mt p�=Ut pL |=Q@ uwvm =D xm |OOa |=yVwQ|=yVwQ 'xDi=}xaUwD OwOLt <=RH= '|RQt u=tr= 'OwOLt <=RH= VwQ %R= OvDQ=@aQw]x@ OwOLt <=RH= VwQ "xOW T=}kt |RQt OwOLt <=RH= VwQ w 'u=tr= uwO@xm G}=Q |=yu=tr= |rw 'CU= xOW xOQ@ Q=mx@ CUmW l}v=mt p�=Ut QO |a}UwOvQ=Ov |@ wN CkO '=yxQiL w =ylQD |m}ORv QO 'OvwQ|t Q=m x@ QwmPt VwQ QOs}y=wNv R=}vOQwt CkO x@ 's}vm O=} R xOwOLt u}= QO =Q =yu=tr= O=OaD Qo = |DL w

"Ov=xOQm P=ND= |Y=N Q}@=OD u=QoWywSB 'xOWu=wva pmWt pL |=Q@ "O}UQpw�Ut xOvU} wv �

"1394 10 2 VQ}PB '1394 9 23 x}LqY= '1394 5 25 Ci=} QO %M} Q=D

"""p�=Ut

QORt=}r}

w|QU

?}=Q[

?}QkD

=RHt CqO=at VwQ |v=@t "2CqO=at?}=Q[T} QD=t =D 'OwW|t xO=iDU= |O}rm Q=R@= 4 R= =RHt CqO=at VwQ QOxDWwv sy R= pkDUt w =RHt CQwYx@ sm =L CqO=at x=oDUO w OwW |Q]k sm =LCW=ov `@=wD "2 &q=@ |x@DQt pmW `@=wD "1 %R= xO=iDU= =@ hOy u}= x@ uO}UQ "OvwW|xrO=at |r=QoDv= sQi O}rwD OvwQ "4 w &R}DQwm - =Wv}rm |Q}op=QoDv= VwQ "3 &hW}@J=RHt Cq=at VwQ R= |}xYqN '=Hv}= QO "CU= xOW =}yt \w@ Qt |xr�Ut Q@ sm =L

%CU= xOW x�=Q= Q} R CQwYx@

sm =L CqO=at G=QNDU= "1"2=OD@= '=RHt CqO=at VwQ QO xr�Ut l} R}i u}vJty w xUOvy |R=UpOt Qw_vt x@C=}YwYN s=tD w ?=NDv= (LCO) |rLt C=YDNt `HQt u=wvax@ x]kv l}|=yRQt \ki 'p=L u}a QO "OvwW|t |@=} RQ= u; x@ C@Uv xr�Ut |m} R}i w |UOvy|Ov@u=tr= xr�Ut |=[i Oa)@ R= QDtm Oa)@ l} =@ |}=yu=tr= R= xO=iDU= =@ xr�UtC=YDNt x=oDUO QO x=wNrO |xUOvy l} C=YNWt '1 pmW j@=]t "OvwW|t?=NDv= x@ xHwD =@ "CU= xOW xO=O u=Wv xOWT=}kt C=YDNt x=oDUO w |rY=|}=yRQt "1 %OvwW|t s}UkD xDUO 2 x@ |rm Qw]x@ xr�Ut |=yRQt '|rLt |=yQwLt"2 w 'OvQ}o|t Q=Qk � |a=aW QwLt |wQ w OvQPo|t LCO R= =yu; O=ODt= xmQO "�OvDUy � |R=wt xm |}=yRQt� OvQPo|tv LCO R= =yu; O=ODt= xm |}=yRQtC=Q}}eD |xOwOLt "OQm |Ov@u=tr= =Q swO wv |=yRQt O}=@ \ki '=RHt CqO=at VwQ

[32]"|Oa)@wO p�=Ut |R=UpOt |xwLv "1 pmW

OwOLt <=RH= VwQ |Q}PBh=]av= '|RQt u=tr= VwQ xm CU= QmPx@ sRq xD@r='OvQ=Ov |R=UxDUUo x@ |R=}v xmv}= x@ xHwD=@ 'u=tr= uwO@ |=yVwQ [10�9]"OQ=Ov =Q'Q}N= |=yp=U QO u}vJty [11]"Ov=xDiQo Q=Qk xHwD OQwt R}v lQD p�=Ut pL |=Q@pL |=Q@ 'OQ=O xm |@U=vt CkO x@ xHwD =@ xOWT=}kt |RQt OwOLt <=RH= VwQ?}mQD =@ QwmPt VwQ "CU= xDiQo Q=Qk xHwD OQwt Q=}U@ CUmW l}v=mt p�=UtVwQ QO "OQ=O |OQi x@ QYLvt |=y|oS} w '|RQt u=tr= w OwOLt <=RH= |=yVwQ|R=UxDUUo xr�Ut RQt \ki |RQt u=tr= x@=Wt xOWT=}kt |RQt OwOLt <=RH=|RQt OwOLt <=RH= VwQ [13�12]"OQ=Ov|U=U= pL x@ |R=}v xm Cw=iD u}= =@ 'OwW|t|Ov@Vt xrtH R= OwOLt <=RH= VwQ QO OwHwt CqmWt R= |Q=}U@ 'xOWT=}ktQO xOWC} wkD XwYNt |=yu=tr= R= xO=iDU= =} w lQD lwv h=Q]= QO R} Q Q=}U@<=RH= VwQ xm CU= xOW xO=O u=Wv |}xar=]t QO "Ovm|t hPL =Q lQD lwv h=Q]=[14]"Ovm x@U=Lt =Q VvD COW ?} Q[ Ov=wD|t |DL=Q x@ xOWT=}kt |RQt OwOLt<=RH= VwQ QO =Q CUmW l}v=mt hrDNt p�=Ut 'u=QoWywSB |NQ@ 'xt=O= QO

[17�15]"Ov=xO=O Q=Qk |@=} RQ= OQwt xOWT=}kt |RQt OwOLtxHwD =@ "CU= syt w |}x}=B s}y=it R= |m} Rt=}r} w |QU 'CUmW l}v=mt QOQ=OlQD \}Lt QO VvD w |}=Hx@=H u=O}t 'lQD lwv QO OwHwt |=y|ov}mD x@�1957� Rt=}r} w \UwD Q=@ u}rw= xm 'OwW|t u=}@ xS} w |y=vDt=v |QU l} CQwYx@VvD COW ?} Q[ x@ \w@ Qt Rt=}r} w |QU R= xrtH u}rw= [18]"CU= xOW O=yvW}B|x@U=Lt w 'J p=QoDv= 'lQD OWQ |x@U=Lt QO VvD COW ?} Q[ [19]"CU=|=Q@ |a}Uw C=ar=]t Qw_vt u}tyx@ w 'OQ=O |v=w=Qi OQ@ Q=m Q=OlQD \}Lt |SQv=\w@ Qt Rt=}r} w |QU R= swO ?} Q[ [20�14w11�2]"CU= xDiQo CQwY u; |x@U=Lt|x@U=Lt w 'lQD |Q=O}=B 'lQD OWQ CyH |x@U=Lt QO T VvD "CU= T VvD x@QO OwHwt |=yC}OwOLt x@ xHwD =@ "OQ=O OQ@ Q=m CUmW l}v=mt QO |otQkJVvD |x@U=Lt QO |QDtm C=ar=]t '|OOa |=yVwQ w |y=oW}=tR; |=yVywSB

[24�21]"CU= xOW s=Hv= VvD COW ?} Q[ x@ C@Uv TswO w pw= |=yxrtH Q@ xwqa xm CU= xO=O u=Wv Q}N= |=yp=U C=ar=]t?}=Q[ OQ@ Q=m u} QDsyt "OvQ=O |}=RUx@ C}ty= R}v u; QDq=@ CqtH 'Rt=}r} w |QUQO [25]"CU= Q=OlQD \}Lt QO xR=Ov= QF= |x@U=Lt x@ \w@ Qt Rt=}r} w |QU QDq=@p=U QO "CU= xOW s=Hv= |QDtm Q=}U@ C=ar=]t 'Rt=}r} w |QU QDq=@?}=Q[|x@U=Lt=@ C}=yv|@ jQw |=Q@ |r}rLD |WwQ 'Rt=}r} w |QU ?}=Q[ |x@U=Lt CyH '2012|QO=v |=yVywSB '|OOa |=yVwQ OQwt QO [26]"CU= xOW O=yvW}B |RmQt lQDVwQ x@ u=wD|t u=}t u}= QO "Ovvm x@U=Lt =Q Rt=}r} w |QU ?}=Q[ xm Ov=xDUv=wDOwOLt <=RH= VwQ w u=tr= uwO@ |@}mQD VwQ [27]'xOWT=}kt |RQt OwOLt <=RH=|@}mQD OwOLt <=RH= VwQ w [30w29]'|@}mQD lQD u=tr= VwQ [28]'xOWT=}kt |RQt"Ov=xOQm x@U=Lt pw= Owt |=Q@ =Q Rt=}r} w |QU pw= |xrtH 5 xm 'OQm xQ=W= [31]'p=Dm =Qiw |H=N \UwD xm 'CU=RHt CqO=at VwQ 'O}OH ,=D@Uv |=yVwQ R= |m}"CU= xDiQo Q=Qk|UQQ@ OQwt|irDNtp�=UtpL |=Q@ w 'CU= xOW O=yvW}B u=Q=mty[32]'p}Uv=DB p�=Ut pL x@ u=wD|t =RHt CqO=at VwQ |=yOQw=DUO u} QDsyt R="OQm xQ=W= [37]'|a=HDQ= G=wt= Q=WDv= w [36w35]'l}t=v}OwDUq= [34w33]'l}D=DU=wDUq=Rt=}r} w |QU uOR?} QkD |=Q@ =RHt CqO=at VwQ |xaUwD 'Q[=L VywSB R= hOyQO &=RHt CqO=at |r}rLDxt}v VwQ R= |}xYqN '2VN@ QO Qw_vt u}= |=Q@ "CU=G=QNDU= |xwLv '4 VN@ QO &Q=OlQD |=y\}Lt QO QwmPt VwQ |xaUwD '3 VN@|HvUQ=@Da= '5 VN@ QO &=RHt CqO=at VwQ R= xO=iDU= =@ Rt=}r} w |QU ?}=Q[Q=DWwv |=yOQw=DUO '6 VN@ QO s=Hv=QU w &p=Ft 2 R= xO=iDU= =@ Q[=L VwQ

"CU= xOW x�=Q= Q[=L

88

2/2|

xQ=tW'

33�2

|xQwO'

�1396

u=DU@=D�

h}QW

u=Qta|

UOvyt %O};|t CUOx@ 9 |x]@=Q R= w CU= p=kDv= |@ wm =S T} QD=t J(�; �) 'u; QO xm

J(�; �) =

"x;�(�; �) y;�(�; �)

x;�(�; �) y;�(�; �)

#(9)

10 |x]@=Q CQwYx@ 4 w 3 \@=wQ R= xO=iDU= =@ R}v =yRQt |wQ |@ wm =S T} QD=t%OwW|t x@U=Lt

J(�) =

"x(�) y(�)

x;�(�) y;�(�)

#(10)

%OwW|th} QaD 11 |x]@=Q CQwYx@ [L] jDWt QwD=QB=T} QD=t '|Oa)@ wO p�=Ut QO

[L] =

"@@x 0 @

@y

0 @@y

@@x

#T(11)

=@ |rLt w |rm C=YDNt x=oDUO wO QO QwmPt Q=OQ@ |=yjDWt |x]@=Q 'u}=Q@=v@%OwW|t u=}@ 12 |x]@=Q R= xO=iDU=

[L] = [b1(�)] @@�

+ [b2(�)]1�@@�

(12)

%Ov};|t CUOx@ 14 w 13 |=yx]@=Q R= [b2(�)] w [b1(�)] u; QO xm

[b1(�)] = 1jJ(�)j

2664 y(�);� 00 �x(�);�

�x(�);� y(�);�

3775 (13)

[b2(�)] = 1jJ(�)j

2664�y(�) 00 x(�)

x(�) �y(�)

3775 (14)

O=ODt= u; QO p=tQv Q=OQ@ uDUv=O x@ R=}v 'O=ODt= Qy QO uWmQD |x@U=Lt Qw_vtx@15 |x]@=Q CQwYx@ xr�Ut |=yRQt |wQ Q@ fng K]U Q@ Owta p=tQv Q=OQ@ "CU=

%OvwW|t h} QaDfng = 1

jr~xj r~x (15)

CQwYx@ u=wD|t =Q � w � CyH wO |=Q@ 15 |x]@=Q '4 w 3 |=yx]@=Q R= xO=iDU= =@%CWwv 17 w 16 |=yx]@=Q

[n�(�)] = 1�����( y;�(�)

�x;�(�)

)�����2664 y;�(�) 0

0 �x;�(�)

�x;�(�) y;�(�)

3775 (16)

[n�(�)] = 1�����(�y(�)

x(�)

)�����2664�y(�) 0

0 x(�)

x(�) �y(�)

3775 (17)

QO xm 'OwW|t xO=iDU= |Y=N |=y|oS} w =@ pmW `@=wD R= =RHt CqO=at VwQ QOxO=iDU= =@ =yRQt |wQ Q@ `@=wD |@=}uwQO "OvwW|t xO=O u=Wv [N ] =@ |rm Cr=L|}xQo \=kv QO "1 %OvQ=O syt |oS} w wO xm 'OQ}o|t s=Hv= QwmPt pmW `@=wD R=

C=Q}}eD 'OwOLt p�=Ut QO "CU= �1 � � � +1 CQwYx@ � |U=tt QwLt"CU= �=yRQt |wQ Q@� l} w �LCO QO� QiY u}@ � |a=aW QwLt

=@ |rm C=YDNt QO xr�Ut |xRwL uwQO |x]kv Qy C=YDNt 'QwmPt VwQ QO=@ R}v xr�Ut |=yRQt R= x]kv Qy C=YDNt xm |r=L QO 'OwW|t XNWt (x; y)

x@ (x; y) |rm C=YDNt R= xr�Ut |xUOvy p=kDv= Qw_vtx@ "OwW|t u}}aD (x; y)

|q=@ |x@DQt |=y|}xrtHOvJ wv R= xm CW=ov `@=wD R= '(�; �) |rLt C=YDNt|=yRQt |wQ x]kv Qy C=YDNt u}=Q@=v@ "OwW|t xO=iDU= 'OvDUy [�(�)] hW}@JOy=wN x@U=Lt p@=k 2 w 1 |=yx]@=Q CQwYx@ CW=ov `@=wD R= xO=iDU= =@ xr�Ut

%Ow@

x(�) =n�+1Xi=1

'i(�)xi (1)

y(�) =n�+1Xi=1

'i(�) yi (2)

O=OaD n� w |rm C=YDNt x=oDUO QO RQt |wQ \=kv C=YDNt y w x '=yu; QO xmx]kv Qy C=YDNt '=RHt CqO=at VwQ QO "OvDUy RQt |wQ |=yu=tr= |}xQo \=kv

%OwW|t x@U=Lt 4 w 3 \@=wQ R= xO=iDU= =@ xr�Ut |xRwL uwQO

x(�; �) = � x(�) = �n�Xi=1

�i(�)xi (3)

y(�; �) = � y(�) = �n�Xi=1

�i(�) yi (4)

hW}@J |=y|}xrtHOvJ R= xO=iDU= =@ '|}xQo n� +1 u=tr= l} |=Q@ CW=ov `@=D%OwW|t u}}aD 5 |x]@=Q CQwYx@

'i(�) = 2n�

n�Xn=0

1c(i�1)cn

Tn(�(i�1))Tn(�) (5)

|=Q@ u}vJty "CU= n |x@DQt R= pw= wv hW}@J |}xrtHOvJ Tn(�) 'u; QO xmQ=Okt n = 0 w n� Q}O=kt |=Q@ w CU= cn = 1 Q=Okt '0 < n < n� Q}O=ktQm}vwQm |=DrO C}Y=N 'xOt; CUOx@ CW=ov `@=wD '?}DQD u}= x@ "CU= cn = 2

%�6 |x]@=Q� OvQ=O =yxQo R= l} Qy QO'�(��) = ��� (6)

=@ xm OvDUy |iW}@J \=kv 'CU= n = 0; :::; n� =yv; QO xm �n |}xQo \=kv%Ov};|t CUOx@ 7 |x]@=Q R= xO=iDU=

�n = � cos(n�n�

); n = 0; :::; n� (7)

|}=y|}xrtHOvJ =yu; R= xO=iDU= =@ xm 'CU= xQo n� +1 |=Q=O u=tr= Qy 'u; QO xmsm =L \@=wQ G=QNDU= Qw_vt x@ "OvwW|t O}rwD xUOvy |@=}uwQO |=Q@ n� |xHQO R=C=YDNt QO u=tr= K]U <RH "CU= R=}v x}=B\@=wQ R= |NQ@ x@ '|rLt C=YDNt QOCQwYx@ |}x]@=Q (d�d�) |rLt C=YDNt QO u=tr= K]U <RH =@ (dx dy) |rm

%OQ=O 8 |x]@=Qd = dx dy =

����J(�; �)���� d� d� = �

����J(�)���� d� d� (8)

89

"""p�=Ut

QORt=}r}

w|QU

?}=Q[

?}QkDx@ =} '|r}rLD |=yVwQ R= xO=iDU= =@ w |wk sQi x@ u=wD|t =Q 27 sm =L

|=Q@ xOWx�=Q= VwQ |=v@t "OQm pL |OOa CQwYx@ w Q=OuRw |xOv=t}k=@ sQis}y=wN =Q 28 |x]@=Q u}=Q@=v@ 'CU= Q=OuRw |=yxOv=t}k=@ VwQ 'p�=Ut pL

%CW=OZw (�ij;j + fi) d = 0 (28)

QO pO=aD |xrO=at 'Q=OuRw |=yxOv=t}k=@ VwQ R= xO=iDU= =@ '28 |x]@=Q pL =@[33]%OwW|t G=QNDU= 29 |x]@=Q CQwYx@ =RHt CqO=at VwQ

�[D0]fu(�)g;�� + [D1]fu(�)g;� + �fF b(�)g = f0g (29)

h} QaD 32 |r= 30 \@=wQ CQwYx@ OwHwt Q=OQ@ w ?}=Q[ |=yT} QD=t 'u; QO xm%OvwW|t

[D0] =Z +1

�1[B1(�)]T [D][B1(�)]jJ(�)jd� (30)

[D1] =Z +1

�1[B2(�)]T [D][B1(�)];�jJ(�)jd� (31)

fF b(�)g =Z +1

�1[N(�)]T ffb(�; �)gjJ(�)jd� (32)

"CU=yxQo QO |tHL |=ywQ}v Q=OQ@ fF bg = [F bx F by ]T 'Q}N= |x]@=Q QO xmVwQ R= xO=iDU= =@ C@=F ?}=Q[ |=yT} QD=t �29 |x]@=Q� p}Uv=Qi}O |xrO=at QO|Q}op=QoDv= VwQ u}= R= xO=iDU= "CU= xOW x@U=Lt T}DQwm - =Wvrm |Q}opQoDv=?}=Q[ |=yT} QD=t O}rwD ?@U xOW|iQat |xS} w CW=ov w pmW `@=wD x=Qty x@

%�34 w 33 \@=wQ� OwW|t |Q]kD0ij = 2�ijwi[B1(�i)]T [D][B1(�i)]jJ(�i)j (33)

D1ij = 2�ijwi[B2(�i)]T [D][B1(�i)];�jJ(�i)j (34)

29 |x]@=Q Q}oQO CqO=at x=oDUO u}=Q@=v@ "CU= Qm}vwQm |=DrO �ij '=yu; QO xm%CWwv i |O=R; xHQO Qy <=R= x@ 35 |x]@=Q CQwYx@ u=wD|t =Q

�D0ii ui;��(�) +D1

ii ui;�(�) + � F bi (�) = 0 (35)

sm =L |xrO=at pL OvwQ "2"2Q=Qk |Q=PoQ=@ CLD xm |}=yxQo |=Q@ \ki 35 |x]@=Q 'QwmPt VwQ CUNv s=o QO|O=R; C=HQO R= l} Qy |=Q@ VvD C=Q}}eD 'swO s=o QO "OwW|t x@U=Lt 'OvQ=O=@ TBU "OwW|t u}}aD 26 |x]@=Q R= xO=iDU= =@ � QwLt O=ODt= QO QmPr=jwiQO xQo Qy =@ \@DQt RmQtDt |rN=O |=ywQ}v |=yxir-wt Q=Okt 'pO=aD \@=wQ R= xO=iDU=|x]@=Q T=U=Q@ LCO |x]kv QO |rN=O VvD u=R}t u}vJty w � QwLt O=ODt=

%OwW|t x@U=Lt 36

f�LCOg =nXi=1f�LCOig (36)

R= LCO QO |rN=O VvD u}= VNB R=@ x@ xHwD =@ =yxQo R= l} Qy syU 'u; QO xm%O};|t CUOx@ 37 |x]@=Q�

�LCOi

= D0ii

nPj=1

D0jj

��LCO

(37)

|=yQwLt x@ C@Uv =yu; pw= jDWt "2 w &OvQ=O Qm}vwQm |=DrO C}Y=N '=yu=tr=%�19 w 18 |=yx]@=Q� CU= QiY Q@=Q@ =yxQo s=tD QO |U=tt |rLt

N�(��) = ��� (18)

N�;�(��) = 0 (19)

|x@DQt R= |}xrtHOvJ l} '|}xQo n� + 1 u=tr= l} |=Q@ |O=yvW}B pmW `@=wD%OQ=O pwyHt QDt=Q=B 2n� xm 'CU= 20 |x]@=Q CQwYx@ 2n� + 1

Ni(�) =2n�+1Xm=0

am�m = Ni(�) = a0 + a1� + a2�2 + a3�3 + :::

+a2n� + 1�2n�+1 (20)

"OvwW|t u}}aD 19 w 18 |=yx]@=Q \}=QW p=ta= =@ '20 |x]@=Q QO C@=F ?}=Q['(�; �) C=YDNt x@ Q_vOQwt |xr�Ut |=[i R= x]kv Qy QO u=mtQ}}eD |=yxir-wt=@ 'OvwW|t h} QaD fu(�; �)g = [ux(�; �) uy(�; �)]T |x]@=Q CQwYx@ xm=@ RQt |wQ |=yu=tr= Q@ `k=w |=yxQo u=mtQ}}eD ?UL Q@ pmW `@=wD R= xO=iDU=

%OvwW|t x@U=Lt 21 |x]@=Q R= xO=iDU=

fu(�; �)g = [N(�)] fu(�)g = [N(�)]�ux(�) uy(�)

�T(21)

xr�Ut |=[i QO (�; �) |x]kv QO VvQm |=yxir-wt '14 w 13 \@=wQ R= xO=iDU= =@%OvwW|t u=}@ 22 |x]@=Q CQwYx@

f"(�; �)g =�"x(x; y) "y(x; y) xy(x; y)

�T =�B1(�)

��u(�)

;�

+1��B2(�)

��u(�)

(22)

%Ov};|t CUOx@ 24 w 23 |=yx]@=Q R= [B2(�)] w [B1(�)] u; QO xm�B1(�)

�=�b1(�)

��N(�)

�(23)�

B2(�)�

=�b2(�)

��N(�)

�;� (24)

C=YDNt x@ x]kv Qy QO VvD |=yxir-wt OQwt QO 'lwy uwv=k R= xO=iDU= =@ u}vJty%OQm xO=iDU= 26 w 25 |=yx]@=Q R= u=wD|t '(�; �)

f�(�; �)g = [D]f"(�; �)gf�(�; �)g = [D]

�[b1(�)][N(�)]fu(�)g;� (25)

+1�

[b2(�)][N(�)];�fu(�)g�

(26)

"CU= Kr=Yt X=wN T} QD=t Qou=}@ [D] '=yu; QO xmu=}@ 27 |x]@=Q CQwYx@ |Oa)@ wO l}D=DU=wDUq= p�=Ut Q@ sm =L pO=aD |xrO=at

%OwW|t�ij;j + fi = 0 (27)

|=yxir-wt R}v fi w CU= |Oa)@ wO VvD QwUv=D <=RH= Qou=}@ �ij 'u; QO xmQO xm CU= QmP x@ sRq "OvDUy xr�Ut |=[i Q@ xOWp=ta= |tHL |=ywQ}v|xrO=at "OvDUy j = x; y w i = x; y 'l}D=DU=wDUq= p�=Ut |Oa)@ wO Cr=L

90

2/2|

xQ=tW'

33�2

|xQwO'

�1396

u=DU@=D�

h}QW

u=Qta|

UOvyt

"|@]k C=YDNt w xOWT=}kt C=YDNt x=oDUO u}@ |x]@=Q "2 pmW

%CWwv u=wD|t u}=Q@=v@r = �r�(�) (41)

%CW=O s}y=wN =Q 42 |x]@=Q 2 pmW j@] 'Qo}O hQ] R=

� = tan�1�yx

�= tan�1

�y(�)x(�)

�(42)

?UL Q@ =RHt CqO=at VwQ QO VvD w |}=Hx@=H 'CU= xOW QmP ,q@k xm Qw]u=tyR= xO=iDU= =@ u}=Q@=v@ "OvwW|t u=}@ (r; �) ?UL Q@ CUmW l}v=mt QO w (�; �)

|@]k C=YDNt x=oDUO QO =Q =yQDt=Q=B u}= u=wD|t |DL=Qx@ 42 w 41 |=yx]@=Q"OQm p}O@D Qo}Om} x@ |rLt w

lQD l} R}i |R=UpOt "2"3l}D=DU=wDUq= x@=WtlQDCUOQwO\=kv QOVvD u=O}t `} RwD 'Q=OlQD |=y\}Lt QOO}OW C=Q}}eD 'lQD lwv QOC}=yv|@VvD OwHwCrax@ lQD lwvl}ORv|L=wv QO wR= xO=iDU= =@ =Q |}=Hx@=H w VvD u=O}t u=wD|t u}=Q@=v@ "OQ=O OwHw VvD u=O}t QOl}v=mt p�=Ut QO |DQ=mO C=YDNt =} w |@]k C=YDNt QO Rt=}r} w |QU \U@xm OvwW|t x@U=Lt |Qw] |QU ?}=Q[ 'Rt=}r} w |QU QO [19w5]"CWwv CUmWl} R}i |xaUwD |=Q@ u}=Q@=v@ "OW=@ u=}@ p@=k |r}rLD CQwYx@ lQD lwv |ov}mDxO=iDU=|tHL |=ywQ}v|}xQo Q=OQ@ R= O}OH sQil} R= '=RHt CqO=at VwQ QO lQDCU= xOW O=yvW}B u=vJ |tHL |=ywQ}v |}xQo Q=OQ@ R= O}OH sQi u}= "CU= xOWsRq "Ovm u=}@ =Q CUmW l}v=mt |xr�Ut |ak=w l} R}i |ov}mD |xHQO Ov=wD@ xmfbi (�) xQo Qy |=Q@ |tHL |=ywQ}v `@=D 'l}D=DU=wDUq= p�=Ut QO xm CU= QmPx@|=Q@ "(fbi (�) = ai� + bi) Ovm|t Q}}eD � = 1 =D � = 0 R= |]N CQwYx@� = 1 w � = 0 QO uWmQD |RQt \}=QW R= bi w ai pwyHt wO uOQw; CUOx@Q@=Q@ lQD lwv QO VvD 'CUmW l}v=mt p�=Ut QO xm |}=Hv; R= "OwW|t xO=iDU=sQi l} '=RHt CqO=at VwQ =@ lQD |xr�Ut l} R}i u=}@ |=Q@ 'CU= C}=yv|@

%CU= xOW x�=Q= 43 |x]@=Q j@=]t fbi (�) R= O}OHfbi (�) = aip

�+ bi�p�

(43)

lwv QO |ov}mD '=ywQ}v Q=OQ@ R= O}OH sQi u}= |x�=Q= =@ '=RHt CqO=at VwQ QO|x]@=Q u=wD|t jwi C=L}[wD x@ xHwD =@ "OwW|t <=[Q= |m} R}i CQwYx@ lQD|ov}mD xm |}=Hv; R= %OQm x}HwD xvwou}= |m} R}i ^=Lr R= =Q |O=yvW}B |[=} Q

|tHL |rN=O |wQ}v uDiQo Q_vQO =@ 'sm =L |xrO=at Qo}O Q=@ '|Oa@ s=o QO|x]@=Q � O=ODt= QO |tHL Q=@ u=wvax@ LCO |rN=O VvD T=U=Q@ xOWx@U=LtO=ODt= QO xQo Qy x@ \w@ Qt u=mtQ}}eD |=yxir-wt w pL |O=R; xHQO Qy <=R= x@ '38

%OwW|t x@U=Lt � QwLtffbi (�)g = [n�]T f�i(�)g (38)

O=ODt= QO xQo Qy |=Q@ u=mtQ}}eD |=yxir-wt `@=D uOWXNWt =@ w '|v=}=B s=o QOu}vJty "OwW|t |@=}uwQO pmW `@=wD R= xO=iDU= =@ '\=kv Q}=U |=Q@ MU=B '� QwLt"OwW|t u}}aD 26 |x]@=Q R= xO=iDU= =@ R}v xr�Ut |xRwL R= x]kv Qy QO VvD u=R}t

CUmW l}v=mt QO =RHt CqO=at VwQ |xaUwD "3Qy |=Q@ =RHt CqO=at VwQ |xaUwD |=Q@ '2 VN@ QO xOWx�=Q= ?r=]t j@]=RHt CqO=at VwQ QO \w@ Qt |xr�Ut l} R}i w xUOvy xm CU= sRq |}xr�UtCUmWl}v=mtp�=Ut |=Q@VN@ Q} R wO 'xt=O= QO Qw_vt u}ty |=Q@ "OwW G=QNDU=

%CU= xOW x�=Q=

lQD |xUOvy |R=UpOt "1"3CqO=at VwQ |v=@t T=U=Q@ lQD |xr�Ut u=}@ |=Q@ xOWs=Hv= C=ar=]t x@ xHwD =@x@ xr�Ut |xUOvy QO lQD u=}@ CyH ZQi 3 R= CUmW l}v=mt s}y=it w =RHt

%CU= xOW xO=iDU= CQwY u}='OwOLt <=RH= |=yVwQ x@=Wt %xr�Ut |xRwL QO lQD |xUOvy |R=UpOt �xr�Ut |xRwL QO lQD |xUOvy 'xOWT=}kt |RQt OwOLt <=RH= w |RQt u=tr=

"OwW|t |R=UpOt |mJwm Q=}U@ |r=N |=[i CQwYx@C=}YwYN|t=tD '=RHtCqO=atVwQ QO %lQDlwv QO `HQt|x]kvuDiQo Q_vQO �xHwD =@ "OwW|t u=}@ `HQt C=YDNt x=oDUO l} QO xr�Ut |m} R}i w |UOvypLt QO lQD lwv 'lQD |xr�Ut h} QaD |=Q@ `HQt x]kv h} QaD \@=w[ x@Q=m u}= s=Hv= p}qO u} QDsyt R= |m} "CU= xOW xDiQo Q_v QO LCO |x]kvR= xO=iDU= =@ lQD |m} R}i C=}YwYN |xty 'CUmW l}v=mt QO xm CU= u;LCO |x]kv =RHt CqO=at VwQ QO xm |}=Hv; R= w 'OvwW|t u=}@ lQD lwvw =RHt CqO=at VwQ u}@ \=@DQ= O=H}= |=Q@ u}=Q@=v@ 'OQ=O =Q Vkv u}ty R}v

"CU= xOW xO=iDU= QwmPt ZQi R= 'CUmW l}v=mt |x} Q_v|}=Hv; R= %|@]k C=YDNt w xOWT=}kt C=YDNt x=oDUO u}@ x]@=Q G=QNDU= �QO xm |DqO=at CU= sRq 'OwW|t u=}@ |@]k C=YDNt QO lQD |xr�Ut xmu=}@ |@]k C=YDNt QO C}=yv QO R}v 'OvwW|t G=QNDU= =RHt CqO=at VwQ

"OvwWxOWT=}kt C=YDNt x=oDUO w |@]k C=YDNt u}@ |x]@=Q '2 pmW x@ xHwD =@

%OwW|t xDWwv 39 |x]@=Q CQwYx@ =RHt CqO=at VwQ QOr2 = x2 + y2 (39)

w �4 w 3 \@=wQ� =RHt CqO=at VwQ QO xOWT=}kt C=YDNt h} QaD x@ xHwD =@%O};|t CUOx@ 40 |x]@=Q '39 |x]@=Q QO =yv; |Q=Po|=H

r2 = x2 + y2 =��x(�)

�2 +��y(�)

�2 = �2�x(�)2 + y(�)2�= �2r2

�(�) (40)

91

"""p�=Ut

QORt=}r}

w|QU

?}=Q[

?}QkD�a2

n

��2 + n

2 � (�1)n�

sin�n2 � 1

��

��n2 � 1

�sin�n2 � 3

����

(46)

�yy(r; �) =1Xn=1

n2 r

(n=2)�1�a1n

��2� n

2 � (�1)n�

cos�n2 � 1

�� +

�n2 � 1

�cos�n2 � 3

���

�a2n

��2� n

2 + (�1)n�

sin�n2 � 1

��

+�n2 � 1

�sin�n2 � 3

����

(47)

�xy(r; �) =1Xn=1

n2 r

(n=2)�1�a1n

���n2 + (�1)n

�sin�n2 � 1

�� +

�n2 � 1

�sin�n2 � 3

���

�a2n

��n2 � (�1)n

�cos�n2 � 1

��

��n2 � 1

�cos�n2 � 3

����

(48)

ux(r; �) =1Xn=1

rn=2

2G�a1n

���+ n

2 + (�1)n�

cos n2 � �n2 cos�

n2 � 2

���� a2

n

���+ n

2 � (�1)n�

sin n2 �

�n2 sin�n2 � 2

����

(49)

uy(r; �) =1Xn=1

rn=2

2G�a1n

���� n

2 � (�1)n�

sin n2 �

+n2 sin

�n2 � 2

���

+ a2n

���� n

2 + (�1)n�

cos n2 � + n2 cos

�n2 � 2

����

(50)

|=Q@ � = (3 � �)=(1 + �) 'uwU=wB ?} Q[ � '|WQ@ pwOt G '=yu; QO xmu}vJty "CU= K]Ut VvQm p�=Ut |=Q@ � = (3� 4�) w K]Ut VvD p�=Uta2n w CUmW pw= Owt w uQ=kDt xS} w `@=D |=Q@ Rt=}r} w |QU ?}=Q[ a1

n ?} Q[T=U=Q@ w 'OvDUy CUmW swO Owt w uQ=kDt=v |xS} w `@=D |=Q@ Rt=}r} w |QU ?}=Q[Rt=}r} w |QU '4 pmW j@=]t "OvwW|t x@U=Lt xr�Ut |Q=PoQ=@ w |RQt \}=QWOwOLt=v jQw |=Q@ Rt=}r} w |QU |}=Qosy `=aW w OwW|t xO=O \U@ lQD lwv pwL"CU= �lQD pw] |xR=Ov= x@� r = 2a Q@=Q@ CU=Q CtU lwv QO |RmQt lQD =@"OQ=Ov =Q xOwOLt u}= R= QD=Qi QO |}=Qosy |}=v=wD Qw@ Rt |QU xm CU= QmP x@ sRqw |}=Hx@=H MU=B 'CU= |r}rLDxt}v VwQ l} =RHt CqO=at VwQ xm |}=Hv; R=|OOa CQwYx@ =yRQt QO MU=B w |r}rLD CQwYx@ u=O}t R= x]kv Qy QO VvD'`HQt |x]kv R= xQo Qy |=DU=Q QO '3 pmW j@=]t u}=Q@=v@ "Ot; Oy=wN CUOx@|r}rLD CQwYx@ 26 |x]@=Q j@=]t VvD w 45 |x]@=Q j@=]t |}=Hx@=H MU=Bpw= Owt |=Q@ Rt=}r} w |QU ?}=Q[ |x@U=Lt |=Q@ "O};|t CUOx@ � R= |a@=D w45 |x]@=Q j@=]t |}=Hx@=H u=O}t '42 w 41 |=yx]@=Q R= xO=iDU= =@ CUmW|QU TBU w p}O@D |@]k C=YDNt x=oDUO x@ xOWT=}kt C=YDNt x=oDUO R=

"OwW|t xOR ?} QkD pw= |xrtH 4 R= xO=iDU= =@ Rt=}r} w

"=RHt CqO=at VwQ QO lQD |xr�Ut u=}@ "3 pmW

u}=Q@=v@ 'OwW|t u=}@ r(n2 �1) ?} Q[ =@ Rt=}r} w |QU QO lQD lwv QO VvDxOW O=yvW}B R}v |tHL |=ywQ}v |}xQo Q=OQ@ '|ov}mD sQi u}ty R= xO=iDU= =@\ki |tHL |}xQo |=ywQ}v Q=OQ@ '=RHt CqO=at VwQ QO uwJ |iQ] R= "CU=xOW xO=iDU= |O=yvW}B `@=D QO xrtH wO R= \ki 'OW=@ xDW=O pwyHt wO Ov=wD|t

"CU=R= bi w ai ?}=Q[ �3 pmW� xr�Ut |RQt \}=QW p=ta= =@ 'Qo}O p�=Ut x@=Wtu=wD|t p=L "OwW|t x@U=Lt � = 1 w � = 0 QO (Tb; T ) uWm =QD |RQt \}=QWj@=]t =Q =RHt CqO=at VwQ QO xQo Qy |}=yv |}=Hx@=H |=Q@ p}Uv=Qi}O |xrO=at

%OQm |U} wvR=@ 44 |x]@=Q j@=]t Q=OlQD |=yu=O}t |=Q@ 35 |x]@=Q

� D0ii ui;��(�) +D1

iiui;�(�) + aip� + bip

�= 0 (44)

i |O=R; xHQO Qy x@ \w@ Qt MU=B '�44 |x]@=Q� Q}N= p}Uv=Qi}O |xrO=at pL =@%Ot; Oy=wN CUOx@ 45 |x]@=Q j@=]t

ui(�) = Ai�

�D0ii�D1

iiD0ii

�+ Bi

(D0ii �D1

ii)� 4ai

(D0ii � 2D1

ii)p�

� 4bi3(D0

ii + 2D1ii)�p� (45)

QO |Q=PoQ=@ w lQD lwv QO |}=Hx@=H� xr�Ut |RQt \}=QW R= Bi w Ai 'u; QO xm"OvwW|t x@U=Lt |O=R; xHQO Qy |=Q@ �=yxQo

Rt=}r}w |QU ?}=Q[ G=QNDU= "4Q}O=kt w `@=wD O=yvW}B w xS} w \U@ VwQ R= xO=iDU= =@ Rt=}r} w '1957 p=U QO|r=wL QO u}mD VvD w |}=Hx@=H u=O}t '|Q}*= VvD `@=D QO l}vwtQ=y|=@ |xS} w|Oa)@ wO |a=HDQ= |=y\}Lt |=Q@ u}=Q@=v@ "CU= xO=O O=yvW}B =Q Q=OlQD |=y\}Lt|r= 46 |=yx]@=Q CQwYx@ Rt=}r} w |QU j@=]t |}=Hx@=H w VvD u=O}t 'Q=OlQD

[18]%CU= x@U=Lt p@=k 50

�xx(r; �) =1Xn=1

n2 r

(n=2)�1�a1n

��2 + n

2 + (�1)n�

cos�n2 � 1

�� �

�n2 � 1

�cos�n2 � 3

���

92

2/2|

xQ=tW'

33�2

|xQwO'

�1396

u=DU@=D�

h}QW

u=Qta|

UOvyt

"pw= p=Ft QO |UQQ@ OQwt |Oa)@wO |xr�Ut "5 pmW

CyH `H=Qt Q}=U QO xOWx@U=Lt Rt=}r}w |QU pw= ?} Q[ 5 |x@U=Lt "1 pwOH"pw= p=Ft |HvUQ=@Da=

`HQt Rt=}r}w |QU ?}=Q[[29] [22] [22]

0 7665 0 7670 0 7680 a11�0 2779 �0 2760 �0 2777 a120 1915 0 1880 0 1866 a13�0 0018 �0 0033 �0 0030 a14�0 0235 �0 0320 �0 0279 a15

"C}=yv|@ jQw QO CU=Q CtU lQD lwv QO Rt=}r}w |QU |}=Qosy `=aW "4 pmW

=RHt CqO=at VwQ |HvUQ=@Da= "5VwQ R= xO=iDU= =@ VvD u=O}t w |}=Hx@=H u=O}t 'Rt=}r} w |QU ?} QkD CyH"CU= xOW xO=iDU= MATLAB Q=Ri=sQv |U} wvxt=vQ@ \}Lt R= =RHt CqO=atu; G}=Dv w ?=NDv= Q[=L VwQ |HvUQ=@Da= CyH uwtR; p=Ft 2 'Qw_vt u}tyx@

"CU= xOW xU}=kt OwHwt `H=Qt Q}=U =@

|RmQt lQD =@ OwOLt jQw "1"5j@=]t |RmQt lQD =@ OL=w b � 2h O=a@= =@ OwOLt |kQw x@ \w@ Qt 'pw= p=FtpwOt '� = 1 Q@=Q@ u; Q=Okt w s�=k CyH QO |Q=PoQ=@ "CU= �hr=5� pmW=@ "CU= a=b = 0 25 C@Uv w '� = 0 25 uwU=wB ?} Q[ 'E = 1 |a=HDQ==@ jQw sQ=yJ l} \ki uQ=kD Crax@ 'p=Ft u}= QO K]Ut VvD Q=DiQ uDiQo Q_vQOxOW |R=UpOt �G w ?5� |=ypmW j@=]t �|O=R; xHQO 42� |}xQoxU u=tr= 10j@=]t [29w22]'`H=Qt |OOa G}=Dv R= =RHt CqO=at VwQ |HvUCLY |=Q@ "CU=w x@U=Lt |}=Hx@=H u=O}t '45 |x]@=Q R= xO=iDU= =@ "CU= xOW xO=iDU= 1 pwOH"CU= xOW x@U=Lt 2 pwOH j@=]t Rt=}r} w |QU pw= |xrtH 4 ?}=Q[ 'xt=O= QO51 |x]@=Q R= 2 pwOH QO xOWx�=Q= |=]N |x@U=Lt |=Q@ xm CU= QmP x@ sRq

%CU= xOW xO=iDU=Error(�) = aave � aDEM

aave� 100 (51)

VwQ QO xOWx@U=Lt Q=Okt aDem w Qo}O `H=Qt Q}O=kt u}ov=}t aave 'u; QO xm|=yu=O}t 'xOW xOR?} QkD Rt=}r} w |QU R= xO=iDU= =@ TBU "CU=RHt CqO=atu=v}t]= pwYL CyH "CU= xOW x�=Q= 7 w 6 |=ypmW j@=]t |}=Hx@=H w VvDQO VvD w |}=Hx@=H C=Q}}eD OvwQ 9 w 8 |=ypmW QO xOt; CUOx@ |=yQwDv=m R=|m} 'Rt=}r} w |QU ?}=Q[ Q@ xwqa "CU= xOW xU}=kt =yVywSB Q}=U =@ y |=DU=QCU= Q=OlQD |=y\}Lt QO |}=Hx@=H wVvD u=O}t?} QkD 'Q[=L Q=DWwv h=Oy= R=uO=O u=Wv '7 w 6 |=ypmW R= hOy u}=Q@=v@ "�CUOQwO \=kv QO w lQD l}ORv QO�9 w 8 |=ypmW \UwD u; CkO xm 'CU= Q=OlQD \}Lt |}=Hx@=H w VvD u=O}tQO VvD COW xm OwW|t xOy=Wt 8 pmW QO "CU= xDiQo Q=Qk |HvUQ=@Da= OQwtG}=Dv x@=Wt w xOW sm VvD Q=Okt QDQwO |=yu=O}t QO w xOw@ O=} R Q=}U@ lQD h=Q]=

"CU= xOW l}D=DU=wDUq=G}=Dv Q}=U =@ 2 pwOH QO =RHt CqO=at VwQ R= pY=L G}=Dv |xU}=kt =@u}ov=}t u}@ hqDN= T=U=Q@� xOWO=H}= |=]N |x@U=Lt w 1 pwOH QO |OOa

93

"""p�=Ut

QORt=}r}

w|QU

?}=Q[

?}QkD

QO |Q=PoQ=@ pLt =D lQD lwv R= jQw |rw] |=DU=Q QO uy |}=Hx@=H "9 pmW"pw= p=Ft

"xRwL pm QO VvD u=O}t QO Rt=}r}w |QU CqtH O=OaD QF= "10 pmW

CqO=at VwQ '9 w 8 |=ypmW j@=]t "OQ=O xrtH 4 \ki 'Ovm|t pL |r}rLDu}=Q@=v@ "CU= xOR ?} QkD Q=OlQD \}Lt QO =Q VvD w |}=Hx@=H |DUQOx@ =RHt|QU QO OwHwt |xrtH C}=yv|@ '=RHt CqO=at VwQ xm CiQo xH}Dv u=wD|tQO Rt=}r} w |QU QDq=@ CqtH C=QF= w xOR ?} QkD |DUQOx@ xrtH 4 QO =Q Rt=}r} w10 pmW R= xm Qw]u=ty "CU= xOW xOv=Hvo =RHt CqO=at VwQ swO |xrtH'OwW ?=NDv= Rt=}r} w |QU ?}=Q[ R= |tm O=OaD xm |DQwY QO 'CU= XNWtQDQwO \=kv QO w OwW|t u=}@ |DUQO x@ lQD lwv |o}=Uty QO \ki VvD u=O}tCqO=at VwQ QO xm CU= |Qw;O=} x@ sRq "CU= QwO C}ak=w R= G}=Dv lQD R=CUOx@ CUmW l}v=mt |RQt \}=QW R= swU w pw= |=yxrtH xm |}=Hv; R= '=RHtswO |=y?} Q[ "OvDUy ?U=vt CkO =@ lQD |m}ORv QO =yMU=B u}=Q@=v@ 'Ov=xOt;|yHwD p@=k |=]N xm 'Ov=xOt; CUOx@ l}D=DU=wDUq= |RQt \}=QW R= sQ=yJ wVwQ xmv}= x@ xHwD =@ xm OQm \=@vDU= u}vJ 10 pmW j@=]t u=wD|t u}=Q@=v@ "OvQ=O?}=Q[ O}=@ 'xRwL pm R= j}kO G}=Dv pwYL |=Q@ 'OQ=O xrtH 4 \ki =RHt CqO=atxOv=Hvo QDu}}=B |x@DQt ?}=Q[ QO |wLvx@ lQD CUOQwO |L=wv QF= w QDq=@ |x@DQt'xOWQmP ?r=]t j@=]t u}=Q@=v@ "CU= xO=Di= j=iD= swO |xrtH QO Qt= u}= "OvwWu=}@ =Q lQD R= QwO j]=vt w QDq=@ CqtH C=QF= 'Rt=}r} w |QU swO |xrtH ?} Q[xOW Q_vhQY swO |xrtH |=Q@ =]N |x�=Q= R= 2 pwOH QO p}rO u}ty x@ "Ovm|t

"CU=

|}x@r lQD =@ OwOLt jQw "2"5j@=]t |}x@r lQD =@ OL=w b � 2h O=a@= =@ OwOLt |kQw x@ \w@ Qt swO p=Ft

VwQ R= xO=iDU= =@ xOWx@U=Lt Rt=}r}w |QU pw= ?} Q[ 4 |x@U=Lt "2 pwOH"pw= p=Ft QO =RHt CqO=at

��� =]N u}ov=}t Q[=L VwQ Rt=}r}w |QU ?}=Q[1 42 0 7812 a11

� 0 2280 a123 77 0 1901 a13

0 0 a14

"pw= p=Ft QO =RHt CqO=at VwQ QO xOWx@U=Lt �yy VvD u=O}t "6 pmW

"pw= p=Ft QO =RHt CqO=at VwQ QO xOWx@U=Lt uy |}=Hx@=H u=O}t "7 pmW

"pw= p=Ft QO |Q=PoQ=@ pLt =D lQD lwv R= jQw |rw] |=DU=Q QO �yy "8 pmW

CkO Q[=L VwQ xm OwW|t xOy=Wt '�=RHt CqO=at VwQ =@ `H=Qt Q}=U G}=DvQ_vQO CQwY QO "OQ=O sQ=yJ w 'swU 'pw= ?}=Q[ |x@U=Lt |=Q@ |@U=vt Q=}U@|}q=@ ,=D@Uv |=]N 'Rt=}r} w |QU swO ?} Q[ 'Rt=}r} w |QU ?}=Q[ Ctqa uDiQou=O}t CU= XNWt xm Qw]u=ty %OQm x}HwD u}vJ u=wD|t =Q u; Cra xm 'OQ=O�OQ=O xrtHC}=yv|@ xm � Rt=}r} w |QU T=U=Q@ Q=OlQD |=y\}LtVvD w |}=Hx@=HCQwYx@ =Q xr�Ut |xvt=O xm '=RHt CqO=at VwQ xm |r=L QO 'OwW|t x@U=Lt

94

2/2|

xQ=tW'

33�2

|xQwO'

�1396

u=DU@=D�

h}QW

u=Qta|

UOvyt CyH `H=Qt Q}=U QO xOWx@U=Lt Rt=}r}w |QU pw= ?} Q[ 5 |x@U=Lt "3 pwOH

"swO p=Ft |HvUQ=@Da=`HQt Rt=}r} w |QU ?}=Q[

[27] [30] [31]

1 0596 1 0597 1 0585 a11�0 1506 �0 1503 �0 1513 a120 0794 0 0820 0 08150 a13�0 0476 �0 0486 �0 0478 a140 0056 0 0063 0 0055 a15

VwQ R= xO=iDU= =@ xOWx@U=Lt Rt=}r} w |QU pw= ?} Q[ 4 |x@U=Lt "4 pwOH"swO p=Ft QO =RHt CqO=at

��� =]N u}ov=}t Q[=L VwQ Rt=}r} w |QU ?}=Q[0 16 1 0973 a11

� 0 2951 a126 37 0 0781 a13

0 0 a14

"swO p=Ft QO =RHt CqO=at VwQ QO xOWx@U=Lt �yy VvD u=O}t "12 pmW

"swO p=Ft QO =RHt CqO=at VwQ QO xOWx@U=Lt uy |}=Hx@=H u=O}t "13 pmW

=RHt CqO=at VwQ R= pY=L G}=Dv |xU}=kt =@ "CU= xOW xU}=kt =yVywSB'pw= p=Ft Ovv=ty R}v swO p=Ft QO xm OwW|t xOy=Wt |OOa G}=Dv Q}=U =@p=Ft Ovv=ty =Q u; Cra u=wD|t xm 'OQ=O =]N |tm Rt=}r} w |QU swO ?} Q[OQwt QO Rt=}r} w |QU sQ=yJ ?} Q[ xm CU= QmP x@ sRq xD@r= "OQm x}HwD pw=hqDN= sy =@ �50 OwOL QO R}v [29w22]'u=QoWywSB Q}=U C=ar=]t QO swO p=Ft

"OvQ=O

"swO p=Ft QO |UQQ@ OQwt |Oa)@ wO |xr�Ut "11 pmW

p=Ft x@=Wt 'swO p=Ft QO Kr=Yt C=YNWt w |Q=PoQ=@ "CU= �hr=11� pmWG}=Dv R= =RHt CqO=at VwQ |HvUCLY |=Q@ "CU= xOW xDiQo Q_vQO pw=Q_vQO =@ "CU= xOW xO=iDU= 3 pwOH j@=]t [31w30'27]=yVywSB Q}=U |OOa13 =@ jQw hYv \ki uQ=kD Crax@ 'p=Ft u}= QO K]Ut VvD Q=DiQ uDiQo|R=UpOt �G w ?� 11 pmW j@=]t |O=R; xHQO 26 =@ wtHt QO w xQo4 ?}=Q[ 'xt=O= QO w |}=Hx@=H u=O}t '45 |x]@=Q R= xO=iDU= =@ "CU= xOWxO=iDU= =@ TBU "CU= xOW x@U=Lt 4 pwOH j@=]t Rt=}r} w |QU pw= |xrtH12 |=ypmW j@=]t VvD w |}=Hx@=H u=O}t 'xOW xOR?} QkD Rt=}r} w |QU R=QO xOt; CUOx@ |=yQwDv=m R= u=v}t]= pwYL CyH "CU= xOW x�=Q= 13 wQ}=U =@ y |=DU=Q QO VvD w |}=Hx@=H C=Q}}eD OvwQ '15 w 14 |=ypmW

95

"""p�=Ut

QORt=}r}

w|QU

?}=Q[

?}QkD

|Q}oxH}Dv "6?} Q[ |x@U=Lt QO OQ@ Q=m Cra x@ CUmW l}v=mt |x} Q_v QO Rt=}r} w |QU|xrtH ?} Q[� T � VvD ?} Q[ '�Rt=}r} w |QU pw= |xrtH ?} Q[� VvD COW|x@U=Lt w �Rt=}r} w |QU QO q=@ |x@DQt?}=Q[� lQD |xR=Ov= QF= '�Rt=}r} w |QU swO,=Q}N= "OQ=O |v=w=Qi C}ty= 'Q=OlQD |=y\}Lt uwt=Q}B QO VvD w |}=Hx@=H u=O}txDUw}B |=y\}LthrDNt p�=Ut pL |=Q@ =RHt CqO=at VwQ s=v x@ |O}OH VwQ

"CU= xOW xO=O xaUwD|xUOvy |R=UpOt 'lQD lwv QO `HQt |x]kv uDiQoQ=Qk =@ 'Q[=L Q=DWwv QOVwQ QO CUmW l}v=mt |xr�Ut '|tHL |}xQo |=ywQ}v R= |O}OH sQi w lQDp�=Ut QO Rt=}r} w |QU pw= |xrtH 4 'xt=O= QO "CU= xOW h} QaD =RHt CqO=at2 pL =@ 'C}=yv QO "CU= xOW xO=O xaUwD CUmW pw= |=yOwt |=Q@ |Oa)@ wO|m =L G}=Dv "CU= xDiQo Q=Qk |HvUCLY OQwt |O=yvW}B VwQ '|OOa p=FtRt=}r} w |y=vDt=v |QU |DUQOx@ =RHt CqO=at VwQ QO xrtH 4 xm CU= u; R=|x@U=Lt |=Q@ |@U=vt CkO '=RHt CqO=at VwQ u}vJty "CU= xOR ?} QkD =Q?} Q[ xm CU= |r=L QO u}= w OQ=O Rt=}r} w |QU R= sQ=yJ w 'swU 'pw= |=yxrtHQDq=@ CqtH `}tHD u=wD|t =Q u; Cra xm 'OQ=O |O=} R |=]N 'Rt=}r} w |QU swOVwQ QO lQD CUOQwO |L=wv QF= u}vJty w 'swO |xrtH QO Rt=}r} w |y=vDt=v |QU

"CUv=O =RHt CqO=at

"swO p=Ft QO |Q=PoQ=@ pLt =D lQD lwv R= jQw |rw] |=DU=Q QO �yy "14 pmW

QO |Q=PoQ=@ pLt =D lQD lwv R= jQw |rw] |=DU=Q QO uy |}=Hx@=H "15 pmW"swO p=Ft

�References� `@=vt1. Henshell, R.D. and Shaw, K.G. \Crack tip �nite ele-

ments are unnecessary", International Journal for Nu-merical Methods in Engineering, 9(3), pp. 495-507(1975).

2. Barsoum, R.S. \On the use of isoparametric �nite ele-ments in linear fracture mechanics", International Jour-nal for Numerical Methods in Engineering, 10(1), pp.25-37 (1976).

3. Liebowitz, H. and Moyer Jr, E.T. \Finite element meth-ods in fracture mechanics", Computers and Structures,31(1), pp. 1-9 (1989).

4. Belytschko, T. and Black, T. \Elastic crack growth in�nite elements with minimal remeshing", InternationalJournal for Numerical Methods in Engineering, 45(5),pp. 601-620 (1999).

5. Mohammadi, S., Extended Finite Element Method: forFracture Analysis of Sructures, John Wiley & Sons, 280p. (2008).

6. Ghasemi Ghalebahman, A. and Salavati, S. \Utilizingthe extended fnite element method for determining crackstress intensity factors and higher order terms coef-�cients", Modares Mechanical Engineering, 15(2), pp.135-146 (In Persian) (2015).

7. Aliabadi, M.H., The Boundary Element Method, Volume2: Applications in Solids and Structures, Wiley, 598 p.(2002).

8. Cruse, T.A. and Vanburen, W. \Three-dimensional elas-tic stress analysis of a fracture specimen with an edge

crack", International Journal of Fracture Mechanics,7(1), pp. 1-15 (1971).

9. Cruse, T.A. \Numerical evaluation of elastic stress inten-sity factors by the boundary-integral equation method,Surface Crack", Physical Problems and ComputationalSolutions, ASME, New York, pp. 153-170 (1972).

10. Portela, A., Aliabadi, M.H. and Rooke, D.P. \E�cientboundary element analysis of sharp notched plates", In-ternational Journal for Numerical Methods in Engineer-ing, 32(3), pp. 445-470 (1991).

11. Belytschko, T., Lu, Y.Y. and Gu, L. \Element-freeGalerkin methods", International Journal for Numeri-cal Methods in Engineering, 37(2), pp. 229-256 (1994).

12. Song, C. and Wolf, J.P. \Scaled boundary �nite-elementmethod - a primer: Derivations", Computers and Struc-tures, 78(1), pp. 191-210 (2000).

13. Wolf, J.P. and Song, C. \Scaled boundary �nite-elementmethod - a primer: Solution procedures", Computers andStructures, 78(1), pp. 211-225 (2000).

14. Song, C. and Wolf, J.P. \Semi-analytical representationof stress singularities as occurring in cracks in anisotropicmulti-materials with the scaled boundary �nite-elementmethod", Computers and Structures, 80(2), pp. 183-197(2002).

15. Yang, Z. \Fully automatic modelling of mixed-modecrack propagation using scaled boundary �nite elementmethod", Engineering Fracture Mechanics, 73(12), pp.1711-1731 (2006).

16. Yang, Z. \Application of scaled boundary �nite elementmethod in static and dynamic fracture problems", ActaMechanica Sinica, 22(3), pp. 243-256 (2006).

96

2/2|

xQ=tW'

33�2

|xQwO'

�1396

u=DU@=D�

h}QW

u=Qta|

UOvyt

17. Song, C. and Vrcelj, Z. \Evaluation of dynamic stressintensity factors and T-stress using the scaled boundary�nite-element method", Engineering Fracture Mechan-ics, 75(8), pp. 1960-1980 (2008).

18. Williams, M.L. \On the stress distribution at the baseof a stationary crack", Journal of Applied Mechanics,24(1), pp. 109-114 (1957).

19. Anderson, T.L., Fracture Mechanics: Fundamental andApplications, Taylor & Francis (2005).

20. Bird, G.E., Trevelyan, J. and Augarde, C.E. \A cou-pled BEM/scaled boundary FEM formulation for accu-rate computations in linear elastic fracture mechanics",Engineering Analysis with Boundary Elements, 34(6),pp. 599-610 (2010).

21. Sham, T.L. \The determination of the elastic T-term us-ing higher order weight functions", International Journalof Fracture, 48(2), pp. 81-102 (1991).

22. Fett, T. \T-stresses in rectangular plates and circulardisks", Engineering Fracture Mechanics, 60(5), pp. 631-652 (1998).

23. Chen, C.S., Krause, R., Pettit, R.G., Banks-Sills, L. andIngra�ea, A.R. \Numerical assessment of T-stress com-putation using a p-version �nite element method", Inter-national Journal of Fracture, 107(2), pp. 177-199 (2001).

24. Song, C. \Evaluation of power-logarithmic singularities,T-stresses and higher order terms of in-plane singu-lar stress �elds at cracks and multi-material corners",Engineering Fracture Mechanics, 72(10), pp. 1498-1530(2005).

25. Karihaloo, B.L. \Size e�ect in shallow and deep notchedquasi-brittle structures", International Journal of Frac-ture, 95(1), pp. 379-390 (1999).

26. Hello, G., Ben Tahar, M. and Roelandt, J.M. \Analyt-ical determination of coe�cients in crack-tip stress ex-pansions for a �nite crack in an in�nite plane medium",International Journal of Solids and Structures, 49(3),pp. 556-566 (2012).

27. Chidgzey, S.R. and Deeks, A.J. \Determination of co-e�cients of crack tip asymptotic �elds using the scaledboundary �nite element method", Engineering FractureMechanics, 72(13), pp. 2019-2036 (2005).

28. He, Y.Q., Yang, H.T. and Deeks, A.J. \Determination ofcoe�cients of crack tip asymptotic �elds by an element-free Galerkin scaled boundary method", Fatigue and

Fracture of Engineering Materials and Structures, 35(8),pp. 767-785 (2012).

29. Xiao, Q.Z., Karihaloo, B.L. and Liu, X.Y. \Direct deter-mination of SIF and higher order terms of mixed modecracks by a hybrid crack element", International Journalof Fracture, 125(3-4), pp. 207-225 (2004).

30. Karihaloo, B.L. and Xiao, Q.Z. \Accurate determinationof the coe�cients of elastic crack tip asymptotic �eld bya hybrid crack element with p-adaptivity", EngineeringFracture Mechanics, 68(15), pp. 1609-1630 (2001).

31. Su, R.K.L. and Fok, S.L. \Determination of coe�cientsof the crack tip asymptotic �eld by fractal hybrid �niteelements", Engineering Fracture Mechanics, 74(10), pp.1649-1664 (2007).

32. Khaji, N. and Khodakarami, M.I. \A new semi-analytical method with diagonal coe�cient matrices forpotential problems", Engineering Analysis with Bound-ary Elements, 35(6), pp. 845-854 (2011).

33. Khodakarami, M. and Khaji, N. \Analysis of elastostaticproblems using a semi-analytical method with diagonalcoe�cient matrices", Engineering Analysis with Bound-ary Elements, 35(12), pp. 1288-1296 (2011).

34. Khaji, N. and Khodakarami, M.I. \A semi-analyticalmethod with a system of decoupled ordinary di�erentialequations for three-dimensional elastostatic problems",International Journal of Solids and Structures, 49(18),pp. 2528-2546 (2012).

35. Khodakarami, M.I., Khaji, N. and Ahmadi, M.T. \Mod-eling transient elastodynamic problems using a novelsemi-analytical method yielding decoupled partial dif-ferential equations", Computer Methods in AppliedMechanics and Engineering, 213-216(0), pp. 183-195(2012).

36. Mirzajani, M. and Khaji, N. \Decoupled equationsmethod for solving two-dimensional elastodynamic prob-lems in frequency domain", Sharif Journal of Civil En-gineering, 30(3), pp. 65-74 (In Persian)(2014).

37. Khodakarami, M.I. and Khaji, N. \Wave propagationin semi-in�nite media with topographical irregularitiesusing Decoupled Equations Method", Soil Dynamics andEarthquake Engineering, 60, pp. 102-112 (2014).

97