VUMAT for Fabric Reinforced Composites
|
|
- Lisa Reeves
- 5 years ago
- Views:
Transcription
1 VUMAT or Fabric Reinorce Composites. Introuction This ocument escribes a constitutive mo or abric reinorce composites that was introuce in Abaqus/Exicit 6.8. The mo has been imemente as a built-in VUMAT user subroutine. It can be accesse by naming your material such that it begins with the string ABQ_PLY_FABRIC, e.g. ABQ_PLY_FABRIC_. The mo is currently supporte or ane-stress ements; this inclues shl (S4R an S3R, continuum shl (SC6R an SC8R, ane stress (CPS amily an membrane (M3D amily ements. User materials are currently not supporte in ABAQUS/Exicit or small-strain shl ements (S4RS. When user-eine materials are emoye to eine the material response o shl ements, ABAQUS/Exicit cannot calculate a eault value or the transverse shear stiness o the ement. Hence, you must manually eine the ement's transverse shear stiness. See Shl section behavior, Section o the Version 6.8 ABAQUS Analysis User s Manual, or guiines on choosing this stiness. This ocument escribes the basic equations o the constitutive mo an provies etaile inormation o the user interace to the VUMAT imementation in ABAQUS.. Continuum amage mo or abric reinorce composites A schematic representation o the geometry o the woven abric reinorcement consiere in the constitutive mo is shown in Figure. The iber irections are assume to be orthogonal. Figure : Schematic representation o woven abric. Fibers are aligne with irections an o a local coorinate system. The constitutive stress-strain rations are ormulate in a local Cartesian coorinate system with base vectors aligne with the iber irections, as shown in Figure. The abric-reinorce y is moe as a homogeneous orthotropic astic material with the potential to sustain progressive stiness egraation ue to iber/matrix cracking, an astic eormation uner shear loaing. The ierent aspects o the mo are iscusse next. Copyright Dassault Systèmes, 8 Page o 9
2 VUMAT or Fabric Reinorce Composites.. Elastic stress-strain rations It is assume that the astic stress-strain rations are given by orthotropic amage asticity. Reerre to a local coorinate system aligne with the iber irections (Figure the astic rations take the orm: ν ( E E σ ν σ E ( E σ ( ( The amage variables an are associate with iber racture along the an irections respectivy, whereas is rate to matrix micro-cracking ue to shear eormation. The mo ierentiates between tensile an compressive iber ailure moes by activating the corresponing amage variable epening on the stress state in the iber irections. Thus: < σ < σ < σ < σ + + ; + ( σ σ σ σ In orer to incorporate ierent initial (unamage stiness in tension an compression, the values o the astic constants E, E an ν are assume to take their tensile or compressive values epening on the sign o tr ( +... Fiber response The material response along the iber irections is characterize with amage asticity. It is assume that the iber amage variables are a unction o the corresponing eective stress, that is: The eective stresses are eine as: ( ( < σ ( < σ ( ( ( < σ ( < σ ( (3 (4 In orer to simiy notation, an inex will be use in subsequent iscussions, such that it takes the value (+/, (+/, epening on the sign o the corresponing stresses. Thus, the above our equations or the amage variables are rewritten as: (5 (
3 VUMAT or Fabric Reinorce Composites It is note that the eective stresses are irectly rate to the thermoynamic orces, Y, that are work conjugate to the amage variables, through the rationship E Y. Thereore the above equation states that the iber amage variables epen only on the corresponing thermoynamic orce. At any given time the astic omain is eine in terms o the amage activation unctions, F, as F φ r (6 The unctions φ provie a criterion or iber ailure an are assume to take the orm σ φ ; ( +,, +, X (7 where X are the tensile/compressive strengths or uniaxial loaing along the iber irections. The amage threshols, r, are initially set to one. Ater amage activation ( φ they increase with increasing amage accoring to: r ( t φ ( τ (8 τ t The einition ensures that the amage threshols are non-ecreasing quantities ( r ( t. The amage threshols are assume to obey the Kuhn-Tucker comementary conitions: an the consistency conition: F r r F (9 r F ( Note that the ormulation can be easily enhance to take the eects o amage upon loa reversal into account. For exame, compressive amage will usually egrae the tensile response i the loaing is reverse rom compression to tension. On the other han, tensile cracks close uner compressive loaing an have little eect on the compressive response. The evolution o the amage variables are a unction o the amage threshols an the racture energy per unit area uner uniaxial tensile/compressive loaing,. The ormulation o the amage evolution law ensures that the amage variables are monotonically increasing quantities. It also ensures that the correct amount o energy is issipate when the lamina is subjecte to uniaxial loaing conitions along the iber irections. The evolution o the amage variables is given by the equation: exp( A ( r ;, ( r where g L A ( L c g c 3
4 VUMAT or Fabric Reinorce Composites Here, L c is the characteristic length o the ement, tensile/compressive loaing, an amage initiation: is the racture energy per unit area uner uniaxial g is the astic energy ensity (i.e. per unit volume at the point o X g (3 E The ormulation o the amage evolution law ensures that the amage variables are monotonically increasing quantities. It also ensures that the correct amount o energy is issipate when the lamina is subjecte to uniaxial loaing conitions along the iber irections. For instance, uner uniaxial tensile loaing in the iber irection, the issipate energy per unit area is equal to the racture energy. This hols true provie that: (4 g g Lc Lc < L The ormulation thereore imposes a restriction on the imum ement size that can be use to accuraty capture the right amount o energy issipation uring racture. I the characteristic ement size o the FE mesh is greater than L, the analysis will over-preict the energy issipation. Note that, rom Abaqus 6.-EF onwars, the critical ement length or each o the ABQ_PLY_FABRIC materials an a representative list o ements that excee this criterion are printe in the.sta ile..3. Shear response As mentione earlier, the shear response is ominate by the non-linear behavior o the matrix, which inclues both asticity an stiness egraation ue to matrix microcracking. The main ingreients o the shear response are iscusse bow..3. Elasticity The astic rations give the eective (unamage stress in terms o astic strain:.3. Plasticity Yi unction: σ ( (5 ( F σ σ ( (6 The harening unction is assume to be o the orm: ~ + Flow rule: Assuming associate low, then p ( σ y C ( (7 F ( ~ ~ sign σ (8 σ 4
5 VUMAT or Fabric Reinorce Composites The evolution o the astic work uring yiing is given as U σ ( ( (9.3.3 Damage The astic omain is eine in terms o the amage activation unction, F, as F φ r ( The unction φ provies the criterion or initiation o shear amage o the matrix, which is assume to be o the orm ~σ φ ( S Here σ σ /( is the eective shear stress, an S is the shear stress or initiation o matrix amage. The amage threshol, r, is initially set to one an increases ater amage activation ( φ accoring to ( τ τ t r t φ ( ( Finally, base on [] it is assume that the shear amage variable increases with the logarithm o r until a imum value is reache. Thus: min( ln( r, (3 where, an are material properties..4. Element Detion The VUMAT provies two options to ete ements: The ement is ete when any one tensile/compressive amage variable along the iber irections reaches a imum speciie value, or, or when the astic strain ue to shear eormation reaches a imum speciie value,. This option is activate by setting the lag ldflag. The ement is ete when the amage variables along both iber irections reach a imum speciie value,, or when the astic strain ue to shear eormation reaches a imum speciie value,. This option is activate by setting the lag ldflag. These two options can be combine with a eormation-base ement etion criterion base on the values o the imum ( ˆ an minimum ( ˆmin < principal logarithmic strains that the ement can sustain beore it gets ete. 5
6 VUMAT or Fabric Reinorce Composites 3. Calibration proceure The astic constants an the iber tension/compression strengths, X, are easily measure rom stanar coupon tests in uniaxial tension/compression loaing o /9 laminates. The calibration o amage evolution in the iber ailure moes is base on the racture energy per unit area o the material,, which can be measure experimentally. The shear response is usually calibrate with a cyclic tensile test on a ±45 laminate, where the strains along the iber irections can be neglecte. Figure shows the typical shear response o a abric reinorce composite. It is note that the unloaing/roaing paths in this igure correspon to an iealization o the actual response, which usually exhibits hysteretic behavior. The igure will serve as the starting point or the iscussion o a general calibration proceure or the parameters that enter the amage an asticity equations. σ σ y S ( Figure : Schematic representation o typical shear response o a abric reinorce composite. The lev o amage can be measure rom the ratio o the unloaing stiness to the initial (unamage astic stiness. This allows us to compute pairs o stress-amage values, ( σ,, or each unloaing curve. This ata can be represente in the space o versus ln(, where σ /(. A linear it o the ata provies the values o (slope o the line an S (intersection with the horizontal axis as shown in Figure 3. Sometimes the amage ata shows inication o a saturation value, which woul be use to etermine. Otherwise a value o shoul be use. Figure 3: Calibration o the shear amage parameters an S ln(s ln( 6
7 VUMAT or Fabric Reinorce Composites Finally, or each unloaing curve in Figure, the astic strain at the onset o unloaing is etermine rom the value o resiual eormation in the unloae state. The values o ~ ( σ, at the onset o unloaing are then use to it the parameters o the harening curve, as illustrate in Figure 4. ~σ ~ σ y + p y C ( Figure 4: Calibration o the shear harening curve 4. User interace To activate the mo or abric reinorce composites the material name must start with the string ABQ_PLY_FABRIC, e.g. ABQ_PLY_FABRIC_. A synopsis o the interace is shown bow. The number o solution epenent variables (uner keywor option *DEPVAR is 6, an the DELETE parameter is equal to 6. Reer to Table or a etaile escription o each material constant speciie in the keywor interace. *MATERIAL, NAME ABQ_PLY_FABRIC *DENSITY ρ *USER MATERIAL, CONSTANTS4 ** Line : E, E, ν +,, E, E, ν ** Line : X, X, X, X, S ** Line 3:,,, ** Line 4: y, C, p ** Line 5: ldflag,, *DEPVAR, DELETE6 6,,, ˆ, min ˆ Table : User material constants or the abric reinorce composite material mo. 7
8 VUMAT or Fabric Reinorce Composites LINE: Pos. Symbol Description 3 + E Young s moulus along iber irection when tr ( E Young s moulus along iber irection when tr ( ν Poisson ratio when tr ( 4 Shear moulus E Young s moulus along iber irection when tr ( < E Young s moulus along iber irection when tr ( < ν Poisson ratio when tr ( < 8 Not use LINE: Damage initiation coeicients Pos. Symbol Description X Tensile strength along iber irection X Compressive strength along iber irection 3 X Tensile strength along iber irection 4 X Compressive strength along iber irection 5 S Shear stress at the onset o shear amage 6-8 Not use LINE: 3 Damage evolution coeicients Pos. Symbol Description 3 4 Energy per unit area or tensile racture along iber irection Energy per unit area or compressive racture along iber irection Energy per unit area or tensile racture along iber irection Energy per unit area or compressive racture along iber irection 5 Parameter in the equation o shear amage 6 Maximum shear amage 7-8 Not use 8
9 VUMAT or Fabric Reinorce Composites LINE: 4 Shear asticity coeicients: ( + C ( Pos. Symbol Description y y Initial eective shear yi stress C Coeicient in harening equation 3 p Power term in harening equation 4-8 Not use p LINE: 5 Controls or material point ailure Pos. Symbol Description ldflag Element etion lag: ldflag: Element is not ete (eault ldflag: Element is ete when either iber ails, or, or when ldflag: Element is ete when both ibers ail,, or when. Maximum value o amage variable use in ement etion criterion Maximum value o equivalent astic strain or ement etion criterion. 3 (A value o zero means that is not use as criterion or ement etion Maximum (positive principal logarithmic strain beyon which the 4 ˆ ement will get ete. Ignore i zero, not speciie, or ldflag. Minimum (negative principal logarithmic strain beyon which the 5 min ˆ ement will get ete. Ignore i zero, not speciie, or ldflag. 6-8 Not use 9
10 VUMAT or Fabric Reinorce Composites 5. Output In aition to the stanar (material-inepenent output variables in Abaqus/Exicit or stressisacement ements (such as stress, S, strain, LE, ement STATUS, etc. the ollowing output variables have a special meaning or the user material or abric-reinorce composites: Output Variable Symbol Description SDV Tensile amage along iber irection SDV Compressive amage along iber irection SDV3 Tensile amage along iber irection SDV4 Compressive amage along iber irection SDV5 Shear amage SDV6 r Tensile amage threshol along iber irection SDV7 r Compressive amage threshol along iber irection SDV8 r Tensile amage threshol along iber irection SDV9 SDV r r Compressive amage threshol along iber irection Shear amage threshol SDV SDV SDV3 SDV4 SDV5 SDV6 Equivalent astic strain Elastic strain component Elastic strain component Not use Elastic strain component MpStatus Material point status: i active, i aile. 6. Reerences [] Alastair F. Johnson an Jose Simon, Moing Fabric Reinorce Composites uner Impact Loas. In EUROMECH 4: Impact an Damage Tolerance o Composite Materials an Structures. Imperial College o Science Technology & Meicine, Lonon 7-9 September 999.
Strength Analysis of CFRP Composite Material Considering Multiple Fracture Modes
5--XXXX Strength Analysis of CFRP Composite Material Consiering Multiple Fracture Moes Author, co-author (Do NOT enter this information. It will be pulle from participant tab in MyTechZone) Affiliation
More 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 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 informationMICROMECHANICAL FAILURE ANALYSIS OF UNIDIRECTIONAL FIBER-REINFORCED COMPOSITES UNDER IN-PLANE AND TRANSVERSE SHEAR
THE 19 TH INTERNATIONAL CONFERENCE ON COMPOSITE MATERIALS MICROMECHANICAL FAILURE ANALYSIS OF UNIDIRECTIONAL FIBER-REINFORCED COMPOSITES UNDER IN-PLANE AND TRANSVERSE SHEAR Lei Yang*, Ying Yan, Zhiguo
More informationOn interface property characterization and performance of fiber-reinforced cementitious composites
Concrete Science an Engineering, Vol. 1, September 1999, pp 173-174 On interace property characterization an perormance o iber-reinorce cementitious composites Z. Lin 1, T. Kana an V. C. Li 1 (1 Avance
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 informationDamage in fiber bundle models
Eur Phys J B 17, 269 279 (2) THE EUROPEAN PHYSICAL JOURNAL B c EDP Sciences Società Italiana i isica Springer-Verlag 2 Damage in iber bunle moels Kun 1,2,a, S Zapperi 3,4, an HJ Herrmann 1,3 1 ICA 1, University
More 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 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 informationA DYNAMIC MODEL OF POLYELECTROLYTE GELS
A DYNAMIC MODEL OF POLYELECTROLYTE GELS M.CARME CALDERER, HAORAN CHEN, CATHERINE MICEK, AND Y. MORI Abstract. We erive a moel o the couple mechanical an electrochemical eects o polyelectrolyte gels. We
More informationRecommendations: Part 7: Transient Creep for service and accident conditions
Materials an Structures/Matériaux et Constructions, Vol. 31, June 1998, pp 290-295 RILEM TECHNICAL COMMITTEES RILEM TC 129-MHT: TEST METHODS FOR MECHANICAL PROPERTIES OF CONCRETE AT HIGH TEMPERATURES Recommenations:
More informationA Short Note on Self-Similar Solution to Unconfined Flow in an Aquifer with Accretion
Open Journal o Flui Dynamics, 5, 5, 5-57 Publishe Online March 5 in SciRes. http://www.scirp.org/journal/oj http://x.oi.org/.46/oj.5.57 A Short Note on Sel-Similar Solution to Unconine Flow in an Aquier
More 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 informationS. Srinivasan, Technip Offshore, Inc., Houston, TX
9 th ASCE Specialty Conerence on Probabilistic Mechanics and Structural Reliability PROBABILISTIC FAILURE PREDICTION OF FILAMENT-WOUND GLASS-FIBER Abstract REINFORCED COMPOSITE TUBES UNDER BIAXIAL LOADING
More informationON AN ESTIMATION OF THE DAMPING PROPERTIES OF WOVEN FABRIC COMPOSITES
ON AN ESTIMATION OF THE DAMPING PROPERTIES OF WOVEN FABRIC COMPOSITES Masaru Zao 1, Tetsusei Kurashii 1, Yasumasa Naanishi 2 and Kin ya Matsumoto 2 1 Department o Management o Industry and Technology,
More informationAXIALLY LOADED FRP CONFINED REINFORCED CONCRETE CROSS-SECTIONS
AXIALLY LOADED FRP CONFINED REINFORCED CONCRETE CROSS-SECTIONS Bernát Csuka Budapest University o Technology and Economics Department o Mechanics Materials and Structures Supervisor: László P. Kollár 1.
More informationOptimal design of reinforced concrete T-sections in bending
Engineering Structures 5 (003) 951 964 www.elsevier.com/locate/engstruct Optimal esign o reinorce concrete T-sections in ening C.C. Ferreira a,, M.H.F.M. Barros a, A.F.M. Barros a Department o Civil Engineering,
More informationDelamination Prediction and Non-local Averaging using a Composite Micro-Mechanical Model
Delamination Preiction an Non-local Averaging using a Composite Micro-Mechanical Moel Ala Tabiei *, Saneep Meikona * * Department o Mechanical an Materials Engineering, University o Cincinnati, Cincinnati,
More informationTHREE-DIMENSIONAL THERMO-PORO-MECHANICAL MODELING OF RESERVOIR STIMULATION AND INDUCED MICROSEISMICITY IN GEOTHERMAL RESERVOIR
PROCEEDINGS, Thirty-Sixth Workshop on Geothermal Reservoir Engineering Stanor University, Stanor, Caliornia, January 31 - February 3, 2011 SGP-TR-191 THREE-DIMENSIONAL THERMO-PORO-MECHANICAL MODELING OF
More 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 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 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 informationFATIGUE DURABILITY OF CONCRETE EXTERNALLY STRENGTHENED WITH FRP SHEETS
FATIGUE DURABILITY OF CONCRETE EXTERNALLY STRENGTHENED WITH FRP SHEETS H. Diab () and Zhishen Wu () () Department o Urban and Civil Engineering, Ibaraki University, Japan Abstract A primary concern o the
More information'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 informationDUCTILE FRACTURE PREDICTION IN THIN METALLIC PLATES SUBJECTED TO EXPLOSIONS
DUCTLE FRACTURE PREDCTON N TN METALLC PLATES SUBJECTED TO EXPLOSONS Łukasz Pyrzowski 1, Paweł Kłosowski, Krzyszto Woznica 3, Olivier Pennetier 4 1 Gansk University o Technology, Narutowicza str. 11/1,
More informationConservation Laws. Chapter Conservation of Energy
20 Chapter 3 Conservation Laws In orer to check the physical consistency of the above set of equations governing Maxwell-Lorentz electroynamics [(2.10) an (2.12) or (1.65) an (1.68)], we examine the action
More informationResearch on Thin Film Thickness Uniformity for Deposition of Rectangular Planar Sputtering Target
Available online at www.scienceirect.com Physics Proceia 3 (01 ) 903 913 18 th International Vacuum Congress (IVC-18) Research on Thin Film Thickness Uniormity or Deposition o Rectangular Planar puttering
More informationCharacterization of Internal State Variable for fiber fracture in UD Composite
Characterization o Internal State ariable or iber racture in UD Composite Modris Megnis,, Povl Brondsted 3, Sai A. Rehman 4, Tanveer Ahmad 5 Summary The continuum damage mechanics is used to describe the
More informationA Formulation for Hyperelastic Damaged Materials
Paper 0123456789 Civil-Comp Press, 2015 Proceeings of the Fifteenth International Conference on Civil, Structural an Environmental Engineering Computing, J. Kruis, Y. Tsompanakis an B.H.V. Topping, (Eitors),
More informationensembles When working with density operators, we can use this connection to define a generalized Bloch vector: v x Tr x, v y Tr y
Ph195a lecture notes, 1/3/01 Density operators for spin- 1 ensembles So far in our iscussion of spin- 1 systems, we have restricte our attention to the case of pure states an Hamiltonian evolution. Toay
More informationOn the control of the load increments for a proper description of multiple delamination in a domain decomposition framework
On the control of the loa increments for a proper escription of multiple elamination in a omain ecomposition framework Olivier Allix, Pierre Kerfrien, Pierre Gosselet To cite this version: Olivier Allix,
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 informationLongitudinal Waves in a Rotating Solid Cylinder Immersed in an Inviscid Fluid Using Chebyshev Polynomial
Longituinal Waves in a Rotating Soli Cyliner Immerse in an Invisci Flui Using Chebyshev Polynomial R. Selvamani *1, P.Ponnusamy 1 Department o Mathematics, Karunya University, Coimbatore, TamilNau, Inia
More informationInfluence of Interfacial Delamination on Channel Cracking of Elastic Thin Films. Haixia Mei, Yaoyu Pang, Rui Huang
Inluence o Interacial Delamination on Channel Cracking o Elastic Thin Films Haixia Mei, Yaoyu Pang, Rui Huang Department o Aerospace Engineering an Engineering Mechanics, University o Texas, Austin, TX
More informationDesign criteria for Fiber Reinforced Rubber Bearings
Design criteria or Fiber Reinorced Rubber Bearings J. M. Kelly Earthquake Engineering Research Center University o Caliornia, Berkeley A. Calabrese & G. Serino Department o Structural Engineering University
More information1 Introuction In the past few years there has been renewe interest in the nerson impurity moel. This moel was originally propose by nerson [2], for a
Theory of the nerson impurity moel: The Schrieer{Wol transformation re{examine Stefan K. Kehrein 1 an nreas Mielke 2 Institut fur Theoretische Physik, uprecht{karls{universitat, D{69120 Heielberg, F..
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 informationDelamination Prediction of Uni-Directional Composite Laminates using Shell Elements and a Strain Rate Dependent Micro-mechanical Material Model
14 th International LS-DYNA Users Conerence Session: Coposites Delaination Preiction o Uni-Directional Coposite Lainates using Shell Eleents an a Strain Rate Depenent Micro-echanical Material Moel Saneep
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 informationA finite element model of ultrasonically assisted drilling in carbon/epoxy composites
Available online at wwwsciencedirectcom Procedia CIRP 8 (13 ) 141 146 14 th CIRP Conerence on Modeling o Machining Operations (CIRP CMMO) A inite element model o ultrasonically assisted drilling in carbon/epoxy
More informationComputed Tomography Notes, Part 1. The equation that governs the image intensity in projection imaging is:
Noll 3 CT Notes : Page Compute Tomograph Notes Part Challenges with Projection X-ra Sstems The equation that governs the image intensit in projection imaging is: z I I ep µ z Projection -ra sstems are
More informationComputed Tomography Notes, Part 1. The equation that governs the image intensity in projection imaging is:
Noll 6 CT Notes : Page Compute Tomograph Notes Part Challenges with Projection X-ra Sstems The equation that governs the image intensit in projection imaging is: z I I ep μ z Projection -ra sstems are
More informationMomentum and Energy. Chapter Conservation Principles
Chapter 2 Momentum an Energy In this chapter we present some funamental results of continuum mechanics. The formulation is base on the principles of conservation of mass, momentum, angular momentum, an
More informationThermal runaway during blocking
Thermal runaway uring blocking CES_stable CES ICES_stable ICES k 6.5 ma 13 6. 12 5.5 11 5. 1 4.5 9 4. 8 3.5 7 3. 6 2.5 5 2. 4 1.5 3 1. 2.5 1. 6 12 18 24 3 36 s Thermal runaway uring blocking Application
More information1. Filling an initially porous tube under a constant head imposed at x =0
Notes on Moving Bounary problems, Voller U o M, volle00@umn.eu. Filling an initially porous tube uner a constant hea impose at x =0 Governing equation is base on calculating the water volume lux by the
More 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 informationDesign and assessment of structures
Design an assessment o structures Limit state esign requires the structure to satisy in two principal criterias: the ultimate limit state (ULS) an the serviceability limit state (SLS) In Europe, the Limit
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 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 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 informationDuctility and Failure Modes of Single Reinforced Concrete Columns. Hiromichi Yoshikawa 1 and Toshiaki Miyagi 2
Ductility an Failure Moes of Single Reinforce Concrete Columns Hiromichi Yoshikawa 1 an Toshiaki Miyagi 2 Key Wors: seismic capacity esign, reinforce concrete column, failure moes, eformational uctility,
More informationDesign Method for RC Building Structure Controlled by Elasto-Plastic Dampers Using Performance Curve
Design Metho or RC Biling Strctre Controlle by Elasto-Plastic Dampers Using Perormance Crve W. P Whan University o Technology, China K. Kasai Tokyo Institte o Technology, Japan SUMMARY: This paper proposes
More informationModeling of vegetated rivers for inbank and overbank ows
Loughborough University Institutional Repository Moeling o vegetate rivers or inbank an overbank ows This item was submitte to Loughborough University's Institutional Repository by the/an author. Citation:
More informationLife Prediction Under Multiaxial Fatigue
Lie Prediction Under Multiaxial Fatigue D. Ramesh and M.M. Mayuram Department o Mechanical Engineering Indian Institute o Technology, Madras Chennai-600 036 (India) e-mail: mayuram@iitm.ac.in ABSTRACT
More informationX-ray Diffraction from Materials
X-ray Diffraction from Materials 8 Spring Semester Lecturer; Yang Mo Koo Monay an Wenesay 4:45~6: Plastically ll eforme metal -Macro strain: uniform elastic strain over relatively large istance iffraction
More informationMath 342 Partial Differential Equations «Viktor Grigoryan
Math 342 Partial Differential Equations «Viktor Grigoryan 6 Wave equation: solution In this lecture we will solve the wave equation on the entire real line x R. This correspons to a string of infinite
More informationSteady Load Failure Theories
tea Loa Failure Theories Ductile Materials Uniaxial tress/train Fiel Maximum-Normal-tress Maximum-Normal-train Maximum-hear-tress Distortion-Energ hear-energ Von Mises-Henck Octaheral-hear-tress Internal-Friction
More informationAn Approach for Design of Multi-element USBL Systems
An Approach for Design of Multi-element USBL Systems MIKHAIL ARKHIPOV Department of Postgrauate Stuies Technological University of the Mixteca Carretera a Acatlima Km. 2.5 Huajuapan e Leon Oaxaca 69000
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 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 informationModule 2. DC Circuit. Version 2 EE IIT, Kharagpur
Moule 2 DC Circuit Lesson 9 Analysis of c resistive network in presence of one non-linear element Objectives To unerstan the volt (V ) ampere ( A ) characteristics of linear an nonlinear elements. Concept
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 informationFiber / Toughened Epoxy Composites under
Fatigue Behavior o Impact-Damaged Carbon Fiber / Toughened Epoxy Composites under Compressive loading Toshio Ogasawara, Hisaya Katoh, unao ugimoto, and Takashi Ishikawa Advanced d Composite Technology
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 informationFinite element modeling incorporating nonlinearity of material behavior based on the fib Model Code 2010
Peer-reviewed & Open access journal www.academicpublishingplatorms.com Finite element modeling incorporating non-linearity o material behavior ATI - Applied Technologies & Innovations Volume 5 Issue November
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 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 informationBasic Thermoelasticity
Basic hermoelasticity Biswajit Banerjee November 15, 2006 Contents 1 Governing Equations 1 1.1 Balance Laws.............................................. 2 1.2 he Clausius-Duhem Inequality....................................
More informationANALYSIS OF SECTION. Behaviour of Beam in Bending
ANALYSIS OF SECTION Behaviour o Beam in Bening Conier a imply upporte eam ujecte to graually increaing loa. The loa caue the eam to en an eert a ening moment a hown in igure elow. The top urace o the eam
More informationFinite element analysis of the dynamic behavior of concrete-filled double skin steel tubes (CFDST) under lateral impact with fixed-simply support
Structures Uner Shock an Impact XIII 383 Finite element analysis of the ynamic behavior of concrete-fille ouble skin steel tubes (CFDST) uner lateral impact with fixe-simply support M. Liu & R. Wang Architecture
More informationA FINITE ELEMENT MODEL TO PREDICT MULTI- AXIAL STRESS-STRAIN RESPONSE OF CERAMIC MATRIX COMPOSITES WITH STRAIN INDUCED DAMAGE
A FINITE ELEMENT MODEL TO PREDICT MULTI- AXIAL STRESS-STRAIN RESPONSE OF CERAMIC MATRIX COMPOSITES WITH STRAIN INDUCED DAMAGE Daxu Zhang and D. R. Hayhurst School of Mechanical, Aerospace and Civil Engineering,
More informationAn Upper Bound on the Minimum Distance of Serially Concatenated Convolutional Codes
SUBMITTED TO IEEE TRANSACTIONS ON INFORMATION THEORY, FEBRUARY 2004 1 An Upper Boun on the Minimum Distance o Serially Concatenate Convolutional Coes Alberto Perotti, Member, IEEE, an Sergio Beneetto,
More informationConvective heat transfer
CHAPTER VIII Convective heat transfer The previous two chapters on issipative fluis were evote to flows ominate either by viscous effects (Chap. VI) or by convective motion (Chap. VII). In either case,
More information3.2 Shot peening - modeling 3 PROCEEDINGS
3.2 Shot peening - moeling 3 PROCEEDINGS Computer assiste coverage simulation François-Xavier Abaie a, b a FROHN, Germany, fx.abaie@frohn.com. b PEENING ACCESSORIES, Switzerlan, info@peening.ch Keywors:
More 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 informationEC5555 Economics Masters Refresher Course in Mathematics September 2013
EC5555 Economics Masters Reresher Course in Mathematics September 3 Lecture 5 Unconstraine Optimization an Quaratic Forms Francesco Feri We consier the unconstraine optimization or the case o unctions
More informationComposite models 30 and 131: Ply types 0 and 8 calibration
Model calibration Composite Bi-Phase models 30 and 3 for elastic, damage and failure PAM-CRASH material model 30 is for solid and 3 for multi-layered shell elements. Within these models different ply types
More informationMASSACHUSETTS INSTITUTE OF TECHNOLOGY DEPARTMENT OF MATERIALS SCIENCE AND ENGINEERING CAMBRIDGE, MASSACHUSETTS 02139
MASSACHUSTTS INSTITUT OF TCHNOLOGY DPARTMNT OF MATRIALS SCINC AND NGINRING CAMBRIDG, MASSACHUSTTS 0239 322 MCHANICAL PROPRTIS OF MATRIALS PROBLM ST 4 SOLUTIONS Consider a 500 nm thick aluminum ilm on a
More informationThe derivative of a function f(x) is another function, defined in terms of a limiting expression: f(x + δx) f(x)
Y. D. Chong (2016) MH2801: Complex Methos for the Sciences 1. Derivatives The erivative of a function f(x) is another function, efine in terms of a limiting expression: f (x) f (x) lim x δx 0 f(x + δx)
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 informationELECTRON DIFFRACTION
ELECTRON DIFFRACTION Electrons : wave or quanta? Measurement of wavelength an momentum of electrons. Introuction Electrons isplay both wave an particle properties. What is the relationship between the
More informationDOUBLE PENDULUM VIBRATION MOTION IN FLUID FLOW
ENGINEERING FOR RURA DEEOPMENT Jelgava,.-.5.5. DOUBE PENDUUM IBRATION MOTION IN FUID FOW Janis iba, Maris Eiuks, Martins Irbe Riga Technical University, atvia janis.viba@rtu.lv, maris.eiuks@rtu.lv, martins.irbe@rtu.lv
More informationA Simple Model for the Calculation of Plasma Impedance in Atmospheric Radio Frequency Discharges
Plasma Science an Technology, Vol.16, No.1, Oct. 214 A Simple Moel for the Calculation of Plasma Impeance in Atmospheric Raio Frequency Discharges GE Lei ( ) an ZHANG Yuantao ( ) Shanong Provincial Key
More informationNon-Superconducting Fault Current Limiter
Proceeings o the 5th WEA International Conerence on Applications o Electrical Engineering, Prague, Czech Republic, March -, 6 (pp-7) Non-uperconucting Fault Current Limiter A. CHARMIN*, M. TARAFDAR HAQUE,
More information6. Friction and viscosity in gasses
IR2 6. Friction an viscosity in gasses 6.1 Introuction Similar to fluis, also for laminar flowing gases Newtons s friction law hols true (see experiment IR1). Using Newton s law the viscosity of air uner
More informationWork and Kinetic Energy
Work Work an Kinetic Energy Work (W) the prouct of the force eerte on an object an the istance the object moes in the irection of the force (constant force only). W = " = cos" (N " m = J)! is the angle
More informationBasic Strategy for Card-Counters: An Analytic Approach. by N. Richard Werthamer
Basic Strategy or Car-Counters: An Analytic Approach by N. Richar Werthamer Abstract We evelop an eicient metho, similar to that o Marcus (7) who calls it Counter Basic Strategy, or casino blackack car
More informationSources and Sinks of Available Potential Energy in a Moist Atmosphere. Olivier Pauluis 1. Courant Institute of Mathematical Sciences
Sources an Sinks of Available Potential Energy in a Moist Atmosphere Olivier Pauluis 1 Courant Institute of Mathematical Sciences New York University Submitte to the Journal of the Atmospheric Sciences
More informationSeparation of Variables
Physics 342 Lecture 1 Separation of Variables Lecture 1 Physics 342 Quantum Mechanics I Monay, January 25th, 2010 There are three basic mathematical tools we nee, an then we can begin working on the physical
More informationBoth the ASME B and the draft VDI/VDE 2617 have strengths and
Choosing Test Positions for Laser Tracker Evaluation an Future Stanars Development ala Muralikrishnan 1, Daniel Sawyer 1, Christopher lackburn 1, Steven Phillips 1, Craig Shakarji 1, E Morse 2, an Robert
More informationQUESTION BANK Composite Materials
QUESTION BANK Composite Materials 1. Define composite material. 2. What is the need for composite material? 3. Mention important characterits of composite material 4. Give examples for fiber material 5.
More informationProblems Governed by PDE. Shlomo Ta'asan. Carnegie Mellon University. and. Abstract
Pseuo-Time Methos for Constraine Optimization Problems Governe by PDE Shlomo Ta'asan Carnegie Mellon University an Institute for Computer Applications in Science an Engineering Abstract In this paper we
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 informationLQG FLUTTER CONTROL OF WIND TUNNEL MODEL USING PIEZO-CERAMIC ACTUATOR
5 TH INTERNATIONAL CONGRESS OF THE AERONAUTICAL SCIENCES LQG FLUTTER CONTROL OF WIND TUNNEL MODEL USING PIEZO-CERAMIC ACTUATOR Tatsunori Kaneko* an Yasuto Asano* * Department of Mechanical Engineering,
More informationSYNTHESIS ON THE ASSESSMENT OF CHIPS CONTRACTION COEFFICIENT C d
SYNTHESIS ON THE ASSESSMENT OF HIPS ONTRATION OEFFIIENT Marius MILEA Technical University Gheorghe Asachi of Iași, Machine Manufacturing an Inustrial Management Faculty, B-ul D Mangeron, nr 59, A, 700050,
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 informationGravity wave propagation in a nonisothermal atmosphere with height varying background wind
GEOPHYSICAL RESEARCH LETTERS, VOL. 34, L2383, oi:.29/27gl36, 27 Gravity wave propaation in a nonisothermal atmosphere with heiht varyin backroun win Qihou Zhou an Yu T. Morton Receive 2 July 27; revise
More informationSchöck Isokorb type KST
Schöck Isokorb type Schöck Isokorb type Contents Page Element arrangements/connection layouts 288-289 Views/Dimensions 290-293 Design an capacity table 294 Torsion spring strength/notes on calculations
More informationProduct and Quotient Rules and Higher-Order Derivatives. The Product Rule
330_003.q 11/3/0 :3 PM Page 119 SECTION.3 Prouct an Quotient Rules an Higher-Orer Derivatives 119 Section.3 Prouct an Quotient Rules an Higher-Orer Derivatives Fin the erivative o a unction using the Prouct
More informationEvaporating droplets tracking by holographic high speed video in turbulent flow
Evaporating roplets tracking by holographic high spee vieo in turbulent flow Loïc Méès 1*, Thibaut Tronchin 1, Nathalie Grosjean 1, Jean-Louis Marié 1 an Corinne Fournier 1: Laboratoire e Mécanique es
More information