On interface property characterization and performance of fiber-reinforced cementitious composites

Size: px
Start display at page:

Download "On interface property characterization and performance of fiber-reinforced cementitious composites"

Transcription

1 Concrete Science an Engineering, Vol. 1, September 1999, pp On interace property characterization an perormance o iber-reinorce cementitious composites Z. Lin 1, T. Kana an V. C. Li 1 (1 Avance Civil Engineering Materials Research Laboratory, Department o Civil an Environmental Engineering, University o Michigan, Ann Arbor, MI , USA ( Kajima Technical Research Institute, Tobitakyu, Chou-shi, Tokyo 18, Japan RESEARCH PAPERS ABSTRACT It has been well recognize that iber/matrix interace properties have signiicant eects on the perormance o a iber-reinorce cement composite, incluing its racture toughness, tensile an lexural strength an uctility. Proper characterization o the interace properties via micromechanics moels can lea to eective tools or esigning high perormance an cost-eective cement-base materials. In this paper, a micromechanics moel is evelope to characterize the interace properties at single iber pullout level. In the moel, interacial racture toughness, rictional bon strength an post-eboning slip-harening coeicient are explicitly accounte or. Fiber rupture an iber strength reuction ue to incline iber pullout are also consiere. The complete composite briging stress versus crack opening curve ( B - δ relation is erive analytically. Implications o the present moel on various composite properties, incluing composite tensile strength, racture energy an complementary energy (a measure o uctility, are iscusse along with an example o PVA iber reinorce cement composites. 1. INTRODUCTION Introuction o ibers in a cement-base brittle matrix can signiicantly increase the ailure strain an racture toughness o the composite by orers o magnitue [1-4]. In orer to achieve high perormance an cost-eectiveness o such composites, quantitative materials esign tools are neee. Micromechanics moels have been proven to be very eective in high perormance an costeective cement-base composites evelopment [5-7]. A generation o super-uctile short-iber cement composites so-calle Engineere Cementitious Composites (ECC have been evelope an start to gain momentum in real-worl applications [8-1]. Presente in this paper is a generic micromechanics moel that can be use or almost all the iber types, incluing steel ibers an synthetic (polymeric ibers with both hyrophobic an hyrophilic nature. The current moel extens the work by Maalej et al. [3] to inclue chemical boning an slipepenent interacial properties. It provies not only a more realistic representation o iber-matrix interacial behaviors in many cases but also signiicantly broaens the range o viable iber types in achieving both high perormance an cost-eectiveness, especially or those ibers with strong chemical bon to cement matrix such as Poly-Vinyl-Alchohol (PVA iber, which will be investigate in etail. Fig. 1 shows the scope o the present work within the perormance riven esign approach (PDDA [11]. The remainer o this paper is organize as ollows. A theoretical single iber eboning an pullout moel is irst presente with three most important, physically meaningul parameters: (1 Chemical bon strength quantiie by interacial racture toughness, ( constant rictional bon strength or small sliing an (3 slipharening coeicient that characterizes the increasing eective rictional bon uring large sliing (pullout stage [7, 1]. Then, a composite briging stress versus crack opening ( B - δ relation is erive in close-orm base on the single iber pullout moel an ranom iber istribution assumption. As a unamental material Fig. 1 Scope o the present research within the Perormance Driven Design Approach. ISSN en cours/99 RILEM Publications S.A.R.L. 173

2 Concrete Science an Engineering, Vol. 1, September 1999 property, the B - δ relation has been well recognize as a crucial link between composite constituents an overall composite properties. Most irect an important implications o this moel incluing tensile strength an uctility are presente in the ollowing section along with an example o PVA iber-cement composite. Conclusions rawn rom the current research are mae in the en.. INTERFACE CHARACTERIZATION A commonly-use technique to investigate ibermatrix interacial behaviors is single iber pullout. Fig. shows schematically the set-up or a typical single iber pullout test. The specimens are abricate accoring to the technique or microibers escribe in [13] in orer to ensure accurate alignment. The tensile loa on the iber is measure with a miniature loa cell o 1-Newton capacity. Fig. 3 is a typical pullout curve or a PVA iber o 14.8 µm iameter. The iber embement length L e is chosen in such a way that iber rupture is avoie uring the whole process an both eboning an rictional pullout behaviors can be completely capture. As can be seen rom Fig. 3, there are three stages associate with the loa-isplacement curve: initial elastic stretching o the iber ree length (the portion not embee, ollowe by eboning stage, which is simulate in the present moel by a moe-ii tunneling crack avance with non-zero crack-tip racture toughness. The eboning stage continues until reaching the maximum loa an a istinct loa rop occurs. This loa rop is an inication o chemical boning because it woul not appear i the interace is rictionally bone only. Physically, the loa rop represents the transition rom both chemical bon an rictional bon controlle eboning stage to the pullout stage with rictional bon only. Ater complete eboning o the interace, chemical bon oes not exist but rictional bon coul eectively increase ue to abrillation o iber surace sliing against surrouning matrix. A concave upwar portion o the curve inicates this so-calle slip-harening behavior, which has been investigate elsewhere [7, 1]. A theoretical single iber eboning an pullout moel is erive, base on simple stress analysis an energy balance principle (see Appenix -A. The main assumptions mae in this moel are: (1 Fibers are o high aspect ratio (> 1 so that the en eect on the total eboning loa is negligible. This assumption is generally satisie or most available ibers, an it greatly simpliies the analysis without losing accuracy. Fig. Schematic o a single iber pullout test set-up. Fig. 3 Single iber pullout curve o a PVA iber: test result vs. moel preiction ( During the eboning stage, the slip-epenent eect is negligible since relative slippage between the iber an the matrix in the ebone portion is small. Hence, the rictional stress within the ebone zone remains at a constant τ. (3 Poisson s eect is negligible. For lexible iber/cement systems, Poisson s eect is usually iminishe ue to inevitable slight misalignment an surace roughness o the iber (no alignment at all in a ranom, short iber composite [1]. (4 Elastic stretch o the iber ater complete eboning is negligible, compare with slip magnitue. Following the etaile erivation in Appenix A, we have the theoretical pullout loa P versus pullout isplacement δ (P-δ relation: 3 3 πτe 1 η π GE P = δ, δ δ or the eboning stage, an: ( ( < < P = π τ 1 β δ δ / L δ δ, δ δ L e e (1 ( or the pullout stage, where L e is the iber embement length an δ correspons to the isplacement at which ull-eboning is complete. It is given by: 174

3 Lin, Kana, Li τ Le 1 η 8GL e 1 η δ (3 = E E where η = V E /V m E m, V an V m are the volume ractions o iber an matrix respectively. In terms o iber eboning length L, the eboning loa P can be also expresse as: 3 / P = π τ L π G E (4 At ull-eboning, L = L e ; hence, the maximum eboning loa is given by: Pa = Pb π GE 3 / (5 where P b = π τ L e is the initial rictional pullout loa. Equation (5 can be conveniently use to calibrate the chemical bon strength G an rictional bon strength τ rom the maximum eboning loa P a an the initial rictional pullout loa P b. The slip-harening coeicient β is obtaine rom ( by best-itting the rictional pullout portion o the P-δ curve. The ashe curve in Fig. 3 is the moel result with three interacial parameters: G = 6. J/m, τ = 3. MPa an β =.5. To urther veriy the moel, the maximum eboning loa P a as a unction o iber embement length L e or the same PVA iber is preicte using (5 with the same set o parameters. As can be seen rom Fig. 4, goo agreement between moel preiction an experimental result is oun. The simple P-δ relations in equations (1 to (3 are the builing blocks in constructing the composite briging law. It is also important to stuy the iber alignment eect on the pullout loa an iber in-situ strength because in an actual short iber composite ibers are ranomly oriente. The eect o iber alignment on the pullout loa, so-calle snubbing eect, was investigate in [14, 15]. An empirical relation is given by: = P P e φ φ (6 where is the snubbing coeicient. Kana an Li [16] investigate the eect o iber alignment on iber in-situ strength. Base on the experimental results or PVA ibers, a similar exponential expression was oun to characterize this eect well: = φ u φ u e (7 where u is the in-situ iber strength an is the iber strength reuction coeicient. Fig. 5 shows such a iber strength reuction eect on PVA ibers. The parameters introuce in this section, such as the slip-harening coeicient β, snubbing eect coeicient, an iber strength reuction coeicient, are empirical, curve-itting parameters. They are etermine by single Fig. 4 Maximum iber eboning stress as a unction o iber embement length or PVA (RMU iber: experiment vs. moel preiction. Fig. 5 Eect o iber inclination on in-situ iber strength o PVA (RMU ibers (ater Kana an Li [16]. iber pullout (straight or incline tests as iscusse herein an in [7, 15 an 16]. The range o β or polymeric ibers such as Polypropylene, Polyethylene an PVA ibers is between.5 to.5. Snubbing coeicient ranges rom.5 to.9. The iber strength reuction coeicient is about.3 or the PVA ibers, an it has not been teste or other iber types. 3. COMPOSITE BRIDGING STRESS- CRACK OPENING RELATION Crack briging law is a unamental material property o iber reinorce composites. It is cast in terms o the crack briging stress vs. the crack opening relation ( B -δ relation. Base on the P-δ relation o a single iber pullout, the composite B -δ curve can be obtaine by averaging over the contributions o only those iniviual ibers that 175

4 Concrete Science an Engineering, Vol. 1, September 1999 cross the matrix crack plane as the crack opens up [17]: B 4V π / L / cosφ φ ( δ = P( δ e p( φ p( z zφ π φ= z = (8 where p(φ an p(z are the probability ensity unctions o the orientation angle φ an the centroial istance o a iber rom the crack plane z, respectively. The relation between ecreasing iber iameter an increasing number o ibers has been consiere in the above equation. For a 3-D ranom istribution, p(φ = sinφ an p(z = /L. It shoul be pointe out that in the crack briging case, the crack opening δ results rom two-sie iber eboning an one-sie pullout (shorter embement sie. Thereore, P(δ in equation (8 shoul be moiie rom the one-sie eboning result (equation (1 as ollows: 3 3 P = π Eτ δ 1 η / 4 π EG / (9 By changing variable, equation (8 can be re-cast into the ollowing orm[3]: V π / 1 φ B( δ = P( δφ e ( φ φ (1 π φ=,, l sin l l = where is normalize iniviual iber embement length (= /(L /. In orer to consier the inluence o iber rupture, when carrying out the integration (1, we nee to keep iscounting the contribution o broken ibers as the crack opening δ increases. To illustrate, a typical integration omain at an arbitrary crack opening δ is plotte in Fig. 6. There are several important non-imensional variables associate with this integral omain, whose physical meanings are explaine in the ollowing: Critical iber embement length (normalize by L / or iniviual ibers: ( φ l c = L c e / (11 u where L c = is the normalize critical embement τl length o an aligne iber (φ = in absence o chemical τl bon [3] an = ( τ δ 8G /, with δ = is E ( 1 η a relative measure o chemical bon strength to the rictional bon strength. Equation (11 is obtaine rom equations (4, (6 an (7 with P set to be the critical value P = π u (φ / 4. It eines the potential iber rupture subomain: ibers that locate above this curve will eventually rupture as a matrix crack opens up. From (11, we can eine the critical iber embement length or aligne iber (φ = with both chemical bon an rictional bon: L c = Lc / (1 Fibers at all orientation angles have the potential to rupture i its hal-length L / > L c (i.e. Lc < 1. Clearly, chemical bon tens to promote iber rupture. Also, by setting φ = π/ in (11, another limiting iber length can be eine: ( π / Lr = Lce (13 I L r > L (i.e. L r >, no iber rupture will occur regarless o orientation angle. Also in Fig. 6, the variable (δ = c (φ = φ c eines the critical embement length at which eboning o the iber has complete at crack opening δ an just survive rupture. For iber length between L r an L c, some ibers will rupture, epening on their orientation angles. The limiting iber rupture angle φ b : φ b = 1 ln 1 L L c c, φb π/ (14 eines the potential iber rupture space in terms o orientation angle: only those ibers orient at an angle higher than φ b will eventually rupture. Current iber rupture angle φ c (δ: Fig. 6 Integration omain at an arbitrary crack opening δ or normal chemical bon strength case: L c > 1. φ δ δ c = 1 ln (15 ( L c δ 4 L c It eines the iber orientation angle above which ibers must have rupture at a given crack opening (the riving orce or iber rupture. This quantity is boune between φ b an π/. As crack opening δ increases, φ c starts rom π/ an approaches φ b. Note that within the potential iber rupture space, whether a speciic iber will rupture or not still epens on its embement length. 176

5 Lin, Kana, Li Equation (15 is generally not applicable to uctile metal ibers ue to the elasticity assumption in its erivation. For uctile ibers, i not too ar beyon yieling is expecte, (15 coul be a goo estimation. Base on Fig. 6, it is not iicult to evaluate the lethan sie o (1 analytically by taking only the shae omains marke as eboning an pullout ibers in the igure. Fibers that have broken or have been pulle-out (embement length less that crack opening are eliminate. Close-orm expressions or the crack briging are inclue in Appenix B, which orm the theoretical basis or iscussions in the ollowing sections. For more etaile erivation, one may reer to the previous works [3] an [16]. It is worth consiering a special case o extremely strong chemical boning ( >>1. In this case, part or all o ibers will rupture either uring the eboning process or even beore eboning (recall that there is a threshol or initiation o eboning governe by chemical bon strength (re. equation (1. Fig. 7 shows such a case. Fibers belonging to φ r < φ < π/ an < < 1 will rupture uring the eboning initiation process, where φ r is given by: Fig. 7 Integration omain at an arbitrary crack opening δ or extreme chemical bon strength case: L c < 1, L r <. φr = 1 ln (16 L c The maximum briging stress is attaine at the initial stage (δ =, an it can be estimate by summing the contributions rom eboning an non-eboning (rupture beore eboning ibers: V φr 1 φ B = e sinφl φ (17 V π / 1 φ φ φ φ ue sin l r where = an = 1 / Vτ ( L /. It results in: V Fig. 8 Pre-ully-eboning briging stress - crack opening curves or a system with low tenency o iber rupture. ( r π φ V / r u e e sinφr cosφ B = 4 φr e ( sinφr cosφr 4 (18 Figs. 8 to 1 show graphically the non-imensional crack briging stress - crack opening curves erive analytically. Figs. 8 an 9 are the pre-ully-eboning curves or iber/matrix systems with low tenency o iber rupture (L c =. an high tenency o iber rupture (L c =.5. In both Figs. 8 an 9, the parameter / L c (rom. to.9 inicates the relative chemical bon strength, Fig. 9 Pre-ully-eboning briging stress - crack opening curves or a system with high tenency o iber rupture. given other iber an interace characteristics such as iber length, iameter, strength an rictional bon strength. It can be seen rom these two igures, chemical 177

6 Concrete Science an Engineering, Vol. 1, September 1999 tenency is low, chemical bon increases the maximum briging stress much more ramatically than the case with high iber rupture tenency (L c =.5. Also, at the en o eboning stage, because o less severe iber rupture, the briging stress remains higher at a larger corresponing crack opening in L c =. case. This also contributes to larger post-eboning racture energy - area uner the escening part o the B -δ curve (see Figs. 1 an 11. In the case o slip-harening (Fig. 1, both briging stress an racture energy can be urther increase, epening on the slip-harening coeicient. Fig. 1 Post-ully-eboning briging stress - crack opening curves or a system with low tenency o iber rupture (continuation o the curves in Fig MODEL IMPLICATIONS ON MATERIALS DESIGN FOR DESIRED PROPERTIES With the analytical tools evelope, it is convenient to preict the composite overall properties rom its constituents - iber, matrix an interace. In this section, some irect implications rom the current micromechanics moel are presente, incluing tensile strength, racture energy an complementary energy (a measure o material uctility. Tensile Strength, Fracture Energy, Ductility Fig. 11 Post-ully-eboning briging stress - crack opening curves or a system with high tenency o iber rupture (continuation o the curves in Fig. 9. The tensile strength preicte by the moel is essentially the maximum briging stress or a iber/cement composite which shows multiple cracking (re. [18]. Fig. 13 shows the epenency o tensile strength on iber length in a nonimensional orm. Given iber volume raction V, iameter, iber strength u an rictional bon strength τ, the composite tensile strength cu monotonically increases with iber length in zero chemical boning case (the lowest curve in Fig. 13. This is consistent with the result erive by Fig. 1 Eect o slip-harening interace on the post-ullyeboning B -δ relation. bon tens to increase initial briging stress an reuce the complementary area to the let o the B -δ curve up to the maximum briging stress point (this area is reerre as the complementary energy [4]. When iber rupture Fig. 13 Normalize composite tensile strength as a unction o normalize iber length an relative chemical bon strength. 178

7 Lin, Kana, Li Fig. 14 Normalize composite racture energy as a unction o normalize iber length an relative chemical bon strength. Maalej et al. [3]. As chemical bon increases, tensile strength or short iber length can be ramatically increase. The plateaus at small iber lengths inicate that composite tensile strength is controlle by the initial eboning stress which is in turn governe by chemical bon. As iber becomes longer, subsequent eboning ater initial eboning results in higher briging stress an composite tensile strength. From this igure, it is clear that when iber is too long or chemical bon is too strong, composite tensile strength becomes almost inepenent o iber length. Physically, those two limiting cases lea to severe iber rupture an composite tensile strength shoul be mainly controlle by iber volume raction an iber strength. Shown in Fig. 14 is composite racture energy as a unction o iber length an the relative chemical bon strength. This racture energy is the energy absorbe by briging ibers in the crack plane an it is evaluate by integrating the area uner the post-ully-eboning portion o a briging stress - crack opening curve (Figs The peak shown in each curve in Fig. 14 inicates that there exists an optimal iber length to maximize the racture energy o a iber composite. It is also clear that chemical bon tens to reuce racture energy because it promotes more iber rupture. The concept o multiple cracking in iber reinorce brittle matrix composites has been extensively stuie. Marshall an Cox [18] propose a simple means o etermining the conition or steay state cracking, require or multiple cracking. Base on a J-integral analysis, the conition can be written in terms o the complementary energy G c an the crack tip toughness J : Gc J (19 where G c is eine in terms o iber briging property via the B -δ curve: G δ c cu cu B = cu δ δ δ ( an cu an δ cu are the maximum briging stress an the corresponing crack opening. Graphically, G c is simply the complementary area to the let o the B -δ curve up to peak stress. The above iscussion suggests that high uctility o a short iber composite requires large complementary energy. Equation (19 becomes a very important esign criterion or uctile iber composites. Using the generic micromechanics moel evelope in this research, one can easily compute this quantity with ierent combinations o iber, matrix an interace properties an choose the one that satisies conition (19 to guarantee the uctility o such engineere iber composites. Shown in Fig. 15 are the normalize complementary energy as a unction o iber length an relative chemical bon strength. Fig. 15 Normalize composite complementary energy as a unction o normalize iber length an relative chemical bon strength. An Example - PVA iber/cement composites Recently, Polyvinyl Alcohol (PVA iber has been recognize or its high potential as reinorcement in high perormance cementitious composites ue to its high strength. This type o iber can also orm strong chemical bon to the cement matrix. In orer to maximize the potential o PVA ibers as viable reinorcement or cement base composites which have both high strength an high uctility, analytical material esign tools are certainly neee. To emonstrate the useulness o the present micromechanics moel, two types o cement composites mae rom two types o PVA ibers (RMU an RK with 14 an 4 µm iameter respectively are investigate. Table 1 lists the iber an interace parameters or both PVA ibers. The interacial properties were erive rom the simple iber pullout tests iscusse previously. PVA-RMU composite contains 1.5% volume raction o RMU ibers while PVA-RK composite has % RK ibers. 179

8 Concrete Science an Engineering, Vol. 1, September 1999 Table 1 Fiber an interacial parameters or PVA-RMU an RK ibers Fiber E L τ G u, Type (GPa (mm (µm (MPa (Nm (MPa RMU ,.3 RK ,.3 Table Matrix mix proportion Cement San Water Super Viscous Plasticizer Agent The matrix mix proportion or the two composites are liste in Table. They were teste ater seven ays water curing an one ay air rying. Stanar og-bone shape specimens were use. Fig. 16 shows the uniaxial tension test results or both composites. As can be seen rom this igure, PVA-RMU has higher tensile strength but much lower uctility than the PVA-RK composite. This is because RMU iber is much thinner in iameter. Although it is a higher iber strength than the RK iber, uner similar interacial conitions, RMU ibers suere more severe iber rupture, which leas to lower uctility. The contrast in uctility between these two composites can be explaine more clearly using the complementary concept later on. Fig. 17 shows the comparison between moel preictions an tensile test results in terms o tensile strength. The moel preicts higher tensile strength or PVA-RMU than or PVA-RK, consistent with experimental observations. Note however, that the PVA-RMU oes not satisy the multiple cracking conition (equation (19 (see iscussion to ollow an the measure tensile strength cannot be expecte to be the same as the calculate maximum briging stress. Also, the current tensile strengths or both composites are very close to their limit values an reuction in iber length will not cause much loss in tensile strength. From the racture energy preictions (Fig. 18, the current iber lengths (1 mm or RK, 6 mm or RMU are well over the optimal values or maximum racture energy. The ierence in tensile strength as L approaches zero is ue to the iber strength ierence. Fig. 19 illustrates the complementary energy as a unction o iber length or both composites. Uner current composite esign, PVA-RK has a G c value o 3.4 J/m while PVA-RMU has only.8 J/m. Compare to the racture toughness o cement paste, which is about J/m, the PVA-RK composite satisies the multiple cracking conition (19 an this is the reason behin the contrast in tensile uctility o the two composites. Fig. 16 Uniaxial tensile test results or PVA-RMU an PVA-RK composites. Fig. 17 Composite tensile strength as a unction o iber length or both PVA-RMU an RK composites: experiments vs moel preictions. Fig. 18 Preicte composite racture energy as a unction o iber length or both PVA-RMU an RK composites. Fig. shows the eect o chemical bon strength G on both racture energy an tensile strength o the PVA- RMU composite. Clearly, chemical bon ecreases the racture energy ramatically while increases the tensile 18

9 Lin, Kana, Li Fig. 19 Preicte composite complementary energy as a unction o iber length or both PVA-RMU an RK composites. (3 Via parametric stuy presente in this paper, it is oun that in iber rupture cases, interacial racture toughness (chemical bon ten to increase composite tensile strength but ecreases composite racture energy an complementary energy as well as uctility. Single iber pullout problem is a complicate soli mechanics problem. The erivation inclue in the present stuy is only an approximate solution to this complicate problem, base on simpliie orce an energy balance. The attempt has been mae to capture the physical essence o the interacial parameters in a simple, easy-to-use ormulation without getting into complex mathematics. In the current micromechanics moel, the iber ispersion eect an the statistical istribution o iber strength are not inclue. The comparison between the two iber/cement composites PVA-RMU an PVA-RK is meant to emonstrate the useulness o the moel in iber section or esire properties (say, uctility rather than just a comparision, which has been the central theme o this work. ACKNOWLEDGMENTS Support o this work has been provie by the National Science Founation (CMS-9616 to the University o Michigan. Fig. Preicte composite racture energy an tensile strength as a unction o chemical bon strength or the PVA-RMU composite. strength to the limit set by iber strength an volume raction. The kink in the tensile strength curve is ue to the transition rom normal chemical boning to extreme chemical boning case as iscusse previously (Re. equation (18. It is important to tailor the interacial properties incluing chemical bon to achieve both high strength an high uctility. 4. CONCLUSIONS AND DISCUSSIONS (1 Several important interacial behaviors characterize: chemical bon, rictional bon, an slip-harening. These interacial parameters can be etermine rom a single iber pullout test. ( An analytical briging stress-crack opening relation is erive, which can be use as a materials esign tool or esire properties an cost-eectiveness via intelligent iber selection (type, strength, length, content,, matrix moiication an iber/matrix interace tailoring. REFERENCES [1]Krenchel, H. an Hansen, S., New recipes an new prouction techniques or high perormance FRC-materials, in High Perormance Fiber Reinorce Cement Composites, Reinhart H. W., Naaman, A. E., Es., E& FN Spon, Lonon, 199, pp [] Naaman, A. E., SIFCON: Tailore properties or structural perormance, in High Perormance Fiber Reinorce Cement Composites, Reinhart H. W, Naaman, A. E., Es., E& FN Spon, Lonon, 199; pp [3] Maalej M., Li, V. C. an Hashia, T., Design an Structural Applications o Stress-Crack With Relations in Fiber Reinorce Concrete, J. o Engineering Mechanics 8 ( [4] Li, V. C., Kana, T. an Lin, Z., The Inluence o Fiber/Matrix Interace Properties on Complementary Energy an Composite Damage Tolerance, in Key Engineering Materials: Proc. 3r Con. on Frac. & Strength o Solis, Hong Kong, 1997, pp [5] Li, V. C. an Leung, C. K. Y., Steay state an multiple cracking o short ranom iber composites, ASCE J. Engng Mech. 118 ( [6] Li, V. C. an Wu, H. C., Conitions or pseuo strain-harening in Fiber Reinorce Brittle Matrix Composites, J. Appl. Mech. Rev. 45 (8 ( [7] Lin, Z. an Li, V. C., Crack briging in iber reinorce cementitious composites with slip-harening interaces, J. 181

10 Concrete Science an Engineering, Vol. 1, September 1999 Mech. Phys. Solis 45 (5 ( [8] Li, V. C. an Kana, T., Engineere cementitious composites or structural applications, ASCE J. Mater. Civil Eng. 1 ( ( [1] Kabele, P., Li, V. C., Horii, H., Kana, T. an Takeuchi, S., Use o BMC or Ductile Structural Members, in Proc. o 5th Int. Sym. on Brittle Matrix Composites, Warsaw, Polan, 1997; pp [1] Kana, T., Watanabe, S. an Li, V. C., Application o Pseuo Strain Harening Cementitious Composites to Shear Resistant Structural Elements, in Proc. Fracture Mechanics o Concrete Structures-3, AEDIFICATIO Publishers, Freiburg, Germany, Oct., 1998, pp [11] Li, V. C., From micromechanics to structural engineering - the esign o cementitious composites or civil engineering applications, JSCE J. Struct. Mech. Earthquake Eng. 1 ( ( [1] Wang, Y., Li, V. C. an Backer, S., Moeling o iber pull-out rom a cement matrix, Int l. J. Cement Comp. Lightweight Concrete 1 (3 ( [13] Katz, A., Li, V. C., A Special technique or etermining the bon strength o carbon ibers in cement matrix by pullout test, J. Materials Science Letters 15 ( [14] Morton, J. an Groves, G. W., The eect o metal wires on the racture o a brittle matrix composite, J. Mater. Sci. 11 ( [15] Li, V. C., Post-crack scaling relations or iber reinorce cementitious composites, ASCE J. Mater. Civil Engng. 4 (1 ( [16] Kana, T. an Li, V. C., Interace property an apparent strength o a high strength hyrophilic iber in cement matrix, ASCE J. Mater. Civil Engng. 1 (1 ( [17] Li, V. C., Wang, Y. an Backer, S., A micromechanical moel o tension-sotening an briging toughening o short ranom iber reinorce brittle matrix composites, J. Mech. Phys. Solis 39 (5 ( [18] Marshall, D. B. an Cox, B. N., A J-integral metho or calculating steay-state matrix cracking stress in composites, Mech. Mater. 7 ( [19] Gao, Y. C., Mai, Y.-W. an Cotterell, B., Fracture o iberreinorce materials, J. Appl. Math. Phys. 39 ( Appenix A APPENDIX Derivation o P-δ relation with combine chemical bon an rictional bon Consier a iber with length L e embee in a matrix (Fig. A1. A portion o the iber/matrix interace is ebone. The ebone zone length is L, within which a constant rictional stress τ is assume to exist. Chemical bon strength or the bone region is G. The task here is to erive the relation between iber pullout loa P or stress an iber pullout isplacement δ. From equilibrium requirements, the one-imensional Fig. A1. stress istribution in the ebone region is given by: z ( z = L (A-1 m z 1 z/l m (A- where an m are the normel stresses in the iber an the matrix respectively at z =, an they are given by: 4τL = (A-3 m (A-4 where V an V m are the volume ractions o iber an matrix respectively. Deine the relative isplacement between iber an matrix in the ebone region as: = m z u z u z Then, = ( 4τLV = V m = = (A-5 z u z um z z m z (A-6 z z z E Em Substituting (A-1 through (A-4 an applying the bounary conition (z = =, we have: ( z z Lz z = 4τ τ 1 η 1 η E (A-7 E E where η= VE. The iber pullout istance is then VmEm given by δ τ η = ( L = L L 1 (A-8 E E Now, we nee to in = (L, G, τ,. An energy base eboning criterion is employe in establishing such a relation [19]. For any ininitesimal avance o the ebone zone, A, the energy conservation requires 18

11 Lin, Kana, Li Pu = U W GA (A-9 It ollows that where u is the axial isplacement at the pulle iber en (work conjugate o applie loa P = π τ η η / 4, U is the = L ( GE ( change in strain energy in the system, W is the change in (A-16 energy issipate by riction at the interace, an G A is or: the energy consume in the avance o ebone zone. On 3 the other han, i the ebone zone is consiere as a π GE ( 1 η P= π/ 4= πτl( 1 η stress bounary o the system, it can be shown that [19]: (A-17 U = 1 In a single iber pullout specimen, η can be neglecte ( Pu W (A-1 without losing accuracy. Thereore, equation (5 in the From (A-9 an (A-1, we have: main text results rom (A-17. Also, when L = L e, maximum eboning loa is reache. The corresponing is- GA = 1 ( Pu W (A-11 placement is given by (A-16 an (A-8: where u an W can be obtaine rom: τl e( 1 η 8GL e 1 η L ( z u (A-1 E z V = δ = ( (A-18 ( E L e L E E c Combining (A-17 an (A-8, we have: an 3 3 πτ L E W = π ( 1 η π GE ( z z (A-13 P = δ (A-19 In equation (A-1, I η =, (A-19 reuces to equation ( in the main Ec = VE VmEm (A-14 text. By using relevant equations alreay erive, (A-11 becomes: Ater eboning, the interace is controlle by rictional bon only. By taking S = δ-δ in equation (1 in the τ η τ η η main text an rom orce balance, it is reaily shown that 8 L( 1 16 L GE( = (A-15 or iber pullout stage, P= π τ 1 β δ δ / L δ δ (A- ( ( e Appenix B Analytical crack briging relation For L r < L < L c B δ δ φ π g 41 ( ( 1 u = G φ ( c, 41 ( ( 1 Lc φc δ δ A(, L c A φ φ δ φ π c, A ( 1 ( c, u = 4 ( φ, ( L ( φ, δ G b 1 c A b Lc A( φb, δ Lr A φ ( b, ( 1 c δ δ δ G( b, ( 1 L c B( b, a, L c φ δ φ φ δ B( φb, φa, δ B φ ( b, φa, L ( 1 cδ r 1 183

12 Concrete Science an Engineering, Vol. 1, September 1999 B where For L > L c δ δ φ π g 41 ( ( 1 u = G φ ( c, 41 ( ( 1 LcA φc δ δ (, L ca φ φ δ φ π ( c, 1 A c, u = u φ = 4 π L G L δ G π c, c, G π, ( 1 c ( = L r δ δu φ δ L cg( a, L φ c δ G( φa, G δ φa, ( 1 c Lr L c 1 βl Vτ L /, c = = τ L = 8G / τδ, δ, δ = δ 1, δ / L / E 1 η L u c =, = Lc/ ( L / = Lc Lc, L Lr / L / L ce τ L r φ φ φ u c b a = Lc e ( 1 ( φ δ φ L c = 1 δ ln ( c Lc c 4L, = / 1 c 1 1 = ln ( L c L c = 1 ln L L = c c = π / = = For 3-D istribution: [ ] = φ G φ, 1 e sinφ cosφ 4 g G π /, = = [ ] φ π A φ, 1 / e cosφ sinφ e 4 φ e sin φ cos φ B( φ, φ, = φ1 e sinφ cosφ ( 1 1 For -D istribution: = φ G φ, 4 e cosφ sinφ π 1 g G π /, = = [ ] [ ] π φ A φ, 4 / e e cosφ sinφ π( 1 φ e ( cosφ sinφ B( φ, φ, = 4 1 φ π( 1 1 e ( cosφ1 sinφ 1 184

VUMAT for Fabric Reinforced Composites

VUMAT for Fabric Reinforced Composites VUMAT or Fabric Reinorce Composites. Introuction This ocument escribes a constitutive mo or abric reinorce composites that was introuce in Abaqus/Exicit 6.8. The mo has been imemente as a built-in VUMAT

More information

Debond crack growth in fatigue along fiber in UD composite with broken fibers

Debond crack growth in fatigue along fiber in UD composite with broken fibers Debon crack growth in atigue along iber in UD composite with broken ibers Johannes Eitzenberger an Janis arna Dept o Applie Physics an Mechanical Engineering Lulea University o Technology, Sween SE 97

More information

3-D FEM Modeling of fiber/matrix interface debonding in UD composites including surface effects

3-D FEM Modeling of fiber/matrix interface debonding in UD composites including surface effects IOP Conference Series: Materials Science an Engineering 3-D FEM Moeling of fiber/matrix interface eboning in UD composites incluing surface effects To cite this article: A Pupurs an J Varna 2012 IOP Conf.

More information

Influence of Interfacial Delamination on Channel Cracking of Elastic Thin Films. Haixia Mei, Yaoyu Pang, Rui Huang

Influence of Interfacial Delamination on Channel Cracking of Elastic Thin Films. Haixia Mei, Yaoyu Pang, Rui Huang Inluence o Interacial Delamination on Channel Cracking o Elastic Thin Films Haixia Mei, Yaoyu Pang, Rui Huang Department o Aerospace Engineering an Engineering Mechanics, University o Texas, Austin, TX

More information

Thermal conductivity of graded composites: Numerical simulations and an effective medium approximation

Thermal conductivity of graded composites: Numerical simulations and an effective medium approximation JOURNAL OF MATERIALS SCIENCE 34 (999)5497 5503 Thermal conuctivity of grae composites: Numerical simulations an an effective meium approximation P. M. HUI Department of Physics, The Chinese University

More information

Damage in fiber bundle models

Damage in fiber bundle models Eur Phys J B 17, 269 279 (2) THE EUROPEAN PHYSICAL JOURNAL B c EDP Sciences Società Italiana i isica Springer-Verlag 2 Damage in iber bunle moels Kun 1,2,a, S Zapperi 3,4, an HJ Herrmann 1,3 1 ICA 1, University

More information

Research on Thin Film Thickness Uniformity for Deposition of Rectangular Planar Sputtering Target

Research on Thin Film Thickness Uniformity for Deposition of Rectangular Planar Sputtering Target Available online at www.scienceirect.com Physics Proceia 3 (01 ) 903 913 18 th International Vacuum Congress (IVC-18) Research on Thin Film Thickness Uniormity or Deposition o Rectangular Planar puttering

More information

Modeling of vegetated rivers for inbank and overbank ows

Modeling of vegetated rivers for inbank and overbank ows Loughborough University Institutional Repository Moeling o vegetate rivers or inbank an overbank ows This item was submitte to Loughborough University's Institutional Repository by the/an author. Citation:

More information

1. Filling an initially porous tube under a constant head imposed at x =0

1. Filling an initially porous tube under a constant head imposed at x =0 Notes on Moving Bounary problems, Voller U o M, volle00@umn.eu. Filling an initially porous tube uner a constant hea impose at x =0 Governing equation is base on calculating the water volume lux by the

More information

Theoretical Investigation Heat Transfer Mechanisms in Nanofluids and the Effects of Clustering on Thermal Conductivity

Theoretical Investigation Heat Transfer Mechanisms in Nanofluids and the Effects of Clustering on Thermal Conductivity International Journal o Bioscience, Biochemistry an Bioinormatics, Vol., No., March 01 Theoretical Investigation Heat Transer Mechanisms in Nanoluis an the Eects o Clustering on Thermal Conuctivity Mohamma

More information

A DYNAMIC MODEL OF POLYELECTROLYTE GELS

A DYNAMIC MODEL OF POLYELECTROLYTE GELS A DYNAMIC MODEL OF POLYELECTROLYTE GELS M.CARME CALDERER, HAORAN CHEN, CATHERINE MICEK, AND Y. MORI Abstract. We erive a moel o the couple mechanical an electrochemical eects o polyelectrolyte gels. We

More information

Analytic Scaling Formulas for Crossed Laser Acceleration in Vacuum

Analytic Scaling Formulas for Crossed Laser Acceleration in Vacuum October 6, 4 ARDB Note Analytic Scaling Formulas for Crosse Laser Acceleration in Vacuum Robert J. Noble Stanfor Linear Accelerator Center, Stanfor University 575 San Hill Roa, Menlo Park, California 945

More information

inflow outflow Part I. Regular tasks for MAE598/494 Task 1

inflow outflow Part I. Regular tasks for MAE598/494 Task 1 MAE 494/598, Fall 2016 Project #1 (Regular tasks = 20 points) Har copy of report is ue at the start of class on the ue ate. The rules on collaboration will be release separately. Please always follow the

More information

Modelling dowel action of discrete reinforcing bars in cracked concrete structures

Modelling dowel action of discrete reinforcing bars in cracked concrete structures Title Moelling owel action of iscrete reinforcing bars in cracke concrete structures Author(s) Kwan, AKH; Ng, PL; Lam, JYK Citation The 2n International Symposium on Computational Mechanics an the 12th

More information

MULTISCALE FRICTION MODELING FOR SHEET METAL FORMING

MULTISCALE FRICTION MODELING FOR SHEET METAL FORMING MULTISCALE FRICTION MODELING FOR SHEET METAL FORMING Authors J. HOL 1, M.V. CID ALFARO 2, M.B. DE ROOIJ 3 AND T. MEINDERS 4 1 Materials innovation institute (M2i) 2 Corus Research Centre 3 University of

More information

THREE-DIMENSIONAL THERMO-PORO-MECHANICAL MODELING OF RESERVOIR STIMULATION AND INDUCED MICROSEISMICITY IN GEOTHERMAL RESERVOIR

THREE-DIMENSIONAL THERMO-PORO-MECHANICAL MODELING OF RESERVOIR STIMULATION AND INDUCED MICROSEISMICITY IN GEOTHERMAL RESERVOIR PROCEEDINGS, Thirty-Sixth Workshop on Geothermal Reservoir Engineering Stanor University, Stanor, Caliornia, January 31 - February 3, 2011 SGP-TR-191 THREE-DIMENSIONAL THERMO-PORO-MECHANICAL MODELING OF

More information

EVALUATION OF LIQUEFACTION RESISTANCE AND LIQUEFACTION INDUCED SETTLEMENT FOR RECLAIMED SOIL

EVALUATION OF LIQUEFACTION RESISTANCE AND LIQUEFACTION INDUCED SETTLEMENT FOR RECLAIMED SOIL 386 EVALUATION OF LIQUEFACTION RESISTANCE AND LIQUEFACTION INDUCED SETTLEMENT FOR RECLAIMED SOIL Lien-Kwei CHIEN 1, Yan-Nam OH 2 An Chih-Hsin CHANG 3 SUMMARY In this stuy, the fille material in Yun-Lin

More information

Optimal design of reinforced concrete T-sections in bending

Optimal design of reinforced concrete T-sections in bending Engineering Structures 5 (003) 951 964 www.elsevier.com/locate/engstruct Optimal esign o reinorce concrete T-sections in ening C.C. Ferreira a,, M.H.F.M. Barros a, A.F.M. Barros a Department o Civil Engineering,

More information

Strength Analysis of CFRP Composite Material Considering Multiple Fracture Modes

Strength Analysis of CFRP Composite Material Considering Multiple Fracture Modes 5--XXXX Strength Analysis of CFRP Composite Material Consiering Multiple Fracture Moes Author, co-author (Do NOT enter this information. It will be pulle from participant tab in MyTechZone) Affiliation

More information

Least-Squares Regression on Sparse Spaces

Least-Squares Regression on Sparse Spaces Least-Squares Regression on Sparse Spaces Yuri Grinberg, Mahi Milani Far, Joelle Pineau School of Computer Science McGill University Montreal, Canaa {ygrinb,mmilan1,jpineau}@cs.mcgill.ca 1 Introuction

More information

'HVLJQ &RQVLGHUDWLRQ LQ 0DWHULDO 6HOHFWLRQ 'HVLJQ 6HQVLWLYLW\,1752'8&7,21

'HVLJQ &RQVLGHUDWLRQ LQ 0DWHULDO 6HOHFWLRQ 'HVLJQ 6HQVLWLYLW\,1752'8&7,21 Large amping in a structural material may be either esirable or unesirable, epening on the engineering application at han. For example, amping is a esirable property to the esigner concerne with limiting

More information

The Phenomenon of Anomalous Rock Embrittlement

The Phenomenon of Anomalous Rock Embrittlement https://papers.acg.uwa.eu.au/p/574_29_tarasov/ B.G. Tarasov A.V. Dyskin School of Civil an Resource Engineering The University of Western Australia The paper analyses a phenomenon of rock behaviour - the

More information

EXPERIMENTAL AND NUMERICAL STUDY ON FREE SPANNING PIPELINES SUBJECTED TO EARTHQUAKES

EXPERIMENTAL AND NUMERICAL STUDY ON FREE SPANNING PIPELINES SUBJECTED TO EARTHQUAKES The th Worl Conerence on Earthquake Enineerin October -7, 8, Beijin, China EXPERIMENTAL AND NUMERICAL STUDY ON FREE SPANNING PIPELINES SUBJECTED TO EARTHQUAKES X. Li, Q. Jin, R.B. Don 3 an J. Zhou Associate

More information

Optimum design of tuned mass damper systems for seismic structures

Optimum design of tuned mass damper systems for seismic structures Earthquake Resistant Engineering Structures VII 175 Optimum esign of tune mass amper systems for seismic structures I. Abulsalam, M. Al-Janabi & M. G. Al-Taweel Department of Civil Engineering, Faculty

More information

A note on the Mooney-Rivlin material model

A note on the Mooney-Rivlin material model A note on the Mooney-Rivlin material moel I-Shih Liu Instituto e Matemática Universiae Feeral o Rio e Janeiro 2945-97, Rio e Janeiro, Brasil Abstract In finite elasticity, the Mooney-Rivlin material moel

More information

Qubit channels that achieve capacity with two states

Qubit channels that achieve capacity with two states Qubit channels that achieve capacity with two states Dominic W. Berry Department of Physics, The University of Queenslan, Brisbane, Queenslan 4072, Australia Receive 22 December 2004; publishe 22 March

More information

International Conference on Advances in Energy, Environment and Chemical Engineering (AEECE-2015)

International Conference on Advances in Energy, Environment and Chemical Engineering (AEECE-2015) International Conference on Avances in Energy, Environment an Chemical Engineering (AEECE-2015) Stuy on Damage Characteristic of Unergroun Cavern Blasting Excavation Base on Dynamic Damage Constitutive

More information

OF CHS. associated. indicate. the need. Rio de Janeiro, Brazil. a) Footbridge Rio. d) Maria Lenk. CHS K joints

OF CHS. associated. indicate. the need. Rio de Janeiro, Brazil. a) Footbridge Rio. d) Maria Lenk. CHS K joints EUROSTEEL 2, August 3 September 2, 2, Buapest, Hungary A NUMERICAL EVALUATION OF CHS T JOINTS UNDER AXIAL LOADS Raphael S. a Silva a, Luciano R. O. e Lima b, Pero C. G. a S. Vellasco b, José G. S. a Silva

More information

Crack onset assessment near the sharp material inclusion tip by means of modified maximum tangential stress criterion

Crack onset assessment near the sharp material inclusion tip by means of modified maximum tangential stress criterion Focuse on Mechanical Fatigue of Metals Crack onset assessment near the sharp material inclusion tip by means of moifie maximum tangential stress criterion Onřej Krepl, Jan Klusák CEITEC IPM, Institute

More information

Discrete Time Linear Quadratic Optimal Control for an Electrical Servo Drive System

Discrete Time Linear Quadratic Optimal Control for an Electrical Servo Drive System Discrete ime Linear Quaratic Optimal Control or an Electrical Servo Drive System CORNELIU BOAN, FLORIN OSAFI an ALEXANDRU ONEA Dept. o Automatic Control an Inustrial Inormatics Gh. Asachi echnical University

More information

Bond strength model for interfaces between nearsurface mounted (NSM) CFRP strips and concrete

Bond strength model for interfaces between nearsurface mounted (NSM) CFRP strips and concrete University o Wollongong Research Online Faculty o Engineering and Inormation Sciences - Papers: Part A Faculty o Engineering and Inormation Sciences 2014 Bond strength model or interaces between nearsurace

More information

Application of the homotopy perturbation method to a magneto-elastico-viscous fluid along a semi-infinite plate

Application of the homotopy perturbation method to a magneto-elastico-viscous fluid along a semi-infinite plate Freun Publishing House Lt., International Journal of Nonlinear Sciences & Numerical Simulation, (9), -, 9 Application of the homotopy perturbation metho to a magneto-elastico-viscous flui along a semi-infinite

More information

ensembles When working with density operators, we can use this connection to define a generalized Bloch vector: v x Tr x, v y Tr y

ensembles When working with density operators, we can use this connection to define a generalized Bloch vector: v x Tr x, v y Tr y Ph195a lecture notes, 1/3/01 Density operators for spin- 1 ensembles So far in our iscussion of spin- 1 systems, we have restricte our attention to the case of pure states an Hamiltonian evolution. Toay

More information

6. Friction and viscosity in gasses

6. Friction and viscosity in gasses IR2 6. Friction an viscosity in gasses 6.1 Introuction Similar to fluis, also for laminar flowing gases Newtons s friction law hols true (see experiment IR1). Using Newton s law the viscosity of air uner

More information

05 The Continuum Limit and the Wave Equation

05 The Continuum Limit and the Wave Equation Utah State University DigitalCommons@USU Founations of Wave Phenomena Physics, Department of 1-1-2004 05 The Continuum Limit an the Wave Equation Charles G. Torre Department of Physics, Utah State University,

More information

MICROMECHANICAL FAILURE ANALYSIS OF UNIDIRECTIONAL FIBER-REINFORCED COMPOSITES UNDER IN-PLANE AND TRANSVERSE SHEAR

MICROMECHANICAL FAILURE ANALYSIS OF UNIDIRECTIONAL FIBER-REINFORCED COMPOSITES UNDER IN-PLANE AND TRANSVERSE SHEAR THE 19 TH INTERNATIONAL CONFERENCE ON COMPOSITE MATERIALS MICROMECHANICAL FAILURE ANALYSIS OF UNIDIRECTIONAL FIBER-REINFORCED COMPOSITES UNDER IN-PLANE AND TRANSVERSE SHEAR Lei Yang*, Ying Yan, Zhiguo

More information

Study on aero-acoustic structural interactions in fan-ducted system

Study on aero-acoustic structural interactions in fan-ducted system Stuy on aero-acoustic structural interactions in fan-ucte system Yan-kei CHIANG 1 ; Yat-sze CHOY ; Li CHENG 3 ; Shiu-keung TANG 4 1,, 3 Department of Mechanical Engineering, The Hong Kong Polytechnic University,

More information

Lectures - Week 10 Introduction to Ordinary Differential Equations (ODES) First Order Linear ODEs

Lectures - Week 10 Introduction to Ordinary Differential Equations (ODES) First Order Linear ODEs Lectures - Week 10 Introuction to Orinary Differential Equations (ODES) First Orer Linear ODEs When stuying ODEs we are consiering functions of one inepenent variable, e.g., f(x), where x is the inepenent

More information

Harmonic Modelling of Thyristor Bridges using a Simplified Time Domain Method

Harmonic Modelling of Thyristor Bridges using a Simplified Time Domain Method 1 Harmonic Moelling of Thyristor Briges using a Simplifie Time Domain Metho P. W. Lehn, Senior Member IEEE, an G. Ebner Abstract The paper presents time omain methos for harmonic analysis of a 6-pulse

More information

Model for predicting the UHPFRC tensile hardening response

Model for predicting the UHPFRC tensile hardening response UHPC-2008: The Second International Symposium on Ultra High Performance Concrete, March 05-07, 2008, Kassel, Germany. John Wuest Dr., Civil Engineer Ecole Polytechnique Fédérale de Lausanne Lausanne, Switzerland

More information

3.2 Shot peening - modeling 3 PROCEEDINGS

3.2 Shot peening - modeling 3 PROCEEDINGS 3.2 Shot peening - moeling 3 PROCEEDINGS Computer assiste coverage simulation François-Xavier Abaie a, b a FROHN, Germany, fx.abaie@frohn.com. b PEENING ACCESSORIES, Switzerlan, info@peening.ch Keywors:

More information

Ductility and Failure Modes of Single Reinforced Concrete Columns. Hiromichi Yoshikawa 1 and Toshiaki Miyagi 2

Ductility and Failure Modes of Single Reinforced Concrete Columns. Hiromichi Yoshikawa 1 and Toshiaki Miyagi 2 Ductility an Failure Moes of Single Reinforce Concrete Columns Hiromichi Yoshikawa 1 an Toshiaki Miyagi 2 Key Wors: seismic capacity esign, reinforce concrete column, failure moes, eformational uctility,

More information

Longitudinal Waves in a Rotating Solid Cylinder Immersed in an Inviscid Fluid Using Chebyshev Polynomial

Longitudinal Waves in a Rotating Solid Cylinder Immersed in an Inviscid Fluid Using Chebyshev Polynomial Longituinal Waves in a Rotating Soli Cyliner Immerse in an Invisci Flui Using Chebyshev Polynomial R. Selvamani *1, P.Ponnusamy 1 Department o Mathematics, Karunya University, Coimbatore, TamilNau, Inia

More information

Estimation of hardness by nanoindentation of rough surfaces

Estimation of hardness by nanoindentation of rough surfaces Journal of MATERIALS RESEARCH Welcome Comments Help Estimation of harness by nanoinentation of rough surfaces M. S. Bobji an S. K. Biswas Department of Mechanical Engineering, Inian Institute of Science,

More information

Approaches for Predicting Collection Efficiency of Fibrous Filters

Approaches for Predicting Collection Efficiency of Fibrous Filters Volume 5, Issue, Summer006 Approaches for Preicting Collection Efficiency of Fibrous Filters Q. Wang, B. Maze, H. Vahei Tafreshi, an B. Poureyhimi Nonwovens Cooperative esearch Center, North Carolina State

More information

FATIGUE DURABILITY OF CONCRETE EXTERNALLY STRENGTHENED WITH FRP SHEETS

FATIGUE DURABILITY OF CONCRETE EXTERNALLY STRENGTHENED WITH FRP SHEETS FATIGUE DURABILITY OF CONCRETE EXTERNALLY STRENGTHENED WITH FRP SHEETS H. Diab () and Zhishen Wu () () Department o Urban and Civil Engineering, Ibaraki University, Japan Abstract A primary concern o the

More information

arxiv: v1 [physics.flu-dyn] 8 May 2014

arxiv: v1 [physics.flu-dyn] 8 May 2014 Energetics of a flui uner the Boussinesq approximation arxiv:1405.1921v1 [physics.flu-yn] 8 May 2014 Kiyoshi Maruyama Department of Earth an Ocean Sciences, National Defense Acaemy, Yokosuka, Kanagawa

More information

A Short Note on Self-Similar Solution to Unconfined Flow in an Aquifer with Accretion

A Short Note on Self-Similar Solution to Unconfined Flow in an Aquifer with Accretion Open Journal o Flui Dynamics, 5, 5, 5-57 Publishe Online March 5 in SciRes. http://www.scirp.org/journal/oj http://x.oi.org/.46/oj.5.57 A Short Note on Sel-Similar Solution to Unconine Flow in an Aquier

More information

12.11 Laplace s Equation in Cylindrical and

12.11 Laplace s Equation in Cylindrical and SEC. 2. Laplace s Equation in Cylinrical an Spherical Coorinates. Potential 593 2. Laplace s Equation in Cylinrical an Spherical Coorinates. Potential One of the most important PDEs in physics an engineering

More information

Thermal runaway during blocking

Thermal runaway during blocking Thermal runaway uring blocking CES_stable CES ICES_stable ICES k 6.5 ma 13 6. 12 5.5 11 5. 1 4.5 9 4. 8 3.5 7 3. 6 2.5 5 2. 4 1.5 3 1. 2.5 1. 6 12 18 24 3 36 s Thermal runaway uring blocking Application

More information

Time-of-Arrival Estimation in Non-Line-Of-Sight Environments

Time-of-Arrival Estimation in Non-Line-Of-Sight Environments 2 Conference on Information Sciences an Systems, The Johns Hopkins University, March 2, 2 Time-of-Arrival Estimation in Non-Line-Of-Sight Environments Sinan Gezici, Hisashi Kobayashi an H. Vincent Poor

More information

An Upper Bound on the Minimum Distance of Serially Concatenated Convolutional Codes

An Upper Bound on the Minimum Distance of Serially Concatenated Convolutional Codes SUBMITTED TO IEEE TRANSACTIONS ON INFORMATION THEORY, FEBRUARY 2004 1 An Upper Boun on the Minimum Distance o Serially Concatenate Convolutional Coes Alberto Perotti, Member, IEEE, an Sergio Beneetto,

More information

A Simple Model for the Calculation of Plasma Impedance in Atmospheric Radio Frequency Discharges

A Simple Model for the Calculation of Plasma Impedance in Atmospheric Radio Frequency Discharges Plasma Science an Technology, Vol.16, No.1, Oct. 214 A Simple Moel for the Calculation of Plasma Impeance in Atmospheric Raio Frequency Discharges GE Lei ( ) an ZHANG Yuantao ( ) Shanong Provincial Key

More information

The Principle of Least Action

The Principle of Least Action Chapter 7. The Principle of Least Action 7.1 Force Methos vs. Energy Methos We have so far stuie two istinct ways of analyzing physics problems: force methos, basically consisting of the application of

More information

INFLUENCE OF SURFACE ROUGHNESS THROUGH A SERIES OF FLOW FACTORS ON THE PERFORMANCE OF A LONGITUDINALLY ROUGH FINITE SLIDER BEARING

INFLUENCE OF SURFACE ROUGHNESS THROUGH A SERIES OF FLOW FACTORS ON THE PERFORMANCE OF A LONGITUDINALLY ROUGH FINITE SLIDER BEARING ANNALS of Faculty Engineering Huneoara International Journal of Engineering Tome XIV [2016] Fascicule 2 [May] ISSN: 1584-2665 [print; online] ISSN: 1584-2673 [CD-Rom; online] a free-accessmultiisciplinarypublication

More information

An Anisotropic Hardening Model for Springback Prediction

An Anisotropic Hardening Model for Springback Prediction An Anisotropic Harening Moel for Springback Preiction Danielle Zeng an Z. Ceric Xia Scientific Research Laboratories For Motor Company Dearborn, MI 48 Abstract. As more Avance High-Strength Steels (AHSS

More information

Finite element analysis of the dynamic behavior of concrete-filled double skin steel tubes (CFDST) under lateral impact with fixed-simply support

Finite element analysis of the dynamic behavior of concrete-filled double skin steel tubes (CFDST) under lateral impact with fixed-simply support Structures Uner Shock an Impact XIII 383 Finite element analysis of the ynamic behavior of concrete-fille ouble skin steel tubes (CFDST) uner lateral impact with fixe-simply support M. Liu & R. Wang Architecture

More information

Design criteria for Fiber Reinforced Rubber Bearings

Design criteria for Fiber Reinforced Rubber Bearings Design criteria or Fiber Reinorced Rubber Bearings J. M. Kelly Earthquake Engineering Research Center University o Caliornia, Berkeley A. Calabrese & G. Serino Department o Structural Engineering University

More information

Crack-tip stress evaluation of multi-scale Griffith crack subjected to

Crack-tip stress evaluation of multi-scale Griffith crack subjected to Crack-tip stress evaluation of multi-scale Griffith crack subjecte to tensile loaing by using periynamics Xiao-Wei Jiang, Hai Wang* School of Aeronautics an Astronautics, Shanghai Jiao Tong University,

More information

Conservation laws a simple application to the telegraph equation

Conservation laws a simple application to the telegraph equation J Comput Electron 2008 7: 47 51 DOI 10.1007/s10825-008-0250-2 Conservation laws a simple application to the telegraph equation Uwe Norbrock Reinhol Kienzler Publishe online: 1 May 2008 Springer Scienceusiness

More information

Chapter 6: Energy-Momentum Tensors

Chapter 6: Energy-Momentum Tensors 49 Chapter 6: Energy-Momentum Tensors This chapter outlines the general theory of energy an momentum conservation in terms of energy-momentum tensors, then applies these ieas to the case of Bohm's moel.

More information

EXPONENTIAL FOURIER INTEGRAL TRANSFORM METHOD FOR STRESS ANALYSIS OF BOUNDARY LOAD ON SOIL

EXPONENTIAL FOURIER INTEGRAL TRANSFORM METHOD FOR STRESS ANALYSIS OF BOUNDARY LOAD ON SOIL Tome XVI [18] Fascicule 3 [August] 1. Charles Chinwuba IKE EXPONENTIAL FOURIER INTEGRAL TRANSFORM METHOD FOR STRESS ANALYSIS OF BOUNDARY LOAD ON SOIL 1. Department of Civil Engineering, Enugu State University

More information

Michaelis-Menten Plot

Michaelis-Menten Plot NAM: HAITI CHMISTRY, SPRING, 17 Circle Section Number: 1 11 amination, April 15, 17 Answer each question in the space provie; use back o page i etra space is neee. Answer questions so the graer can RADILY

More information

Chapter 9 Method of Weighted Residuals

Chapter 9 Method of Weighted Residuals Chapter 9 Metho of Weighte Resiuals 9- Introuction Metho of Weighte Resiuals (MWR) is an approimate technique for solving bounary value problems. It utilizes a trial functions satisfying the prescribe

More information

arxiv: v2 [cond-mat.stat-mech] 11 Nov 2016

arxiv: v2 [cond-mat.stat-mech] 11 Nov 2016 Noname manuscript No. (will be inserte by the eitor) Scaling properties of the number of ranom sequential asorption iterations neee to generate saturate ranom packing arxiv:607.06668v2 [con-mat.stat-mech]

More information

Recommendations: Part 7: Transient Creep for service and accident conditions

Recommendations: Part 7: Transient Creep for service and accident conditions Materials an Structures/Matériaux et Constructions, Vol. 31, June 1998, pp 290-295 RILEM TECHNICAL COMMITTEES RILEM TC 129-MHT: TEST METHODS FOR MECHANICAL PROPERTIES OF CONCRETE AT HIGH TEMPERATURES Recommenations:

More information

Elastic nucleon-deuteron scattering and breakup with chiral forces

Elastic nucleon-deuteron scattering and breakup with chiral forces Elastic nucleon-euteron scattering an breakup with chiral orces H.Witała, J.Golak, R.Skibiński, K.Topolnicki JAGIELLONIAN UNIVERSITY LENPIC Collaboration Jagiellonian University, Kraków Ruhr-Universität,

More information

Basic Strategy for Card-Counters: An Analytic Approach. by N. Richard Werthamer

Basic Strategy for Card-Counters: An Analytic Approach. by N. Richard Werthamer Basic Strategy or Car-Counters: An Analytic Approach by N. Richar Werthamer Abstract We evelop an eicient metho, similar to that o Marcus (7) who calls it Counter Basic Strategy, or casino blackack car

More information

Adhesive Wear Theory of Micromechanical Surface Contact

Adhesive Wear Theory of Micromechanical Surface Contact International Journal Of Computational Engineering esearch ijceronline.com Vol. Issue. hesive Wear Theory of Micromechanical Surface Contact iswajit era Department of Mechanical Engineering National Institute

More information

APPROXIMATE SOLUTION FOR TRANSIENT HEAT TRANSFER IN STATIC TURBULENT HE II. B. Baudouy. CEA/Saclay, DSM/DAPNIA/STCM Gif-sur-Yvette Cedex, France

APPROXIMATE SOLUTION FOR TRANSIENT HEAT TRANSFER IN STATIC TURBULENT HE II. B. Baudouy. CEA/Saclay, DSM/DAPNIA/STCM Gif-sur-Yvette Cedex, France APPROXIMAE SOLUION FOR RANSIEN HEA RANSFER IN SAIC URBULEN HE II B. Bauouy CEA/Saclay, DSM/DAPNIA/SCM 91191 Gif-sur-Yvette Ceex, France ABSRAC Analytical solution in one imension of the heat iffusion equation

More information

DUCTILE FRACTURE PREDICTION IN THIN METALLIC PLATES SUBJECTED TO EXPLOSIONS

DUCTILE FRACTURE PREDICTION IN THIN METALLIC PLATES SUBJECTED TO EXPLOSIONS DUCTLE FRACTURE PREDCTON N TN METALLC PLATES SUBJECTED TO EXPLOSONS Łukasz Pyrzowski 1, Paweł Kłosowski, Krzyszto Woznica 3, Olivier Pennetier 4 1 Gansk University o Technology, Narutowicza str. 11/1,

More information

Transmission Line Matrix (TLM) network analogues of reversible trapping processes Part B: scaling and consistency

Transmission Line Matrix (TLM) network analogues of reversible trapping processes Part B: scaling and consistency Transmission Line Matrix (TLM network analogues of reversible trapping processes Part B: scaling an consistency Donar e Cogan * ANC Eucation, 308-310.A. De Mel Mawatha, Colombo 3, Sri Lanka * onarecogan@gmail.com

More information

Experimental Robustness Study of a Second-Order Sliding Mode Controller

Experimental Robustness Study of a Second-Order Sliding Mode Controller Experimental Robustness Stuy of a Secon-Orer Sliing Moe Controller Anré Blom, Bram e Jager Einhoven University of Technology Department of Mechanical Engineering P.O. Box 513, 5600 MB Einhoven, The Netherlans

More information

Size Effect on Strength of Bimaterial Joints: Computational Approach versus Analysis and Experiment

Size Effect on Strength of Bimaterial Joints: Computational Approach versus Analysis and Experiment Proc. (CD), 12th Intern. Con. on Fracture (ICF-12) (held in Ottawa, Canada, July 2009), Robert Bell, ed., Carleton University, pp. 1 10. Size Eect on Strength o Bimaterial Joints: Computational Approach

More information

THE VAN KAMPEN EXPANSION FOR LINKED DUFFING LINEAR OSCILLATORS EXCITED BY COLORED NOISE

THE VAN KAMPEN EXPANSION FOR LINKED DUFFING LINEAR OSCILLATORS EXCITED BY COLORED NOISE Journal of Soun an Vibration (1996) 191(3), 397 414 THE VAN KAMPEN EXPANSION FOR LINKED DUFFING LINEAR OSCILLATORS EXCITED BY COLORED NOISE E. M. WEINSTEIN Galaxy Scientific Corporation, 2500 English Creek

More information

A Numerical Analysis on a Compact Heat Exchanger in Aluminum Foam

A Numerical Analysis on a Compact Heat Exchanger in Aluminum Foam Journal o Physics: Conerence Series PAPER OPEN ACCESS A Numerical Analysis on a Compact Heat Exchanger in Aluminum Foam To cite this article: B Buonomo et al 2016 J. Phys.: Con. Ser. 745 032141 View the

More information

A SIMPLE ENGINEERING MODEL FOR SPRINKLER SPRAY INTERACTION WITH FIRE PRODUCTS

A SIMPLE ENGINEERING MODEL FOR SPRINKLER SPRAY INTERACTION WITH FIRE PRODUCTS International Journal on Engineering Performance-Base Fire Coes, Volume 4, Number 3, p.95-3, A SIMPLE ENGINEERING MOEL FOR SPRINKLER SPRAY INTERACTION WITH FIRE PROCTS V. Novozhilov School of Mechanical

More information

water adding dye partial mixing homogenization time

water adding dye partial mixing homogenization time iffusion iffusion is a process of mass transport that involves the movement of one atomic species into another. It occurs by ranom atomic jumps from one position to another an takes place in the gaseous,

More information

A 3D SEDIMENT TRANSPORT MODEL FOR COMBINED WAVE-CURRENT FLOWS

A 3D SEDIMENT TRANSPORT MODEL FOR COMBINED WAVE-CURRENT FLOWS A 3D SEDIMENT TRANSPORT MODEL FOR COMBINED WAVE-CURRENT FLOWS Peifeng Ma 1 an Ole Secher Masen Accurate preiction of current velocity an bottom shear stress, which both can be significantly influence by

More information

EVOLUTION OF PARTICLE SIZE DISTRIBUTION IN AIR IN THE RAINFALL PROCESS VIA THE MOMENT METHOD

EVOLUTION OF PARTICLE SIZE DISTRIBUTION IN AIR IN THE RAINFALL PROCESS VIA THE MOMENT METHOD 137 THERMAL SCIENCE, Year 1, Vol. 16, No. 5, pp. 137-1376 EVOLUTION OF PARTICLE SIZE DISTRIBUTION IN AIR IN THE RAINFALL PROCESS VIA THE MOMENT METHOD by Fu-Jun GAN a an Jian-Zhong LIN a,b * a Department

More information

Assessment of the Buckling Behavior of Square Composite Plates with Circular Cutout Subjected to In-Plane Shear

Assessment of the Buckling Behavior of Square Composite Plates with Circular Cutout Subjected to In-Plane Shear Assessment of the Buckling Behavior of Square Composite Plates with Circular Cutout Sujecte to In-Plane Shear Husam Al Qalan 1)*, Hasan Katkhua 1) an Hazim Dwairi 1) 1) Assistant Professor, Civil Engineering

More information

Non-Superconducting Fault Current Limiter

Non-Superconducting Fault Current Limiter Proceeings o the 5th WEA International Conerence on Applications o Electrical Engineering, Prague, Czech Republic, March -, 6 (pp-7) Non-uperconucting Fault Current Limiter A. CHARMIN*, M. TARAFDAR HAQUE,

More information

Math 342 Partial Differential Equations «Viktor Grigoryan

Math 342 Partial Differential Equations «Viktor Grigoryan Math 342 Partial Differential Equations «Viktor Grigoryan 6 Wave equation: solution In this lecture we will solve the wave equation on the entire real line x R. This correspons to a string of infinite

More information

AXIALLY LOADED FRP CONFINED REINFORCED CONCRETE CROSS-SECTIONS

AXIALLY LOADED FRP CONFINED REINFORCED CONCRETE CROSS-SECTIONS AXIALLY LOADED FRP CONFINED REINFORCED CONCRETE CROSS-SECTIONS Bernát Csuka Budapest University o Technology and Economics Department o Mechanics Materials and Structures Supervisor: László P. Kollár 1.

More information

Situation awareness of power system based on static voltage security region

Situation awareness of power system based on static voltage security region The 6th International Conference on Renewable Power Generation (RPG) 19 20 October 2017 Situation awareness of power system base on static voltage security region Fei Xiao, Zi-Qing Jiang, Qian Ai, Ran

More information

PCCP PAPER. 1 Introduction. A. Nenning,* A. K. Opitz, T. M. Huber and J. Fleig. View Article Online View Journal View Issue

PCCP PAPER. 1 Introduction. A. Nenning,* A. K. Opitz, T. M. Huber and J. Fleig. View Article Online View Journal View Issue PAPER View Article Online View Journal View Issue Cite this: Phys. Chem. Chem. Phys., 2014, 16, 22321 Receive 4th June 2014, Accepte 3r September 2014 DOI: 10.1039/c4cp02467b www.rsc.org/pccp 1 Introuction

More information

IERCU. Institute of Economic Research, Chuo University 50th Anniversary Special Issues. Discussion Paper No.210

IERCU. Institute of Economic Research, Chuo University 50th Anniversary Special Issues. Discussion Paper No.210 IERCU Institute of Economic Research, Chuo University 50th Anniversary Special Issues Discussion Paper No.210 Discrete an Continuous Dynamics in Nonlinear Monopolies Akio Matsumoto Chuo University Ferenc

More information

To understand how scrubbers work, we must first define some terms.

To understand how scrubbers work, we must first define some terms. SRUBBERS FOR PARTIE OETION Backgroun To unerstan how scrubbers work, we must first efine some terms. Single roplet efficiency, η, is similar to single fiber efficiency. It is the fraction of particles

More information

Prediction of ECC tensile stress-strain curves based on modified fiber bridging relations considering fiber distribution characteristics

Prediction of ECC tensile stress-strain curves based on modified fiber bridging relations considering fiber distribution characteristics Computers and Concrete, Vol. 7, No. 5 (010) 455-468 455 Prediction of ECC tensile stress-strain curves based on modified fiber bridging relations considering fiber distribution characteristics Bang Yeon

More information

3.5 Analysis of Members under Flexure (Part IV)

3.5 Analysis of Members under Flexure (Part IV) 3.5 Analysis o Members under Flexure (Part IV) This section covers the ollowing topics. Analysis o a Flanged Section 3.5.1 Analysis o a Flanged Section Introduction A beam can have langes or lexural eiciency.

More information

The effect of nonvertical shear on turbulence in a stably stratified medium

The effect of nonvertical shear on turbulence in a stably stratified medium The effect of nonvertical shear on turbulence in a stably stratifie meium Frank G. Jacobitz an Sutanu Sarkar Citation: Physics of Fluis (1994-present) 10, 1158 (1998); oi: 10.1063/1.869640 View online:

More information

Assignment 1. g i (x 1,..., x n ) dx i = 0. i=1

Assignment 1. g i (x 1,..., x n ) dx i = 0. i=1 Assignment 1 Golstein 1.4 The equations of motion for the rolling isk are special cases of general linear ifferential equations of constraint of the form g i (x 1,..., x n x i = 0. i=1 A constraint conition

More information

arxiv:physics/ v2 [physics.ed-ph] 23 Sep 2003

arxiv:physics/ v2 [physics.ed-ph] 23 Sep 2003 Mass reistribution in variable mass systems Célia A. e Sousa an Vítor H. Rorigues Departamento e Física a Universiae e Coimbra, P-3004-516 Coimbra, Portugal arxiv:physics/0211075v2 [physics.e-ph] 23 Sep

More information

TOWARDS THERMOELASTICITY OF FRACTAL MEDIA

TOWARDS THERMOELASTICITY OF FRACTAL MEDIA ownloae By: [University of Illinois] At: 21:04 17 August 2007 Journal of Thermal Stresses, 30: 889 896, 2007 Copyright Taylor & Francis Group, LLC ISSN: 0149-5739 print/1521-074x online OI: 10.1080/01495730701495618

More information

This section outlines the methodology used to calculate the wave load and wave wind load values.

This section outlines the methodology used to calculate the wave load and wave wind load values. COMPUTERS AND STRUCTURES, INC., JUNE 2014 AUTOMATIC WAVE LOADS TECHNICAL NOTE CALCULATION O WAVE LOAD VALUES This section outlines the methoology use to calculate the wave loa an wave win loa values. Overview

More information

P. A. Martin b) Department of Mathematics, University of Manchester, Manchester M13 9PL, United Kingdom

P. A. Martin b) Department of Mathematics, University of Manchester, Manchester M13 9PL, United Kingdom Time-harmonic torsional waves in a composite cyliner with an imperfect interface J. R. Berger a) Division of Engineering, Colorao School of Mines, Golen, Colorao 80401 P. A. Martin b) Department of Mathematics,

More information

S. Srinivasan, Technip Offshore, Inc., Houston, TX

S. Srinivasan, Technip Offshore, Inc., Houston, TX 9 th ASCE Specialty Conerence on Probabilistic Mechanics and Structural Reliability PROBABILISTIC FAILURE PREDICTION OF FILAMENT-WOUND GLASS-FIBER Abstract REINFORCED COMPOSITE TUBES UNDER BIAXIAL LOADING

More information

Model for Dopant and Impurity Segregation During Vapor Phase Growth

Model for Dopant and Impurity Segregation During Vapor Phase Growth Mat. Res. Soc. Symp. Proc. Vol. 648, P3.11.1-7 2001 Materials Research Society Moel for Dopant an Impurity Segregation During Vapor Phase Growth Craig B. Arnol an Michael J. Aziz Division of Engineering

More information

Survey Sampling. 1 Design-based Inference. Kosuke Imai Department of Politics, Princeton University. February 19, 2013

Survey Sampling. 1 Design-based Inference. Kosuke Imai Department of Politics, Princeton University. February 19, 2013 Survey Sampling Kosuke Imai Department of Politics, Princeton University February 19, 2013 Survey sampling is one of the most commonly use ata collection methos for social scientists. We begin by escribing

More information

Module 5 Couplings. Version 2 ME, IIT Kharagpur

Module 5 Couplings. Version 2 ME, IIT Kharagpur Moule 5 Couplings Version ME, IIT Kharagpur Lesson Design proceures for rigi an flexible rubber-bushe couplings Version ME, IIT Kharagpur Instructional Objectives At the en of this lesson, the stuents

More information