Coupled Structural/Thermal Analysis: University of Alberta
Coupled Structural/Thermal Analysis: University of Alberta
Coupled Structural/Thermal Analysis: University of Alberta
Design Optimization
Coupled Structural/Thermal Analysis
Substructuring
Index Loads will not be applied to the link, only a temperature change of 75 degrees Celsius. The link is steel with a
modulus of elasticity of 200 GPa, a thermal conductivity of 60.5 W/m*K and a thermal expansion coefficient of
12e-6 /K.
Contributions
Comments
MecE 563
Mechanical Engineering
University of Alberta
ANSYS Inc.
Copyright © 2001
University of Alberta
Thus, each different physics environment must be constructed seperately so they can be used to determine the
coupled physics solution. However, it is important to note that a single set of nodes will exist for the entire
model. By creating the geometry in the first physical environment, and using it with any following coupled
environments, the geometry is kept constant. For our case, we will create the geometry in the Thermal
Environment, where the thermal effects will be applied.
Although the geometry must remain constant, the element types can change. For instance, thermal elements are
required for a thermal analysis while structural elements are required to deterime the stress in the link. It is
important to note, however that only certain combinations of elements can be used for a coupled physics
analysis. For a listing, see Chapter 2 of the ANSYS Coupled-Field Guide located in the help file.
The process requires the user to create all the necessary environments, which are basically the preprocessing
portions for each environment, and write them to memory. Then in the solution phase they can be combined to
solve the coupled analysis.
www.mece.ualberta.ca/…/Coupled.html 1/7
5/10/2010 U of A ANSYS Tutorials - Coupled Stru…
Thermal Environment - Create Geometry and Define Thermal Properties
3. Define Keypoints
Preprocessor > Modeling > Create > Keypoints > In Active CS...
K,#,x,y,z
We are going to define 2 keypoints for this link as given in the following table:
4. Create Lines
Preprocessor > Modeling > Create > Lines > Lines > In Active Coord
L,1,2
For this problem we will use the LINK33 (Thermal Mass Link 3D conduction) element. This
element is a uniaxial element with the ability to conduct heat between its nodes.
In the 'Real Constants for LINK33' window, enter the following geometric properties:
Preprocessor > Material Props > Material Models > Thermal > Conductivity > Isotropic
In the window that appears, enter the following geometric properties for steel:
i. KXX: 60.5
Preprocessor > Meshing > Size Cntrls > ManualSize > Lines > All Lines...
For this example we will use an element edge length of 0.1 meters.
Preprocessor > Meshing > Mesh > Lines > click 'Pick All'
The thermal environment (the geometry and thermal properties) is now fully described and can be
written to memory to be used at a later time.
Preprocessor > Physics > Environment > Write
In the window that appears, enter the TITLE Thermal and click OK.
Doing this clears all the information prescribed for the geometry, such as the element type, material
properties, etc. It does not clear the geometry however, so it can be used in the next stage, which
is defining the structural environment.
Since the geometry of the problem has already been defined in the previous steps, all that is required is to detail
the structural variables.
This will switch to the complimentary structural element automatically. In this case it is LINK 8.
For more information on this element, see the help file. A warning saying you should modify the
new element as necessary will pop up. In this case, only the material properties need to be
modified as the geometry is staying the same.
Preprocessor > Material Props > Material Models > Structural > Linear > Elastic > Isotropic
In the window that appears, enter the following geometric properties for steel:
Preprocessor > Material Props > Material Models > Structural > Thermal Expansion Coef >
Isotropic
i. ALPX: 12e-6
3 W it E i t
www.mece.ualberta.ca/…/Coupled.html 3/7
5/10/2010 U of A ANSYS Tutorials - Coupled Stru…
3. Write Environment
If the Physics option is not available under Solution, click Unabridged Menu at the bottom of the
Solution menu. This should make it visible.
3. Apply Constraints
Solution > Define Loads > Apply > Thermal > Temperature > On Keypoints
It is very important to click Finish as it closes that environment and allows a new one to be
opened without contamination. If this is not done, you will get error messages.
The thermal solution has now been obtained. If you plot the steady-state temperature on the link, you will
see it is a uniform 348 K, as expected. This information is saved in a file labelled Jobname.rth, were
.rth is the thermal results file. Since the jobname wasn't changed at the beginning of the analysis, this data
can be found as file rth We will use these results in determing the structural effects
www.mece.ualberta.ca/…/Coupled.html 4/7
5/10/2010 U of A ANSYS Tutorials - Coupled Stru…
can be found as file.rth. We will use these results in determing the structural effects.
7. Apply Constraints
Solution > Define Loads > Apply > Structural > Displacement > On Keypoints
Solution > Define Loads > Apply > Structural > Temperature > From Therm Analy
As shown below, enter the file name File.rth. This couples the results from the solution of the
thermal environment to the information prescribed in the structural environment and uses it during
the analysis.
Preprocessor > Loads > Define Loads > Settings > Reference Temp
For this example set the reference temperature to 273 degrees Kelvin.
Hand calculations were performed to verify the solution found using ANSYS:
www.mece.ualberta.ca/…/Coupled.html 5/7
5/10/2010 U of A ANSYS Tutorials - Coupled Stru…
As shown, the stress in the link should be a uniform 180 MPa in compression.
Since the element is only a line, the stress can't be listed in the normal way. Instead, an element
table must be created first.
General Postproc > Element Table > Define Table > Add
Fill in the window as shown below. [CompStr > By Sequence Num > LS > LS,1
ETABLE,CompStress,LS,1
General Postproc > Element Table > List Elem Table > COMPSTR > OK
PRETAB,CompStr
The following list should appear. Note the stress in each element: -0.180e9 Pa, or 180 MPa in
compression as expected.
www.mece.ualberta.ca/…/Coupled.html 6/7
5/10/2010 U of A ANSYS Tutorials - Coupled Stru…
www.mece.ualberta.ca/…/Coupled.html 7/7