Abstract
Skin meshing facilitates the greater expansion of donor skin through patterned slits and is widely used for treating extensive burn injuries. However, the actual expansion often falls below manufacturers' claims. Previous computational analyses using the isotropic Yeoh model have shown that Langer's line orientation and slit direction significantly affect induced stress and meshing ratios, yet the use of nonlinear anisotropic models that represent collagen fiber alignment corresponding to Langer's lines remains unexplored. This study employs a nonlinear anisotropic Gasser-Ogden-Holzapfel (GOH) model with slit orientations of 0°, 45°, and 90°, consistent with geometries reported in the literature, to quantify induced stress in skin meshing by incorporating collagen fibers within the dermis layer. The GOH parameters were calibrated to human back skin data uniaxially stretched parallel and perpendicular to Langer's lines using Levenberg-Marquardt optimization in the GIBBON toolbox (MATLAB R2023a) coupled with FEBio v4.0, achieving a standard deviation of 3% relative to experimental data. The GOH model predicted the highest induced stress at 100% strain for the 45° slit parallel to Langer's lines and the lowest for the 90° slit perpendicular, exceeding 40 MPa due to absence of damage and rupture modeling but accurately representing anisotropic mesh behavior.