We present a systematic derivation of hydrodynamic theories for nonhomogeneous nematic liquid crystal polymers (LCPs) by approximating the molecules as rigid ellipsoids, which can be either uniaxial molecules (spheroids) or biaxial ones. The short range interaction is assumed to be dominated by the excluded volume effect. Additional molecular properties with ellipsoidal molecules, e.g., a dipole-dipole interaction in extended nematics and chiral molecular structure in cholesterics, are accounted for through additional intermolecular potentials. Long-range molecular interaction is implemented through an averaged mean-field potential characterized by interaction functions. The extra elastic stress tensor is calculated using an extended virtual work principle consistent with conservation of angular momentum on the material volume, whereas the extra viscous stress is obtained by Batchelor's volume averaging method. In the isothermal case, the theories are shown to satisfy the second law of thermodynamics, i.e., they admit positive production of entropy or energy dissipation. In the case of cholesterics, the kinetic theory reduces to the Leslie-Ericksen theory in the limit of weak translational diffusion, weak long range interaction, and weak flow. (C) 2004 Elsevier B.V. All rights reserved.