<p>We present a fully coupled hydro-thermomechanical framework for multiphase model of porous media under extreme thermomechanical conditions, possibly involving extremely large deformation, heat conduction, phase transition and internal flow. The proposed computational framework combines the hot optimal transportation meshfree (HOTM) method and a meshfree formulation of Darcy’s law. In specific, the optimal transportation theory is introduced for temporal discretization, while material-point sampling method is employed for spatial discretization of the porous media. The linear momentum and energy conservation are formulated in the Lagrangian configuration and jointly solved in the HOTM framework to predict the solid-skeleton deformation and the temperature evolution in porous media. Meanwhile, the mass conservation and Darcy’s law are formulated in the current configuration at the material-point level and solved via a weighted residual method to predict the internal fluid flow and the porosity distribution in porous media. The detailed formulation of the computational framework is stated and validated. Then, we focus on a particular application: the hot-forming process of resin-based friction composites. In this simulation, the resin-based matrix is modeled as a continuous porous medium, while particles and fibers are modeled as explicit spheres and cylinders embedded in the porous matrix. Simulations with various loading conditions are conducted to investigate the effects of loading parameters on the internal fluid flow and porosity in the product. The sensitivity of product’s porosity on loading conditions including pressure and temperature is further studied.</p>