I agree with Qasim's recommendation - *MAT110. You can find an LS-DYNA user conference paper on the application of this model (in addition to parameters) for glass, which may be a good starting point....
If you have a license of COMSOL all you need to do is to type in the equation that defines the constitutive law or perturbs the linear one. If you intention is to capture material nonlinearity on a mainly linear material because of microcracking then you best bet is to weaken the elements of the Hooke's tensor as a function of strain. You can do this in both a the usual volumetric sense by weakening the diagonal elements of the tensor (that would be the most reasonable for a ceramic material) or in the distortional sense by weakening the non-diagonal terms. In COMSOL this is trivial. In other codes like ANSYS and LS-DYNA and ABAQUS you need to select a material with damage (isotropic) with a value of the damage parameter very near 0.0
Hi John, what do you mean with "you need to select a material with damage (isotropic) with a value of the damage parameter very near 0.0". I am trying to define an elastic perfectly brittle material with Ansys but I don't know how. Can you please explain more? Thank you.