Optimization as a tool for the inverse identification of parameters of nonlinear material models
|
|
- Phebe Webster
- 5 years ago
- Views:
Transcription
1 Optimization as a tool for the inverse identification of parameters of nonlinear material models FILIP HOKES, JIRI KALA, ONDREJ KRNAVEK Institute of Structural Mechanics Brno University of Technology, Faculty of Civil Engineering Veveri 331/95, Brno, CZECH REPUBLIC hokes.f@fce.vutbr.cz, kala.j@fce.vutbr.cz, krnavek.o@fce.vutbr.cz Abstract: - At present, optimization methods are finding application in many technical as well as nontechnical fields. One of the possible applications of optimization procedures is in the identification of formerly unknown parameters of nonlinear constitutive relationships of construction materials which are used in numerical simulations of the behaviour of built structures. In this respect, optimization constitutes a certain counterweight to other soft-computing methods which also include methods of identification based on the exercise of artificial neural networks. The topic of the following contribution is the description of the possible connection of a computational system based on the finite element method with an optimization tool, and an analysis of the identification of the mechanico-physical and fracture mechanical parameters of a material model for concrete based on the Willam-Warnke yield surface from an experimentally-measured loading curve. The individual steps of this inverse identification process are described in individual paragraphs, including a description of the model used, all with the aim of creating a comprehensive procedure which can be utilized in further research or in the practical design of efficient built structures. Key-Words: - FEA, nonlinear material model, inverse identification, optimization, concrete, ANSYS 1 Introduction The numerical solution of tasks in the area of the design of structures has been very popular for many years. This popularity is a result of both the already very good existing knowledge of numerical methods of mathematical modelling and, among other things, the fact that the numerical simulation of certain complex and atypical structural elements in computational systems such as ANSYS [1] or LS- Dyna [2] is significantly cheaper than the experimental testing of real specimens [3,4]. Despite the level of development of current computational tools, efforts are constantly underway to improve them still further, which is most true of all for the area of nonlinear mechanics. Research into the description of the nonlinear behaviour of materials and the testing of the usability of derived constitutive relationships in the calculations needed for specific civil engineering and mechanical engineering tasks is of high importance. With the correct use of an exact model, and when the given prerequisites are fulfilled, the closer approximation of a simulation to the real behaviour of a structure can be achieved with numerical calculations. One result of such an approximation may be the design of still more economical and safer buildings. An example of such a calculation can be found in [5]. With regard to the efforts mentioned earlier, a database of elastoplastic material models, multiplas, has been developed for nonlinear simulations in the ANSYS system. However, the use of this and similar tools results in a problem in real life in the form of the large amount of parameters which need to be known for the selected material model before the launch of the numerical simulation itself. Unfortunately, these parameters, which can be both mechanico-physical and fracture-mechanical in nature, may not be well known in advance. In the above-mentioned situation it is possible to make advantageous use of experimentally obtained load displacement curves measured during the testing of test specimens fabricated from selected materials and reverse identify needed parameters with the aid of inverse identification. The inverse identification process is currently executed by exercising artificial neural networks or with the help of optimization algorithms implemented in academic or commercial optimization software tools. In the case of inverse identification using optimization, the optimization task can be defined as the minimization of the ISBN:
2 difference between the experimental load displacement curve and a curve obtained as the output of a numerical simulation in the selected computational system. In such a task, the design vector is created by previously unknown parameters of the nonlinear material models while the delimiting conditions are derived from the prerequisites from which the used material model was derived, or from the requirement for the smooth convergence of the solution. An example of a similar task can be found in conference contribution [5]. The presented article deals with the abovementioned inverse identification of the parameters of a nonlinear material model from the multiplas library with the help of an optimization algorithm implemented in the optislang programme. These parameters were established from an L-d load displacement curve which was obtained from a three-point bending test on a notched concrete beam. The ANSYS computational system was selected for the numerical simulation of this task. The aim of the article is thus to describe the individual steps for the above-mentioned inverse identification method and provide a comprehensive procedure which is useable in further research and in the design of larger structural units. 2 Definition of the problem The inverse identification of parameters for the Menétrey-Willam material model from the multiplas library was performed on a loaddisplacement curve obtained from an experimental investigation published in [6]. This curve was the output of a three-point bending test carried out on a notched concrete beam with the dimensions 360 x 120 x 58 mm. This is one of several typical tests used to identify the material and fracture parameters of quasi-brittle materials [6]. In order to achieve savings with regard to the computation time needed for the nonlinear simulation, the complexity of the plane stress task was reduced from 3D to 2D. With regard to the simplification mentioned above, boundary conditions were also simplified. In real experiments, support and the application of load both take place using steel cylinders. When the numerical simulations were carried out, only zero vertical displacements were prescribed at the supports. The application of load was also realized in the form of vertical displacement: in order to maintain solvability and obtain the convergence of the solution, horizontal displacement was prevented at the load input location. The final appearance of the idealization of the solved task is depicted in Fig. 1. Fig.1 Idealization of the testing configuration 2.1 Analysis of the input data The cited publication [6] presented tests on four sets of concrete samples of strength class C25/30 according to [7]. Within these four sets, concrete mixtures with slump values of F45 and F70 were combined and then tested at the age of 28 and 170 days. Each set was composed of 5 test specimens. As previously mentioned, the concrete samples were subjected to three-point bending tests, the output of which was a diagram of the dependence of loading force L (measured on a testing machine) and displacement d (measured on the test specimen at midspan). For the purposes of the presented article, only one load-displacement curve was chosen for one of the test specimens with a slump value of F45 and an age of 28 days. The form of this L-d curve is clearly depicted in Fig. 2. A maximum deformation value d max of m was derived from this reference curve. The stated load value enabled the post-critical behaviour of the concrete sample (softening curve) to be expressed within the numerical simulation. Fig.2 Reference L-d curve from a three-point bending test ISBN:
3 3 Description of the process of inverse identification using optimization methods The whole process of inverse identification in the optislang programme was divided into two main parts: sensitivity analysis and global optimization. The individual numerical simulations were dealt with using the Newton-Raphson nonlinear equation solver implemented in ANSYS software. 3.1 Nonlinear numerical analysis A numerical analysis of the three-point bending test was conducted via the finite element method within the ANSYS 15.0 computational system. The nonlinear behaviour of the computational model used in the simulation was determined using the Menétrey-Willam material model from the multiplas library of elasto-plastic material models. Due to task automation requirements, the whole calculation was controlled by a programmed macro which set up the geometry of the computational model and the parameters of the chosen constitutive law, input the solver settings, solved the calculation and finally exported the file containing the points for the L-d diagrams The computational model A computational model of a concrete specimen was put together with previously-defined dimensions from a total of 4800 PLANE182 4-node planar elements with 3 mm long edges and a thickness of 58 mm. At the location of the notch, two parallel lines were created with a shared node at the top of the notch. With regard to the fact that the monitored area with strongly nonlinear behaviour was located above the top of the notch, and with regard to the occurrence of local peaks of tension in the area above the supports, a linear material model was assigned for the elements at both ends of the computational models within a 30 mm-wide area on both sides of the supports The final appearance of the computational model is depicted in Fig The material model A material model from the multiplas library called the Menétrey-Willam model was selected for the numerical solution of a three-point bending test simulation in the ANSYS system. This material model belongs to a group of material models of concrete which do not consider the influence of the strain rate on stress and in which the decomposition of the total plastic strain vector, ε tot, into an elastic component, ε el, and a plastic part, ε pl, is thus expected [8]. ε = ε + ε (1) tot el pl From the perspective of plastic flow, the model used ranks among the group of models with nonassociated plastic flow. The Menétrey-Willam material model [9] is based on the Willam-Warnke yield surface [10], which is, in contrast with the Drucker-Prager surface, a function not only of the first and second but also of the third invariant of the stress deviator (known as the Lode angle ). This modification achieves the softening of the corners of deviatoric planes of the yield surface which (in addition) do not lie at a constant distance from the hydrostatic axis in Haigh-Westergaard space. This aspect is the basic difference in relation to the modified Drucker- Prager model, and it thus expands the possibilities for the utilization of the plasticity model for other types of loading conditions. From the point of view of the use of FEM the chosen material model utilizes the smeared crack concept [11]. The given problem was solved with the aid of a softening function based on the dissipation of specific fracture energy G ft, which thus acts as one of the sought parameters. With regard to the need to remove the negative dependence of the solution on the size of the mesh of finite elements, the nonlinear Menétrey-Willam model makes use of Bažant s Crack Band concept. [12]. In order for this material model to function correctly, a further 11 parameters had to be defined in addition to the above-mentioned specific fracture energy. Basic description of these parameters can be found in Tab. 1. Fig.3 The computational model ISBN:
4 Table 1 Parameters of the selected material model Par. Unit Description E [Pa] Young s modulus of elasticity ν [-] Poisson s ratio f c [Pa] Uniaxial compression strength f t [Pa] Uniaxial tension strength k [-] Ratio between biaxial compressive strength and uniaxial compressive strength ψ [ ] Dilatancy angle (friction angle) ε ml [-] Plastic strain corresponding to the maximum load G fc [Nm/m 2 Specific fracture energy in ] compression Relative stress level at the start of Ω ci [-] nonlinear hardening in compression Ω cr [-] Residual relative stress level in compression G ft [Nm/m 2 ] Specific fracture energy in tension Ω tr [-] Residual relative stress level in tension Within the sensitivity analysis, a total of 300 random realizations of design vectors were executed. Sufficient coverage of the design space was ensured using the ALHS method [14]. The form of the boundary curves and their position with regard to the reference L-d curve is depicted in Fig. 4 along with sufficient coverage of the design space via random realizations. Sensitivity analysis proved that the following material model parameters were most used to find the corresponding form of the L-d curve: Young s elasticity modulus E, Poisson s ratio ν, specific fracture energy in tension G ft, and the relative value of residual tensile strength Ω tr. Another interim result of the sensitivity analysis was the first generation of design vectors used for the subsequent optimization process. These realizations were selected on the basis of the calculated value of the target function. 3.2 Inverse identification The inverse analysis itself was performed in the environment of the optislang optimization system in two stages. The aim of the sensitivity analysis was to reveal the level of sensitivity of the individual parameters to the resultant form of the L- d curve, and thus limit the number of parameters of the material model to the necessary minimum and possibly modify the intervals of the values of these parameters. Within the second stage, the optimization itself was performed with a reduced vector of the designed parameters using a suitable optimization algorithm Sensitivity analysis Sensitivity analysis is basically a task which seeks the level to which output data uncertainties are influenced by the variability of input data [13]. Before the beginning of the sensitivity analysis itself, a set of purely empirical testing simulations was peformed whose aim was to obtain boundary curves. The parameters which correspond to these curves were then used as limit boundaries of the design intervals for the sensitivity analysis. The secondary aim of these testing calculations was to prove the solvability of the calculations for the boundary values of parameters in the computational system, which also implied their solvability within the selected intervals as a result. Fig.4 Realizations of L-d curves in sensitivity analysis Global optimization According to [15], the term optimization can be defined as an effort to obtain the best possible result under a given set of conditions. Optimization is used for the minimization of needed effort or the maximization of an effect. The described inverse identification task can also be defined in the stated intentions. This is namely the effort to obtain an L-d curve from nonlinear numerical simulations that would best correspond to the sample curve obtained from the experiment. In other words, the aim of the executed optimization is to minimize the difference between the reference and calculated curve by varying the input parameters of the constitutive relation. The global optimization process was carried out using an algorithm which ranks among what are known as population methods. These are stochastic ISBN:
5 optimization methods based on the principle of the imitation of biological processes [16]. Two types of population method were implemented in the optislang program, (a) genetic algorithms (GA) and (b) evolution strategies (ES); a combination of these methods could be used within the given software for the solution of the described problem. The objective function for the given optimization problem had the following form: n i= 1 * i y i 2 ERROR = ( y ), (2) where y i * represented the value of the force originating from the numerical simulation and y i was the value of the force obtained from an experimentally obtained L-d curve. With regard to the results of the previous sensitivity analysis, the reduced design vector in the form given below was selected for the optimization task in order to reduce the time required for the calculation: red { E, ν, f, ψ,, Ω, G } T X = ε (3) t ml tr The initial population of design vectors for the optimization was composed of the 10 best sensitivity analysis realizations. In order to find an optimum solution, a total of 25 generations of the design vectors were created and thus a total of 250 numerical simulations were carried out in the ANSYS system. An optimum result was achieved with a final objective function value ERROR = A comparison of the curve for the final values of the parameters and the reference L-d curve is shown in Fig. 5. More accurate results can be achieved with more generations of design vectors and further improvement of the results can be increased by implementing a local optimization process, which is included in the future plans of the authors. Fig.5 L-d curve for optimized parameters ft 4 Results The form of the resultant curve for the identified parameters of the material model has already been presented in Fig. 5. The sizes of the individual values of the parameters of the material model are documented further in Tab. 2. Apart from the final values, the table also shows the values of the preoptimized parameters which originate from the sensitivity analysis, including the value of the target function for this set of parameters. Table 2 Results final values of parameters Par. Units Sensitivity Global optimization ERROR [N 2 ] E [Pa] ν [-] f c [Pa] f t [Pa] k [-] ψ [ ] 10 9 ε ml [-] G fc [Nm/m 2 ] Ω ci [-] Ω cr [-] G ft [Nm/m 2 ] Ω tr [-] Conclusion The results stated in the previous paragraphs document the usability of optimization methods in the area of the inverse identification of the parameters of nonlinear material models. It became apparent, though, that when solving the calculations for a specific task such as a simulated three-point bending test, the only parameters that can successfully be defined are those which are utilized in the given loading method. In order to completely obtain all the parameters of the material model, parameters also need to be identified from experimental data originating from different tests. An undeniable advantage of the described inverse identification method is the fact that every numerically obtained curve for a certain set of parameters represents a truly converged numerical solution. This demonstrates the functionality and wide range of applications of the selected material model. A possible refinement of the optimum can be achieved by increasing the number of design vector generations or via the use of other, different ISBN:
6 optimization procedures, the particularities involved in whose application options are currently the focus of attention of the authors of this contribution. 6 Acknowledgement This contribution was created with the financial aid of project GACR S Aspects of the use of complex nonlinear material models provided by the Czech Science Foundation and under the project No. LO1408 "AdMaS UP - Advanced Materials, Structures and Technologies", supported by Ministry of Education, Youth and Sports under the National Sustainability Programme I".. References: [1] ANSYS, Inc.: ANSYS Mechanical Theory Reference, Release 15.0, [2] Livermore Software Technology Corporation (LSTC). LS-DYNA Theory Manual, [3] F. Hokes, Different Approaches to Numerical Simulations of Prestressed Concrete Structural Elements, Applied Research in Materials and Mechanics Engineering, Vol. 621, 2014, ISSN: , pp [4] P. Hradil, J. Kala, Analysis of the Shear Failure of a Reinforced Concrete Wall, Applied Research in Materials and Mechanics Engineering, Vol. 621, 2014, ISSN: , pp [5] T. Most, Identification of the parameters of complex constitutive models: Least squares minimization vs. Bayesian updating, Reliability Conference in München, [6] A. Strauss, T. Zimmermann, D. Lehký, D. Novák, Z. Keršner, Stochastic fracturemechanical parameters for the performance based design of concrete structures. Structural Concrete, Vol. 15, 2014, pp [7] European Committee for Standardization: Concrete Part 1: Specification, performance, production and conformity, EN 206-1, [8] Dynardo, GmbH.: Multiplas, User s Manual Release for ANSYS 15.0, Weimar, [9] P. G. Menétrey, Numerical Analysis of Punching Failure in Reinforced Concrete Structures, PhD Thesis, École polytechnique fédérale de Lausane EPFL, Lausanne, [10] K. J. Willam, E. P. Warnke, Constitutive Models of Triaxial Behavior of Concrete. In Proceedings of the International Association for Bridge and Structural Engineering, Bergamo, vol. 19, pp [11] R. Pölling, Eine praxisnahe, schädigungsorientierte Materialbeschreibung von Stahlbeton für Strukturanalysen. PhD Thesis, Ruhr-Universität Bochum, Bochum, [12] Z. P. Bazant, B. H. Oh, Crack band theory for fracture of concrete. Material and Structures, Rilem, vol. 16, 1983, pp [13] A. Saltelli, M. Ratto, T. Andres, F. Campolongo, J. Cariboni, D. Gatelli, M. Saisana, S. Tarantola, Global Sensitivity Analysis: The primer. Wiley, London, [14] D. E. Huntington, C. S. Lyrintzis, Improvement to limitations of Latin hypercube sampling. Probabilistic Engineering Mechanics, Elsevier, Amsterdam, vol. 13, 1998, pp [15] S. S. Rao, Engineering Optimization: Theory and Practice, 4th edn. Wiley, London, [16] Dynardo GmbH.: Method for multi-disciplinary optimization and robustness analysis. Germany, [17] Kala, J., Hradil, P., Bajer, M. Reinforced concrete wall under shear load Experimental and nonlinear simulation, International Journal of Mechanics, 9, pp , ISBN:
Numerical Simulation of Reinforced Concrete Beam with Utilization of Elasto-plastic Material Model of Concrete
Numerical Simulation of Reinforced Concrete Beam with Utilization of Elasto-plastic Material Model of Concrete FILIP HOKEŠ, MARTIN HUŠEK, PETR KRÁL, JIŘÍ KALA Faculty of Civil Engineering, Institute of
More informationAvailable online at ScienceDirect. Procedia Engineering 161 (2016 )
Available online at www.sciencedirect.com ScienceDirect Procedia Engineering 161 (2016 ) 892 897 World Multidisciplinary Civil Engineering-Architecture-Urban Planning Symposium 2016, WMCAUS 2016 Parameter
More informationInverse identification of the material parameters of a nonlinear concrete constitutive model based on the triaxial compression strength testing
Focussed on Modelling in Mechanics Inverse identification of the material parameters of a nonlinear concrete constitutive model based on the triaxial compression strength testing Petr Král, Petr Hradil,
More informationAiMT Advances in Military Technology
AiMT Advances in Military Technology Vol. 13, No. 1 (2018), pp. 33-45 ISSN 1802-2308, eissn 2533-4123 DOI 10.3849/aimt.01213 Identification of Concrete Material Model Parameters P. Král 1*, M. Hušek 1,
More informationUncertainty modelling using software FReET
Uncertainty modelling using software FReET D. Novak, M. Vorechovsky, R. Rusina Brno University of Technology Brno, Czech Republic 1/30 Outline Introduction Methods and main features Software FReET Selected
More informationVeveri 331/95, Brno, Czech Republic
Advanced Materials Research Online: 2014-06-24 ISSN: 1662-8985, Vol. 969, pp 298-301 doi:10.4028/www.scientific.net/amr.969.298 2014 Trans Tech Publications, Switzerland Probabilistic Simulation of the
More informationTHE BEHAVIOUR OF REINFORCED CONCRETE AS DEPICTED IN FINITE ELEMENT ANALYSIS.
THE BEHAVIOUR OF REINFORCED CONCRETE AS DEPICTED IN FINITE ELEMENT ANALYSIS. THE CASE OF A TERRACE UNIT. John N Karadelis 1. INTRODUCTION. Aim to replicate the behaviour of reinforced concrete in a multi-scale
More informationALGORITHM FOR NON-PROPORTIONAL LOADING IN SEQUENTIALLY LINEAR ANALYSIS
9th International Conference on Fracture Mechanics of Concrete and Concrete Structures FraMCoS-9 Chenjie Yu, P.C.J. Hoogenboom and J.G. Rots DOI 10.21012/FC9.288 ALGORITHM FOR NON-PROPORTIONAL LOADING
More informationNonlinear FE Analysis of Reinforced Concrete Structures Using a Tresca-Type Yield Surface
Transaction A: Civil Engineering Vol. 16, No. 6, pp. 512{519 c Sharif University of Technology, December 2009 Research Note Nonlinear FE Analysis of Reinforced Concrete Structures Using a Tresca-Type Yield
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 informationParameter identification of damage parameters of LS-DYNA GURSON material model from a tensile test. Lectures. Johannes Will
Lectures Parameter identification of damage parameters of LS-DYNA GURSON material model from a tensile test Johannes Will presented at the Weimar Optimization and Stochastic Days 2009 Source: www.dynardo.de/en/library
More information6. NON-LINEAR PSEUDO-STATIC ANALYSIS OF ADOBE WALLS
6. NON-LINEAR PSEUDO-STATIC ANALYSIS OF ADOBE WALLS Blondet et al. [25] carried out a cyclic test on an adobe wall to reproduce its seismic response and damage pattern under in-plane loads. The displacement
More informationSENSITIVITY ANALYSIS OF LATERAL BUCKLING STABILITY PROBLEMS OF HOT-ROLLED STEEL BEAMS
2005/2 PAGES 9 14 RECEIVED 18.10.2004 ACCEPTED 18.4.2005 Z. KALA, J. KALA SENSITIVITY ANALYSIS OF LATERAL BUCKLING STABILITY PROBLEMS OF HOT-ROLLED STEEL BEAMS ABSTRACT Doc. Ing. Zdeněk KALA, PhD. Brno
More informationCard Variable MID RO E PR ECC QH0 FT FC. Type A8 F F F F F F F. Default none none none 0.2 AUTO 0.3 none none
Note: This is an extended description of MAT_273 input provided by Peter Grassl It contains additional guidance on the choice of input parameters beyond the description in the official LS-DYNA manual Last
More informationLecture #8: Ductile Fracture (Theory & Experiments)
Lecture #8: Ductile Fracture (Theory & Experiments) by Dirk Mohr ETH Zurich, Department of Mechanical and Process Engineering, Chair of Computational Modeling of Materials in Manufacturing 2015 1 1 1 Ductile
More informationNonlinear Analysis of Reinforced Concrete Shells Subjected to Impact Loads
Transactions of the 7 th International Conference on Structural Mechanics in Reactor Technology (SMiRT 7) Prague, Czech Republic, August 7, 00 Paper # J0- Nonlinear Analysis of Reinforced Concrete Shells
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 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 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 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 informationEQUIVALENT FRACTURE ENERGY CONCEPT FOR DYNAMIC RESPONSE ANALYSIS OF PROTOTYPE RC GIRDERS
EQUIVALENT FRACTURE ENERGY CONCEPT FOR DYNAMIC RESPONSE ANALYSIS OF PROTOTYPE RC GIRDERS Abdul Qadir Bhatti 1, Norimitsu Kishi 2 and Khaliq U Rehman Shad 3 1 Assistant Professor, Dept. of Structural Engineering,
More information5 ADVANCED FRACTURE MODELS
Essentially, all models are wrong, but some are useful George E.P. Box, (Box and Draper, 1987) 5 ADVANCED FRACTURE MODELS In the previous chapter it was shown that the MOR parameter cannot be relied upon
More informationSIMPLIFIED CONCRETE MODELING WITH *MAT_CONCRET_DAMAGE_REL3
SIMPLIFIED CONCRETE MODELING WITH *MAT_CONCRET_DAMAGE_REL3 Leonard E Schwer Schwer Engineering & Consulting Services, Windsor CA, USA and L. Javier Malvar Karagozian & Case Structural Engineers, Burbank
More information3D Finite Element analysis of stud anchors with large head and embedment depth
3D Finite Element analysis of stud anchors with large head and embedment depth G. Periškić, J. Ožbolt & R. Eligehausen Institute for Construction Materials, University of Stuttgart, Stuttgart, Germany
More information1 Introduction. Abstract
Abstract This paper presents a three-dimensional numerical model for analysing via finite element method (FEM) the mechanized tunneling in urban areas. The numerical model is meant to represent the typical
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 informationAnálisis Computacional del Comportamiento de Falla de Hormigón Reforzado con Fibras Metálicas
San Miguel de Tucuman, Argentina September 14 th, 2011 Seminary on Análisis Computacional del Comportamiento de Falla de Hormigón Reforzado con Fibras Metálicas Antonio Caggiano 1, Guillermo Etse 2, Enzo
More informationFracture Mechanics of Non-Shear Reinforced R/C Beams
Irina Kerelezova Thomas Hansen M. P. Nielsen Fracture Mechanics of Non-Shear Reinforced R/C Beams DANMARKS TEKNISKE UNIVERSITET Report BYG DTU R-154 27 ISSN 161-2917 ISBN 97887-7877-226-5 Fracture Mechanics
More informationNumerical Modelling of Blockwork Prisms Tested in Compression Using Finite Element Method with Interface Behaviour
13 th International Brick and Block Masonry Conference Amsterdam, July 4-7, 2004 Numerical Modelling of Blockwork Prisms Tested in Compression Using Finite Element Method with Interface Behaviour H. R.
More informationAPPLICATION OF ACOUSTIC EMISSION METHOD DURING CYCLIC LOADING OF CONCRETE BEAM
More info about this article: http://www.ndt.net/?id=21866 Abstract IX th NDT in PROGRESS October 9 11, 2017, Prague, Czech Republic APPLICATION OF ACOUSTIC EMISSION METHOD DURING CYCLIC LOADING OF CONCRETE
More information*MAT_PAPER and *MAT_COHESIVE_PAPER: Two New Models for Paperboard Materials
14 th International LS-DYNA Users Conference Session: Constitutive Modeling *MAT_PAPER and *MAT_COHESIVE_PAPER: Two New Models for Paperboard Materials Jesper Karlsson 1, Mikael Schill 1, Johan Tryding
More informationDetermination of statistical material parameters of concrete using fracture test and inverse analysis based on FraMePID-3PB tool
Determination of statistical material parameters of concrete using fracture test and inverse analysis based on FraMePID-3PB tool David Lehký, Zbyněk Keršner, Drahomír Novák Brno University of Technology,
More informationAuthor(s) Okajima, Kenji; Tanaka, Tadatsugu; Symposium on Backwards Problem in G.
Title Backwards Analysis for Retaining Wa based upon ateral Wall Displacemen Author(s) Okajima, Kenji; Tanaka, Tadatsugu; Proceeding of TC302 Symposium Osaka Citation Symposium on Backwards Problem in
More informationEnhancing Prediction Accuracy In Sift Theory
18 TH INTERNATIONAL CONFERENCE ON COMPOSITE MATERIALS Enhancing Prediction Accuracy In Sift Theory J. Wang 1 *, W. K. Chiu 1 Defence Science and Technology Organisation, Fishermans Bend, Australia, Department
More informationCritical applied stresses for a crack initiation from a sharp V-notch
Focussed on: Fracture and Structural Integrity related Issues Critical applied stresses for a crack initiation from a sharp V-notch L. Náhlík, P. Hutař Institute of Physics of Materials, Academy of Sciences
More informationStochastic Simulation of Aircraft Fuselage Assembly Considering Manufacturing Uncertainties
9th European LS-DYNA Conference 0 Stochastic Simulation of Aircraft Fuselage Assembly Considering Manufacturing Uncertainties Dietmar C. Vogt, Sönke Klostermann EADS Innovation Works, Germany Introduction
More informationTHE S-S CURVE APPROXIMATION FOR GFRP TENSILE SPECIMEN BY USING INVERSE METHOD AND OPTIMIZATION
THE S-S CURVE APPROXIMATION FOR GFRP TENSILE SPECIMEN BY USING INVERSE METHOD AND OPTIMIZATION D. S. Shin and E. S. Jeon Department of Mechanical Engineering, Kongju National University, Korea E-Mail:
More informationMODELING OF THE WEDGE SPLITTING TEST USING AN EXTENDED CRACKED HINGE MODEL
Engineering MECHANICS, Vol. 21, 2014, No. 1, p. 67 72 67 MODELING OF THE WEDGE SPLITTING TEST USING AN EXTENDED CRACKED HINGE MODEL Tomáš Pail, Petr Frantík* The present paper describes a semi-analytical
More informationNon-linear and time-dependent material models in Mentat & MARC. Tutorial with Background and Exercises
Non-linear and time-dependent material models in Mentat & MARC Tutorial with Background and Exercises Eindhoven University of Technology Department of Mechanical Engineering Piet Schreurs July 7, 2009
More informationMODELING OF CONCRETE MATERIALS AND STRUCTURES. Kaspar Willam
MODELING OF CONCRETE MATERIALS AND STRUCTURES Class Meeting #1: Fundamentals Kaspar Willam University of Colorado at Boulder Notation: Direct and indicial tensor formulations Fundamentals: Stress and Strain
More informationChapter. Materials. 1.1 Notations Used in This Chapter
Chapter 1 Materials 1.1 Notations Used in This Chapter A Area of concrete cross-section C s Constant depending on the type of curing C t Creep coefficient (C t = ε sp /ε i ) C u Ultimate creep coefficient
More informationMicroplane Model formulation ANSYS, Inc. All rights reserved. 1 ANSYS, Inc. Proprietary
Microplane Model formulation 2010 ANSYS, Inc. All rights reserved. 1 ANSYS, Inc. Proprietary Table of Content Engineering relevance Theory Material model input in ANSYS Difference with current concrete
More informationDynamic Analysis of a Reinforced Concrete Structure Using Plasticity and Interface Damage Models
Dynamic Analysis of a Reinforced Concrete Structure Using Plasticity and Interface Damage Models I. Rhee, K.J. Willam, B.P. Shing, University of Colorado at Boulder ABSTRACT: This paper examines the global
More informationMASONRY MICRO-MODELLING ADOPTING A DISCONTINUOUS FRAMEWORK
MASONRY MICRO-MODELLING ADOPTING A DISCONTINUOUS FRAMEWORK J. Pina-Henriques and Paulo B. Lourenço School of Engineering, University of Minho, Guimarães, Portugal Abstract Several continuous and discontinuous
More information/tvsb Transactions of the VŠB Technical University of Ostrava Civil Engineering Series,Vol. 16, No.
10.1515/tvsb-2016-0020 Transactions of the VŠB Technical University of Ostrava Civil Engineering Series,Vol. 16, No. 2, 2016 paper #20 Josef KVĚTOŇ 1, Jan ELIÁŠ 2 DISCRETE MODEL FOR CONCRETE FRACTURE:
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 informationDesign issues of thermal induced effects and temperature dependent material properties in Abaqus
Materials Characterisation VII 343 Design issues of thermal induced effects and temperature dependent material properties in Abaqus I. Both 1, F. Wald 1 & R. Zaharia 2 1 Department of Steel and Timber
More informationGlobal Sensitivity Analysis in Industrial Application with LS-OPT
9. LS-DYNA Forum, Bamberg 2010 Optimierung II Global Sensitivity Analysis in Industrial Application with LS-OPT Björn Hohage*, Anja Förderer**, Gordon Geißler**, Heiner Müllerschön** *Audi AG, Ingolstadt,
More informationLoading capacity of yielding connections used in steel arch roadway supports
Ground Support 2013 Y. Potvin and B. Brady (eds) 2013 Australian Centre for Geomechanics, Perth, ISBN 978-0-9806154-7-0 https://papers.acg.uwa.edu.au/p/1304_31_horyl/ Loading capacity of yielding connections
More informationFRACTURE IN HIGH PERFORMANCE FIBRE REINFORCED CONCRETE PAVEMENT MATERIALS
FRACTURE IN HIGH PERFORMANCE FIBRE REINFORCED CONCRETE PAVEMENT MATERIALS ERIK DENNEMAN A thesis submitted in partial fulfilment of the requirements for the degree of PHILOSOPHIAE DOCTOR (ENGINEERING)
More informationAn orthotropic damage model for crash simulation of composites
High Performance Structures and Materials III 511 An orthotropic damage model for crash simulation of composites W. Wang 1, F. H. M. Swartjes 1 & M. D. Gan 1 BU Automotive Centre of Lightweight Structures
More informationLimit analysis of brick masonry shear walls with openings under later loads by rigid block modeling
Limit analysis of brick masonry shear walls with openings under later loads by rigid block modeling F. Portioli, L. Cascini, R. Landolfo University of Naples Federico II, Italy P. Foraboschi IUAV University,
More informationGEO E1050 Finite Element Method Mohr-Coulomb and other constitutive models. Wojciech Sołowski
GEO E050 Finite Element Method Mohr-Coulomb and other constitutive models Wojciech Sołowski To learn today. Reminder elasticity 2. Elastic perfectly plastic theory: concept 3. Specific elastic-perfectly
More informationMATERIAL MECHANICS, SE2126 COMPUTER LAB 2 PLASTICITY
MATERIAL MECHANICS, SE2126 COMPUTER LAB 2 PLASTICITY PART A INTEGRATED CIRCUIT An integrated circuit can be thought of as a very complex maze of electronic components and metallic connectors. These connectors
More informationCDPM2: A damage-plasticity approach to modelling the failure of concrete
CDPM2: A damage-plasticity approach to modelling the failure of concrete Peter Grassl 1, Dimitrios Xenos 1, Ulrika Nyström 2, Rasmus Rempling 2, Kent Gylltoft 2 arxiv:1307.6998v1 [cond-mat.mtrl-sci] 26
More informationUSER S MANUAL 1D Seismic Site Response Analysis Example University of California: San Diego August 30, 2017
USER S MANUAL 1D Seismic Site Response Analysis Example http://www.soilquake.net/ucsdsoilmodels/ University of California: San Diego August 30, 2017 Table of Contents USER'S MANUAL TABLE OF CONTENTS Page
More informationAnisotropic modeling of short fibers reinforced thermoplastics materials with LS-DYNA
Anisotropic modeling of short fibers reinforced thermoplastics materials with LS-DYNA Alexandre Hatt 1 1 Faurecia Automotive Seating, Simplified Limited Liability Company 1 Abstract / Summary Polymer thermoplastics
More informationCHAPTER 3 EXPERIMENTAL STUDY
Experimental Study 42 CHAPTER 3 EXPERIMENTAL STUDY 3.1. INTRODUCTION The experimental study that has been carried out in this thesis has two main objectives: 1. Characterise the concrete behaviour in mode
More informationNumerical modeling of standard rock mechanics laboratory tests using a finite/discrete element approach
Numerical modeling of standard rock mechanics laboratory tests using a finite/discrete element approach S. Stefanizzi GEODATA SpA, Turin, Italy G. Barla Department of Structural and Geotechnical Engineering,
More informationDiscrete Element Modelling of a Reinforced Concrete Structure
Discrete Element Modelling of a Reinforced Concrete Structure S. Hentz, L. Daudeville, F.-V. Donzé Laboratoire Sols, Solides, Structures, Domaine Universitaire, BP 38041 Grenoble Cedex 9 France sebastian.hentz@inpg.fr
More informationCracking in Quasi-Brittle Materials Using Isotropic Damage Mechanics
Cracking in Quasi-Brittle Materials Using Isotropic Damage Mechanics Tobias Gasch, PhD Student Co-author: Prof. Anders Ansell Comsol Conference 2016 Munich 2016-10-12 Contents Introduction Isotropic damage
More informationNonlinear Finite Element Modeling of Nano- Indentation Group Members: Shuaifang Zhang, Kangning Su. ME 563: Nonlinear Finite Element Analysis.
ME 563: Nonlinear Finite Element Analysis Spring 2016 Nonlinear Finite Element Modeling of Nano- Indentation Group Members: Shuaifang Zhang, Kangning Su Department of Mechanical and Nuclear Engineering,
More informationModelling of ductile failure in metal forming
Modelling of ductile failure in metal forming H.H. Wisselink, J. Huetink Materials Innovation Institute (M2i) / University of Twente, Enschede, The Netherlands Summary: Damage and fracture are important
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 informationA FINITE ELEMENT MODEL FOR SIZE EFFECT AND HETEROGENEITY IN CONCRETE STRUCTURES
A FINITE ELEMENT MODEL FOR SIZE EFFECT AND HETEROGENEITY IN CONCRETE STRUCTURES Roque Luiz Pitangueira 1 and Raul Rosas e Silva 2 1 Department of Structural Engineering -Federal University of Minas Gerais
More informationOPTIMAL SENSOR LOCATION FOR PARAMETER IDENTIFICATION IN SOFT CLAY
2 th International Conference on the Application of Computer Science and Mathematics in Architecture and Civil Engineering K. Gürlebeck and T. Lahmer (eds.) Weimar, Germany, 2-22 July 215 OPTIMAL SENSOR
More informationThe design of a formula student front impact attenuator
The design of a formula student front impact attenuator J.M.J. Schormans MT 10.06 Supervisors: dr. ir. Varvara Kouznetsova Eindhoven University of Technology Department of Mechanical Engineering Computational
More informationUSER S MANUAL 1D Seismic Site Response Analysis Example University of California: San Diego August 30, 2017
USER S MANUAL 1D Seismic Site Response Analysis Example http://www.soilquake.net/ucsdsoilmodels/ University of California: San Diego August 30, 2017 Table of Contents USER'S MANUAL TABLE OF CONTENTS Page
More informationStatistical Analysis of Stresses in Rigid Pavement
Statistical Analysis of Stresses in Rigid Pavement Aleš Florian, Lenka Ševelová, and Rudolf Hela Abstract Complex statistical analysis of stresses in concrete slab of the real type of rigid pavement is
More information3D MATERIAL MODEL FOR EPS RESPONSE SIMULATION
3D MATERIAL MODEL FOR EPS RESPONSE SIMULATION A.E. Swart 1, W.T. van Bijsterveld 2, M. Duškov 3 and A. Scarpas 4 ABSTRACT In a country like the Netherlands, construction on weak and quite often wet soils
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 informationDeterministic and stochastic investigation of welded, high strength steel T joint
Szakály Ferenc Student, MSc in Computational Structural Engineering, BUTE (M.Sc.) Deterministic and stochastic investigation of welded, high strength steel T joint Scientific Student s Association Thesis
More informationNumerical Characterization of Concrete Heterogeneity
Vol. Materials 5, No. Research, 3, 2002Vol. 5, No. 3, Statistical 309-314, 2002. Characterization of the Concrete Numerical Modeling of Size Effect In Heterogeneity 2002 309 Numerical Characterization
More informationA Constitutive Framework for the Numerical Analysis of Organic Soils and Directionally Dependent Materials
Dublin, October 2010 A Constitutive Framework for the Numerical Analysis of Organic Soils and Directionally Dependent Materials FracMan Technology Group Dr Mark Cottrell Presentation Outline Some Physical
More informationMicro-meso draping modelling of non-crimp fabrics
Micro-meso draping modelling of non-crimp fabrics Oleksandr Vorobiov 1, Dr. Th. Bischoff 1, Dr. A. Tulke 1 1 FTA Forschungsgesellschaft für Textiltechnik mbh 1 Introduction Non-crimp fabrics (NCFs) are
More informationUsing the Abaqus CDP Model in Impact Simulations
SAFIR project (http://safir2018.vtt.fi/) The Finnish Research Programme on Nuclear Power Plant Safety 2015-2018 Using the Abaqus CDP Model in Impact Simulations Alexis Fedoroff Technical reseach centre
More informationModel for predicting the UHPFRC tensile hardening response
UHPC-2008: The Second International Symposium on Ultra High Performance Concrete, March 05-07, 2008, Kassel, Germany. John Wuest Dr., Civil Engineer Ecole Polytechnique Fédérale de Lausanne Lausanne, Switzerland
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 informationA rate-dependent Hosford-Coulomb model for predicting ductile fracture at high strain rates
EPJ Web of Conferences 94, 01080 (2015) DOI: 10.1051/epjconf/20159401080 c Owned by the authors, published by EDP Sciences, 2015 A rate-dependent Hosford-Coulomb model for predicting ductile fracture at
More informationFinite element analysis of diagonal tension failure in RC beams
Finite element analysis of diagonal tension failure in RC beams T. Hasegawa Institute of Technology, Shimizu Corporation, Tokyo, Japan ABSTRACT: Finite element analysis of diagonal tension failure in a
More informationMODELING 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 informationPlasticity R. Chandramouli Associate Dean-Research SASTRA University, Thanjavur
Plasticity R. Chandramouli Associate Dean-Research SASTRA University, Thanjavur-613 401 Joint Initiative of IITs and IISc Funded by MHRD Page 1 of 9 Table of Contents 1. Plasticity:... 3 1.1 Plastic Deformation,
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 informationIMPACT RESPONSE ANALYSIS OF LARGE SCALE RC GIRDER WITH SAND CUSHION
-Technical Paper- IMPACT RESPONSE ANALYSIS OF LARGE SCALE RC GIRDER WITH SAND CUSHION Abdul Qadir BHATTI *1, Norimitsu KISHI *2, Shin-ya OKADA *3 and Hisashi KONNO *4 ABSTRACT In order to establish a proper
More informationMESOSCOPIC MODELLING OF MASONRY USING GFEM: A COMPARISON OF STRONG AND WEAK DISCONTINUITY MODELS B. Vandoren 1,2, K. De Proft 2
Blucher Mechanical Engineering Proceedings May 2014, vol. 1, num. 1 www.proceedings.blucher.com.br/evento/10wccm MESOSCOPIC MODELLING OF MASONRY USING GFEM: A COMPARISON OF STRONG AND WEAK DISCONTINUITY
More informationExperiments and Numerical Simulations on Stress-State-Dependence of Ductile Damage Criteria
Experiments and Numerical Simulations on Stress-State-Dependence of Ductile Damage Criteria Michael Brünig, Steffen Gerke and Daniel Brenner Abstract The paper deals with a series of new experiments and
More informationANSYS Mechanical Basic Structural Nonlinearities
Lecture 4 Rate Independent Plasticity ANSYS Mechanical Basic Structural Nonlinearities 1 Chapter Overview The following will be covered in this Chapter: A. Background Elasticity/Plasticity B. Yield Criteria
More informationEffect Of The In-Situ Stress Field On Casing Failure *
Effect Of The In-Situ Stress Field On Casing Failure * Tang Bo Southwest Petroleum Institute, People's Republic of China Lian Zhanghua Southwest Petroleum Institute, People's Republic of China Abstract
More informationME 2570 MECHANICS OF MATERIALS
ME 2570 MECHANICS OF MATERIALS Chapter III. Mechanical Properties of Materials 1 Tension and Compression Test The strength of a material depends on its ability to sustain a load without undue deformation
More informationPredicting Fatigue Life with ANSYS Workbench
Predicting Fatigue Life with ANSYS Workbench How To Design Products That Meet Their Intended Design Life Requirements Raymond L. Browell, P. E. Product Manager New Technologies ANSYS, Inc. Al Hancq Development
More informationUSER S MANUAL. 1D Seismic Site Response Analysis Example. University of California: San Diego.
USER S MANUAL 1D Seismic Site Response Analysis Example http://www.soilquake.net/ucsdsoilmodels/ University of California: San Diego August 2, 2017 Table of Contents USER'S MANUAL TABLE OF CONTENTS Page
More informationA Performance Modeling Strategy based on Multifiber Beams to Estimate Crack Openings ESTIMATE in Concrete Structures CRACK
A Performance Modeling Strategy based on Multifiber Beams to Estimate Crack Openings ESTIMATE in Concrete Structures CRACK A. Medjahed, M. Matallah, S. Ghezali, M. Djafour RiSAM, RisK Assessment and Management,
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 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 informationFEM model of pneumatic spring assembly
FEM model of pneumatic spring assembly Tien Tran Xuan 1, David Cirkl 2 Department of Applied Mechanics, Faculty of Mechanical Engineering, Technical University of Liberec, Liberec, Czech Republic 1 Corresponding
More informationA Semianalytical Model for the Simulation of Polymers
A Semianalytical Model for the Simulation of Polymers Paul Du Bois 3, Stefan Kolling 1, Markus Feucht 1 & André Haufe 2 1 DaimlerChrysler AG, Sindelfingen, Germany 2 Dynamore GmbH, Stuttgart, Germany 3
More informationCONTENTS. Lecture 1 Introduction. Lecture 2 Physical Testing. Lecture 3 Constitutive Models
CONTENTS Lecture 1 Introduction Introduction.......................................... L1.2 Classical and Modern Design Approaches................... L1.3 Some Cases for Numerical (Finite Element) Analysis..........
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 informationEXPERIMENTAL IDENTIFICATION OF HYPERELASTIC MATERIAL PARAMETERS FOR CALCULATIONS BY THE FINITE ELEMENT METHOD
Journal of KONES Powertrain and Transport, Vol. 7, No. EXPERIMENTAL IDENTIFICATION OF HYPERELASTIC MATERIAL PARAMETERS FOR CALCULATIONS BY THE FINITE ELEMENT METHOD Robert Czabanowski Wroclaw University
More informationMODELING GEOMATERIALS ACROSS SCALES
MODELING GEOMATERIALS ACROSS SCALES JOSÉ E. ANDRADE DEPARTMENT OF CIVIL AND ENVIRONMENTAL ENGINEERING AFOSR WORKSHOP ON PARTICULATE MECHANICS JANUARY 2008 COLLABORATORS: DR XUXIN TU AND MR KIRK ELLISON
More information