Trabecular bone microdamage and microstructural stresses under uniaxial compression

The balance between local remodeling and accumulation of trabecular bone microdamage is believed to play an important role in the maintenance of skeletal integrity. However, the local mechanical parameters associated with microdamage initiation are not well understood. Using histological damage labeling, micro-CT imaging, and image-based finite element analysis, regions of trabecular bone microdamage were detected and registered to estimated microstructural von Mises effective stresses and strains, maximum principal stresses and strains, and strain energy density (SED). Bovine tibial trabecular bone cores underwent a stepwise uniaxial compression routine in which specimens were micro-CT imaged following each compression step. The results indicate that the mode of trabecular failure observed by micro-CT imaging agreed well with the polarity and distribution of stresses within an individual trabecula. Analysis of on-axis subsections within specimens provided significant positive relationships between microdamage and each estimated tissue stress, strain and SED parameter. In a more localized analysis, individual microdamaged and undamaged trabeculae were extracted from specimens loaded within the elastic region and to the apparent yield point. As expected, damaged trabeculae in both groups possessed significantly higher local stresses and strains than undamaged trabeculae. The results also indicated that microdamage initiation occurred prior to apparent yield at local principal stresses in the range of 88-121 MPa for compression and 35-43 MPa for tension and local principal strains of 0.46-0.63% in compression and 0.18-0.24% in tension. These data provide an important step towards understanding factors contributing to microdamage initiation and establishing local failure criteria for normal and diseased trabecular bone.     

Molecular dynamics simulations showing 1-palmitoyl-2-oleoyl-phosphatidylcholine (POPC) membrane mechanoporation damage under different strain paths

Continuum finite element material models used for traumatic brain injury lack local injury parameters necessitating nanoscale mechanical injury mechanisms be incorporated. One such mechanism is membrane mechanoporation, which can occur during physical insults and can be devastating to cells, depending on the level of disruption. The current study investigates the strain state dependence of phospholipid bilayer mechanoporation and failure. Using molecular dynamics, a simplified membrane, consisting of 72 1-palmitoyl-2-oleoyl-phosphatidylcholine (POPC) phospholipids, was subjected to equibiaxial, 2:1 non-equibiaxial, 4:1 non-equibiaxial, strip biaxial, and uniaxial tensile deformations at a von Mises strain rate of 5.45 × 10⁸ s^?1, resulting in velocities in the range of 1 to 4.6 m·s^?1. A water bridge forming through both phospholipid bilayer leaflets was used to determine structural failure. The stress magnitude, failure strain, headgroup clustering, and damage responses were found to be strain state-dependent. The strain state order of detrimentality in descending order was equibiaxial, 2:1 non-equibiaxial, 4:1 non-equibiaxial, strip biaxial, and uniaxial. The phospholipid bilayer failed at von Mises strains of .46, .47, .53, .77, and 1.67 during these respective strain path simulations. Additionally, a Membrane Failure Limit Diagram (MFLD) was created using the pore nucleation, growth, and failure strains to demonstrate safe and unsafe membrane deformation regions. This MFLD allowed representative equations to be derived to predict membrane failure from in-plane strains. These results provide the basis to implement a more accurate mechano-physiological internal state variable continuum model that captures lower length scale damage and will aid in developing higher fidelity injury models.     

Subject-specific finite element models implementing a maximum principal strain criterion are able to estimate failure risk and fracture location on human femurs tested in vitro

No agreement on the choice of the failure criterion to adopt for the bone tissue can be found in the literature among the finite element studies aiming at predicting fracture risk of bones. The use of stress-based criteria seems to prevail on strain-based ones, while basic bone biomechanics suggest using strain parameters to describe failure. The aim of the present combined experimental-numerical study was to verify, using subject-specific finite element models able to accurately predict strains, if a strain-based failure criterion could identify the failure patterns of bones. Three cadaver femurs were CT-scanned and subsequently fractured in a clinically relevant single-stance loading scenario. Load-displacement curves and high-speed movies were acquired to define the failure load and the location of fracture onset, respectively. Subject-specific finite element models of the three femurs were built from CT data following a validated procedure. A maximum principal strain criterion was implemented in the finite element models, and two stress-based criteria selected for comparison. The failure loads measured were applied to the models, and the computed risks of fracture were compared to the results of the experimental tests. The proposed principal strain criterion managed to correctly identify the level of failure risk and the location of fracture onset in all the modelled specimens, while Von Mises or maximum principal stress criteria did not give significant information. A maximum principal strain criterion can thus be defined a suitable candidate for the in vivo risk factor assessment on long bones.     

