地震引致儲液槽內的流體動力分析---假頻譜矩陣元素數值模式建立與平行化(I)

全文

(1)

(2) !"$# %& ' () *,+ -,.0/ 1 2 3 4 4 5 6 798 :,; <=$> ?@ A0B CD E$F G 4 4 H HHHHHHHHHHHHHHHHHHHHHHHH IKJMLON P Q RSNSTUIKJ VXWUYSTUIKJ IKJMZ\[ P ]_â`cbed f ghgjikg f l f dmihi f dminb f oUprqtsu vew x wzy { wm| }~ vew x wh { | } K ur\. SrMUX¡ ¢¤£ u ¥\¦§M¨K©cªX«c¬®°¯¤± rr²¤³ ¥\¦X´XµM¶K·©cªX«c¬®°¯¤± rr²¤³ ¥¸©º¹§¡»O¼¤½¿¾ ÀÁ¤ÂrÃrÄUÅÇÆXÈ\ÉÊÌËÍÎ¤Ï ÐÑ¡ÒÔÓÇÕÇÖØ×cÙKÚÑMÛÖØ×cÃrÄÜÇÎ¤Ï ÝUÞÌß¡à áãâMäUåærçè¡é\êXëÇìí î ï ð ñ òôó õ ö ÷ øúù û.

(3)

(4) !#" $&%&')(&*,+ - . / 01325476&8#9:<;>=7?,@BA>CDEF GIHKJMLKNPORQTSVUXWYLKZ\[]SVU ^`_TacbKadbKegfihjbKkVlMbKmgnjfoeVpPqT_YrMstbKavuYbKw\x<fzyTs{a}|~{ fXf ^^Kst^Tas{lVmVoas{wtf ^f b]j^Ksus{eMomVst^ { ¢¡{£ ¤ ¥{¥§¦¨¥ ¤ © ¤ £z¦{¦ ¤ £z¦ª¡ «¬¯®±° ¡{£ ² ¥ ³ ¦ íµ ¡{£ ² ¶ ³ ·§¦ ´ ¸¹»º I¼¯½¾ ¿ÁÀÃÂÅÄ «¬ÆÇÉÈ¢ÊÌË¯ÍÃÎ ÏÐÑÅÒ¯ÓÅÔÌÆÕ#Ö ×ØÚÙÛ º,Ü VÝÅÞ&ß W=0/(%!!- &6- 50O1910(+,- $7"!')> #1')(035C'J& 1'J0- 0!8 àáãâIä}å\æ +,5+,01& $%(0:?C(%* (&6%@- 0!8@& 27011+,')- $%(%*!')>#* &SMH ç}èjéjêìëîíðï}ñRòôó}õðöî÷ìøðùûúðüýÿþ Q2M9:AB- * *!?7C'=9:%;6*;=5'X10!'#&6(0!- 0!83;6')(%$& 1')!;

