Keratoconus prediction using a finite element model of the cornea with local biomechanical properties

PURPOSE: The ability to predict and understand which biomechanical properties of the cornea are responsible for the stability or progression of keratoconus may be an important clinical and surgical tool for the eye-care professional. We have developed a finite element model of the cornea, that tries to predicts keratoconus-like behavior and its evolution based on material properties of the corneal tissue. METHODS: Corneal material properties were modeled using bibliographic data and corneal topography was based on literature values from a schematic eye model. Commercial software was used to simulate mechanical and surface properties when the cornea was subject to different local parameters, such as elasticity. RESULTS: The simulation has shown that, depending on the corneal initial surface shape, changes in local material properties and also different intraocular pressures values induce a localized protuberance and increase in curvature when compared to the remaining portion of the cornea. CONCLUSIONS: This technique provides a quantitative and accurate approach to the problem of understanding the biomechanical nature of keratoconus. The implemented model has shown that changes in local material properties of the cornea and intraocular pressure are intrinsically related to keratoconus pathology and its shape/curvature.

Stress distribution in the cervical region of an upper central incisor in a 3D finite element model

The aim of this study was to evaluate the stress distribution in the cervical region of a sound upper central incisor in two clinical situations, standard and maximum masticatory forces, by means of a 3D model with the highest possible level of fidelity to the anatomic dimensions. Two models with 331,887 linear tetrahedral elements that represent a sound upper central incisor with periodontal ligament, cortical and trabecular bones were loaded at 45º in relation to the tooth's long axis. All structures were considered to be homogeneous and isotropic, with the exception of the enamel (anisotropic). A standard masticatory force (100 N) was simulated on one of the models, while on the other one a maximum masticatory force was simulated (235.9 N). The software used were: PATRAN for pre- and post-processing and Nastran for processing. In the cementoenamel junction area, tensile forces reached 14.7 MPa in the 100 N model, and 40.2 MPa in the 235.9 N model, exceeding the enamel's tensile strength (16.7 MPa). The fact that the stress concentration in the amelodentinal junction exceeded the enamel's tensile strength under simulated conditions of maximum masticatory force suggests the possibility of the occurrence of non-carious cervical lesions such as abfractions.

A Study on Unconditionally Stable FDTD Methods for the Modeling of Metamaterials

We assess the performance of three unconditionally stable finite-difference time-domain (FDTD) methods for the modeling of doubly dispersive metamaterials: 1) locally one-dimensional FDTD; 2) locally one-dimensional FDTD with Strang splitting; and (3) alternating direction implicit FDTD. We use both double-negative media and zero-index media as benchmarks.

A discussion on volume change in the plastic phase

