U.S. Department of Transportation
Federal Highway Administration
1200 New Jersey Avenue, SE
Washington, DC 20590
202-366-4000
Federal Highway Administration Research and Technology
Coordinating, Developing, and Delivering Highway Transportation Innovations
This report is an archived publication and may contain dated technical, contact, and link information |
|
Publication Number: FHWA-HRT-05-062
Date: May 2007 |
|||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Users Manual for LS-DYNA Concrete Material Model 159PDF Version (1.49 KB)
PDF files can be viewed with the Acrobat® Reader® ForewordThis report documents a concrete material model that has been implemented into the dynamic finite element code, LS-DYNA, beginning with version 971. This model is in keyword format as MAT_CSCM for Continuous Surface Cap Model. This material model was developed to predict the dynamic performance-both elastic deformation and failure-of concrete used in roadside safety structures when involved in a collision with a motor vehicle. An example of a roadside safety structure is a concrete safety barrier that divides opposing lanes of traffic on a roadway. Default input parameters for concrete are stored in the model and can be accessed for use. This material model only replicates the concrete aggregate. Appropriate reinforcement bars or rods must be included in the structure model separately. The Users Manual for LS-DYNA Concrete Material Model 159 is the first of two reports that completely document this material model. This report documents the theoretical basis, the required input format, and includes limited hypothetical problems for the user. The second report, Evaluation of LS-DYNA Concrete Material Model 159 (FHWA-HRT-05-063), documents the testing performed to document the model's performance and accuracy of results. This report will be of interest to research engineers who are associated with the evaluation and crashworthy performance of roadside safety structures, particularly engineers responsible for predicting the crash response of such structures when using the finite element code, LS-DYNA. Michael Trentacoste Director, Office of Safety R&D Notice This document is disseminated under the sponsorship of the U.S. Department of Transportation in the interest of information exchange. The U.S. Government assumes no liability for the use of the information contained in this document. The U.S. Government does not endorse products or manufacturers. Trademarks or manufacturers' names appear in this report only because they are considered essential to the objective of the document. Quality Assurance Statement The Federal Highway Administration (FHWA) provides high-quality information to serve Government, industry, and the public in a manner that promotes public understanding. Standards and policies are used to ensure and maximize the quality, objectivity, utility, and integrity of its information. FHWA periodically reviews quality issues and adjusts its programs and processes to ensure continuous quality improvement. Technical Report Documentation Page
Form DOT F 1700.7 (8-72) Reproduction of completed page authorized SI (Modern Metric) Conversion FactorsTable of ContentsChapter 2. Theoretical Manual Chapter 3. Users Manual Appendix A. Modeling Softening Appendix C. Single Element Input File Appendix D. CEB Specification for Rate Effects List of FiguresFigure 1. Graph. Example concrete data from Mills and Zimmermann plotted in the meridian plane.(8) Figure 3. Graph. Example plots of the failure surfaces of LS-DYNA Model 159 in the meridian plane Figure 6. Graph. Variation of concrete softening response with confinement. Source: Joy and Moxley.(12) Figure 15. Illustration. General shape of the concrete model yield surface in three dimensions Figure 17. Equation. Stress invariant J 1, J′2, and J′3 Figure 18. Equation. Yield function f. Figure 19. Equation. Shear failure surface function Ff. Figure 20. Graph. Schematic of shear surface Figure 21. Graph. Schematic of two-part cap function Figure 22. Graph. Schematic of multiplicative formulation of the shear and cap surfaces Figure 23. Equation. Cap failure surface function Fc Figure 24. Equation. L of kappa Figure 25. Equation. Simple cap failure surface function Fc Figure 26. Equation. X as a function of kappa. Figure 27. Equation. Plastic volume strain ε pv Figure 29. Equation. Angle beta hat in the deviatoric plane Figure 30. Equation. Relationship between beta hat and J hat Figure 31. Equation. Rubin scaling function ℜ Figure 32. Equation. Most general form for Q1 and Q2 Figure 33. Equation. Mohr-Coulomb form for Q1, Q2 Figure 34. Equation. Willam-Warnke form for Q1 Figure 35. Equation. Damaged stress σ dij Figure 36. Graph. This cap model simulation demonstrates strain softening and modulus reduction Figure 37. Equation. Brittle damage threshold τb Figure 38. Equation. Ductile damage threshold τd Figure 39. Equation. Viscoplastic damage threshold r0 Figure 40. Equation. Incremental damage threshold, small rn+1 Figure 41. Equation. Brittle damage small d of tau Figure 42. Equation. Ductile damage small d of tau Figure 43. Equation. Variation of dmax with stress invariant ratio Figure 44. Equation. Variation of dmax with rate effects Figure 45. Schematic representation of four stress paths and their stress invariant ratios Figure 46. Equation. Reduction of A with confinement Figure 47. Equation. Fracture energy integral for Gf Figure 48. Equation. Brittle damage fracture energy Gf Figure 49. Equation. Brittle damage threshold difference τ minus small r0b Figure 50. Equation. Brittle softening parameter C Figure 51. Equation. Ductile damage fracture energy Gf Figure 52. Equation. Ductile damage threshold difference τ − r0d Figure 53. Equation. Ductile softening parameter A Figure 54. Equation. Brittle damage threshold Gf Brittle Figure 55. Equation. Ductile damage threshold Gf Ductile Figure 56. Equation. The fracture energy with rate effects, Gvpf Figure 57. Equation. Default damage recovery of d of τt Figure 58. Equation. Optional damage recovery of d of τt Figure 59. Equation. Viscoplastic stress update for σvpij Figure 60. Equation. Two-parameter η Figure 61. Equation. Dynamic strengths, f ′T dynamic, and f ′C dynamic Figure 62. Equation. Effective strain rate Figure 63. Equation. Variation of fluidity parameter η in tension Figure 64. Equation. Variation of fluidity parameter η in compression Figure 65. Equation. Effective fluidity parameters, ηt, ηc, and η s Figure 66. Equation. Overstress limit of η Figure 67. Equation. Back stress α ij n + 1 Figure 68. Equation. Updated stress with hardening, σPij n+1 Figure 69. Equation. Incremental back stress, Δαij Figure 70. Equation. Brittle rate of translation CH Brittle Figure 71. Equation. Ductile rate of translation CH Ductile Figure 72. Equation. The limiting function Gα. Figure 73. Equation. Modified shear failure surface, Ff Figure 74. Equation. Default Young's modulus E Figure 75. Equation. Shear and bulk moduli, G and K Figure 76. Equation. ACI Young's modulus, Ec Figure 77. Equation. Reduced ACI Young's modulus, Ec Figure 78. Equation. TXC Strength Figure 79. Equation. Interpolation parameter P Figure 80. Equation. Most general form for Q1, Q2 Figure 82. Equation. The default fracture energy GF Figure 84. Illustration. General shape of the concrete model yield surface in two dimensions Figure 85. Equation. Three stress invariants, J1, J′2, J′3 Figure 86. Equation. Plasticity yield function f Figure 87. Equation. Shear surface function Ff Figure 88. Equation. Most general form for scaling functions Q1, Q2 Figure 89. Equation. Cap surface function, Fc Figure 90. Equation. Definition of L of kappa Figure 91. Equation. Pressure invariant X as a function of kappa Figure 92. Equation. Plastic volume strain hardening rule, ε pv Figure 93. Equation. Transformation of viscoplastic stress to damaged stress, σdij. Figure 94. Equation. Ductile damage accumulation, τ d Figure 95. Equation. Brittle damage accumulation, τb Figure 96. Equation. Brittle damage, d of tb Figure 97. Equation. Ductile damage, d of τd Figure 98. Equation. Reduction of A with confinement Figure 99. Equation. Brittle and ductile damage thresholds, Gf Figure 100. Equation. Viscoplastic stress, σvpij Figure 101. Equation. Variation of the fluidity parameter η in tension and compression Figure 102. Definition of effective strain rate Figure 103. Equation. Overstress limit of η Figure 104. Equation. Fracture energy with rate effects Figure 108. Equation. Old generic damage, small d of τ Figure 109. Equation. New generic damage, small d of τ Figure 110. Graph. Behavior of the original softening function Figure 111. Graph. Behavior of the updated softening function Figure 114. Equation. CEB tensile strength dynamic increase factor, DIFten Figure 115. Equation. CEB compressive strength dynamic increase factor, DIFcomp Figure 116. Graph. Dynamic increase factors specified in CEB List of TablesTable 3. Approximate strength measurements used to set default TXC yield surface parameters. Table 4. TXC yield surface input parameters as a function of unconfined compression strength. Table 6. TOR yield surface input parameters as a function of unconfined compression strength. Table 7. TXE yield surface input parameters as a function of unconfined compression strength. Table 9. Coefficients for the fracture energy equation. Table 10. Tensile fracture energies tabulated in CEB as a function of concrete strength. Table 11. Example load curve for modeling rebar strain rate effects with LS-DYNA Material Model #24.
|