December 16, 2020 at 12:27 pmPaulvhSubscriber
Currently I'm working on a static structural analysis with two models that describe the failure of individual material and the interface in a multi-material (PC) part. The individual material failure is modelled with the use of EKILL based on the ultimate plastic strain theory. The interface failure is modelled by a separation-distance based debonding (CZM) on contact elements.
The problem that I encounter in this simulation is the high occurence of element HDST violations, i.e. element distortions. To a degree this behavior is expected as when elements are killed via EKILL, their stiffness is reduced significantly, subjecting them to high degrees of deformation. I have however tried many approaches to minimize this effect, including:
- Model setup: decreased loading increments and increased amount of substeps (1mm displacement is loaded linearly ramped over 100 loadsteps, with each having minimum 5 and maximum 1*10^5 (!) substeps), varied CZM artificial damping factor, ensured contact pinball radius is at least the length of the contact gap at the completion of debonding.
- Mesh quality: varied between all mesh methods (multizone/quadrilateral/triangles), varied mesh size (ranging from 0.01mm to 1mm), introduced initial biased aspect ratios of elements that are to be stretched.
Unfortunately, none of the approaches minimizes the occurence of element distortions to an extend where the simulation is able to completely converge in a robust manner. Observing the force convergence solution output and the Newton-Raphson residual force plot (see figures), I can only conclude the simulation is prematurely halted because of the element distortion (highlighted in green in the NR plot).January 7, 2021 at 11:06 amAshish KhemkaAnsys EmployeennWhat is the number of equilibrium iterations you are using? If you have not changed the same then can you try increasing the number of equilibrium iterations using the command in the Setup branch:nneqit,100 ! This increases the equlibrium iterationsnnRegards,nAshish KhemkanJanuary 7, 2021 at 2:02 pmPaulvhSubscriberHi Akhemka,nI have come across this suggestion on another forum as well and attempted with increased neqit,50. Unfortunately I noticed no changes. Thanks for thinking along though!nIn the meantime a promising change I have made was changing the discplacement boundary condition to a remote displacement boundary condition, allowing me to change its behaviour from rigid to deformable. Switching between these settings drastically affected the robustness of the model with limited effect on the results. nKind regards,nPaulnViewing 2 reply threads
Ansys Innovation Space
- You must be logged in to reply to this topic.
Simulation World 2022
Check out more than 70 different sessions now available on demand. Get inspired as you hear from visionary companies, leading researchers and educators from around the globe on a variety of topics from life-saving improvements in healthcare, to bold new realities of space travel. Take a leap of certainty and check out a session today here.
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.
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.Trending discussions
- How to calculate the residual stress on a coating by Vickers indentation?
- The solver engine was unable to converge on a solution for the nonlinear problem as constrained.
- Command Object Snippet
- whether have the difference between using contact and target bodies
- Errors – Reinforced Concrete Beam
- Massive amount of memory (RAM) required for solve
- An Unknown error occurred during solution. Check the Solver Output…..
- Large deflection
- Transient Structural analysis?
- Revolute joint (Static structural) Very high deformation
Top Rated Tags