This technical note discusses the possibility of using a more simplified scheme to estimate the plastic multiplier when some material shows volume changes, e.g. soil, balsa wood foam and other similar materials. Two procedures regarding volume changes during the plastic phase are discussed here. The first one is the classic procedure applied to non-associative plasticity, for which a Drucker-Prager-like surface is adopted to represent the plastic potential. For the second procedure, the plastic potential is not explicitly known, however, its orthogonal direction is chosen respecting a plastic volume change parameter similar to Poisson`s ratio. Copyright (C) 2007 John Wiley & Sons, Ltd.

GENERALIZED FINITE ELEMENT METHOD FOR NONLINEAR THREE-DIMENSIONAL ANALYSIS OF SOLIDS

The Generalized Finite Element Method (GFEM) is employed in this paper for the numerical analysis of three-dimensional solids tinder nonlinear behavior. A brief summary of the GFEM as well as a description of the formulation of the hexahedral element based oil the proposed enrichment strategy are initially presented. Next, in order to introduce the nonlinear analysis of solids, two constitutive models are briefly reviewed: Lemaitre`s model, in which damage and plasticity are coupled, and Mazars`s damage model suitable for concrete tinder increased loading. Both models are employed in the framework of a nonlocal approach to ensure solution objectivity. In the numerical analyses carried out, a selective enrichment of approximation at regions of concern in the domain (mainly those with high strain and damage gradients) is exploited. Such a possibility makes the three-dimensional analysis less expensive and practicable since re-meshing resources, characteristic of h-adaptivity, can be minimized. Moreover, a combination of three-dimensional analysis and the selective enrichment presents a valuable good tool for a better description of both damage and plastic strain scatterings.

On the finite element modeling of fatigue crack growth in pressurized cylindrical shells

A methodology for the computational modeling of the fatigue crack growth in pressurized shell structures, based on the finite element method and concepts of Linear Elastic Fracture Mechanics, is presented. This methodology is based on that developed by Potyondy [Potyondy D, Wawrzynek PA, Ingraffea, AR. Discrete crack growth analysis methodology for through crack in pressurized fuselage structures. Int J Numer Methods Eng 1995;38:1633-1644], which consists of using four stress intensity factors, computed from the modified crack integral method, to predict the fatigue propagation life as well as the crack trajectory, which is computed as part of the numerical simulation. Some issues not presented in the study of Potyondy are investigated herein such as the influence of the crack increment size and the number of nodes per element (4 or 9 nodes) on the simulation results by means of a fatigue crack propagation simulation of a Boeing 737 airplane fuselage. The results of this simulation are compared with experimental results and those obtained by Potyondy [1]. (C) 2008 Elsevier Ltd. All rights reserved.

Galerkin finite element method for incompressible thermofluid flows framed within the bond graph theory

In this paper a bond graph methodology is used to model incompressible fluid flows with viscous and thermal effects. The distinctive characteristic of these flows is the role of pressure, which does not behave as a state variable but as a function that must act in such a way that the resulting velocity field has divergence zero. Velocity and entropy per unit volume are used as independent variables for a single-phase, single-component flow. Time-dependent nodal values and interpolation functions are introduced to represent the flow field, from which nodal vectors of velocity and entropy are defined as state variables. The system for momentum and continuity equations is coincident with the one obtained by using the Galerkin method for the weak formulation of the problem in finite elements. The integral incompressibility constraint is derived based on the integral conservation of mechanical energy. The weak formulation for thermal energy equation is modeled with true bond graph elements in terms of nodal vectors of temperature and entropy rates, resulting a Petrov-Galerkin method. The resulting bond graph shows the coupling between mechanical and thermal energy domains through the viscous dissipation term. All kind of boundary conditions are handled consistently and can be represented as generalized effort or flow sources. A procedure for causality assignment is derived for the resulting graph, satisfying the Second principle of Thermodynamics. (C) 2007 Elsevier B.V. All rights reserved.

Coarse-grid higher order finite-difference time-domain algorithm with low dispersion errors

Higher order (2,4) FDTD schemes used for numerical solutions of Maxwell`s equations are focused on diminishing the truncation errors caused by the Taylor series expansion of the spatial derivatives. These schemes use a larger computational stencil, which generally makes use of the two constant coefficients, C-1 and C-2, for the four-point central-difference operators. In this paper we propose a novel way to diminish these truncation errors, in order to obtain more accurate numerical solutions of Maxwell`s equations. For such purpose, we present a method to individually optimize the pair of coefficients, C-1 and C-2, based on any desired grid size resolution and size of time step. Particularly, we are interested in using coarser grid discretizations to be able to simulate electrically large domains. The results of our optimization algorithm show a significant reduction in dispersion error and numerical anisotropy for all modeled grid size resolutions. Numerical simulations of free-space propagation verifies the very promising theoretical results. The model is also shown to perform well in more complex, realistic scenarios.

Three-Dimensional Stress Distribution in the Human Periodontal Ligament in Masticatory, Parafunctional, and Trauma Loads: Finite Element Analysis

Background: The presence of the periodontal ligament (PDL) makes it possible to absorb and distribute loads produced during masticatory function and other tooth contacts into the alveolar process via the alveolar bone proper. However, several factors affect the integrity of periodontal structures causing the destruction of the connective matrix and cells, the loss of fibrous attachment, and the resorption of alveolar bone. Methods: The purpose of this study was to evaluate the stress distribution by finite element analysis in a PDL in three-dimensional models of the upper central incisor under three different load conditions: 100 N occlusal loading at 45 degrees (model 1: masticatory load); 500 N at the incisal edge at 45 degrees (model 2: parafunctional habit); and 800 N at the buccal surface at 90 degrees (model 3: trauma case). The models were built from computed tomography scans. Results: The stress distribution was quite different among the models. The most significant values (harmful) of tensile and compressive stresses were observed in models 2 and 3, with similarly distinct patterns of stress distributions along the PDL. Tensile stresses were observed along the internal and external aspects of the PDL, mostly at the cervical and middle thirds. Conclusions: The stress generation in these models may affect the integrity of periodontal structures. A better understanding of the biomechanical behavior of the PDL under physiologic and traumatic loading conditions might enhance the understanding of the biologic reaction of the PDL in health and disease. J Periodontol 2009;80:1859-1867.

A comparison of methods to identify open-lung PEEP

Purpose: Many methods exist in the literature for identifying PEEP to set in ARDS patients following a lung recruitment maneuver (RM). We compared ten published parameters for setting PEEP following a RM. Methods: Lung injury was induced by bilateral lung lavage in 14 female Dorset sheep, yielding a PaO(2) 100-150 mmHg at F(I)O(2) 1.0 and PEEP 5 cmH(2)O. A quasi-static P-V curve was then performed using the supersyringe method; PEEP was set to 20 cmH(2)O and a RM performed with pressure control ventilation (inspiratory pressure set to 40-50 cmH(2)O), until PaO(2) + PaCO(2) > 400 mmHg. Following the RM, a decremental PEEP trial was performed. The PEEP was decreased in 1 cmH(2)O steps every 5 min until 15 cmH(2)O was reached. Parameters measured during the decremental PEEP trial were compared with parameters obtained from the P-V curve. Results: For setting PEEP, maximum dynamic tidal respiratory compliance, maximum PaO(2), maximum PaO(2) + PaCO(2), and minimum shunt calculated during the decremental PEEP trial, and the lower Pflex and point of maximal compliance increase on the inflation limb of the P-V curve (Pmci,i) were statistically indistinguishable. The PEEP value obtained using the deflation upper Pflex and the point of maximal compliance decrease on the deflation limb were significantly higher, and the true inflection point on the inflation limb and minimum PaCO(2) were significantly lower than the other variables. Conclusion: In this animal model of ARDS, dynamic tidal respiratory compliance, maximum PaO(2), maximum PaO(2) + PaCO(2), minimum shunt, inflation lower Pflex and Pmci,i yield similar values for PEEP following a recruitment maneuver.

Effectiveness of the CAGE questionnaire, gamma-glutamyltransferase and mean corpuscular volume of red blood cells as markers for alcohol-related problems in the workplace

Objective: To evaluate the usefulness of gamma-glutamyltransferase (GGT) and mean corpuscular volume (MCV), as well as that of the CAGE questionnaire, in workplace screening for alcohol abuse/dependence. Methods: A total of 183 male employees were submitted to structured interviews (Structured Clinical Interview for DSM-IV 2.0 and CAGE questionnaire). Blood samples were collected. Diagnostic accuracy and odds ratio were determined for the CAGE, GGT and MCV. Results: The CAGE questionnaire presented the best sensitivity for alcohol dependence (91%; specificity, 87.8%) and for alcohol abuse (87.5%, specificity, 80.9%), which increased when the questionnaire was used in combination with GGT (sensitivity, 100% and 87.5%, respectively; specificity, 68% and 61.5, respectively). CAGE positive results and/or alterations in GGT were less likely to occur among employees not presenting alcohol abuse/ dependence than among those presenting such abuse (OR for CAGE = 13, p < 0.05; OR for CAGE-GGT = 11, p < 0.05) or dependence (OR for CAGE = 76, p < 0.0 1; OR for GGT = 5, p < 0.0 1). Employees not presenting alcohol abuse/dependence were also several times more likely to present negative CAGE or GGT results. Conclusions: The use short, simple questionnaires, combined with that of low-cost biochemical markers, such as GGT, can serve as an initial screening for alcohol-related problems, especially for employees in hazardous occupations. The data provided can serve to corroborate clinical findings. (C) 2008 Elsevier Ltd. All rights reserved.

Urodynamic evaluation of patients with lower urinary tract symptoms and small prostate volume

Introduction: The relevance of prostate size in the pathophysiology of lower urinary tract symptoms (LUTS) is controversial. We evaluated the urodynamic findings in patients with LUTS and small prostate volumes. Materials and Methods: 84 patients aged >= 50 years with LUTS and prostates < 40 ml were evaluated. All had an International Prostate Symptom Score (IPSS) >= 8. Average age was 62.0 +/- 8.1 years. We evaluated the impact of bladder outlet obstruction (BOO) and detrusor overactivity (DO) on the voiding symptoms and urodynamic findings. Results: Mean prostate volume was 29.2 +/- 7.2 ml and mean IPSS was 13.5 +/- 4.6. BOO was the main finding, affecting 42 (50.0%) patients, followed by detrusor underactivity (DU) in 41 (48.8%) and DO in 28 (33.3%) patients. Patients without BOO were significantly older than the obstructed (64.0 +/- 8.8 and 60.1 +/- 6.9 years, respectively; p = 0.026) and had an increased prevalence of DU (76.2 and 21.4%, respectively; p < 0.001). Comparison of patients with and without DO showed reduced bladder capacity and compliance in the DO group (p < 0.001). No other comparisons were significant. Conclusion: Half of the patients with LUTS and small prostates are not obstructed and may have DO or decreased detrusor contractility as the basis for their voiding symptoms. Our results emphasize the value of urodynamics in this population, especially when invasive treatments are being considered. Copyright (c) 2008 S. Karger AG, Basel.

A Comparison Between Pulse Pressure Variation and Right End Diastolic Volume Index as Guides to Resuscitation in a Model of Hemorrhagic Shock in Pigs

Background: Different hemodynamic parameters including static indicators of cardiac preload as right ventricular end-diastolic volume index (RVEDVI) and dynamic parameters as pulse pressure variation (PPV) have been used in the decision-making process regarding volume expansion in critically ill patients. The objective of this study was to compare fluid resuscitation guided by either PPV or RVEDVI after experimentally induced hemorrhagic shock. Methods: Twenty-six anesthetized and mechanically ventilated pigs were allocated into control (group I), PPV (group II), or RVEDVI (group III) group. Hemorrhagic shock was induced by blood withdrawal to target mean arterial pressure of 40 mm Hg, maintained for 60 minutes. Parameters were measured at baseline, time of shock, 60 minutes after shock, immediately after resuscitation with hydroxyethyl starch 6% (130/0.4), 1 hour and 2 hours thereafter. The endpoint of fluid resuscitation was determined as the baseline values of PPV and RVEDVI. Statistical analysis of data was based on analysis of variance for repeated measures followed by the Bonferroni test (p < 0.05). Results: Volume and time to resuscitation were higher in group III than in group II (group III = 1,305 +/- 331 mL and group II = 965 +/- 245 mL, p < 0.05; and group III = 24.8 +/- 4.7 minutes and group II = 8.8 +/- 1.3 minutes, p < 0.05, respectively). All static and dynamic parameters and biomarkers of tissue oxygenation were affected by hemorrhagic shock and nearly all parameters were restored after resuscitation in both groups. Conclusion: In the proposed model of hemorrhagic shock, resuscitation to the established endpoints was achieved within a smaller amount of time and with less volume when guided by PPV than when guided by pulmonary artery catheter-derived RVEDVI.

Small volume resuscitation with 3% hypertonic saline solution decrease inflammatory response and attenuates end organ damage after controlled hemorrhagic shock

BACKGROUND: Recently, studies have been conducted examining the efficacy of 3% hypertonic saline solution (HS) for the treatment of traumatic brain injury; however, few studies have analyzed the effects of 3% hemorrhagic shock during hemorrhagic shock. The aim of this study was to test the potential immunomodulatory benefits of 3% hemorrhagic shock resuscitation over standard fluid resuscitation. METHODS: Wistar rats were bled to a mean arterial pressure of 35 mm Hg and then randomized into 3 groups: those treated with lactated Ringer`s solution (LR; 33 mL/kg, n = 7), 3% HS (10 mL/kg, n = 7), and 7.5% HS (4 mL/kg, n = 7). Half of the extracted blood was reinfused after fluid resuscitation. Animals that did not undergo shock served as controls (n = 5). Four hours after hemorrhagic shock, blood was collected for the evaluation of tumor necrosis factor-a and interleukin-6 by enzyme immunoassay. Lung and intestinal samples were obtained for histopathologic analysis. RESULTS: Animals in the HS groups had significantly higher mean arterial pressure than those in the LR group 1 hour after treatment. Osmolarity and sodium levels were markedly elevated in the HS groups. Tumor necrosis factor-alpha and interleukin-6 levels were similar between the control and HS groups but significantly higher in the LR group (P < .05). The lung injury score was significantly higher in the LR group compared with the 7.5% HS and 3% HS groups (5.7 +/- 0.7, 2.1 +/- 0.4, and 2.7 +/- 0.5, respectively). Intestinal injury was attenuated in the 7.5% HS and 3% HS groups compared with the LR group (2.0 +/- 0.6, 2.3 +/- 0.4, and 5.9 +/- 0.6, respectively). CONCLUSIONS: A small-volume resuscitation strategy modulates the inflammatory response and decreases end-organ damage after HS. Three percent HS provides immunomodulatory and metabolic effects similar to those observed with conventional concentrations of HS. (C) 2009 Elsevier Inc. All rights reserved.

