Analysis of Interfacial Damage and Debonding Life Estimation of Cold Forging Tool Coated with Hard Film

Size: px
Start display at page:

Download "Analysis of Interfacial Damage and Debonding Life Estimation of Cold Forging Tool Coated with Hard Film"

Transcription

1 Materials Trasactios, Vol. 45, No. 9 (24) pp to 2837 #24 The Japa Society for Techology of Plasticity Aalysis of Iterfacial Damage ad Debodig Life Estimatio of Cold Forgig Coated with Hard Kuio Hayakawa 1, Tamotsu Nakamura 1, Shigekazu Taaka 1 ad Kuiaki Harada 2 1 Departmet of Mechaical Egieerig, Shizuoka Uiversity, Hamamatsu , Japa 2 NTN Corporatio, Iwata , Japa The behavior of the iterface of a cold forgig tool coated with a hard film is aalyzed by the fiite elemet method (FEM). The mechaical properties of the iterface betwee a hard film ad a tool material are modeled uder the framework of the irreversible thermodyamic costitutive equatio theory with cotiuum damage mechaics. The iterface damage variable ad its evolutio law are itroduced, which express the degradatio of the iterface. The hemispherical tool of SUJ2 i JIS coated with a electroplated coatig of chromium idets ito a cylidrical workpiece of S25C i JIS, which is performed to examie the validity of the proposed model. I FE aalysis, the proposed costitutive equatios are implemeted by the oliear sprigs betwee opposite odes of the iterface. Calculated results show that the maximum iterfacial ormal displacemet is observed close to the frot of the lip of idetatio, whereas maximum iterfacial tagetial displacemet is observed close to the back of the lip of idetatio. The regio where debodig will first occur as idicated by the calculatio results, is similar to the regio of the behavior durig a actual forgig operatio. Fially, a method of evaluatig the iterfacial debodig life of coated tools is proposed. (Received October 6, 23; Accepted Jue 1, 24) Keywords: cold forgig, hard-film-coated tool, iterface, fiite elemet method, oliear sprig model, iterface damage mechaics, iterfacial separatio, debodig life 1. Itroductio Recetly, forgig tools coated with hard films such as electroplated coatigs of chromium, titaium itride (TiN) ad titaium carbide (TiC), are used with the aim of icreasig the service life of tools by about 2 to 1 times. 1 3) I the practical use of coated tools, however, we ofte ecouter debodig at the iterface of the film ad the tool material. It is importat to grasp the behavior of the iterface of a hard film ad a tool material to accurately evaluate iterfacial stregth. Numerical simulatios of the mechaical behavior of materials coated with a hard film have bee performed. 4,5) The objects of the simulatios are flat hard films o a flat substrate ideted with a ideter similar to that used i a scratch test or a idetatio test. 6) Moreover, complete adhesio at the iterface is postulated. However, the deformatio of the film i a scratch test or a idetatio test does ot correspod to that of a actual hard film coated o forgig tools sice the film itself is ideted. Furthermore, because complete adhesio at the iterface is postulated, the deterioratio of the iterface is ot cosidered. I the preset paper, the fiite elemet aalysis of the mechaical behavior of the hard film - tool material iterface of forgig tools uder coditios closer to actual forgig situatios tha those for covetioal scratch or idetatio tests was performed by iterface damage mechaics (IDM). First, the examiatio model of idetatio of hemispherical coated tool was proposed. Next, after itroducig a oliear iterfacial sprig model ad defiig the iterfacial damage variable, the costitutive ad damage evolutio equatios at the iterface betwee the hard film ad the tool material are formulated uder the framework of the irreversible thermodyamic costitutive theory. 8,9,11) To apply the above equatios to the umerical aalysis of the iterface by FEM, the cocept of the oliear iterfacial sprig model is employed. 7) O the basis of the calculated results o the deformatio behavior of the iterface, the probable debodig mechaisms of the preset forgig tool coated with the hard film are discussed. Fially, a method of evaluatig the iterfacial debodig life is proposed. Moreover, the calculated results o the debodig life of the preset forgig tool are show. 2. Examied Model for Iterfacial Behavior of Forgig Coated with Hard I the scratch or idetatio tests, the hard film o the substrate is ideted by a slider or a ideter. The, the aspect of the deformatio of the film does ot correspod to that of actual forgig tools. It is favorable to evaluate the behavior of the iterface at coditios closer to those of actual forgig tha to coditios of these tests. I the preset study, we aalyze the behavior of the iterface betwee the hard film ad the tool material by the examiatio model as show i Fig. 1. I this model, the hemispherical coated tool of 5 mm radius is ideted to the R Thickess t Fig. 1 Schematics of examiatio model of idetatio of hemispherical tool coated with hard film.

2 Aalysis of Iterfacial Damage ad Debodig Life Estimatio of Cold Forgig Coated with Hard 2833 Iterfacial regio h Microdefects workpiece of 2 mm i diameter ad 1 mm i height. The deformatio of the bottom ad the periphery of the workpiece are costraied. The proposed model is itroduced as the simplificatio of the actual forgig. The effect of the plastic flow of the workpiece ad the frictio coditio betwee the film ad the workpiece o the behavior of the iterface ca be take ito accout more effectively. 3. Modelig of Iterface by Noliear Iterfacial Sprig Model 3.1 Noliear iterfacial sprig model Let us postulate the iterfacial behavior that occurs i the iterfacial regio show i Fig. 2(a). I the iterfacial regio, microdefects develop whe the exteral load is subjected, ad adhesio stregth deteriorates. I the preset study, o calculatig the iterfacial behavior by FEM, the adhesio of the iterfacial regio betwee the hard film ad the tool material is modeled usig distributed oliear sprigs proposed by Ma ad Kishimoto, 7) as show i Fig. 2(b). I this model, the displacemet of the iterfacial regio ca be expressed by the differece i displacemet betwee the iitial poits o the opposite sides of the film ad the tool material P ad P across the iterfacial regio, as show i Fig. 2(b). Whe exteral load is applied, poit P will move to poit Q, ad these poits will produce the iterfacial displacemet U ¼ P Q. The iterfacial displacemet vector U ca be expressed by the compoets of the ormal directio, tagetial directio t ad aother directio b accordig to a right-had Cartesia coordiate system i Fig. 2, expressed as U, U t ad U b. Moreover, we defie the iterfacial relative displacemet u as u ¼ðu ; u t ; u b Þ T ¼ 1 h ðu ; U t ; U b Þ T where h is the thickess of the iterfacial regio. Moreover, the superscript T meas traspositio. 3.2 Iterfacial damage variable Whe the iterface is subject to a exteral load, the iterface deteriorates gradually because of the iitiatio ad growth of microscopic cracks ad cavities. Figure 3 shows the odamaged ad damaged states of the iterface. The elemetal area at the udamaged state da is U P U t Q Fig. 2 Schematics of Iterfacial regio with microdefects ad equivalet iterfacial sprig model. (a) Iterfacial regio (b) Equivalet sprig model P b t ð1þ da No-damaged state of iterface deteriorated by the exteral load, ad the elemetal area da becomes d ~a. The, the iterfacial damage variable D is defied as d ~a ¼ð1 D Þda ð3þ where is the uit ormal vector to the iterface. I the preset study, we express the deterioratio of the iterface usig the iterfacial damage variable D. 4. Modelig of Iterfacial Costitutive ad Damage Evolutio Equatios 4.1 Iterfacial stress iterfacial relative displacemet relatioship with iterfacial damage The irreversible thermodyamic costitutive equatio theory 8) is oe of the most ratioal frameworks for the formulatio of the costitutive equatios with irreversible chages of iteral states such as damage ad plastic deformatio. I the preset study, we apply this theory to the formulatio of the costitutive ad damage evolutio equatios of the iterface. If the thickess h is small ad the strai throughout the thickess is uiform, iterfacial strai tesor " is defied as " ¼ 1 2 ðu þ u Þ ð3þ Whe we restrict the adiabatic, costat temperature coditio, Clausius Duhem iequality i the iterfacial regio ca be express as : _" _ ð4þ where ad are the iterfacial stress ad Helmholtz free eergy i the iterfacial regio. Helmholtz eergy is a fuctio of iterfacial relative displacemet u ad iterfacial damage variable D. By use of the iterfacial stress vector T ¼, eq. (4) ca be rewritte as To hold for eq. _D ð5þ For the Helmholtz free eergy, we defie da Damaged state of iterface Fig. 3 Schematics of o-damaged ad damaged states of elemetal area i hard film-tool material iterface. ð6þ

