6.6.2 Coupled thermal-electrical analysis

**Products: **ABAQUS/Standard ABAQUS/CAE

Coupled thermal-electrical problems:

are those in which coupling between the electrical potential and temperature fields make it necessary to solve both fields simultaneously;

require the use of coupled thermal-electrical elements, although pure heat transfer elements can also be used in the model;

can include a specification of the fraction of electrical energy that will be released as heat;

can include thermal interactions such as gap radiation, gap conductance, and heat generation between surfaces (see “Thermal contact properties,” Section 30.2.1);

can include cavity radiation effects (see “Cavity radiation,” Section 32.1.1);

can include electrical interactions such as electrical current flowing across surfaces (see “Electrical contact properties,” Section 30.3.1);

allow for transient or steady-state thermal solutions and for steady-state electrical solutions; and

can be linear or nonlinear.

Coupled thermal-electrical analysis

Joule heating arises when the energy dissipated by an electrical current flowing through a conductor is converted into thermal energy. ABAQUS/Standard provides a fully coupled thermal-electrical procedure for analyzing this type of problem: the coupled thermal-electrical equations are solved simultaneously for both temperature and electrical potential at the nodes.

The capability includes the analysis of the electrical problem, the thermal problem, and the coupling between the two problems. Coupling arises from two sources: temperature-dependent electrical conductivity and internal heat generation, which is a function of the electrical current density. The thermal part of the problem can include heat conduction and heat storage (“Thermal properties: overview,” Section 20.2.1) as well as cavity radiation effects (“Cavity radiation,” Section 32.1.1). Forced convection caused by fluid flowing through the mesh is not considered.

The thermal-electrical equations are unsymmetric; therefore, the unsymmetric solver is invoked automatically if you request coupled thermal-electrical analysis. For problems where coupling between the thermal and electrical solutions is weak or where a pure electrical conduction analysis is required for the entire model, the unsymmetric terms resulting from the interfield coupling may be small or zero. In these problems you can invoke the less costly symmetric storage and solution scheme by solving the thermal and electrical equations separately. The separated technique uses the symmetric solver by default. The thermal-electrical solution schemes are discussed below.

The theoretical basis of coupled thermal-electrical analysis is described in detail in “Coupled thermal-electrical analysis,” Section 2.12.1 of the ABAQUS Theory Manual.

Governing electric field equation

The electric field in a conducting material is governed by Maxwell's equation of conservation of charge. Assuming steady-state direct current, the equation reduces to

where

The flow of electrical current is described by Ohm's law:

where

is the electrical field intensity, defined as the negative of the gradient of the electrical potential ,

is the electrical potential,

)

is the electrical conductivity matrix,

is the temperature, and

are predefined field variables.

Using Ohm's law in the conservation equation, written in variational form, provides the governing equation of the finite element model:

where is the current density entering the control volume across

The electrical conductivity, , can be isotropic, orthotropic, or fully anisotropic (see “Electrical conductivity,” Section 20.6.1). Ohm's law assumes that the electrical conductivity is independent of the electrical field, . The coupled thermal-electrical problem is nonlinear when the electrical conductivity depends on temperature.

Specifying the amount of thermal energy generated due to electrical current

Joule's law describes the rate of electrical energy, , dissipated by current flowing through a conductor as

The amount of this energy released as internal heat within the body is , where is an energy conversion factor. You specify in the material definition. It is assumed that all the electrical energy is converted into heat () if you do not include the joule heat fraction in the material description. The fraction given can include a unit conversion factor, if required.

Input File Usage: | *JOULE HEAT FRACTION |

ABAQUS/CAE Usage: | Property module: material editor: |

Steady-state analysis

Steady-state analysis provides the steady-state solution directly. Steady-state thermal analysis means that the internal energy term (the specific heat term) in the governing heat transfer equation is omitted. Only direct current is considered in the electrical problem, and it is assumed that the system has negligible capacitance. (Electrical transient effects are so rapid that they can be neglected.)

Input File Usage: | *COUPLED THERMAL-ELECTRICAL, STEADY STATE |

ABAQUS/CAE Usage: | Step module: |

