The elastoplastic constitutive equation with the damage is formulated incorporating the subloading surface model. Further, it is extended to describe the unilateral damage phenomenon by formulating the actual Young’s modulus tensor as the function of the signs of the principal actual damaged stresses and the damage variable. Here, we may perform the ordinary deformation analysis simply in the virtual undamaged state
