Description: Coupled temperature–pore pressure elements are similar to existing pore pressure elements with the important difference that they have temperature as a nodal degree of freedom in addition to pore pressure and displacement. The element formulation solves the heat transfer equation in addition to and in a fully coupled manner with mechanical equilibrium and continuity equations. The formulation is important in situations where there is a strong coupling between pore fluid flow and heat transfer and the temperature field has a relatively strong dependence on the stresses and the pore pressure. For example, changes in the pore pressure may affect the flow rates, which in turn may affect the temperature distribution due to the modified convection rates. Although these elements assume by default that nodal temperature is an active degree of freedom, they can be used in a pure consolidation analysis without heat transfer effects by turning off the temperature degrees of freedom. The element library includes first-order brick, second-order modified tetrahedron, and first-order axisymmetric formulations. The thermal material properties that govern the heat transfer behavior can be independently specified for both the solid and the pore fluid phases.