3 2834 K. Hayakawa, T. Nakamura, S. Taaka ad K. Harada ¼ 1 2 ð1 D Þ u K u where K is the iitial iterfacial rigidity tesor at the itact state. Here, we will adopt followig simple form: K 1 t K ¼ B Kb A ð8þ K By usig eqs. (6) ad (7), the relatioship betwee iterfacial stress vector ad iterfacial relative displacemet ca be obtaied as T ¼ð1 D ÞK u ð9þ O the other had, we defie the thermodyamic cojugate force of the iterfacial damage variable R ¼ 1 2 ðk t u2 t þ K b u2 b þ K u2 Þ ð7þ ð1þ 4.2 Evolutio equatio of iterfacial damage variable Uder the restrictio of the framework of the irreversible thermodyamic costitutive equatio, the evolutio equatio of the iterfacial damage variable _D ca be derived usig the ormality rule, so that _D ca develop to the outer ormal directio of the damage dissipatio potetial defied i the space of the cojugate force of the iterfacial damage F d ðr Þ _D ¼ _ d ð _ Þ where _ is a udetermied multiplier that prescribes the magitude of damage developmet. I the preset study, we postulate a simple form for the damage dissipatio potetial F d ðr Þ as F d ðr Þ¼R ð12þ For the multiplier _, we give the followig form so that the iterfacial fatigue behavior uder the cyclic loadig ca be expressed as _ ¼ R R D 1 h _R i ð13þ S S Usig eqs. (11) ad (12), eq. (11) is expressed as _D ¼ R R D 1 h _R i ð14þ S S where, S ad R are material costats. Moreover, the brackets hi are the Macauley brackets. As observed i eq. (13), the iterfacial damage variable D develops whe the cojugate force of the iterfacial damage R exceeds the threshold value R ad the derivative of R with respect to time icreases. Debodig is assumed to occur whe the iterfacial damage variable D reaches the threshold damage variable D cr. Whe debodig occurs, the iterfacial rigidity is reduced to zero. Adhesio stregth ca be evaluated by iterfacial stress i eq. (9). 4.3 Iitial iterfacial rigidity The deformatio of the iterface regio depeds o the hard film ad tool materials. The material that ca deform more easily tha others will give a stroger effect o the deformatio property of the iterfacial regio. Therefore, the iitial iterfacial rigidities Kb, K t ad K are postulated to be expressed as Kt ¼ Kb ¼ 2G f G t =ðg f þ G t Þ ð15þ K ¼ 2E f E t =ðe f þ E t Þ ð16þ where E f ad E t are Youg s moduli of the hard film ad the tool materials, ad G f ad G t are the their shear moduli, respectively. 5. Fiite Elemet Aalysis of Iterfacial Behavior Figure 4 shows the discretized model, boudary coditios ad the equivalet iterfacial sprig of the iterface for aalysis. I the preset study, the axisymmetric model is used. The geometries of the hemispherical coated tool ad the workpiece are the same as show i Fig. 1. For the boudary coditios of the hemispherical tool, the odes o the axisymmetric axis are restricted o the x-coordiate, ad the idetatio depth s o the y-coordiate is give to the odes of the upper edge of the tool. The maximum idetatio depth s max is 4.5 mm. For the parameter idicatig the positio of the iterface, we adopt the agle ( 9 ), as show i Fig. 4. For the boudary coditio of the workpiece, we give the costrait of the x-directio to the odes o the axisymmetric axis ad the outer edge, ad the costrait of the y-directio to the odes o the bottom edge. The commercial fiite elemet code MSC. Marc2 is used i the preset study. The oliear iterfacial sprigs betwee the two material poits o opposite sides of the iterface, as show i Fig. 4, are implemeted usig the user subrouties equipped i the code. For fiite elemets, we use axisymmetric 4-ode isoparametric quadrilateral elemets. The umbers of the elemets x Idetatio y depth s θ Coated tool a =.13mm Equivalet iterfacial sprig Fig. 4 Aalyzed model of idetatio of hemispherical tool coated with hard film. t

4 Aalysis of Iterfacial Damage ad Debodig Life Estimatio of Cold Forgig Coated with Hard 2835 Table 1 Elastic ad plastic properties materials used i FEM. Material Youg s modulus /GPa Poisso s ratio of the hard film, the tool material ad the workpiece are 18, 175 ad 64, respectively. Characteristic legth h is defied as h ¼ :1 mm. The depedecy of h o the umerical results is ot cosidered here because we aim at the validatio of the proposed costitutive equatios of the iterface. For frictio coditios, a coulomb frictio coefficiet of ¼ :3 is employed by assumig coditios of o lubricatio to promote the developmet of the iterfacial damage. Table 1 shows the elastic properties of the materials used i the aalysis. As materials, we employ a chromium electroplatig, a ball-bearig steel SUJ2 (i JIS) ad a lowcarbo steel S25C (i JIS) for the hard film, tool ad the workpiece, respectively. Usig the values i Table 1, the iterfacial rigidity of eqs. (17) ad (18) ca be calculated to be K ¼ 237: GPa ad K t ¼ Kb ¼ 93:92 GPa. Figure 5 shows the flow stress - plastic strai relatioship of the tool material ad the workpiece. Although the material costats D, S D ad R of eq. (15) should be determied experimetally, we employ the followig values, for the purpose of determiig the characteristics of the proposed damage evolutio equatio qualitatively: ¼ 2:, S ¼ 5: ad R ¼ :3. As observed i eq. (14), the evolutio of the iterfacial damage _D is attributed to a form proportioal to the term ðr R Þ=S for ¼ 2:. Furthermore, S ca determie the magitude of iterfacial damage evolutio. However, the ifluece of the value of S to the qualitative characteristics of the iterfacial damage is small. 6. Results ad Discussio Yield stress /MPa Behavior of iterface by iterfacial sprig model Figure 6 shows the distributio of equivalet stress i the hard film, tool material ad workpiece at S ¼ 1 mm. The Equivalet stress σ eq / MPa Equivalet plastic strai ε eq p 1. Fig. 5 Flow stress - plastic strai curves of tool ad workpiece materials used for FEM. Equivalet stress / MPa Fig Idetatio depth s = 1 mm Distributio of equivalet stress at idetatio depth of s ¼ 1 mm. plastic deformatio of the tool material does ot yield. I this figure, high equivalet stress i the tool material ad hard film is observed i the viciity of the lip of idetatio. Compressive stress is take ito accout i the calculatio of equivalet stress show i Fig. 6. For debodig, however, the effect of compressive stress ca be cosidered ot to be high. Figures 7(a) ad 7(b) show the relatioship betwee the iterfacial relative displacemets u ad u t, ad the positio of the iterface with the parameter of idetatio depth s.i these figures, the positios of the cotact edge betwee the tool ad workpiece at each idetatio depth are idicated by solid circles. As observed i Fig. 7(a), the iterfacial relative displacemet i the ormal directio u takes a maximum value at ¼ 42. This positio is close to the back of the cotact edge. A similar relatioship ca be observed at ay idetatio depth. O the other had, the iterfacial relative displacemet o tagetial directio u t is distributed moderately. At ay idetatio depths, the maximum u t is i frot of the cotact edge. Let us ivestigate i more detail the positio of the maximum iterfacial relative displacemets u ad u t. Figure 8 shows the positio of the maximum iterfacial displacemet i the ormal ad tagetial directios u ;max ad u t;max at s ¼ 2 mm. As metioed above, u ;max is located at the back of the cotact edge, whereas u t;max is i frot of the cotact edge as observed i Figs. 7(a) ad 7(b). 6.2 Developmet of iterfacial damage variable Figure 9 shows the chage i the iterfacial relative displacemets u ad u t, ad the iterfacial damage variable D at the idetatio depth s at ¼ 45. The iterfacial relative displacemets u ad u t are plotted as left ordiate, whereas the iterfacial damage variable D is plotted as right ordiate. We ca observe the rapid icrease i u at s = approximately 1.15 mm. The iterfacial relative displacemet i the tagetial directio u t starts icreasig at

