Stress-compatible embedded cohesive crack in CST element

Size: px
Start display at page:

Download "Stress-compatible embedded cohesive crack in CST element"

Transcription

1 Fratur Mhani of Conrt an Conrt Strutur - Rnt Avan in Fratur Mhani of Conrt - B. H. Oh, t al.() 2010 Kora Conrt Intitut, Soul, ISBN Str-ompatibl mb ohiv rak in CST lmnt J.F. Oln & P.N. Pouln Thnial Univrity of Dnmark, Dnmark ABSTRACT: A impl lmnt ith an mb trong iontinuity for moling ohiv raking of onrt i prnt. Th lmnt iffr from prviou lmnt of th mb typ, in that a onitnt tr fil i obtain by irt nformnt of tr ontinuity aro th rak. Th iplamnt iontinuity i mol in an XFEM fahion; hovr, th iontinuou iplamnt fil i pial, alloing for th irt nformnt of tr ontinuity. Thi in turn allo for limination of xtra gr of from nary for ribing th rak formation, thu th lmnt ha th am numbr of from a it ontinuou bai: CST. Th goo prforman of th lmnt i montrat by it ability to imulat thrpoint bning of a noth onrt bam. Th avantag of th lmnt i it impliity an th traightforar implmntation of it. Hanling ituation ith multipl rak ill not rquir furthr vlopmnt, although thi i not montrat hr. 1 INTRODUCTION Th 2D moling of onrt raking in th framork of FEM ha alray rah an avan lvl, an vlopmnt in thi fil ovr th pat a hav bn haratriz a XFEM, th Xtn Finit Elmnt Mtho, (Blythko & Blak 1999), (Moë t al. 1999) an (Moë & Blythko 2002). For rfinmnt of th mtho appli to ohiv raking alloing for partly rak lmnt, (Afrg t al. 2007) an (Mougaar t al. 2009). Although th mtho ar vry fftiv, th rabak, hovr, i that it tak a lot of bookkping, at th ytm lvl a ll a at th lmnt lvl, to nur onformity an to mol paration of th iontinuou part of an lmnt auratly, i.. in uh a ay that it allo for fully inpnnt formation of th to parat part. Th lmnt hih i propo hr i a 2D lmnt of th mb iontinuity typ ith a trong iontinuity; thi mb iontinuity typ a originally introu by (Dvorkin, E.N. t al. 1990). Furthr vlopmnt r prnt by (Olivr 1996); thi ork i ba on th aum nhan train approah an iffr from th prnt ork in that it o not nfor tr ontinuity irtly, but in a variational mannr. Th avantag of th mb lmnt typ i that it i impl to implmnt an that it o not ompliat th builing of th FE ytm Equation. Th rabak of it i that it lak th ability to rat ompltly inpnnt lmnt paration an that iplamnt ontinuity i violat on lmnt i ro by a rak; hovr, in many a th fiini might b onir to b of minor importan. Th prnt lmnt tak avantag of th ytmati ontrution of iplamnt fil in XFEM, xpt that th iontinuou hap funtion ar not th am a th ontinuou. Bi from bing of th mb typ th prnt lmnt i haratriz by having a onitnt tr fil hih i ontinuou aro th iontinuity. Enforing thi tr ontinuity in th lmnt obtain th poibility of liminating th xtra gr of from (DOF) nary for ribing th iontinuity; hrby an lmnt i raliz ith no xtra DOF ompar to th ontinuou vrion of th lmnt. Th prnt lmnt allo for ohiv raking ith a tr ritrion for rak initiation, i.. a rak initiat if th uni-axial tnil trngth i x by th firt prinipl tr an no tr intnity i aum at th rak tip, auming mooth rak lour. 2 KINEMATICS Th moling of a trong iontinuity, uh a a rak, i ba on th tratgy of th xtn finit lmnt mtho for th approximation of a iplamnt fil ith a iontinuity. Thu, th iplamnt fil approximation in an lmnt ith a iontinuity i tablih by ombining th iplamnt fil orrponing to th ontinuou lmnt ith th iplamnt fil orrponing

