Thermal stress is a type of stress induced by a temperature change, which can lead to fracture or deformation of an object. In the following model a steel beam is fixed at the left and heated by a constant inward heat flux on the top surface. The left and the lower surfaces are kept at while the right end is assumed to be thermally insulated.
During the heating process the temperature on the top surface will be gradually increased. The resulting temperature gradient will then generate thermal stress and bend the beam downwards.
The evolution of the temperature field and the thermally induced deformation are simulated by sequentially the heat transfer and the structural mechanics models. The location and values of horizontal and vertical displacement will then be calculated.
The symbols and corresponding units used here are summarized in the Nomenclature section.
Please refer to the information provided in "Heat Transfer" for a more general theoretical background for heat transfer analysis.
Load the finite element package.
There are two physical domains involved in this application: the heat transfer model and the structural mechanics model. The coupling between two models is one-way since only the structural deformation depends on the temperature but the temperature field does not depend on the deformation of the object in this model. This type of problem is considered as a sequential multiphysics model, and will be solved in two steps.
First the heat transfer model is built to simulate the temperature field of the beam. The structural mechanics model is then constructed and uses the temperature field previously computed to show the thermally induced deformation.
Heat Transfer Model
The heat equation (1) is used to solve for the temperature field in a heat transfer model:
For a transient heat transfer model without sources, the source term in (2) is set to zero. Since a solid is modeled, any internal velocity also vanishes and the heat equation simplifies to:
Set up a 2D transient heat transfer model.
Structural Mechanics Model
In structural mechanics, the plane stress relation (3) describes the deformation of thin objects. Here, "thin" means thin relative to the other dimensions of the object. The two dependent variables , denote the deformation in and directions, respectively. As material data, Young's modulus and the Poisson ratio need to be specified.
The equilibrium equation (4), which denotes the force balance in the and directions, is then used as the governing PDE of the structural mechanics model. On the right hand side and model any external force that may act on the object.
Define a transient plane stress operator.
Set up material parameters for the steel beam.
is the density , is the heat capacity , is the thermal conductivity , is the Young's modulus , is the thermal expansion coefficient , is the Poisson ratio.
Set up the model parameters.
The steel beam to be modeled has a length of and a width of .
Define the 2D domain .
In order to get a good result, a finer than the default grid is applied for the mesh generation. Here the maximum grid size is set to so that about one thousand elements are used.
Discretize the domain with a prescribed mesh size.
Solve the PDE Model
In the following section the heat transfer model will be solved first to simulate the temperature field of the beam. The structural mechanics model is then constructed to show the thermally induced deformation.
Heat Transfer Model
At the beginning of the simulation, the temperature of the steel beam is set to the temperature .
Set up the initial temperature of the beam.
The left and bottom surfaces are kept at the initial temperature at all time.