(5) ! & 27&6(01G H. " #%$'&)(*+,.-0/12*34&657 8*3 9 :; YXZ1[!\^]<_Fà b ;6')%T#1');=(%$%>O1%('J& 2E5(G>O"!#"%$& ')(%* <%= >%?A@BCD? &FEGH2*I9 J KLM ?ON +,& 2! @ 7 8 D? 9P+%Q BC?ARST &UEG2PVW( *+ BC2T0X &ZY2[(*+\] I >%?^N_@ H S cedgfihkjmlon T 9P+%Q BC2`2ab c degf0h0ijk l mno p q m%r%s t.u2vwxi.jkAyz2{ | c}d}~2 v xr}~ T r m e tAu u0Zno p0q vxpryqtzs uwv xzzy|t{R m rt20|%2e rt ¡£¢^¤ ¦¥§ ¨ª©r«k¬®¯i§ok°± ²zt³kí ¡ Ff0h¢m£¤2¥.¦ §¨2©0ª « ¬0 µi¶¸· ±t¹k²z³k´ º »r¼ª²½§¿¾À · ±t¹Áº ®%¯° ±0²³´2µ.¶ ©¸· ²½Â»¦ÃrÄioÅwÀÆkÇwÈtÉÊËtÌ ¢XÍÎ ºk § · ²ÔLÕ×ÖkØioÒwÓt¾À · ÖzÑ ¹ º »2¼ ±2½ ¾¿ À Á2ÂÃÄ ÅÆ Ç±È.É Ïº±w±rÙzÐÚ ÑroÛ¿ÒwÜtÓ ¯ Ý ²z³kí æ ÊË Ì Í Ë ÎÐÏÑÒ ½ ÓÔÕÖ×ÅA± ÈËØ ÙÚ· çtèré êwë¦ìítîÞ ïk¢ß¦ðàwñákò¿âóõãô÷ärö¦økå ùr¹kúz ÀÁ ÂÑÒÖ×Å¿Û%±ÝÜ Þ ßAà0áãâåäZæç2èé ÿ

(6) û

(7) ûküýö¿þ ê2ë ìíî ï2ð2ñ òóô õOö ô"!#$%&zñ('ï)iö+*,-/.wû 10 2+34 , 5 6798):);<5 6

(8) =?>A@"BCñ ÷ø%ù úüûý íþÿ ÿ . ï 3 D 8 EGFHJILK MONQPSR T)KAUWVVXY1Z[ P]\ ^.

(9) . !#"%$& ')(%*!+,(& ')- ./%* +,01&+,& 2!3-4"!')5"% &6/- +3* (&672- 852* 9:050* - 0%('<;=* - 3910(+,- $7')>#"50 - 05$ %@?19,%-4#+,- $7AB(C>D- 0/(* - E5- /#&65')(F8!7&6(01G H I %-4#+,- $7AB(C>D$%(7C- ?!')(&6- 50D!;(* - E5- /#&65')(F8! &6(01G/;=- '#&JHLKD- E5- D$%501&6(%- 0%3(&& 2* - E5- @&6(01G,* #2 57&6/#& 'J$& 1')(%*!C- ?!')(&6- 50MHLKD- E5- @+,5&6- 50D$%(% ?19,%-4#+,- $7AB(C>"!')!5$%910(+,- $7"!')>N#1')50:& 2 AB(%* *!; & 2* - E5- @&6(01G/(0O1;61'J& 2'J+,5')>O1($%("- P>- 0!8 +,5+,01&JHRQ2910(+,- $7"!')>N#1')7')>#* &SD- 0:?!$G* - 0!83!; & 2T#2%* *SHLU:0$%>O& 272595')!910(+,- $7')>#"5072(>D( C'=9385')%(& - 0;=* 0$%- 0/;6')(%$& 1')!; & 27&6(01G HRQ2 %-4#+,- $72595')!910(+,- $7')>#"50T#25* @?- 01C>#&6- 8!(&6% &63(C!- /#$2/;6')(%$& 1')5HRQ272- 852/5')!' "!#"%$& ')(%*!+,(& ')- ./%* +,01&+,& 2!3-41&6- * - P>% ?%$%(7& 2- 05$%%/%-4#+,- $;=* - @+,5&6- 50/-4C'=9 2- 852* 9:050* - 0%('VH. _àbdcefghbijklmn)oqpr`s tuv mwyx v{z| }~)v{+I +~ v)+) `{ E d9

(10)

(11) ¡¢9£¤¥ ¦¨§ª©¦ ¤«/¬¨¯® ©¦ ¤{°±²/´³µ« ¬

(12) ¶·¸9¹º»

(13) ¼½¾¿²

(14) ÀÁ+ ²ÂÃ Ä«/¬ GÅ?°±² Æ¸)¹º» ¼{¾ ©{Ç©{Ç Ã ÄÈÉ²Ê+Ë ³GÌ¥ ¦ ÍÎ)Ï Ð{ÑÒÓ{ÔÕÖ×(Ø Ï ÐÓ)ÙÚ{ÛÝÜ/Þß Îàáâ ãqäªåæçèéê Ñë)ì íïî]ð ñòó)ô Ïõ Ñö÷ ÒÓøùÑòó)ôú û àáü9ý Ó/þÿØ

(15) !#" $&%'

(16) (*)+, - .0/1#2 3 405&60708 9 ;: < =>?&@ 1A&B CD&E01&FG =>#HI JK LM'N

(17) #( 1: < => OPRQTSMUVWVXY Z0[ \]^\_ ` Uacb&d e fgPh0VNi&jNkl m no gpP q&rtsu.

(18)

(19) "!#"$%&(')* ,+-./0-. !,#1324 54-637!8#- 9$%,5* :; <=?>A@CB -D* *EF 8&G H EF . !8# 4I H >KJCL M N*OPQSR T UV*W X YZ X [7\,]S^ _Sàbc d7e f-ghjilk*mnod-bce f-gh prqts Ttutvxwzy {?|t}~ efr ) * 8¡¢¤£¥¦ §©¨-ª«-¬ ®,¯ °± ²-³´µ¶·¸t¹ º-»¼½ ¢3£¾¸À¿ÁÂÃÄ ÅÆÇ-È ÉÊ Ë Ì3ÍÎÏÐ Ñ*ÒÓ3ÔCÕ8Ö×ØÙÚ ÛÛÜ ÝÞ"ß ÖàáÜâã-äåæ*ç8è*æ*é8êë-ì-í?îðï ÕÖ*× ñòóÖà ë-ìô-õ÷öðøÖù"úûüëý7Ü âøþÿ- ç

(20) ë Sÿ ý .

(21) !#"$%&'(*)*+-,/.021354 67289;:<=>!(@?A BC :D2E*F;G H*I5J K L2M!N

(22) O P($%;& '*Q5RSTVUW5XYZ [ \^]_5à b ` ced f!g2h\5i

(23) hjk

(24) l ]m npo qrs!tûwvx2yz5{|} ]!~ oq} ]# u c; |!#} ]`!*; yz c@^