Biomechanical failure properties and microstructural content of ruptured and unruptured abdominal aortic aneurysms

PurposeTo test the hypothesis that ruptured abdominal aortic aneurysms (AAA) are globally weaker than unruptured ones.MethodsFour ruptured and seven unruptured AAA specimens were harvested whole from fresh cadavers during autopsies performed over an 18-month period. Multiple regionally distributed longitudinally oriented rectangular strips were cut from each AAA specimen for a total of 77 specimen strips. Strips were subjected to uniaxial extension until failure. Sections from approximately the strongest and weakest specimen strips were studied histologically and histochemically. From the load-extension data, failure tension, failure stress and failure strain were calculated. Rupture site characteristics such as location, arc length of rupture and orientation of rupture were also documented.ResultsThe failure tension, a measure of the tissue mechanical caliber was remarkably similar between ruptured and unruptured AAA (group mean±standard deviation of within-subject means: 11.2±2.3 versus 11.6±3.6 N/cm; p=0.866 by mixed model ANOVA). In post-hoc analysis, there was little difference between the groups in other measures of tissue mechanical caliber as well such as failure stress (95±28 versus 98±23 N/cm²; p=0.870), failure strain (0.39±0.09 versus 0.36±0.09; p=0.705), wall thickness (1.7±0.4 versus 1.5±0.4 mm; p=0.470) , and % coverage of collagen within tissue cross section (49.6±12.9% versus 60.8±9.6%; p=0.133). In the four ruptured AAA, primary rupture sites were on the lateral quadrants (two on left; one on left-posterior; one on right). Remarkably, all rupture lines had a longitudinal orientation and ranged from 1 to 6 cm in length.ConclusionThe findings are not consistent with the hypothesis that ruptured aortic aneurysms are globally weaker than unruptured ones.     

Computational analysis of biomechanical contributors to possible endovascular graft failure

This paper evaluates numerically coupled blood flow and wall structure interactions in a representative stented abdominal aortic aneurysm (AAA) model, leading potentially to endovascular graft (EVG) failure. A total of 12 biomechanical contributors to possible EVG migration were considered. The results show that after EVG insertion for the given model, the peak AAA sac-pressure was reduced to 14.2 mmHg (11.8% of p_lumen), and hence the maximum von Mises wall stress and wall deformation dropped by factors of 20 and 10, respectively. Thus, an EVG can significantly reduce sac pressure, mechanical stress, pulsatile wall motion, and the maximum diameter in AAAs and hence prevent AAA rupture effectively. In the absence of endoleaks, elevated sac-pressure can still be caused by fluid-structure interactions between the EVG, stagnant blood, and AAA wall. EVG migration forces vary from 1.4 to 7 N for different EVG geometries, material properties, and hemodynamic conditions. AAA-neck angle, iliac bifurcation angle, neck aorta-to-iliac diameter ratio, EVG size, aorto-uni-iliac EVG, and hypertension play important roles in generating forces potentially leading to EVG migration.     

Regional distribution of wall thickness and failure properties of human abdominal aortic aneurysm

The regional distribution of wall thickness and failure properties in human abdominal aortic aneurysm (AAA) was explored. Three unruptured and one ruptured AAA were harvested as a whole during necropsy. Thickness was measured at about every 1.5 cm² wall surface area for an average of 100 measurement sites per AAA. Multiple longitudinally oriented rectangular specimen strips were cut at various locations from each AAA for a total of 48 strips. The strips were subjected to uniaxial extension until failure. Wall thickness varied regionally and between AAA from as low as 0.23 mm at a rupture site to 4.26 mm at a calcified site (median=1.48 mm). Wall thickness was slightly lower in the posterior and right regions. The failure tension (ultimate) of specimen strips varied regionally and between AAA from 5.5 N/cm close to a blister site in the ruptured AAA to 42.3 N/cm at the undilated neck of a 4 cm diameter unruptured AAA (median=14.8 N/cm). Failure stress (ultimate) varied from 33.6 to 235.1 N/cm² (median=126.6 N/cm²). There was no perceptible pattern in failure properties along the circumference. Failure tension of specimen strips at or close to blisters was mostly low. The rupture site in the ruptured aneurysm had the lowest recorded wall thickness of 0.23 mm with only slightly higher readings within a 1 cm radius. The failure tension of the specimen strip close to the rupture site was low (11.1 N/cm) compared to its neighborhood in the ruptured aneurysm.     

