Crack onset assessment near the sharp material inclusion tip by means of modified maximum tangential stress criterion
|
|
- Shanna Fowler
- 5 years ago
- Views:
Transcription
1 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 of Physics of Materials AS CR, Žižkova 22, Brno , Czech Republic ABSTRACT. In the case of particle reinforce composites, where the particles are in a form of sharp material inclusions, singular stress concentration exists on each tip of each inclusion. This is ue to the geometric an material iscontinuities between matrix an particle. These points of stress concentration are susceptible of crack initiation an thus often responsible for failure of the whole structure. The moifie maximum tangential stress criterion is employe in orer to preict crack onset conitions. KEYWORDS. Sharp Material Inclusion; Singular an Non-singular Stress Terms; Generalize Fracture Mechanics. Citation: Krepl, O., Klusák, J., Crack onset assessment near the sharp material inclusion tip by means of moifie maximum tangential stress criterion, Frattura e Integrità Strutturale, 42 (217) Receive: Accepte: Publishe: Copyright: 217 This is an open access article uner the terms of the CC-BY 4., which permits unrestricte use, istribution, an reprouction in any meium, provie the original author an source are creite. INTRODUCTION C omposites fin application in a variety of engineering structures. Overall excellent properties, which can only be achieve by combination of two or more homogenous materials, bring solution for high emans given by contemporary avance technologies. owever, the very nature of composites also causes a ifficulty in their assessment in terms of fracture mechanics. Common composite types, i.e. fibre reinforce composites or particle reinforce composites give rise to points of singular stress concentration. These points are foun in loci where fibre ens are embee in matrix or at each tip of each sharp particle an are referre as General Singular Stress Concentrators (GSSCs). This paper eals with the case of particle embee in matrix which is in a form of Sharp Material Inclusion (SMI). The strength of singularity in this case is ifferent an lower than in the case of a crack. owever, the GSSCs are points often susceptible for crack initiation an thus they can be responsible for failure of the whole structure. A thorough unerstaning of crack initiation conitions in GSSC leas to improve estimation of the critical failure loa. Stress, energy or couple criteria are commonly use in fracture mechanics analyses of many GSSC types such as sharp notches. These criteria can be employe in orer to preict crack onset conitions to the case of SMI. The GSSC characterize by weak singularity often require consieration of higher orer non-singular terms of the asymptotic stress series to escribe the stress state an thus to receive consistent critical loa preictions. Suggestion of moifie maximum tangential stress stability criterion for the SMI with consieration of non-singular terms is the main objective of this paper. 66
2 FRACTURE MECANICS OF SARP MATERIAL INCLUSION T he SMI is moelle in 2D as a special case of multi-material junction, the bi-material junction [1, 2] as shown in Fig. 1. Both material regions are consiere to consist of linear elastic material an fully escribe by Young s mouli E 1 an E 2 an Poisson s ratios ν 1 an ν 2. The geometry is characterize by the opening angle α. Perfect boning (traction an isplacement continuity) is assume at both interfaces Γ an Γ 1. The problem is consiere either in a state of plane stress or plane strain. The stress fiel in the vicinity of bi-material junction tip, i.e. when r, is escribe by following asymptotic series: (1) ij r fij 1m r fij 2m r fij 3m where k is the kth Generalize Stress Intensity Factor (GSIF) which correspons to the kth eigenvalue k, that forms the stress singularity exponent (1 k ). The terms of the series can be either singular, when < ( k) < 1 or non-singular when 1 < ( k ). As r the singular terms become unboune, while the non-singular terms vanish. f ijkm (θ) is the imensionless angular eigenfunction constructe for ijth component of the stress tensor, kth eigenvalue an mth material as in [2, 1]. Note that the series above is in sake of simplicity written for real eigenvalues an real GSIFs. Such form of the series provies satisfactory escription for most of the cases. Figure 1: Sharp material inclusion moel. In majority of the fracture mechanics analyses of GSSCs only the singular terms are use for escription of stress fiel [3, 4] an following etermination of crack onset conitions. In [5, 6] the effect of first non-singular stress term on stress istribution in case of sharp V-notch is stuie. The effect of the first non-singular term in case of bi-material notches is stuie in [7, 8]. In [9] Klusák et al. have shown the significance of consieration of the first non-singular term in the case of sharp bi-material orthotropic plate. A stuy which has shown the effect of non-singular terms in the cases of SMI has been conucte in [1]. In orer to ientify singular an non-singular terms for given configuration the epenence of the eigenvalues 1, 2 an 3 of chosen geometric SMI cases on Young s mouli ratio E 1/E 2 is presente in Fig. 2. ere, material region 1 an Young s moulus E 1 belongs to inclusion an the region 2 with moulus E 2 correspons to the matrix. In the case of rectangular SMI characterize by = 9 two singular an one non-singular terms for all E 1 /E 2 ratios are foun. In the case of sharper SMI with = 6 there are two singular terms together with one non-singular terms for cases where inclusion is more compliant than matrix, i.e. E 1 /E 2 < 1 an one singular term together with two non-singular terms for cases of inclusion stiffer than matrix, i.e. E 1 /E 2 >1. The SMI with a blunt opening angle of = 12 shows reverse tren, as the cases of more compliant inclusion have one singular term together with two non-singular terms an cases with stiffer inclusion have two singular together with one non-singular term. 67
3 Figure 2: Depenence of 1, 2 an 3 on E 1 /E 2 shown for 3 SMI geometries = 6, 9, 12. To preict crack initiation conitions, a moifie maximum tangential stress criterion is use. The tangential stress epens on raial istance from the singular concentrator tip. To mitigate the raial istance epenence the mean value can be calculate as: 1 m r, r (2) The parameter is relate to the fracture mechanism, e.g. in case of cleavage fracture it can be set as = 2-5 grain size of the material [13]. The criterion states that the crack will initiate in the irection θ,m of maximum value of mean tangential stress. The extreme value is foun as: 2 ; 2, m, m (3) Let s consier the first 3 terms (singular an non-singular) of the stress series (1). By averaging it over specific istance as in (2) an substituting it into Eq. (3) we obtain: f f f m 2m 3m (4) The first GSIF 1 is factore out an the consequent ratios are enote as k1 = k/ 1. Equation has now the only unknown, the crack initiation irection θ,m. f f f m 2m 3m (5) 68
4 In general, global an local extremes can be foun in bi-material problem. The crack will be initiate in the irection of maximum of an corresponing minimum value of generalize fracture toughness [14]. Note that the crack initiation irection oes not epen on the absolute values of GSIFs but rather on their ratios k1. The general form of the equation for n terms of the stress series is: n k1 k k1 k f km (6) The crack initiation conitions of the SMI are escribe by means of critical stress quantity. In this way, Knésl in [13] assesses a stability of V-notch. e proposes comparison of such critical quantity of the crack an the V-notch. Since in both cases (crack tip an the V-notch tip) the stress is singular, he supposes that crack initiation mechanism in V-notch will be the same as the crack propagation mechanism. Then crack initiation conitions in case of GSSCs can be assesse by means of critical average stress ascertaine for a crack. For a crack in homogeneous material an loae in normal moe I, the critical mean tangential stress value is: IC C 2K (7) 2 Whereas for the GSIF with consieration of first three singular an non-singular terms of the series (1), it is: C, m C, m 11 f 1m 21 f 2m 31 f 3m (8) By comparison of the two relations above the generalize fracture toughness 1C,m is: 2 IC, m 1C, m K f ( ) f ( ) f ( ) 11 1m, m 21 2m, m 31 3m, m (9) which in general form for n terms of the stress series is written: 2K IC, m 1C, m 1 k 2 n k1k1 fkm, m k ( ) The generalize fracture toughness 1C,m epens on fracture toughness K IC,m of the material m. It is obvious that in case of SMI two materials m = 1 an 2 shoul be consiere an possibility of crack initiation into the matrix an into the inclusion shoul be evaluate. If the value 1C,1 is lower than 1C,2 crack initiation is expecte into the inclusion, otherwise it occurs into matrix. In sake of completeness, the case of crack initiation into the interface shoul be also evaluate. In that case, the value of 1C,interface has to be calculate base on fracture toughness of the interface K IC,interface. Note that for all the critical values 1C,1, 1C,2, an 1C,interface the shape functions f,m shoul contain corresponing angle of potential crack initiation,m (m = 1, 2, interface). The crack initiation occurs if the following stability criterion is violate (11): min,, 1 1C,1 1C,2 1C,interface Crack is not initiate in SMI tip if the value of GSIF 1 is lower than its critical value. Finally, the critical applie loa σ crit can be calculate as: (1) (11) 69
5 crit appl 1 appl min 1C,1,1, 1C,2,2, 1C,interface 1,2 ) (12) NUMERICAL EXAMPLE T he SMI is moelle by FEM in 2D. The moel geometry an bounary conitions are shown in Fig. 3. The rectangular inclusion with α = 9 in a plane strain conitions is stuie. The moelle specimen size is mm. The moel is loae with unit tension of σ appl = 1 MPa. In this numerical example the inclusion more compliant than matrix is consiere, which represents sanstone inclusion in cement paste of various stiffness. The moelle bimaterial configurations are liste in Tab. 1. The Poisson s ratios are ientical for both materials, ν 1 = ν 2 =.2. The eigenvalues λ k are etermine as a solution of an eigenvalue problem [2, 1]. The GSIFs k are calculate semi-analytically by means of overeterministic metho [2, 1], with resulting values as shown in Tab. 1. The iameter which provies inputs (noal isplacements) for GSIFs calculation is chosen as r 1 = 1 mm. The characteristic length which correspons to the fracture mechanism is for the case of cement paste composites chosen as = 1 mm. Fig. 4 shows the tangential stress istribution for the case E 1 /E 2 =.5, nevertheless istribution of this kin is characteristic to the cases of inclusion more compliant than matrix in general. The istribution isposes of two extremes, global which is foun in the matrix θ,2 = 18 an local foun in the inclusion θ,1 =. The values of crack initiation angles θ,m are given by the symmetry of the problem. For the same reason the o term, the skew relate GSIF 2 = MPa m 1-λ2. It is necessary to asses potential crack initiation in all possible ways, i.e. in the irection of the global maximum, local maximum an the interface. The crack initiation angle in the latter case is θ,interface = 45. In Tab. 1 the generalize fracture toughness is calculate for global function extremes, the local one an the interface. The fracture toughness value of matrix, inclusion an the interface is chosen equally, K IC,m = 1 MPa m 1/2 for m = {1, 2, interface}. Generalize fracture toughness for unit K IC,m can be unerstoo as normalize value an in Tab. 1 is enote as. * 1C,m E 1 [GPa] E 2 [GPa] E 1 /E 2 1 [MPa m 1-λ1 ] 3 [MPa m 1-λ3 ] K IC [MPa m 1/2 ] * 1C,2 [MPa m 1-λ1 ] * 1C,1 [MPa m 1-λ1 ] * 1C,interface [MPa m 1-λ1 ] Table 1: Bi-material configurations with resulting GSIFs an values of generalize fracture toughness. Figure 3: The numerical moel. 7
6 1, 3 It is evient that the stress escription (Fig. 4) by means of the first an the thir stress terms ( ) is more precise than 1 the escription by the first stress term only ( ). It can be seen especially in irections θ = 9 an 27 corresponing to regions of lower tangential stresses. The stability criterion employing critical value 1C following from all three stress terms (9), (1) is use. Fig. 5 shows results of the values 1 an the critical values * 1C corresponing to crack initiation into matrix, into the inclusion an into the interface. The values 1 follow from FEM numerical solution of moel with unit loaing 1 MPa. The critical values * 1C,m are ascertaine for unit fracture toughness K IC,m = 1 MPa m 1/2 for m = {1, 2, interface}. Particular values 1C,m for given fracture toughness of the matrix, the inclusion an the interface are: K (13) * 1C, m 1C,m IC, m Figure 4: Tangential stress istribution for the case E 1 /E 2 =.5. Yellow an magenta curve stay for the analytical tangential stress solution on r = 1 mm obtaine by one term 1 an by two terms 1 an 3 respectively. The cyan curve is the mean tangential stress obtaine by two terms 1 an 3. Black ots represent the FEM solution on r = 1 mm. The yellow vertical lines enote the interfaces an black vertical lines the function local an global extreme. The results show that crack initiation conitions epen on the ratio of materials Young's mouli, an on the fracture toughness of material components (matrix, inclusion) an the interface. In the numerical example E 1 correspons to Young's moulus of the inclusion. The inclusion here is consiere to be a sanstone aggregate with E 1 = 2 GPa. Corresponing fracture toughness of sanstone is between.28 an.52 MPa m 1/2. On the other han, E 2 correspons to Young's moulus of matrix. In the numerical example it varies from 26.6 to 8 GPa. The fracture toughness of common harene cement paste is between.1 an.8 MPa m 1/2, where higher values of E 2 usually match to higher values of K IC,2. The fracture toughness of interface can wiely vary. In case of silicate base composites it is epenent on the manufacturing process an evelopment of interfacial transition zone. The values 1C,interface are higher than 1C,1 an 1C,2 thus it seems that the crack kink to the interface is not probable, but min that usually K IC,interface is lower than K IC,matrix an K IC,inclusion. When consiering particular values of K IC,m for particular ratios E 1/E 2, the curves of the critical values 1C,m for matrix, the inclusion an the interface will change by multiplying as inicate in (13). Thus every particular case will show if crack initiation occurs to matrix, to the inclusion or the interface, i.e. min{ 1C,1, 1C,2, 1C,interface}. Let us note that stuie material moel is universal. Although it supposes sharp concave inclusion embee in matrix, at the same time it can escribe convex corner of stiffer inclusion fille with matrix. In this case E 1 woul match to matrix an E 2 to the inclusion. Principally, the approaches will be the same an the results very similar. 71
7 Figure 5: GSIF 1 for unit applie loa an the critical values 1C neee for crack initiation to matrix an to the inclusion. CONCLUSIONS T he article eals with generalize fracture mechanics of a crack initiation in a tip of polygonal material inclusion embee in matrix. Stability criterion for etermination of crack initiation conitions is propose by means of average value of tangential stress. The stability criterion is written for stress escribe by singular an non-singular (higher orer) terms. The criterion employing higher orer terms is necessary for the assessment of sharp material inclusion, because singular stress terms only o not escribe stress state satisfactorily in some material an geometrical configurations. The propose criterion allows inicating whether crack initiation occurs into matrix, into the inclusion or into the interface. Its general form allows such assessment for particular known fracture toughness of the matrix, inclusion or the interface. Description of stability of stress concentrators of this kin contribute to better unerstaning of toughening mechanisms in particle composites. This kin of unerstaning leas to composite esign optimisation or composite structure failure prevention. ACKNOWLEDGEMENTS T his research has been financially supporte by the Ministry of Eucation, Youth an Sports of the Czech Republic uner the project CEITEC 22 (LQ161). The authors woul like to thank the Czech Science Founation for financial support through the Grant 16/1872S. REFERENCES [1] Paggi, M., Carpinteri, A., On the Stress Singularities at Multimaterial Interfaces an Relate Analogies With Flui Dynamics an Diffusion, Applie Mechanics Reviews, 61 (28) [2] Krepl, O., Klusák, J., Reconstruction of a 2D stress fiel aroun the tip of a sharp material inclusion, Proceia Structural Integrity, 2 (216)
8 [3] Williams, M.L., Stress singularities resulting from various bounary conitions in angular corners of plates in extension. J Appl Mech,19 (1952) [4] Seweryn, A., Molski, K., Elastic stress singularities an corresponing generalize stress intensity factors for angular corners uner various bounary conitions, Engineering Fracture Mechanics, 55(4) (1996) [5] Kim, J.K., Cho, S.B., Effect of secon non-singular term of moe I near the tip of a V-notche crack, Fatigue an Fracture of Engineering Materials an Structures, 32 (29) [6] Ayatollahi, M.R., Dehghany, M., Nejati, M., Fracture analysis of V-notche components Effects of first non-singular stress term, International Journal of Solis an Structures, 48 (211) [7] Ayatollahi, M.R., Mirsayar, M.M., Nejati, M., Evaluation of first non-singular stress term in bi-material notches, Computational Materials Science, 5 (21) [8] Mirsayar, M.M., Aliha, M.R.M., Samaei, A.T., On fracture initiation angle near bi-material notches Effects of first non-singular stress term, Engineering Fracture Mechanics, 119 (214) [9] Klusák, J., rstka, M., Profant, T., Krepl, O., Ševeček, O., Kotoul M., The influence of the first non-singular stress terms on crack initiation irection in an orthotropic bi-material plate, Theoretical an Applie Fracture Mechanics, 71 (214) [1] Krepl, O., Klusák, J., The influence of non-singular terms on the precision of stress escription near a sharp material inclusion tip, Theoretical an Applie Fracture Mechanics, in press. DOI: 1.116/j.tafmec [11] Sih, G.C., A special theory of crack propagation, in: G.C. Sih (E.), Mechanics of Fracture Methos of analysis an solutions of crack problems, Noorhoff International Publishing, Leyen, The Netherlans, (1973) XXI-XLV. [12] Sih, G.C., Mechanics of Fracture Initiation an Propagation, Kluwer Acaemic Publishing, Boston, (1991). [13] Knésl, Z., A criterion of V-notch stability, Int J of Fracture, 48 (1991) [14] Klusák, J. Profant, T., Knésl, Z, Kotoul M., The influence of iscontinuity an orthotropy of fracture toughness on conitions of fracture initiation in singular stress concentrators, Engineering Fracture Mechanics, 11 (213)
Critical applied stresses for a crack initiation from a sharp V-notch
Focussed on: Fracture and Structural Integrity related Issues Critical applied stresses for a crack initiation from a sharp V-notch L. Náhlík, P. Hutař Institute of Physics of Materials, Academy of Sciences
More information3-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 informationCRACK INITIATION CRITERIA FOR SINGULAR STRESS CONCENTRATIONS Part I: A Universal Assessment of Singular Stress Concentrations
Engineering MECHANICS, Vol. 14, 2007, No. 6, p. 399 408 399 CRACK INITIATION CRITERIA FOR SINGULAR STRESS CONCENTRATIONS Part I: A Universal Assessment of Singular Stress Concentrations Zdeněk Knésl, Jan
More informationCrack-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'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 informationStrength 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 informationThermal 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 information05 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 informationOF 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 informationChapter 4. Electrostatics of Macroscopic Media
Chapter 4. Electrostatics of Macroscopic Meia 4.1 Multipole Expansion Approximate potentials at large istances 3 x' x' (x') x x' x x Fig 4.1 We consier the potential in the far-fiel region (see Fig. 4.1
More informationP. 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 informationThe thin plate theory assumes the following deformational kinematics:
MEG6007 (Fall, 2017) - Solutions of Homework # 8 (ue on Tuesay, 29 November 2017) 8.0 Variational Derivation of Thin Plate Theory The thin plate theory assumes the following eformational kinematics: u
More informationAssessment 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 informationSimulation of Angle Beam Ultrasonic Testing with a Personal Computer
Key Engineering Materials Online: 4-8-5 I: 66-9795, Vols. 7-73, pp 38-33 oi:.48/www.scientific.net/kem.7-73.38 4 rans ech ublications, witzerlan Citation & Copyright (to be inserte by the publisher imulation
More informationThe 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 informationEngineering Solid Mechanics
Engineering Solid Mechanics (3) 49-53 Contents lists available at GrowingScience Engineering Solid Mechanics homepage: www.growingscience.com/esm Calculation of stress intensity factors for an interfacial
More informationA simple model for the small-strain behaviour of soils
A simple moel for the small-strain behaviour of soils José Jorge Naer Department of Structural an Geotechnical ngineering, Polytechnic School, University of São Paulo 05508-900, São Paulo, Brazil, e-mail:
More informationCode_Aster. Detection of the singularities and computation of a card of size of elements
Titre : Détection es singularités et calcul une carte [...] Date : 0/0/0 Page : /6 Responsable : Josselin DLMAS Clé : R4.0.04 Révision : 9755 Detection of the singularities an computation of a car of size
More informationTOWARDS 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 informationVectors in two dimensions
Vectors in two imensions Until now, we have been working in one imension only The main reason for this is to become familiar with the main physical ieas like Newton s secon law, without the aitional complication
More informationCharacterization of lead zirconate titanate piezoceramic using high frequency ultrasonic spectroscopy
JOURNAL OF APPLIED PHYSICS VOLUME 85, NUMBER 1 15 JUNE 1999 Characterization of lea zirconate titanate piezoceramic using high frequency ultrasonic spectroscopy Haifeng Wang, Wenhua Jiang, a) an Wenwu
More informationVERIFICATION OF BRITTLE FRACTURE CRITERIA FOR BIMATERIAL STRUCTURES
VERIFICATION OF BRITTLE FRACTURE CRITERIA FOR BIMATERIAL STRUCTURES Grzegorz MIECZKOWSKI *, Krzysztof MOLSKI * * Faculty of Mechanical Engineering, Białystok University of Technology, ul. Wiejska 45C,
More informationChapter 2 Governing Equations
Chapter 2 Governing Equations In the present an the subsequent chapters, we shall, either irectly or inirectly, be concerne with the bounary-layer flow of an incompressible viscous flui without any involvement
More informationThe Three-dimensional Schödinger Equation
The Three-imensional Schöinger Equation R. L. Herman November 7, 016 Schröinger Equation in Spherical Coorinates We seek to solve the Schröinger equation with spherical symmetry using the metho of separation
More informationInternational 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 informationCode_Aster. Detection of the singularities and calculation of a map of size of elements
Titre : Détection es singularités et calcul une carte [...] Date : 0/0/0 Page : /6 Responsable : DLMAS Josselin Clé : R4.0.04 Révision : Detection of the singularities an calculation of a map of size of
More informationModelling 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 informationEXPONENTIAL 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 informationSYNCHRONOUS SEQUENTIAL CIRCUITS
CHAPTER SYNCHRONOUS SEUENTIAL CIRCUITS Registers an counters, two very common synchronous sequential circuits, are introuce in this chapter. Register is a igital circuit for storing information. Contents
More informationDebond 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 informationarxiv: 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 information5-4 Electrostatic Boundary Value Problems
11/8/4 Section 54 Electrostatic Bounary Value Problems blank 1/ 5-4 Electrostatic Bounary Value Problems Reaing Assignment: pp. 149-157 Q: A: We must solve ifferential equations, an apply bounary conitions
More informationSensors & Transducers 2015 by IFSA Publishing, S. L.
Sensors & Transucers, Vol. 184, Issue 1, January 15, pp. 53-59 Sensors & Transucers 15 by IFSA Publishing, S. L. http://www.sensorsportal.com Non-invasive an Locally Resolve Measurement of Soun Velocity
More informationLecture XII. where Φ is called the potential function. Let us introduce spherical coordinates defined through the relations
Lecture XII Abstract We introuce the Laplace equation in spherical coorinates an apply the metho of separation of variables to solve it. This will generate three linear orinary secon orer ifferential equations:
More informationTime-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 informationThe total derivative. Chapter Lagrangian and Eulerian approaches
Chapter 5 The total erivative 51 Lagrangian an Eulerian approaches The representation of a flui through scalar or vector fiels means that each physical quantity uner consieration is escribe as a function
More informationChapter 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 informationInfluence the Nozzle Shape on Local Heat Transfer in Impinging Jet
Issue 6, Volume 6, 12 Influence the Nozzle Shape on Local Heat Transfer in Impinging Jet M. Attalla an M. S. Ahme 1 Abstract The local Nusselt number istributions of circular nozzle on a heate flat plate
More informationOptimization of Geometries by Energy Minimization
Optimization of Geometries by Energy Minimization by Tracy P. Hamilton Department of Chemistry University of Alabama at Birmingham Birmingham, AL 3594-140 hamilton@uab.eu Copyright Tracy P. Hamilton, 1997.
More information(2012) , ISBN
Cregan, V. an O'Brien, Stephen B.G. an McKee, Sean (2012) Asymptotics of a small liqui rop on a cone an plate rheometer. In: Progress in Inustrial Mathematics at ECMI 2010. Mathematics in Inustry: Progress
More informationApplication 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 informationinflow 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 informationAnalytic 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 informationIntroduction to the Vlasov-Poisson system
Introuction to the Vlasov-Poisson system Simone Calogero 1 The Vlasov equation Consier a particle with mass m > 0. Let x(t) R 3 enote the position of the particle at time t R an v(t) = ẋ(t) = x(t)/t its
More informationStudy 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 informationLectures - 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θ x = f ( x,t) could be written as
9. Higher orer PDEs as systems of first-orer PDEs. Hyperbolic systems. For PDEs, as for ODEs, we may reuce the orer by efining new epenent variables. For example, in the case of the wave equation, (1)
More informationAdhesive 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 informationAn 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 informationOptimal LQR Control of Structures using Linear Modal Model
Optimal LQR Control of Structures using Linear Moal Moel I. Halperin,2, G. Agranovich an Y. Ribakov 2 Department of Electrical an Electronics Engineering 2 Department of Civil Engineering Faculty of Engineering,
More informationMULTISCALE 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 informationResilient Modulus Prediction Model for Fine-Grained Soils in Ohio: Preliminary Study
Resilient Moulus Preiction Moel for Fine-Graine Soils in Ohio: Preliminary Stuy by Teruhisa Masaa: Associate Professor, Civil Engineering Department Ohio University, Athens, OH 4570 Tel: (740) 59-474 Fax:
More informationTMA 4195 Matematisk modellering Exam Tuesday December 16, :00 13:00 Problems and solution with additional comments
Problem F U L W D g m 3 2 s 2 0 0 0 0 2 kg 0 0 0 0 0 0 Table : Dimension matrix TMA 495 Matematisk moellering Exam Tuesay December 6, 2008 09:00 3:00 Problems an solution with aitional comments The necessary
More informationINTEGRATED NUMERICAL ANALYSIS FOR COMPOSITE DAMAGE TOLERANCE
ICAS CONGRESS INTEGRATED NUMERICAL ANALSIS FOR COMPOSITE DAMAGE TOLERANCE Guinar Stéphane, Thévenet Pascal, Vinet Alain, Maison-Le Poëc Serge EADS-CCR Keywors: impact, amage tolerance, laminate composites
More informationDAMAGE DETECTIONS IN NONLINEAR VIBRATING THERMALLY LOADED STRUCTURES 1
11 th National Congress on Theoretical an Applie Mechanics, 2-5 Sept. 2009, Borovets, Bulgaria DAMAGE DETECTIONS IN NONLINEAR VIBRATING THERMALLY LOADED STRUCTURES 1 E. MANOACH Institute of Mechanics,
More informationLecture 2 Lagrangian formulation of classical mechanics Mechanics
Lecture Lagrangian formulation of classical mechanics 70.00 Mechanics Principle of stationary action MATH-GA To specify a motion uniquely in classical mechanics, it suffices to give, at some time t 0,
More informationAPPROXIMATE 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 informationNON-SMOOTH DYNAMICS USING DIFFERENTIAL-ALGEBRAIC EQUATIONS PERSPECTIVE: MODELING AND NUMERICAL SOLUTIONS. A Thesis PRIYANKA GOTIKA
NON-SMOOTH DYNAMICS USING DIFFERENTIAL-ALGEBRAIC EQUATIONS PERSPECTIVE: MODELING AND NUMERICAL SOLUTIONS A Thesis by PRIYANKA GOTIKA Submitte to the Office of Grauate Stuies of Texas A&M University in
More informationLie symmetry and Mei conservation law of continuum system
Chin. Phys. B Vol. 20, No. 2 20 020 Lie symmetry an Mei conservation law of continuum system Shi Shen-Yang an Fu Jing-Li Department of Physics, Zhejiang Sci-Tech University, Hangzhou 3008, China Receive
More informationThe Principle of Least Action and Designing Fiber Optics
University of Southampton Department of Physics & Astronomy Year 2 Theory Labs The Principle of Least Action an Designing Fiber Optics 1 Purpose of this Moule We will be intereste in esigning fiber optic
More informationComputing Exact Confidence Coefficients of Simultaneous Confidence Intervals for Multinomial Proportions and their Functions
Working Paper 2013:5 Department of Statistics Computing Exact Confience Coefficients of Simultaneous Confience Intervals for Multinomial Proportions an their Functions Shaobo Jin Working Paper 2013:5
More informationV-notched elements under mode II loading conditions
Structural Engineering and Mechanics, Vol. 49, No. 4 (2014) 499-508 DOI: http://dx.doi.org/10.12989/sem.2014.49.4.499 499 V-notched elements under mode loading conditions Alberto Sapora, Pietro Cornetti
More informationAsymptotics of a Small Liquid Drop on a Cone and Plate Rheometer
Asymptotics of a Small Liqui Drop on a Cone an Plate Rheometer Vincent Cregan, Stephen B.G. O Brien, an Sean McKee Abstract A cone an a plate rheometer is a laboratory apparatus use to measure the viscosity
More informationNonPAS: A Program for Nonlinear Analysis of Flexible Pavements
International Journal of Integrate Engineering, Vol. 7 No. (05) pp. -8 NonPAS: A Program for Nonlinear Analysis of Flexible Pavements Ali eza Ghanizaeh,*, Arash Ziaie, Department of Civil Engineering,
More informationNonlinear Dielectric Response of Periodic Composite Materials
onlinear Dielectric Response of Perioic Composite aterials A.G. KOLPAKOV 3, Bl.95, 9 th ovember str., ovosibirsk, 639 Russia the corresponing author e-mail: agk@neic.nsk.su, algk@ngs.ru A. K.TAGATSEV Ceramics
More informationThe 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 informationImpurities in inelastic Maxwell models
Impurities in inelastic Maxwell moels Vicente Garzó Departamento e Física, Universia e Extremaura, E-671-Baajoz, Spain Abstract. Transport properties of impurities immerse in a granular gas unergoing homogenous
More informationMonitoring of Mechanical Properties of Gypsum Using the Impulse Excitation Method
Monitoring of Mechanical Properties of Gypsum Using the Impulse Excitation Metho RICHARD ŤOUPEK, TOMÁŠ PLACHÝ, MICHAL POLÁK, PAVEL TESÁREK Department of Mechanics, Faculty of Civil Engineering Czech Technical
More informationHomework 7 Due 18 November at 6:00 pm
Homework 7 Due 18 November at 6:00 pm 1. Maxwell s Equations Quasi-statics o a An air core, N turn, cylinrical solenoi of length an raius a, carries a current I Io cos t. a. Using Ampere s Law, etermine
More informationShear stiffness of granular material at small strains: does it depend on grain size?
Title Shear stiffness of granular material at small strains: oes it epen on grain size? Author(s) Yang, J; Gu, X Citation Géotechnique, 213, v. 63 n. 2, p. 165-179 Issue Date 213 URL http://hl.hanle.net/1722/18899
More information2Algebraic ONLINE PAGE PROOFS. foundations
Algebraic founations. Kick off with CAS. Algebraic skills.3 Pascal s triangle an binomial expansions.4 The binomial theorem.5 Sets of real numbers.6 Surs.7 Review . Kick off with CAS Playing lotto Using
More informationMagnetic field generated by current filaments
Journal of Phsics: Conference Series OPEN ACCESS Magnetic fiel generate b current filaments To cite this article: Y Kimura 2014 J. Phs.: Conf. Ser. 544 012004 View the article online for upates an enhancements.
More informationTHE ACCURATE ELEMENT METHOD: A NEW PARADIGM FOR NUMERICAL SOLUTION OF ORDINARY DIFFERENTIAL EQUATIONS
THE PUBISHING HOUSE PROCEEDINGS O THE ROMANIAN ACADEMY, Series A, O THE ROMANIAN ACADEMY Volume, Number /, pp. 6 THE ACCURATE EEMENT METHOD: A NEW PARADIGM OR NUMERICA SOUTION O ORDINARY DIERENTIA EQUATIONS
More informationTopological Sensitivity Analysis for Three-dimensional Linear Elasticity Problem
Topological Sensitivity Analysis for Three-imensional Linear Elasticity Problem A.A. Novotny, R.A. Feijóo, E. Taroco Laboratório Nacional e Computação Científica LNCC/MCT, Av. Getúlio Vargas 333, 25651-075
More informationA new identification method of the supply hole discharge coefficient of gas bearings
Tribology an Design 95 A new ientification metho of the supply hole ischarge coefficient of gas bearings G. Belforte, F. Colombo, T. Raparelli, A. Trivella & V. Viktorov Department of Mechanics, Politecnico
More informationA 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 informationProf. Dr. Ibraheem Nasser electric_charhe 9/22/2017 ELECTRIC CHARGE
ELECTRIC CHARGE Introuction: Orinary matter consists of atoms. Each atom consists of a nucleus, consisting of protons an neutrons, surroune by a number of electrons. In electricity, the electric charge
More informationApproaches 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 informationEXACT TRAVELING WAVE SOLUTIONS FOR A NEW NON-LINEAR HEAT TRANSFER EQUATION
THERMAL SCIENCE, Year 017, Vol. 1, No. 4, pp. 1833-1838 1833 EXACT TRAVELING WAVE SOLUTIONS FOR A NEW NON-LINEAR HEAT TRANSFER EQUATION by Feng GAO a,b, Xiao-Jun YANG a,b,* c, an Yu-Feng ZHANG a School
More informationPARALLEL-PLATE CAPACITATOR
Physics Department Electric an Magnetism Laboratory PARALLEL-PLATE CAPACITATOR 1. Goal. The goal of this practice is the stuy of the electric fiel an electric potential insie a parallelplate capacitor.
More informationThe 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 informationGeneralization of the persistent random walk to dimensions greater than 1
PHYSICAL REVIEW E VOLUME 58, NUMBER 6 DECEMBER 1998 Generalization of the persistent ranom walk to imensions greater than 1 Marián Boguñá, Josep M. Porrà, an Jaume Masoliver Departament e Física Fonamental,
More informationIn the usual geometric derivation of Bragg s Law one assumes that crystalline
Diffraction Principles In the usual geometric erivation of ragg s Law one assumes that crystalline arrays of atoms iffract X-rays just as the regularly etche lines of a grating iffract light. While this
More informationDusty Plasma Void Dynamics in Unmoving and Moving Flows
7 TH EUROPEAN CONFERENCE FOR AERONAUTICS AND SPACE SCIENCES (EUCASS) Dusty Plasma Voi Dynamics in Unmoving an Moving Flows O.V. Kravchenko*, O.A. Azarova**, an T.A. Lapushkina*** *Scientific an Technological
More informationPower Generation and Distribution via Distributed Coordination Control
Power Generation an Distribution via Distribute Coorination Control Byeong-Yeon Kim, Kwang-Kyo Oh, an Hyo-Sung Ahn arxiv:407.4870v [math.oc] 8 Jul 204 Abstract This paper presents power coorination, power
More informationDiagonalization of Matrices Dr. E. Jacobs
Diagonalization of Matrices Dr. E. Jacobs One of the very interesting lessons in this course is how certain algebraic techniques can be use to solve ifferential equations. The purpose of these notes is
More informationExtinction, σ/area. Energy (ev) D = 20 nm. t = 1.5 t 0. t = t 0
Extinction, σ/area 1.5 1.0 t = t 0 t = 0.7 t 0 t = t 0 t = 1.3 t 0 t = 1.5 t 0 0.7 0.9 1.1 Energy (ev) = 20 nm t 1.3 Supplementary Figure 1: Plasmon epenence on isk thickness. We show classical calculations
More informationOn the use of leaky modes in open waveguides for the sound propagation modeling in street canyons
On the use of leaky moes in open waveguies for the soun propagation moeling in street canyons Arien Pelat, a Simon Félix, an Vincent Pagneux LAUM, CNRS, Université u Maine, avenue Olivier Messiaen, 7285
More informationAssignment 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 informationEVALUATION 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 informationSingle Arm, Centrifugal, Water Turbine for Low Head and Low Flow Application: Part 1- Theory and Design
Energy an Power 2018, 8(2): 51-55 DOI: 10.5923/j.ep.20180802.03 Single Arm, Centrifugal, Water Turbine for Low ea an Low Flow Application: Part 1- Theory an Design Kiplangat C. Kononen 1, Augustine B.
More informationStatic analysis of beams on elastic foundation by the method of discrete singular convolution
International Journal of Engineering & Applie Sciences (IJEAS) Beir Agöz, Kair Mercan, Çiğem Demir, Ömer Civale Vol.8, Issue 3(016) 67-73 Static analysis of beams on elastic founation by the metho of iscrete
More informationTensile behaviour of anti-symmetric CFRP composite
Available online at www.sciencedirect.com Procedia Engineering 1 (211) 1865 187 ICM11 Tensile behaviour of anti-symmetric CFRP composite K. J. Wong a,b, *, X. J. Gong a, S. Aivazzadeh a, M. N. Tamin b
More informationDamage detection of shear building structure based on FRF response variation
, pp.18-5 http://x.oi.org/10.1457/astl.013.3.05 Damage etection of shear builing structure base on FRF response variation Hee-Chang Eun 1,*, Su-Yong Park 1, Rae-Jung im 1 1 Dept. of Architectural Engineering,
More informationSchrödinger s equation.
Physics 342 Lecture 5 Schröinger s Equation Lecture 5 Physics 342 Quantum Mechanics I Wenesay, February 3r, 2010 Toay we iscuss Schröinger s equation an show that it supports the basic interpretation of
More informationSPE Copyright 1999, Society of Petroleum Engineers Inc.
SPE 664 Effect of Flow Through a Choke Valve on Emulsion Stability M.J. van er Zane, SPE, K.R. van Heuven, J.H. Muntinga, SPE, an W.M.G.T. van en Broek, SPE, Delft University of Technology Copyright 1999,
More informationINTRODUCTION & PHASE SYSTEM
INTRODUCTION & PHASE SYSTEM Dr. Professor of Civil Engineering S. J. College of Engineering, Mysore 1.1 Geotechnical Engineering Why? 1. We are unable to buil castles in air (yet)! 2. Almost every structure
More informationSemiclassical analysis of long-wavelength multiphoton processes: The Rydberg atom
PHYSICAL REVIEW A 69, 063409 (2004) Semiclassical analysis of long-wavelength multiphoton processes: The Ryberg atom Luz V. Vela-Arevalo* an Ronal F. Fox Center for Nonlinear Sciences an School of Physics,
More informationPlacement and tuning of resonance dampers on footbridges
Downloae from orbit.tu.k on: Jan 17, 19 Placement an tuning of resonance ampers on footbriges Krenk, Steen; Brønen, Aners; Kristensen, Aners Publishe in: Footbrige 5 Publication ate: 5 Document Version
More informationAngles-Only Orbit Determination Copyright 2006 Michel Santos Page 1
Angles-Only Orbit Determination Copyright 6 Michel Santos Page 1 Abstract This ocument presents a re-erivation of the Gauss an Laplace Angles-Only Methos for Initial Orbit Determination. It keeps close
More information