Analysis of composites with multiple rigid-line reinforcements by the BEM
|
|
- Julie Cannon
- 5 years ago
- Views:
Transcription
1 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, Gliwice, Polad Abstract This paper presets the boudary elemet method (BEM) for aalysis of plates cotaiig multiple rigid-lie reiforcemets. The accuracy of the method is verified by comparig computed stresses i the eighbourhood of a rigid fibre i a ifiite plate subjected to tesio with aalytical solutios. The proposed approach is applied for aalysis of composites reiforced by multiple rigid fibres. The ifluece of the distace betwee the fibres o stress distributio i the matrix ad the effective Youg modulus is aalyzed. Keywords: boudary elemet methods, composites, microstructures, homogeizatio, elasticity 1. Itroductio The aim of usig fibres i composite materials is to icrease stregth ad stiffess of structures. It is kow that these two properties are difficult to improve simultaeously [3]. Stiffess depeds o overall properties of the structure ad stregth o microstructure. If the stiffess of fibers is much greater tha the matrix the the fibers ca be modelled as rigid stiffeers i a deformable matrix. The examples of such atural composites are: silk, acre ad eamel [6]. The aalysis of composites with may fibers requires computatioal methods. Li ad Tig [4] aalyzed a lie iclusio i a aisotropic elastic ifiite plate subjected to uiform loadig at ifiity. They used the Stroh formalism to calculate the displacemet ad stress fields for the rigid ad elastic iclusio. Pigle et al. [6] used the duality priciple to derive rigid lie iclusio solutios from the crack solutios. They computed stress fields aroud a rigid lie iclusio ad derived a compliace cotributio tesor for a sigle ad multiple lie iclusios. Gorbatikh et al. [3] determied the relatio betwee stress itesity factors at the tips of rigid lie iclusios ad effective compliace of the material. Liu et al. [5] applied the fast multipole BEM to aalyze composites reiforced by aotubes. They calculated effective material properties of composites modelled as three-dimesioal structures havig large umber of degrees of freedom. The formulatio ad applicatio of the BEM for a sigle rigid fibre i a fiite plate was preseted by Fedeliski [2]. I the preset work the method is exteded ad applied to solids with multiple rigid-lie reiforcemets. 2. Boudary elemet method for plates with rigid fibres 2.1. Boudary itegral equatios for a plate with fibres Cosider a plate made of homogeous, isotropic ad liearelastic material. The boudary of the plate is deoted by ad its domai by Ω (Fig. 1). The plate is statically loaded alog the exteral boudary by boudary tractios t j ad the domai Ω by body forces f j. The relatio betwee the loadig of the plate ad its displacemets ca be expressed by the Somigliaa idetity [1] cij ( x') u j ( x') + Tij ( x', x) u j( x) d ( x) =, (1) Uij ( x', x) t j ( x) d ( x) + Uij( x', X ) f j( X ) dω( X ) Ω where: x is the collocatio poit, for which the above itegral equatio is applied, x is the boudary poit, X is the domai poit ad c ij is a costat, which depeds o the positio of the poit x, U ij ad T ij are the Kelvi fudametal solutios of elastostatics. I equatios the Eistei summatio covetio is used ad the idices for two-dimesioal problems have values i,j=1,2. Figure 1. Elastic plate - loadig ad displacemets Assume that the plate is reiforced by straight, thi ad rigid fibres, which are perfectly boded to the matrix. If the plate is deformed the forces of iteractio occur alog the lies of attachmet of fibers (Fig. 2). These forces of iteractio ca be treated as particular body forces actig alog the lies i the domai of the body. The boudary itegral equatio (1) for the plate loaded by boudary tractios ad forces of iteractio of fibres has the form [2] *Ackowledgmet: The scietific research is fiaced by the Miistry of Sciece ad Higher Educatio of Polad i years
2 cij ( x') u j ( x') + Tij ( x', x) u j( x) d ( x) = N, (2) Uij ( x', x) t j( x) d ( x) + Uij ( x', X ) t j ( X ) d( X ) = 1 where: N is the umber of fibres, is the lie of attachmet ad t j are the forces of iteractio. The cosidered structure ad therefore each fibre is i equilibrium. The forces actig o each fibre (Fig. 4) should satisfied the followig equilibrium equatios t1 ( x ) d ( x ) = 0, (5) t2 ( x ) d ( x ) = 0, (6) [ t1 ( x)si α + t2 ( x)cos α ] r( x) d ( x) = 0, (7) The last equatio (7) is the equatio of momets of forces with respect of the fibre tip x o. Figure 2. Elastic plate with fibres 2.2. Displacemets ad equatios of equilibrium for stiff fibres Figure 3. Displacemets of the rigid fibre 2.3. Numerical implemetatio of the method The boudary of the plate ad the fibres are divided ito boudary elemets (Fig. 5). I the developed computer code 3- ode quadratic boudary elemets are used. Alog the exteral boudary of the plate the variatios of coordiates, displacemets ad tractios, ad alog the fibres the variatios of iteractio forces are iterpolated. The boudary itegral equatios (2) are used for odes alog the exteral boudary ad the fibres. The displacemets of fibre odes ca be expressed by the displacemets of fibre tips ad their agles of rotatio, by usig Eqs (3) ad (4). These equatios ca writte i the followig matrix form u= Iu f, (8) where the matrix u cotais the compoets of displacemets of fibre odes, the matrix I depeds o the positio of odes ad the matrix u f cotais compoets of displacemets of fibre tips ad their agles of rotatio. Figure 4. Forces actig o the fibre Deformatios of the matrix ifluece displacemets of fibres. The displacemet of a arbitrary poit x of the fibre ca be expressed by the displacemets of the fibre tip x o ad the agle of rotatio of the fibre ϕ (Fig. 3). For small agles of rotatio of fibres, the compoets of displacemets of a arbitrary poit of the fibre are expressed i the form u1( x) = u1( xo) ϕr( x)siα, (3) u2( x) = u2( xo ) + ϕr( x)cosα, (4) Figure 5. Discretizatio of the matrix ad fibres by quadratic boudary elemets The equilibrium equatios for fibres (5), (6) ad (7) ca be writte i the matrix form Et f = 0, (9) where the matrix E depeds o the positio of fibre odes ad the matrix t f cotais odal values of compoets of tractios i fibres. The matrix E is obtaied by itegratio of expressios i Eqs (5), (6) i (7), by assumig quadratic variatios of forces alog the fibres. Because the equilibrium equatios have very simple forms, the itegrals are computed aalytically. where: α is the iitial agle betwee the fibre ad the axis x 1 of the global coordiate system, r is the distace betwee the poit x ad the fibre tip x o.
3 Boudary itegral equatios (2), supplied with Eqs (8) ad (9) ca be writte i the matrix form H 0 Gee G ee ef ue te H fe I = G fe G ff, (10) u 0 0 f t 0 f E where the submatrices with the idex e are related to the exteral boudary ad the submatrices deoted by the idex f are related to the fibres. The submatrices H ad G deped o boudary itegrals of fudametal solutios, shape fuctios ad are itegrated umerically by usig the Gauss quadrature. Next the system of algebraic equatios is rearraged. The ukow quatities are o the left side of the equatio ad the kow quatities o the right side of the equatio. The first modificatio refers to the ukow iteractio forces t f Hee Gef 0 u e G ee H fe G ff I t f = G fe [ te]. (11) E u f I the fial modificatio the kow ad ukow boudary coditios are rearraged. The modified matrix equatio is solved ad the ukow displacemets ad tractios alog the exteral boudary ad the displacemets ad iteractio forces for fibres are obtaied. 3. Numerical examples To demostrate the accuracy of the method ad possible applicatios three umerical examples are solved Rigid fibre i a ifiite plate A ifiite plate with a rigid fibre of legth 2l is subjected to the parallel loadig q 1 or to the perpedicular loadig q 2, as show i Fig. 6. The structure is modelled as a fiite square plate with a rigid fibre ad the dimesios of the plate are 10 times larger tha the fibre. The material of the plate has the Poisso ratio ν=0.25 ad is i plae stress state. The rigid fibre is divided ito 20 boudary elemets ad the square plate ito 160 boudary elemets. Stresses are computed at 245 iteral poits i the eighbourhood of the fibre. This field is marked as a grey square i Fig. 6. The cotour plot of ormalized stresses σ 11 /q 1 for the parallel loadig q 1 ad ormalized stresses σ 22 /q 2 for the perpedicular loadig q 2 are show i Fig. 7. Figure 7: Normal stresses i the viciity of the rigid fibre: parallel loadig ormalized stress σ 11 /q 1, perpedicular loadig ormalized stress σ 22 /q 2 For the parallel loadig stresses σ 11 have the smallest values alog the fibre ad the largest values alog the extesio of the fibre. For the perpedicular loadig the stresses σ 22 are more uiformly distributed i the aalyzed field. The maximum values of stresses exceed the applied tesio by about 2%. The distributio of stresses agree very well with stresses computed aalytically ad preseted by Pigle et al. [6] Two rigid fibres i a ifiite plate Two parallel rigid fibres i a ifiite plate are subjected to tesio q 1 (Fig. 8). The legth of the fibres is 2l ad the distace betwee the fibres is 2d 2. The same material properties ad discretizatio are assumed as i the previous example 3.1. The ifluece of the distace betwee the fibres o stress distributio is studied. The ormalized stresses σ 11 /q 1, i the marked field show i Fig. 8, for the distaces d 2 /l=0.5 ad d 2 /l=1.0 are preseted i Fig. 9. For the cosidered distaces the iteractio of fibres has small ifluece o stress distributio. Figure 6: Rigid fibre i a ifiite plate dimesios ad loadig Figure 8. Two rigid fibres i a ifiite plate dimesios ad loadig
4 Figure 11. Rigid fibres i a rectagular plate iitial shape (dashed lie) ad deformed shape (cotiuous lie) Figure 9. Ifluece of the distace betwee two fibres o the stress distributio σ 11 /q 1 : d 2 /l=0.5, d 2 /l= Rigid fibres i a rectagular plate A rectagular plate of legth 2b ad height 2h cotais 13 rigid fibres of legth 2l, as show i Fig. 10. The horizotal distace betwee cetres of eighbourig fibres is d 1 ad the vertical distace is d 2. The ratios of dimesios are: b/l=5, h/l=4, d 1 /l=3 ad d 2 /l=1.6. The material of the plate has the Poisso ratio ν=0.3 ad is i plae strai state. The plate is subjected to the horizotal loadig q 1. Each rigid fibre is divided ito 8 boudary elemets ad the exteral boudary ito 72 boudary elemets. The iitial shape ad the deformed shape of the plate are show i Fig. 11. The stress distributios i the matrix, i the marked field i Fig. 10, are show i Fig. 12. The relative effective Youg modulus is computed as E c /E m =1.345, where E c ad E m are the Youg modulus of the composite ad the matrix, respectively. The effective Poisso ratio of the composite is the same as for the matrix. c) Figure 12. Stress distributios i the matrix: σ 11 /q 1, σ 22 /q 1, c) σ 12 /q 1 Figure 10: Rigid fibres i a rectagular plate dimesios ad loadig
5 The ifluece of the horizotal ad vertical distace betwee the fibres o the effective Youg modulus is aalyzed. The dimesios of fibres are costat ad the dimesios of the plate are chaged proportioally to the distace betwee the fibres. The depedece of the relative Youg modulus o the horizotal ad vertical distace is preseted i Fig. 13. If the horizotal or vertical distace are smaller tha half of the legth of the fibre the stiffess of the composite sigificatly icreases. For the icreasig distace betwee the fibres the stiffess of the composite teds to the stiffess of the matrix. Refereces [1] Becker A.A., The boudary elemet method i egieerig. A complete course. McGraw-Hill Book Compay, Lodo, [2] Fedeliski, P., Boudary elemet method for aalysis of elastic structures with rigid fibers, Eg. Modell., 32, pp , 2006 (i Polish). [3] Gorbatikh, L., Lomov S.V. ad Verpoest, I., Relatio betwee elastic properties ad stress itesity factors for composites with rigid-lie reiforcemets, It. J. Fract., 161, pp , [4] Li, Q. ad Tig, T.C.T., Lie iclusios i aisotropic elastic solids, Tras. ASME, J. Appl. Mech., 56, pp , [5] Liu, Y., Nishimura, N. ad Otai, Y., Large-scale modelig of carbo-aotube composites by a fast multipole boudary elemet method, Comp. Mater. Sc., 34, pp , [6] Pigle, P., Sherwood, J. ad Gorbatikh, L., Properties of rigid-lie iclusios as buildig blocks of aturally occurrig composites, Composites Sc. Techol., 68, pp , Figure 13. Ifluece of the distace betwee the fibres o the relative effective Youg modulus E c /E m : horizotal distace d 1 /2l, vertical distace d 2 /2l 4. Coclusios The boudary elemet method is a very efficiet method for modellig elastic composites with rigid iclusios, particularly fibres. The modellig of composites is sigificatly simplified i compariso to domai methods, for example the fiite elemet method, because odes are situated oly alog the exteral boudary ad fibers. Positios, shapes ad dimesios of fibers ad their umber ca be very simply modified. The compariso of stresses for the fibre i the ifiite plate obtaied by the preset method ad the aalytical method shows that the proposed approach gives very accurate results. For distaces betwee the fibres smaller tha the half legth of fibres a ifluece of iteractio of fibres o stress distributio i the matrix ad sigificat icrease of stiffess ca be observed.
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 informationMETHOD 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 informationDiscrete Orthogonal Moment Features Using Chebyshev Polynomials
Discrete Orthogoal Momet Features Usig Chebyshev Polyomials R. Mukuda, 1 S.H.Og ad P.A. Lee 3 1 Faculty of Iformatio Sciece ad Techology, Multimedia Uiversity 75450 Malacca, Malaysia. Istitute of Mathematical
More informationMATHEMATICAL 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 informationANALYSIS OF THERMOELASTIC STRESSES IN LAYERED PLATES
6t Iteratioal DAAAM Baltic Coferece "INDUSTRIAL ENGINEERING" 4-6 April 8, Talli, Estoia ANALYSIS OF THERMOELASTIC STRESSES IN LAYERED PLATES Kõo, J. & Valgur, J. Abstract: Tis paper deals wit termoelastic
More informationDYNAMIC 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 informationNumerical 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 informationEVALUATION 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 informationCO-LOCATED DIFFUSE APPROXIMATION METHOD FOR TWO DIMENSIONAL INCOMPRESSIBLE CHANNEL FLOWS
CO-LOCATED DIFFUSE APPROXIMATION METHOD FOR TWO DIMENSIONAL INCOMPRESSIBLE CHANNEL FLOWS C.PRAX ad H.SADAT Laboratoire d'etudes Thermiques,URA CNRS 403 40, Aveue du Recteur Pieau 86022 Poitiers Cedex,
More information2C09 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 informationDeterministic Model of Multipath Fading for Circular and Parabolic Reflector Patterns
To appear i the Proceedigs of the 5 IEEE outheastco, (Ft. Lauderdale, FL), April 5 Determiistic Model of Multipath Fadig for Circular ad Parabolic Reflector Patters Dwight K. Hutcheso dhutche@clemso.edu
More informationSECTION 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 informationANALYSIS 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 informationAdded 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 informationChapter 9: Numerical Differentiation
178 Chapter 9: Numerical Differetiatio Numerical Differetiatio Formulatio of equatios for physical problems ofte ivolve derivatives (rate-of-chage quatities, such as velocity ad acceleratio). Numerical
More informationA numerical Technique Finite Volume Method for Solving Diffusion 2D Problem
The Iteratioal Joural Of Egieerig d Sciece (IJES) Volume 4 Issue 10 Pages PP -35-41 2015 ISSN (e): 2319 1813 ISSN (p): 2319 1805 umerical Techique Fiite Volume Method for Solvig Diffusio 2D Problem 1 Mohammed
More informationFailure 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 informationPILOT 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 informationGROUND 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 information6.1 Analysis of frequency selective surfaces
6.1 Aalysis of frequecy selective surfaces Basic theory I this paragraph, reflectio coefficiet ad trasmissio coefficiet are computed for a ifiite periodic frequecy selective surface. The attetio is tured
More informationPrinciple 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 informationDynamical stability and parametrical vibrations of the laminated plates with complex shape
(3) 75 88 Dyamical stability ad parametrical vibratios of the lamiated plates with comple shape Abstract The problem of oliear vibratios ad stability aalysis for the symmetric lamiated plates with comple
More informationStreamfunction-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 informationMONITORING 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 informationStochastic finite element method study on temperature sensitivity of the steel telecommunications towers
CMM-11 Computer Methods i Mechaics 9 1 May 11, Warsaw, Polad Stochastic fiite elemet method study o temperature sesitivity of the steel telecommuicatios towers Marci M. Kamińsi 1 ad Jace Szafra 1 Professor,
More informationSize, 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 informationTHEORETICAL RESEARCH REGARDING ANY STABILITY THEOREMS WITH APPLICATIONS. Marcel Migdalovici 1 and Daniela Baran 2
ICSV4 Cairs Australia 9- July, 007 THEORETICAL RESEARCH REGARDING ANY STABILITY THEOREMS WITH APPLICATIONS Marcel Migdalovici ad Daiela Bara Istitute of Solid Mechaics, INCAS Elie Carafoli, 5 C-ti Mille
More informationTaylor polynomial solution of difference equation with constant coefficients via time scales calculus
TMSCI 3, o 3, 129-135 (2015) 129 ew Treds i Mathematical Scieces http://wwwtmscicom Taylor polyomial solutio of differece equatio with costat coefficiets via time scales calculus Veysel Fuat Hatipoglu
More informationThe 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 informationNumerical Solution of the Two Point Boundary Value Problems By Using Wavelet Bases of Hermite Cubic Spline Wavelets
Australia Joural of Basic ad Applied Scieces, 5(): 98-5, ISSN 99-878 Numerical Solutio of the Two Poit Boudary Value Problems By Usig Wavelet Bases of Hermite Cubic Splie Wavelets Mehdi Yousefi, Hesam-Aldie
More informationDECOMPOSITION METHOD FOR SOLVING A SYSTEM OF THIRD-ORDER BOUNDARY VALUE PROBLEMS. Park Road, Islamabad, Pakistan
Mathematical ad Computatioal Applicatios, Vol. 9, No. 3, pp. 30-40, 04 DECOMPOSITION METHOD FOR SOLVING A SYSTEM OF THIRD-ORDER BOUNDARY VALUE PROBLEMS Muhammad Aslam Noor, Khalida Iayat Noor ad Asif Waheed
More informationElastic 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 informationInelastic 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 informationNumerical Methods in Fourier Series Applications
Numerical Methods i Fourier Series Applicatios Recall that the basic relatios i usig the Trigoometric Fourier Series represetatio were give by f ( x) a o ( a x cos b x si ) () where the Fourier coefficiets
More informationA 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 informationPOWER SERIES SOLUTION OF FIRST ORDER MATRIX DIFFERENTIAL EQUATIONS
Joural of Applied Mathematics ad Computatioal Mechaics 4 3(3) 3-8 POWER SERIES SOLUION OF FIRS ORDER MARIX DIFFERENIAL EQUAIONS Staisław Kukla Izabela Zamorska Istitute of Mathematics Czestochowa Uiversity
More informationOrthogonal Gaussian Filters for Signal Processing
Orthogoal Gaussia Filters for Sigal Processig Mark Mackezie ad Kiet Tieu Mechaical Egieerig Uiversity of Wollogog.S.W. Australia Abstract A Gaussia filter usig the Hermite orthoormal series of fuctios
More informationThe Method of Least Squares. To understand least squares fitting of data.
The Method of Least Squares KEY WORDS Curve fittig, least square GOAL To uderstad least squares fittig of data To uderstad the least squares solutio of icosistet systems of liear equatios 1 Motivatio Curve
More informationCHAPTER 8 SYSTEMS OF PARTICLES
CHAPTER 8 SYSTES OF PARTICLES CHAPTER 8 COLLISIONS 45 8. CENTER OF ASS The ceter of mass of a system of particles or a rigid body is the poit at which all of the mass are cosidered to be cocetrated there
More informationThe target reliability and design working life
Safety ad Security Egieerig IV 161 The target reliability ad desig workig life M. Holický Kloker Istitute, CTU i Prague, Czech Republic Abstract Desig workig life ad target reliability levels recommeded
More informationwavelet collocation method for solving integro-differential equation.
IOSR Joural of Egieerig (IOSRJEN) ISSN (e): 5-3, ISSN (p): 78-879 Vol. 5, Issue 3 (arch. 5), V3 PP -7 www.iosrje.org wavelet collocatio method for solvig itegro-differetial equatio. Asmaa Abdalelah Abdalrehma
More informationBoundary 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 informationOn Strongly Consistent Finite Dierence Approximations
D.Michels et al. (KAUST,JINR,SSU) Strogly Cosistet Approximatios 17 April 2018 1 / 25 O Strogly Cosistet Fiite Dierece Approximatios Domiik Michels 1, Vladimir Gerdt 2, Dmitry Lyakhov 1, ad Yuri Blikov
More informationA 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 informationDYNAMIC STUDY OF THE CONTACT ANGLE HYSTERESIS IN THE PRESENCE OF PERIODIC DEFECTS *
11 th Natioal Cogress o Theoretical ad Applied Mechaics, 2-5 Sept. 2009, Borovets, Bulgaria DYNAMIC STUDY OF THE CONTACT ANGLE HYSTERESIS IN THE PRESENCE OF PERIODIC DEFECTS * STANIMIR ILIEV, NINA PESHEVA
More informationNumerical Conformal Mapping via a Fredholm Integral Equation using Fourier Method ABSTRACT INTRODUCTION
alaysia Joural of athematical Scieces 3(1): 83-93 (9) umerical Coformal appig via a Fredholm Itegral Equatio usig Fourier ethod 1 Ali Hassa ohamed urid ad Teh Yua Yig 1, Departmet of athematics, Faculty
More informationa 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 informationFREE 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 informationPAijpam.eu ON TENSOR PRODUCT DECOMPOSITION
Iteratioal Joural of Pure ad Applied Mathematics Volume 103 No 3 2015, 537-545 ISSN: 1311-8080 (prited versio); ISSN: 1314-3395 (o-lie versio) url: http://wwwijpameu doi: http://dxdoiorg/1012732/ijpamv103i314
More informationUsing 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 informationNUMERICAL METHOD FOR SINGULARLY PERTURBED DELAY PARABOLIC PARTIAL DIFFERENTIAL EQUATIONS
THERMAL SCIENCE, Year 07, Vol., No. 4, pp. 595-599 595 NUMERICAL METHOD FOR SINGULARLY PERTURBED DELAY PARABOLIC PARTIAL DIFFERENTIAL EQUATIONS by Yula WANG *, Da TIAN, ad Zhiyua LI Departmet of Mathematics,
More informationAvoidance of numerical singularities in free vibration analysis of Euler-Bernoulli beams using Green functions
WSEAS TRASACTIOS o APPLIED ad THEORETICAL MECHAICS Goraka Štimac Ročević, Braimir Ročević, Ate Skoblar, Saji Braut Avoidace of umerical sigularities i free vibratio aalysis of Euler-Beroulli beams usig
More informationCEE 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 informationTR Use of the Implicit HHT-I3 and the Explicit ADAMS Methods with the Absolute Nodal Coordinate Formulation
R-007-05 Use of the Implicit HH-I3 ad the Explicit ADAMS Methods with the Absolute Nodal Coordiate Formulatio Bassam Hussei, Da Negrut, Ahmed A. Shabaa August 007 Abstract his ivestigatio is cocered with
More informationDETERMINATION OF HIGHER ORDER COEFFICIENTS FOR A CRACK PARALLEL TO AN INTERFACE
Clemso Uiversity TigerPrits All Theses Theses 8- DETERMINATION OF HIGHER ORDER COEFFICIENTS FOR A CRACK PARALLEL TO AN INTERFACE Vekata lakshma Kolluru Clemso Uiversity, lakshma.kvl@gmail.com Follow this
More informationTMA4205 Numerical Linear Algebra. The Poisson problem in R 2 : diagonalization methods
TMA4205 Numerical Liear Algebra The Poisso problem i R 2 : diagoalizatio methods September 3, 2007 c Eiar M Røquist Departmet of Mathematical Scieces NTNU, N-749 Trodheim, Norway All rights reserved A
More informationTHE SYSTEMATIC AND THE RANDOM. ERRORS - DUE TO ELEMENT TOLERANCES OF ELECTRICAL NETWORKS
R775 Philips Res. Repts 26,414-423, 1971' THE SYSTEMATIC AND THE RANDOM. ERRORS - DUE TO ELEMENT TOLERANCES OF ELECTRICAL NETWORKS by H. W. HANNEMAN Abstract Usig the law of propagatio of errors, approximated
More informationOn Random Line Segments in the Unit Square
O Radom Lie Segmets i the Uit Square Thomas A. Courtade Departmet of Electrical Egieerig Uiversity of Califoria Los Ageles, Califoria 90095 Email: tacourta@ee.ucla.edu I. INTRODUCTION Let Q = [0, 1] [0,
More informationFor 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 informationModeling the Radiation Field of Normal-Beam Ultrasonic Transducers by Multi-Gaussian Beam Technique
Modelig the Radiatio Field of Normal-Beam Ultrasoic Trasducers by Multi-Gaussia Beam Techique Ami Yaghootia, Farhag Hoarvar, ad Mehdi Zeighami NDE Lab., Faculty of Mechaical Egieerig, K. N. Toosi Uiversity
More informationImproving the Localization of Eigenvalues for Complex Matrices
Applied Mathematical Scieces, Vol. 5, 011, o. 8, 1857-1864 Improvig the Localizatio of Eigevalues for Complex Matrices P. Sargolzaei 1, R. Rakhshaipur Departmet of Mathematics, Uiversity of Sista ad Baluchesta
More information3 Balance equations ME338A CONTINUUM MECHANICS
ME338A CONTINUUM MECHANICS lecture otes 1 thursy, may 1, 28 Basic ideas util ow: kiematics, i.e., geeral statemets that characterize deformatio of a material body B without studyig its physical cause ow:
More informationAdvanced Analysis. Min Yan Department of Mathematics Hong Kong University of Science and Technology
Advaced Aalysis Mi Ya Departmet of Mathematics Hog Kog Uiversity of Sciece ad Techology September 3, 009 Cotets Limit ad Cotiuity 7 Limit of Sequece 8 Defiitio 8 Property 3 3 Ifiity ad Ifiitesimal 8 4
More informationBoundary 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 informationInverse Nodal Problems for Differential Equation on the Half-line
Australia Joural of Basic ad Applied Scieces, 3(4): 4498-4502, 2009 ISSN 1991-8178 Iverse Nodal Problems for Differetial Equatio o the Half-lie 1 2 3 A. Dabbaghia, A. Nematy ad Sh. Akbarpoor 1 Islamic
More informationU8L1: Sec Equations of Lines in R 2
MCVU U8L: Sec. 8.9. Equatios of Lies i R Review of Equatios of a Straight Lie (-D) Cosider the lie passig through A (-,) with slope, as show i the diagram below. I poit slope form, the equatio of the lie
More informationProblem Cosider the curve give parametrically as x = si t ad y = + cos t for» t» ß: (a) Describe the path this traverses: Where does it start (whe t =
Mathematics Summer Wilso Fial Exam August 8, ANSWERS Problem 1 (a) Fid the solutio to y +x y = e x x that satisfies y() = 5 : This is already i the form we used for a first order liear differetial equatio,
More informationTESTING OF THE FORCES IN CABLE OF SUSPENSION STRUCTURE AND BRIDGES
TSTING OF TH FORCS IN CABL OF SUSPNSION STRUCTUR AND BRIDGS Zhou, M. 1, Liu, Z. ad Liu, J. 1 College of the Muicipal Techology, Guagzhou Uiversity, Guagzhou. Guagzhou Muicipal ad Ladscape gieerig Quality
More informationSHAFTS: 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 informationComparison 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 information577. 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 informationBalancing. Rotating Components Examples of rotating components in a mechanism or a machine. (a)
alacig NOT COMPLETE Rotatig Compoets Examples of rotatig compoets i a mechaism or a machie. Figure 1: Examples of rotatig compoets: camshaft; crakshaft Sigle-Plae (Static) alace Cosider a rotatig shaft
More informationThe Poisson Process *
OpeStax-CNX module: m11255 1 The Poisso Process * Do Johso This work is produced by OpeStax-CNX ad licesed uder the Creative Commos Attributio Licese 1.0 Some sigals have o waveform. Cosider the measuremet
More information17 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 informationSCIENCE CHINA Technological Sciences. Transverse vibration of pre-tensioned nonlocal nanobeams with precise internal axial loads
SCIENCE CHINA Techological Scieces RESEARCH PAPER August Vol5 No8: 7 3 doi: 7/s3--79-9 Trasverse vibratio of pre-tesioed olocal aobeams with precise iteral axial loads I Cheg,,3, IM C W,3*, YU Jii,3 &
More informationLecture 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 informationFAILURE 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 informationModeling of Finite Gradient Elasticity Using Optimal Transportation Method
9 th Iteratioal Cogress of Croatia Society of Mechaics 18-22 September 2018 Split, Croatia Modelig of Fiite Gradiet Elasticity Usig Optimal Trasportatio Method Boris JALUŠIĆ *, Christia WEIßENFELS +, Jurica
More informationComparison of Minimum Initial Capital with Investment and Non-investment Discrete Time Surplus Processes
The 22 d Aual Meetig i Mathematics (AMM 207) Departmet of Mathematics, Faculty of Sciece Chiag Mai Uiversity, Chiag Mai, Thailad Compariso of Miimum Iitial Capital with Ivestmet ad -ivestmet Discrete Time
More informationNumerical 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 informationSimilarity Solutions to Unsteady Pseudoplastic. Flow Near a Moving Wall
Iteratioal Mathematical Forum, Vol. 9, 04, o. 3, 465-475 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/0.988/imf.04.48 Similarity Solutios to Usteady Pseudoplastic Flow Near a Movig Wall W. Robi Egieerig
More informationA theoretical and computational framework for studying creep crack growth
It J Fract DOI.7/s74-7-23-2 IUTAM BALTIMORE A theoretical ad computatioal framework for studyig creep crack growth Elsiddig Elmukashfi Ala C. F. Cocks Received: 2 March 27 / Accepted: 29 Jue 27 The Author(s)
More informationSimilarity 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 informationTHE FRACTAL TECHNIQUES APPLIED IN NUMERICAL ITERATIVE METHODS FRAME
6 th INTERNATIONAL MULTIDISCIPLINARY CONFERENCE THE FRACTAL TECHNIQUES APPLIED IN NUMERICAL ITERATIVE METHODS FRAME George MAHALU, Radu PENTIUC Uiversity of Suceava - mahalu@eed.usv.ro Abstract: May methods
More information(c) Write, but do not evaluate, an integral expression for the volume of the solid generated when R is
Calculus BC Fial Review Name: Revised 7 EXAM Date: Tuesday, May 9 Remiders:. Put ew batteries i your calculator. Make sure your calculator is i RADIAN mode.. Get a good ight s sleep. Eat breakfast. Brig:
More informationMathematical 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 informationNotes on iteration and Newton s method. Iteration
Notes o iteratio ad Newto s method Iteratio Iteratio meas doig somethig over ad over. I our cotet, a iteratio is a sequece of umbers, vectors, fuctios, etc. geerated by a iteratio rule of the type 1 f
More informationA collocation method for singular integral equations with cosecant kernel via Semi-trigonometric interpolation
Iteratioal Joural of Mathematics Research. ISSN 0976-5840 Volume 9 Number 1 (017) pp. 45-51 Iteratioal Research Publicatio House http://www.irphouse.com A collocatio method for sigular itegral equatios
More informationThe 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 informationTHE 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 informationModern Applied Science
www.ccseet.org/joural.html Vol., No. 4 November 7 Performace Aalysis o Equivalet Elasticity of D 4-directioal Braided omposites Dog he, Li he, Yig Su Key Laboratory of Advaced Textile omposite Material
More informationPoincaré Problem for Nonlinear Elliptic Equations of Second Order in Unbounded Domains
Advaces i Pure Mathematics 23 3 72-77 http://dxdoiorg/4236/apm233a24 Published Olie Jauary 23 (http://wwwscirporg/oural/apm) Poicaré Problem for Noliear Elliptic Equatios of Secod Order i Ubouded Domais
More informationA COMPLEX VARIABLE BOUNDARY ELEMENT METHOD FOR AN ELLIPTIC PARTIAL DIFFERENTIAL EQUATION WITH VARIABLE COEFFICIENTS
A COMPLEX VARIABLE BOUNDARY ELEMENT METHOD FOR AN ELLIPTIC PARTIAL DIFFERENTIAL EQUATION WITH VARIABLE COEFFICIENTS Y. S. PARK a adw.t.ang b Computatioal ad Mathematical Scieces Group Faculty of Iformatio
More informationHigher-order iterative methods by using Householder's method for solving certain nonlinear equations
Math Sci Lett, No, 7- ( 7 Mathematical Sciece Letters A Iteratioal Joural http://dxdoiorg/785/msl/5 Higher-order iterative methods by usig Householder's method for solvig certai oliear equatios Waseem
More informationIIT JAM Mathematical Statistics (MS) 2006 SECTION A
IIT JAM Mathematical Statistics (MS) 6 SECTION A. If a > for ad lim a / L >, the which of the followig series is ot coverget? (a) (b) (c) (d) (d) = = a = a = a a + / a lim a a / + = lim a / a / + = lim
More informationAnalysis of Dual Boundary Element in Shape Optimization
IOSR Joural of Mechaical ad Civil Egieerig (IOSR-JMCE) e-issn: 78-1684,p-ISSN: 30-334X, Volume 13, Issue 3 Ver. IV (May- Ju. 016), PP 17-141 www.iosrjourals.org Aalysis of Dual Boudary Elemet i Shape Optimizatio
More informationAnalytic Continuation
Aalytic Cotiuatio The stadard example of this is give by Example Let h (z) = 1 + z + z 2 + z 3 +... kow to coverge oly for z < 1. I fact h (z) = 1/ (1 z) for such z. Yet H (z) = 1/ (1 z) is defied for
More informationINF-GEO Solutions, Geometrical Optics, Part 1
INF-GEO430 20 Solutios, Geometrical Optics, Part Reflectio by a symmetric triagular prism Let be the agle betwee the two faces of a symmetric triagular prism. Let the edge A where the two faces meet be
More informationFailure 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 informationSequences, Mathematical Induction, and Recursion. CSE 2353 Discrete Computational Structures Spring 2018
CSE 353 Discrete Computatioal Structures Sprig 08 Sequeces, Mathematical Iductio, ad Recursio (Chapter 5, Epp) Note: some course slides adopted from publisher-provided material Overview May mathematical
More information