We present and discuss a set of boundary conditions (BCs) to use in three-dimensional, mesoscopic, finite element models of mid-infrared pulsed laser ablation of brittle materials. These models allow the study of the transient displacement and stress fields generated at micrometer scales during and after one laser pulse, where using conventional BCs may lead to some results without physical significance that can be considered an artefact of the calculations. The proposed BCs are tested and applied to a micrometer-scale continuous model of human dental enamel under CO2 radiation (10.6 mu m, 0.35 mu s pulse, sub-ablative fluence), giving rise to the following results: the highest stress is obtained at the irradiated surface of the model, at the end of the laser pulse, but afterwards it decreases rapidly until it becomes significantly lower than the stress in a region 2.5 mu m deep in the model; a thermally induced vibration in the material is predicted. This non-intuitive dynamics in stress and displacement distribution cannot be neglected and has to be considered in dynamic laser ablation models, since it may have serious implications in the mechanisms of ablation. (c) 2005 Elsevier B.V. All rights reserved.