Measuring and modeling patient-specific distributions of material properties in abdominal aortic aneurysm wall

Both the clinically established diameter criterion and novel approaches of computational finite element (FE) analyses for rupture risk stratification of abdominal aortic aneurysms (AAA) are based on assumptions of population-averaged, uniform material properties for the AAA wall. The presence of inter-patient and intra-patient variations in material properties is known, but has so far not been addressed sufficiently. In order to enable the preoperative estimation of patient-specific AAA wall properties in the future, we investigated the relationship between non-invasively assessable clinical parameters and experimentally measured AAA wall properties. We harvested n = 163 AAA wall specimens (n = 50 patients) during open surgery and recorded the exact excision sites. Specimens were tested for their thickness, elastic properties, and failure loads using uniaxial tensile tests. In addition, 43 non-invasively assessable patient-specific or specimen-specific parameters were obtained from recordings made during surgery and patient charts. Experimental results were correlated with the non-invasively assessable parameters and simple regression models were created to mathematically describe the relationships. Wall thickness was most significantly correlated with the metabolic activity at the excision site assessed by PET/CT (ρ = 0.499, P = 4 × 10^?7) and to thrombocyte counts from laboratory blood analyses (ρ = 0.445, P = 3 × 10^?9). Wall thickness was increased in patients suffering from diabetes mellitus, while it was significantly thinner in patients suffering from chronic kidney disease (CKD). Elastic AAA wall properties had significant correlations with the metabolic activity at the excision site (PET/CT), with existent calcifications, and with the diameter of the non-dilated aorta proximal to the AAA. Failure properties (wall strength and failure tension) had correlations with the patient’s medical history and with results from laboratory blood analyses. Interestingly, AAA wall failure tension was significantly reduced for patients with CKD and elevated blood levels of potassium and urea, respectively, both of which are associated with kidney disease. This study is a first step to a future preoperative estimation of AAA wall properties. Results can be conveyed to both the diameter criterion and FE analyses to refine rupture risk prediction. The fact that AAA wall from patients suffering from CKD featured reduced failure tension implies an increased AAA rupture risk for this patient group at comparably smaller AAA diameters.     

A methodology for developing anisotropic AAA phantoms via additive manufacturing

An Abdominal Aortic Aneurysm (AAA) is a permanent focal dilatation of the abdominal aorta at least 1.5 times its normal diameter. The criterion of maximum diameter is still used in clinical practice, although numerical studies have demonstrated the importance of biomechanical factors for rupture risk assessment. AAA phantoms could be used for experimental validation of the numerical studies and for pre-intervention testing of endovascular grafts. We have applied multi-material 3D printing technology to manufacture idealized AAA phantoms with anisotropic mechanical behavior. Different composites were fabricated and the phantom specimens were characterized by biaxial tensile tests while using a constitutive model to fit the experimental data. One composite was chosen to manufacture the phantom based on having the same mechanical properties as those reported in the literature for human AAA tissue; the strain energy and anisotropic index were compared to make this choice. The materials for the matrix and fibers of the selected composite are, respectively, the digital materials FLX9940 and FLX9960 developed by Stratasys. The fiber proportion for the composite is equal to 0.15. The differences between the composite behavior and the AAA tissue are small, with a small difference in the strain energy (0.4%) and a maximum difference of 12.4% in the peak Green strain ratio. This work represents a step forward in the application of 3D printing technology for the manufacturing of AAA phantoms with anisotropic mechanical behavior.     

A model of growth and rupture of abdominal aortic aneurysm