A steady-state analysis has no intrinsic physically meaningful time scale. Nevertheless, you can assign a “time” scale to the analysis step, which is often convenient for output identification and for specifying prescribed temperatures, electrical potential, and fluxes (heat flux and current density) with varying magnitudes. Thus, when steady-state analysis is chosen, you specify a “time” increment and a “time” period for the step; ABAQUS/Standard then increments through the step accordingly.

Any fluxes or boundary condition changes to be applied during a steady-state step should be given using appropriate amplitude references to specify their “time” variations (“Amplitude curves,” Section 27.1.2). If fluxes and boundary conditions are specified for the step without amplitude references, they are assumed to change linearly with “time” during the step—from their magnitudes at the end of the previous step (or zero, if this is the beginning of the analysis) to their newly specified magnitudes at the end of this step (see “Procedures: overview,” Section 6.1.1).

Transient analysis

Alternatively, the thermal portion of the coupled thermal-electrical problem can be considered transient. As in steady-state analysis, electrical transient effects are neglected. See “Uncoupled heat transfer analysis,” Section 6.5.2, for a more detailed description of the heat transfer capability in ABAQUS/Standard.

Input File Usage: | *COUPLED THERMAL-ELECTRICAL |

ABAQUS/CAE Usage: | Step module: |

Time integration in the transient heat transfer problem is done with the same backward Euler method used in uncoupled heat transfer analysis. This method is unconditionally stable for linear problems. You can specify the time increments directly, or ABAQUS can select them automatically based on a user-prescribed maximum nodal temperature change in an increment. Automatic time incrementation is generally preferred.

The time increment size can be selected automatically based on a user-prescribed maximum allowable nodal temperature change in an increment, . ABAQUS/Standard will restrict the time increments to ensure that these values are not exceeded at any node (except nodes with boundary conditions) during any increment of the analysis (see “Time integration accuracy in transient problems,” Section 7.2.4).

Input File Usage: | *COUPLED THERMAL-ELECTRICAL, DELTMX= |

ABAQUS/CAE Usage: | Step module: |

If you select fixed time incrementation and do not specify , fixed time increments equal to the user-specified initial time increment, , will then be used throughout the analysis.

Input File Usage: | *COUPLED THERMAL-ELECTRICAL |

ABAQUS/CAE Usage: | Step module: |

In transient heat transfer analysis with second-order elements there is a relationship between the minimum usable time increment and the element size. A simple guideline is

where is the time increment, is the density,

There is no upper limit on the time increment size (the integration procedure is unconditionally stable) unless nonlinearities cause convergence problems.

By default, a transient analysis will end when the specified time period has been completed. Alternatively, you can specify that the analysis should continue until steady-state conditions are reached. Steady state is defined by the temperature change rate; when the temperature changes at a rate that is less than the user-specified rate (given as part of the step definition), the analysis terminates.

Input File Usage: | Use the following option to end the analysis when the time period is reached: |

*COUPLED THERMAL-ELECTRICAL, END=PERIOD (default) Use the following option to end the analysis based on the temperature change rate: *COUPLED THERMAL-ELECTRICAL, END=SS |

ABAQUS/CAE Usage: | Step module: |

Fully coupled solution schemes

ABAQUS/Standard offers an exact as well as an approximate implementation of Newton's method for coupled thermal-electrical analysis.

An exact implementation of Newton's method involves a nonsymmetric Jacobian matrix as is illustrated in the following matrix representation of the coupled equations:

where and are the respective corrections to the incremental electrical potential and temperature, are submatrices of the fully coupled Jacobian matrix, and and are the electrical and thermal residual vectors, respectively.

Solving this system of equations requires the use of the unsymmetric matrix storage and solution scheme. Furthermore, the electrical and thermal equations must be solved simultaneously. The method provides quadratic convergence when the solution estimate is within the radius of convergence of the algorithm. The exact implementation is used by default.

Some problems require a fully coupled analysis in the sense that the electrical and thermal solutions evolve simultaneously, but with a weak coupling between the two solutions. In other words, the components in the off-diagonal submatrices are small compared to the components in the diagonal submatrices . For these problems a less costly solution may be obtained by setting the off-diagonal submatrices to zero, so that we obtain an approximate set of equations:

