Chapter 4 Multifield Surface Bone Remodeling

Similar documents
A Study of the Solutions of the Lotka Volterra. Prey Predator System Using Perturbation. Technique

Life Science Journal 2014;11(9) An Investigation of the longitudinal fluctuations of viscoelastic cores

16.512, Rocket Propulsion Prof. Manuel Martinez-Sanchez Lecture 3: Ideal Nozzle Fluid Mechanics

Jonathan Turner Exam 2-10/28/03

Math 266, Practice Midterm Exam 2

Single Correct Type. cos z + k, then the value of k equals. dx = 2 dz. (a) 1 (b) 0 (c)1 (d) 2 (code-v2t3paq10) l (c) ( l ) x.

UNSTEADY HEAT TRANSFER

CS 541 Algorithms and Programs. Exam 2 Solutions. Jonathan Turner 11/8/01

1 Finite Automata and Regular Expressions

The model proposed by Vasicek in 1977 is a yield-based one-factor equilibrium model given by the dynamic

Mathcad Lecture #4 In-class Worksheet Vectors and Matrices 1 (Basics)

Library Support. Netlist Conditioning. Observe Point Assessment. Vector Generation/Simulation. Vector Compression. Vector Writing

Revisiting what you have learned in Advanced Mathematical Analysis

3.4 Repeated Roots; Reduction of Order

The Laplace Transform

Inventory Management Model with Quadratic Demand, Variable Holding Cost with Salvage value

Chapter 3. The Fourier Series

Decline Curves. Exponential decline (constant fractional decline) Harmonic decline, and Hyperbolic decline.

Midterm. Answer Key. 1. Give a short explanation of the following terms.

Boyce/DiPrima 9 th ed, Ch 7.8: Repeated Eigenvalues

Global Solutions of the SKT Model in Population Dynamics

More on FT. Lecture 10 4CT.5 3CT.3-5,7,8. BME 333 Biomedical Signals and Systems - J.Schesser

Transfer function and the Laplace transformation

Lecture 21 : Graphene Bandstructure

Laplace Transform. National Chiao Tung University Chun-Jen Tsai 10/19/2011