We present here a coupled mathematical model of growth and failure of the abdominal aortic aneurysm (AAA). The failure portion of the model is based on the constitutive theory of softening hyperelasticity where the classical hyperelastic law is enhanced with a new constant indicating the maximum energy that an infinitesimal material volume can accumulate without failure. The new constant controls material failure and it can be interpreted as the average energy of molecular bonds from the microstructural standpoint. The constitutive model is compared to the data from uniaxial tension tests providing an excellent fit to the experiment. The AAA failure model is coupled with a phenomenological theory of soft tissue growth. The unified theory includes both momentum and mass balance laws coupled with the help of the constitutive equations. The microstructural alterations in the production of elastin and remodeling of collagen are reflected in the changing macroscopic parameters characterizing tissue stiffness, strength and density. The coupled theory is used to simulate growth and rupture of an idealized spherical AAA. The results of the simulation showing possible AAA ruptures in growth are reasonable qualitatively while the quantitative calibration of the model will require further clinical observations and in vitro tests. The presented model is the first where growth and rupture are coupled.     

Mechanical stresses in abdominal aortic aneurysms: influence of diameter, asymmetry, and material anisotropy

Biomechanical studies suggest that one determinant of abdominal aortic aneurysm (AAA) rupture is related to the stress in the wall. In this regard, a reliable and accurate stress analysis of an in vivo AAA requires a suitable 3D constitutive model. To date, stress analysis conducted on AAA is mainly driven by isotropic tissue models. However, recent biaxial tensile tests performed on AAA tissue samples demonstrate the anisotropic nature of this tissue. The purpose of this work is to study the influence of geometry and material anisotropy on the magnitude and distribution of the peak wall stress in AAAs. Three-dimensional computer models of symmetric and asymmetric AAAs were generated in which the maximum diameter and length of the aneurysm were individually controlled. A five parameter exponential type structural strain-energy function was used to model the anisotropic behavior of the AAA tissue. The anisotropy is determined by the orientation of the collagen fibers (one parameter of the model). The results suggest that shorter aneurysms are more critical when asymmetries are present. They show a strong influence of the material anisotropy on the magnitude and distribution of the peak stress. Results confirm that the relative aneurysm length and the degree of aneurysmal asymmetry should be considered in a rupture risk decision criterion for AAAs.     

Strains and stresses in sub-dermal tissues of the buttocks are greater in paraplegics than in healthy during sitting

A pressure-related deep tissue injury (DTI) is a severe pressure ulcer, which initiates in muscle tissue overlying a bony prominence (e.g. the ischial tuberosities, IT) and progresses outwards through fat and skin, unnoticed by the paralyzed patient. We recently showed that internal strains and stresses in muscle and fat of individuals at anatomical sites susceptible to DTI can be evaluated by integrating Open-MRI scans with subject-specific finite element (FE) analyzes (Linder-Ganz et al., Journal of Biomechanics, 2007); however, sub-dermal soft tissue strains/stresses from paraplegics are still missing in literature. We hypothesize that the pathoanatomy of the buttocks in paraplegia increases the internal soft tissue loads under the IT, making these patients inherently susceptible to DTI. We hence compared the strain and stress peaks in the gluteus muscle and fat tissues under the IT of six healthy and six paraplegic patients, using the coupled MRI-FE method. Peak principal compression, principal tension, von Mises and shear strains in the gluteus were 1.2-, 3.1-, 1.4- and 1.4-fold higher in paraplegics than in healthy, respectively (p<0.02). Likewise, peak principal compression, principal tension, von Mises and shear stresses in the gluteus were 1.9-, 2.5-, 2.1- and 1.7-fold higher for the paraplegics (p<0.05). Peak gluteal compression and shear stresses decreased by as much as 70% when the paraplegic patients moved from a sitting to a lying posture, indicating on the effectiveness of recommending such patients to lie down after prolonged periods of sitting. This is the first attempt to compare internal soft tissue loads between paraplegic and healthy subjects, using an objective standardized bioengineering method of analysis. The findings support our hypothesis that internal tissue loads are significantly higher in paraplegics, and that postural changes significantly affect these loads. The method of analysis is useful for quantifying the effectiveness of various interventions to alleviate sub-dermal tissue loads at sites susceptible to pressure ulcers and DTI, including cushions, mattresses, recommendations for posture and postural changes, etc.     

The effects of aneurysm on the biaxial mechanical behavior of human abdominal aorta