As a result of this approximation the electrical and thermal equations can be solved separately, with fewer equations to consider in each subproblem. The savings due to this approximation, measured as solver time per iteration, will be of the order of a factor of two, with similar significant savings in solver storage of the factored stiffness matrix. Further, in situations without strong thermal loading due to cavity radiation, the subproblems may be fully symmetric or approximated as symmetric, so that the less costly symmetric storage and solution scheme can be used. The solver time savings for a symmetric solution is an additional factor of two. Unless you explicitly select the unsymmetric solver for the step (“Procedures: overview,” Section 6.1.1), the symmetric solver will be used with this separated technique.

This modified form of Newton's method does not affect solution accuracy since the fully coupled effect is considered through the residual vector at each increment in time. However, the rate of convergence is no longer quadratic and depends strongly on the magnitude of the coupling effect, so more iterations are generally needed to achieve equilibrium than with the exact implementation of Newton's method. When the coupling is significant, the convergence rate becomes very slow and may prohibit the attainment of a solution. In such cases the exact implementation of Newton's method is required. In cases where it is possible to use this approximation, the convergence in an increment will depend strongly on the quality of the first guess to the incremental solution, which you can control by selecting the extrapolation method used for the step (see “Procedures: overview,” Section 6.1.1).

Input File Usage: | Use the following option to specify a separated solution scheme: |

*SOLUTION TECHNIQUE, TYPE=SEPARATED |

ABAQUS/CAE Usage: | Step module: |

Uncoupled electric conduction and heat transfer analysis

The coupled thermal-electrical procedure can also be used to perform uncoupled electric conduction analysis for the whole model or just part of the model (using coupled thermal-electrical elements). Uncoupled electrical analysis is available by omitting the thermal properties from the material description, in which case only the electric potential degrees of freedom are activated in the element and all heat transfer effects are ignored. If heat transfer effects are ignored in the entire model, you should invoke the separated solution technique described above. Use of this technique will then invoke the symmetric storage and solution scheme, which is an exact representation of a purely electrical problem.

Similarly, coupled thermal-electrical elements can be used in an uncoupled heat transfer analysis (“Uncoupled heat transfer analysis,” Section 6.5.2), in which case all electric conduction effects are ignored. This feature is useful if a thermal-electrical analysis is followed by a pure heat conduction analysis. A typical example is a welding process, where the electric current is applied instantaneously, followed by a cooldown period during which no electrical effects need to be considered. The symmetric solver is activated by default in an uncoupled heat transfer analysis.

Cavity radiation

Cavity radiation can be activated in a coupled thermal-electrical step. When the thermal emissivity is a function of temperature or field variables, you can specify the maximum allowable emissivity change during an increment in addition to the maximum temperature change (discussed above in “Automatic incrementation”) to control the time incrementation. See “Cavity radiation,” Section 32.1.1, for more information.

Input File Usage: | Use the following option in the step definition to activate cavity radiation: |

*RADIATION VIEWFACTOR Use the following option to specify the maximum allowable emissivity change: *COUPLED THERMAL-ELECTRICAL, MXDEM= |

ABAQUS/CAE Usage: | Cavity radiation is not supported in ABAQUS/CAE. However, you can specify the maximum allowable emissivity change for a coupled thermal-electrical step. |

Step module: |

By default, the initial temperature of all nodes is zero. You can specify nonzero initial temperatures or field variables (see “Initial conditions,” Section 27.2.1). Since only steady-state electrical currents are considered, the initial value of the electrical potential is not relevant.

Boundary conditions can be used to prescribe the electrical potential, (degree of freedom 9), and the temperature, (degree of freedom 11), at the nodes. See “Boundary conditions,” Section 27.3.1.

Boundary conditions can be specified as functions of time by referring to amplitude curves (see “Amplitude curves,” Section 27.1.2).

A boundary without any prescribed boundary conditions corresponds to an insulated surface.

Both thermal and electrical loads can be applied in a coupled thermal-electrical analysis.

The following types of thermal loads can be prescribed in a coupled thermal-electrical analysis, as described in “Thermal loads,” Section 27.4.4:

Concentrated heat fluxes

Body fluxes and distributed surface fluxes

Convective film conditions and radiation conditions

The following types of electrical loads can be prescribed:

Concentrated current

Distributed current densities, (see Part VI, “Elements,” for the current density types available with particular elements)

