Right now I'm running a coupled temp-disp analysis in ABAQUS/CAE, however, my job is very slow, 15 minutes have passed and only 1 increment has been made, I tried modifying the step to a time period of 1, maximum number of increments of 50, initial increment size 0.01, maximum 1 and minimum 1e-15. When I increase the maximum number of increments to a 1000 that seems to run faster but it may take days to finish. What can I do?
Check the warning messages printed when the solver is running. The model seems to be underconstrained. There's just one BC in the direction normal to the plate.
Check "Spurious oscillations due to small time increments" in help manual. "In transient heat transfer analysis with second-order elements there is a relationship between the minimum usable time increment and the element size."
Some observations
-1W/mm^3 heat flux defined.
-Check units of material data. Unit of density for INVAR 36 seems wrong.
The warning messages are:
For contact pair (assembly__pickedsurf14-assembly__pickedsurf13), not all the nodes that have been adjusted were printed. Specify *preprint,contact=yes for complete printout.
Surface interaction intprop-1, associated with one or more contact pairs with surfaces having underlying elements with temperature and displacement degrees of freedom, does not have any thermal interaction defined. No thermal interaction will occur.
Where is the heat flux defined? Should I change it?
Concerning the units of the INVAR 36, I got the density from this link https://f.hubspotusercontent20.net/hubfs/7407324/M...
Since my model is in mm I converted the .291 lb/in3 to tonne/mm3, I got 8.05e-6 tonne/mm3
Apart from the initial temperature, there is a body heat flux load with a unit magnitude in your file. Maybe that's what causes problems, remove it if you haven't defined it intentionally.
Also, enable parallelization in job settings and specify the number of CPUs that your machine has. It will run much faster then.
After enabling parallelization and changing the density the job ran much faster but the job aborted at the 50th increment, which was my maximum, should I increase the maximum number of increments? The only warnings I got were the warnings that I already talked about.
If there is no any other error that is related to convergence, you need to definitely increase the maximum increments as analysis requires more increments to complete. You can use 1000 increments as maximum increments if you not sure about how much is required.
Please refrain from double posting. You can use this post to ask more questions related to this (same) post in the same post.
For the heat flux I just realized that I did add one. I'm trying to mimick an autoclave process. For this I needed the model to pass through different temperatures at different times, to achieve this I added an amplitude, with the time being in seconds and the amplitude representing the temperature in celsius. Is this correct?
Yes, you can do it this way. Just remember to specify a magnitude of 1 for the BC or load using this amplitude and to stick with Celsius units for the rest of the analysis setup.
My job aborted after 207 increments, the error messages were:
Time increment required is less than the minimum specified
The analysis has been terminated due to previous errors. All output requests have been written for the last converged increment.
The warnings were:
For contact pair (assembly__pickedsurf14-assembly__pickedsurf13), not all the nodes that have been adjusted were printed. Specify *preprint,contact=yes for complete printout.
Surface interaction intprop-1, associated with one or more contact pairs with surfaces having underlying elements with temperature and displacement degrees of freedom, does not have any thermal interaction defined. No thermal interaction will occur.
In the message file, the following warnings were present:
***WARNING: THE SYSTEM MATRIX HAS 3 NEGATIVE EIGENVALUES.
EXPLANATIONS ARE SUGGESTED AFTER THE FIRST OCCURRENCE OF THIS MESSAGE.
***WARNING: SOLVER PROBLEM. NUMERICAL SINGULARITY WHEN PROCESSING NODE
COMPOSITE-1.27 D.O.F. 1 RATIO = 32.7548E+09 .
***WARNING: SOLVER PROBLEM. NUMERICAL SINGULARITY WHEN PROCESSING NODE
COMPOSITE-1.27 D.O.F. 2 RATIO = 33.0342E+09 .
My minimum time increment is 1e-15 so I don't think the problem is there, it is already pretty low. What can I do?
No, how do I eliminate that? Right now the only BC that I have is a Displacement/Rotation with U3=0 in the bottom face of the larger body. In step time it went till 0.000155 on the last increment.
Try with also U1 and U2 fixed for this face. If it doesn’t help, replace contact with tie constraint. If it still doesn’t work, we will know that the problem is not caused by an underconstraint. Then the thermal aspects will be the likely cause.
With the U1 and U2 fixed, the singularities disappear but I notice that the top surface moves away from the bottom, I don't know if it's a problem with the contact. I tried ticking the "tie ajusted surfaces" box but that creates excessive distortion warnings.
Try using tie constraint (available from the Create constraint dialog in the Interaction module). You just may have to disable adjustment and refine the mesh a bit.
I added contact control and that seems to keep the top surface in place. However, the results were excessive with the von mises results being up to e+4, I showed my boss at my internship and he said that I should change the body heat flux to a surface heat flux for more accuracy and he also made me change the boundary conditions to:
1. Displacement/Rotation BC at the bottom layer with U3=0
2. Displacement/Rotation BC at a point with U2=0
3. Displacement/Rotation BC at a point with U1=0
I ran the job and currently I'm receiving extreme distortion warnings and also the bottom surface is moving. I will upload the updated model.
Learn methods and guidelines for using stereolithography (SLA) 3D printed molds in the injection molding process to lower costs and lead time. Discover how this hybrid manufacturing process enables on-demand mold fabrication to quickly produce small batches of thermoplastic parts. Download Now
Examine how the principles of DfAM upend many of the long-standing rules around manufacturability - allowing engineers and designers to place a part’s function at the center of their design considerations. Download Now
Metal 3D printing has rapidly emerged as a key technology in modern design and manufacturing, so it’s critical educational institutions include it in their curricula to avoid leaving students at a disadvantage as they enter the workforce. Download Now
This ebook covers tips for creating and managing workflows, security best practices and protection of intellectual property, Cloud vs. on-premise software solutions, CAD file management, compliance, and more. Download Now
RE: ABAQUS/CAE JOB IS RUNNING SLOW
RE: ABAQUS/CAE JOB IS RUNNING SLOW
Some observations
-1W/mm^3 heat flux defined.
-Check units of material data. Unit of density for INVAR 36 seems wrong.
RE: ABAQUS/CAE JOB IS RUNNING SLOW
For contact pair (assembly__pickedsurf14-assembly__pickedsurf13), not all the nodes that have been adjusted were printed. Specify *preprint,contact=yes for complete printout.
Surface interaction intprop-1, associated with one or more contact pairs with surfaces having underlying elements with temperature and displacement degrees of freedom, does not have any thermal interaction defined. No thermal interaction will occur.
Where is the heat flux defined? Should I change it?
Concerning the units of the INVAR 36, I got the density from this link https://f.hubspotusercontent20.net/hubfs/7407324/M...
Since my model is in mm I converted the .291 lb/in3 to tonne/mm3, I got 8.05e-6 tonne/mm3
RE: ABAQUS/CAE JOB IS RUNNING SLOW
Do you use parallelization (multiple CPUs) for this analysis ?
Heat flux is your thermal load. Just make sure that the value is correct. Units in thermal analyses can be really tricky.
It should be 8.05e-9 after conversion.
RE: ABAQUS/CAE JOB IS RUNNING SLOW
I didn't use parallelization.
I don't remember adding heat flux, I added an initial temperature of 20 Celsius.
Thanks for telling me, I will change the density.
RE: ABAQUS/CAE JOB IS RUNNING SLOW
Also, enable parallelization in job settings and specify the number of CPUs that your machine has. It will run much faster then.
RE: ABAQUS/CAE JOB IS RUNNING SLOW
RE: ABAQUS/CAE JOB IS RUNNING SLOW
Please refrain from double posting. You can use this post to ask more questions related to this (same) post in the same post.
RE: ABAQUS/CAE JOB IS RUNNING SLOW
Here is the image of the amplitude.
RE: ABAQUS/CAE JOB IS RUNNING SLOW
RE: ABAQUS/CAE JOB IS RUNNING SLOW
Time increment required is less than the minimum specified
The analysis has been terminated due to previous errors. All output requests have been written for the last converged increment.
The warnings were:
For contact pair (assembly__pickedsurf14-assembly__pickedsurf13), not all the nodes that have been adjusted were printed. Specify *preprint,contact=yes for complete printout.
Surface interaction intprop-1, associated with one or more contact pairs with surfaces having underlying elements with temperature and displacement degrees of freedom, does not have any thermal interaction defined. No thermal interaction will occur.
In the message file, the following warnings were present:
***WARNING: THE SYSTEM MATRIX HAS 3 NEGATIVE EIGENVALUES.
EXPLANATIONS ARE SUGGESTED AFTER THE FIRST OCCURRENCE OF THIS MESSAGE.
***WARNING: SOLVER PROBLEM. NUMERICAL SINGULARITY WHEN PROCESSING NODE
COMPOSITE-1.27 D.O.F. 1 RATIO = 32.7548E+09 .
***WARNING: SOLVER PROBLEM. NUMERICAL SINGULARITY WHEN PROCESSING NODE
COMPOSITE-1.27 D.O.F. 2 RATIO = 33.0342E+09 .
My minimum time increment is 1e-15 so I don't think the problem is there, it is already pretty low. What can I do?
RE: ABAQUS/CAE JOB IS RUNNING SLOW
Those are the singularities I’ve mentioned before. Have you changed the boundary condition to eliminate rigid body motions ?
RE: ABAQUS/CAE JOB IS RUNNING SLOW
RE: ABAQUS/CAE JOB IS RUNNING SLOW
RE: ABAQUS/CAE JOB IS RUNNING SLOW
RE: ABAQUS/CAE JOB IS RUNNING SLOW
RE: ABAQUS/CAE JOB IS RUNNING SLOW
1. Displacement/Rotation BC at the bottom layer with U3=0
2. Displacement/Rotation BC at a point with U2=0
3. Displacement/Rotation BC at a point with U1=0
I ran the job and currently I'm receiving extreme distortion warnings and also the bottom surface is moving. I will upload the updated model.