The biomechanical response of normal and pathologic human abdominal aortic tissue to uniaxial loading conditions is insufficient for the characterization of its three-dimensional (3D) mechanical behavior. Planar biaxial mechanical evaluation allows for 3D constitutive modeling of nearly incompressible tissues, as well as the investigation of the nature of mechanical anisotropy. In the current study, 26 abdominal aortic aneurysm (AAA) tissue samples and 8 age-matched (> 60 years of age) nonaneurysmal abdominal aortic (AA) tissue samples were obtained and tested using a tension-controlled biaxial testing protocol. Graphical response functions (Sun et al., 2003. J. Biomech. Eng. 125, 372-380) were used as a guide to describe the pseudo-elastic response of AA and AAA. Based on the observed pseudo-elastic response, a four-parameter exponential strain energy function developed by Vito (1990. J. Biomech. Eng. 112, 153-159) was used from which both an individual specimen and group material parameter sets were determined for both AA and AAA. Peak Green strain values in the circumferential (Ethetatheta,max) and longitudinal (ELL,max) directions under an equibiaxial tension of 120 N/m were also compared. The strain energy function fit all of the individual specimens well with an average R2 of 0.95 +/- 0.02 and 0.90 +/- 0.02 (mean +/- SEM) for the AA and AAA groups, respectively. The average Ethetatheta,max at 200 N/m equibiaxial tension was found to be significantly smaller for AAAs as compared to AAs (0.07 +/- 0.01 versus 0.13 +/- 0.03, respectively; p < 0.01). There was also a pronounced increase in the circumferential stiffness for AAA tissue as compared to AA tissue, indicating a larger degree of anisotropy for this tissue as compared to age-matched AA tissue. We also observed that the four-parameter Fung-elastic model was not able to fit the AAA tissue mechanical response using physically realistic material parameter values. It was concluded that aneurysmal degeneration of the abdominal aorta is associated with an increase in mechanical anisotropy, with preferential stiffening in the circumferential direction.     

Age-related changes in human trabecular bone: Relationship between microstructural stress and strain and damage morphology

Accumulation of microdamage in aging and disease can cause skeletal fragility and is one of several factors contributing to osteoporotic fractures. To better understand the role of microdamage in fragility fracture, the mechanisms of bone failure must be elucidated on a tissue-level scale where interactions between bone matrix properties, the local biomechanical environment, and bone architecture are concurrently examined for their contributions to microdamage formation. A technique combining histological damage assessment of individual trabeculae with linear finite element solutions of trabecular von Mises and principal stress and strain was used to compare the damage initiation threshold between pre-menopausal (32-37 years, n=3 donors) and post-menopausal (71-80 years, n=3 donors) femoral cadaveric bone. Strong associations between damage morphology and stress and strain parameters were observed in both groups, and an age-related decrease in undamaged trabecular von Mises stress was detected. In trabeculae from younger donors, the 95% CI for von Mises stress on undamaged regions ranged from 50.7-67.9MPa, whereas in trabeculae from older donors, stresses were significantly lower (38.7-50.2, p<0.01). Local microarchitectural analysis indicated that thinner, rod-like trabeculae oriented along the loading axis are more susceptible to severe microdamage formation in older individuals, while only rod-like architecture was associated with severe damage in younger individuals. This study therefore provides insight into how damage initiation and morphology relate to local trabecular microstructure and the associated stresses and strains under loading. Furthermore, by comparison of samples from pre- and post-menopausal women, the results suggest that trabeculae from younger individuals can sustain higher stresses prior to microdamage initiation.     

Femoral head apparent density distribution predicted from bone stresses

A new theory relating bone morphology to applied stress is used to predict the apparent density distribution in the femoral head and neck. Cancellous bone is modeled as a self-optimizing material and cortical bone as a saturated (maximum possible bone density) response to stress in the bone tissue. Three different approaches are implemented relating bone apparent density to: (1) the von Mises stress, (2) the strain energy density in the mineralized tissue and (3) a defined closed effective stress (spherical stress). An iterative nonlinear three-dimensional finite element model is used to predict the apparent density distribution in the femoral head and neck for each of the three approaches. It is shown that the von Mises stress (an open effective stress) cannot accurately predict bone apparent density. It is shown that strain energy density and the defined closed effective stress can predict apparent density and that they give predictions consistent with the observed density pattern in the femoral head and neck.     

In vitro analysis of localized aneurysm rupture

