TAGGED: #mechanical#workbench, microplanemodel, softening


February 10, 2023 at 3:43 pmBabaloluwa PopoolaSubscriber
Hello, I was able to define nonlinear material properties for concrete using APDL commands provided by Professor Brock D. Hedegaard as a base in Ansys Mechanical. The parameters are for a "DruckerPrager Microplane plasticity " model and attached below is the code:
/PREP7! Step 1: Define elastic properties of materialE = 30000 ! Young's Modulus, MPanu = 0.2 ! Poisson's Ratio, unitless! Step 2: Define microplane model properties!  Approximations and typical ranges per intuition... your concrete may be different!fc = 80 ! Uniaxial compressive strength, MPafb = 92 ! Biaxial compressive strength, MPa (estimate as 1.15*fc if unknown)ft = 2 ! Tensile strength, MPa (estimate as 0.1*fc if unknown)Rt = 1 ! Tension cap hardening constant, unitless (typically = 1)D = 20e4 ! Compressive hardening constant, MPa^2 (typically between 1e4 and 50e4)sigma_cv = 32 ! Int. of comp. cap and DP yield, MPa (always negative, approximately 2/3*fb or lower)Rc = 2 ! Compression cap ratio constant, unitless (typically = 2)gamma_t = 0 ! Tensile damage threshhold, unitless (typically = 0)gamma_c = 5e5 ! Compressive damage threshhold, unitless (typically between 1e5 to 10e5)beta_t = 6000 ! Tension damage evolution, unitless (typically 1.5*beta_c)beta_c = 4000 ! Compression damage evolution, unitless (typically between 1000 and 10000)! Step 3: Define nonlocal parameters c and mc = 650 ! Nonlocal range parameter, mm^2 (element size should be < 0.5*sqrt(c))m = 2 ! Overnonlocal parameter, unitless (typically 1 to 3)! Step 4: Select element typeelement = 215 ! Element type number (Linear = 215, Quadratic = 216)! ===================================================================! You should not need to change code below this point, unless! you are debugging, modifying, or otherwise know what you are doing! ===================================================================! Assign Elastic PropertiesMP,EX,conc,E ! Assign Young's ModulusMP,NUXY,conc,nu ! Assign Poisson's Ratio! Assign microplane properties for coupled damagedplasticity (MPLANEDPC)TB,MPLANE,conc,,,DPCTBDATA,1,fc,fb,ft,Rt,D,sigma_cvTBDATA,7,Rc,gamma_t,gamma_c,beta_t,beta_c! Assign nonlocal properties to microplane model (MPLANENLOCAL)TB,MPLANE,conc,,,NLOCALTBDATA,1,c,m! Define coupled elements (this code assumes ITYPE = matid, which it should be in almost all cases)itype = concet,itype,elementkeyopt,itype,18,2 ! Element needs 2 nonlocal parameters for keyopt(18)! Switch element type of those that are assigned Concreteesel,s,mat,,concemodif,all,type,itype! To check that everything worked out, print the element listallseletlist,all! Now proceed to the solution and output all results/SOLUoutres,all,allTrying to validate the defined parameters, I modelled a generic beam and subjected it to uniaxial and flexural deflection schemes.The uniaxial loading conditions gave the expected results, with both compressive and tensile stress values peaking at the defined strength parameters before softening occurred. I have attached images below:Fig.1 Peak compression
Fig.2 Peak tension
In contrast, flexural loading conditions give results for peak normal stress which significantly exceed defined parameters, even after varying the applied displacement on the beam.
Fig .3 Peak compressive and tensile stresses for the flexural loading condition.
It is my understanding that the normal stress in the zdirection is composed of Tensile and compressive stresses beneath and above the neutral axis respectively. I would appreciate feedback on possible errors in my assumptions or potential reasons for the disparities in results.
Thanks in Advance.

February 14, 2023 at 2:19 pmAshish KhemkaAnsys Employee
Concrete Material Model (ansys.com)
Query regarding Validation of Technology Showcase Example 54 by ANSYS?
Regards,
Ashish Khemka


February 16, 2023 at 7:55 pmBabaloluwa PopoolaSubscriber
Thanks for the provided links and response Ashish
The apdl commands used in the youtube video served as the base for the one I used. The displacement applied to the beam in the video induced a peak stress which was below the defined yield stress, whereas in my analysis I'm trying to investigate the behaviour of the material when damage/softening occurs.
I also went through the other links and found the information very useful, however, none of the changes I made seem to address the large variations in peak stress when comparing both systems, even after I varied damage parameters.
Thanks for the provided links and response Ashish
The apdl commands used in the youtube video served as the base for the one I used. The displacement applied to the beam in the video induced a peak stress which was below the defined yield stress, whereas in my analysis I'm trying to investigate the behaviour of the material when damage/softening occurs.
I also went through the other links and found the information very useful, however, none of the changes I made seem to address the large variations in peak stress when comparing both systems, even after I varied damage parameters.
 You must be logged in to reply to this topic.
Boost Ansys Fluent Simulations with AWS
Computational Fluid Dynamics (CFD) helps engineers design products in which the flow of fluid components is a significant challenge. These different use cases often require large complex models to solve on a traditional workstation. Click here to join this event to learn how to leverage Ansys Fluids on the cloud, thanks to Ansys Gateway powered by AWS.
Earth Rescue – An Ansys Online Series
The climate crisis is here. But so is the human ingenuity to fight it. Earth Rescue reveals what visionary companies are doing today to engineer radical new ideas in the fight against climate change. Click here to watch the first episode.
Ansys Blog
Subscribe to the Ansys Blog to get great new content about the power of simulation delivered right to your email on a weekly basis. With content from Ansys experts, partners and customers you will learn about product development advances, thought leadership and trends and tips to better use Ansys tools. Sign up here.
 Saving & sharing of Working project files in .wbpz format
 Solver Pivot Warning in Beam Element Model
 Understanding Force Convergence Solution Output
 An Unknown error occurred during solution. Check the Solver Output…..
 What is the difference between bonded contact region and fixed joint
 whether have the difference between using contact and target bodies
 The solver engine was unable to converge on a solution for the nonlinear problem as constrained.
 User manual
 Colors and Mesh Display
 material damping and modal analysis

3744

2573

1803

1236

594
© 2023 Copyright ANSYS, Inc. All rights reserved.