2 to th iontinuity. For th linar intrpolation of th ontinuou iplamnt fil onir a thr-no triangular lmnt, th CST. Th lmnt i givn th poibility of a trong iontinuity along a traight lin roing th lmnt. Th nformnt of tr ompatibility btn th tr in th ontinuou part of th lmnt an th tr briging th iontinuity lin i th ntial ia of th prnt ork. Str ompatibility i nur by maning that th orr of variation of th iontinuou iplamnt fil along th iontinuity lin mut math th orr of variation of th ontinuou tr fil, an that th iontinuou iplamnt fil prou qual tr on oppoit i of th iontinuity. In th a of a CST thi la to a iplamnt fil ith a ontant jump in th iplamnt along th iontinuity lin, hih ill prou ontant briging tr along th iontinuity mathing th ontant tr fil of th CST. Furthr, thi iplamnt fil mut prou qual graint on ithr i of th iontinuity. A hap funtion that allo for a ontant jump aro a traight lin an at th am tim only introu qual an ontant train on ithr i of th iontinuity lin i hon in Figur 1. Th hap funtion may b rittn in trm of th natural oorinat ζ aoiat ith th lmnt vrtx hih i alo a vrtx of th triangular ub-omain : ζ 1 in N (1) ζ in Thi typ of hap funtion ha prviouly bn appli by (Olivr 1996). Only to DOF ar n to rib th iontinuou formation. W ho th global omponnt of th iontinuity vtor an ollt th DOF in th vtor T V [ ΔV1, ΔV2]. Eah DOF i aoiat ith a iontinuou hap funtion a givn in (1), an th iontinuou iplamnt fil vtor u i introu a u N 0 NV, N 0 N (2) By ombining u ith th ontinuou iplamnt fil vtor u N V for th normal CST lmnt arriv at th total formation fil: u u u N V N V (3) hr N i th uual linar intrpolation matrix an V th uual of vtor. J ) D ( h, T h Th atr ontnt an b xpr a of th vaporabl atr (apillary a Figur 1. Top: Straight vapor, iontinuity an aorb lin fining atr) lmnt an ubomain, an th aoiat (hmially natural oorinat boun) ζ. atr n (Mil th non- Pantazopoulo & Mill 1995). It i ra Bottom: Illutration aum of that ontant th vaporabl iontinuity atr i a fu hap funtion. rlativ humiity, h, gr of hyration Aopting th Voigt gr notation of ilia an fum auming ration, a lin- ar train maur, th ag-pnnt train in th ontinuou orption/orption part, i.. of th lmnt may (Norling b rittn Mjonll a th 1997). vtor Unr thi aum by ubtituting Equation 1 into Equati ε BV BV obtain (4) hr B i th uual ontant train intrpolation matrix an B ( D h) i th ontant train intrpolation & & h matrix riv from th iontinuou hap fun-tion. Not that B i th am on both i of th iontinuity. hr / i th lop of th orption/ iothrm (alo all moitur apa govrning quation (Equation 3) mut b 3 VARIATIONAL by appropriat AND FEM bounary FORMULATIONS an initial oniti Th rlation btn th amount of Conir a linar atr lati an boy rlativ in a humiity tat of plan i all tr or plan train. iothrm Th boy if maur i ituat in ith a plan inraing Cartian oorinat humiity ytm an ( x orption 1, x 2). Aum iothrm that a in th rak ha form a. in th Nglting boy, Figur thir iffrn 2. Th unit (Xi t al. tangnt vtor to th th folloing, urvilinar orption rak path iothrm -ilnot by, an rfrn th normal to to both th orption rak path an i orption b not by n. At any By point th ay, along th if th rak hytri ( n, ) -ofin a loal right-han iothrm oul oorinat b takn ytm. into aount, Th to th ngativ i of rlation, th rak vaporabl i fin atr a v th rlativ i humi hr n oini b u ith aoring th outar to th normal ign to of th th varia rak fa. Th tr rlativity tranfr humiity. ovr th rak Th i hap givn of th a th normal an iothrm tangntial for HPC tration, i influn σ n an by τ n, many p rptivly. Th pially ar th tho gnraliz that influn tr xtnt in an th rak, an thir hmial ork-onjugat ration gnraliz an, in train turn, trm ar th normal opning trutur of an th por rak, iz Δ uitribution n un un,(atr- ratio, lip mnt of hmial th rak ompoition, fa, SF an th tangntial Δ u u u. Th uring gnraliz tim an tr mtho, an tmpratur, train mix in th rak ar ollt t.). In in th vtor litratur aoring variou to formulatio foun to rib th orption iothrm σ r n onrt u n, Δ (Xi t al. 1994). Hovr, in th σ papr τ u th mi-mpirial xprion (5) pro n Δu Norling Mjornll (1997) i aopt b Th proportionality offiint D(h,T) moitur prmability ζ an it i a nonlina of th rlativ humiity h an tmpratur & Najjar 1972). Th moitur ma balan that th variation in tim of th atr ma volum 1 of onrt (atr ontnt ) b q ivrgn of th moitur flux J t J Proing of FraMCoS-7, May 23-28, 2010

3 hr J D ( h, ut ) hignifi a jump in th iplamnt (1) fil u rlat to th loal ( n, ) oorinat ytm. Th proportionality offiint D(h,T) i all moitur Th tr prmability briging an th it rak i a ar nonlinar typially funtion funtion of th of rlativ th rak humiity opning h an an tmpratur liing a ll T (Bažant a of th & Najjar hitory 1972). of opning Th moitur an liing. ma At balan thi point rquir aum that th that variation th briging in tim tr of th atr may b ma rittn pr a unit a funtion volum of of onrt u : (atr ontnt ) b qual to th ivrgn of th moitur flux J r r σ σ ( u ) (6) J I (2) Th t intrnal virtual ork δw an th xtrnal E virtual ork δ W may b tat a Th atr ontnt an b xpr a th um of I T T r δw th vaporabl δεσ δ S atr u σ (apillary atr, atr (7) S vapor, an aorb atr) an th non-vaporabl E T T (hmially δw δu f boun) δ Γ u t atr (8) Γ n (Mill 1966, t Pantazopoulo & Mill 1995). It i raonabl to aum that th vaporabl atr i a funtion of hr σ an ε ar th ork onjugat tr an rlativ humiity, h, gr of hyration, train vtor, f i th omain loa vtor,, an t i gr of ilia fum ration, th prrib bounary tration, i.. vtor, an (h, th, prfix δ not th variation of th ubqunt fil. ) ag-pnnt orption/orption iothrm (Norling Mjonll 1997). Unr thi aumption an Th un-rak bulk matrial i aum to b linar lati an ontitutiv Equation i rittn in by ubtituting Equation 1 into Equation 2 on obtain th form σ D ε (9) ( D h) & & & n (3) h hr D i th appropriat matrial tiffn matrix. hr / i th lop of th orption/orption iothrm Str ompatibility (alo all aro moitur th iontinuity apaity). Th i r nur govrning by quation maning (Equation th tration 3) mut t b on omplt th rak fa by appropriat to qual th bounary briging an tr initial onition. Th rlation btn th amount of vaporabl r r r (10) T atr t an σ rlativ humiity i all aorption iothrm if maur ith inraing rlativity r hr humiity T an i orption a tranformation iothrm matrix in pning th oppoit on th a. orintation Nglting of thir rak. iffrn Th tration (Xi t al. may 1994), b trmin th folloing, from orption iothrm ill b u ith in rfrn to both orption an orption onition. By th ay, if n th 1 0 hytri n r 2 of th moitur iothrm t mσ, oul m (11) b takn into 0 n2 n aount, to iffrnt 1 rlation, vaporabl atr v rlativ humiity, mut b u aoring to th ign of th variation of th hr m i a matrix ba on th lmnt of th rlativity humiity. Th hap of th orption normal vtor of th rak fa n ( n1, n2). Introuing into (6), (10) an (11) th lmnt fil ap- iothrm for HPC i influn by many paramtr, pially tho that influn xtnt an rat of th proximation (4) an (9), an utilizing that hmial r ration an, in turn, trmin por u T V th tr ompatibility rquirmnt trutur an por iz itribution (atr-to-mnt may b xpr at th lmnt lvl a ratio, mnt hmial ompoition, SF ontnt, T uring r md tim B V an B mtho, V tmpratur, σ r ( T r V mix aitiv, t.). In th litratur variou formulation ) an (12) b foun to rib th orption iothrm of normal onrt Th natur (Xi t of al. thi 1994). Equation Hovr, pn in on th (6), prnt an in papr gnral th it mi-mpirial i nonlinar in xprion V. Hovr, propo it may b by olv Norling for Mjornll V at th (1997) lmnt i lvl, aopt thu bau alloing it for xpliitly th limination aount of for th th DOF volution ribing of hyration th iontinuity. ration an Th SF variation ontnt. of Thi V may orption b tablih iothrm by ra taking th variation of (12): 1 ( h,, ) G (, ) ( g ) h 1 Figur hr 2. th Global firt oorinat trm (gl ytm iothrm) ( x1, x2) an rprnt loal rak th oorinat phyially ytm boun ( n, )(aorb). atr an th on trm (apillary iothrm) rprnt th apillary atr. Thi xprion i vali r r δ δ σ only u for lo ontnt T md B V B V δv of SF. Th offiint G (13) 1 rprnt u V th amount of atr pr unit volum hl in th gl por at 100% rlativ r W introu humiity, D an a it an th b gnraliz xpr tangntial (Norling tiffn Mjornll of 1997) th rak a briging, i.. r r G (, σ) k k 1 vg vg D u (5) (14) hr k vg an k vg ar matrial paramtr. From th r maximum Ralizing amount that of uatr V pr T unit, th volum variation that an of V fill all may por b iolat (both apillary from (13) por in th an form gl por), on an alulat K 1 a on obtain δv ZδV (15) 10 g h G ( r ) r r 1 K r, 1 10 g h 1 1 hr 10( g ) h K (, ) (4) (6) Z D T T mdb T mdb (16) Th approximation to th intrnal virtual ork in an lmnt Th matrial may no paramtr b xpr k vg an by k vg an g 1 an b alibrat by fitting xprimntal ata rlvant to fr I T T T T T rt r δw δ V (vaporabl) [ ] S B Z B atr σ ontnt Z T in σ onrt (17) at S variou ag (Di Luzio & Cuati 2009b). hr all quantiti no rlat to an lmnt. Th tr 2.2 Tmpratur that ntr volution into thi xprion ar trmin Not that, a: at arly ag, in th hmial ration aoiat ith mnt hyration an SF ration σ ar DBV xothrmi, BV th tmpratur fil i not uniform (18) for r r r σ non-aiabati σ ( T V ytm vn if th nvironmntal ) (19) tmpratur i ontant. Hat onution an b rib in onrt, at lat for tmpratur not xing Equation 100 C (17) ontitut (Bažant th & lmnt Kaplan noal 1996), for by vtor Fourir q la, uh hih that ra may rit δw I δv T q. Th iffrntial form of th intrnal virtual ork ra q λ T (7) I T T r W δ δε σ δ S u σ (20) hr q i th hat S flux, T i th abolut tmpratur, an λ i th hat onutivity; in thi Proing of FraMCoS-7, May 23-28, 2010

4 In an lmnt thi may b xpr olly in trm of th uual (ontinuou) DOF vtor I T T T T W δ δv B Z B D B B Z T rt r r S ZT DT Z S V V : (21) Thi Equation ontitut th lmnt tiffn matrix kt uh that may rit δw I δv T ktv. Th tangnt tiffn matrix i only ymmtri if th r matrix D i ymmtri, hih i not narily r th a. Hovr, for itrativ purpo D may r rt b rpla by th ymmtri matrix ½ D D. 4 NUMERICS On th global lvl th virtual ork Equation furnih th irt quilibrium Equation: QV ( ) R (22) hr Q i th um of lmnt noal for tablih through Equation (17)-(19), V i th global DOF vtor an R i th global noal loa vtor trmin from (8). Equation (22) i a t of nonlinar Equation hih ar olv itrativly applying th linar inrmntal rlation K V R (23) T hr K T i th global tangnt tiffn matrix ambl from lmnt matri tablih through (21), V i th global inrmntal DOF vtor an R i global inrmntal loa vtor. Th lmnt prnt in th prviou tion i ba on th CST, hovr, it allo for th formation of a iplamnt iontinuity or a rak. Thu hav nam th lmnt CST. Th CST ha thr no an ix DOF, to at ah no ribing th iplamnt vtor. Th atual valu of th iontinuity vtor i alulat at th lmnt lvl an no global DOF ar n to rprnt th vtor lmnt. A rak i form if th prinipl tr in an lmnt x th uni-axial tnil trngth, an th normal to th iontinuity lin i paralll to th prinipl tr vtor. If a nighboring lmnt i in th rak tat, th rak in th atual lmnt i for to onnt to th nighboring rak; othri it i for to pa through th ntr of th lmnt. Th nonlinar quilibrium Equation may b olv by tanar FEM prour; hr hav appli th orthogonal riual algorithm (Krnk 1995). Convrgn i nur at vry loa tp in trm of an nrgy J ritrion D ( h, T ) rlat h to th lati nrgy of th initial lati loa tp. Th proportionality offiint D(h,T) moitur prmability an it i a nonlina 5 EXAMPLE of th rlativ humiity h an tmpratur & Najjar 1972). Th moitur ma balan Th apabiliti that of th th CST variation lmnt in tim ar of montrat through th volum moling of onrt of th (atr RILEM ontnt thr-) b q th atr ma point bning tt ivrgn (Vanrvall of th 2000) moitur for trmining fratur mhanial proprti of onrt, flux J hih ha bom a bnhmark tt alo for th numrial moling of ohiv J fratur propagation. t A i vi of th tt tup i hon in Figur 3. Th bam ha a quar Th atr ro-tion ontnt an an a 25 b xpr mm a p noth ut prpniular of th vaporabl to th atr bottom fa in (apillary a th mi-tion plan; vapor, it an i imply aorb upport atr) an at th th non- n an loa (hmially by a loa ating boun) vrtially atr onar at mi-pan. n (Mil Pantazopoulo & Mill 1995). It i ra aum that th vaporabl atr i a fu rlativ humiity, h, gr of hyration gr of ilia fum ration,, i.. ag-pnnt orption/orption (Norling Mjonll 1997). Unr thi aum by ubtituting Equation 1 into Equati obtain Figur 3. Gomtry of th RILEM tt ( D bam h) ith a & 25 mm & h noth. Th bam i aum hr to / b at i th in lop onrt of th ith orption/ haratriti aoring iothrm to Tabl (alo 1. all Th tnil moitur oftning of th onrt govrning i aum quation to (Equation oby a linar 3) mut b apa ohiv la a by pit appropriat in Figur bounary 4. Th an initial rak oniti propagation in thi tt Th i rlation ominat btn by th th opning amount of of th rak, hn atr th an mix rlativ mo humiity haratriti i all of th onrt ar iothrm not onir. if maur Hovr, ith liing inraing of th rak fa humiity mut prou an orption om tr iothrm in orr in th to avoi unrtrit a. Nglting rigi boy thir motion iffrn of th (Xi t al. rak lmnt th part folloing, hih i only orption attah iothrm to on ill b no. Thrfor, rfrn liing i to mol both orption latially an orption ith an infrior tiffn. By No th intration ay, if btn th hytri th to of th rak formation iothrm i onir. oul b takn into aount, to rlation, vaporabl atr v rlativ humi Tabl 1. Conrt matrial paramtr. b u aoring to th ign of th varia Paramtr rlativity humiity. Valu Th hap of th Young moulu, Eiothrm for HPC i 37.4 influn GPa by many p pially tho that influn xtnt an Poion ratio, ν 0.2 hmial ration an, in turn, trm Tnil trngth, f t trutur an por iz 3.5 MPa itribution (atrratio, mnt hmial 160 J/mompoition, 2 SF Fratur nrgy, G f uring tim an mtho, tmpratur, mix Th CST lmnt t.). i In tt th litratur againt th variou prforman of an intrfa foun lmnt to rib uppli th ith orption th om- iothrm formulatio mrially availabl onrt FEM (Xi o t al. DIANA 1994). Hovr, (Diana in th 2003). In th DIANA papr rfrn th mi-mpirial mol th xprion bulk of pro th bam i mol Norling ith Mjornll CST lmnt, (1997) an i intr- aopt b Proing of FraMCoS-7, May 23-28, 2010

5 fa J D ( lmnt h, T ) h ar pr-loat to mol a rak path (1) onfin to th mi-pan tion of th bam. Th proportionality offiint D(h,T) i all moitur prmability an it i a nonlinar funtion of th rlativ humiity h an tmpratur T (Bažant & Najjar 1972). Th moitur ma balan rquir that th variation in tim of th atr ma pr unit volum of onrt (atr ontnt ) b qual to th ivrgn of th moitur flux J Figur 4. Linar J tnion oftning urv. t (2) Th atr ontnt an b xpr a th um of th vaporabl atr (apillary atr, atr vapor, an aorb atr) an th non-vaporabl (hmially boun) atr n (Mill 1966, Pantazopoulo & Mill 1995). It i raonabl to aum that th vaporabl atr i a funtion of rlativ humiity, h, gr of hyration,, an gr of ilia fum ration,, i.. (h,, ) ag-pnnt orption/orption iothrm (Norling Mjonll 1997). Unr thi aumption an by ubtituting Equation 1 into Equation 2 on obtain Figur 5. Loa-fltion urv of th RILEM tt bam ith ( D h) & & & n (3) a 25 mm h noth. Comparion h btn a DIANA mol ith 48 intrfa lmnt ovr th ligamnt hight, an to CST mol ith untrutur mh having 14 an 28 lmnt ovr hr th ligamnt / i hight, th lop rptivly. of th orption/orption iothrm (alo all moitur apaity). Th govrning To mak quation a fair omparion (Equation 3) btn mut omplt th CST mol by appropriat an th bounary rfrn an mol, initial th onition. CST mol i allo Th rlation to vlop btn on rak th amount only, although of vaporabl it ha th atr apability an rlativ of moling humiity multi i rak all vlopmnt. aorption Th iothrm on-rak if maur rtrittion ith i obtain inraing by only rlativity alloing humiity an lmnt an orption to rak if it iothrm ha a nighbor in th hih oppoit i rak. a. Nglting Hovr, thir to gt iffrn tarting, (Xi th t firt al. 1994), lmnt in to th rak folloing, i xmpt orption from iothrm thi rul. ill b u ith rfrn Th rfrn to both orption bam a an mol orption in onition. a rgular mh By th ith ay, 48 lmnt if th ovr hytri th hight of th of th moitur bam. Untrutur iothrm oul mh b takn r into u aount, for to th iffrnt CSTmol rlation, ith vaporabl to niti atr v furnihing rlativ humiity, 14 an 28 mut lmnt b u ovr aoring th ligamnt to th hight, ign of rptivly. th variation of th rlativity Simulat humiity. loa-fltion Th hap urv of ar th hon orption in Figur iothrm 5. for Th HPC fltion i influn i maur by many at paramtr, mi-pan rlativ pially to th tho man that valu influn of th xtnt formation an rat of th bam hmial mi-plan ration at th an, upport. in turn, Rult trmin for both por mh trutur ompar an por nily iz itribution ith th rult (atr-to-mnt for th rfrn ratio, mol, mnt although hmial th CST ompoition, mol lightly SF ontnt, ovrtimat uring tim th an pak mtho, loa. tmpratur, It houl b mix mphaiz, aitiv, though, t.). In that th th litratur rak path variou ha formulation not bn prfin an b in foun th CST to rib mol. th In orption Figur 6 iothrm lmnt of mh normal lo onrt to th (Xi noth t al. an 1994). rak Hovr, path at in pak th loa prnt ar hon, papr th rptivly mi-mpirial for th xprion to mh niti propo appli. Norling Only Mjornll th mil (1997) part i of aopt th bam bau (±25mm by it xpliitly aount for th volution of hyration ration an SF ontnt. Thi orption iothrm ra 1 ( h,, ) G (, ) ( g ) h 1 10( g ) h K (, ) (4) hr th firt trm (gl iothrm) rprnt th phyially boun (aorb) atr an th on trm (apillary iothrm) rprnt th apillary atr. Thi xprion i vali only for lo ontnt of SF. Th offiint G 1 rprnt th amount of atr pr unit volum hl in th gl por at 100% rlativ humiity, an it an b xpr (Norling Mjornll 1997) a G (, ) k k 1 vg vg (5) hr k vg an k vg ar matrial paramtr. From th Figur maximum 6. Elmnt amount mh of for atr CST pr mol unit lo volum to th that noth, an an fill all rak por path (both at pak apillary loa. Th por mil an part gl of por), th bam on (±25mm an alulat from th K 1 mi-tion) a on obtain i hon for to mh niti. from th mi-tion) i hon. Th 10 gining h 1 of th CST rak 0 i of our G a 1 onqun of th oar (, ontant ) train fil of th lmnt. Hovr, th ovrall irtion 10 g of th h rak i in lin (6) K 1 1 ith th mi-tion plan. 1 Th matrial paramtr k vg an k vg an g 1 an 6 b CONCLUSION alibrat by fitting xprimntal ata rlvant to fr (vaporabl) atr ontnt in onrt at A variou impl ag lmnt (Di Luzio all & Cuati CST ith 2009b). an mb iontinuity for moling onrt raking ha bn 2.2 Tmpratur prnt. Th volution lmnt iffr from prviou trong iontinuity lmnt of th mb or XFEM Not that, typ, at arly in that ag, th in tr th fil hmial i onitnt ration ithin aoiat th lmnt ith mnt nuring hyration tr ontinuity an SF ration aro th ar rak. xothrmi, Th irt th tmpratur nformnt fil of tr i not ompatibility for non-aiabati ha allo ytm for th limination vn if th of nvironmntal xtra DOF uniform for tmpratur ribing i th ontant. rak opning Hat an onution liing, thu an th b CST rib ha th in am onrt, numbr at lat of DOF for tmpratur a th CST. not xing Th goo 100 C prforman (Bažant of & th Kaplan CST 1996), ha bn by montrat Fourir la, in hih th ra Mo I bnhmark tt: th thr-point bning of a noth onrt bam. q λ T Th avantag of th CST lmnt i th impliity of it. Implmntation i traightforar an (7) th hr hanling q i of th a multipl hat flux, rak T ituation i th ill abolut not rquir tmpratur, furthr an vlopmnt. λ i th hat onutivity; in thi Proing of FraMCoS-7, May 23-28, 2010

6 REFERENCES Afrg, J.L., Pouln, P.N. an Niln, L.O A onitnt partly rak XFEM lmnt for ohiv rak groth. Int. J. Numr. Mth. Engng. Vol. 72: Blythko, T. an Blak, T Elati rak groth in finit lmnt ith minimal rmhing. Int. J. Numr. Mth. Engng. Vol. 45: Diana Ur Manual DIANA Ur Manual, Elmnt Library, (Eition 8.1). TNO Builing an Contrution Rarh, Dlft, Th Nthrlan. Dvorkin, E.N., Cuitiño, A.M. an Gioia, G Finit lmnt ith iplamnt intrpolat mb loalization lin innitiv to mh iz an itortion. Comp. Mth. in Appl. Mh. an Engng. Vol. 90: Krnk, S An orthogonal riual prour for nonlinar finit lmnt quation. Int. J. Numr. Mth. Engng. Vol. 38: Moë, N., Dolbo, J. J an DBlythko, ( h, T ) h T A finit lmnt mtho for rak groth ithout rmhing. Int. J. Numr. Mth. Engng. Vol. 46: Moë, N., an Blythko, Th T. proportionality Extn finit offiint lmnt D(h,T) mtho for ohiv moitur rak prmability groth. Engng. an Frat. it i Mh. a nonlina Vol. 63: of th rlativ humiity h an tmpratur Mougaar, J.F., Pouln, & Najjar P.N. an 1972). Niln, Th L.O. moitur A partly ma balan an fully rak XFEM lmnt ba on highr orr that th variation in tim of th atr ma polynomial hap funtion for moling ohiv fratur. Submitt for publiation. volum of onrt (atr ontnt ) b q Olivr, J Molling ivrgn trong of iontinuiti th moitur in flux oli J mhani via train oftning ontitutiv quation. Part 1: Funamntal. Part 2: Numrial imulation. Int. J. Numr. Mth. Engng. Vol. 39: J Vanrvall, L Tt t an ign mtho for fibr rinfor onrt. Rommnation for bning tt. Matr. & Strut. Vol. 33: 3-5. Th atr ontnt an b xpr a of th vaporabl atr (apillary a vapor, an aorb atr) an th non- (hmially boun) atr n (Mil Pantazopoulo & Mill 1995). It i ra aum that th vaporabl atr i a fu rlativ humiity, h, gr of hyration gr of ilia fum ration,, i.. ag-pnnt orption/orption (Norling Mjonll 1997). Unr thi aum by ubtituting Equation 1 into Equati obtain ( D h) h & & hr / i th lop of th orption/ iothrm (alo all moitur apa govrning quation (Equation 3) mut b by appropriat bounary an initial oniti Th rlation btn th amount of atr an rlativ humiity i all iothrm if maur ith inraing humiity an orption iothrm in th a. Nglting thir iffrn (Xi t al. th folloing, orption iothrm ill b rfrn to both orption an orption By th ay, if th hytri of th iothrm oul b takn into aount, to rlation, vaporabl atr v rlativ humi b u aoring to th ign of th varia rlativity humiity. Th hap of th iothrm for HPC i influn by many p pially tho that influn xtnt an hmial ration an, in turn, trm trutur an por iz itribution (atrratio, mnt hmial ompoition, SF uring tim an mtho, tmpratur, mix t.). In th litratur variou formulatio foun to rib th orption iothrm onrt (Xi t al. 1994). Hovr, in th papr th mi-mpirial xprion pro Norling Mjornll (1997) i aopt b Proing of FraMCoS-7, May 23-28, 2010

Fracture simulation of fiber reinforced concrete by visco-elasto-plastic suspension element method

Fracture simulation of fiber reinforced concrete by visco-elasto-plastic suspension element method Fratur Mhani Conrt Conrt Strutur - High Prforman, Fibr Rinford Conrt, Spial Loading Strutural Appliation- B. H. Oh, t al. (d) 200 Kora Conrt Intitut, ISBN 978-89-5708-82-2 Fratur imulation fibr rinford

More information

Fracture properties of high-strength steel fiber concrete

Fracture properties of high-strength steel fiber concrete Fratur Mhani of Conrt and Conrt Strutur - High Prforman, Fibr Rinford Conrt, Spial Loading and Strutural Appliation- B. H. Oh, t al. (d) 200 Kora Conrt Intitut, ISBN 978-89-5708-82-2 Fratur proprti of

More information

Verification of wet and dry packing methods with experimental data

Verification of wet and dry packing methods with experimental data Fratur Mhani Conrt Conrt Strutur - High Prforman, Fibr Rinfor Conrt, Spial Loaing Strutural Appliation- B. H. Oh, t al. () Kora Conrt Intitut, ISBN 978-89-578-8- Vrifiation wt ry paking with xprimntal

More information

Behaviors of FRP sheet reinforced concrete to impact and static loading

Behaviors of FRP sheet reinforced concrete to impact and static loading Fratur Mhani of Conrt Conrt Strutur - High Prforman, Fibr Rinford Conrt, Spial Loading Strutural Appliation- B. H. Oh, t al. (d) 21 Kora Conrt Intitut, ISBN 978-89-578-182-2 Bhavior of FRP ht rinford onrt

More information

A local bond stress-slip model for reinforcing bars in self-compacting concrete

A local bond stress-slip model for reinforcing bars in self-compacting concrete Fratur Mhani of Conrt Conrt Strutur - Amnt, Durability, Monitoring Rtrofitting of Conrt Strutur- B. H. Oh, t al. (d) 200 Kora Conrt Intitut, Soul, ISBN 978-89-5708-8-5 A loal bond tr-lip modl for rinforing

More information

Chloride diffusion in the cracked concrete

Chloride diffusion in the cracked concrete Fratur Mhani of Conrt and Conrt Strutur - Amnt, Durability, Monitoring and Rtrofitting of Conrt Strutur- B. H. Oh, t al. (d) Kora Conrt Intitut, Soul, ISBN 97-9-57--5 Chlorid diffuion in th rakd onrt W.L.

More information

Bond analysis model of deformed bars to concrete

Bond analysis model of deformed bars to concrete Fratur Mhani of Conrt Conrt Strutur - Amnt, Durability, Monitoring Rtrofitting of Conrt Strutur- B. H. Oh, t al. (d) Kora Conrt Intitut, Soul, ISBN 978-89-578-8-5 Bond analyi modl of dformd bar to onrt

More information

Fiber reinforced concrete characterization through round panel test - Part II: analytical and numerical study

Fiber reinforced concrete characterization through round panel test - Part II: analytical and numerical study Fratur Mhani Conrt Conrt Strutur - High rforman Fibr Rinford Conrt Spial Loading Strutural Appliation- B. H. Oh t al. (d) Kora Conrt Intitut ISBN 978-89-578-8- Fibr rinford onrt haratrization through round

More information

Determination of fracture parameters of concrete interfaces using DIC

Determination of fracture parameters of concrete interfaces using DIC Fratur Mhani of Conrt Conrt Strutur - Amnt, Durability, Monitoring Rtrofitting of Conrt Strutur- B. H. Oh, t al. (d) 2010 Kora Conrt Intitut, Soul, ISBN 978-89-5708-181-5 Dtrmination of fratur paramtr

More information

Engineered cementitious composites with low volume of cementitious materials

Engineered cementitious composites with low volume of cementitious materials Fratur Mhani Conrt Conrt Strutur - High Prforman, Fibr Rinford Conrt, Spial Loading Strutural Appliation- B. H. Oh, t al. (d) 2 Kora Conrt Intitut, ISBN 978-89-578-82-2 Enginrd mntitiou ompoit ith lo volum

More information

An analytical study on the stress-strain relation of PVA-ECC under tensile fatigue

An analytical study on the stress-strain relation of PVA-ECC under tensile fatigue Fratur Mhani Conrt Conrt Strutur - High Prforman, Fibr Rinford Conrt, Spial Loading Strutural Appliation- B. H. Oh, t al. (d) Kora Conrt Intitut, ISBN 978-89-578-8- An analytial tudy on th tr-train rlation

More information

Fiber reinforced concrete characterization through round panel test - part I: experimental study

Fiber reinforced concrete characterization through round panel test - part I: experimental study Fratur Mhani of Conrt and Conrt Strutur - High Prforman, Fibr Rinford Conrt, Spial Loading and Strutural Appliation- B. H. Oh, t al. (d) 2 Kora Conrt Intitut, ISBN 978-89-78-82-2 Fibr rinford onrt haratrization

More information

Crack width control of reinforced concrete one-way slabs utilizing expansive strain-hardening cement-based composites (SHCCs)

Crack width control of reinforced concrete one-way slabs utilizing expansive strain-hardening cement-based composites (SHCCs) Fratur Mhani Conrt Conrt Strutur - High Prforman, Fibr Rinford Conrt, Spial Loading Strutural Appliation- B. H. Oh, t al. (d) 2 Kora Conrt Intitut, ISBN 978-89-578-82-2 Crak width ontrol rinford onrt on-way

More information

Measuring crack width and spacing in reinforced concrete members

Measuring crack width and spacing in reinforced concrete members Fratur Mhani Conrt and Conrt Strutur - Rnt Advan in Fratur Mhani Conrt - B. H. Oh, t al.(d) 200 Kora Conrt Intitut, Soul, ISBN 978-89-5708-80-8 Mauring rak idth and paing in rinford onrt mmbr S. Yair Alam,

More information

Fuzzy Logic Model of Fiber Concrete

Fuzzy Logic Model of Fiber Concrete Fratur Mhani Conrt Conrt Strutur - High Prforman, Fibr Rinford Conrt, Spial Loading Strutural Appliation- B. H. Oh, t al. (d) 2 Kora Conrt Intitut, ISBN 978-89-78-82-2 Fuzzy Logi Modl Fibr Conrt A. Kohoutková,

More information

Static and fatigue failure simulation of concrete material by discrete analysis

Static and fatigue failure simulation of concrete material by discrete analysis Fratur Mhani Conrt Conrt Strutur - Rnt Advan in Fratur Mhani Conrt - B. H. Oh, t al.(d) 2010 Kora Conrt Intitut, Soul, ISBN 978-89-5708-180-8 Stati fatigu failur imulation onrt matrial dirt analyi K. Nagai

More information

Experimental investigation of compressive concrete elements confined with shape memory Ni-Ti wires

Experimental investigation of compressive concrete elements confined with shape memory Ni-Ti wires Fratur Mhani Conrt Conrt Strutur - Amnt, Durability, Monitoring Rtritting Conrt Strutur- B. H. Oh, t al. (d) 200 Kora Conrt Intitut, Soul, ISBN 978-89-5708-8-5 Exprimntal invtigation ompriv onrt lmnt onfind

More information

Experimental study on the flexural behaviour of fibre reinforced concretes strengthened with steel and macro-synthetic fibres

Experimental study on the flexural behaviour of fibre reinforced concretes strengthened with steel and macro-synthetic fibres Fratur Mhani o Conrt Conrt Strutur - Amnt, Durability, Monitorg Rtroittg o Conrt Strutur- B. H. Oh, t al. (d) 2 Kora Conrt Intitut, Soul, ISBN 978-89-578-8-5 Exprimntal tudy on th lxural bhaviour o ibr

More information

Behavior of concrete members constructed with SHCC/GFRP permanent formwork

Behavior of concrete members constructed with SHCC/GFRP permanent formwork Fratur Mhani Conrt Conrt Strutur - High Prforman, Fibr Rinford Conrt, Spial Loading Strutural Appliation- B. H. Oh, t al. (d) 200 Kora Conrt Intitut, ISBN 978-89-5708-82-2 Bhavior onrt mmbr ontrutd ith

More information

O. Omikrine-Metalssi & V.-D. Le Université Paris-Est, Paris, France

O. Omikrine-Metalssi & V.-D. Le Université Paris-Est, Paris, France Fratur Mhani Conrt Conrt Strutur - Rnt Advan in Fratur Mhani Conrt - B. H. Oh, t al.(d) 2 Kora Conrt Intitut, Soul, ISBN 978-89-578-8-8 Intgration ontat lmnt in RGIB-modul th finit lmnt tar CESAR-LCPC

More information

Blast loading response of ultra high performance concrete and reactive powder concrete slabs

Blast loading response of ultra high performance concrete and reactive powder concrete slabs Fratur Mhani of Conrt and Conrt Strutur - High Prforman, Fibr Rinford Conrt, Spial Loading and Strutural Appliation- B. H. Oh, t al. (d) 2 Kora Conrt Intitut, ISBN 978-89-578-82-2 Blat loading rpon of

More information

Durability performance of UFC sakata-mira footbridge under sea environment

Durability performance of UFC sakata-mira footbridge under sea environment Fratur Mhani of Conrt and Conrt Strutur - High Prforman, Fibr Rinford Conrt, Spial Loading and Strutural Appliation- B. H. Oh, t al. (d) 2 Kora Conrt Intitut, ISBN 978-89-578-82-2 Durability prforman of

More information

Fracture energy of high performance mortar subjected to high temperature

Fracture energy of high performance mortar subjected to high temperature Fratur Mhani Conrt Conrt Strutur - Rnt Advan in Fratur Mhani Conrt - B. H. Oh, t al.(d) 2 Kora Conrt Intitut, Soul, ISBN 978-89-578-8-8 Fratur nrgy high prforman mortar ubjtd high tmpratur S. Djaknoun

More information

Toughness indices of fiber reinforced concrete subjected to mode II loading

Toughness indices of fiber reinforced concrete subjected to mode II loading Fratur Mhani of Conrt Conrt Strutur - Rnt Advan in Fratur Mhani of Conrt - B. H. Oh, t al.(d) 200 Kora Conrt Intitut, Soul, ISBN 978-89-5708-80-8 Toughn indi of fibr rinford onrt ubjtd to mod II loading

More information

Cracking analysis of brick masonry arch bridge

Cracking analysis of brick masonry arch bridge Fratur Mhani of Conrt and Conrt Strutur - High Prforman, Fibr Rinford Conrt, Spial Loading and Strutural Appliation- B. H. Oh, t al. (d) 200 Kora Conrt Intitut, ISBN 978-89-5708-82-2 Craking analyi of

More information

Crack formation and tensile stress-crack opening behavior of fiber reinforced cementitious composites (FRCC)

Crack formation and tensile stress-crack opening behavior of fiber reinforced cementitious composites (FRCC) Fratur Mhani Conrt Conrt Strutur - High Prforman, Fibr Rinford Conrt, Spial Loading Strutural Appliation- B. H. Oh, t al. (d) 200 Kora Conrt Intitut, ISBN 978-89-5708-82-2 Crak formation tnil tr-rak opning

More information

Effect of loading condition, specimen geometry, size-effect and softening function on double-k fracture parameters of concrete

Effect of loading condition, specimen geometry, size-effect and softening function on double-k fracture parameters of concrete Fratur Mhani onrt onrt Strutur - Rnt Advan Fratur Mhani onrt - B. H. Oh, t al.(d) 2010 Kora onrt Intitut, Soul, ISBN 978-89-5708-180-8 Efft loadg ondition, pimn gomtry, iz-fft tng ftion on doubl-k fratur

More information

Simulation of tensile performance of fiber reinforced cementitious composite with fracture mechanics model

Simulation of tensile performance of fiber reinforced cementitious composite with fracture mechanics model Fratur Mhani Conrt Conrt Strutur - High Prorman, Fibr Rinord Conrt, Spial Loading Strutural Appliation- B. H. Oh, t al. (d) Kora Conrt Intitut, ISBN 978-89-578-8- Simulation tnil prorman ibr rinord mntitiou

More information

Study of the effect of alkali-silica reaction on properties of concrete by means of conventional test methods and non-destructive test methods

Study of the effect of alkali-silica reaction on properties of concrete by means of conventional test methods and non-destructive test methods Fratur Mhani of Conrt and Conrt Strutur - Amnt, Durability, Monitoring and Rtrofitting of Conrt Strutur- B. H. Oh, t al. (d) 2010 Kora Conrt Intitut, Soul, ISBN 978-89-5708-181-5 Study of th fft of alkali-ilia

More information

Effect of short fibres on fracture behaviour of textile reinforced concrete

Effect of short fibres on fracture behaviour of textile reinforced concrete Fratur Mhani Conrt and Conrt Strutur - High Prforman, Fibr Rinford Conrt, Spial Loading and Strutural Appliation- B. H. Oh, t al. (d) 2010 Kora Conrt Intitut, ISBN 978-89-5708-182-2 Efft hort fibr on fratur

More information

Degradation of reinforced concrete structures under atmospheric corrosion

Degradation of reinforced concrete structures under atmospheric corrosion Fratur Mhani Conrt and Conrt Strutur - Amnt, Durability, Monitoring and Rtritting Conrt Strutur- B. H. Oh, t al. (d) 200 Kora Conrt Intitut, Soul, ISBN 978-89-5708-8-5 Dgradation rinford onrt trutur undr

More information

Analysis of balanced double-lap joints with a bi-linear softening adhesive

Analysis of balanced double-lap joints with a bi-linear softening adhesive Fratur Mhani Conrt and Conrt Strutur - Ant, Durability, Monitoring and Rtritting Conrt Strutur- B. H. Oh, t al. (d 00 Kora Conrt ntitut, Soul, SBN 978-89-5708-8-5 Analyi baland doubl-lap joint with a bi-linar

More information

Pre and post-cracking behavior of steel-concrete composite deck subjected to high cycle load

Pre and post-cracking behavior of steel-concrete composite deck subjected to high cycle load Fratr Mhani Conrt Conrt Strtr - High Prforman, Fibr Rinford Conrt, Spial Loading Strtral Appliation- B. H. Oh, t al. (d) 2 Kora Conrt Intitt, ISBN 978-89-578-82-2 Pr pot-raking bhavior tl-onrt ompoit dk

More information

Size-scale effects on minimum flexural reinforcement in RC beams: application of the cohesive crack model

Size-scale effects on minimum flexural reinforcement in RC beams: application of the cohesive crack model Fratur Mhani Conrt Conrt Strutur - Rnt Advan Fratur Mhani Conrt - B. H. Oh, t al.(d) 2010 Kora Conrt Intitut, Soul, ISBN 978-89-5708-180-8 Siz-al fft on mimum flxural rformnt RC am: appliation ohiv rak

More information

Quantified estimation of rebar corrosion by means of acoustic emission technique

Quantified estimation of rebar corrosion by means of acoustic emission technique Fratur Mani Conrt Conrt Strutur - Amnt Durability Monitoring Rtritting Conrt Strutur- B. H. O t al. (d) 200 Kora Conrt Intitut Soul ISBN 978-89-5708-8-5 Quantifid timation rbar orroion by man aouti mion

More information

Influence of temperature and composition upon drying of concretes

Influence of temperature and composition upon drying of concretes Fratur Mhani Conrt Conrt Strutur - Amnt, Durability, Monirg Rtrittg Conrt Strutur- B. H. Oh, t al. (d) 2 Kora Conrt Intitut, Soul, ISBN 978-89-578-8-5 Influn tmpratur ompoition upon dryg onrt F. Bru, Y.

More information

Cover cracking in RC columns subjected to reinforcement corrosion under sustained load

Cover cracking in RC columns subjected to reinforcement corrosion under sustained load Fratur Mhani Conrt Conrt Strutur - Amnt, Durability, Monirg Rtrittg Conrt Strutur- B. H. Oh, t al. (d) 200 Kora Conrt Intitut, Soul, ISBN 978-89-5708-8-5 Covr rakg RC olumn ubjtd rformnt orroion undr utad

More information

Solutions to Supplementary Problems

Solutions to Supplementary Problems Solution to Supplmntary Problm Chaptr 5 Solution 5. Failur of th tiff clay occur, hn th ffctiv prur at th bottom of th layr bcom ro. Initially Total ovrburn prur at X : = 9 5 + = 7 kn/m Por atr prur at

More information

Relating tensile properties with flexural properties in SHCC

Relating tensile properties with flexural properties in SHCC Fratr Mhani of Conrt and Conrt Strtr - High Prforman, Fibr Rinford Conrt, Spial Loading and Strtral Appliation- B. H. Oh, t al. (d) Kora Conrt Intitt, ISBN 978-89-78-8- Rlating tnil proprti ith flxral

More information

Multiple Short Term Infusion Homework # 5 PHA 5127

Multiple Short Term Infusion Homework # 5 PHA 5127 Multipl Short rm Infusion Homwork # 5 PHA 527 A rug is aministr as a short trm infusion. h avrag pharmacokintic paramtrs for this rug ar: k 0.40 hr - V 28 L his rug follows a on-compartmnt boy mol. A 300

More information

Rebar bond slip in diagonal tension failure of reinforced concrete beams

Rebar bond slip in diagonal tension failure of reinforced concrete beams Rbar bond lip in diagonal nion failur of rinford onr bam T. Hagaa Iniu of Tnology Simizu Corporaion Tokyo apan Rbar bond lip in diagonal nion failur of rinford onr bam T. Hagaa Iniu of Tnology Simizu Corporaion

More information

Carriers Concentration in Semiconductors - VI. Prof.P. Ravindran, Department of Physics, Central University of Tamil Nadu, India

Carriers Concentration in Semiconductors - VI. Prof.P. Ravindran, Department of Physics, Central University of Tamil Nadu, India Carrirs Conntration in Smionutors - VI 1 Prof.P. Ravinran, Dpartmnt of Pysis, Cntral Univrsity of Tamil au, Inia ttp://folk.uio.no/ravi/smi01 P.Ravinran, PHY0 Smionutor Pysis, 17 January 014 : Carrirs

More information

A Hybrid Model of Rubber Elasticity in Simple Extension

A Hybrid Model of Rubber Elasticity in Simple Extension A Hybri Mol of Rubbr Elatiity in Simpl Extnion Arkay I. Lonov Dpartmnt of Polymr Enginring, Th Univrity of Akron, Akron, OH 4435-31 E-mail ar: lonov@uakron.u (A.I. Lonov) Abtrat A thrmoynamially rlat mol

More information

MA1506 Tutorial 2 Solutions. Question 1. (1a) 1 ) y x. e x. 1 exp (in general, Integrating factor is. ye dx. So ) (1b) e e. e c.

MA1506 Tutorial 2 Solutions. Question 1. (1a) 1 ) y x. e x. 1 exp (in general, Integrating factor is. ye dx. So ) (1b) e e. e c. MA56 utorial Solutions Qustion a Intgrating fator is ln p p in gnral, multipl b p So b ln p p sin his kin is all a Brnoulli quation -- st Sin w fin Y, Y Y, Y Y p Qustion Dfin v / hn our quation is v μ

More information

WEEK 3 Effective Stress and Pore Water Pressure Changes

WEEK 3 Effective Stress and Pore Water Pressure Changes WEEK 3 Effctiv Str and Por Watr Prur Chang 5. Effctiv tr ath undr undraind condition 5-1. Dfinition of ffctiv tr: A rvi A you mut hav larnt that th ffctiv tr, σ, in oil i dfind a σ σ u Whr σ i th total

More information

Punching of flat slabs: Design example

Punching of flat slabs: Design example fi Mol Co 010 Pnhing of flat la: Dign apl Stfan Lip, Arlio Mttoni, Migl Frnánz Riz Eol Polthniq Fééral Laann, Switzrlan, 0.11.010 Lip / Mttoni / Frnánz Riz / Eol Polthniq Fééral Laann, Switzrlan 1 1 Bai

More information

Characteristic Equations and Boundary Conditions

Characteristic Equations and Boundary Conditions Charatriti Equation and Boundary Condition Øytin Li-Svndn, Viggo H. Hantn, & Andrw MMurry Novmbr 4, Introdution On of th mot diffiult problm on i onfrontd with In numrial modlling oftn li in tting th boundary

More information

Case Study Vancomycin Answers Provided by Jeffrey Stark, Graduate Student

Case Study Vancomycin Answers Provided by Jeffrey Stark, Graduate Student Cas Stuy Vancomycin Answrs Provi by Jffry Stark, Grauat Stunt h antibiotic Vancomycin is liminat almost ntirly by glomrular filtration. For a patint with normal rnal function, th half-lif is about 6 hours.

More information

Numerical methods, Mixed exercise 10

Numerical methods, Mixed exercise 10 Numrial mthos, Mi ris a f ( ) 6 f ( ) 6 6 6 a = 6, b = f ( ) So. 6 b n a n 6 7.67... 6.99....67... 6.68....99... 6.667....68... To.p., th valus ar =.68, =.99, =.68, =.67. f (.6).6 6.6... f (.6).6 6.6.7...

More information

Performance prediction of eccentrically loaded RC columns wrapped with FRP

Performance prediction of eccentrically loaded RC columns wrapped with FRP Fourt Intrnational Conrn on FR Compoit in Civil Enginring (CICE008-4Jul 008, Zuri, Switzrlan rorman prition o ntriall loa RC olumn wrapp wit FR T. El Maaaw Unit rab Emirat Univrit, l-in, bu Dabi, Unit

More information

Experimental study on the ultimate strength of R/C curved beam

Experimental study on the ultimate strength of R/C curved beam Fraur Mhani Conr and Conr Sruur - High Prforman, Fibr Rinford Conr, Spial Loading and Sruural Appliaion- B. H. Oh, al. (d) 2 Kora Conr Iniu, ISBN 978-89-578-82-2 Exprimnal udy on h ulima rngh R/C urvd

More information

Digital Signal Processing, Fall 2006

Digital Signal Processing, Fall 2006 Digital Signal Prossing, Fall 006 Ltur 7: Filtr Dsign Zhng-ua an Dpartmnt of Eltroni Systms Aalborg Univrsity, Dnmar t@om.aau. Cours at a glan MM Disrt-tim signals an systms Systm MM Fourir-omain rprsntation

More information

INTRODUCTION TO AUTOMATIC CONTROLS INDEX LAPLACE TRANSFORMS

INTRODUCTION TO AUTOMATIC CONTROLS INDEX LAPLACE TRANSFORMS adjoint...6 block diagram...4 clod loop ytm... 5, 0 E()...6 (t)...6 rror tady tat tracking...6 tracking...6...6 gloary... 0 impul function...3 input...5 invr Laplac tranform, INTRODUCTION TO AUTOMATIC

More information

DESIGN SPECTRA REDUCTION COEFFICIENTS FOR SYSTEMS WITH SEISMIC ENERGY DISSIPATING DEVICES LOCATED ON FIRM GROUND

DESIGN SPECTRA REDUCTION COEFFICIENTS FOR SYSTEMS WITH SEISMIC ENERGY DISSIPATING DEVICES LOCATED ON FIRM GROUND Th 14 th Worl Confrn on Earthquak Enginring DESIGN SPECTRA REDUCTION COEFFICIENTS FOR SYSTEMS WITH SEISMIC ENERGY DISSIPATING DEVICES LOCATED ON FIRM GROUND S. E. Ruiz 1, J. P. H. Toxqui 2 an J. L. Rivra

More information

Finite Element Analysis

Finite Element Analysis Finit Elmnt Analysis L4 D Shap Functions, an Gauss Quaratur FEA Formulation Dr. Wiong Wu EGR 54 Finit Elmnt Analysis Roamap for Dvlopmnt of FE Strong form: govrning DE an BCs EGR 54 Finit Elmnt Analysis

More information

2. SIMPLE SOIL PROPETIES

2. SIMPLE SOIL PROPETIES 2. SIMPLE SOIL PROPETIES 2.1 EIGHT-OLUME RELATIONSHIPS It i oft rquir of th gotchical gir to collct, claify a ivtigat oil ampl. B it for ig of fouatio or i calculatio of arthork volum, trmiatio of oil

More information

ES 240 Solid Mechanics

ES 240 Solid Mechanics ES4 Solid Mhani Fall ES 4 Solid Mhani Joot J Vlaak Shool of Enginring and Applid Sin Harvard Univrity Th not ar largly bad on our not put togthr by Prof Suo whn h taught ES 4 in 6 and on our not dvlopd

More information

Characteristics of beam-electron cloud interaction

Characteristics of beam-electron cloud interaction Charatriti of bam-ltron loud intration Tun hift and intabilit K. Ohmi KEK Int. Workhop on Two-tram Intabiliti in Partil Alrator and Storag Ring @ KEK Tukuba Japan Bam-ltron intration Bam partil ar loalid

More information

Outline. Heat Exchangers. Heat Exchangers. Compact Heat Exchangers. Compact Heat Exchangers II. Heat Exchangers April 18, ME 375 Heat Transfer 1

Outline. Heat Exchangers. Heat Exchangers. Compact Heat Exchangers. Compact Heat Exchangers II. Heat Exchangers April 18, ME 375 Heat Transfer 1 Hat Exangr April 8, 007 Hat Exangr Larry artt Manial Engrg 375 Hat ranfr April 8, 007 Outl Bai ida f at xangr Ovrall at tranfr ffiint Lg-man tmpratur diffrn mtd Efftivn NU mtd ratial nidratin Hat Exangr

More information

Hybrid Inversion technique for predicting geometrical parameters of Porous Materials

Hybrid Inversion technique for predicting geometrical parameters of Porous Materials Aoustis 8 Paris Hybri Invrsion thniqu for priting gomtrial paramtrs of Porous Matrials P. hravag, P. Bonfiglio F. Pompoli Dipartimnto i Inggnria - Univrsity of Frrara, Via aragat 1, 441 Frrara, Italy franso.pompoli@unif.it

More information

Recent advances on self healing of concrete

Recent advances on self healing of concrete Fratur Mhani Conrt Conrt Strutur - Rnt Advan Fratur Mhani Conrt - B. H. Oh, t al.(d) 21 Kora Conrt Intitut, Soul, ISBN 978-89-578-18-8 Rnt advan on lf onrt E. Shlangn, H. Jonkr, S. Qian & A. Garia Dlft

More information

Modeling with first order equations (Sect. 2.3).

Modeling with first order equations (Sect. 2.3). Moling with first orr quations (Sct. 2.3. Main xampl: Salt in a watr tank. Th xprimntal vic. Th main quations. Analysis of th mathmatical mol. Prictions for particular situations. Salt in a watr tank.

More information

Engineering Differential Equations Practice Final Exam Solutions Fall 2011

Engineering Differential Equations Practice Final Exam Solutions Fall 2011 9.6 Enginring Diffrntial Equation Practic Final Exam Solution Fall 0 Problm. (0 pt.) Solv th following initial valu problm: x y = xy, y() = 4. Thi i a linar d.. bcau y and y appar only to th firt powr.

More information

(1) Then we could wave our hands over this and it would become:

(1) Then we could wave our hands over this and it would become: MAT* K285 Spring 28 Anthony Bnoit 4/17/28 Wk 12: Laplac Tranform Rading: Kohlr & Johnon, Chaptr 5 to p. 35 HW: 5.1: 3, 7, 1*, 19 5.2: 1, 5*, 13*, 19, 45* 5.3: 1, 11*, 19 * Pla writ-up th problm natly and

More information

Nonlinear reduced Braginskii equations with ion thermal dynamics in toroidal plasma

Nonlinear reduced Braginskii equations with ion thermal dynamics in toroidal plasma Nonlinar ruc Braginkii quation with ion thrmal ynamic in toroial plama A. Zilr a) Max-Planck-Intitut für Plamaphyik, EURATOM Aoation, 85748 Garching, Grmany J. F. Drak an B. Rogr Intitut for Plama Rarch,

More information

Uncertainty in non-linear long-term behavior and buckling of. shallow concrete-filled steel tubular arches

Uncertainty in non-linear long-term behavior and buckling of. shallow concrete-filled steel tubular arches CCM14 8-3 th July, Cambridg, England Unrtainty in non-linar long-trm bhavior and bukling of shallow onrt-filld stl tubular arhs *X. Shi¹, W. Gao¹, Y.L. Pi¹ 1 Shool of Civil and Environmnt Enginring, Th

More information

COMSACO INC. NORFOLK, VA 23502

COMSACO INC. NORFOLK, VA 23502 YMOL 9. / 9. / 9. / 9. YMOL 9. / 9. OT:. THI RIG VLOP ROM MIL--/ MIL-TL-H, TYP II, L ITH VITIO OLLO:. UPO RULT O HOK TTIG, HOK MOUT (ITM, HT ) HV IR ROM 0.0 THIK TO 0.090 THIK LLO Y MIL-TL-H, PRGRPH...

More information

First order differential equation Linear equation; Method of integrating factors

First order differential equation Linear equation; Method of integrating factors First orr iffrntial quation Linar quation; Mtho of intgrating factors Exampl 1: Rwrit th lft han si as th rivativ of th prouct of y an som function by prouct rul irctly. Solving th first orr iffrntial

More information

Solution 3.1 Prove the following: γ d. (a) Start with fundamental definitions: V = (b) e = 1 n. wg e S =

Solution 3.1 Prove the following: γ d. (a) Start with fundamental definitions: V = (b) e = 1 n. wg e S = Solution. Prove the folloing: (a) G + e Start ith funamental efinition: W ; W V G ; V V V V G G V (l + e) l + e (l + e) ; ubitute for W an V (b) e n n S G e G ( n) n Solution.2 Dr D r relative enity hich

More information

A Simple Method of Tuning PI Controllers for Interval Plant of Cold Rolling Mill

A Simple Method of Tuning PI Controllers for Interval Plant of Cold Rolling Mill ntrnational Journal of Rnt Trns in Enginring, Vol. 1, No. 4, May 009 A Simpl Mtho of Tuning P Controllrs for ntrval Plant of Col Rolling Mill S.Umamahswari 1, V.Palanisamy, M.Chiambaram 3, 1 Dpartmnt of

More information

3) Use the average steady-state equation to determine the dose. Note that only 100 mg tablets of aminophylline are available here.

3) Use the average steady-state equation to determine the dose. Note that only 100 mg tablets of aminophylline are available here. PHA 5127 Dsigning A Dosing Rgimn Answrs provi by Jry Stark Mr. JM is to b start on aminophyllin or th tratmnt o asthma. H is a non-smokr an wighs 60 kg. Dsign an oral osing rgimn or this patint such that

More information

a 1and x is any real number.

a 1and x is any real number. Calcls Nots Eponnts an Logarithms Eponntial Fnction: Has th form y a, whr a 0, a an is any ral nmbr. Graph y, Graph y ln y y Th Natral Bas (Elr s nmbr): An irrational nmbr, symboliz by th lttr, appars

More information

Calculation of electromotive force induced by the slot harmonics and parameters of the linear generator

Calculation of electromotive force induced by the slot harmonics and parameters of the linear generator Calculation of lctromotiv forc inducd by th lot harmonic and paramtr of th linar gnrator (*)Hui-juan IU (**)Yi-huang ZHANG (*)School of Elctrical Enginring, Bijing Jiaotong Univrity, Bijing,China 8++58483,

More information

Finite element discretization of Laplace and Poisson equations

Finite element discretization of Laplace and Poisson equations Finit lmnt discrtization of Laplac and Poisson quations Yashwanth Tummala Tutor: Prof S.Mittal 1 Outlin Finit Elmnt Mthod for 1D Introduction to Poisson s and Laplac s Equations Finit Elmnt Mthod for 2D-Discrtization

More information

AS 5850 Finite Element Analysis

AS 5850 Finite Element Analysis AS 5850 Finit Elmnt Analysis Two-Dimnsional Linar Elasticity Instructor Prof. IIT Madras Equations of Plan Elasticity - 1 displacmnt fild strain- displacmnt rlations (infinitsimal strain) in matrix form

More information

A Quadratic Serendipity Plane Stress Rectangular Element

A Quadratic Serendipity Plane Stress Rectangular Element MAE 323: Chaptr 5 Putting It All Togthr A Quadratic Srndipity Plan Str Rctangular Elmnt In Chaptr 2, w larnd two diffrnt nrgy-bad mthod of: 1. Turning diffrntial quation into intgral (or nrgy) quation

More information

Weighted graphs -- reminder. Data Structures LECTURE 15. Shortest paths algorithms. Example: weighted graph. Two basic properties of shortest paths

Weighted graphs -- reminder. Data Structures LECTURE 15. Shortest paths algorithms. Example: weighted graph. Two basic properties of shortest paths Dt Strutur LECTURE Shortt pth lgorithm Proprti of hortt pth Bllmn-For lgorithm Dijktr lgorithm Chptr in th txtook (pp ). Wight grph -- rminr A wight grph i grph in whih g hv wight (ot) w(v i, v j ) >.

More information

The Exponential, Gaussian and Uniform Truncated Discrete Density Functions for Discrete Time Systems Analysis

The Exponential, Gaussian and Uniform Truncated Discrete Density Functions for Discrete Time Systems Analysis tvan Brbr Th Exponntial, Gauian an Uniform Truncat Dicrt Dnity Function for Dicrt Tim ytm nalyi TEVN BERBER Elctrical an Computr Enginring Dpartmnt Th Univrity of ucklan Princ trt, 00 ucklan Cntral NEW

More information

ESCI 341 Atmospheric Thermodynamics Lesson 16 Pseudoadiabatic Processes Dr. DeCaria

ESCI 341 Atmospheric Thermodynamics Lesson 16 Pseudoadiabatic Processes Dr. DeCaria ESCI 34 Atmohi hmoynami on 6 Puoaiabati Po D DCaia fn: Man, A an FE obitaill, 97: A omaion of th uialnt otntial tmatu an th tati ngy, J Atmo Si, 7, 37-39 Btt, AK, 974: Futh ommnt on A omaion of th uialnt

More information

Y 0. Standing Wave Interference between the incident & reflected waves Standing wave. A string with one end fixed on a wall

Y 0. Standing Wave Interference between the incident & reflected waves Standing wave. A string with one end fixed on a wall Staning Wav Intrfrnc btwn th incint & rflct wavs Staning wav A string with on n fix on a wall Incint: y, t) Y cos( t ) 1( Y 1 ( ) Y (St th incint wav s phas to b, i.., Y + ral & positiv.) Rflct: y, t)

More information

Equil. Properties of Reacting Gas Mixtures. So far have looked at Statistical Mechanics results for a single (pure) perfect gas

Equil. Properties of Reacting Gas Mixtures. So far have looked at Statistical Mechanics results for a single (pure) perfect gas Shool of roa Engnrng Equl. Prort of Ratng Ga Mxtur So far hav lookd at Stattal Mhan rult for a ngl (ur) rft ga hown how to gt ga rort (,, h, v,,, ) from artton funton () For nonratng rft ga mxtur, gt mxtur

More information

ELEC 372 LECTURE NOTES, WEEK 11 Dr. Amir G. Aghdam Concordia University

ELEC 372 LECTURE NOTES, WEEK 11 Dr. Amir G. Aghdam Concordia University ELEC 37 LECTURE NOTES, WEE Dr Amir Aghdam Cncrdia Univrity Part f th nt ar adaptd frm th matrial in th fllwing rfrnc: Mdrn Cntrl Sytm by Richard C Drf and Rbrt H Bihp, Prntic Hall Fdback Cntrl f Dynamic

More information

are given in the table below. t (hours)

are given in the table below. t (hours) CALCULUS WORKSHEET ON INTEGRATION WITH DATA Work th following on notbook papr. Giv dcimal answrs corrct to thr dcimal placs.. A tank contains gallons of oil at tim t = hours. Oil is bing pumpd into th

More information

2. Finite Impulse Response Filters (FIR)

2. Finite Impulse Response Filters (FIR) .. Mthos for FIR filtrs implmntation. Finit Impuls Rspons Filtrs (FIR. Th winow mtho.. Frquncy charactristic uniform sampling. 3. Maximum rror minimizing. 4. Last-squars rror minimizing.. Mthos for FIR

More information

STRIPLINES. A stripline is a planar type transmission line which is well suited for microwave integrated circuitry and photolithographic fabrication.

STRIPLINES. A stripline is a planar type transmission line which is well suited for microwave integrated circuitry and photolithographic fabrication. STIPLINES A tiplin i a plana typ tanmiion lin hih i ll uitd fo mioav intgatd iuity and photolithogaphi faiation. It i uually ontutd y thing th nt onduto of idth, on a utat of thikn and thn oving ith anoth

More information

Utilizing exact and Monte Carlo methods to investigate properties of the Blume Capel Model applied to a nine site lattice.

Utilizing exact and Monte Carlo methods to investigate properties of the Blume Capel Model applied to a nine site lattice. Utilizing xat and Mont Carlo mthods to invstigat proprtis of th Blum Capl Modl applid to a nin sit latti Nik Franios Writing various xat and Mont Carlo omputr algorithms in C languag, I usd th Blum Capl

More information

Fracture analysis of strain hardening cementitious composites by means of discrete modeling of short fibers

Fracture analysis of strain hardening cementitious composites by means of discrete modeling of short fibers Fraur Mhani Conr and Conr Sruur - Rn Advan in Fraur Mhani Conr - B. H. Oh, al.(d) 200 Kora Conr Iniu, Soul, ISBN 978-89-5708-80-8 Fraur analyi rain hardning mniiou ompoi by man dir modling hor fibr M.

More information

Case Study 4 PHA 5127 Aminoglycosides Answers provided by Jeffrey Stark Graduate Student

Case Study 4 PHA 5127 Aminoglycosides Answers provided by Jeffrey Stark Graduate Student Cas Stuy 4 PHA 527 Aminoglycosis Answrs provi by Jffry Stark Grauat Stunt Backgroun Gntamicin is us to trat a wi varity of infctions. Howvr, u to its toxicity, its us must b rstrict to th thrapy of lif-thratning

More information

Thomas Whitham Sixth Form

Thomas Whitham Sixth Form Thomas Whitham Sith Form Pur Mathmatics Unit C Algbra Trigonomtr Gomtr Calculus Vctor gomtr Pag Algbra Molus functions graphs, quations an inqualitis Graph of f () Draw f () an rflct an part of th curv

More information

A L A BA M A L A W R E V IE W

A L A BA M A L A W R E V IE W A L A BA M A L A W R E V IE W Volume 52 Fall 2000 Number 1 B E F O R E D I S A B I L I T Y C I V I L R I G HT S : C I V I L W A R P E N S I O N S A N D TH E P O L I T I C S O F D I S A B I L I T Y I N

More information

Fracture mechanics of early-age concrete

Fracture mechanics of early-age concrete Fraur Mhani Conr Conr Sruur - Rn Advan in Fraur Mhani Conr - B. H. Oh, al.(d Kora Conr Iniu, Soul, ISBN 978-89-578-8-8 Fraur mhani arly-ag onr V.T.N. Dao, P.H. Morri & P.F. Dux Shool Civil Enginring, Th

More information

ENGR 7181 LECTURE NOTES WEEK 5 Dr. Amir G. Aghdam Concordia University

ENGR 7181 LECTURE NOTES WEEK 5 Dr. Amir G. Aghdam Concordia University ENGR 78 LETURE NOTES WEEK 5 r. mir G. dam onordia Univrity ilinar Tranformation - W will now introdu anotr mtod of tranformation from -plan to t - plan and vi vra. - Ti tranformation i bad on t trapoidal

More information

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

Instructors Solution for Assignment 3 Chapter 3: Time Domain Analysis of LTIC Systems Inrucor Soluion for Aignmn Chapr : Tim Domain Anali of LTIC Sm Problm i a 8 x x wih x u,, an Zro-inpu rpon of h m: Th characriic quaion of h LTIC m i i 8, which ha roo a ± j Th zro-inpu rpon i givn b zi

More information

ENGR 7181 LECTURE NOTES WEEK 3 Dr. Amir G. Aghdam Concordia University

ENGR 7181 LECTURE NOTES WEEK 3 Dr. Amir G. Aghdam Concordia University ENGR 78 LECTURE NOTES WEEK Dr. ir G. ga Concoria Univrity DT Equivalnt Tranfr Function for SSO Syt - So far w av tui DT quivalnt tat pac ol uing tp-invariant tranforation. n t ca of SSO yt on can u t following

More information

EE750 Advanced Engineering Electromagnetics Lecture 17

EE750 Advanced Engineering Electromagnetics Lecture 17 EE75 Avan Engnrng Eltromagnt Ltur 7 D EM W onr a D ffrntal quaton of th form α α β f ut to p on Γ α α. n γ q on Γ whr Γ Γ Γ th ontour nlong th oman an n th unt outwar normal ot that th ounar onton ma a

More information

Acoustic Analysis with Consideration of Damping Effects of Air Viscosity in Sound Pathway

Acoustic Analysis with Consideration of Damping Effects of Air Viscosity in Sound Pathway Aousti Analysis with Consiration of Damping Effts of Air Visosity in Soun Pathway M. Sasajima, M. Watana, T. Yamaguhi, Y. Kurosawa, an Y. Koi Astrat Soun pathways in th nlosurs of small arphons ar vry

More information

CS553 Lecture Register Allocation I 3

CS553 Lecture Register Allocation I 3 Low-Lvl Issus Last ltur Intrproural analysis Toay Start low-lvl issus Rgistr alloation Latr Mor rgistr alloation Instrution shuling CS553 Ltur Rgistr Alloation I 2 Rgistr Alloation Prolm Assign an unoun

More information

Evans, Lipson, Wallace, Greenwood

Evans, Lipson, Wallace, Greenwood Camrig Snior Mathmatial Mthos AC/VCE Units 1& Chaptr Quaratis: Skillsht C 1 Solv ah o th ollowing or x: a (x )(x + 1) = 0 x(5x 1) = 0 x(1 x) = 0 x = 9x Solv ah o th ollowing or x: a x + x 10 = 0 x 8x +

More information

Additional Math (4047) Paper 2 (100 marks) y x. 2 d. d d

Additional Math (4047) Paper 2 (100 marks) y x. 2 d. d d Aitional Math (07) Prpar b Mr Ang, Nov 07 Fin th valu of th constant k for which is a solution of th quation k. [7] Givn that, Givn that k, Thrfor, k Topic : Papr (00 marks) Tim : hours 0 mins Nam : Aitional

More information