In this study, bulge inflation tests were used to characterize the failure response of 15 layers of human ascending thoracic aortic aneurysms (ATAA). Full field displacement data were collected during each of the mechanical tests using a digital image stereo-correlation (DIS-C) system. Using the collected displacement data, the local stress fields at burst were derived and the thickness evolution was estimated during the inflation tests. It was shown that rupture of the ATAA does not systematically occur at the location of maximum stress, but in a weakened zone of the tissue where the measured fields show strain localization and localized thinning of the wall. Our results are the first to show the existence of weakened zones in the aneurysmal tissue when rupture is imminent. An understanding these local rupture mechanics is necessary to improve clinical assessments of aneurysm rupture risk. Further studies must be performed to determine if these weakened zones can be detected in vivo using non-invasive techniques.     

A simple method of estimating the stress acting on a bilaterally symmetric abdominal aortic aneurysm

We extended a method of estimating the stress acting on an axisymmetric abdominal aortic aneurysm (AAA) under a load in vivo (Elger, D. F., Blackketter, D. M., Budwig, R. S., Johansen, K. H. (1996) The influence of shape on the stresses in model abdominal aortic aneurysms, Journal of Biomechanical Engineering, 118, pp. 326-32.) to bilaterally-symmetric AAAs, which are symmetric about the sagittal plane. Stresses were calculated along the anterior and posterior median lines of the AAA wall. Of the two force equilibrium equations, the Laplace equation held in this study. The longitudinal equilibrium was extended to hold by approximating the meridional tension and the directional cosine of the wall surface as constants along the circumference. The estimated stresses were compared with the results of a finite element analysis. Comparisons showed that the maximal principal stress, usually the circumferential stress or sometimes the meridional stress depending on location, sufficiently represented the wall stress. The proposed method provides a reasonable index for evaluating the rupture risk using the peak value of the maximal principal stress and its location without using the stress-free geometry and constitutive equation.     

A comparison of uniaxial and biaxial mechanical properties of the annulus fibrosus: a porcine model

The annulus fibrosus of the intervertebral disk experiences multidirectional tension in vivo, yet the majority of mechanical property testing has been uniaxial. Therefore, our understanding of how this complex multilayered tissue responds to loading may be deficient. This study aimed to determine the mechanical properties of porcine annular samples under uniaxial and biaxial tensile loading. Two-layer annulus samples were isolated from porcine disks from four locations: anterior superficial, anterior deep, posterior superficial, and posterior deep. These tissues were then subjected to three deformation conditions each to a maximal stretch ratio of 1.23: uniaxial, constrained uniaxial, and biaxial. Uniaxial deformation was applied in the circumferential direction, while biaxial deformation was applied simultaneously in the circumferential and compressive directions. Constrained uniaxial consisted of a stretch ratio of 1.23 in the circumferential direction while holding the tissue stationary in the axial direction. The maximal stress and stress-stretch ratio (S-S) moduli determined from the biaxial tests were significantly higher than those observed during both the uniaxial tests (maximal stress, 97.1% higher during biaxial; p=0.002; S-S moduli, 117.9% higher during biaxial; p=0.0004) and the constrained uniaxial tests (maximal stress, 46.8% higher during biaxial; S-S moduli, 82.9% higher during biaxial). These findings suggest that the annulus is subjected to higher stresses in vivo when under multidirectional tension.     

Influence of bicortical techniques in internal connection placed in premaxillary area by 3D finite element analysis

The aim of study was to evaluate the stress distribution in implant-supported prostheses and peri-implant bone using internal hexagon (IH) implants in the premaxillary area, varying surgical techniques (conventional, bicortical and bicortical in association with nasal floor elevation), and loading directions (0°, 30° and 60°) by three-dimensional (3D) finite element analysis. Three models were designed with Invesalius, Rhinoceros 3D and Solidworks software. Each model contained a bone block of the premaxillary area including an implant (IH, Ø4 × 10 mm) supporting a metal-ceramic crown. 178 N was applied in different inclinations (0°, 30°, 60°). The results were analyzed by von Mises, maximum principal stress, microstrain and displacement maps including ANOVA statistical test for some situations. Von Mises maps of implant, screws and abutment showed increase of stress concentration as increased loading inclination. Bicortical techniques showed reduction in implant apical area and in the head of fixation screws. Bicortical techniques showed slight increase stress in cortical bone in the maximum principal stress and microstrain maps under 60° loading. No differences in bone tissue regarding surgical techniques were observed. As conclusion, non-axial loads increased stress concentration in all maps. Bicortical techniques showed lower stress for implant and screw; however, there was slightly higher stress on cortical bone only under loads of higher inclinations (60°).