This work summarizes and generalizes the results of long term researches of mathematic modeling of stress-strain state (SSS) of shell molds (SM) considering their morphological structure and types of hydromechanical and temperature impacts on them. A problem for determination of SSS of multilayer SM with internal porous layer during filling this mold with melt and mold hardening considering the phase changes was set and solved. At the same time the calculation methodology described in the early works of authors, namely constitutive equation system, numerical method, numerical scheme and solution algorithm was used. The results of numerical solution of the problem, namely diagrams of stresses arising in the SM during casting hardening are given here.