- Research article
- Open Access
Finite element analysis of posterior acetabular column plate and posterior acetabular wall prostheses in treating posterior acetabular fractures
Journal of Orthopaedic Surgery and Research volume 18, Article number: 94 (2023)
The purpose of this study was to investigate the mechanical stability of the posterior acetabular column plate and different posterior acetabular wall prostheses used in treating posterior acetabular fractures with or without comminution.
The unilateral normal ilium was reconstructed, and a model of posterior acetabular wall fracture was established on this basis. The fracture fragment accounted for approximately 40% of the posterior acetabular wall. The posterior acetabular column plate and different posterior acetabular wall prostheses were also designed. Using static and dynamic analysis methods, we observed and compared the changes in the stress and displacement values of different models at different hip joint flexion angles under external forces.
At different hip flexion angles, the stress of each model mainly fluctuated between 37.98 MPa and 1129.00 MPa, and the displacement mainly fluctuated between 0.076 and 6.955 mm. In the dynamic analysis, the nodal stress‒time curves of the models were nonlinear, and the stress changed sharply during the action time. Most of the nodal displacement‒time curves of the models were relatively smooth, with no dramatic changes in displacement during the action time; additionally, most of the curves were relatively consistent in shape.
For simple posterior acetabular wall fractures, we recommend using a posterior acetabular column plate. In the case of comminuted posterior acetabular fractures, we recommend the use of a nonflanked posterior acetabular prosthesis or a biflanked posterior acetabular prosthesis. Regarding the method of acetabular prosthesis design, we propose the concept of “Break up to Make up” as a guide.
Posterior acetabular fractures are the most common type of acetabular fracture, accounting for approximately 15–47% of acetabular fractures [1,2,3,4], and are commonly caused by car accidents or falls from heights . Currently, reconstructing the integrity of the posterior wall of the acetabulum is difficult. Many complications may occur after posterior acetabular wall fracture, with the most common complications including traumatic hip arthritis, femoral head necrosis, hip instability and redislocation [6, 7], which can greatly increase the financial burden on patients. At present, for common posterior acetabular wall fractures, if the fracture is a simple posterior wall fracture without fragmentation, the original fragment is used for reduction. If the fracture fragment is comminuted, then autologous ilium reconstruction is applied; however, because this method cannot achieve complete reconstruction of the acetabular articular surface, the probability of postoperative hip osteoarthritis or femoral head necrosis is relatively high.
To solve the problem of reconstructing the posterior wall of the acetabulum, we designed a plate for the posterior column of the acetabulum, which is mainly intended for use in treating simple posterior wall fractures of the acetabulum. To solve the problems related to comminuted fractures of the posterior acetabular wall, we designed three types of posterior acetabular prostheses, with the aim of reconstructing the articular surface of the posterior acetabular wall and maintaining the anatomical reduction of the posterior acetabular wall to the greatest extent possible. The purpose of this study was to understand and compare the mechanical response of the above different methods for the internal fixation of posterior acetabular fractures and to select the best internal fixator according to the type of posterior acetabular fracture. We use the finite element method to compare the forces in different models at different hip flexion angles using static and dynamic analysis to understand the mechanical properties of the different models.
Computed tomography (CT) data of a volunteer’s hip were obtained and imported into Mimics 20.0 software (Materialise, Leuven, Belgium) to reconstruct a normal ilium. The acetabular fossa of the ilium was divided into 12 equal parts (clock model) and marked as 12 points, with the 12th point corresponding to the intersection of the line connecting the anteroinferior iliac spine and the ischial tuberosity with the top of the acetabulum, and the sixth point symmetrical to the 12th point symmetrical across the centre of the acetabulum. Each point was named according to the corresponding point on a clock. The posterior wall ranged from the 7th point to the 11th point, accounting for 120°; in this ilium model, the diameter of the acetabular fossa was 55.23 mm. Then, a fragment was created in the posterior wall of the acetabulum, accounting for approximately 40% of the posterior wall, i.e., 48°, as shown in Fig. 1. After the posterior acetabular wall fracture model was established, with or without comminution, different models for fixing the posterior wall of the acetabulum were established, including the following prosthesis models: a common reconstruction plate (CRP) model, consisting of two medial and lateral reconstruction plates for posterior wall fracture fixation; a posterior acetabular column plate (PACP) model, consisting of a main plate and two small upper and lower “T” plates on the same side; a nonflanked posterior acetabular wall prosthesis (NPAWP) model; a biflanked posterior acetabular wall prosthesis (BPAWP) model, consisting of two (upper and lower) lateral wings; and a triflanked posterior acetabular wall prosthesis (TPAWP) model, consisting of three (upper, middle and lower) lateral wings. The model of the ilium was imported into 3-matic to construct the models of the steel plates and prostheses. To reduce the complexity of the analysis, the model was simplified by ignoring the soft tissue around the hip joint and the femoral head. The final result is shown in Fig. 2. The above five models were then imported into HyperMesh 13.0 software (Altair Company, Troy, MI, USA) for meshing using four-node tetrahedral elements (C3D4). Then, the separate fracture model was imported into Mimics 20.0 software for material assignment. We applied titanium alloy material properties for the steel plates and prostheses and smoothly polished the hip joint surface of the prostheses. Finally, all models were imported into Abaqus 2021 software (Dassault; France) for analysis. Static and dynamic analyses were carried out to obtain numerical values for stress and displacement. In the dynamic analysis, plots of stress versus time and displacement versus time were also obtained.
The plates, screws and prostheses in this study were all assigned the properties of titanium alloy, with a Young modulus of 110 GPa and Poisson ratio of 0.316 . It was assumed that the materials of the internal fixators, prostheses and ilium were all isotropic. When assigning material properties to the iliac bone model, we followed our previous method , based on the CT value. According to the formula, different CT values correspond to different Young’s moduli; thus, the true material properties of the ilium can be approximated. In the static analysis, surface‒surface contact, hard contact tangent behaviour, a penalty contact attribute and a friction coefficient between bone and bone of 0.40 were applied ; surface‒surface contact, hard contact tangent behaviour, a penalty contact property and a friction coefficient between the internal fixators or prostheses and the ilium of 0.45 were applied . As boundary conditions, the pubic symphysis, articular surface of the ilium and sacroiliac joint were fixed. In the dynamic analysis, general contact (display dynamic analysis), hard contact tangent behaviour and a penalty contact attribute were applied. The mechanical analysis allowed only one type of contact to exist; thus, only the friction between the internal fixators or prostheses and the ilium was set to 0.45, and the friction between bone and bone was ignored.
Boundaries and loads
Loading was applied in two ways. In the static analysis, the five models were simulated in the standing position and at 30°, 60°, 90° and 120° of hip joint flexion and were subjected to a force of 2300 N ; the stress distribution of the internal fixators and the displacement of the models were determined under all these conditions. In the dynamic analysis, at 90° of hip flexion, the femoral head was simulated to impact the posterior wall of the acetabulum with the same force. The force increased from 0 to 2300 N in the first 0.5 ms and decreased to 0 in the last 0.5 ms. The total duration of force application was 1 ms , representing the mechanical response of contact with the posterior wall of the acetabulum. The relationships between the maximum stress and displacement nodes and times for the five model plates and screws or prostheses and screws were obtained.
To compare whether the stress and displacement of different models were significantly different, the Kruskal‒Wallis test method was used in this study, and P < 0.05 was considered statistically significant. If a significant difference between the models was calculated, pairwise comparisons were performed using a multivariate ANOVA method.
A CRP model (592,718 units, 133,777 nodes), an APCP (450,007 units, 100,879 nodes), an NPAWP model (466,968 units, 101,596 nodes), a BPAWP model (440,890 units, 97,965 nodes) and a TPAWP model (488,085 units, 108,440 nodes) were established.
In the static and dynamic analyses, the CRP, APCP, NPAWP, BPAWP and TPAWP models were placed in different positions, and the maximum stress in each position is shown in Table 1. The maximum stress of each plate and prosthesis model fluctuated between 37.98 and 1129.00 MPa. The maximum displacement is shown in Table 2. In different positions and under different conditions, the displacement of the models fluctuated between 0.076 mm and 6.955 mm. The static displacement results are similar to those of the in vitro biomechanical study reported by Altun et al. , supporting the credibility of the experimental results of this study.
The stress distributions of different models in the static and dynamic analyses are shown in Figs. 3 and 4, respectively. In this study, the stress distribution was affected by the degree of hip flexion as well as the force. As the hip flexed, the distribution of the maximum stress was mostly vertical, along the internal fixator or prosthesis around the acetabular joint. Additionally, stress was commonly found at the junction of the plate and screw hole, plate and bone surface, “T” plate and main plate, prosthesis and screw and prosthesis and bone surface.
The relationship between stress and action time during the dynamic analysis is shown in Fig. 5. The relationship between stress and time was irregular, with the relationship for each model showing a unique and irregular shape. However, many of the model stress‒time plots exhibited a sawtooth-like shape.
The relationship between displacement and action time during the dynamic analysis is shown in Fig. 6. The graphs of the displacement–time relationship for each model were essentially the same, indicating nonlinear relationships.
Statistical analysis using the Kruskal‒Wallis test showed a significant difference in stress among the models (P = 0.005 < 0.05) but not displacement (P = 0.792 > 0.05). When the displacement data were subjected to multivariate analysis of variance, the pairwise comparisons between different models indicated a statistically significant difference between the CRP model and the TPAWP model (P = 0.02 < 0.05). There were no significant differences between the other models.
The results of this study suggest that the stress of the plate and prosthesis models ranges between 37.98 and 1129.00 MPa. According to the statistical analysis results, the magnitude of the stress varied significantly among the models, mainly between the CRP and TPAWP models. There were other differences between the models. The stress values of the TPAWP model in various positions were larger than those of the other models, indicating that the design of that prosthesis was more likely to cause stress concentration. If this prosthesis was subjected to the complex mechanical environment in vivo, the design may lead to an increased risk of internal fixation failure. In this study, titanium material properties were applied, and the yield strength of titanium is between 800 and 1000 MPa . Applied forces exceeding this range will likely cause internal fixation fracture. According to our research results, when the hip joint was flexed at 30°, 60° and 120°, the stress of the acetabular prosthesis models fell within this range or exceeded this range, indicating that when the external force is 2300 N, the acting time of the force is 1 ms, and the hip joint is flexed at one of the above angles, loosening of or damage to the prosthesis may easily occur. We recommend avoiding strenuous movements with the hip in flexion when using this type of prosthesis. The research of the static and dynamic analyses of the hip joint showed that under different degrees of hip flexion, the displacement fluctuated from 0.076 to 6.953 mm, and statistical analysis revealed no significant differences. This finding indicates that the models in this study preliminarily meet the requirements of internal fixation stability under the applied load type.
The maximum stress of the models was mainly observed at the edge of the plate, at the edge of the screw holes, at the junction between the “T”-shaped components and the main plate and at the junction between the flank and the bone surface. The common features of areas of stress concentration were sudden changes in the geometric shape and cuts in the material of the components . On this basis, repeated loading during activities of daily life may cause fatigue damage . Therefore, to improve the design of our plates and prostheses, some additional processing can be performed in places where the shape suddenly changes, such as the application of some circular chamfering structures to serve as a transitional structure between different shapes and mitigate stress concentration.
Most of the stress‒time curves exhibited a sawtooth-like structure, indicating that after the node was acted on by an external force, the stress of the node changed violently and nonlinearly. The reason for this may be that the dynamic behaviour of the same material is unstable and varies with the loading parameters . Moreover, the output of this study is represented in terms of changes at the nodes of the models, and these changes may be more obvious at the node level than at the model component level. In the displacement‒time relationship diagrams, the changes at the nodes were relatively smooth. Except for the curve of the prosthesis model (e.g., the “CRP model”) being close to linear, the curves of other models were not very different and were relatively smooth, showing an “S” shape.
The stress and displacement results indicate that the stress of the prostheses is relatively concentrated relative to the stress of the plate and that the stress of the TPAWP model is larger than that of the other models, with a significant difference compared with the CRP model. We preliminarily conclude that for internal fixation stability and surgical success, other internal fixation methods should be given priority. For example, if there is a noncomminuted fracture of the posterior acetabular wall, the PACP should be preferred because the plate is anatomically designed to avoid the need to prebend the plate during the operation, which would require additional operation time and increase the risk of the operation. In addition, the PACP features two “T”-shaped plates to reinforce the posterior wall fracture fragments, and each “T”-shaped plate has screw holes (the actual product has 3 screw holes; the model was simplified in this study) and can be fixed with screws as needed. Both the NPAWP and the BPAWP can be used to reconstruct the anatomy of the posterior acetabular wall. However, the maximum stress demonstrated by the NPAWP in this study was 1050.35 MPa, while the maximum stress demonstrated by the BPAWP was 1028.63 MPa, indicating that the BPAWP may have better mechanical stability than NPAWP. Structurally, the NPAWP has a simple design and allows a good mechanical effect to be achieved with only two screws. Under the premise of avoiding strenuous exercise after surgery, to save operation time, the NPAWP is recommended.
In actual clinical practice, to create a posterior acetabular prosthesis that fits the shape of the fractured posterior wall of the acetabulum to the greatest extent, we recommend individualizing the design of the prosthesis and printing it; additionally, an osteotomy guide plate should be designed at the same time to appropriately remove a part of the bone such that the prosthesis fits the posterior wall of the acetabulum to the greatest extent. Furthermore, the prosthesis should be designed with a porous structure , and bone tissue engineering strategies can be applied to create a scaffold . The purpose of these strategies is to maximize the preservation of bone mass and facilitate future hip replacement. Under the premise of strength, the pore size needs to be maximized for autologous bone grafting to retain sufficient bone mass. At the same time, to improve the design of the prosthesis, absorbable metals can be used  in the fabrication process. However, it is necessary to ensure that absorbable material remains in the body for a sufficiently long period. When reconstructing the acetabular articular surface, to prevent the prosthesis from rubbing against the femoral head, the articular surface of the prosthesis needs to be polished and the edges chamfered. Alternatively, on the premise of ensuring the strength of the prosthesis, the articular surface of the prosthesis can be designed with a thin layer of titanium and a baffle at the rear of the prosthesis to support the prosthesis and prevent the femoral head from affecting the prosthesis under physiological conditions. Such a baffle structure can dissipate dangerous kinetic energy . The posterior part of the prosthesis can consist of a frame structure, which is filled in with grafted bone. In this case, several baffle-like structures are needed, and screws can be passed through the flank structures to fix the prosthesis. The main advantages of these prosthesis design concepts are as follows: anatomically reconstructing the posterior wall of the acetabulum; maintaining the stability of the acetabulum; maintaining the bone volume of the posterior wall via the frame structure; and avoiding wear of the femoral head by polishing the articular surface. Li et al.  claimed that reconstructing an acetabular with large defects was a challenge during surgery. In some complex cases, the key to success was how to attach the fixations to the remaining autologous bone. The advantage of individualized acetabular components was that they could stably cover large acetabular defects and, at the same time, they matched the implant site. The prosthesis in this study is equivalent to a partial acetabular component; maximizing the preservation of the remaining bone after placement is key. Based on Lagoa et al. , they developed a 3D printed repair material for bone defects that were elastic to avoid stress occlusion, had a supporting effect, were biocompatible and could be partially biodescended. Combined with our research, these materials can be used to reconstruct the posterior acetabular wall and maintain bone mass in the posterior acetabular wall, and the properties of partially degradable materials can also be appropriate for patients who will need an artificial hip replacement in the future. Thus, the new concept of posterior acetabular wall reconstruction proposed in this study can be summarized as follows: Break up to Make up.
The prognosis for posterior acetabular fractures is unsatisfactory. Currently, it is common to use autologous iliac bone grafting to rebuild the posterior wall of the acetabulum. This procedure increases the surgical incisions, the amount of intraoperative bleeding and the operation time. Autologous iliac bone grafting sometimes fails to completely dissect and reconstruct the posterior acetabular wall, resulting in joint surface spaces or steps. Uneven joint surfaces can lead to postoperative traumatic arthritis , and eventually the patient undergoes artificial hip replacement or other secondary surgery. In a retrospective study by Iselin et al. , 10 of 33 patients (30%) with acetabular fracture involving the posterior wall developed symptoms of osteoarthritis, of whom 6 (18%) required total hip replacement or hip arthropexy. The personalized acetabular posterior wall prosthesis, polished articular surface and degradable material properties used in this study have the potential to be another ideal treatment modality.
In this study, calculations were performed only on a hemipelvic model, and analysis of the whole pelvis was not performed. This study assumed that the properties of bone materials were isotropic, ignoring other properties of bone, which will inevitably result in differences from actual values. In addition, the acetabulum of the human body is a complex component, including ligaments, the joint capsule, synovial fluid and surrounding soft tissue structures, which will have a certain impact on the force of the acetabulum. This study consisted of only a numerical simulation, which failed to simulate the friction of the acetabular prosthesis on the femoral head. Thus, the wear caused by long-term contact between the acetabular prosthesis and femoral head needs to be further verified by animal experiments.
In this study, the use of different methods for the internal fixation of posterior acetabular wall fractures was analysed by finite element analysis. According to the results, a PACP is recommended for simple posterior acetabular fractures. In the case of a comminuted posterior acetabular wall fracture, either an NPAWP or BPAWP can be used, depending on the patient’s condition. The TPAWP in this study is not recommended due to the areas of high stress. Regarding the method of acetabular prosthesis design, we propose the concept of “Break up to Make up” as a guide.
Availability of data and materials
The datasets used and/or analysed during the current study are available from the corresponding author on reasonable request.
Common reconstruction plate
Posterior acetabular column plate
Nonflanked posterior acetabular wall prosthesis
Biflanked posterior acetabular wall prosthesis
Triflanked posterior acetabular wall prosthesis
Chen MJ, Hollyer I, Wadhwa H, et al. Management of the posterior wall fracture in associated both column fractures of the acetabulum. Eur J Orthop Surg Traumatol. 2021;31:1047–54. https://doi.org/10.1007/s00590-020-02850-9.
Kim HT, Ahn J-M, Hur J-O, Lee J-S, Cheon S-J. Reconstruction of acetabular posterior wall fractures. Clin Orthop Surg. 2011;3:114–20. https://doi.org/10.4055/cios.2011.3.2.114.
Moed BR. Pearls: how to reduce and fix comminuted posterior acetabular wall fractures. Clin Orthop Relat Res. 2017;475:39–42. https://doi.org/10.1007/s11999-016-5104-9.
Ebraheim NA, Patil V, Liu J, Sanford CG Jr, Haman SP. Reconstruction of comminuted posterior wall fractures using the buttress technique: a review of 32 fractures. Int Orthop (SICOT). 2007;31:671–5. https://doi.org/10.1007/s00264-006-0246-0.
Perdue PW Jr, Tainter D, Toney C, Lee C. Evaluation and management of posterior wall acetabulum fractures. J Am Acad Orthop Surg. 2021;29(21):e1057–67. https://doi.org/10.5435/JAAOS-D-20-01301.
Lian S, Yang Z, Hu Z, et al. A combination of the K-L and S-P approaches for treating acetabular posterior wall factures accompanied by femoral head fractures with open reduction and internal fixation. BMC Surg. 2022;22:165. https://doi.org/10.1186/s12893-022-01597-w.
Sharma V, Bansal H, Mittal S, et al. Tricortical iliac crest graft as a salvageable option in the reconstruction of comminuted posterior wall acetabular fractures: our experience from a level 1 trauma centre. Arch Orthop Trauma Surg. 2021. https://doi.org/10.1007/s00402-021-04064-y.
Huang G, Chen K, Wang Y, Guo X. Finite element analysis of the Union Plate in treating elderly acetabular fracture patients. J Orthop Surg Res. 2022;17(1):56. https://doi.org/10.1186/s13018-022-02951-7.
Altun G, Saka G, Demir T, Elibol FKE, Polat MO. Precontoured buttress plate vs reconstruction plate for acetabulum posterior wall fractures: a biomechanical study. World J Orthop. 2019;10(5):219–27. https://doi.org/10.5312/wjo.v10.i5.219.
Majumder S, Roychowdhury A, Pal S. Dynamic response of the pelvis under side impact load—a three-dimensional finite element approach. Int J Crashworthiness. 2004;9:89–103. https://doi.org/10.1533/ijcr.2004.0274.
Ghouse S, Babu S, Nai K, Hooper PA, Jeffers JRT. The influence of laser parameters, scanning strategies and material on the fatigue strength of a stochastic porous structure. Addit Manuf. 2018;22:290–301. https://doi.org/10.1016/j.addma.2018.05.024.
Abbasnia A, Jafari M, Rohani A. A novel method for estimation of stress concentration factor of central cutouts located in orthotropic plate. J Braz Soc Mech Sci Eng. 2021;43:348. https://doi.org/10.1007/s40430-021-03061-x.
Fang H, Iqbal N, Van Staen G, et al. Experimental and numerical investigation of stress concentration at rib-to-crossbeam joint. Int J Steel Struct. 2021;21:360–80. https://doi.org/10.1007/s13296-020-00443-0.
Mikhailova NV, Petrov YV. Effect of impact time parameters on the dynamic strength in spall fracture. Phys Mesomech. 2021;24:9–13. https://doi.org/10.1134/S1029959921010021.
Abar B, Alonso-Calleja A, Kelly A, Kelly C, Gall K, West JL. 3D printing of high-strength, porous, elastomeric structures to promote tissue integration of implants. J Biomed Mater Res A. 2021;109(1):54–63. https://doi.org/10.1002/jbm.a.37006.
Baptista R, Guedes M. Morphological and mechanical characterization of 3D printed PLA scaffolds with controlled porosity for trabecular bone tissue replacement. Mater Sci Eng C Mater Biol Appl. 2021;118:111528. https://doi.org/10.1016/j.msec.2020.111528.
Yang N, Venezuela J, Almathami S, Dargusch M. Zinc-nutrient element based alloys for absorbable wound closure devices fabrication: current status, challenges, and future prospects. Biomaterials. 2022;280:121301. https://doi.org/10.1016/j.biomaterials.2021.121301.
Bi Y, He S, Du Y, et al. Effects of the configuration of a baffle–avalanche wall system on rock avalanches in Tibet Zhangmu: discrete element analysis. Bull Eng Geol Environ. 2019;78:2267–82. https://doi.org/10.1007/s10064-018-1284-8.
Li H, Qu X, Mao Y, et al. Custom acetabular cages offer stable fixation and improved hip scores for revision THA With severe bone defects. Clin Orthop Relat Res. 2016;474:731–40. https://doi.org/10.1007/s11999-015-4587-0.
Lagoa ALC, Wedemeyer C, von Knoch M, et al. A strut graft substitute consisting of a metal core and a polymer surface. J Mater Sci: Mater Med. 2008;19:417–24. https://doi.org/10.1007/s10856-006-0022-0.
Wu X. A biomechanical comparison of different fixation techniques for fractures of the acetabular posterior wall. Int Orthop (SICOT). 2018;42:673–9. https://doi.org/10.1007/s00264-017-3728-3.
Iselin LD, Wahl P, Studer P, et al. Associated lesions in posterior wall acetabular fractures: not a valid predictor of failure. J Orthopaed Traumatol. 2013;14:179–84. https://doi.org/10.1007/s10195-013-0247-x.
The National Natural Science Foundation of China (No.81873999, No.82072446).
Ethics approval and consent to participate
This research had been ethically approved by Huazhong University of Science and Technology Tongji Medical College Medical Ethics Committee (NO. S1060). Informed consent was obtained from all individual participants included in the study. All methods were performed in accordance with the Declarations of Helsinki.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
Huang, G., Wan, Y., Chen, K. et al. Finite element analysis of posterior acetabular column plate and posterior acetabular wall prostheses in treating posterior acetabular fractures. J Orthop Surg Res 18, 94 (2023). https://doi.org/10.1186/s13018-023-03535-9
- Acetabular posterior wall fracture
- Acetabular posterior column plate
- Acetabular posterior wall prosthesis
- Finite element analysis