SPICE is inadequate for this job because the new crystals have much larger nonlinear sensitivities. You need a general purpose model that starts from the material anisotropic stress/strain elasticity tensor and its dependence on a host of other factors such as drive level as basis and goes from there. The usual applications engineering linear time invariant model with fixed components is a joke.