Accurately predicting the forces driving crack propagation is critical for engineering applications where fracture is an important failure criterion. A generalized approach to unifying crack initiation as well as the force available to drive the crack in different directions is through the theory of configurational force. The configurational force provides the energy required to change the surface configuration at a phase interface or crack. We demonstrate the calculation of configurational force within an isogeometric computational framework for fracture modeling. The implementation relies on a recently developed enrichment procedure for multi-material corner singularities, of which cracks are a special case. Enriched displacement and temperature functions applied at crack tips are illustrated and demonstrated through several numerical examples on thermoelastic solids. Examples include (a) crack driving force estimation under adiabatic/isothermal conditions and (b) edge crack propagation simulation.