5 2836 K. Hayakawa, T. Nakamura, S. Taaka ad K. Harada Relative displacemet u / s = 1. mm s = 2. mm s = 3. mm s = 4. mm Positio of cotact edge Positio of iterface θ / (a) 8 Relative displacemet u, u t / θ = 45 Iterfacial damage D u t (Tagetial) u (Normal) Idetatio depth s / mm Fig. 9 Chage i relative displacemet u, u t ad iterfacial damage variable D to idetatio depth s. at iterfacial positio ¼ Iterfacial damage D / 1-6 Relative displacemet u t / s = 1. mm s = 2. mm s = 3. mm s = 4. mm Positio of cotact edge Positio of iterface θ / (b) B 4 5 A Fig. 7 Effects of idetatio depth o distributio of iterfacial relative displacemets u ad u t. (a) Normal directio (b) Tagetial directio Fig. 1 Schematic of debodig mechaism of iterface betwee tool ad workpiece..2mm u t,max u,max mm Fig. 8 Positio of maximum relative displacemets u ;max ad u t;max at idetatio depth s ¼ 2 mm. s = approximately 1.5 mm, ad it takes the maximum at s = approximately 1.4 mm. Iterfacial damage variable D also icreases as the iterfacial relative displacemets u ad u t icrease. However, the developmet of D ceases after the iterfacial relative displacemets u ad u t stop icreasig. 6.3 Mechaism of debodig of hard film Figure 1 shows the schematic of the iterfacial stresses subjected i the viciity of the lip of the idetatio as well as the iterfacial displacemets obtaied by the calculatio. Let us discuss the mechaisms of the debodig of the coated tool from this figure. The hard film is subjected to stress ` ad its reactio stress. Uder such a coditio, the hard film is subjected to compressio by the force of idetatio À. Cosequetly, the hard film will be i the state of buckig. The mechaism cosidered above ca be observed i the practical cold forgig processes usig tools coated with the hard film, particulary a puch. Debodig from the straight lad of the puch ca be ofte experieced. I these cases, the adhesio stregth is smaller tha the stregths of the tool material ad hard film, ad thus debodig by the mechaism cosidered above will occur. The mechaism of debodig by the bucklig of the film has bee reported by Strawbridge ad Evas. 1) They ivestigated such a mechaism usig a flat film o a flat substrate. I the preset study, it is revealed that a similar mechaism may be a cause of the debodig of the iterface of a practical forgig tool with more complicated geometry, i case that the coditios regardig the mechaical properties of the tool material ad hard film are similar to those i the preset calculatio. The iterface is also subject to iterfacial shear stress by frictioal force ˆ. The, debodig by shear ca occur at the

6 Aalysis of Iterfacial Damage ad Debodig Life Estimatio of Cold Forgig Coated with Hard 2837 positio of the maximum iterfacial relative displacemet u t;max, poit B, i Fig. 1. I practical cases, debodig by shear at poit B will be restricted, as large hydrostatic pressure is subjected to poit B. 7. Evaluatio of Iterfacial Debodig Life of Coated Forgig I covetioal forgig, the debodig of the hard film of a tool occurs after some umber of cycles. The cotiuum damage mechaics may be used to evaluate the above pheomea by calculatig iterfacial damage. I the preset chapter, we attempt to evaluate the iterfacial debodig life of a forgig tool coated with a hard film uder the cyclic idetatio show i Fig. 1. A idetatio depth s of 3 mm is employed. The cyclic idetatio is calculated by FEM as follows: First, the sigle idetatio is calculated. Next, the pressure agaist the hard film at the istataeous idetatio depth s obtaied i the calculatio is coverted to the odal exteral force. The, the calculatio of the cyclic idetatio is simulated usig the coverted odal exteral forces as boudary coditios applied to the surface of the coated tool. Let us ote that calculatio after the iterfacial debodig is ot performed i the preset calculatio, sice the evaluatio of chages i the iterfacial damage variable ad iterfacial relative displacemet util the debodig is maily focused o. Figure 11 shows the distributio of the accumulated iterfacial damage variable D util the cyclic idetatio of N ¼ 22. The maximum accumulated damage ca be observed at a iterface positio of ¼ 28. From this result, the iterfacial debodig is cosidered to occur at this positio. Figure 12 shows chages i the ormal relative displacemet u ad iterfacial damage variable D to umber of idetatio at a iterface positio of ¼ 28. As umber of idetatio icreases, iterfacial damage variable ad ormal relative displacemet icrease. Moreover, as observed i Iterfacial damage D Positio of iterface θ / N = 5 N = 1 N = 15 N = 2 N = 22 Fig. 11 Effects of umber of idetatio N o distributio of iterfacial damage variable D. 8 Relative displacemet u / eq. (1), iterfacial damage variable decreases iterfacial rigidity K. 8. Coclusios I the preset paper, the FE aalysis of the iterfacial behavior of the idetatio of a hemispherical tool coated with a chromium electroplatig ito a workpiece is performed i order to elucidate mechaical behavior of the iterface of the cold forgig tool. The results of the FE aalysis show that the iterfacial relative displacemet i the ormal directio takes its maximum close to the back of the lip of the idetatio at ay idetatio depths. O the other had, the maximum iterfacial relative displacemet i the tagetial directio is located i frot of the edge of the idetatio. The possible mechaisms of the debodig of the hard film are iterfacial displacemet i ormal ad tagetial directios. Fially, iterfacial debodig life is evaluated. The aspect of the local accumulatio of iterfacial damage variables with the icrease i umber of idetatio ca be described. REFERENCES θ = 28 Iterfacial damage D Normal relative displacemet u Number of idetatio N Fig. 12 Chage i relative displacemet u ad iterfacial damage variable D to umber of idetatio N at ¼ 28. 1) Y. C. Lee ad F. K. Che: J. Mat. Process. Techol. 15 (2) ) D. C. Ko, D. H. Kim ad B. M. Kim: Wear 252 (22) ) Ch. Hisel, M. Celeghili, U. Egel ad M. Geiger: Proc. 6th ICTP 1 (1999) ) E. R. Kral, K. Komvopoulos ad D. B. Bogy: J. Appl. Mech., Tras. ASME 62 (1995) ) N. Naka, K. Yamamoto, Y. Isoo, T. Taaka ad N. Terayama: J. JSME, A 66 (2) (i Japaese). 6) P. Bejami ad C. Weaver: Proc. R. Soc., Ser. A 254 (196) ) F. S. Ma ad K. Kishimoto: JSME It. J., Ser. A 39 (1996) ) J. Lemaitre ad J. L. Chaboche: Mechaics of Solid Materials, (Cambridge Uiversity Press, 199). 9) K. Hayakawa ad S. Murakami: It. J. Damage Mech. 6 (1997) ) A. Strawbridge ad H. E. Evas: Egg. Failure Aal. 2 (1995) ) K. Hayakawa, T. Nakamura ad S. Taaka: J. JSTP 43 (22) (i Japaese) Iterfacial Damage D

EVALUATION OF GLASS FIBER/EPOXY INTERFACIAL STRENGTH BY THE CRUCIFORM SPECIMEN METHOD

EVALUATION OF GLASS FIBER/EPOXY INTERFACIAL STRENGTH BY THE CRUCIFORM SPECIMEN METHOD EVALUATION OF GLASS FIBER/EPOX INTERFACIAL STRENGTH B THE CRUCIFORM SPECIMEN METHOD Ju KOANAGI, Hajime KATO, Akihiro KASHIMA, uichi IGARASHI, Keichi WATANABE 3, Ichiro UENO 4 ad Shiji OGIHARA 4 Istitute

More information

Analysis of composites with multiple rigid-line reinforcements by the BEM

Analysis of composites with multiple rigid-line reinforcements by the BEM Aalysis of composites with multiple rigid-lie reiforcemets by the BEM Piotr Fedeliski* Departmet of Stregth of Materials ad Computatioal Mechaics, Silesia Uiversity of Techology ul. Koarskiego 18A, 44-100

More information

MATHEMATICAL MODELLING OF ARCH FORMATION IN GRANULAR MATERIALS

MATHEMATICAL MODELLING OF ARCH FORMATION IN GRANULAR MATERIALS 6 th INTERNATIONAL MULTIDISCIPLINARY CONFERENCE MATHEMATICAL MODELLING OF ARCH FORMATION IN GRANULAR MATERIALS Istva eppler SZIE Faculty of Mechaics, H-2103 Gödöllő Páter. 1., Hugary Abstract: The mathematical

More information

Principle Of Superposition

Principle Of Superposition ecture 5: PREIMINRY CONCEP O RUCUR NYI Priciple Of uperpositio Mathematically, the priciple of superpositio is stated as ( a ) G( a ) G( ) G a a or for a liear structural system, the respose at a give

More information

A CONFINEMENT MODEL OF HIGH STRENGTH CONCRETE

A CONFINEMENT MODEL OF HIGH STRENGTH CONCRETE 3 th World Coferece o Earthquake Egieerig Vacouver, B.C., Caada August -6, 24 Paper No. 873 A CONFINEMENT MODEL OF HIGH STRENGTH CONCRETE Nobutaka NAKAZAWA, Kazuhiko KAWASHIMA 2, Gakuho WATANABE 3, Ju-ichi

More information

A statistical method to determine sample size to estimate characteristic value of soil parameters

A statistical method to determine sample size to estimate characteristic value of soil parameters A statistical method to determie sample size to estimate characteristic value of soil parameters Y. Hojo, B. Setiawa 2 ad M. Suzuki 3 Abstract Sample size is a importat factor to be cosidered i determiig

More information

DYNAMIC ANALYSIS OF BEAM-LIKE STRUCTURES SUBJECT TO MOVING LOADS

DYNAMIC ANALYSIS OF BEAM-LIKE STRUCTURES SUBJECT TO MOVING LOADS DYNAMIC ANALYSIS OF BEAM-LIKE STRUCTURES SUBJECT TO MOVING LOADS Ivaa Štimac 1, Ivica Kožar 1 M.Sc,Assistat, Ph.D. Professor 1, Faculty of Civil Egieerig, Uiverity of Rieka, Croatia INTRODUCTION The vehicle-iduced

More information

Failure analysis of single-bolted joint for lightweight composite laminates and metal plate

Failure analysis of single-bolted joint for lightweight composite laminates and metal plate IOP Coferece Series: Materials Sciece ad Egieerig PAPER OPEN ACCESS Failure aalysis of sigle-bolted joit for lightweight composite lamiates ad metal plate To cite this article: Lijie Li et al 018 IOP Cof.

More information

Damped Vibration of a Non-prismatic Beam with a Rotational Spring

Damped Vibration of a Non-prismatic Beam with a Rotational Spring Vibratios i Physical Systems Vol.6 (0) Damped Vibratio of a No-prismatic Beam with a Rotatioal Sprig Wojciech SOCHACK stitute of Mechaics ad Fudametals of Machiery Desig Uiversity of Techology, Czestochowa,

More information

a b c d e f g h Supplementary Information

a b c d e f g h Supplementary Information Supplemetary Iformatio a b c d e f g h Supplemetary Figure S STM images show that Dark patters are frequetly preset ad ted to accumulate. (a) mv, pa, m ; (b) mv, pa, m ; (c) mv, pa, m ; (d) mv, pa, m ;

More information

Numerical Simulation of Thermomechanical Problems in Applied Mechanics: Application to Solidification Problem

Numerical Simulation of Thermomechanical Problems in Applied Mechanics: Application to Solidification Problem Leoardo Joural of Scieces ISSN 1583-0233 Issue 9, July-December 2006 p. 25-32 Numerical Simulatio of Thermomechaical Problems i Applied Mechaics: Applicatio to Solidificatio Problem Vicet Obiajulu OGWUAGWU

More information

APPENDIX A EARLY MODELS OF OXIDE CMP

APPENDIX A EARLY MODELS OF OXIDE CMP APPENDIX A EALY MODELS OF OXIDE CMP Over the past decade ad a half several process models have bee proposed to elucidate the mechaism ad material removal rate i CMP. Each model addresses a specific aspect

More information

ANALYSIS OF STEADY WEAR PROCESSES FOR PERIODIC SLIDING

ANALYSIS OF STEADY WEAR PROCESSES FOR PERIODIC SLIDING Joural of Computatioal ad Applied Mechaics, Vol. 1., No. 2., (215), pp. 231 268 ANALYSIS OF STEADY WEAR PROCESSES FOR PERIODIC SLIDING Istvá Páczelt Istitute of Applied Mechaics, Uiversity of Miskolc H-3515

More information

NUCLEATION 7.1 INTRODUCTION 7.2 HOMOGENEOUS NUCLEATION Embryos and nuclei CHAPTER 7

NUCLEATION 7.1 INTRODUCTION 7.2 HOMOGENEOUS NUCLEATION Embryos and nuclei CHAPTER 7 CHAPER 7 NUCLEAION 7.1 INRODUCION I this text, we focus our attetio o crystallie solids that form from the melt. he process begis with the creatio of a cluster of atoms of crystallie structure, which may

More information

SECTION 2 Electrostatics

SECTION 2 Electrostatics SECTION Electrostatics This sectio, based o Chapter of Griffiths, covers effects of electric fields ad forces i static (timeidepedet) situatios. The topics are: Electric field Gauss s Law Electric potetial

More information

Similarity between quantum mechanics and thermodynamics: Entropy, temperature, and Carnot cycle

Similarity between quantum mechanics and thermodynamics: Entropy, temperature, and Carnot cycle Similarity betwee quatum mechaics ad thermodyamics: Etropy, temperature, ad Carot cycle Sumiyoshi Abe 1,,3 ad Shiji Okuyama 1 1 Departmet of Physical Egieerig, Mie Uiversity, Mie 514-8507, Japa Istitut

More information

11 Correlation and Regression

11 Correlation and Regression 11 Correlatio Regressio 11.1 Multivariate Data Ofte we look at data where several variables are recorded for the same idividuals or samplig uits. For example, at a coastal weather statio, we might record

More information

Plastic deformation in flat-section band clamps

Plastic deformation in flat-section band clamps 93 Plastic deformatio i flat-sectio bad clamps K Shoghi, S M Barras, ad H V Rao Departmet of Mechaical Egieerig, Uiversity of Huddersfield, Huddersfield, UK The mauscript was received o 20 May 2004 ad

More information

Code_Aster. Viscoplastic constitutive law VISC_DRUC_PRAG

Code_Aster. Viscoplastic constitutive law VISC_DRUC_PRAG Titre : Loi de comportemet viscoplastique VISC_DRUC_PRAG Date : 29/05/2013 Page : 1/17 Viscoplastic costitutive law VISC_DRUC_PRAG Summarized: This documet the mod describes viscoplastic costitutive law

More information

MODELLING LOW VELOCITY IMPACT IN C/EPOXY LAMINATES

MODELLING LOW VELOCITY IMPACT IN C/EPOXY LAMINATES 3 RD INTERNATIONAL SYMPOSIUM ON BALLISTICS TARRAGONA, SPAIN 6- APRIL 7 MODELLING LOW VELOCITY IMPACT IN C/EPOXY LAMINATES E. Arévalo, C. Gozález,, F. Gálvez y J. LLorca Departmet of Material Sciece. Polytechic

More information

MONITORING THE STABILITY OF SLOPES BY GPS

MONITORING THE STABILITY OF SLOPES BY GPS MONITORING THE STABILITY OF SLOPES BY GPS Prof. S. Sakurai Costructio Egieerig Research Istitute Foudatio, Japa Prof. N. Shimizu Dept. of Civil Egieerig, Yamaguchi Uiversity, Japa ABSTRACT The stability

More information

The z-transform. 7.1 Introduction. 7.2 The z-transform Derivation of the z-transform: x[n] = z n LTI system, h[n] z = re j

The z-transform. 7.1 Introduction. 7.2 The z-transform Derivation of the z-transform: x[n] = z n LTI system, h[n] z = re j The -Trasform 7. Itroductio Geeralie the complex siusoidal represetatio offered by DTFT to a represetatio of complex expoetial sigals. Obtai more geeral characteristics for discrete-time LTI systems. 7.

More information

Elastic Plastic Behavior of Geomaterials: Modeling and Simulation Issues

Elastic Plastic Behavior of Geomaterials: Modeling and Simulation Issues Elastic Plastic Behavior of Geomaterials: Modelig ad Simulatio Issues Boris Zhaohui Yag (UA), Zhao Cheg (EarthMechaics Ic.), Mahdi Taiebat (UBC) Departmet of Civil ad Evirometal Egieerig Uiversity of Califoria,

More information

Probability, Expectation Value and Uncertainty

Probability, Expectation Value and Uncertainty Chapter 1 Probability, Expectatio Value ad Ucertaity We have see that the physically observable properties of a quatum system are represeted by Hermitea operators (also referred to as observables ) such

More information

Lecture 8: Solving the Heat, Laplace and Wave equations using finite difference methods

Lecture 8: Solving the Heat, Laplace and Wave equations using finite difference methods Itroductory lecture otes o Partial Differetial Equatios - c Athoy Peirce. Not to be copied, used, or revised without explicit writte permissio from the copyright ower. 1 Lecture 8: Solvig the Heat, Laplace

More information

CUMULATIVE DAMAGE ESTIMATION USING WAVELET TRANSFORM OF STRUCTURAL RESPONSE

CUMULATIVE DAMAGE ESTIMATION USING WAVELET TRANSFORM OF STRUCTURAL RESPONSE CUMULATIVE DAMAGE ESTIMATION USING WAVELET TRANSFORM OF STRUCTURAL RESPONSE Ryutaro SEGAWA 1, Shizuo YAMAMOTO, Akira SONE 3 Ad Arata MASUDA 4 SUMMARY Durig a strog earthquake, the respose of a structure

More information

LS-DYNA Simulation of Hemispherical-punch Stamping Process Using an Efficient Algorithm for Continuum Damage Based Elastoplastic Constitutive Equation

LS-DYNA Simulation of Hemispherical-punch Stamping Process Using an Efficient Algorithm for Continuum Damage Based Elastoplastic Constitutive Equation LS-YNA Simulatio of Hemispherical-puch Stampig Process Usig a Efficiet Algorithm for Cotiuum amage Based Elastoplastic Costitutive Equatio Nima Salajegheh*, Nader Abedrabbo ad Farhag Pourboghrat epartmet

More information

REGRESSION (Physics 1210 Notes, Partial Modified Appendix A)

REGRESSION (Physics 1210 Notes, Partial Modified Appendix A) REGRESSION (Physics 0 Notes, Partial Modified Appedix A) HOW TO PERFORM A LINEAR REGRESSION Cosider the followig data poits ad their graph (Table I ad Figure ): X Y 0 3 5 3 7 4 9 5 Table : Example Data

More information

ANALYSIS OF DAMPING EFFECT ON BEAM VIBRATION

ANALYSIS OF DAMPING EFFECT ON BEAM VIBRATION Molecular ad Quatum Acoustics vol. 7, (6) 79 ANALYSIS OF DAMPING EFFECT ON BEAM VIBRATION Jerzy FILIPIAK 1, Lech SOLARZ, Korad ZUBKO 1 Istitute of Electroic ad Cotrol Systems, Techical Uiversity of Czestochowa,

More information

577. Estimation of surface roughness using high frequency vibrations

577. Estimation of surface roughness using high frequency vibrations 577. Estimatio of surface roughess usig high frequecy vibratios V. Augutis, M. Sauoris, Kauas Uiversity of Techology Electroics ad Measuremets Systems Departmet Studetu str. 5-443, LT-5368 Kauas, Lithuaia

More information

THE LEVEL SET METHOD APPLIED TO THREE-DIMENSIONAL DETONATION WAVE PROPAGATION. Wen Shanggang, Sun Chengwei, Zhao Feng, Chen Jun

THE LEVEL SET METHOD APPLIED TO THREE-DIMENSIONAL DETONATION WAVE PROPAGATION. Wen Shanggang, Sun Chengwei, Zhao Feng, Chen Jun THE LEVEL SET METHOD APPLIED TO THREE-DIMENSIONAL DETONATION WAVE PROPAGATION We Shaggag, Su Chegwei, Zhao Feg, Che Ju Laboratory for Shock Wave ad Detoatio Physics Research, Southwest Istitute of Fluid

More information

Failure Theories Des Mach Elem Mech. Eng. Department Chulalongkorn University

Failure Theories Des Mach Elem Mech. Eng. Department Chulalongkorn University Failure Theories Review stress trasformatio Failure theories for ductile materials Maimum-Shear-Stress Theor Distortio-Eerg Theor Coulomb-Mohr Theor Failure theories for brittle materials Maimum-Normal-Stress

More information

Quantum Annealing for Heisenberg Spin Chains

Quantum Annealing for Heisenberg Spin Chains LA-UR # - Quatum Aealig for Heiseberg Spi Chais G.P. Berma, V.N. Gorshkov,, ad V.I.Tsifriovich Theoretical Divisio, Los Alamos Natioal Laboratory, Los Alamos, NM Istitute of Physics, Natioal Academy of

More information

Castiel, Supernatural, Season 6, Episode 18

Castiel, Supernatural, Season 6, Episode 18 13 Differetial Equatios the aswer to your questio ca best be epressed as a series of partial differetial equatios... Castiel, Superatural, Seaso 6, Episode 18 A differetial equatio is a mathematical equatio

More information

Approximating the ruin probability of finite-time surplus process with Adaptive Moving Total Exponential Least Square

Approximating the ruin probability of finite-time surplus process with Adaptive Moving Total Exponential Least Square WSEAS TRANSACTONS o BUSNESS ad ECONOMCS S. Khotama, S. Boothiem, W. Klogdee Approimatig the rui probability of fiite-time surplus process with Adaptive Movig Total Epoetial Least Square S. KHOTAMA, S.

More information

Finite Element Analysis of Rubber Bumper Used in Air-springs

Finite Element Analysis of Rubber Bumper Used in Air-springs Available olie at www.sciecedirect.com Procedia Egieerig 48 (0 ) 388 395 MMaMS 0 Fiite Elemet Aalysis of Rubber Bumper Used i Air-sprigs amas Makovits a *, amas Szabó b a Departmet of Mechaical Egieerig,

More information

Boundary layer problem on conveyor belt. Gabriella Bognár University of Miskolc 3515 Miskolc-Egyetemváros, Hungary

Boundary layer problem on conveyor belt. Gabriella Bognár University of Miskolc 3515 Miskolc-Egyetemváros, Hungary Boudary layer problem o coveyor belt Gabriella Bogár Uiversity of Miskolc 355 Miskolc-Egyetemváros, Hugary e-mail: matvbg@ui-miskolc.hu Abstract: A techologically importat source of the boudary layer pheomeo

More information

A new formulation of internal forces for non-linear hypoelastic constitutive models verifying conservation laws

A new formulation of internal forces for non-linear hypoelastic constitutive models verifying conservation laws A ew formulatio of iteral forces for o-liear hypoelastic costitutive models verifyig coservatio laws Ludovic NOELS 1, Lauret STAINIER 1, Jea-Philippe PONTHOT 1 ad Jérôme BONINI. 1 LTAS-Milieux Cotius et

More information

Boundary Element Method (BEM)

Boundary Element Method (BEM) Boudary Elemet Method BEM Zora Ilievski Wedesday 8 th Jue 006 HG 6.96 TU/e Talk Overview The idea of BEM ad its advatages The D potetial problem Numerical implemetatio Idea of BEM 3 Idea of BEM 4 Advatages

More information

The axial dispersion model for tubular reactors at steady state can be described by the following equations: dc dz R n cn = 0 (1) (2) 1 d 2 c.

The axial dispersion model for tubular reactors at steady state can be described by the following equations: dc dz R n cn = 0 (1) (2) 1 d 2 c. 5.4 Applicatio of Perturbatio Methods to the Dispersio Model for Tubular Reactors The axial dispersio model for tubular reactors at steady state ca be described by the followig equatios: d c Pe dz z =

More information

Lecture III-2: Light propagation in nonmagnetic

Lecture III-2: Light propagation in nonmagnetic A. La Rosa Lecture Notes ALIED OTIC Lecture III2: Light propagatio i omagetic materials 2.1 urface ( ), volume ( ), ad curret ( j ) desities produced by arizatio charges The objective i this sectio is

More information

The Mathematical Model and the Simulation Modelling Algoritm of the Multitiered Mechanical System

The Mathematical Model and the Simulation Modelling Algoritm of the Multitiered Mechanical System The Mathematical Model ad the Simulatio Modellig Algoritm of the Multitiered Mechaical System Demi Aatoliy, Kovalev Iva Dept. of Optical Digital Systems ad Techologies, The St. Petersburg Natioal Research

More information

MECHANICAL INTEGRITY DESIGN FOR FLARE SYSTEM ON FLOW INDUCED VIBRATION

MECHANICAL INTEGRITY DESIGN FOR FLARE SYSTEM ON FLOW INDUCED VIBRATION MECHANICAL INTEGRITY DESIGN FOR FLARE SYSTEM ON FLOW INDUCED VIBRATION Hisao Izuchi, Pricipal Egieerig Cosultat, Egieerig Solutio Uit, ChAS Project Operatios Masato Nishiguchi, Egieerig Solutio Uit, ChAS

More information

FAILURE CRITERIA: MOHR S CIRCLE AND PRINCIPAL STRESSES

FAILURE CRITERIA: MOHR S CIRCLE AND PRINCIPAL STRESSES LECTURE Third Editio FAILURE CRITERIA: MOHR S CIRCLE AND PRINCIPAL STRESSES A. J. Clark School of Egieerig Departmet of Civil ad Evirometal Egieerig Chapter 7.4 b Dr. Ibrahim A. Assakkaf SPRING 3 ENES

More information

Mathematical Modeling of Optimum 3 Step Stress Accelerated Life Testing for Generalized Pareto Distribution

Mathematical Modeling of Optimum 3 Step Stress Accelerated Life Testing for Generalized Pareto Distribution America Joural of Theoretical ad Applied Statistics 05; 4(: 6-69 Published olie May 8, 05 (http://www.sciecepublishiggroup.com/j/ajtas doi: 0.648/j.ajtas.05040. ISSN: 6-8999 (Prit; ISSN: 6-9006 (Olie Mathematical

More information

Performance of Stirling Engines * (Arranging Method of Experimental Results and Performance Prediction) Abstract. 1. Introduction

Performance of Stirling Engines * (Arranging Method of Experimental Results and Performance Prediction) Abstract. 1. Introduction 1 Performace of Stirlig Egies (Arragig Method of Experimetal Results ad Performace Predictio) Shoichi IWAMOTO Koichi HIRATA ad Fujio TODA Key Words: Stirlig Egie, Heat Exchager, Frictio, Egie Desig, Mechaical

More information

DETERMINATION OF MECHANICAL PROPERTIES OF A NON- UNIFORM BEAM USING THE MEASUREMENT OF THE EXCITED LONGITUDINAL ELASTIC VIBRATIONS.

DETERMINATION OF MECHANICAL PROPERTIES OF A NON- UNIFORM BEAM USING THE MEASUREMENT OF THE EXCITED LONGITUDINAL ELASTIC VIBRATIONS. ICSV4 Cairs Australia 9- July 7 DTRMINATION OF MCHANICAL PROPRTIS OF A NON- UNIFORM BAM USING TH MASURMNT OF TH XCITD LONGITUDINAL LASTIC VIBRATIONS Pavel Aokhi ad Vladimir Gordo Departmet of the mathematics

More information

METHOD OF FUNDAMENTAL SOLUTIONS FOR HELMHOLTZ EIGENVALUE PROBLEMS IN ELLIPTICAL DOMAINS

METHOD OF FUNDAMENTAL SOLUTIONS FOR HELMHOLTZ EIGENVALUE PROBLEMS IN ELLIPTICAL DOMAINS Please cite this article as: Staisław Kula, Method of fudametal solutios for Helmholtz eigevalue problems i elliptical domais, Scietific Research of the Istitute of Mathematics ad Computer Sciece, 009,

More information

FREE VIBRATION RESPONSE OF A SYSTEM WITH COULOMB DAMPING

FREE VIBRATION RESPONSE OF A SYSTEM WITH COULOMB DAMPING Mechaical Vibratios FREE VIBRATION RESPONSE OF A SYSTEM WITH COULOMB DAMPING A commo dampig mechaism occurrig i machies is caused by slidig frictio or dry frictio ad is called Coulomb dampig. Coulomb dampig

More information

3. Z Transform. Recall that the Fourier transform (FT) of a DT signal xn [ ] is ( ) [ ] = In order for the FT to exist in the finite magnitude sense,

3. Z Transform. Recall that the Fourier transform (FT) of a DT signal xn [ ] is ( ) [ ] = In order for the FT to exist in the finite magnitude sense, 3. Z Trasform Referece: Etire Chapter 3 of text. Recall that the Fourier trasform (FT) of a DT sigal x [ ] is ω ( ) [ ] X e = j jω k = xe I order for the FT to exist i the fiite magitude sese, S = x [

More information

Holistic Approach to the Periodic System of Elements

Holistic Approach to the Periodic System of Elements Holistic Approach to the Periodic System of Elemets N.N.Truov * D.I.Medeleyev Istitute for Metrology Russia, St.Peterburg. 190005 Moskovsky pr. 19 (Dated: February 20, 2009) Abstract: For studyig the objectivity

More information

Lecture 2: Monte Carlo Simulation

Lecture 2: Monte Carlo Simulation STAT/Q SCI 43: Itroductio to Resamplig ethods Sprig 27 Istructor: Ye-Chi Che Lecture 2: ote Carlo Simulatio 2 ote Carlo Itegratio Assume we wat to evaluate the followig itegratio: e x3 dx What ca we do?

More information

1 Adiabatic and diabatic representations

1 Adiabatic and diabatic representations 1 Adiabatic ad diabatic represetatios 1.1 Bor-Oppeheimer approximatio The time-idepedet Schrödiger equatio for both electroic ad uclear degrees of freedom is Ĥ Ψ(r, R) = E Ψ(r, R), (1) where the full molecular

More information

Chapter 6 Sampling Distributions

Chapter 6 Sampling Distributions Chapter 6 Samplig Distributios 1 I most experimets, we have more tha oe measuremet for ay give variable, each measuremet beig associated with oe radomly selected a member of a populatio. Hece we eed to

More information

11TH INTERNATIONAL SYMPOSIUM ON PARTICLE IMAGE VELOCIMETRY - PIV15 Santa Barbara, California, Sept , 2015

11TH INTERNATIONAL SYMPOSIUM ON PARTICLE IMAGE VELOCIMETRY - PIV15 Santa Barbara, California, Sept , 2015 11TH INTERNATIONAL SYMPOSIUM ON PARTICLE IMAGE VELOCIMETRY - PIV15 Sata Barbara, Califoria, Sept. 14-16, 2015 ABSTRACT HELE-SHAW RHEOMETRY BY MEANS OF PARTICLE IMAGE VELOCIMETRY Sita Drost & Jerry Westerweel

More information

SHAFTS: STATICALLY INDETERMINATE SHAFTS

SHAFTS: STATICALLY INDETERMINATE SHAFTS LECTURE SHAFTS: STATICALLY INDETERMINATE SHAFTS Third Editio A.. Clark School of Egieerig Departmet of Civil ad Eviromet Egieerig 7 Chapter 3.6 by Dr. Ibrahim A. Assakkaf SPRING 2003 ENES 220 Mechaics

More information

Analysis of deep drawing process to predict the forming severity considering inverse finite element and extended strain-based forming limit diagram

Analysis of deep drawing process to predict the forming severity considering inverse finite element and extended strain-based forming limit diagram Aalysis of deep drawig process to predict the formig severity cosiderig iverse fiite elemet ad exteded strai-based formig limit diagram M. Bosta Shiri a, R. Hashemi b ad A. Assempour c,* a School of biomedical

More information

FREE VIBRATIONS OF SIMPLY SUPPORTED BEAMS USING FOURIER SERIES

FREE VIBRATIONS OF SIMPLY SUPPORTED BEAMS USING FOURIER SERIES Abdullah : FREE VIBRATIONS OF SIMPY SUPPORTED BEAMS FREE VIBRATIONS OF SIMPY SUPPORTED BEAMS USING FOURIER SERIES SAWA MUBARAK ABDUAH Assistat ecturer Uiversity of Mosul Abstract Fourier series will be

More information

CEE 522 Autumn Uncertainty Concepts for Geotechnical Engineering

CEE 522 Autumn Uncertainty Concepts for Geotechnical Engineering CEE 5 Autum 005 Ucertaity Cocepts for Geotechical Egieerig Basic Termiology Set A set is a collectio of (mutually exclusive) objects or evets. The sample space is the (collectively exhaustive) collectio

More information

DISPERSION OF THE SEISMIC RAYLEIGH WAVES IN ELASTIC CRUSTAL LAYER ON VISCOELASTIC MANTLE

DISPERSION OF THE SEISMIC RAYLEIGH WAVES IN ELASTIC CRUSTAL LAYER ON VISCOELASTIC MANTLE DISPERSION OF THE SEISMIC RAYLEIGH WAVES IN ELASTIC CRUSTAL LAYER ON VISCOELASTIC MANTLE S.D. Akbarov, ad M. Negi 3 Prof. Dr., Mechaical Eg. Departmet, Yildiz Techical Uiversity, Istabul Istitute of Mathematics

More information

For example suppose we divide the interval [0,2] into 5 equal subintervals of length

For example suppose we divide the interval [0,2] into 5 equal subintervals of length Math 1206 Calculus Sec 1: Estimatig with Fiite Sums Abbreviatios: wrt with respect to! for all! there exists! therefore Def defiitio Th m Theorem sol solutio! perpedicular iff or! if ad oly if pt poit

More information

Streamfunction-Vorticity Formulation

Streamfunction-Vorticity Formulation Streamfuctio-Vorticity Formulatio A. Salih Departmet of Aerospace Egieerig Idia Istitute of Space Sciece ad Techology, Thiruvaathapuram March 2013 The streamfuctio-vorticity formulatio was amog the first

More information

PHYC - 505: Statistical Mechanics Homework Assignment 4 Solutions

PHYC - 505: Statistical Mechanics Homework Assignment 4 Solutions PHYC - 55: Statistical Mechaics Homewor Assigmet 4 Solutios Due February 5, 14 1. Cosider a ifiite classical chai of idetical masses coupled by earest eighbor sprigs with idetical sprig costats. a Write

More information

DISTRIBUTION LAW Okunev I.V.

DISTRIBUTION LAW Okunev I.V. 1 DISTRIBUTION LAW Okuev I.V. Distributio law belogs to a umber of the most complicated theoretical laws of mathematics. But it is also a very importat practical law. Nothig ca help uderstad complicated

More information

2C09 Design for seismic and climate changes

2C09 Design for seismic and climate changes 2C09 Desig for seismic ad climate chages Lecture 02: Dyamic respose of sigle-degree-of-freedom systems I Daiel Grecea, Politehica Uiversity of Timisoara 10/03/2014 Europea Erasmus Mudus Master Course Sustaiable

More information

BALKANTRIB O5 5 th INTERNATIONAL CONFERENCE ON TRIBOLOGY JUNE Kragujevac, Serbia and Montenegro

BALKANTRIB O5 5 th INTERNATIONAL CONFERENCE ON TRIBOLOGY JUNE Kragujevac, Serbia and Montenegro BALKANTRIB O5 5 th INTERNATIONAL CONFERENCE ON TRIBOLOGY JUNE.15-18. 5 Kragujevac, Serbia ad Moteegro A NEW APPROACH FOR DETERMINING THE FRICTION LAW IN METAL FORMING Sergei Alexadrov, Istitute for Problems

More information

INFLUENCE OF VARIABLE HEAT TRANSFER COEFFICIENT OF FIREWORKS AND CRACKERS ON THERMAL EXPLOSION CRITICAL AMBIENT TEMPERATURE AND TTIME TO IGNITION

INFLUENCE OF VARIABLE HEAT TRANSFER COEFFICIENT OF FIREWORKS AND CRACKERS ON THERMAL EXPLOSION CRITICAL AMBIENT TEMPERATURE AND TTIME TO IGNITION Guo, Z., et al.: Ifluece of Variable Heat Trasfer Coefficiet of Fireworks... THERMAL SCIENCE: Year 2016, Vol. 20, No. 5, pp. 1671-1679 1671 INFLUENCE OF VARIABLE HEAT TRANSFER COEFFICIENT OF FIREWORKS

More information

Numerical simulation of two-phase Darcy-Forchheimer flow during CO 2 injection into deep saline aquifers. Andi Zhang Feb. 4, 2013

Numerical simulation of two-phase Darcy-Forchheimer flow during CO 2 injection into deep saline aquifers. Andi Zhang Feb. 4, 2013 Numerical simulatio of two-phase Darcy-Forchheimer flow durig CO 2 ijectio ito deep salie aquifers Adi Zhag Feb. 4, 2013 Darcy flow VS o-darcy flow Darcy flow A liear relatioship betwee volumetric flow

More information

Applied Mathematical Modelling

Applied Mathematical Modelling Applied Mathematical Modellig 37 (013) 1177 1186 Cotets lists available at SciVerse ScieceDirect Applied Mathematical Modellig joural homepage: www.elsevier.com/locate/apm The effective permeability of

More information

Inelastic spherical caps with defects

Inelastic spherical caps with defects Ielastic spherical caps with defects JAAN LLLP RNST TUNGL Istitute of athematics Uiversity of Tartu Liivi str Tartu 549 STONIA aalellep@utee Tartu College Talli Uiversity of Techology Puiestee 78 58 Tartu

More information

THE NUMERICAL SOLUTION OF THE NEWTONIAN FLUIDS FLOW DUE TO A STRETCHING CYLINDER BY SOR ITERATIVE PROCEDURE ABSTRACT

THE NUMERICAL SOLUTION OF THE NEWTONIAN FLUIDS FLOW DUE TO A STRETCHING CYLINDER BY SOR ITERATIVE PROCEDURE ABSTRACT Europea Joural of Egieerig ad Techology Vol. 3 No., 5 ISSN 56-586 THE NUMERICAL SOLUTION OF THE NEWTONIAN FLUIDS FLOW DUE TO A STRETCHING CYLINDER BY SOR ITERATIVE PROCEDURE Atif Nazir, Tahir Mahmood ad

More information

Size, shape and temperature effect on nanomaterials

Size, shape and temperature effect on nanomaterials Idia Joural of Pure & Applied Physics Vol. 53, November 2015, pp. 768-775 Size, shape ad temperature effect o aomaterials G Sharma, S Bhatt, R Kumar & M Kumar* Departmet of Physics, G.B. Pat Uiversity

More information

MINIMUM-WEIGHT DESIGN OF COMPRESSIVELY LOADED COMPOSITE PLATES AND STIFFENED PANELS FOR POSTBUCKLING STRENGTH BY GENETIC ALGORITHM

MINIMUM-WEIGHT DESIGN OF COMPRESSIVELY LOADED COMPOSITE PLATES AND STIFFENED PANELS FOR POSTBUCKLING STRENGTH BY GENETIC ALGORITHM MINIMUM-WEIGHT DESIGN OF COMPRESSIVELY LOADED COMPOSITE PLATES AND STIFFENED PANELS FOR POSTBUCKLING STRENGTH BY GENETIC ALGORITHM Ji-Ho Kag a, Jae-Hu Lee a, Joo-Hyu Ha a, Chu-Go Kim a, ad Dog-Mi Lee b

More information

Microscopic traffic flow modeling

Microscopic traffic flow modeling Chapter 34 Microscopic traffic flow modelig 34.1 Overview Macroscopic modelig looks at traffic flow from a global perspective, whereas microscopic modelig, as the term suggests, gives attetio to the details

More information

8. СОВЕТУВАЊЕ. Охрид, септември ANALYSIS OF NO LOAD APPARENT POWER AND FREQUENCY SPECTRUM OF MAGNETIZING CURRENT FOR DIFFERENT CORE TYPES

8. СОВЕТУВАЊЕ. Охрид, септември ANALYSIS OF NO LOAD APPARENT POWER AND FREQUENCY SPECTRUM OF MAGNETIZING CURRENT FOR DIFFERENT CORE TYPES 8. СОВЕТУВАЊЕ Охрид, 22 24 септември Leoardo Štrac Frajo Keleme Kočar Power Trasformers Ltd. ANALYSS OF NO LOAD APPARENT POWER AND FREQENCY SPECTRM OF MAGNETZNG CRRENT FOR DFFERENT CORE TYPES ABSTRACT

More information

Direction of Arrival Estimation Method in Underdetermined Condition Zhang Youzhi a, Li Weibo b, Wang Hanli c

Direction of Arrival Estimation Method in Underdetermined Condition Zhang Youzhi a, Li Weibo b, Wang Hanli c 4th Iteratioal Coferece o Advaced Materials ad Iformatio Techology Processig (AMITP 06) Directio of Arrival Estimatio Method i Uderdetermied Coditio Zhag Youzhi a, Li eibo b, ag Hali c Naval Aeroautical

More information

A proposed discrete distribution for the statistical modeling of

A proposed discrete distribution for the statistical modeling of It. Statistical Ist.: Proc. 58th World Statistical Cogress, 0, Dubli (Sessio CPS047) p.5059 A proposed discrete distributio for the statistical modelig of Likert data Kidd, Marti Cetre for Statistical

More information

Using An Accelerating Method With The Trapezoidal And Mid-Point Rules To Evaluate The Double Integrals With Continuous Integrands Numerically

Using An Accelerating Method With The Trapezoidal And Mid-Point Rules To Evaluate The Double Integrals With Continuous Integrands Numerically ISSN -50 (Paper) ISSN 5-05 (Olie) Vol.7, No., 017 Usig A Acceleratig Method With The Trapezoidal Ad Mid-Poit Rules To Evaluate The Double Itegrals With Cotiuous Itegrads Numerically Azal Taha Abdul Wahab

More information

GUIDELINES ON REPRESENTATIVE SAMPLING

GUIDELINES ON REPRESENTATIVE SAMPLING DRUGS WORKING GROUP VALIDATION OF THE GUIDELINES ON REPRESENTATIVE SAMPLING DOCUMENT TYPE : REF. CODE: ISSUE NO: ISSUE DATE: VALIDATION REPORT DWG-SGL-001 002 08 DECEMBER 2012 Ref code: DWG-SGL-001 Issue

More information

PILOT STUDY ON THE HORIZONTAL SHEAR BEHAVIOUR OF FRP RUBBER ISOLATORS

PILOT STUDY ON THE HORIZONTAL SHEAR BEHAVIOUR OF FRP RUBBER ISOLATORS Asia-Pacific Coferece o FRP i Structures (APFIS 2007) S.T. Smith (ed) 2007 Iteratioal Istitute for FRP i Costructio PILOT STUDY ON THE HORIZONTAL SHEAR BEHAVIOUR OF FRP RUBBER ISOLATORS T.B. Peg *, J.Z.

More information

Obtaining Constants of Johnson-Cook Material Model Using a Combined Experimental, Numerical Simulation and Optimization Method

Obtaining Constants of Johnson-Cook Material Model Using a Combined Experimental, Numerical Simulation and Optimization Method Vol:0, No:9, 06 Obtaiig Costats of Johso-Cook Material Model Usig a Combied Experimetal, Numerical Simulatio ad Optimizatio Method F. Rahimi Dehgola, M. Behzadi, J. Fathi Sola Iteratioal Sciece Idex, Mechaical

More information

Course 4: Preparation for Calculus Unit 1: Families of Functions

Course 4: Preparation for Calculus Unit 1: Families of Functions Course 4: Preparatio for Calculus Uit 1: Families of Fuctios Review ad exted properties of basic fuctio families ad their uses i mathematical modelig Develop strategies for fidig rules of fuctios whose

More information

Comparison Study of Series Approximation. and Convergence between Chebyshev. and Legendre Series

Comparison Study of Series Approximation. and Convergence between Chebyshev. and Legendre Series Applied Mathematical Scieces, Vol. 7, 03, o. 6, 3-337 HIKARI Ltd, www.m-hikari.com http://d.doi.org/0.988/ams.03.3430 Compariso Study of Series Approimatio ad Covergece betwee Chebyshev ad Legedre Series

More information

Finally, we show how to determine the moments of an impulse response based on the example of the dispersion model.

Finally, we show how to determine the moments of an impulse response based on the example of the dispersion model. 5.3 Determiatio of Momets Fially, we show how to determie the momets of a impulse respose based o the example of the dispersio model. For the dispersio model we have that E θ (θ ) curve is give by eq (4).

More information

Numerical Study on MHD Flow And Heat Transfer With The Effect Of Microrotational Parameter In The Porous Medium

Numerical Study on MHD Flow And Heat Transfer With The Effect Of Microrotational Parameter In The Porous Medium IOSR Joural of Egieerig (IOSRJEN) ISSN (e): 5-3, ISSN (p): 78-879 Vol. 5, Issue 4 (April. 5), V PP 8-7 www.iosrje.org Numerical Study o MHD Flow Ad Heat rasfer With he Effect Of Microrotatioal Parameter

More information

Lecture 9: Diffusion, Electrostatics review, and Capacitors. Context

Lecture 9: Diffusion, Electrostatics review, and Capacitors. Context EECS 5 Sprig 4, Lecture 9 Lecture 9: Diffusio, Electrostatics review, ad Capacitors EECS 5 Sprig 4, Lecture 9 Cotext I the last lecture, we looked at the carriers i a eutral semicoductor, ad drift currets

More information

GROUND MOTION OF NON-CIRCULAR ALLUVIAL VALLEY FOR INCIDENT PLANE SH-WAVE. Hui QI, Yong SHI, Jingfu NAN

GROUND MOTION OF NON-CIRCULAR ALLUVIAL VALLEY FOR INCIDENT PLANE SH-WAVE. Hui QI, Yong SHI, Jingfu NAN The th World Coferece o Earthquake Egieerig October -7, 8, Beiig, Chia GROUND MOTION OF NON-CIRCULAR ALLUVIAL VALLEY FOR INCIDENT PLANE SH-WAVE Hui QI, Yog SHI, Jigfu NAN ABSTRACT : Professor, Dept. of

More information

PSYCHOLOGICAL RESEARCH (PYC 304-C) Lecture 9

PSYCHOLOGICAL RESEARCH (PYC 304-C) Lecture 9 Hypothesis testig PSYCHOLOGICAL RESEARCH (PYC 34-C Lecture 9 Statistical iferece is that brach of Statistics i which oe typically makes a statemet about a populatio based upo the results of a sample. I

More information

SHAFTS: STATICALLY INDETERMINATE SHAFTS

SHAFTS: STATICALLY INDETERMINATE SHAFTS SHAFTS: STATICALLY INDETERMINATE SHAFTS Up to this poit, the resses i a shaft has bee limited to shearig resses. This due to the fact that the selectio of the elemet uder udy was orieted i such a way that

More information

Chapter 1. Complex Numbers. Dr. Pulak Sahoo

Chapter 1. Complex Numbers. Dr. Pulak Sahoo Chapter 1 Complex Numbers BY Dr. Pulak Sahoo Assistat Professor Departmet of Mathematics Uiversity Of Kalyai West Begal, Idia E-mail : sahoopulak1@gmail.com 1 Module-2: Stereographic Projectio 1 Euler

More information

Finite Difference Derivations for Spreadsheet Modeling John C. Walton Modified: November 15, 2007 jcw

Finite Difference Derivations for Spreadsheet Modeling John C. Walton Modified: November 15, 2007 jcw Fiite Differece Derivatios for Spreadsheet Modelig Joh C. Walto Modified: November 15, 2007 jcw Figure 1. Suset with 11 swas o Little Platte Lake, Michiga. Page 1 Modificatio Date: November 15, 2007 Review

More information

Let us give one more example of MLE. Example 3. The uniform distribution U[0, θ] on the interval [0, θ] has p.d.f.

Let us give one more example of MLE. Example 3. The uniform distribution U[0, θ] on the interval [0, θ] has p.d.f. Lecture 5 Let us give oe more example of MLE. Example 3. The uiform distributio U[0, ] o the iterval [0, ] has p.d.f. { 1 f(x =, 0 x, 0, otherwise The likelihood fuctio ϕ( = f(x i = 1 I(X 1,..., X [0,

More information

Added mass estimation of open flat membranes vibrating in still air

Added mass estimation of open flat membranes vibrating in still air Added mass estimatio of ope flat membraes vibratig i still air *Yua-qi Li 1),Yi Zhou ), Akihito YOHIDA 3) ad Yukio TAMURA 4) 1), ) Departmet of Buildig Egieerig, Togji Uiversity, haghai, 9, Chia 1) liyq@togji.edu.c

More information

Supplementary Information: Flexible and tunable silicon photonic circuits on plastic substrates

Supplementary Information: Flexible and tunable silicon photonic circuits on plastic substrates Supplemetary Iformatio: Flexible ad tuable silico photoic circuits o plastic substrates Yu Che 1, Hua Li 1, ad Mo Li 11 1 Departmet of Electrical ad Computer Egieerig, Uiversity of Miesota, Mieapolis,

More information

17 Phonons and conduction electrons in solids (Hiroshi Matsuoka)

17 Phonons and conduction electrons in solids (Hiroshi Matsuoka) 7 Phoos ad coductio electros i solids Hiroshi Matsuoa I this chapter we will discuss a miimal microscopic model for phoos i a solid ad a miimal microscopic model for coductio electros i a simple metal.

More information

MATHEMATICS. The assessment objectives of the Compulsory Part are to test the candidates :

MATHEMATICS. The assessment objectives of the Compulsory Part are to test the candidates : MATHEMATICS INTRODUCTION The public assessmet of this subject is based o the Curriculum ad Assessmet Guide (Secodary 4 6) Mathematics joitly prepared by the Curriculum Developmet Coucil ad the Hog Kog

More information

ANALYSIS OF EXPERIMENTAL ERRORS

ANALYSIS OF EXPERIMENTAL ERRORS ANALYSIS OF EXPERIMENTAL ERRORS All physical measuremets ecoutered i the verificatio of physics theories ad cocepts are subject to ucertaities that deped o the measurig istrumets used ad the coditios uder

More information

If, for instance, we were required to test whether the population mean μ could be equal to a certain value μ

If, for instance, we were required to test whether the population mean μ could be equal to a certain value μ STATISTICAL INFERENCE INTRODUCTION Statistical iferece is that brach of Statistics i which oe typically makes a statemet about a populatio based upo the results of a sample. I oesample testig, we essetially

More information

True Nature of Potential Energy of a Hydrogen Atom

True Nature of Potential Energy of a Hydrogen Atom True Nature of Potetial Eergy of a Hydroge Atom Koshu Suto Key words: Bohr Radius, Potetial Eergy, Rest Mass Eergy, Classical Electro Radius PACS codes: 365Sq, 365-w, 33+p Abstract I cosiderig the potetial

More information

1. Hydrogen Atom: 3p State

1. Hydrogen Atom: 3p State 7633A QUANTUM MECHANICS I - solutio set - autum. Hydroge Atom: 3p State Let us assume that a hydroge atom is i a 3p state. Show that the radial part of its wave fuctio is r u 3(r) = 4 8 6 e r 3 r(6 r).

More information