Three-Dimensional Explicit Parallel Finite Element Analysis of Functionally Graded Solids under Impact Loading. Ganesh Anandakumar and Jeong-Ho Kim
|
|
- Ralph Phelps
- 5 years ago
- Views:
Transcription
1 Three-Dimensional Eplicit Parallel Finite Element Analsis of Functionall Graded Solids under Impact Loading Ganesh Anandaumar and Jeong-Ho Kim Department of Civil and Environmental Engineering, Universit of Connecticut, 6 Glenbroo Rd, U-7, Storrs, CT 669, U.S.A. Januar, 9 Abstract This paper presents two numerical eamples to investigate the behavior of three-dimensional (D) functionall graded (FG) solids under dnamic loading using eplicit parallel finite element method. In the first eample, wave propagation in a D FG bar under transient sinusoidal loading is investigated. Material gradation and thicness (D) effects are seen in the dnamic stress behavior of the FG bar. In the second eample, a three-point bending beam made of epo and glass phases under velocit impact is studied. Bending stress histor for beam with higher values of material properties at the loading edge is consistentl higher than that of the homogeneous beam and the beam with lower values of material properties at the loading edge. Larger bending stresses for the foremost beam ma indicate earlier crac initiation time than the other two beams which was proven b eperiments performed b other researchers. Kewords: functionall graded material (FGM), dnamic analsis, three-dimensional wave propagation, finite element methods (FEM), parallel computing. Introduction Functionall graded material (FGM) is a material solution & concept used for a new advanced class of composites. is characteried b a gradual variation in composition, microstructure and material properties. eperiences severe dnamic thermal and mechanical loadings. tpicall fails through cracing and spallation. therefore requires a detailed dnamic analsis to determine the effect of material gradation and obtain locations of pea values of stresses. Several numerical methods have been used to investigate FGMs, including integral equations (Otur and Erdogan, 997), boundar element methods (Sutradhar et al., ), finite element methods (Santare and Lambros, ), etc. In this stud, we used the displacement-based FEM and graded finite elements (Kim and Paulino, ) to model FGMs using the direct Gaussian integration formulation. Parallel eplicit FEM is used for obtaining the dnamic response of D FG solids using Message Passing Interface (MPI) standard (MPI, 997). The following are the novelties of this stud: Wave propagation analsis of a FG D bar using FEM Dnamic analsis of a FG D beam under velocit impact load using FEM Parallel Eplicit Dnamic FEA using Newmar-β method The steps involved in the eplicit Newmar-β (γ=.5 and β=) method (Newmar, 959) are given below. Corresponding Author, anandg@engr.uconn.edu
2 calculate velocities at time t n : u n = (ü un n +.5 t n + ü n ) calculate displacements at time t n+ : u n+ = u n + t n u n +.5 t nü n compute effective force vector (EFV) at time t n+ : solve for accelerations at time t n+ : ü n+ = M fn+ n+ f = fet n+ fint n+ where u, u, and ü represent displacement, velocit, and acceleration vectors, respectivel, M is the lumped mass matri, fet is the eternal force vector, n is the time step number, fint(=ku) is the internal force vector, and denotes degree of freedom (DOF). Figure shows a flowchart of the steps involved in the parallel eecution of the dnamic FE code using a master-slave approach. Tpicall, the echange of EFV needs to be done onl at processors that share the node through which partition is done. But this leads to man small messages being sent from one processor to another and ma lead to increase in communication time which is undesirable. To overcome this, the EFV at the slaves are sent to the master for assembl and the master returns the assembled EFV bac to the slaves for calculating acceleration vector. Initialie number of processors using MPI_INIT and MPI_COMM_RANK Slaves send EFV to master for assembl using MPI_SEND and MPI_RECV Partition FE mesh in to sub meshes manuall Calculate K & M matrices for local elements Master returns assembled force vector to slaves. Calculate acceleration vector. Loop over number of timesteps for time integration using Newmar beta method Repeat time integration for the rest of the simulation b this echange algorithm Calculate displacement, velocit, and effective force vector (EFV) at slaves & master Eit parallel code using MPI_FINALIZE Numerical Eamples Figure : Flow chart of the parallel eecution of the eplicit FE code. Eample : Wave Propagation in a D FG Bar: Wave propagation in a fied-free bar with graded materials in the direction is simulated to obtain the effect of material gradation. Consider a fied-free square bar (Figure (a)) of length L =m and height H =.5 m under a transient sinusoidal load (Figure (b)) applied at the free end of the bar. W=.5 m H=.5 m L =. m (a) P f(t) f(t) t(µs) (b) Figure : Eample : Wave propagation in FG D bar (a) schematic; (b) sinusoidal load. The material properties var linearl from alumina in the bottom to steel at the top. The material properties of steel and alumina are shown in Table. Due to fine mesh criterion, the D FE mesh is discretied into 5 5 quads, each quad divided into four 5-node wedge elements, totaling 7, elements and 75,86 nodes.
3 Table : Eample : Material properties of steel and alumina considered Material E(GP a) ν ρ(g/m ) C d (m/s) = =.5L =L =.5L =L Midplane σ /σ L σ /σ L Midplane H = L =.5L = L 5 6 t*(c d ) /L 5 6 t*(c d ) /L (a) (b) Figure : Eample : Normalied longitudinal stress histor of 9 points (see the insert) on a graded bar subjected to transient sinusoidal loading at (a) =.5 m and (b) =m. Solid, dashed, and dash-dot lines indicate points at =, =.5L, and = L, respectivel. Thic, intermediate-thic and thin lines indicate -rich side, mid-plane and -rich side points, respectivel. Figure shows the stress histor at 9 locations for the FG bar at (a) =.5 m and (b) =m obtained using the parallel FE code. Despite the gradation, the stress wave remains the same for the three locations at = L, as it is the prescribed initial condition. At other locations ( =.5L and = ), we see that the stress wave gets distorted in time and the magnitude differs considerabl along the direction due to differences in material properties and wave speeds. The alumina side at different locations undergoes higher stresses when compared with the steel side, more so at the fied end than at other locations. It is interesting to see that the stresses at the fied end = is considerabl higher at =m than at =.5m. This is because of the D discretiation and the fied boundar condition at that location of the bar. This outcome is a novelt of this stud. Eample : An FG D Beam Subjected to Velocit Impact: A three point bending beam (TPBB) under velocit impact is studied to understand the influence of material gradation on the beam bending behavior. The beam, which is a real FGM sstem, is made of glass/epo phases. The dnamic fracture eperiments of the linearl graded specimens has been conducted b Rousseau and Tippur (). This stud offers nowledge of the dnamic behavior of this material sstem and understanding of the stress field which helps to predict fracture initiation times in various graded specimens. Consider a TPBB under velocit impact load of m/s applied at the top as shown in Figure (a). Due to the smmetr of the geometr and the loading conditions, onl one-fourth of the beam is modeled for numerical analsis (Figure (b)). The D FE mesh (Figure (c)) has an uniform element sie of 9.5 µm along the loaded edge. Point P (,.W ) is of significance because it corresponds to the location of the crac tip in the dnamic fracture analsis of the beam (Rousseau and Tippur, ) and is also the location at which stress results are obtained in this stud. Three material gradation cases are considered for the dnamic analsis of the TPBB:
4 V = m/s V = m/s V = m/s E E E W = 7 mm W = 7 mm W = 7 mm E L = 5 mm B = 6 mm P.W L = 76 mm E B = mm P.W L = 76 mm (a) (b) (c) (d) E B = mm Figure : Eample : Epo/glass beam subjected to velocit impact; (a) geometr and boundar conditions; (b) line load case; (b) point load case; (d) D FE mesh of one-quarter model with node wedge elements and 87 nodes. Homogeneous beam (Homog, E = E ) Beam stiffer at the impacted surface (StiffT op, E >E ) Beam softer at the impacted surface (StiffBot, E <E ) where subscripts and denote bottom and top surfaces of the beam, respectivel. The material properties for the beam are obtained from Zhang and Paulino (7). Two tpes of smmetric loadings are considered in this stud. The first (Figure (b)) being the load being applied throughout the thicness of the beam (line load) and the second (point load) where the load is applied onl at a central node (Figure (c)). The point load case can onl be analed using D finite elements and is a novelt of this stud. 5 point load σ (MPa) 5 σ (MPa) line load 5 line load point load t*c /W d(avg) t*c /W d(avg) Figure 5: Eample : Stress histories σ and σ at location P(,.W ) for homogeneous and graded beams subjected to impact velocit of m/s as a line load (thic lines) and point load (thin lines). Solid, dashed, and dash-dot lines correspond to Stif f T op, Homog, and Stif f Bot beams, respectivel. Figure 5 shows the comparison of σ and σ at point P for the three beams under impact velocit of V = m/s applied as line load and as point load. We see that the stresses for the latter case is much lower compared to the former. Maimum tensile stress (σ ) is eperienced b the StiffT op beam followed b the Homog beam and the StiffBot beam, for both loading scenarios. This ma indicate earlier crac initiation time for StiffT op beam when compared to other two beams which was proven b eperiments performed b Rousseau and Tippur (). Figure 6 shows the stress σ contour at t = tc d(avg) /W 6 for the three graded beams and two load cases. Since the stress magnitudes for the two load cases are different, we use different ranges for the contour plot. There is quite a different behavior of stress waves in the top quarter region when comparing the si cases.
5 (a) (b) (c) (d) (e) (f) Figure 6: Eample : Stress contour (σ, MPa) along the direction at t 6. StiffBot: (a) line load; (d) point load; Homog: (b) line load; (e) point load; StiffT op: (c) line load; (f) point load. Concluding Remars Dnamic behavior of D FG solids is investigated using eplicit parallel FE formulation. Thicness (D) and gradation effects are seen in the stress behavior of the FG bar. Material gradation considerabl affects the dnamic stress behavior of the beam. Tensile stress is maimum for StiffT op beam at the imaginar crac-tip location indicating that crac initiation will occur earlier for this beam which was verified b Rousseau and Tippur (). References [] J.-H. Kim and G. H. Paulino. Isoparametric graded finite elements for nonhomogeneous isotropic and orthotropic materials. Journal of Applied Mechanics, Transactions ASME, 69():5 5,. [] MPI-:. Etension to the Message Passing Interface. Universit of Tennessee, Knoville, Tennessee, 997. [] N. M. Newmar. A method of computation for structural dnamics. Journal of Engineering Mechanics (ASCE), 85:67 9, 959. [] M. Otur and F. Erdogan. Mode I crac problem in an inhomogeneous orthotropic medium. International Journal of Engineering Science, 5(9):869 88, 997. [5] C.-E. Rousseau and H. V. Tippur. Dnamic fracture of compositionall graded materials with cracs along the elastic gradient: eperiments and analsis. Mechanics of Materials, :, a. [6] M. H. Santare and J. Lambros. Use of graded finite elements to model the behavior of nonhomogeneous materials. Journal of Applied Mechanics, Transactions ASME, 67():89 8,. [7] S. Sutradhar, G. H. Paulino, and L. J. Gra. Transient heat conduction in homogeneous and nonhomogeneous materials b the Laplace transform galerin boundar element method. Engineering Analsis with Boundar Elements, 6():9,. [8] Z. Zhang and G. H. Paulino. Wave propagation and dnamic analsis of smoothl graded heterogeneous continua using graded finite elements. International Journal of Solids and Structures, :6 66, 7. 5
Chapter 11 Three-Dimensional Stress Analysis. Chapter 11 Three-Dimensional Stress Analysis
CIVL 7/87 Chapter - /39 Chapter Learning Objectives To introduce concepts of three-dimensional stress and strain. To develop the tetrahedral solid-element stiffness matri. To describe how bod and surface
More informationMMJ1153 COMPUTATIONAL METHOD IN SOLID MECHANICS PRELIMINARIES TO FEM
B Course Content: A INTRODUCTION AND OVERVIEW Numerical method and Computer-Aided Engineering; Phsical problems; Mathematical models; Finite element method;. B Elements and nodes, natural coordinates,
More informationEshan V. Dave, Secretary of M&FGM2006 (Hawaii) Research Assistant and Ph.D. Candidate. Glaucio H. Paulino, Chairman of M&FGM2006 (Hawaii)
Asphalt Pavement Aging and Temperature Dependent Properties through a Functionally Graded Viscoelastic Model I: Development, Implementation and Verification Eshan V. Dave, Secretary of M&FGM2006 (Hawaii)
More informationThe Plane Stress Problem
. 4 The Plane Stress Problem 4 Chapter 4: THE PLANE STRESS PROBLEM 4 TABLE OF CONTENTS Page 4.. INTRODUCTION 4 3 4... Plate in Plane Stress............... 4 3 4... Mathematical Model.............. 4 4
More informationMixed-Mode Crack Propagation in Functionally Graded Materials
Materials Science Forum Vols. 492-493 (25) pp. 49-414 online at http://www.scientific.net 25 Trans Tech Publications, Switzerland Mixed-Mode Crack Propagation in Functionally Graded Materials Jeong-Ho
More informationChapter 6 2D Elements Plate Elements
Institute of Structural Engineering Page 1 Chapter 6 2D Elements Plate Elements Method of Finite Elements I Institute of Structural Engineering Page 2 Continuum Elements Plane Stress Plane Strain Toda
More informationThe Plane Stress Problem
. 14 The Plane Stress Problem 14 1 Chapter 14: THE PLANE STRESS PROBLEM 14 TABLE OF CONTENTS Page 14.1. Introduction 14 3 14.1.1. Plate in Plane Stress............... 14 3 14.1.. Mathematical Model..............
More informationChapter 2 Overview of the Simulation Methodology
Chapter Overview of the Simulation Methodolog. Introduction This chapter presents an overview of the simulation methodolog that comprises both the art and science involved in simulating phsical phenomena.
More informationThe Plane Stress Problem
14 The Plane Stress Problem IFEM Ch 14 Slide 1 Plate in Plane Stress Thickness dimension or transverse dimension z Top surface Inplane dimensions: in, plane IFEM Ch 14 Slide 2 Mathematical Idealization
More informationDIC to Study Deformations. * Jajam, Tippur, Int. J. Solids/Structures, Strains from DIC
DIC to Stud Deformations * Jajam, Tippur, Int. J. Solids/Structures, Strains from DIC Sub image Size Effect a) b) c) d) Contours of displacement from piel sub-images (a, b) indicate nois data, while piel
More informationINTERFACE CRACK IN ORTHOTROPIC KIRCHHOFF PLATES
Gépészet Budapest 4-5.Ma. G--Section-o ITERFACE CRACK I ORTHOTROPIC KIRCHHOFF PLATES András Szekrénes Budapest Universit of Technolog and Economics Department of Applied Mechanics Budapest Műegetem rkp.
More informationResearch Article Equivalent Elastic Modulus of Asymmetrical Honeycomb
International Scholarl Research Network ISRN Mechanical Engineering Volume, Article ID 57, pages doi:.5//57 Research Article Equivalent Elastic Modulus of Asmmetrical Honecomb Dai-Heng Chen and Kenichi
More informationKINEMATIC RELATIONS IN DEFORMATION OF SOLIDS
Chapter 8 KINEMATIC RELATIONS IN DEFORMATION OF SOLIDS Figure 8.1: 195 196 CHAPTER 8. KINEMATIC RELATIONS IN DEFORMATION OF SOLIDS 8.1 Motivation In Chapter 3, the conservation of linear momentum for a
More informationStability Analysis of a Geometrically Imperfect Structure using a Random Field Model
Stabilit Analsis of a Geometricall Imperfect Structure using a Random Field Model JAN VALEŠ, ZDENĚK KALA Department of Structural Mechanics Brno Universit of Technolog, Facult of Civil Engineering Veveří
More informationBOUNDARY EFFECTS IN STEEL MOMENT CONNECTIONS
BOUNDARY EFFECTS IN STEEL MOMENT CONNECTIONS Koung-Heog LEE 1, Subhash C GOEL 2 And Bozidar STOJADINOVIC 3 SUMMARY Full restrained beam-to-column connections in steel moment resisting frames have been
More informationWilliam G. Buttlar Glaucio H. Paulino Harry H. Hilton Phillip B. Blankenship Hervé DiBenedetto. All colleagues and friends
Ph.D. Final Exam Asphalt Pavement Aging and Temperature Dependent Properties using a Functionally Graded Viscoelastic Model Eshan V. Dave June 5th 29 Department of Civil and Environmental Engineering University
More informationTHE GENERAL ELASTICITY PROBLEM IN SOLIDS
Chapter 10 TH GNRAL LASTICITY PROBLM IN SOLIDS In Chapters 3-5 and 8-9, we have developed equilibrium, kinematic and constitutive equations for a general three-dimensional elastic deformable solid bod.
More informationEVALUATION OF STRESS IN BMI-CARBON FIBER LAMINATE TO DETERMINE THE ONSET OF MICROCRACKING
EVALUATION OF STRESS IN BMI-CARBON FIBER LAMINATE TO DETERMINE THE ONSET OF MICROCRACKING A Thesis b BRENT DURRELL PICKLE Submitted to the Office of Graduate Studies of Teas A&M Universit in partial fulfillment
More informationMAE 323: Chapter 4. Plane Stress and Plane Strain. The Stress Equilibrium Equation
The Stress Equilibrium Equation As we mentioned in Chapter 2, using the Galerkin formulation and a choice of shape functions, we can derive a discretized form of most differential equations. In Structural
More informationINTERNATIONAL JOURNAL OF APPLIED ENGINEERING RESEARCH, DINDIGUL Volume 2, No 1, 2011
Interlaminar failure analysis of FRP cross ply laminate with elliptical cutout Venkateswara Rao.S 1, Sd. Abdul Kalam 1, Srilakshmi.S 1, Bala Krishna Murthy.V 2 1 Mechanical Engineering Department, P. V.
More informationShear Stress and Interlaminar Shear Strength Tests of Cross-laminated Timber Beams
Shear Stress and nterlaminar Shear Strength Tests of Cross-laminated Timber Beams Yao Lu, Wenbo Xie, Zheng Wang,* and Zizhen Gao The interlaminar shear stresses of the three-laer, five-laer, and sevenlaer
More informationConsistent Formulations of the Interaction Integral Method for Fracture of Functionally Graded Materials
Jeong-Ho Kim Glaucio H. Paulino 2 e-mail: paulino@uiuc.edu Department of Civil and Environmental Engineering, Newmark Laboratory, The University of Illinois at Urbana-Champaign, 205 North Mathews Avenue,
More informationStability Analysis of Laminated Composite Thin-Walled Beam Structures
Paper 224 Civil-Comp Press, 2012 Proceedings of the Eleventh International Conference on Computational Structures echnolog, B.H.V. opping, (Editor), Civil-Comp Press, Stirlingshire, Scotland Stabilit nalsis
More informationOutline. Organization. Stresses in Beams
Stresses in Beams B the end of this lesson, ou should be able to: Calculate the maimum stress in a beam undergoing a bending moment 1 Outline Curvature Normal Strain Normal Stress Neutral is Moment of
More informationMECHANICS OF MATERIALS REVIEW
MCHANICS OF MATRIALS RVIW Notation: - normal stress (psi or Pa) - shear stress (psi or Pa) - normal strain (in/in or m/m) - shearing strain (in/in or m/m) I - area moment of inertia (in 4 or m 4 ) J -
More informationVibrational Power Flow Considerations Arising From Multi-Dimensional Isolators. Abstract
Vibrational Power Flow Considerations Arising From Multi-Dimensional Isolators Rajendra Singh and Seungbo Kim The Ohio State Universit Columbus, OH 4321-117, USA Abstract Much of the vibration isolation
More informationGeneric Strategies to Implement Material Grading in Finite Element Methods for Isotropic and Anisotropic Materials
University of Nebraska - Lincoln DigitalCommons@University of Nebraska - Lincoln Engineering Mechanics Dissertations & Theses Mechanical & Materials Engineering, Department of Winter 12-9-2011 Generic
More information6.1 The Linear Elastic Model
Linear lasticit The simplest constitutive law or solid materials is the linear elastic law, which assumes a linear relationship between stress and engineering strain. This assumption turns out to be an
More informationFracture Test & Fracture Parameters of Self Compacting Concrete using ANSYS. Zeel Vashi 1,Megha Thomas 2 I. INTRODUCTION
International Journal of Technical Innovation in Modern Engineering & Science (IJTIMES) Impact Factor: 3.45 (SJIF-2015), e-issn: 2455-2584 Volume 3, Issue 05, May-2017 Fracture Test & Fracture Parameters
More informationStructural Analysis of an Exhaust System for Heavy Trucks
Structural Analsis of an Ehaust Sstem for Heav Trucks Mr N. Vasconcellos, Mr F. dos Anjos and Mr M. Argentino, debis humaitá IT Services Latin America, Brail ABSTRACT The present work describes the use
More informationBending of Shear Deformable Plates Resting on Winkler Foundations According to Trigonometric Plate Theory
J. Appl. Comput. ech., 4(3) (018) 187-01 DOI: 10.055/JAC.017.3057.1148 ISSN: 383-4536 jacm.scu.ac.ir Bending of Shear Deformable Plates Resting on Winkler Foundations According to Trigonometric Plate Theor
More informationME 475 Modal Analysis of a Tapered Beam
ME 475 Modal Analysis of a Tapered Beam Objectives: 1. To find the natural frequencies and mode shapes of a tapered beam using FEA.. To compare the FE solution to analytical solutions of the vibratory
More informationSurvey of Wave Types and Characteristics
Seminar: Vibrations and Structure-Borne Sound in Civil Engineering Theor and Applications Surve of Wave Tpes and Characteristics Xiuu Gao April 1 st, 2006 Abstract Mechanical waves are waves which propagate
More informationModule 2 Stresses in machine elements. Version 2 ME, IIT Kharagpur
Module Stresses in machine elements Version M, IIT Kharagpur Lesson 3 Strain analsis Version M, IIT Kharagpur Instructional Objectives At the end of this lesson, the student should learn Normal and shear
More informationFig. 1. Different locus of failure and crack trajectories observed in mode I testing of adhesively bonded double cantilever beam (DCB) specimens.
a). Cohesive Failure b). Interfacial Failure c). Oscillatory Failure d). Alternating Failure Fig. 1. Different locus of failure and crack trajectories observed in mode I testing of adhesively bonded double
More informationLinear plate bending and laminate theory
Linear plate bending and laminate theor 4M: Design Tools Eindhoven Universit of Technolog Linear plate bending A plate is a bod of which the material is located in a small region around a surface in the
More informationFATIGUE ANALYSIS: THE SUPER-NEUBER TECHNIQUE FOR CORRECTION OF LINEAR ELASTIC FE RESULTS
6 TH INTERNATIONAL CONGRESS OF THE AERONAUTICAL SCIENCES FATIGUE ANALYSIS: THE SUPER-NEUBER TECHNIQUE FOR CORRECTION OF LINEAR ELASTIC FE RESULTS L. Samuelsson * * Volvo Aero Corporation, Sweden Kewords:
More informationTheoretical Manual Theoretical background to the Strand7 finite element analysis system
Theoretical Manual Theoretical background to the Strand7 finite element analysis system Edition 1 January 2005 Strand7 Release 2.3 2004-2005 Strand7 Pty Limited All rights reserved Contents Preface Chapter
More informationChapter 3: BASIC ELEMENTS. solid mechanics)
Chapter 3: BASIC ELEMENTS Section 3.: Preliminaries (review of solid mechanics) Outline Most structural analsis FE codes are displacement based In this chapter we discuss interpolation methods and elements
More informationData Repository Hampel et al., page 1/5
GSA DATA REPOSITORY 2138 Data Repositor Hampel et al., page 1/5 SETUP OF THE FINITE-ELEMENT MODEL The finite-element models were created with the software ABAQUS and consist of a 1-km-thick lithosphere,
More informationFinite elements for plates and shells. Advanced Design for Mechanical System LEC 2008/11/04
Finite elements for plates and shells Sstem LEC 008/11/04 1 Plates and shells A Midplane The mid-surface of a plate is plane, a curved geometr make it a shell The thickness is constant. Ma be isotropic,
More informationCOMPUTATIONAL MODELING AND SIMULATION OF LOW VELOCITY IMPACT ON FIBROUS COMPOSITE PANELS DROP_WEIGHT UN_PARTITIONED MODEL
VOL. 4, O., PRIL 9 ISS 89-8 RP Journal of Engineering and pplied Sciences 6-9 sian Research Publishing etwor (RP). ll rights reserved. COMPUTTIOL MOELIG SIMULTIO OF LOW VELOCITY IMPCT O FIROUS COMPOSITE
More informationA PAPER ON DESIGN AND ANALYSIS OF PRESSURE VESSEL
A PAPER ON DESIGN AND ANALYSIS OF PRESSURE VESSEL P.Palanivelu 1, R.Siva Prasad 2, 1 PG Scholar, Department of Mechanical Engineering, Gojan School of Business and Technology, Redhills, Chennai, India.
More informationThree-dimensional free vibration analysis of functionally graded rectangular plates using the B-spline Ritz method
The 01 World Congress on Advances in Civil, Environmental, and Materials Research (ACEM 1) Seoul, Korea, August 6-30, 01 Three-dimensional free viration analsis of functionall graded rectangular plates
More informationModeling of Interfacial Debonding Induced by IC Crack for Concrete Beam-bonded with CFRP
Proceedings of the World Congress on Engineering 21 Vol II WCE 21, June 2 - July 1, 21, London, U.K. Modeling of Interfacial Debonding Induced by IC Crack for Concrete Beam-bonded with CFRP Lihua Huang,
More informationMulti-Point Constraints
Multi-Point Constraints Multi-Point Constraints Multi-Point Constraints Single point constraint examples Multi-Point constraint examples linear, homogeneous linear, non-homogeneous linear, homogeneous
More informationModeling of Failure along Predefined Planes in Fractured Reservoirs
PROCEEDINGS, Thirt-Ninth Workshop on Geothermal Reservoir Engineering Stanford Universit, Stanford, California, Februar 24-26, 2014 SGP-TR-202 Modeling of Failure along Predefined Planes in Fractured Reservoirs
More informationSecond-Order Linear Differential Equations C 2
C8 APPENDIX C Additional Topics in Differential Equations APPENDIX C. Second-Order Homogeneous Linear Equations Second-Order Linear Differential Equations Higher-Order Linear Differential Equations Application
More informationRESEARCH ON PIEZOELECTRIC QUARTZ UNDER MULTIDIMENSIONAL FORCES
aterials Phsics and echanics (015) 94-100 eceived: October 31, 014 ESEACH ON PIEZOELECTIC QUATZ UNDE ULTIDIENSIONAL OCES ei Wang *, uji Wang, Zongjin en, Zhenuan Jia, Liqi Zhao, Dong Li Ke Laborator for
More informationSupplementary Figures
Fracture Strength (GPa) Supplementary Figures a b 10 R=0.88 mm 1 0.1 Gordon et al Zhu et al Tang et al im et al 5 7 6 4 This work 5 50 500 Si Nanowire Diameter (nm) Supplementary Figure 1: (a) TEM image
More informationMODELLING BAMBOO AS A FUNCTIONALLY GRADED MATERIAL
MODELLING BAMBOO AS A FUNCTIONALLY GRADED MATERIAL Emílio Carlos Nelli Silva Associate Professor Department of ofmechatronics and Mechanical Systems Engineering Escola Politécnica da dauniversidade de
More informationx y plane is the plane in which the stresses act, yy xy xy Figure 3.5.1: non-zero stress components acting in the x y plane
3.5 Plane Stress This section is concerned with a special two-dimensional state of stress called plane stress. It is important for two reasons: () it arises in real components (particularl in thin components
More informationCH.7. PLANE LINEAR ELASTICITY. Multimedia Course on Continuum Mechanics
CH.7. PLANE LINEAR ELASTICITY Multimedia Course on Continuum Mechanics Overview Plane Linear Elasticit Theor Plane Stress Simplifing Hpothesis Strain Field Constitutive Equation Displacement Field The
More informationChapter 14 Truss Analysis Using the Stiffness Method
Chapter 14 Truss Analsis Using the Stiffness Method Structural Mechanics 2 ept of Arch Eng, Ajou Univ Outline undamentals of the stiffness method Member stiffness matri isplacement and force transformation
More informationA Software Tool for Mechanics of Composite Materials*
Int. J. Engng Ed. Vol. 13, No., p. 433±441, 17 04-14X/1 $3.00+0.00 Printed in Great Britain. # 17 TEMPUS Publications. A Software Tool for Mechanics of Composite Materials* AUTAR K. KAW and GARY WILLENBRING
More informationDEVELOPMENT OF TEST GUIDANCE FOR COMPACT TENSION FRACTURE TOUGHNESS SPECIMENS CONTAINING NOTCHES INSTEAD OF FATIGUE PRE-CRACKS
Transactions, SMiRT-23 Division II, Paper ID 287 Fracture Mechanics and Structural Integrity DEVELOPMENT OF TEST GUIDANCE FOR COMPACT TENSION FRACTURE TOUGHNESS SPECIMENS CONTAINING NOTCHES INSTEAD OF
More informationApplications of Gauss-Radau and Gauss-Lobatto Numerical Integrations Over a Four Node Quadrilateral Finite Element
Avaiable online at www.banglaol.info angladesh J. Sci. Ind. Res. (), 77-86, 008 ANGLADESH JOURNAL OF SCIENTIFIC AND INDUSTRIAL RESEARCH CSIR E-mail: bsir07gmail.com Abstract Applications of Gauss-Radau
More informationSimulation of Acoustic and Vibro-Acoustic Problems in LS-DYNA using Boundary Element Method
10 th International LS-DYNA Users Conference Simulation Technolog (2) Simulation of Acoustic and Vibro-Acoustic Problems in LS-DYNA using Boundar Element Method Yun Huang Livermore Software Technolog Corporation
More informationStress and Strain ( , 3.14) MAE 316 Strength of Mechanical Components NC State University Department of Mechanical & Aerospace Engineering
(3.8-3.1, 3.14) MAE 316 Strength of Mechanical Components NC State Universit Department of Mechanical & Aerospace Engineering 1 Introduction MAE 316 is a continuation of MAE 314 (solid mechanics) Review
More information3. Overview of MSC/NASTRAN
3. Overview of MSC/NASTRAN MSC/NASTRAN is a general purpose finite element analysis program used in the field of static, dynamic, nonlinear, thermal, and optimization and is a FORTRAN program containing
More informationCH. 1 FUNDAMENTAL PRINCIPLES OF MECHANICS
446.201 (Solid echanics) Professor Youn, eng Dong CH. 1 FUNDENTL PRINCIPLES OF ECHNICS Ch. 1 Fundamental Principles of echanics 1 / 14 446.201 (Solid echanics) Professor Youn, eng Dong 1.2 Generalied Procedure
More informationREVIEW FOR EXAM II. Dr. Ibrahim A. Assakkaf SPRING 2002
REVIEW FOR EXM II. J. Clark School of Engineering Department of Civil and Environmental Engineering b Dr. Ibrahim. ssakkaf SPRING 00 ENES 0 Mechanics of Materials Department of Civil and Environmental
More informationBASE VECTORS FOR SOLVING OF PARTIAL DIFFERENTIAL EQUATIONS
BASE VECTORS FOR SOLVING OF PARTIAL DIFFERENTIAL EQUATIONS J. Roubal, V. Havlena Department of Control Engineering, Facult of Electrical Engineering, Czech Technical Universit in Prague Abstract The distributed
More information445. Investigation of Thermo-Elastic Damping of Vibrations of Rectangular and Ring-Shaped MEMS Resonators
. Investigation of Thermo-Elastic Damping of Vibrations of Rectangular and Ring-Shaped MEMS Resonators R. Barauskas,a, S. Kausinis,b, H. A. C. Tilmans Kaunas Universit of Technolog, Studentu -7, LT-68,
More informationDISPENSA FEM in MSC. Nastran
DISPENSA FEM in MSC. Nastran preprocessing: mesh generation material definitions definition of loads and boundary conditions solving: solving the (linear) set of equations components postprocessing: visualisation
More informationEffect of Specimen Dimensions on Flexural Modulus in a 3-Point Bending Test
Effect of Specimen Dimensions on Flexural Modulus in a 3-Point Bending Test M. Praveen Kumar 1 and V. Balakrishna Murthy 2* 1 Mechanical Engineering Department, P.V.P. Siddhartha Institute of Technology,
More informationMECHANICS OF MATERIALS
00 The McGraw-Hill Companies, Inc. All rights reserved. T Edition CHAPTER MECHANICS OF MATERIALS Ferdinand P. Beer E. Russell Johnston, Jr. John T. DeWolf Lecture Notes: J. Walt Oler Teas Tech Universit
More informationFinite element simulation of residual stresses in laser heating
IAS-2008-66-546ST Finite element simulation of residual stresses in laser heating G. H. Farrahi 1, M. Sistaninia 2, H. Moeinoddini 3 1,2-School of Mechanical Engineering, Sharif University of Technology,
More information3.032 Problem Set 2 Solutions Fall 2007 Due: Start of Lecture,
3.032 Problem Set 2 Solutions Fall 2007 Due: Start of Lecture, 09.21.07 1. In the beam considered in PS1, steel beams carried the distributed weight of the rooms above. To reduce stress on the beam, it
More informationSpace frames. b) R z φ z. R x. Figure 1 Sign convention: a) Displacements; b) Reactions
Lecture notes: Structural Analsis II Space frames I. asic concepts. The design of a building is generall accomplished b considering the structure as an assemblage of planar frames, each of which is designed
More information1.1 The Equations of Motion
1.1 The Equations of Motion In Book I, balance of forces and moments acting on an component was enforced in order to ensure that the component was in equilibrium. Here, allowance is made for stresses which
More informationExercise solutions: concepts from chapter 5
1) Stud the oöids depicted in Figure 1a and 1b. a) Assume that the thin sections of Figure 1 lie in a principal plane of the deformation. Measure and record the lengths and orientations of the principal
More informationMicrostructural Randomness and Scaling in Mechanics of Materials. Martin Ostoja-Starzewski. University of Illinois at Urbana-Champaign
Microstructural Randomness and Scaling in Mechanics of Materials Martin Ostoja-Starzewski University of Illinois at Urbana-Champaign Contents Preface ix 1. Randomness versus determinism ix 2. Randomness
More informationFEA A Guide to Good Practice. What to expect when you re expecting FEA A guide to good practice
FEA A Guide to Good Practice What to expect when you re expecting FEA A guide to good practice 1. Background Finite Element Analysis (FEA) has transformed design procedures for engineers. Allowing more
More informationMathematical aspects of mechanical systems eigentones
Seminar: Vibrations and Structure-Borne Sound in Civil Engineering Theor and Applications Mathematical aspects of mechanical sstems eigentones Andre Kuzmin April st 6 Abstract Computational methods of
More informationEVALUATION OF THERMAL TRANSPORT PROPERTIES USING A MICRO-CRACKING MODEL FOR WOVEN COMPOSITE LAMINATES
EVALUATION OF THERMAL TRANSPORT PROPERTIES USING A MICRO-CRACKING MODEL FOR WOVEN COMPOSITE LAMINATES C. Luo and P. E. DesJardin* Department of Mechanical and Aerospace Engineering Universit at Buffalo,
More informationThe Control-Volume Finite-Difference Approximation to the Diffusion Equation
The Control-Volume Finite-Difference Approimation to the Diffusion Equation ME 448/548 Notes Gerald Recktenwald Portland State Universit Department of Mechanical Engineering gerr@mepdedu ME 448/548: D
More informationDynamics and control of mechanical systems
JU 18/HL Dnamics and control of mechanical sstems Date Da 1 (3/5) 5/5 Da (7/5) Da 3 (9/5) Da 4 (11/5) Da 5 (14/5) Da 6 (16/5) Content Revie of the basics of mechanics. Kinematics of rigid bodies coordinate
More informationPrediction of Micromechanical Behaviour of Elliptical Frp Composites
Prediction of Micromechanical Behaviour of Elliptical Frp Composites Kiranmayee.Nerusu Dept. of Mechanical Engg. P. V. P. Siddhartha Institute of Technology, Vijayawada 520 007, A.P, India. P. Phani Prasanthi
More informationRANGE CONTROL MPC APPROACH FOR TWO-DIMENSIONAL SYSTEM 1
RANGE CONTROL MPC APPROACH FOR TWO-DIMENSIONAL SYSTEM Jirka Roubal Vladimír Havlena Department of Control Engineering, Facult of Electrical Engineering, Czech Technical Universit in Prague Karlovo náměstí
More informationDeflections and Strains in Cracked Shafts due to Rotating Loads: A Numerical and Experimental Analysis
Rotating Machinery, 10(4): 283 291, 2004 Copyright c Taylor & Francis Inc. ISSN: 1023-621X print / 1542-3034 online DOI: 10.1080/10236210490447728 Deflections and Strains in Cracked Shafts due to Rotating
More informationLATERAL BUCKLING ANALYSIS OF ANGLED FRAMES WITH THIN-WALLED I-BEAMS
Journal of arine Science and J.-D. Technolog, Yau: ateral Vol. Buckling 17, No. Analsis 1, pp. 9-33 of Angled (009) Frames with Thin-Walled I-Beams 9 ATERA BUCKING ANAYSIS OF ANGED FRAES WITH THIN-WAED
More informationME 7502 Lecture 2 Effective Properties of Particulate and Unidirectional Composites
ME 75 Lecture Effective Properties of Particulate and Unidirectional Composites Concepts from Elasticit Theor Statistical Homogeneit, Representative Volume Element, Composite Material Effective Stress-
More informationModeling of the Bending Stiffness of a Bimaterial Beam by the Approximation of One-Dimensional of Laminated Theory
. Flores-Domínguez Int. Journal of Engineering Research and Applications RESEARCH ARTICLE OPEN ACCESS odeling of the Bending Stiffness of a Bimaterial Beam by the Approimation of One-Dimensional of Laminated
More informationEstimation Of Linearised Fluid Film Coefficients In A Rotor Bearing System Subjected To Random Excitation
Estimation Of Linearised Fluid Film Coefficients In A Rotor Bearing Sstem Subjected To Random Ecitation Arshad. Khan and Ahmad A. Khan Department of Mechanical Engineering Z.. College of Engineering &
More informationTHREE DIMENSIONAL STRESS ANALYSIS OF THE T BOLT JOINT
THREE DIMENSIONAL STRESS ANALYSIS OF THE T BOLT JOINT Víctor Martínez 1, Alfredo Güemes 2, Norbert Blanco 1, Josep Costa 1 1 Escola Politècnica Superior. Universitat de Girona. Girona, Spain (17071) 2
More informationAxial force-moment interaction in the LARSA hysteretic beam element
Axial force-moment interaction in the LARSA hsteretic beam element This document briefl discusses the modeling of tri-axial interaction (i.e. between axial force and bending moments) in the LARSA beam
More informationAn accelerated predictor-corrector scheme for 3D crack growth simulations
An accelerated predictor-corrector scheme for 3D crack growth simulations W. Weber 1 and G. Kuhn 2 1,2 1 Institute of Applied Mechanics, University of Erlangen-Nuremberg Egerlandstraße 5, 91058 Erlangen,
More informationDetermination of Young s modulus of glass by Cornu s apparatus
Determination of Young s modulus of glass b Cornu s apparatus Objective To determine Young s modulus and Poisson s ratio of a glass plate using Cornu s method. Theoretical Background Young s modulus, also
More informationCitation Key Engineering Materials, ,
NASITE: Nagasaki Universit's Ac Title Author(s) Interference Analsis between Crack Plate b Bod Force Method Ino, Takuichiro; Ueno, Shohei; Saim Citation Ke Engineering Materials, 577-578, Issue Date 2014
More informationAircraft Structures Structural & Loading Discontinuities
Universit of Liège Aerospace & Mechanical Engineering Aircraft Structures Structural & Loading Discontinuities Ludovic Noels Computational & Multiscale Mechanics of Materials CM3 http://www.ltas-cm3.ulg.ac.be/
More informationEffects of birefringence on Fizeau interferometry that uses polarization phase shifting technique
Effects of birefringence on Fizeau interferometr that uses polarization phase shifting technique Chunu Zhao, Dongel Kang and James H. Burge College of Optical Sciences, the Universit of Arizona 1630 E.
More informationOpen-hole compressive strength prediction of CFRP composite laminates
Open-hole compressive strength prediction of CFRP composite laminates O. İnal 1, A. Ataş 2,* 1 Department of Mechanical Engineering, Balikesir University, Balikesir, 10145, Turkey, inal@balikesir.edu.tr
More informationShear and torsion correction factors of Timoshenko beam model for generic cross sections
Shear and torsion correction factors of Timoshenko beam model for generic cross sections Jouni Freund*, Alp Karakoç Online Publication Date: 15 Oct 2015 URL: http://www.jresm.org/archive/resm2015.19me0827.html
More informationtorsion equations for lateral BucKling ns trahair research report r964 July 2016 issn school of civil engineering
TORSION EQUATIONS FOR ATERA BUCKING NS TRAHAIR RESEARCH REPORT R96 Jul 6 ISSN 8-78 SCHOO OF CIVI ENGINEERING SCHOO OF CIVI ENGINEERING TORSION EQUATIONS FOR ATERA BUCKING RESEARCH REPORT R96 NS TRAHAIR
More informationComposite Structures- Modeling, FEA, Optimization and Diagnostics
Composite Structures- Modeling, FEA, Optimization and Diagnostics Ratan Jha Mechanical and Aeronautical Engineering Clarkson University, Potsdam, NY Composite Laminate Modeling Refined Higher Order Displacement
More informationCRITERIA FOR SELECTION OF FEM MODELS.
CRITERIA FOR SELECTION OF FEM MODELS. Prof. P. C.Vasani,Applied Mechanics Department, L. D. College of Engineering,Ahmedabad- 380015 Ph.(079) 7486320 [R] E-mail:pcv-im@eth.net 1. Criteria for Convergence.
More informationINTERFACIAL STRENGTH EVALUATION IN A GLASS FIBER REINFORCED COMPOSITE USING CRUCIFORM SPECIMEN METHOD
INTERFACIAL STRENGTH EVALUATION IN A GLASS FIBER REINFORCED COMPOSITE USING CRUCIFORM SPECIMEN METHOD Shinji Ogihara, oriaki Sakamoto 2, Hajime Kato 2, Satoshi Kobayashi 3 and Ichiro Ueno Dept. of Mechanical
More informationScienceDirect. Analogy of stress singularities analysis between piezoelectric materials and anisotropic materials
Available online at www.sciencedirect.com Scienceirect Procedia Materials Science ( 4 ) 767 77 th European Conference on Fracture (ECF) Analog of stress singularities analsis between piezoelectric materials
More informationHydraulic fracturing in unconventional shale gas reservoirs. Dr.-Ing. Johannes Will, Dynardo GmbH, Weimar, Germany
Hydraulic fracturing in unconventional shale gas reservoirs Dr.-Ing. Johannes Will, Dynardo GmbH, Weimar, Germany Founded: 2001 (Will, Bucher, CADFEM International) More than 35 employees, offices at Weimar
More information