We develop conservation laws for coupled hydro-mechanical processes in unsaturated porous media using three-phase continuum mixture theory. From the first law of thermodynamics, we identify energy-conjugate variables for constitutive modeling at macroscopic scale. Energy conjugate expressions identified relate a certain measure of effective stress to the deformation of the solid matrix, the degree of saturation to the matrix suction, the pressure in each constituent phase to the corresponding intrinsic volume change of this phase, and the seepage forces to the corresponding pressure gradients. We then develop strong and weak forms of boundary-value problems relevant for 3D finite element modeling of coupled hydro-mechanical processes in unsaturated porous media. The paper highlights a 3D numerical example illustrating the advances in the solution of large-scale coupled finite element systems, as well as the challenges in developing more predictive tools satisfying the basic conservation laws and the observed constitutive responses for unsaturated porous materials.