An implicit integration method for solving non-linear problems in mechanics of structures under complex loading histories
|
|
- Samson Lawson
- 6 years ago
- Views:
Transcription
1 An implicit integration method for solving non-linear problems in mechanics of structures under complex loading histories M. Arzt, W. Brocks & R. Mohr GKSS Forschungszentrum Geesthacht GmbH Institutfur Werkstofforschung Max-Planck-Strafie, Geesthacht, Germany Abstract An implicit integration method is discussed. This method allows for studying nonlinear material behaviour, such as viscoplasticity. The material behaviour is described by so-called internal variables. The model applied here was developed by Lemaitre and Chaboche. The integration method is implemented in a usersubroutine in the commercial Finite-Element-Method package ABAQUS. The results of a numerical analysis of the rotation of an aircraft turbine disk are presented. Cyclic loading conditions are studied, and the results are compared to numerical results found in literature. 1 Introduction Structures under cyclic loading conditions are studied. The inelastic behaviour under cyclic loading is an interesting task since the failure of structures is often caused by repeated cycles of loading, e.g. an aircraft turbine disk, Dambrine & Mascarell [7]. To simulate the material behaviour in case of cyclic loading, multiaxial and out-of-phase loading conditions, non-linear constitutive relations have to be applied, Benallal & Ben Cheikh [3] and Ohno [15]. At present, constitutive equations described by internal variables are the common approach to simulate the phenomena of isotropic and kinematic hardening, Nouailhas [13] and Wang & Ohno [17]. The viscoplastic model applied is based on the works of Chaboche [4], [5] and Lemaitre & Chaboche [11]. These constitutive relations are developed to simulate cyclic loading. The model of Chaboche, implemented in the commercial FEM package ABAQUS is rate-independent, plasticity behaviour, Aravas [1].
2 372 Computational Methods and Experimental Measurements The model applied in this paper is rate-dependent, to simulate viscoplastic behaviour. This model is implemented in a user-subroutine of ABAQUS. The unknowns of the problem are defined in the second paragraph. The constitutive relations are strongly non-linear. Therefore, the numerical analysis of structures is complex and the integration methods are difficult to develop. The user-subroutine is based on an implicit integration method, Hornberger & Stamm [8] and Olschewski [14]. An expansion in a Taylor series of the constitutive equations is carried out. An Euler scheme is applied too. The increments of the unknowns are a superposition of the rates at the begin of the increment and at the end of the increment. A weighting factor permits to chose an explicit integration method or an implicit integration method, and the mid-point rule as well. The solution of the linearized system is determined using a public domain solver, the second task. The third task is to determine the so-called Jacobi matrix, or the tangent operator, in order to reduce the number of iterations, and to increase the time increments. The required processing time and memory for storage, to solve these problems, is still important. Even today, there are few complex computations carried out in industry, although the performances of computers is increasing, and strategies are developed to carry out these computations, Arzt et al. [2], Cognard & Ladeveze [6] and Ladeveze & Rougee [9]. These methods allow to reduce the processing time and to decrease the storage capacities required, Ladeveze [10] and Lesne & Savalle[12]. 2 Material Behaviour Here, the unknowns of the problem are defined, stress, back stress, isotropic hardening and the inelastic strain, Chaboche [4], [5] and Lemaitre & Chaboche [11]. The mechanical behaviour is described by so-called internal variables, the equivalent plastic strain, p, and the kinematic hardening variable, a, in order to simulate cyclic phenomena. The expansion or contraction, and translation of the yield surface are represented by the corresponding variables, the isotropic hardening function, R, and the back stress, x, respectively. Isotropic and kinematic hardening are comined, and kinematic hardening is non-linear. The potential of the free energy, where the density is />, is taken as pv = i De*V 4- I caaa + h(p). (1) 2 o The elastic stiffness or Hooke matrix is D. Kinematic hardening constants are c and a. For the isotropic hardening, the function h(p) is chosen as the potential: The asymptotic value of the isotropic hardening function is Q, eqn. (5), and b is an exponent.
3 Computational Methods and Experimental Measurements 373 The state equations are obtained by differentiating the potential of the free energy (1), using (2): a = Be*', (3) 2 x = - caa, (4) The constitutive equations are derived from the dissipation potential: (5) Where n stands for the hardening exponent, and K is the coefficient of resistance. The yield function /(<r, x, R) is defined as: /(<r,x,fl) = J2(<r-x)-/Z-*, (7) and the second invariant or the equivalent von Mises stress is. where cr*** is the deviatoric part of the stress. The yield strength at zero plastic strain is k. A function, the MacCauley bracket, is defined as: (y) = yify>0 and (j/> = 0 if y < 0. Derivation of the dissipation potential leads to the constitutive equations, Lemaitre & Chaboche [11]: R = b(q-r)p. (10) The second term, the relaxation term, of the kinematic hardening is introduced to decrease the effect of the plastic strain. The function of the relaxation term, < (p), is The constants $00 and w are material parameters. The third term describes the static recovery, the parameters are the constant d and the exponent r. This version of the model is a rather simple one. More kinematic hardening variables can be introduced, in order to describe the phenomenon of hysteresis.
4 374 Computational Methods and Experimental Measurements 3 Implicit Integration Method An implicit integration method, Hornberger & Stamm [8], Olschewski [14] and Qi [16], is developed. The constitutive equations, (8), (9) and (10), are taken in the following expression: x = h(cr, x, R) R = r(<r,x, R). (11) (12) The aim is to determine the increments of the inelastic strain, Ae'", the back stress, Ax, and isotropic hardening function, AA, for a given time increment, A(, by the rates at the begin of the increment, t, and the end of the increment, t + At. A superposition, of the rates at the begin and the end of the increment, is applied: (13) Ac'" = (14) The superposition includes the explicit, 0 = 0, and implicit, 0=1, integration method by Euler, and the mid-point rule, 0 = 1/2, as well, since the parameter can be choosen in 0 < 0 < 1. The rates at the end of the increment are unknown. An expansion of the constitutive equations using a Taylor series at the begin of the increment leads to a linearized form: <9g <9g AcrH da t ^ 9x dh c9h X t + At = h t + ^~ ACTt da dr dr ACT 4 da Ax 4- t Ax-f t Ax + t ^g 9h 5r (15) (16) AB (17) (18) (19) The system of equations to be solved leads to 4 unknowns, the increments of the stress, ACT, back stress, Ax, isotropic hardening function, A/2, and the inelastic strain, Ae'*\ To solve the system, the elastic behaviour and the usual decomposition of the total strain are taken = D(Ac- (20) Therefore, the numerical calculations are reduced to solve a system of the size of 13 unknowns. In case of a three dimensional problem, there are 6 components of the stress, 6 components of the back stress, and one of the isotropic hardening. The system to be solved is written in the following form: MAy = (21)
5 Computational Methods and Experimental Measurements 375 The matrix M is of the following structure M = The components of the matrix are: M21 M22 11*23 ni (22) <9g, M. = {l-,a^j gb ah gb oa The right hand side vector is given i - Ae The components are determined by the constitutive equations, (8), (9) and (10), and the rates of the previous time step, applying the superposition, (14), (15) and (16): 63 (23) 63 = The vector of the unknowns is: H- (1 - ACT Ax (24) The system is solved using subroutines developed by Cleve Moler, University of New Mexico, Argonne National Laboratory. Once, the vector of the unknowns is determined, the inelastic strain increment, eqn. (20), is calculated. Then the stress, back stress, isotropic hardening and the inelastic strain are updated. The last task is to define the Jacobi matrix. Therefore, the system, eqn. (21), is taken and transformed to obtain a form like: ACT = KnAe ^n,0- (25) The matrix, K^, is determined for each integration point, n. And then, the matrix, K, of the global problem is assembled in order to calculate the equilibrium solution. The system, (21), is solved analytically. First, the increments of the back stress, Ax, and isotropic hardening function, AH, are determined, since these are independent of the strain increment, Ac. These increments are a function of the increment of the stress, ACT.
6 376 Computational Methods and Experimental Measurements The Jacobi matrix is of the following structure: K-I = -[MH - MizK^ where the matrices are: Kr = m33m22 - m23m32 and M and the vector is Finally, the Jacobi matrix, Kn, is obtained. 4 Numerical Simulation As a numerical simulation is carried out the problem of a rotating aircraft turbine disk, Lesne & Savalle, [12]. The problem is quasi-static and isothermal conditions are assumed. The mechanical loading is the centrifugal force. The problem is axisymmetric. The diameter of the disk is 300 mm, and the diameter of the whole is 75 mm. The thickness of the disk at the bore is 50 mm, and at the neck the thickness is 15 mm. Infig.1 the mesh of the turbine disk is shown. The number of three-node elements is 571, and the number of nodes is 335. The plane of symmetry is the plane 1-3. The total number of degree of freedom is about 629. Figure 1: mesh of the aircraft turbine disk The temperature field is homogeneous, and the temperature is of 550 C. The material is the nickel-based alloy INCONEL 718, the data is given in tab. 1, Benallal & Ben Cheikh, [3]. The data allows for studying cyclic mean-stress relaxation. Table 1: parameters of Chaboche model E = MPa v 0.3 k = 646 MPa Q = -185 MPa c = 500 d = 0 1/s n b a r K MPa MPa p = 8190 kg/nf
7 Computational Methods and Experimental Measurements 377 A loading of 500 cycles is simulated. The shape of thefirstthree loading cycles is presented infig.2. The minimum rotation speed is 1500 r/min, and the maximum rotation speed is r/min. In contrast to Lesne & Savalle, [12], where the maximum rotation speed is r/min. In fact of the possibility to study cyclic mean-stress relaxation, the maximum rotation speed is higher, in order to obtain important inelastic strains time t in [sec] Figure 2: cycles of loading The processing time on an IBM Rise/System 6000 computer is 7 hs, CPU. The integration parameter is chosen equal to 1/2. The results are compared to data determined by Lesne and Savalle, [12]. In fig. 3 the stress versus inelastic strain for a node at the bore and at the neck is presented. An important variation of the stress is noticed. The inelastic strain magnitude is of 0.6 % at the bore. In fact of the mean-stress relaxation, and the higher level of the maximum rotation speed, the stresses and inelastic strains are of the same magnitude. bore, [12] neck, [12] bore, authors neck, authors inelastic strain 33^ in [%] Figure 3: hoop stress versus inelastic hoop strain, at the bore and at the neck
8 378 Computational Methods and Experimental Measurements The hoop stress at the beginning of maximum rotation speed for thefirstand tenth cycle, and for cycles 100 and 500 are given infigs.4, 5, 6 and 7, respectively. The stress distribution versus the number of cycles is observed. The minimum level of the hoop stress decreases, and the maximum level of the hoop stress increases. The distribution is shown, to give an idea of the zones where plastic flow occurs. Remark a high level of the stress at the bore and at the neck. Important stress gradients are observed at the exterior part of the disk. The stable stress distribution is nearly reached at cycle 100. Figure 4: hoop stress at beginning of maximum rotation speed for cycle [MPa] 969 Figure 5: hoop stress at beginning of maximum rotation speed for cycle 10 Figure 6: hoop stress at beginning of maximum rotation speed for cycle 100 [MPa] Figure 7: hoop stress at beginning of maximum rotation speed for cycle 500
9 Computational Methods and Experimental Measurements Conclusions The performances of the method is shown. Numerical results correspond to the results obtained by Lesne & Savalle [12]. Accuracy is sufficient. The aim is to take into account anisotropy, an isothermal conditions and damage as well. The increment size will be controlled by procedures in order to accelerate convergence and decrease processing time. Acknowledgements The authors would like to thank Rainer Sievert for providing experimental data, which allowed us to verify the subroutine. References 1. Aravas, N., On the Numerical Integration of a Class of Pressure-Dependent Plasticity Models, InternationalJournal for Numerical Methods in Engineering, John Wiley & Sons Ltd, vol. 24, pp , Arzt, M., Cognard, J.-Y. & Ladeveze, P., A Large Time Increment Strategy for the Analysis of Viscoplastic Structures Under Complex Loading Histories, Proceedings of the International Seminar on Multiaxial Plasticity, eds. Benallal A., Billardon R. & Marquis D., Laboratoire de Mecanique et Technologic, Ecole Normale Superieure de Cachan, France, pp. 434^460, Benallal, A. & Ben Cheikh, A., Constitutive Equations for Anisothermal Elasto-Viscoplasticity, Constitutive Laws for Engineering Materials: Theory and Applications, eds. Desai, C.S. et al, Elsevier Science Publishing Co Inc, pp , Chaboche, J.-L., Cyclic Viscoplastic Constitutive Equations, Part I: A Thermodynamically Consistent Formulation, Journal of Applied Mechanics, Transactions of the ASME, New York, vol. 60, pp , december Chaboche, J.-L., Cyclic Viscoplastic Constitutive Equations, Part II: Stored Energy Comparison Between Models and Experiments, Journal of Applied Mechanics, Transactions of the ASME, New York, vol. 60, pp , december Cognard, J.-Y. & Ladeveze, P., A Parallel Computer Implementation for Elastoplastic Calculations with the Large Time Increment Method, Nonlinear Engineering Computations, Swansea, pp. 1-10, 1991
10 380 Computational Methods and Experimental Measurements 7. Dambrine, B. & Mascarell, J. P., About the Interest of Using Unified Viscoplastic Models in Engine Hot Components Life Prediction, Proceedings of the International Seminar on High Temperature Fracture Mechanisms and Mechanics, eds. Bensussan, P. & Mascarell, J. P., Mechanical Engineering Publications, London, pp. 1-15, Hornberger, K. & Stamm, T, An Implicit Integration Algorithm with a Projection Method for Viscoplastic Constitutive Equations, International Journal for Numerical Methods in Engineering, John Wiley & Sons Ltd, vol. 28, pp , Ladeveze, P. & Rougee, P., Plasticite et viscoplasticite sous chargement cyclique : propri6tes et calcul du cycle limite, Comptes Rendus de I'Academic des Sciences, Gauthier-Villars, Paris, tome 301, serie II, no. 13, pp , Ladeveze, P., La methode a grand increment de temps pour r analyse de structures a comportement non lineaire decrit par variables intermes, Comptes Rendus de VAcademic des Sciences, Gauthier-Villars, Paris, tome 309, serie II, pp , Lemaitre, J. & Chaboche, J.-L., Mechanics of Solid Materials, Cambridge University Press, UK, Lesne, M. P. & Savalle, S., An Efficient Cycles Jump Technique for Viscoplastic Structure Calculations Involving Large Number of Cycles, Proceedings of the Second International Conference on Computational Plasticity, eds. Owen, D. R. J.; Onate, E. & Hinton E., Pineridge Press Ltd, Swansea, pp , Nouailhas, D., Unified Modelling of Cyclic Viscoplasticity: Application to Austenitic Stainless Steels, International Journal of Plasticity, Pergamon Press pic, vol. 5, pp , Olschewski, J., Plastizitdtstheorie II, Vorlesung an der Technischen Universitat Berlin, BAM /4, Bundesanstalt fiir Materialforschung und -priifung, November Ohno, N., Constitutive Modelling of Cyclic Plasticity with Emphasis on Ratchetting, InternationalJournal of Mechanical Science, Elsevier Science Ltd, vol. 40, nos. 2-3, pp , Qi, W., Entwicklung eines allgemeinen FE-Bausteins fur viskoplastische Materialmodelle auf der Basis der fur thermisch-mechanische Belastungsvorgdnge erweiterten Integrationsmethode nach Tanaka, Diplomarbeit, Technische Universitat Berlin, Institut fur Mechanik, Berlin, Wang, J.-D. & Ohno, N., Two Equivalent Forms of Non-linear Kinematic Hardening: Application to Non-isothermal Plasticity, International Journal of Plasticity, Pergamon Press pic, vol. 7, pp , 1991
MODELING OF ELASTO-PLASTIC MATERIALS IN FINITE ELEMENT METHOD
MODELING OF ELASTO-PLASTIC MATERIALS IN FINITE ELEMENT METHOD Andrzej Skrzat, Rzeszow University of Technology, Powst. Warszawy 8, Rzeszow, Poland Abstract: User-defined material models which can be used
More informationAn Evaluation of Simplified Methods to Compute the Mechanical Steady State
An Evaluation of Simplified Methods to Compute the Mechanical Steady State T. Herbland a,b, G. Cailletaud a, S. Quilici a, H. Jaffal b, M. Afzali b a Mines Paris Paris Tech, CNRS UMR 7633, BP 87, 91003
More informationDEVELOPMENT OF A CONTINUUM PLASTICITY MODEL FOR THE COMMERCIAL FINITE ELEMENT CODE ABAQUS
DEVELOPMENT OF A CONTINUUM PLASTICITY MODEL FOR THE COMMERCIAL FINITE ELEMENT CODE ABAQUS Mohsen Safaei, Wim De Waele Ghent University, Laboratory Soete, Belgium Abstract The present work relates to the
More informationFE FORMULATIONS FOR PLASTICITY
G These slides are designed based on the book: Finite Elements in Plasticity Theory and Practice, D.R.J. Owen and E. Hinton, 1970, Pineridge Press Ltd., Swansea, UK. 1 Course Content: A INTRODUCTION AND
More informationMODELLING OF THE CYCLIC AND VISCOPLASTIC BEHAVIOR OF A COPPER-BASE ALLOY USING CHABOCHE MODEL
SP2016_3124805 MODELLING OF THE CYCLIC AND VISCOPLASTIC BEHAVIOR OF A COPPER-BASE ALLOY USING CHABOCHE MODEL Wissam BOUAJILA (1), Jörg RICCIUS (1) (1) German Aerospace Centre DLR Lampoldshausen, Langer
More informationHERCULES-2 Project. Deliverable: D4.4. TMF model for new cylinder head. <Final> 28 February March 2018
HERCULES-2 Project Fuel Flexible, Near Zero Emissions, Adaptive Performance Marine Engine Deliverable: D4.4 TMF model for new cylinder head Nature of the Deliverable: Due date of the Deliverable:
More informationAPPLICATION OF DAMAGE MODEL FOR NUMERICAL DETERMINATION OF CARRYING CAPACITY OF LARGE ROLLING BEARINGS
INTERNATIONAL ESIGN CONFERENCE - ESIGN ubrovnik, May 14-17,. APPLICATION OF AMAGE MOEL FOR NUMERICAL ETERMINATION OF CARRYING CAPACITY OF LARGE ROLLING BEARINGS Robert Kunc, Ivan Prebil, Tomaž Rodic and
More informationPLASTICITY AND VISCOPLASTICITY UNDER CYCLIC LOADINGS
ATHENS Course MP06 Nonlinear Computational Mechanics March 16 to 20, 2009 PLASTICITY AND VISCOPLASTICITY UNDER CYCLIC LOADINGS Jean-Louis Chaboche ONERA, 29 av. de la Division Leclerc 92320 Châtillon,
More informationMECHANICS OF MATERIALS. EQUATIONS AND THEOREMS
1 MECHANICS OF MATERIALS. EQUATIONS AND THEOREMS Version 2011-01-14 Stress tensor Definition of traction vector (1) Cauchy theorem (2) Equilibrium (3) Invariants (4) (5) (6) or, written in terms of principal
More informationTangent Modulus in Numerical Integration of Constitutive Relations and its Influence on Convergence of N-R Method
Applied and Computational Mechanics 3 (2009) 27 38 Tangent Modulus in Numerical Integration of Constitutive Relations and its Influence on Convergence of N-R Method R. Halama a,, Z. Poruba a a Faculty
More informationAbstract. 1 Introduction
Contact analysis for the modelling of anchors in concrete structures H. Walter*, L. Baillet** & M. Brunet* *Laboratoire de Mecanique des Solides **Laboratoire de Mecanique des Contacts-CNRS UMR 5514 Institut
More informationFinite Element Formulation for Prediction of Over-speed and burst-margin limits in Aero-engine disc
International Journal of Soft Computing and Engineering (IJSCE) Finite Element Formulation for Prediction of Over-speed and burst-margin limits in Aero-engine disc Maruthi B H, M. Venkatarama Reddy, K.
More informationCHOICE AND CALIBRATION OF CYCLIC PLASTICITY MODEL WITH REGARD TO SUBSEQUENT FATIGUE ANALYSIS
Engineering MECHANICS, Vol. 19, 2012, No. 2/3, p. 87 97 87 CHOICE AND CALIBRATION OF CYCLIC PLASTICITY MODEL WITH REGARD TO SUBSEQUENT FATIGUE ANALYSIS Radim Halama*, Michal Šofer*, František Fojtík* Plasticity
More informationComputational Inelasticity FHLN05. Assignment A non-linear elasto-plastic problem
Computational Inelasticity FHLN05 Assignment 2017 A non-linear elasto-plastic problem General instructions A written report should be submitted to the Division of Solid Mechanics no later than October
More informationChapter 2 Finite Element Formulations
Chapter 2 Finite Element Formulations The governing equations for problems solved by the finite element method are typically formulated by partial differential equations in their original form. These are
More informationConstitutive models: Incremental plasticity Drücker s postulate
Constitutive models: Incremental plasticity Drücker s postulate if consistency condition associated plastic law, associated plasticity - plastic flow law associated with the limit (loading) surface Prager
More informationFatigue Damage Development in a Steel Based MMC
Fatigue Damage Development in a Steel Based MMC V. Tvergaard 1,T.O/ rts Pedersen 1 Abstract: The development of fatigue damage in a toolsteel metal matrix discontinuously reinforced with TiC particulates
More informationPredeformation and frequency-dependence : Experiment and FE analysis
Predeformation and frequency-dependence : Experiment and FE analysis Nidhal Jridi 1,2,*, Michelle Salvia 2, Adel Hamdi 1, Olivier Bareille 2, Makrem Arfaoui 1, Mohammed Ichchou 2, Jalel Ben Abdallah 1
More informationThe Influence of Strain Amplitude, Temperature and Frequency on Complex Shear Moduli of Polymer Materials under Kinematic Harmonic Loading
Mechanics and Mechanical Engineering Vol. 21, No. 1 (2017) 157 170 c Lodz University of Technology The Influence of Strain Amplitude, Temperature and Frequency on Complex Shear Moduli of Polymer Materials
More informationCode_Aster. Reaction of the steel subjected to corrosion
Titre : Comportement de l acier soumis à la corrosion Date : 16/0/01 Page : 1/1 Reaction of the steel subjected to corrosion Summarized: This documentation the model presents steel reaction of subjected
More information3D-FE Implementation of Evolutionary Cyclic Plasticity Model for Fully Mechanistic (non S-N curve) Fatigue Life Evaluation
3D-FE Implementation of Evolutionary Cyclic Plasticity Model for Fully Mechanistic (non S-N curve) Fatigue Life Evaluation Bipul Barua, Subhasish Mohanty 1, Joseph T. Listwan, Saurindranath Majumdar, and
More informationSiping Road 1239, , Shanghai, P.R. China
COMPARISON BETWEEN LINEAR AND NON-LINEAR KINEMATIC HARDENING MODELS TO PREDICT THE MULTIAXIAL BAUSCHINGER EFFECT M.A. Meggiolaro 1), J.T.P. Castro 1), H. Wu 2) 1) Department of Mechanical Engineering,
More informationContinuum Mechanics and Theory of Materials
Peter Haupt Continuum Mechanics and Theory of Materials Translated from German by Joan A. Kurth Second Edition With 91 Figures, Springer Contents Introduction 1 1 Kinematics 7 1. 1 Material Bodies / 7
More informationThe Finite Element method for nonlinear structural problems
Non Linear Computational Mechanics Athens MP06/2013 The Finite Element method for nonlinear structural problems Vincent Chiaruttini vincent.chiaruttini@onera.fr Outline Generalities on computational strategies
More informationOnedimensional ratchetting Twodimensional ratchetting Response of unified models Response of 2M1C model Response of polycrystalline models
Ratchetting Onedimensional ratchetting Twodimensional ratchetting Response of unified models Response of 2M1C model Response of polycrystalline models Still an open problem Ratchetting effects Primary
More informationThe Finite Element Method for Solid and Structural Mechanics
The Finite Element Method for Solid and Structural Mechanics Sixth edition O.C. Zienkiewicz, CBE, FRS UNESCO Professor of Numerical Methods in Engineering International Centre for Numerical Methods in
More informationThe Finite Element Method II
[ 1 The Finite Element Method II Non-Linear finite element Use of Constitutive Relations Xinghong LIU Phd student 02.11.2007 [ 2 Finite element equilibrium equations: kinematic variables Displacement Strain-displacement
More informationPREDICTION OF FATIGUE LIFE OF COLD FORGING TOOLS BY FE SIMULATION AND COMPARISON OF APPLICABILITY OF DIFFERENT DAMAGE MODELS
PREDICTION OF FATIGUE LIFE OF COLD FORGING TOOLS BY FE SIMULATION AND COMPARISON OF APPLICABILITY OF DIFFERENT DAMAGE MODELS M. Meidert and C. Walter Thyssen/Krupp Presta AG Liechtenstein FL-9492 Eschen
More informationScienceDirect. Bauschinger effect during unloading of cold-rolled copper alloy sheet and its influence on springback deformation after U-bending
Available online at www.sciencedirect.com ScienceDirect Procedia Engineering 81 (2014 ) 969 974 11th International Conference on Technology of Plasticity, ICTP 2014, 19-24 October 2014, Nagoya Congress
More informationLecture #6: 3D Rate-independent Plasticity (cont.) Pressure-dependent plasticity
Lecture #6: 3D Rate-independent Plasticity (cont.) Pressure-dependent plasticity by Borja Erice and Dirk Mohr ETH Zurich, Department of Mechanical and Process Engineering, Chair of Computational Modeling
More informationCourse in. Geometric nonlinearity. Nonlinear FEM. Computational Mechanics, AAU, Esbjerg
Course in Nonlinear FEM Geometric nonlinearity Nonlinear FEM Outline Lecture 1 Introduction Lecture 2 Geometric nonlinearity Lecture 3 Material nonlinearity Lecture 4 Material nonlinearity it continued
More informationA Simple and Accurate Elastoplastic Model Dependent on the Third Invariant and Applied to a Wide Range of Stress Triaxiality
A Simple and Accurate Elastoplastic Model Dependent on the Third Invariant and Applied to a Wide Range of Stress Triaxiality Lucival Malcher Department of Mechanical Engineering Faculty of Tecnology, University
More informationELASTOPLASTIC STEEL BEAM BENDING ANALYSIS BY USING ABAQUS
11 ELASTOPLASTIC STEEL BEAM BENDING ANALYSIS BY USING ABAQUS Biswajit Jena * * Corresponding author: biswa.tech88@gmail.com Abstract: In recent years tremendous efforts have been made in development of
More informationMODELING OF CONCRETE MATERIALS AND STRUCTURES. Kaspar Willam. Uniaxial Model: Strain-Driven Format of Elastoplasticity
MODELING OF CONCRETE MATERIALS AND STRUCTURES Kaspar Willam University of Colorado at Boulder Class Meeting #3: Elastoplastic Concrete Models Uniaxial Model: Strain-Driven Format of Elastoplasticity Triaxial
More informationMHA042 - Material mechanics: Duggafrågor
MHA042 - Material mechanics: Duggafrågor 1) For a static uniaxial bar problem at isothermal (Θ const.) conditions, state principle of energy conservation (first law of thermodynamics). On the basis of
More informationUniversity of Sheffield The development of finite elements for 3D structural analysis in fire
The development of finite elements for 3D structural analysis in fire Chaoming Yu, I. W. Burgess, Z. Huang, R. J. Plank Department of Civil and Structural Engineering StiFF 05/09/2006 3D composite structures
More informationA numerical parametric study on the plate inclination angle of steel plate shear walls
Journal of Structural Engineering & Applied Mechanics 219 Volume 1 Issue 4 164-173 https://doi.org/1.31462/jseam.219.1164173 www.goldenlightpublish.com RESEARCH ARTICLE A numerical parametric study on
More informationOutline. Tension test of stainless steel type 301. Forming including the TRIP-effect
Outline Erweiterte Kontaktformulierung, User-Reibmodell und thermische Analyse Simulationsmöglichkeiten zur Werkzeugauslegung mit LS-DYNA Workshop Simulation in der Umformtechnik Institut für Statik und
More informationelastoplastic contact problems D. Martin and M.H. Aliabadi Wessex Institute of Technology, Ashurst Lodge, Ashurst, Southampton, SO40 7AA, UK
Non-conforming BEM elastoplastic contact problems D. Martin and M.H. Aliabadi discretisation in Wessex Institute of Technology, Ashurst Lodge, Ashurst, Southampton, SO40 7AA, UK Abstract In this paper,
More informationFinite element analysis of indentation experiments J.M. Olaf Fraunhofer-Insitut fur Werkstoffmechanik, Wohlerstr. 11, D-W Freiburg, Germany
Finite element analysis of indentation experiments J.M. Olaf Fraunhofer-Insitut fur Werkstoffmechanik, Wohlerstr. 11, D-W- 7800 Freiburg, Germany ABSTRACT There are only a few methods suitable for a quantitative
More informationA return mapping algorithm for unified strength theory model
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING Int. J. Numer. Meth. Engng 2015; 104:749 766 Published online 10 June 2015 in Wiley Online Library (wileyonlinelibrary.com)..4937 A return mapping
More informationTransactions on Modelling and Simulation vol 9, 1995 WIT Press, ISSN X
Elastic-plastic model of crack growth under fatigue using the boundary element method M. Scibetta, O. Pensis LTAS Fracture Mechanics, University ofliege, B-4000 Liege, Belgium Abstract Life of mechanic
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 informationEXTENDED ABSTRACT. Dynamic analysis of elastic solids by the finite element method. Vítor Hugo Amaral Carreiro
EXTENDED ABSTRACT Dynamic analysis of elastic solids by the finite element method Vítor Hugo Amaral Carreiro Supervisor: Professor Fernando Manuel Fernandes Simões June 2009 Summary The finite element
More informationFatigue-Ratcheting Study of Pressurized Piping System under Seismic Load
Fatigue-Ratcheting Study of Pressurized Piping System under Seismic Load A. Ravi Kiran, M. K. Agrawal, G. R. Reddy, R. K. Singh, K. K. Vaze, A. K. Ghosh and H. S. Kushwaha Reactor Safety Division, Bhabha
More informationThe Finite Element Method for the Analysis of Non-Linear and Dynamic Systems. Prof. Dr. Eleni Chatzi Lecture ST1-19 November, 2015
The Finite Element Method for the Analysis of Non-Linear and Dynamic Systems Prof. Dr. Eleni Chatzi Lecture ST1-19 November, 2015 Institute of Structural Engineering Method of Finite Elements II 1 Constitutive
More informationABSTRACT INTRODUCTION
Optimization of soil anchorages M. Teschner, C. Mattheck Kernforschungszentrum Karlsruhe GmbH, Institut fur Materialforschung II, W-7500 Karlsruhe 1, Postfach 3640, Germany ABSTRACT A new optimization
More informationOn the Numerical Modelling of Orthotropic Large Strain Elastoplasticity
63 Advances in 63 On the Numerical Modelling of Orthotropic Large Strain Elastoplasticity I. Karsaj, C. Sansour and J. Soric Summary A constitutive model for orthotropic yield function at large strain
More informationLow-Cycle Fatigue Crack Growth in Ti-6242 at Elevated Temperature Rebecka Brommesson 1,a, Magnus Hörnqvist,2,b, and Magnus Ekh 3,c
Low-Cycle Fatigue Crack Growth in Ti-6242 at Elevated Temperature Rebecka Brommesson 1,a, Magnus Hörnqvist,2,b, and Magnus Ekh 3,c 1,3 Department of Applied Mechanics, Chalmers University of Technology,
More informationC.J. Bennett, W. Sun Department of Mechanical, Materials and Manufacturing Engineering, University of Nottingham, Nottingham NG7 2RD, UK
Optimisation of material properties for the modelling of large deformation manufacturing processes using a finite element model of the Gleeble compression test C.J. Bennett, W. Sun Department of Mechanical,
More informationUniaxial and biaxial ratchetting in piping materials experiments and analysis
Uniaxial and biaxial ratchetting in piping materials experiments and analysis S.C. Kulkarni a, Y.M. Desai a, *, T. Kant a, G.R. Reddy b, P. Prasad b, K.K. Vaze b, C. Gupta c a Department of Civil Engineering,
More informationEXPERIMENTAL DETERMINATION OF MATERIAL PARAMETERS USING STABILIZED CYCLE TESTS TO PREDICT THERMAL RATCHETTING
U..B. Sci. Bull., Series D, Vol. 78, Iss. 2, 2016 ISSN 1454-2358 EXERIMENTAL DETERMINATION OF MATERIAL ARAMETERS USING STABILIZED CYCLE TESTS TO REDICT THERMAL RATCHETTING Mohammad ZEHSAZ 1, Farid Vakili
More informationRheology. October 2013
Rheology Georges Cailletaud Centre des Matériaux MINES ParisTech/CNRS October 2013 Georges Cailletaud Rheology 1/44 Contents 1 Mechanical tests Structures Representative material elements 2 Rheological
More informationNumerical Properties of Spherical and Cubical Representative Volume Elements with Different Boundary Conditions
TECHNISCHE MECHANIK, 33, 2, (2013), 97 103 submitted: December 11, 2012 Numerical Properties of Spherical and Cubical Representative Volume Elements with Different Boundary Conditions R. Glüge, M. Weber
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 informationBRIDGES BETWEEN DAMAGE AND FRACTURE MECHANICS
BRIDGES BETWEEN DAMAGE AND FRACTURE MECHANICS Jacky Mazars and Gilles Pijaudier-Cabot,,. Laboratoire de Mecanique et Technologie - Ecole Normale Supeneure 94235 Cachan-France Abstract Fracture mechanics
More informationA Comparative Analysis of Linear and Nonlinear Kinematic Hardening Rules in Computational Elastoplasticity
TECHNISCHE MECHANIK, 32, 2-5, (212), 164 173 submitted: October 2, 211 A Comparative Analysis of Linear and Nonlinear Kinematic Hardening Rules in Computational Elastoplasticity Fabio De Angelis In this
More informationGeometric nonlinear sensitivity analysis for nonparametric shape optimization with non-zero prescribed displacements
0 th World Congress on Structural and Multidisciplinary Optimization May 9-24, 203, Orlando, Florida, USA Geometric nonlinear sensitivity analysis for nonparametric shape optimization with non-zero prescribed
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 informationImpact and Crash Modeling of Composite Structures: A Challenge for Damage Mechanics
Impact and Crash Modeling of Composite Structures: A Challenge for Damage Mechanics Dr. A. Johnson DLR Dr. A. K. Pickett ESI GmbH EURO-PAM 99 Impact and Crash Modelling of Composite Structures: A Challenge
More informationBIAXIAL STRENGTH INVESTIGATION OF CFRP COMPOSITE LAMINATES BY USING CRUCIFORM SPECIMENS
BIAXIAL STRENGTH INVESTIGATION OF CFRP COMPOSITE LAMINATES BY USING CRUCIFORM SPECIMENS H. Kumazawa and T. Takatoya Airframes and Structures Group, Japan Aerospace Exploration Agency 6-13-1, Ohsawa, Mitaka,
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 informationVirtual tests based on model reduction strategies for fatigue analysis
Proceedings of the 7th GACM Colloquium on Computational Mechanics for Young Scientists from Academia and Industry October 11-13, 217 in Stuttgart, Germany Virtual tests based on model reduction strategies
More informationDAMAGE MODEL FOR CONCRETE INCLUDING RESIDUAL HYSTERETIC LOOPS: APPLICATION TO SEISMIC AND DYNAMIC LOADING
Fracture Mechanics of Concrete Structures Proceedings FRAMCOS-3 AEDIFICA TIO Publishers, D-79104 Frei burg, Germany DAMAGE MODEL FOR CONCRETE INCLUDING RESIDUAL HYSTERETIC LOOPS: APPLICATION TO SEISMIC
More informationUSE OF A LINEAR VISCOELASTIC CONSTITUTIVE LAW TO MODEL THE TIME DEPENDENT BEHAVIOUR OF ORTHOTROPIC GLASSFIBRE REINFORCED PLASTIC
USE OF A LINEAR VISCOELASTIC CONSTITUTIVE LAW TO MODEL THE TIME DEPENDENT BEHAVIOUR OF ORTHOTROPIC GLASSFIBRE REINFORCED PLASTIC Roland. M. Hinterhoelzl 1, Guenther Schullerer 1, Kurt Moser 1 1 Institut
More informationComputational Inelasticity FHLN05. Assignment A non-linear elasto-plastic problem
Computational Inelasticity FHLN05 Assignment 2018 A non-linear elasto-plastic problem General instructions A written report should be submitted to the Division of Solid Mechanics no later than November
More informationFinite element prediction of the ultimate axial load capacity of V-section band clamps
Journal of Physics: Conference Series Finite element prediction of the ultimate axial load capacity of V-section band clamps To cite this article: S M Barrans and M Muller 2009 J. Phys.: Conf. Ser. 181
More informationEDEM DISCRETIZATION (Phase II) Normal Direction Structure Idealization Tangential Direction Pore spring Contact spring SPRING TYPES Inner edge Inner d
Institute of Industrial Science, University of Tokyo Bulletin of ERS, No. 48 (5) A TWO-PHASE SIMPLIFIED COLLAPSE ANALYSIS OF RC BUILDINGS PHASE : SPRING NETWORK PHASE Shanthanu RAJASEKHARAN, Muneyoshi
More informationMODELING GEOMATERIALS ACROSS SCALES JOSÉ E. ANDRADE DEPARTMENT OF CIVIL AND ENVIRONMENTAL ENGINEERING EPS SEMINAR SERIES MARCH 2008
MODELING GEOMATERIALS ACROSS SCALES JOSÉ E. ANDRADE DEPARTMENT OF CIVIL AND ENVIRONMENTAL ENGINEERING EPS SEMINAR SERIES MARCH 2008 COLLABORATORS: DR XUXIN TU AND MR KIRK ELLISON THE ROADMAP MOTIVATION
More informationFEM for elastic-plastic problems
FEM for elastic-plastic problems Jerzy Pamin e-mail: JPamin@L5.pk.edu.pl With thanks to: P. Mika, A. Winnicki, A. Wosatko TNO DIANA http://www.tnodiana.com FEAP http://www.ce.berkeley.edu/feap Lecture
More informationCyclic and Tangential Plasticity Effects for the Buckling Behavior of a Thin Wall Pier under Multiaxial and Non-proportional Loading Conditions
Cyclic and Tangential Plasticity Effects for the Buckling Behavior of a Thin Wall Pier under Multiaxial and Non-proportional Loading Conditions MOMII Hideto*, TSUTSUMI Seiichiro** and FINCATO Riccardo***
More informationVORONOI APPLIED ELEMENT METHOD FOR STRUCTURAL ANALYSIS: THEORY AND APPLICATION FOR LINEAR AND NON-LINEAR MATERIALS
The 4 th World Conference on Earthquake Engineering October -7, 008, Beijing, China VORONOI APPLIED ELEMENT METHOD FOR STRUCTURAL ANALYSIS: THEORY AND APPLICATION FOR LINEAR AND NON-LINEAR MATERIALS K.
More informationPrediction of Elastic-Plastic Behaviour of Structures at Notches
Prediction of Elastic-Plastic Behaviour of Structures at Notches TANWEER HUSSAIN*, MUJEEBUDDIN MEMON**, AND ZEESHAN ALI MEMON*** RECEIVED ON 01.05.2012 ACCEPTED ON 21.06.2012 ABSTRACT Under the condition
More informationNUMERICAL SIMULATION OF THE INELASTIC SEISMIC RESPONSE OF RC STRUCTURES WITH ENERGY DISSIPATORS
NUMERICAL SIMULATION OF THE INELASTIC SEISMIC RESPONSE OF RC STRUCTURES WITH ENERGY DISSIPATORS ABSTRACT : P Mata1, AH Barbat1, S Oller1, R Boroschek2 1 Technical University of Catalonia, Civil Engineering
More informationStresses Analysis of Petroleum Pipe Finite Element under Internal Pressure
ISSN : 48-96, Vol. 6, Issue 8, ( Part -4 August 06, pp.3-38 RESEARCH ARTICLE Stresses Analysis of Petroleum Pipe Finite Element under Internal Pressure Dr.Ragbe.M.Abdusslam Eng. Khaled.S.Bagar ABSTRACT
More informationNUMERICAL SIMULATIONS OF CORNERS IN RC FRAMES USING STRUT-AND-TIE METHOD AND CDP MODEL
Numerical simulations of corners in RC frames using Strut-and-Tie Method and CDP model XIII International Conference on Computational Plasticity. Fundamentals and Applications COMPLAS XIII E. Oñate, D.R.J.
More informationFinite Element Method
Finite Element Method Finite Element Method (ENGC 6321) Syllabus Objectives Understand the basic theory of the FEM Know the behaviour and usage of each type of elements covered in this course one dimensional
More informationICM11. Simulation of debonding in Al/epoxy T-peel joints using a potential-based cohesive zone model
Available online at www.sciencedirect.com Procedia Engineering 10 (2011) 1760 1765 ICM11 Simulation of debonding in Al/epoxy T-peel joints using a potential-based cohesive zone model Marco Alfano a,, Franco
More informationnumerical implementation and application for life prediction of rocket combustors Tel: +49 (0)
2nd Workshop on Structural Analsysis of Lightweight Structures. 30 th May 2012, Natters, Austria Continuum damage mechanics with ANSYS USERMAT: numerical implementation and application for life prediction
More informationChaboche Nonlinear Kinematic Hardening Model
Model Sheldon Imaoka Memo Number: STI85A ANSYS Release:.. May 4, 8 Introduction The Chaboche nonlinear kinematic hardening model was added in ANSYS 5.6 to complement the existing isotropic and kinematic
More informationENGN 2290: Plasticity Computational plasticity in Abaqus
ENGN 229: Plasticity Computational plasticity in Abaqus The purpose of these exercises is to build a familiarity with using user-material subroutines (UMATs) in Abaqus/Standard. Abaqus/Standard is a finite-element
More informationRELIABILITY STUDY OF POWER MODULE BY STOCHASTIC UNCERTAINTY METHOD FOR AN AERONAUTICAL APPLICATION
26 TH INTERNATIONAL CONGRESS OF THE AERONAUTICAL SCIENCES RELIABILITY STUDY OF POWER MODULE BY STOCHASTIC UNCERTAINTY METHOD FOR AN A. Micol a, T. Lhommeau b, R. Meuret c, C. Martin a, M. Mermet-Guyennet
More informationA novel approach to predict the growth rate of short cracks under multiaxial loadings
A novel approach to predict the growth rate of short cracks under multiaxial loadings F. Brugier 1&2, S. Pommier 1, R. de Moura Pinho 2, C. Mary 2 and D. Soria 2 1 LMT-Cachan, ENS Cachan / CNRS / UPMC
More informationConstitutive models: Incremental (Hypoelastic) Stress- Strain relations. and
Constitutive models: Incremental (Hypoelastic) Stress- Strain relations Example 5: an incremental relation based on hyperelasticity strain energy density function and 14.11.2007 1 Constitutive models:
More informationIntro to Soil Mechanics: the what, why & how. José E. Andrade, Caltech
Intro to Soil Mechanics: the what, why & how José E. Andrade, Caltech The What? What is Soil Mechanics? erdbaumechanik The application of the laws of mechanics (physics) to soils as engineering materials
More informationPROGRAMMING THE TRANSIENT EXPLICIT FINITE ELEMENT ANALYSIS WITH MATLAB
U.P.B. Sci. Bull., Series D, Vol. 75, Iss. 2, 2013 ISSN 1454-2358 PROGRAMMING THE TRANSIENT EXPLICIT FINITE ELEMENT ANALYSIS WITH MATLAB Andrei Dragoş Mircea SÎRBU 1, László FARKAS 2 Modern research in
More informationNigerian Journal of Technology, Vol. 26, No. 2, June 2007 Edelugo 37
Nigerian Journal of Technology, Vol. 26, No. 2, June 2007 Edelugo 37 APPLICATION OF THE REISSNERS PLATE THEORY IN THE DELAMINATION ANALYSIS OF A THREE-DIMENSIONAL, TIME- DEPENDENT, NON-LINEAR, UNI-DIRECTIONAL
More informationEngineering Solid Mechanics
}} Engineering Solid Mechanics 1 (2013) 1-8 Contents lists available at GrowingScience Engineering Solid Mechanics homepage: www.growingscience.com/esm Impact damage simulation in elastic and viscoelastic
More informationA Krylov Subspace Accelerated Newton Algorithm
A Krylov Subspace Accelerated Newton Algorithm Michael H. Scott 1 and Gregory L. Fenves 2 ABSTRACT In the nonlinear analysis of structural systems, the Newton-Raphson (NR) method is the most common method
More informationKuo-Long LEE, Wen-Fung PAN and Chien-Min HSU
212 Experimental and Theoretical Evaluations of the Effect between Diameter-to-Thickness Ratio and Curvature-Rate on the Stability of Circular Tubes under Cyclic Bending Kuo-Long LEE, Wen-Fung PAN and
More informationComparison of Models for Finite Plasticity
Comparison of Models for Finite Plasticity A numerical study Patrizio Neff and Christian Wieners California Institute of Technology (Universität Darmstadt) Universität Augsburg (Universität Heidelberg)
More information3-D Finite Element Analysis of Instrumented Indentation of Transversely Isotropic Materials
3-D Finite Element Analysis of Instrumented Indentation of Transversely Isotropic Materials Abstract: Talapady S. Bhat and T. A. Venkatesh Department of Material Science and Engineering Stony Brook University,
More informationThe Assessment of Elastic Follow-Up Effects on Cyclic Accumulation of Inelastic Strain Under Displacement- Control Loading
The Assessment of Elastic Follow-Up Effects on Cyclic Accumulation of Inelastic Strain Under Displacement- Control Loading Zabihi Ferezqi, H, Shariati, M & Moud, SH Author post-print (accepted) deposited
More informationF7. Characteristic behavior of solids
F7. Characteristic behavior of solids F7a: Deformation and failure phenomena: Elasticity, inelasticity, creep, fatigue. à Choice of constitutive model: Issues to be considered è Relevance? Physical effect
More informationON THE INTEGRATION OF EQUATIONS OF MOTION: FEM AND MOLECULAR DYNAMICS PROBLEMS
8th International Congress on Computational Mechanics, Volos, 1-15 July 015 ON THE INTEGRATION OF EQUATIONS OF MOTION: FEM AND MOLECULAR DYNAMICS PROBLEMS E.G. Kakouris, V.K. Koumousis Institute of Structural
More informationLecture #7: Basic Notions of Fracture Mechanics Ductile Fracture
Lecture #7: Basic Notions of Fracture Mechanics Ductile Fracture by Dirk Mohr ETH Zurich, Department of Mechanical and Process Engineering, Chair of Computational Modeling of Materials in Manufacturing
More information6 th Pipeline Technology Conference 2011
6 th Pipeline Technology Conference 2011 Mechanical behavior of metal loss in heat affected zone of welded joints in low carbon steel pipes. M. J. Fernández C. 1, J. L. González V. 2, G. Jarvio C. 3, J.
More informationINVERSE METHOD FOR FLOW STRESS PARAMETERS IDENTIFICATION OF TUBE BULGE HYDROFORMING CONSIDERING ANISOTROPY
7 th EUROMECH Solid Mechanics Conference J. Ambrósio et.al. (eds.) Lisbon, Portugal, September 7-11, 2009 INVERSE METHOD FOR FLOW STRESS PARAMETERS IDENTIFICATION OF TUBE BULGE HYDROFORMING CONSIDERING
More informationTolerance Ring Improvement for Reducing Metal Scratch
International Journal of Scientific and Research Publications, Volume 2, Issue 11, November 2012 1 Tolerance Ring Improvement for Reducing Metal Scratch Pattaraweerin Woraratsoontorn*, Pitikhate Sooraksa**
More informationTIME-DEPENDENT BEHAVIOR OF PILE UNDER LATERAL LOAD USING THE BOUNDING SURFACE MODEL
TIME-DEPENDENT BEHAVIOR OF PILE UNDER LATERAL LOAD USING THE BOUNDING SURFACE MODEL Qassun S. Mohammed Shafiqu and Maarib M. Ahmed Al-Sammaraey Department of Civil Engineering, Nahrain University, Iraq
More information