Prognostic value of left atrial volume in patients who underwent dobutamine stress echocardiography for known or suspected coronary artery disease

Background Left atrial volume indexed (LAVI) has been reported as a predictor of cardiovascular events. We sought to determine the prognostic value of LAVI for predicting the outcome of patients who underwent dobutamine stress echocardiography (DSE) for known or suspected coronary artery disease (CAD). Methods From January 2000 to July 2005, we studied 981 patients who underwent DSE and off-line measurements of LAVI. The value of DSE over clinical and LAVI data was examined using a stepwise log-rank test. Results During a median follow-up of 24 months, 56 (6%) events occurred. By univariate analysis, predictors of events were male sex, diabetes mellitus, previous myocardial infarction, left ventricular ejection fraction (LVEF), left atrial diameter indexed, LAVI, and abnormal DSE. By multivariate analysis, independent predictors were LVEF (relative risk [RR] = 0.98, 95% CI 0.95-1.00), LAVI (RR = 1.04, 95% CI 1.02-1.05), and abnormal DSE (RR = 2.70, 95% CI 1.28-5.69). In an incremental multivariate model, LAVI was additional to clinical data for predicting events (chi(2) 36.8, P < .001). The addition of DSE to clinical and LAVI yielded incremental information (chi(2) 55.3, P < .001). The 3-year event-free survival in patients with normal DSE and LAVI <= 33 mL/m(2) was 96%; with abnormal DSE and LAVI <= 33 mL/m(2), 91%; with normal DSE and LAVI >34 mL/m(2), 83%; and with abnormal DSE and LAVI >34 mL/m(2) 51%. Conclusion Left atrial volume indexed provides independent prognostic information in patients who underwent DSE for known or suspected CAD. Among patients with normal DSE, those with larger LAVI had worse outcome, and among patients with abnormal DSE, LAVI was still predictive. (Am Heart J 2008; 156:1110-6.)