For wedge splitting test specimens, the stress and displacement fields both in the vicinity and also in larger distance from the crack tip are investigated by means of numerical methods. Several variants of boundary conditions were modeled. The stress intensity factor K, T-stress and even higher-order terms of William series were determined and subsequently utilized for analytical approximation of the stress field. A good fit between the analytical and numerical solution in dependence on the distance from the crack tip was shown, compared and discussed. Presented approach is considered as suitable for estimation of the fracture process zone extent in silicate composite materials.