March 14, 2021 at 6:31 pmAnsuman_SahuSubscriber
I am performing a transient thermal analysis of fused deposition modelling process using ANSYS APDL (Ansys 2020 R2 - Student version). I am trying to simulate the deposition of a single extruded strand using element birth and death technique. In each load step one element is activated. The initial temperature of the activated element (in the current load step) is to be the nozzle temperature. The initial temperatures of the previously activated elements is obtained from the results (temperature field) of the previous load step.
- The steps adopted in the model are as follows:
- Creating & meshing the 3D model (of a single extruded strand of length=50mm)
- Deactivating all elements (using EKILL)
- Setting up the solution
- Activating one element
- Specifying initial condition on nodes
- Specifying boundary conditions (Convection load on exterior node)
- Specifying the necessary load step options and solving.
- Repeating steps 4 to 7 for each activated element (using *DO loop).
I am facing the following issues:
- The analysis runs properly only for the first load step, in that the decrease in temperature of the element is in accordance with the convective loading applied on it. But from the second and subsequent load steps onwards, the temperature distribution obtained is not proper. When the second element is activated, the face of the element exposed to ambient air has a constant temperature = nozzle temperature even though convective load is applied on this element. Element numbers 3 and beyond, remain entirely at the constant temperature = nozzle temperature. Note: If however, each load step is run in a separate file (with all other parameters constant and by reading in the solution of the previous load step as the initial condition), the results of each load step are in accordance with the boundary conditions and loads assigned to them. I am unable to identify what's causing these improper nodal solutions when solving all the load steps in a single file. Am I missing out on any crucial step in my code? Why is there a difference in nodal solution between these two approaches (solving all load steps in one file v/s solving each load step in a separate file?
- When assigning initial conditions, the temperature distribution from the previous load steps can be read (using LDREAD) for the previously activated elements. However, I am unable to separately provide the initial temperature to the element activated in the current load step. The following warning message gets generated : Initial Conditions not allowed after first Load Step. Currently I have managed to avoid this issue by assigning a uniform temperature (equal to the nozzle temperature) at the beginning of solution for all the nodes (of activated as well as deactivated elements), but I am not sure of it. Is there a proper way to separately provide the initial condition for element activated in the current step?
- Also, is LDREAD command the proper method for using the results of the previous load step as the initial condition for the current step?
Since I am new to ANSYS, any suggestions w.r.t these two issues and help in identifying the mistakes in my code/approach will be greatly appreciated.
The relevant screenshots from APDL as well as the code have been attached for reference.March 19, 2021 at 1:25 pmDave LoomanAnsys EmployeeYour description of the issue is correct. It is difficult to make the initial temperature of an element what you want when it is activated. When elements are inactive (dead) they are still there in the finite element model, they just have watered down material properties. As a result, they can change temperature through conduction from the adjacent active elements. One approach to giving them a desired initial temperature is to constrain their nodes to that temperature while they are inactive and then delete those constraints when they are activated. This method has the drawback that the inactive elements share nodes with the adjacent active elements and those nodes can't be constrained. A thin layer of inactive elements at the interface might help this. Another method is to apply heat to the elements when they are activated to bring them up to the desired temperature. nApril 1, 2021 at 2:43 pmAnsuman_SahuSubscriberThank you very much for these suggestions. I am sorry for replying late. nI tried the first suggestion of constraining the nodes of the inactive elements at the initial temperature, and got a reasonable spatial temperature distribution and thermal histories. However, the drawback, which you have mentioned was quite clear from the contour plots. Although I have managed to resolve it to some extent by using a finer mesh, but the computation time is certainly quite high now and also because I am using a student version, there's the limitation on the number of nodes/elements. nI will definitely work on the other two approaches which you have suggested. nViewing 2 reply threads
Ansys Innovation Space
- 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.
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
- Saving & sharing of Working project files in .wbpz format
- Understanding Force Convergence Solution Output
- An Unknown error occurred during solution. Check the Solver Output…..
- Solver Pivot Warning in Beam Element Model
- Colors and Mesh Display
- whether have the difference between using contact and target bodies
- How to calculate the residual stress on a coating by Vickers indentation?
- What is the difference between bonded contact region and fixed joint
- The solver engine was unable to converge on a solution for the nonlinear problem as constrained.
- User manual
Top Rated Tags
© 2023 Copyright ANSYS, Inc. All rights reserved.Ansys does not support the usage of unauthorized Ansys software. Please visit www.ansys.com to obtain an official distribution.