The magnitude of a concentrated current or distributed current density can be controlled by referring to an amplitude curve (see “Amplitude curves,” Section 27.1.2). If different magnitude variations are needed for different currents, the current definitions can be repeated, with each referring to its own amplitude curve.

Predefined temperature fields are not allowed in coupled thermal-electrical analyses. Boundary conditions should be used instead to specify temperatures, as described above.

Other predefined field variables can be specified in a coupled thermal-electrical analysis. These values affect only field-variable-dependent material properties, if any. See “Predefined fields,” Section 27.6.1.

Both thermal and electrical properties are active in coupled thermal-electrical analyses. If thermal properties are omitted, an uncoupled electrical analysis will be performed.

All mechanical behavior material models (such as elasticity and plasticity) are ignored in a coupled thermal-electrical analysis.

For the heat transfer portion of the analysis, the thermal conductivity must be defined (see “Conductivity,” Section 20.2.2). The specific heat must also be defined for transient heat transfer problems (see “Specific heat,” Section 20.2.3). If changes in internal energy due to phase changes are important, latent heat can be defined (see “Latent heat,” Section 20.2.4). Thermal expansion coefficients (“Thermal expansion,” Section 20.1.2) are not meaningful in a coupled thermal-electrical analysis since deformation of the structure is not considered. Internal heat generation can be specified (see “Uncoupled heat transfer analysis,” Section 6.5.2).

For the electrical portion of the analysis, the electrical conductivity must be defined (see “Electrical conductivity,” Section 20.6.1). The electrical conductivity can be a function of temperature and user-defined field variables. The fraction of electrical energy dissipated as heat can also be defined, as explained above.

The simultaneous solution in a coupled thermal-electrical analysis requires the use of elements that have both temperature (degree of freedom 11) and electrical potential (degree of freedom 9) as nodal variables. The finite element model can also include pure heat transfer elements (so that a pure heat transfer analysis is provided for that part of the model) and coupled thermal-electrical elements for which no thermal properties are given (so that a pure electrical conduction solution is provided for that part of the model).

Coupled thermal-electrical elements are available in ABAQUS/Standard in one dimension, two dimensions (planar and axisymmetric), and three dimensions. See “Choosing the appropriate element for an analysis type,” Section 21.1.3.

The following output variables can be used to request output relating to the electric conduction solution:

Element integration point variables:

EPG | Magnitude and components of the electrical potential gradient vector, . |

EPGM | Magnitude of the electrical potential gradient vector. |

EPG n | Component |

ECD | Magnitude and components of the electrical current density vector, |

JENER | Electrical energy dissipated due to flow of current, . |

Whole element variables:

ECURS | Distributed applied electrical current. |

NCURS | Electrical current at nodes due to electric conduction. |

ELJD | Total electrical energy dissipated due to flow of current, . |

Nodal variables:

EPOT | Electrical potential, . |

RECUR | Reactive electrical current. |

CECUR | Concentrated applied electrical current. |

Whole model variables:

ALLJD | Electrical energy summed over the model. |

Surface interaction variables (see “Electrical contact properties,” Section 30.3.1):

ECD | Electrical current density. |

ECDA | ECD multiplied by area. |

ECDT | Time integrated ECD. |

ECDTA | Time integrated ECDA. |

SJD | Heat flux per unit area generated by the electrical current. |

SJDA | SJD multiplied by area. |

SJDT | Time integrated SJD. |

SJDTA | Time integrated SJDA. |

WEIGHT | Heat distribution between interface surfaces, |

*HEADING … *MATERIAL, NAME=mat1*CONDUCTIVITYData lines to define thermal conductivity*ELECTRICAL CONDUCTIVITYData lines to define electrical conductivity*JOULE HEAT FRACTIONData lines to define the fraction of electric energy released as heat** *STEP *COUPLED THERMAL-ELECTRICALData line to define incrementation and steady state*BOUNDARYData lines to define boundary conditions on electrical potential and temperature degrees of freedom*CECURRENTData lines to define concentrated currents*DECURRENT and/or *DSECURRENTData lines to define distributed current densities*CFLUX and/or *DFLUX and/or *DSFLUXData lines to define thermal loading*FILM and/or *SFILM and/or *RADIATE and/or *SRADIATEData lines to define convective film and radiation conditions… *CONTACT PRINT or *CONTACT FILEData lines to request output of surface interaction variables*END STEP