(25) !2 | ]2g ¡c v 2¢

(26) £ ¤ ¥]*¦ y§o ¨ª©«;¬® ¯±°³² ´

(27) µ¶·!¸ ¹º »¼½V¾¿

(28) À2Á ÂÃÄÅ

(29) ÆÇVÈ*É;Ê5Ë ÌÍ Î Ï Ð Ñ ² Ò±ÓÔÕ!Ö×;ØÙ!ÚÛÜÝ*ÚÞß2à

(30) á2â ãäåæâ@ç èêéêëíìîðïòñóî ôöõ±ëôöî ÷ø÷úùüûýõ±ûýî ùüëõóþ ÿ

(31) "!$#

(32) %'&)( *,+ -./0)123$4,56 7,89:;<3$4,5= @ @ @ ? ? ∂φ ∂φ ∂φ ∂φ > C B @ + B ∂B + B @ ∂θ @ + ∂A @ = ∂ 8EDF= 9 bcd φ GH

(33) I JLKNMPO QSRPT U VXWYQZU[M]\Û[T _PO ` a efgih'jkl b,mnop

(34) d eqsr

(35) tu d eEv nwxyz {|} k ~ r$ k. ∂φ = ∂ tE k d e nwxyz

(36) ,,k n. i$, $¡¢£¤¥,¦

(37) § ¨©ª «

(38) ¬® ¯° ± ³ ∂η ∂φ ∂φ ∂η ² ∂φ ∂η ± −´ ¶ = µ − ´ ´ ∂¸¹º ∂ » ¬∂¼∂® ¯° ∂θ ∂θ · ² ¾ Â ³ Â Ã ∂φ À ² Á =½ ¶ + η + ¿ (∇φ ⋅ ∇φ ) + Á ⋅ ¾ ¶ Ä∂ η ÅiÆÇ$È,ÉÊ$ËÌÍ$ÎÏÑÐÒÓ ÔÕ,ÖÊ ×ÙØ ÅÛÚÜÞÝàßâáãäåÆæÇiÈçÉEèé êLëíìïî î ð ñïòZñïó ô ñïõPö]÷Ñø[÷ùõPõûúýü^÷ùøþöPÿP õ ò [ô òZñ

(39) . !#"%$'& ( )+*,.-/+012 3#45 )76#8'9:;<')%=>?7@A BDC#E FG#H (JI KDL ,MNO#P QSRUT V BWX7OL,#Y[Z \%] Y^_à%)b'c d ] ^#e f'ghiIjL',Y T k L',M%)#< lnmporq's'tu%vw t oxy s'tz {}|~

(40)

(41) D' x.% %7+ mo = ρ  ∂φ − (∇φ ⋅ ∇φ ) −    ∂  ' x7 ' vD' ¡ ¢#£ ¤¥¦ § ¨ ª ´ η π ±³² ¬¯® « « ° θ θ =∫∫ © © ª η π « « µ ±³² = ∫ ∫ ° ¬¯® θ θ ° © © ¶¸·º¹ σ »½¼ ¾ σ ¿ tranformationÀ ÁÂ#Ã%ÄÆÅÈÇrÉËÊÇDÌÍ ¨ Î ÏÐ7Ñ ÒÓÔ ÕÖ#× ÅÙØÚÛ Ü Ý Þrß%à áâã7äåæ7ç è#é#ê ë'ìíîDï éê}ðãäñ èDò.ó ë'ì ç%ô#õö'÷ σ øúùû7üýþ ÿ îï Þ

(42) ð σ ë%ì+í $ó #% &'(*) Philips "!. + , σ*- .0/$1*23547658"9:<;= ' σ>?@A D +C σ = B E FHG I I V W /"XLYZ [5A\7JL] KLM<N7OP7Q0R S T η U -*. `b b b ` `b b b e g de g d θ ` ` ` e b c αg a αc ηg d e αc d b c βθ a βc ηθ σ b β b ` b _ a e αb c g b a d βb c b α f e c d c c b σb ^ σ θ σ hHi s = r n q + lm + k p o j t v ∂η v ∂η u α = −(σ + ) | wyx{z β = −(σ + ) θ }*~7∂" 5$7<∂  ∂Φ α ∂Φ  ∂Φ γ ∂Φ + ∂σ + η +  ∂ + ∂σ  ∂  ∂Φ β ∂Φ   ∂Φ  +  + σ=  +   + ⋅ ∂σ   ∂θ  ∂σ  H ∂η γ = −(σ + ) ¡*¢7£¤$¥ ¦ §"¨ª©5«$∂¬7<®¯° ∂ Φ + ∂.  +  +. ∂  ∂. ∂Φ + ∂. −. ∂ Φ + ∂ ∂. ∂ Φ ∂. ∂  + ∂ . +.  ∂ Φ +  ∂ ∂  . ∂  ∂.    +   . −. ∂   ∂Φ  ∂  ∂.   + . ∂ Φ  = ∂ .

(43) . . 9999=9=ma ¡ ¢ £¤m¥¦§¨©¤m¥ª «¬®¯ °²±´³¶µ·g¸¹º » ¼ ½e¾ ¿ÀÁ Ä Å = ËÍÌÏÎ  Ã Â π  Æ Ã = Ê Æ É ÆÈÇ Ç Ç Æ Â Ð.Ñ N+1 ÒgÓÔ ÁeÕgÖ×ÙØÚÛ ÜÝÞ ßáàâaãgäåæâqçèêé?ëÙìí î ïð ñ ò ó ôõö ÷ ø ù ñ ú ù û üýþ ÿ õ þ . . a b ¶h ¶F ¶h =

(44) -

(45) -

(46) ¶ Ł ¶q ł ¶s ¶F ¶h ¶F ¶h s= - . - ¶ ¶ ¶q ¶q "!$#%'&(*)+,& -."/0,1 2436587$9:<;>=@?A:B@C<DE$FGH$IKJL@M<N O ?PN . 3RQSHA: -$T.UVPWX$Y$Z\[^] Q_a` b@c Vd ? 0egfhijk,lAmonpoXSYSZ ?<q U 3srtu$vgw@x$y h z|{R}~86 $$ A$ a,KSa<ga< ¡£¢ ¤ ¥a¦§$ö©@ªa«¬A® a¯°±@²´³¶µ·¸$¹ º» ¼½¾A¿À½ ÁÃÂÃÄÆÅÈÇÊÉ ËÍÌÈÎÐÏÑÅÊÉ ÒÔÓÕÔÂÍÖ ×Ø<ÙÛÚÝÜ@Þ ß ¿À½>àP×a¿Àáaâ ÁÃÂÃÄÆÅÈÇÊÉ ËÍÌÈÎÑÌÈÇÈÇÊãäËÍÌÈÇæåÖ × çè$é êëíì<îïSðíñ@òóoôAõgö@÷ øùúû'üý<þoÿ üýAÿ.

(47) ! " #%$%&(')+*,$%&('-&. /01 23547698;: <>=?A@CB D ? EF ú,û AG H7I JKML NPO QR ST UMNPO V W(XZY\[^]`_acb%dè f-a(Y ghji U kAl;m>nporqPstu v f-_wje [%w x y{z [%W|(Y z [%} ~(f+ |acf-e XY Q9R A9 AAMP ¡; ¢ £ ¤¥;¦ §¨M© ªr« ¬®{¯°±P²³M´ µ·¶¹¸»º·¼¾½ ¡F¯ ¿,À-ÁÂjÃ Ä%Â Å Æ{ÇcÄ%ÈÉ(ÊËÇcÄ%Ì Í À+Î ÉÏcÀ-ÐÃ ÑÎÊ ÒPÓ Ô9ÕCÖ× ØÙÚÜÛ íî7ï ÝðÞñßò7à ó á. âãäCå æ çcè é êë ç ì. Fourier Pseudospectral Matrix u(x) 0 π. !"#. ü ý = ù û π û = ø ÷ö ô ô ô ô ö ú − õ ô þpÿ ú

(48) %) %+.

(49) "!# + % * & () ' % $ = -/. +, 6 1 34 2 576 34 2 5 1 0 PacbeM Q d8R fh9;BTg :=Si=UW<jl>@VYX"?@k=Z Am B nC D EF o@GC H I D p=J\EF G H [I qrR J] ^ m D K _/Ls ` t uwM ^ ND v J\x [ O y

(50) z{|}~z{ " = ( ( )) op T" ¬ £ § ¥ ¦ « ª ¢ ª § ¡ (© ¡ ¤ (¨ ¦ )) ¯®°±r²³´¦µ ¶ ·"¸e¹ Á ¿ » À ½¾ ¼ » ¿= º ¡ Á Â Ã/Ä Ë Æ ÈÉ Ç = Ê\Ë ÈÉ Ç Ê Æ Å íî ÌÎÍ Ï ÐÒÑÓÕÔ×Ö ÑØÙÚÑÛÜÖ ÝÞÖÕßáàâãäåçæ èêéTëeì ö óõò = ö óôò ö ðñò ï ÷ø"ù"úûü ýÿþ ü =. −. =. −.

(51)

(52) !#"$#% *),+.-/.0214365*7 σ &'( -8:9; <=>@?BA*C,D*E,FGH IKJLMN OM PQR SUTWVYXWZ\[^]:_ FGY`#H a b

(53) Ocdeb

(54) Of PQR SUT VYXWZ6gihj[lkYmn.oqpr.sut#v F*w.x XWZ6y:V z {#|Y}~*, *.,:q.q ~#W.# :W.#{|.: q:qU ¢¡U£¤W¥¦#§#¨ *:©Uª «u¬K ®:¯q°± ¥ #²³q ¥´ }µ¶·ª ¤¸º¹*» °q± ¼ ½.¾Y¿WÀÁÃÂ Ä#ÅYÆÇÈ*É#Ê#¿#ËWÌ*ÍÏÎ. '( & ) * $ = ,.#/ A B C D 3 4 -FFEGH # Fourier 0123 4 57698;:=<?> @ N RV Q X P Ð@Ñ^Ò#ÓÔqÕÖ [ X Q7MSP = Z L U ∑NO M K Y J U W I TW V U NO ×,ØqÙÚ6ÛYÜ#Ý.Þ È ß6àqáßqâ6ã#ä å æ ç6è#éê ìí.îYïqð,ñ*òYó.ôYõ÷öùøú:ûñ*ò:õ.ü ãý*þ \.] ^ _a`bcedgfh i j kml n (-1) o jqp ë ÿ ôYü

(55) Yü ã ò

(56) Yó r stu h j v ô:õ ö !"!$#&% '(*)+,-*./0132 y x z{ | } z{ y * ( 7 5 9 6 3 8 : ; <=> &? @A% B 4CDEFGHI 4 y | w z{ = x~ C* H JK3LNMPO3CDEFGH > S Z []\_^ Y Q X = R TVU W (ωQ ) ` . = acb d ge f9hji ω ekml7nogpqrstcuwvmxzy (1). −. =−. − π.

(57)

(58) 2. Canuto, C., Yousuff, H.M., Quarteroni, A. "!$#&%'( and Zang, T.A. 1988 Spectral Methods in Fluid * Dynamics. Berlin-Heidelberg:Springer-Verlag. 4 4 w + = [ + 5 ,.-0/21 ( + 3 )]) 6 3. Chen, B.F. and Chiang, H.W. 1999 Complete " 4 + = p 7 89;: <=;>?@BAC&DEFGH w 2D and fully nonlinear analysis of ideal fluid in P G H ? S R E B F T W U S V X " Y A tanks. Journal of Engineering Mechanics 125, IJBKLNM;O L wQ 70-78. ZP[N\ 4. Chen, W. 1996 Large amplitude liquid ]^B_àbWcBd e >W?@"fSgGWHShij MkO sloshing in seismically excited tanks. Earthquake GHglA&mnoSgqpBrsA;tuSEF \wvSx y"z Engineering and Structure Dynamics 25, 653-669. M.J. 1999 3-D Pseudospectral Model of { abS|}~B YSFGHSh w 5. Free Chern, Surface and Viscous Flows. D.Phil. thesis, w q ¡q¢£&¤¥ ¦ University of Oxford, Oxford. §¨B©&§Bªq«¬S® °¯±S q ² S³´ © µ 6. Chern, M.J., Borthwick, A.G.L., and Eatock ¶q·¸¹qº¼»2½¾¿ÀÂÁkÃºÂÄÅÆÈÇÉ¾¿·¸ R. 1999 A pseudospectral ¹qºÊË;ÌSÍÎ ÌBÏÐ»Ñ¶Ò Ó ÔÕÖ×kÌSÍØÙ&Ú Taylor, sigma-transformation model of 2D nonlinear Û Å ÜÝSÞÄ ·¸¹qºÚßáàâãÜ·¸¹qºÂÄ;Ì Journal of Fluids and Structures 13(5), Ïä åÎæçè»é×kÌSÍØÙ&ê ëì w í îïð wñ òÈóô waves. 607-630 ¿ õö÷qøù ÷

(59) 7. Chern, M.J., Borthwick, A.G.L. and Eatock w îï wñ ô ò úSó ô ûõüöýP÷qþ øÿ ù Taylor, R. 2000 Simulation of free surface motions in a cylinder tank using Chebyshev-Fourier !"#

(60) $%&')( *

(61) ( * + ! ,.-/01243546 pseudospectral elements. Accepted by 7489:;< 5=>@?BADCE4FHGJI5=>@?HK0 International Journal for Numerical Methods in 1LMONPLHQLHRSAUTV6W < 5=>@?XZY#[UI Fluids. Faltinsen, O.M. 1978 A numerical nonlinear >@?H\J]@^)_ÀbaJcdeHfgHh3546 ij5kl 8. method os sloshing in tanks with / NPmn <oJpqJr<s KHt 8 >@?u[E^vZw x two-dimensional flow. Journal of Ship Research , y4zJ{S|~}J Py@|~} 193-202. @HH#} 4J4J 9. Gropp, W. and Lusk, E. 1999 User’s Guide for HH#b @¡J}H@¢J£¤¥¦ mpich, a Portable Implementation of MPI version §| 1.2.0. Argonne National Laboratory, Rev B. ¨ © ªJ«J¬® ¯° w ±²´³ x wµ ¶¸·¹ º 10. ANL/MCS-TM-ANL-96/6 Gupta, R.K. and Hutchinson, G.L. 1989 »¼¾½¿ÁÀÃÂÄ¾ÅÆ interaction in liquid storage tanks. ÚBÛ º »@¼½ ×@ØHÜÝÞ¶ ßÇHàÈ ÖÉ ÊÌÙ Ë ÍÏÚHÎ Ð á ÑÓÒÕÞ4Ô Ö)â×ãJØ äåÙ Solid-water Journal of Sound and Vibration 135, 357-374. æçZèPéOêPë§Û 11. Haroun, M.A. and Abou0Izzeddine, W. 1992 Parametric study of seismic soil-tank interaction. ìîíðïJñ I: Horizontal Excitation. Journal of Structural 118, 783-797. òóõôö÷øùúHûüJý)þ ÿ 12. Engineering Huang, Z. and Hsiung, C.C. 1996 Nonlinear

(62) !#"$% &('*),+.- /0+ shallow-water flow on deck. Journal of Ship 12 35467 8 9 : ; < = > ? @BA /0C,DFE A /0G 1 2 3 Research 40, 303-315. 4 6 7 H I J K ? L M N O 1 P Q R S T U V W X 13. Hussain, B. 1997 Seismic Design and Retrofit of Y H Z [ \ ] ^ 1 P _ ` a b c d e f g h i j k Cylindrical Liquid Storage Tanks. Ph.D. thesis, ` k l m n oqp r s g t u v f w x y z { b c d e University of California, Irvine. g h i j | } ~ v m d 14. Jones, A.F. and Hulme, A. 1987 The hydrodynamics of water on deck. Journal of Ship g h i j l m } g d g t w x g h u v Research 31, 125-135. ¡ ¢ £ ¤ ¥ ¦ § ¨ 15. Ku, H.C. and Hatziavramidis, D. 1985 Solutions of the two-dimensional Navier-Stokes © ª « ¬ ® ¥ ¯ ° ± ² ³ ´ µ ¶ ¨ © equations by Chebyshev expansion methods. Computers & Fluids 13, 99-113. ·¹¸»º¼½ ¾ 16. Ku, H.C., Hirsh, R.S., Taylor, T.D. and 0. 1.. Barton, D.C. and Parker, J.V. 1987 Finite element analysis of the seismic response of archored and unanchored liquid storage tanks. Earthquake Engineering and Structure Dynamics 15, 299-322. 17.. Rosenberg, A. 1989 A pseudospectral matrix element method for solution of three-dimensional incompressible flows and its parallel implementation. Journal of Computational Physics 83, 260-291. Lui, A.P.C. and Lou, J.Y.K. 1990 Dynamic.

(63) coupling of a liquid-tank system under transient excitations. Ocean Engineering 17, 263-277. 18. Ma, Q.W. and Wu, G.X. 1995 Second order transient waves around a vertical cylinder in a tank. Journal of Hydrodynamics Ser. B 4, 72-81. 19. Manos, G.C. and Clough, R.W. 1985 Tank damage during the way 1983 Coalinga earthquake. Earthquake engineering and structural dynamics 13, 449-466. 20. Mellor, G.L., Blumberg, A.F. Modelling vertical and horizontal diffusivities with the sigma coordinate system. Monthl Weather Review 113, 1379-1383. 21. Philips, N.A. 1957 A coordinate system having some special advantages for numerical forecasting. Journal of Meteorology 14, 184-185. 22. Solaas, F. and Faltinsen, O.M. 1997 Combined numerical and analytical solution for sloshing in two-dimensional tanks of general shape. Journal of Ship Research 41, 118-129. 23. Tang, Y. 1993 Dynamic response of tank containing two liquids. Journal of Engineering Mechanics 119, 531-549. 24. Veletsos, S.A. and Tang, Y. 1987 Rocking response of liquid storage tanks. Journal of Engineering Mechanics 113, 1774-1795. 25. Veletsos, A.S. and Shivakumar, P. 1993 Sloshing response of layered liquids in rigid tanks. Earthquake Engineering and Structural Dynamics 22, 801-821. 26. Veletsos, A.S. and Shivakumar, P. 1995 Hydrodynamics effects in rigid tanks containing layered liquids. Earthquake Engineering and Structural Dynamics 24, 835-860. 27. Wu, G.X., Ma, Q.W. and Eatock Taylor, R. 1998 Numerical simulation of sloshing waves in a 3D tank based on a finite element method. Applied Ocean Research 20, 337-355. 28. 1988.

(64) 29. 1997 $&%'()+*+,-

(65) .0/ 132 4 6!5070"# 8!9: ;<+=3>@?3A0BDCFE 30. K 1997 K L MON&P

(66) QORJSUT+VXWDY K[Z \] G=HJ^&IF _. ` eJf+g K 1999a KbUhcJi

(67) d j?3kA0l3BDmCFn3E oprq@st

(68) u0v w xzy&{ |}+~&

(69) 3@30 DF. 3 . h. . s. . ¡ £¢ s ¦ ¤ ¥ w=0.1w0 w=0.4w0 w=0.5w0 w=0.6w0. 1 0.75 0.5 0.25. h /a. 0. -0.25 -0.5 -0.75 -1. § ¨ª© « ¬ q ® ¯±° ²³+´µ¶ · ¸¹ 0. 10. 20. Nondimensional time. 30.

(70) Z. t*=3 0.04. w=0.1w0 w=0.4w0 w=0.5w0 w=0.6w0. 0.03. X. Y. 0.02 0.5. 0. 0. -0.01. -0.5 -1. -0.02. -1 -0.5. -0.03. -0.5. Y. 0. 0. 0. 10. 20.

(71) Nondimensional time. 30. 1. 1. F. M ed orc. on oti. # $ % & w ' (*) w+ , -/.0'21 3. Z. t*=4. 1.5. w=0.1w0 w=0.4w0 w=0.5w0 w=0.6w0. 1. X. Y. 0.5. Fx/gh2Ro. 0.5. 0. 0. -0.5. -0.5 -1. -1 -0.5. -1. -0.5. Y. 0. 0. 0. 10. 20. !" Nondimensional time. 30. X. 0.5. 0.5. -1.5. h /a. -0.04. X. 0.5. 0.5. 1. 1. F. M ed orc. on oti. # $ 4 & w ' (*) w+ , -/.0' 3. Z. t*=5 X. Y. 0.5. 0. -0.5 -1. -1 -0.5. -0.5. Y. 0. 0. X. 0.5. 0.5 1. 1. F. # $ 1 & w ' (*) w+ , -/.0' 3. M ed orc. on oti. h /a. h /a. h /a. 0.01.

(72) Z. Z. t*=9 Y. X. Y. 0.5. 0.5. 0. 0. h /a. X. -0.5 -1. -0.5. -1 -0.5. -1. -0.5. Y. 0. 0 0.5. 0.5 1. 1. -1 -0.5. -0.5. Y. X. on oti dM e c r Fo. 0. 0 0.5. 0.5 1.

(73) w w . 1. F. X. M ed orc. on oti.

(74) w w . Z. t*=10. Z. t*=7. h /a. t*=6. X X. Y. Y. 0.5. h /a. 0 0. h /a. 0.5. -0.5 -1. -0.5 -1. -0.5. -1 -0.5. -0.5. Y. 0. 0 0.5. 0.5 1. 1. Y. dM rce Fo. on oti. X. Y. 0. -0.5 -1 -0.5. Y. X. 0.5. 0.5 1. 1. F.

(75) w w . M ed orc. on oti. h /a. 0.5. -1. 1. F. X. M ed orc.

(76) w w ! . Z. t*=8. 0. 0.5 1.

(77) w w . 0. -0.5 0. 0 0.5. X. -0.5. -1. on oti.

(78)