RUTH. land_of_israel: the *country *which God gave to his people in the *Old_Testament. [*map # 2]

ktmuwii INDEPENDENT IN Al.t THINCIS. NEUTRAL IN NOTHING* Sold at Cast. AI.GE" IS DKVI). Lowell's Bright Boy Stricken With Small Pox at Manila.

Fourier Series and Parseval s Relation Çağatay Candan Dec. 22, 2013

CSE 421 Algorithms. Warmup. Dijkstra s Algorithm. Single Source Shortest Path Problem. Construct Shortest Path Tree from s

Instructors Solution for Assignment 3 Chapter 3: Time Domain Analysis of LTIC Systems

The Procedure Abstraction Part II: Symbol Tables and Activation Records

Math 3301 Homework Set 6 Solutions 10 Points. = +. The guess for the particular P ( ) ( ) ( ) ( ) ( ) ( ) ( ) cos 2 t : 4D= 2

Final Exam : Solutions

International Journal on Recent and Innovation Trends in Computing and Communication ISSN: Volume: 5 Issue:

UNSTEADY STATE HEAT CONDUCTION

Control Systems. Modelling Physical Systems. Assoc.Prof. Haluk Görgün. Gears DC Motors. Lecture #5. Control Systems. 10 March 2013

a dt a dt a dt dt If 1, then the poles in the transfer function are complex conjugates. Let s look at f t H t f s / s. So, for a 2 nd order system:

where: u: input y: output x: state vector A, B, C, D are const matrices

Chapter 5 The Laplace Transform. x(t) input y(t) output Dynamic System

Advanced Engineering Mathematics, K.A. Stroud, Dexter J. Booth Engineering Mathematics, H.K. Dass Higher Engineering Mathematics, Dr. B.S.

Analytical Solution of A Differential Equation that Predicts the Weather Condition by Lorenz Equations Using Homotopy Perturbation Method

P a g e 5 1 of R e p o r t P B 4 / 0 9

UNIT #5 EXPONENTIAL AND LOGARITHMIC FUNCTIONS

3+<6,&6([DP. September 29, SID (last 5 digits): --

LINEAR 2 nd ORDER DIFFERENTIAL EQUATIONS WITH CONSTANT COEFFICIENTS

K x,y f x dx is called the integral transform of f(x). The function

IJRET: International Journal of Research in Engineering and Technology eissn: pissn:

HIGHER ORDER DIFFERENTIAL EQUATIONS

Chapter4 Time Domain Analysis of Control System

DEPARTMENT OF ELECTRICAL &ELECTRONICS ENGINEERING SIGNALS AND SYSTEMS. Assoc. Prof. Dr. Burak Kelleci. Spring 2018

A Tutorial of The Context Tree Weighting Method: Basic Properties

Section 2: The Z-Transform

Poisson process Markov process

Pupil / Class Record We can assume a word has been learned when it has been either tested or used correctly at least three times.

Shortest Paths. CSE 421 Algorithms. Bottleneck Shortest Path. Negative Cost Edge Preview. Compute the bottleneck shortest paths

PHA Second Exam. Fall 2007

Week 06 Discussion Suppose a discrete random variable X has the following probability distribution: f ( 0 ) = 8

2. The Laplace Transform

Fourier. Continuous time. Review. with period T, x t. Inverse Fourier F Transform. x t. Transform. j t

PFC Predictive Functional Control

( ) ( ) + = ( ) + ( )

Applications of these ideas. CS514: Intermediate Course in Operating Systems. Problem: Pictorial version. 2PC is a good match! Issues?

EE Control Systems LECTURE 11

Charging of capacitor through inductor and resistor

Last time: introduced our first computational model the DFA.

Systems of First Order Linear Differential Equations

Introduction to Laplace Transforms October 25, 2017

Equations and Boundary Value Problems

Stability and Optimal Harvesting of Modified Leslie-Gower Predator-Prey Model

Price Dependent Quadratic Demand Inventory Models with Variable Holding Cost and Inflation Rate

REPETITION before the exam PART 2, Transform Methods. Laplace transforms: τ dτ. L1. Derive the formulas : L2. Find the Laplace transform F(s) if.

Boyce/DiPrima 9 th ed, Ch 2.1: Linear Equations; Method of Integrating Factors

Systems of First Order Linear Differential Equations

FL/VAL ~RA1::1. Professor INTERVI of. Professor It Fr recru. sor Social,, first of all, was. Sys SDC? Yes, as a. was a. assumee.

Combinatorial Optimization

Copyright A.Milenin, 2017, AGH University of Science and Technology

Flow Networks Alon Efrat Slides courtesy of Charles Leiserson with small changes by Carola Wenk. Flow networks. Flow networks CS 445

Engine Thrust. From momentum conservation

Elementary Differential Equations and Boundary Value Problems

CIVL 8/ D Boundary Value Problems - Rectangular Elements 1/7

Relation between Fourier Series and Transform

rank Additionally system of equation only independent atfect Gawp (A) possible ( Alb ) easily process form rang A. Proposition with Definition

14.02 Principles of Macroeconomics Fall 2005 Quiz 3 Solutions

Erlkönig. t t.! t t. t t t tj "tt. tj t tj ttt!t t. e t Jt e t t t e t Jt

Solutions to assignment 3

x, x, e are not periodic. Properties of periodic function: 1. For any integer n,

Midterm exam 2, April 7, 2009 (solutions)

Dark Solitons in Gravitational Wave and Pulsar Plasma Interaction

PHA Second Exam. Fall On my honor, I have neither given nor received unauthorized aid in doing this assignment.

Right Angle Trigonometry

(2) If we multiplied a row of B by λ, then the value is also multiplied by λ(here lambda could be 0). namely

Partial Fraction Expansion

2 T. or T. DSP First, 2/e. This Lecture: Lecture 7C Fourier Series Examples: Appendix C, Section C-2 Various Fourier Series

1K21 LED GR N +33V 604R VR? 1K0 -33V -33V 0R0 MUTE SWTH? JA? T1 T2 RL? +33V 100R A17 CB? 1N N RB? 2K0 QBI? OU T JE182 4K75 RB? 1N914 D?

CHAPTER 9 Compressible Flow

Instructions for Section 1

e t dt e t dt = lim e t dt T (1 e T ) = 1

Vidyalankar. 1. (a) Y = a cos dy d = a 3 cos2 ( sin ) x = a sin dx d = a 3 sin2 cos slope = dy dx. dx = y. cos. sin. 3a sin cos = cot at = 4 = 1

PHA First Exam Fall On my honor, I have neither given nor received unauthorized aid in doing this assignment.

INTEGRALS. Chapter 7. d dx. 7.1 Overview Let d dx F (x) = f (x). Then, we write f ( x)

Transcription:

hr Mulifild Surf on Rmodling In hr, h horil nd numril rul of inrnl on rmodling wr rnd. Exnion o mulifild urf on rmodling i diud in hi hr. horil rdiion of urf on rmodling in h dihyi of h long on undr vriou xrnl lod r md wihin h frmwork of div li hory. h lod inlud xrnl lrl rur, lri nd hrml lod. wo oluion r rnd for nlyzing hrmolroli rolm of urf on rmodling. h nlyil oluion h giv xlii formulion i l of modling homognou on mril, nd h mi-nlyil oluion i uil for nlyzing inhomognou. umril rul r rnd o vrify h rood formulion nd o how h ff of mhnil, hrml nd lri lod on h urf on rmodling ro. In hi hr, h dvlomn in [-] r rifly ummrizd.. Soluion of urf modling for homognou hollow irulr ylindril on In hi ion, nlyil oluion for hrmolroli rolm of urf on rmodling, d on h hory of div liiy [], r rnd o udy h ff of mhnil, hrml, nd lri lod on h urf on rmodling ro. h oluion i ud o invig h urf on rmodling ro on h i of uming homognou on mril [5]... R quion for urf on rmodling h quion for h hory of div liiy. n xndd o inlud izolri ff y dding om nw rm low [,]

U n, Q ij Q, Q, n, Q E Q, E Q, ij zz ij zz rz rz i E r r i E z z i. whr zzzz urf rmodling offiin. E E rz rz r r z z, i r.. Diffrnil fild quion for urf rmodling W now onidr gin h hollow irulr ylindr of on hown in Fig... h hollow irulr ylindr i ujd o n xrnl mrur hng, qui-i xil lod P, n xrnl rur, nd n lri onil lod or/nd hown in Fig... h oundry ondiion r dfind y Eq.5 nd.6. h oluion of dilmn u r, u z nd lri onil o h rolm ov wr givn y Eq.9-.6 in hr. Suiuing Eq.-.6 ino Eq. yild U ln ln U ln ln whr..

* * F F zz zz. * * F F zz.5 * P F zz.6 r zr 5.7 * * F F zz zz.8 * * F F zz.9 * P F zz.

5 zr r... nd h uri or ui nd rfr o riol nd ndol, rivly. For h ron d in Sion.., Eq. n ud hr, h i U d, d d U. d hu, Eq. nd. n wrin d d ln ln.5 d d ln ln.6 whr... Aroximion for Smll hng in Rdii I i rn h Eq.5 nd.6 r non-linr nd nno, in gnrl, olvd nlyilly. owvr, h quion n roximly linrizd whn hy r lid o olv rolm wih

mll hng in rdii. In h on urf rmodling ro, Qin l [] umd h h rdii of h innr nd our urf of h on hng vry lil omrd o hir originl vlu, w don in Sion... hi mn h h hng in nd r mll. hi i ronl umion from h viwoin of h hyi of h rolm. o inrodu h roximion, wo nondimnionl rmr [,],.7 r dod in h following lulion. A rul, nd n wrin,,,.8 Sin oh nd r fr mllr hn on, hir qur n ignord from h quion. onqunly, w n hv h following roximion:.9. ln ln...

ln. whr,, ln,.5 hu, Eq.5 nd.6 n roximly nd in rm of nd follow:, d d d d.6 whr.7.8.9...

.. Anlyil Soluion of Surf Rmodling An nlyil oluion of Eq.6 n oind if homognou rory i umd for on mril. In h, h inhomognou linr diffrnil quion ym.6 n onvrd ino h following homognou on d d, d d. y inroduing wo nw vril uh h,., d M d M.5 M.6 d M.7 h oluion of Eq., uj o h iniil ondiion h nd, n xrd in four oil form h fulfill h hyi of h rolm, i.., whn, nd mu limid qunii, nd h oluion mu l. h form of h oluion dnd on h roo of h following qudri quion: whr rm d M.8 rm.9 All h horilly oil oluion r hown follow:

A: whn, nd, Eq.8 h wo diffrn roo, nd, oh of whih r rl nd diin. hn h oluion of h quion r,. whih n lo wrin,. h formul for h vriion of h rdii, i.., nd, wih im n oind y uiuing Eq. ino Eq.7. hu,. h finl rdii of h ylindr r hn lim, lim. : whn, nd, Eq.8 h wo qul roo,. h oluion of h quion r

. whih n lo wrin,.5 h formul for h vriion of nd wih im n oind y uiuing Eq.5 ino Eq.7,.6 h finl rdii of h ylindr r hn lim, lim.7 : whn nd, h oluion of h quion r

whih n lo wrin,.8,.9 h formul for h vriion of nd wih im n oind y uiuing Eq.9 ino Eq.7, follow:,.5 h finl rdii of h ylindr r hn, lim lim.5 D: whn nd, h oluion of h quion r whih n lo wrin,.5,.5 h formul for h vriion of h rdii wih im n oind y uiuing Eq.5 ino Eq.7. hu.5 h finl rdii of h ylindr r hn, lim lim.55

All h ov oluion r horilly vlid. owvr, h fir i h mo likly oluion o h rolm, i i hyilly oil whn []. hrfor i n ud o lul h on urf rmodling.. Aliion of Smi-Anlyil Soluion o Surf Rmodling of Inhomognou on h mi-nlyil oluion rnd in Sion. n ud o lul rin nd r ny oin on h on urf. h rul form h i for urf on rmodling nlyi. hi ion rn liion of oluion.58 o h nlyi of urf rmodling hvior in inhomognou on. [ I i nod h urf on rmodling i im-dndn ro. h hng in h rdii or n hrfor luld y uing h rngulr lgorihm of ingrl Fig... h rodur i drid hr. Fir, l h ring im nd h lngh of im o onidrd, nd divid h im domin ino m qul inrvl / m. A h im, lul h rin nd lri fild uing Eq.-.6. h rul r hn uiud ino Eq. o drmin h norml r of h urf on rmodling. Auming h i uffiinly mll, w n rl U wih i mn vlu U h im inrvl [, + ]. h hng in h rdii or im n hu drmind uing h rul of urf vloiy. Aordingly, h rin nd lri fild r udd y onidring h hng in h rdii. h udd rin nd lri fild r in urn ud o lul h norml urf vloiy h nx im inrvl. hi ro i rd u o h l im inrvl m, ]. Fig.. how h rngulr lgorihm of ingrl whn w rl U wih i iniil vlu hn i mn vlu U h im inrvl [, + ]. U rhr

Fig.. Illurion of h rngulr lgorihm. Surf Rmodling Equion Modifid y n Inring Mdullr Pin Suiuing Eq.65 ino Eq.5 nd.6 yild M d d ln ln.56 M M d d ln ln.57 whr F M zz.58

F M zz.59 M.6 I n n h Eq.56 nd.57 r imilr o Eq.5 nd.6. onidring Eq.7, Eq.56 nd.57 n imlifid, Y Y Y d d Y Y Y d d.6 whr 5 M Y.6 5 M Y.6 M Y.6 6 7 5 M M Y.65 6 7 5 M M Y.66 6 Y M M.67

wih.68 F 5 F.69.7 6 F.7 F 7 F.7 Eq.6 i imilr o Eq. nd n hu olvd y following h oluion rodur drid in Sion.... umril Exml for hrmoizolri on onidr gin h fmur ud in Sion.7. h gomril nd mril offiin of h fmur r h m ho ud in Sion.6 x h h volum frion hng i now kn o zro hr. In ddiion, h urf rmodling r offiin r umd o

9.6 m /, 7. m /, 5. m /, 8. m /,.6 m /,.8 m / zr zz 9.6 m /, zr m /, nd r zz, 9 V m /.87 m/,.58 m/ nd, r umd. In h following, numril rul r rovidd o how h ff of mrur nd xrnl lri lod on h urf on rmodling ro. h rul for h ff of mhnil loding, inrd in, nd mril inhomogniy on h urf rmodling hvior r omid hr; hy n found in []. W diinguih fiv loding. n MP n=.8,.9,,. nd., P=5, nd no ohr y of lod i lid. h xndd rul for hi loding r hown in Fig.. o udy h ff of xrnl rur on h on rmodling ro. Fig.. illur h h rnvr nd xil lod hv h m ff on h on. h innr rdiu of h on dr whil h our rdiu inr h xrnl rur inr. hi rul in n inr in h on ro ionl r nd, onqunly, hikr nd rongr on. Whn h xrnl rur dr, h innr rdiu inr nd h our rdiu dr, whih mn h h on om hinnr nd wkr. On h ohr hnd, g rur n inr h r of on urf rmodling, whih n lr h rovry of injurd on. Fig.. rn n inring hng of rdiu gin im. h our rdiu of h on inr fir h rnvr rur inr. I gin o dr fr fw dy nd finlly onvrg o l vlu h i g hn i iniil vlu. A imilr rul w rnd y owin nd Vn ukirk []. h rlir nlyi w d on modl of urf rmodling nd w riud

o h uounding on mril oming mor or l iff rhr hn ing du o urf movmn..5..... -. -. -. -. 5 6 7 dy MP.9MP.MP.8MP.MP... MP.MP.MP.9MP.8MP -. -. -.6 -.8 6 8 dy Fig.. Vriion of nd wih im P5, ohr. vru ; vru

9.5, 9.8,,.,.5, V, MP, P=5 Fig.. how h ff of mrur hng on on urf rmodling. In gnrl, h rdii of h on dr whn h mrur inr nd hy inr whn h mrur dr. I n lo n from Fig.. h nd r lmo h m. Sin, h hng of h our urf rdiu i normlly g hn h of h innr on. h r of h on ro ion dr h mrur inr. hi lo ugg h lowr mrur i likly o indu hikr on ruur, whr wrmr nvironmn my imrov h rmodling ro wih l hik on ruur. hi rul m o oinid wih ul f. hikr nd rongr on rh mk ron living in Rui look rongr hn on living in Vinm. I hould mniond hr h how hi hng my ff h on rmodling ro i ill n on quion. A n iniil invigion, h uro of hi udy i o how how on my rond o hrml lod nd o rovid informion for h oil u of imod xrnl mrur fild in mdil rmn nd in onrolling h hling ro of injurd on..8.6... -. -. -.6 -.8 5 5 dy 9.5 o 9.8 o o. o.5 o

.8.6... -. -. -.6 -.8 5 5 dy 9.5 o 9.8 o o. o.5 o Fig.. Vriion of nd wih im for vrl mrur hng. vru ; vru -6V, -V, V, nd 6V, MP, P= 5, nd =. Fig.. how h vriion of nd wih im for vriou vlu of lri onil diffrn. I n n h h ff of h lri onil i xly ooi o h of mrur. A dr in h inniy of h lri fild rul in dr of h innr nd our urf rdii of h on y lmo h m mgniud. horilly, h rul ugg h h rmodling ro migh imrovd y xoing on o n lri fild. lrly, furhr horil nd xrimnl udi r ndd o invig h imliion of hi xour for mdil ri.

.6... -. V 6V V -V -6V -. -.6 5 5 dy... V 6V V -V -6V -. -. -.6 5 5 dy Fig.. Vriion of nd wih im for vrl onil diffrn. vru ; vru P=5, nd inrnl rur i rodud y inring rigid in who rdiu * i g hn. h vlu of nd gin for * -=.mm,.mm, nd.5mm r hown in Fig..5. I n n h oh h innr nd h our urf rdii of h on r rdud fr h in h n inrd, whih will inr h ighn of fi. A ighly

fid in n lo inr h r of on urf rmodling, whih n lr h rovry of h injurd on. On h ohr hnd, h rdiu of h in inr, h our urf rdiu of h on dr mor ignifinly hn do h innr urf rdiu. hi rul in dr of h on ro-ionl r nd, hn, in hinnr nd wkr on ruur. Furhrmor, if h fi i oo igh, h rur on h inrf will u dmg o h on ruur. hu h rdiu of h in hould k wihin n rori rng... -. -. -. -. -.5 -.6 -.7 -.8 = m = m = m 5 5 dy.. -. -. -. -. -.5 -.6 -.7 -.8 5 5 dy =m =m =5m Fig..5 Vriion of nd wih im indud y olid in. vru ; vru 5 A hollow, inhomognou irulr ylindril on ujd o xrnl lod.

h xrnl lod r =MP, P=5, =, nd =V. h gomril nd mril rmr of hi rolm r h m ho ud in h ov, x h ll mril onn h ginning of hi ion r now modifid y mulilir [---r/-], whr nd n rng rduion of iffn h innr urf of h on. I i worh mnioning h y uing h mi-nlyil roh, h form of iffn vriion in h rdil dirion n riy. Fig..6 how h rul of nd h ouid urf of h on for =.995,.99 nd.98. In gnrl, h rmodling r dlin h iniil iffn of h innr on urf dr. h inhomogniy of h on h ignifin ff on h on urf rmodling. A highr urf rmodling r i lwy rld o highr lvl of inhomogniy...8.995.99.98.6... 6 8 dy..8.6 =.995 =.99 =.98... 6 8 dy

Fig..6 Vriion of nd wih im for inhomognou on ujd o ouling lod. vru ; vru.5 Exnion o hrmomgnolroli Solid In h of hrmomgnolroli olid, h oniuiv quion i dfind y Eq.7 nd h rmodling r quion Eq. mu ugmnd y om ddiionl rm rld o mgni fild follow: U ij n, Q ij Q ij Q i Ei Q Ei Q G i i Q i Q zzzz rzrz r Er zez Gr r Gz z.7 whr G i i=r,z r urf rmodling offiin, nd E E G G.7 zz zz rz rz r r z z r r z z W now onidr gin hollow irulr ylindr of on, ujd o n xrnl mrur hng, qui-i xil lod P, n xrnl rur, n lri onil lod or/nd nd mgni onil lod nd/or. h oundry ondiion r dfind y Eq.6 nd.77. h govrning quion r dfind y Eq.7,.7-.76. hn h oluion of U nd U o h rolm ov hv h m form ho of Eq. nd. x h h offiin nd r rld y 5 G 5 zr r zr r.75 5 G 5 zr r zr r.76

hrfor, h oluion o nd hv h m form ho drid in Sion... A numril illurion of h nlyil oluion drid in hi ion, w onidr gin fmur diud in Sion.7 x h h volum frion hng i now kn o zro hr. In ddiion, h urf rmodling r offiin givn in Sion. r ud. h urf rmodling r offiin for mgni fild i G r m / Amr. W diinguih hr loding :.5, -.,,.,.5, V, A, =MP, P 5..8.6... -. -. -.6 -.8 5 5 dy -.5 o -. o o. o.5 o

.8.6... -. -. -.6 -.8 5 5 dy -.5 o -. o o. o Fig..7 Vriion of nd wih im for vrl mrur hng. vru ; vru Similr o h hown in Fig.., Fig.7 how h ff of mrur hng on on urf rmodling undr h ominion loding of mhnil, lri, nd mgni fild. h rul oind i imilr o ho of Fig.. nd h nlyi for Fig. li hr. -6V, -V, V, nd 6V, MP, P 5, A nd.5 o.6... 6V V V -V -6V -. -. -.6 5 5 dy

... 6V V V -V -6V -. -. -.6 5 5 dy Fig..8 Vriion of nd wih im for vrl lri onil diffrn. vru ; vru Fig.8 how h vriion of nd wih im for vriou vlu of mgni-fild diffrn. Similr o h hown in Fig.., i n lo n h h ff of h mgni-fild i xly ooi o h of mrur. A dr in h inniy of h mgnifild rul in dr of h innr nd our urf rdii of h on y lmo h m mgniud. -A, -A, A, nd A, MP, P 5, V nd..5..5a.a A -.5A -.A -.5 -. 5 5 dy

dy.a.5a A -.5A -.A Fig..9 Vriion of nd wih im for vrl mgni onil diffrn. vru ; vru Fig.9 how h vriion of nd wih im for vriou vlu of mgni onil diffrn. h hng of h our nd innr urf of h on du o h mgni influn r imilr o h in []. Rfrn [] Qin Q, Qu Y, Y JQ: hrmo lroli oluion for urf on rmodling undr xil nd rnvr lod. iomril 6, 6798-68 5 [] Qu Y, Qin Q, Kng Y: horil rdiion of urf on rmodling undr xil nd rnvr lod, 7-8. In: Pro. of 9h In onf on Inion, Aril, Rir & Minnn of Sruur, Rn WX, Ong KG, nd n JSY d.. I-Prmir PY D, 5 [] Qin Q-: Muli-Fild on Rmodling Undr Axil nd rnvr od. In: oomingon DR d. w Rrh on iomril.. 9-9 ov Sin Pulihr, w York 7 [] owin S, Vn ukirk W: Surf on rmodling indud y mdullry in. Journl of iomhni, 69-76 979

[5] Qin Q, Y JQ: hrmolroli oluion for inrnl on rmodling undr xil nd rnvr lod. Inrnionl Journl of Solid nd Sruur 9-, 7-6