INTRODUCTION
Major innovations have taken place in implant manufacturing technology for use in cosmetic surgery, followed by an increase in breast cancer reconstructions and the need for aesthetic breast surgery [
1]. Although the shape of the breast changes over time during the postoperative period, it is not clear whether these changes depend on the shape of the implant initially implanted [
2].
In this study, we investigated how the dimensions of the breast implant, including its projection and height, affect the shape of the breast in the long term. We also investigated the ideal implant dimensions for optimizing breast shape from an aesthetic perspective. In addition to the effects of the dimensions of the breast implant on the morphology of the breast, we analyzed the internal stress on the breast during static and dynamic states after implantation.
In order to analyze the postoperative external shape and internal stress on the breast, we simulated an
in situ operative breast and a postoperative breast two-dimensionally using a finite element model, which is a widely used method of computational analysis [
3].
IDEA
Finite element modeling
For two-dimensional finite element analysis (FEA) of the external and internal characteristics of the breast with respect to the initial shape of the breast and its shape after implantation, commercial FEA software was utilized (ABAQUS version 6.12; Dassault Systèmes, Vélizy-Villacoublay, France).
First, as shown in
Fig. 1, we characterized the initial shape of a patient’s breast and implant in terms of shape parameters such as the projection, top-point, and height. Among the shape parameters, the top-point was selected as the main design parameter because the projection and the arc length change depending on the position of the top-point. Two-dimensional geometric models of the patient’s breast and implant were prepared. Additionally, the rigid breast plate, which pushes the implant forward into the breast, was introduced to simulate the process of implant surgery.
Based on the geometric model, a 3-node linear plane stress triangle element (CPS3) was assigned to both the breast and the implant. The finite element model of the implant comprised 430 elements and 245 nodes. For the breast, the finite element model comprised 958 elements and 552 nodes. A 2-node two-dimensional linear rigid link element (R2D2) was assigned to the rigid breast plate. The finite element model of the rigid breast plate comprised 16 elements and 17 nodes.
The implant consisted of a shell and filler, made of polydimethylsiloxane (PDMS) and cohesive silicone gel, respectively. Anatomically, the breast mainly consists of fat tissue and fibroglandular tissue. Generally, compared with implants, breast tissue is much softer and more stretchable, like a rubber material. Thus, the constitutive equations for the breast and implant were derived based on the Mooney-Rivlin [
4] model and elasticity model [
5], respectively.
Loading and boundary conditions
Static analysis of the breast with the implant
Initially, the implant and the breast were separated from each other, so that the bottom surfaces of the implant were aligned in a direction parallel to the bottom surfaces of the breast without surface overlap. The rigid plate was also separated and placed at the rear of the implant, and aligned parallel to the bottom surface of the implant, as shown in
Fig. 1.
Dynamic analysis of the breast with the implant
Following the static analysis for implantation, the dynamic response of the postoperative breast was estimated by simulating the periodic vibration of the implant-embedded breast, as shown in
Fig. 2. For this simulation, the loading and boundary conditions were modified subsequent to implantation.
Static analysis
From the nonlinear static analysis, the extrinsic and intrinsic changes of the breast were analyzed numerically in scenarios involving implants with different top-points. As shown in
Fig. 3, the central cross-sectional location of the top-point of the implant-embedded breast tended to increase upward as the top-point of the implant became lower than the initial top-point of the breast. In contrast, as the top-point became higher, a downward migration was observed.
Dynamic analysis
The maximum stress alternately appeared in the lower chest and the upper chest during up-and-down dynamic behavior. However, the value of the maximum stress was generally about 1.6 to 2 times greater than the stress level in the static state, as shown in
Fig. 4.
DISCUSSION
Implant-based breast augmentation surgery was first performed more than 40 years ago [
6]. Since then, silicone implant technology has been continuously developing, resulting in innovations such as better materials and designs, and the ability to manufacture implants with various dimensions. Although implants with diverse dimensions and designs are used in everyday practice, there have not been any long-term follow-up studies or simulations of morphological changes with different implant shapes.
This study quantitatively estimated the external changes affecting postoperative breast aesthetics according to the height of the breast, the breast implant, and the relative position of the protrusion. In particular, internal characteristics such as mechanical stress showed significant quantitative differences depending on the shape of the breast implant and its location in relation to the chest. In particular, stress was maximized at the lower part of the chest, which is causally related to the clinical manifestation of inframammary fold tear, which most often occurs in the lower chest.
In the present study, it was shown that the mechanical stress applied to the lower chest region varied according to the maximal projection point of the breast implant. The maximum stress value received by the lower chest area was 1.6 to 2 times larger during dynamic motion than in static conditions.
The study has several limitations. First, it was performed using two-dimensional analysis, rather than a three-dimensional simulation, which would have been helpful for a more precise analysis. However, the most intriguing aspect of this study is the change in the central cross-sectional shape between the breast implant and the soft tissue envelope of the breast in the sagittal section after insertion of the implant. Based on the fact that the compression under the human chest after implant insertion (referred to as stress in this study) is also greatest near the bottom of the central section of the breast implant, we hypothesize that a two-dimensional study in a central cross-sectional plane may be more accurate for a static and dynamic analysis of pressure measurements in the central section.
Additionally, the location of the inframammary fold was fixed at a single point, rather than migrating in the downward direction, as in actual practice. However, because the goal of this study was to establish a boundary using a reference point for rigid plate alignment in terms of FEA, and to perform quantitative measurements using stress analysis on this boundary, it was not also possible to measure changes in the position of the inframammary fold, for which additional, separately designed research would be required.