Laser-generated elastic waves have been the subject of numerous experimental, theoretical, and numerical studies to describe the opto-acoustic generation process, involving electromagnetic, thermal, and elastic fields and their couplings in matter. Among the numerical methods for solving this multiphysical problem, the semi-analytic approach is one of the most relevant for obtaining fast and accurate results, when analytic solutions exist. In this paper, a multilayer model is proposed to successively solve electromagnetic, thermal, and elastodynamic problems. The optical penetration of the laser line source, as well as thermal conduction and convection, are accounted for. Optical, thermal, and mechanical coupling conditions are considered between the upper and lower media of the multilayer. The simulation of laser-generated ultrasounds in multilayer structures is of interest for the development of nondestructive evaluation methods of complex structures, such as bonded assemblies in aeronautics [as discussed in Hodé et al., J. Acoust. Soc. Am. 150, 2076 (2021)]. The developed Python code is provided for free at https://doi.org/10.5281/zenodo.4301720.