GENERALIZED FINITE ELEMENT METHOD ON NONCONVENTIONAL HYBRID-MIXED FORMULATION

The generalized finite element method (GFEM) is applied to a nonconventional hybrid-mixed stress formulation (HMSF) for plane analysis. In the HMSF, three approximation fields are involved: stresses and displacements in the domain and displacement fields on the static boundary. The GFEM-HMSF shape functions are then generated by the product of a partition of unity associated to each field and the polynomials enrichment functions. In principle, the enrichment can be conducted independently over each of the HMSF approximation fields. However, stability and convergence features of the resulting numerical method can be affected mainly by spurious modes generated when enrichment is arbitrarily applied to the displacement fields. With the aim to efficiently explore the enrichment possibilities, an extension to GFEM-HMSF of the conventional Zienkiewicz-Patch-Test is proposed as a necessary condition to ensure numerical stability. Finally, once the extended Patch-Test is satisfied, some numerical analyses focusing on the selective enrichment over distorted meshes formed by bilinear quadrilateral finite elements are presented, thus showing the performance of the GFEM-HMSF combination.

Standard surgical treatment for benign prostatic hyperplasia is safe for patients over 75 years: Analysis of 100 cases from a high-volume urologic center

OBJECTIVES: In this study, we aimed to determine the complications of standard surgical treatments among patients over 75 years in a high-volume urologic center. METHODS: We analyzed 100 consecutive patients older than 75 years who had undergone transurethral prostatic resection of the prostate or open prostatectomy for treatment of benign prostatic hyperplasia from January 2008 to March 2010. We analyzed patient age, prostate volume, prostate-specific antigen level, international prostatic symptom score, quality of life score, urinary retention, co-morbidities, surgical technique and satisfaction with treatment. RESULTS: Median age was 79 years. Forty-eight patients had undergone transurethral prostatic resection of the prostate, and 52 had undergone open prostatectomy. The median International Prostatic Symptom Score was 20, the median prostate volume was 83 g, 51% were using an indwelling bladder catheter, and the median prostate-specific antigen level was 5.0 ng/ml. The most common comorbidities were hypertension, diabetes and coronary disease. After a median follow-up period of 17 months, most patients were satisfied. Complications were present in 20% of cases. The most common urological complication was urethral stenosis, followed by bladder neck sclerosis, urinary fistula, late macroscopic hematuria and persistent urinary incontinence. The most common clinical complication was myocardial infarction, followed by acute renal failure requiring dialysis. Incidental carcinoma of the prostate was present in 6% of cases. One case had urothelial bladder cancer. CONCLUSIONS: Standard surgical treatments for benign prostatic hyperplasia are safe and satisfactory among the elderly. Complications are infrequent, and urethral stenosis is the most common. No clinical variable is associated with the occurrence of complications.

Numerical assessment of stability of interface discontinuous finite element pressure spaces

The stability of two recently developed pressure spaces has been assessed numerically: The space proposed by Ausas et al. [R.F. Ausas, F.S. Sousa, G.C. Buscaglia, An improved finite element space for discontinuous pressures, Comput. Methods Appl. Mech. Engrg. 199 (2010) 1019-1031], which is capable of representing discontinuous pressures, and the space proposed by Coppola-Owen and Codina [A.H. Coppola-Owen, R. Codina, Improving Eulerian two-phase flow finite element approximation with discontinuous gradient pressure shape functions, Int. J. Numer. Methods Fluids, 49 (2005) 1287-1304], which can represent discontinuities in pressure gradients. We assess the stability of these spaces by numerically computing the inf-sup constants of several meshes. The inf-sup constant results as the solution of a generalized eigenvalue problems. Both spaces are in this way confirmed to be stable in their original form. An application of the same numerical assessment tool to the stabilized equal-order P-1/P-1 formulation is then reported. An interesting finding is that the stabilization coefficient can be safely set to zero in an arbitrary band of elements without compromising the formulation's stability. An analogous result is also reported for the mini-element P-1(+)/P-1 when the velocity bubbles are removed in an arbitrary band of elements. (C) 2012 Elsevier B.V. All rights reserved.

A new enrichment space for the treatment of discontinuous pressures in multi-fluid flows

In this work, a new enrichment space to accommodate jumps in the pressure field at immersed interfaces in finite element formulations, is proposed. The new enrichment adds two degrees of freedom per element that can be eliminated by means of static condensation. The new space is tested and compared with the classical P1 space and to the space proposed by Ausas et al (Comp. Meth. Appl. Mech. Eng., Vol. 199, 10191031, 2010) in several problems involving jumps in the viscosity and/or the presence of singular forces at interfaces not conforming with the element edges. The combination of this enrichment space with another enrichment that accommodates discontinuities in the pressure gradient has also been explored, exhibiting excellent results in problems involving jumps in the density or the volume forces. Copyright (c) 2011 John Wiley & Sons, Ltd.

Effects of volume resuscitation on splanchnic perfusion in canine model of severe sepsis induced by live Escherichia coli infusion

Abstract Introduction We conducted the present study to investigate whether early large-volume crystalloid infusion can restore gut mucosal blood flow and mesenteric oxygen metabolism in severe sepsis. Methods Anesthetized and mechanically ventilated male mongrel dogs were challenged with intravenous injection of live Escherichia coli (6 × 109 colony-forming units/ml per kg over 15 min). After 90 min they were randomly assigned to one of two groups – control (no fluids; n = 13) or lactated Ringer's solution (32 ml/kg per hour; n = 14) – and followed for 60 min. Cardiac index, mesenteric blood flow, mean arterial pressure, systemic and mesenteric oxygen-derived variables, blood lactate and gastric carbon dioxide tension (PCO2; by gas tonometry) were assessed throughout the study. Results E. coli infusion significantly decreased arterial pressure, cardiac index, mesenteric blood flow, and systemic and mesenteric oxygen delivery, and increased arterial and portal lactate, intramucosal PCO2, PCO2 gap (the difference between gastric mucosal and arterial PCO2), and systemic and mesenteric oxygen extraction ratio in both groups. The Ringer's solution group had significantly higher cardiac index and systemic oxygen delivery, and lower oxygen extraction ratio and PCO2 gap at 165 min as compared with control animals. However, infusion of lactated Ringer's solution was unable to restore the PCO2 gap. There were no significant differences between groups in mesenteric oxygen delivery, oxygen extraction ratio, or portal lactate at the end of study. Conclusion Significant disturbances occur in the systemic and mesenteric beds during bacteremic severe sepsis. Although large-volume infusion of lactated Ringer's solution restored systemic hemodynamic parameters, it was unable to correct gut mucosal PCO2 gap.

Small volume of hypertonic saline as the initial fluid replacement in experimental hypodynamic sepsis

Abstract Introduction We conducted the present study to examine the effects of hypertonic saline solution (7.5%) on cardiovascular function and splanchnic perfusion in experimental sepsis. Methods Anesthetized and mechanically ventilated mongrel dogs received an intravenous infusion of live Escherichia coli over 30 minutes. After 30 minutes, they were randomized to receive lactated Ringer's solution 32 ml/kg (LR; n = 7) over 30 minutes or 7.5% hypertonic saline solution 4 ml/kg (HS; n = 8) over 5 minutes. They were observed without additional interventions for 120 minutes. Cardiac output (CO), mean arterial pressure (MAP), portal and renal blood flow (PBF and RBF, respectively), gastric partial pressure of CO2 (pCO2; gas tonometry), blood gases and lactate levels were assessed. Results E. coli infusion promoted significant reductions in CO, MAP, PBF and RBF (approximately 45%, 12%, 45% and 25%, respectively) accompanied by an increase in lactate levels and systemic and mesenteric oxygen extraction (sO2ER and mO2ER). Widening of venous-arterial (approximately 15 mmHg), portal-arterial (approximately 18 mmHg) and gastric mucosal-arterial (approximately 55 mmHg) pCO2 gradients were also observed. LR and HS infusion transiently improved systemic and regional blood flow. However, HS infusion was associated with a significant and sustained reduction of systemic (18 ± 2.6 versus 38 ± 5.9%) and mesenteric oxygen extraction (18.5 ± 1.9 versus 36.5 ± 5.4%), without worsening other perfusional markers. Conclusion A large volume of LR or a small volume of HS promoted similar transient hemodynamic benefits in this sepsis model. However, a single bolus of HS did promote sustained reduction of systemic and mesenteric oxygen extraction, suggesting that hypertonic saline solution could be used as a salutary intervention during fluid resuscitation in septic patients.

Morse taper implants at different bone levels: a finite element analysis of stress distribution

AIM: To explore the biomechanical effects of the different implantation bone levels of Morse taper implants, employing a finite element analysis (FEA). METHODS: Dental implants (TitamaxCM) with 4x13 mm and 4x11 mm, and their respective abutments with 3.5 mm height, simulating a screwed premolar metal-ceramic crown, had their design performed using the software AnsysWorkbench 10.0. They were positioned in bone blocks, covered by 2.5 mm thickness of mucosa. The cortical bone was designed with 1.5 mm thickness and the trabecular bone completed the bone block. Four groups were formed: group 11CBL (11 mm implant length on cortical bone level), group 11TBL (11 mm implant length on trabecular bone level), group 13CBL (13mm implant length on cortical bone level) and group 13TBL (13 mm implant length on trabecular bone level). Oblique 200 N loads were applied. Von Mises equivalent stresses in cortical and trabecular bones were evaluated with the same design program. RESULTS: The results were shown qualitatively and quantitatively by standard scales for each type of bone. By the results obtained, it can be suggested that positioning the implant completely in trabecular bone brings harm with respect to the generated stresses. Its implantation in the cortical bone has advantages with respect to better anchoring and locking, reflecting a better dissipation of the stresses along the implant/bone interfaces. In addition, the search for anchoring the implant in its apical region in cortical bone is of great value to improve stabilization and consequently better stress distribution. CONCLUSIONS: The implant position slightly below the bone in relation to the bone crest brings advantages as the best long-term predictability with respect to the expected neck bone loss.

Effective properties evaluation for smart composite materials

The purpose of this article is to present a method which consists in the development of unit cell numerical models for smart composite materials with piezoelectric fibers made of PZT embedded in a non-piezoelectric matrix (epoxy resin). This method evaluates a globally homogeneous medium equivalent to the original composite, using a representative volume element (RVE). The suitable boundary conditions allow the simulation of all modes of the overall deformation arising from any arbitrary combination of mechanical and electrical loading. In the first instance, the unit cell is applied to predict the effective material coefficients of the transversely isotropic piezoelectric composite with circular cross section fibers. The numerical results are compared to other methods reported in the literature and also to results previously published, in order to evaluate the method proposal. In the second step, the method is applied to calculate the equivalent properties for smart composite materials with square cross section fibers. Results of comparison between different combinations of circular and square fiber geometries, observing the influence of the boundary conditions and arrangements are presented.

Standard surgical treatment for benign prostatic hyperplasia is safe for patients over 75 years: analysis of 100 cases from a high-volume urologic center

OBJECTIVES: In this study, we aimed to determine the complications of standard surgical treatments among patients over 75 years in a high-volume urologic center. METHODS: We analyzed 100 consecutive patients older than 75 years who had undergone transurethral prostatic resection of the prostate or open prostatectomy for treatment of benign prostatic hyperplasia from January 2008 to March 2010. We analyzed patient age, prostate volume, prostate-specific antigen level, international prostatic symptom score, quality of life score, urinary retention, co-morbidities, surgical technique and satisfaction with treatment. RESULTS: Median age was 79 years. Forty-eight patients had undergone transurethral prostatic resection of the prostate, and 52 had undergone open prostatectomy. The median International Prostatic Symptom Score was 20, the median prostate volume was 83 g, 51% were using an indwelling bladder catheter, and the median prostatespecific antigen level was 5.0 ng/ml. The most common comorbidities were hypertension, diabetes and coronary disease. After a median follow-up period of 17 months, most patients were satisfied. Complications were present in 20% of cases. The most common urological complication was urethral stenosis, followed by bladder neck sclerosis, urinary fistula, late macroscopic hematuria and persistent urinary incontinence. The most common clinical complication was myocardial infarction, followed by acute renal failure requiring dialysis. Incidental carcinoma of the prostate was present in 6% of cases. One case had urothelial bladder cancer. CONCLUSIONS: Standard surgical treatments for benign prostatic hyperplasia are safe and satisfactory among the elderly. Complications are infrequent, and urethral stenosis is the most common. No clinical variable is associated with the occurrence of complications.

Longitudinal bar buckling behavior

Reinforced concrete columns might fail because of buckling of the longitudinal reinforcing bar when exposed to earthquake motions. Depending on the hoop stiffness and the length-over-diameter ratio, the instability can be local (in between two subsequent hoops) or global (the buckling length comprises several hoop spacings). To get insight into the topic, an extensive literary research of 19 existing models has been carried out including different approaches and assumptions which yield different results. Finite element fiberanalysis was carried out to study the local buckling behavior with varying length-over-diameter and initial imperfection-over-diameter ratios. The comparison of the analytical results with some experimental results shows good agreement before the post buckling behavior undergoes large deformation. Furthermore, different global buckling analysis cases were run considering the influence of different parameters; for certain hoop stiffnesses and length-over-diameter ratios local buckling was encountered. A parametric study yields an adimensional critical stress in function of a stiffness ratio characterized by the reinforcement configuration. Colonne in cemento armato possono collassare per via dell’instabilità dell’armatura longitudinale se sottoposte all’azione di un sisma. In funzione della rigidezza dei ferri trasversali e del rapporto lunghezza d’inflessione-diametro, l’instabilità può essere locale (fra due staffe adiacenti) o globale (la lunghezza d’instabilità comprende alcune staffe). Per introdurre alla materia, è proposta un’esauriente ricerca bibliografica di 19 modelli esistenti che include approcci e ipotesi differenti che portano a risultati distinti. Tramite un’analisi a fibre e elementi finiti si è studiata l’instabilità locale con vari rapporti lunghezza d’inflessione-diametro e imperfezione iniziale-diametro. Il confronto dei risultati analitici con quelli sperimentali mostra una buona coincidenza fino al raggiungimento di grandi spostamenti. Inoltre, il caso d’instabilità globale è stato simulato valutando l’influenza di vari parametri; per certe configurazioni di rigidezza delle staffe e lunghezza d’inflessione-diametro si hanno ottenuto casi di instabilità locale. Uno studio parametrico ha permesso di ottenere un carico critico adimensionale in funzione del rapporto di rigidezza dato dalle caratteristiche dell’armatura.

Models and methods for power monolithic microwave integrated circuits

Computer aided design of Monolithic Microwave Integrated Circuits (MMICs) depends critically on active device models that are accurate, computationally efficient, and easily extracted from measurements or device simulators. Empirical models of active electron devices, which are based on actual device measurements, do not provide a detailed description of the electron device physics. However they are numerically efficient and quite accurate. These characteristics make them very suitable for MMIC design in the framework of commercially available CAD tools. In the empirical model formulation it is very important to separate linear memory effects (parasitic effects) from the nonlinear effects (intrinsic effects). Thus an empirical active device model is generally described by an extrinsic linear part which accounts for the parasitic passive structures connecting the nonlinear intrinsic electron device to the external world. An important task circuit designers deal with is evaluating the ultimate potential of a device for specific applications. In fact once the technology has been selected, the designer would choose the best device for the particular application and the best device for the different blocks composing the overall MMIC. Thus in order to accurately reproducing the behaviour of different-in-size devices, good scalability properties of the model are necessarily required. Another important aspect of empirical modelling of electron devices is the mathematical (or equivalent circuit) description of the nonlinearities inherently associated with the intrinsic device. Once the model has been defined, the proper measurements for the characterization of the device are performed in order to identify the model. Hence, the correct measurement of the device nonlinear characteristics (in the device characterization phase) and their reconstruction (in the identification or even simulation phase) are two of the more important aspects of empirical modelling. This thesis presents an original contribution to nonlinear electron device empirical modelling treating the issues of model scalability and reconstruction of the device nonlinear characteristics. The scalability of an empirical model strictly depends on the scalability of the linear extrinsic parasitic network, which should possibly maintain the link between technological process parameters and the corresponding device electrical response. Since lumped parasitic networks, together with simple linear scaling rules, cannot provide accurate scalable models, either complicate technology-dependent scaling rules or computationally inefficient distributed models are available in literature. This thesis shows how the above mentioned problems can be avoided through the use of commercially available electromagnetic (EM) simulators. They enable the actual device geometry and material stratification, as well as losses in the dielectrics and electrodes, to be taken into account for any given device structure and size, providing an accurate description of the parasitic effects which occur in the device passive structure. It is shown how the electron device behaviour can be described as an equivalent two-port intrinsic nonlinear block connected to a linear distributed four-port passive parasitic network, which is identified by means of the EM simulation of the device layout, allowing for better frequency extrapolation and scalability properties than conventional empirical models. Concerning the issue of the reconstruction of the nonlinear electron device characteristics, a data approximation algorithm has been developed for the exploitation in the framework of empirical table look-up nonlinear models. Such an approach is based on the strong analogy between timedomain signal reconstruction from a set of samples and the continuous approximation of device nonlinear characteristics on the basis of a finite grid of measurements. According to this criterion, nonlinear empirical device modelling can be carried out by using, in the sampled voltage domain, typical methods of the time-domain sampling theory.

Safety factors in sheet pile wall design: considerations by using traditional methods and FEM analysis

The purpose of the work is: define and calculate a factor of collapse related to traditional method to design sheet pile walls. Furthermore, we tried to find the parameters that most influence a finite element model representative of this problem. The text is structured in this way: from chapter 1 to 5, we analyzed a series of arguments which are usefull to understanding the problem, while the considerations mainly related to the purpose of the text are reported in the chapters from 6 to 10. In the first part of the document the following arguments are shown: what is a sheet pile wall, what are the codes to be followed for the design of these structures and what they say, how can be formulated a mathematical model of the soil, some fundamentals of finite element analysis, and finally, what are the traditional methods that support the design of sheet pile walls. In the chapter 6 we performed a parametric analysis, giving an answer to the second part of the purpose of the work. Comparing the results from a laboratory test for a cantilever sheet pile wall in a sandy soil, with those provided by a finite element model of the same problem, we concluded that:in modelling a sandy soil we should pay attention to the value of cohesion that we insert in the model (some programs, like Abaqus, don’t accept a null value for this parameter), friction angle and elastic modulus of the soil, they influence significantly the behavior of the system (structure-soil), others parameters, like the dilatancy angle or the Poisson’s ratio, they don’t seem influence it. The logical path that we followed in the second part of the text is reported here. We analyzed two different structures, the first is able to support an excavation of 4 m, while the second an excavation of 7 m. Both structures are first designed by using the traditional method, then these structures are implemented in a finite element program (Abaqus), and they are pushed to collapse by decreasing the friction angle of the soil. The factor of collapse is the ratio between tangents of the initial friction angle and of the friction angle at collapse. At the end, we performed a more detailed analysis of the first structure, observing that, the value of the factor of collapse is influenced by a wide range of parameters including: the value of the coefficients assumed in the traditional method and by the relative stiffness of the structure-soil system. In the majority of cases, we found that the value of the factor of collapse is between and 1.25 and 2. With some considerations, reported in the text, we can compare the values so far found, with the value of the safety factor proposed by the code (linked to the friction angle of the soil).

Ensemble wind forecasting based on the HARMONIE model and adaptive finite elements in complex orography

[EN]Ensemble forecasting is a methodology to deal with uncertainties in the numerical wind prediction. In this work we propose to apply ensemble methods to the adaptive wind forecasting model presented in. The wind field forecasting is based on a mass-consistent model and a log-linear wind profile using as input data the resulting forecast wind from Harmonie, a Non-Hydrostatic Dynamic model used experimentally at AEMET with promising results. The mass-consistent model parameters are estimated by using genetic algorithms. The mesh is generated using the meccano method and adapted to the geometry…

A new insight in stress recovery tecniques and hessian reconstruction methods

Stress recovery techniques have been an active research topic in the last few years since, in 1987, Zienkiewicz and Zhu proposed a procedure called Superconvergent Patch Recovery (SPR). This procedure is a last-squares fit of stresses at super-convergent points over patches of elements and it leads to enhanced stress fields that can be used for evaluating finite element discretization errors. In subsequent years, numerous improved forms of this procedure have been proposed attempting to add equilibrium constraints to improve its performances. Later, another superconvergent technique, called Recovery by Equilibrium in Patches (REP), has been proposed. In this case the idea is to impose equilibrium in a weak form over patches and solve the resultant equations by a last-square scheme. In recent years another procedure, based on minimization of complementary energy, called Recovery by Compatibility in Patches (RCP) has been proposed in. This procedure, in many ways, can be seen as the dual form of REP as it substantially imposes compatibility in a weak form among a set of self-equilibrated stress fields. In this thesis a new insight in RCP is presented and the procedure is improved aiming at obtaining convergent second order derivatives of the stress resultants. In order to achieve this result, two different strategies and their combination have been tested. The first one is to consider larger patches in the spirit of what proposed in [4] and the second one is to perform a second recovery on the recovered stresses. Some numerical tests in plane stress conditions are presented, showing the effectiveness of these procedures. Afterwards, a new recovery technique called Last Square Displacements (LSD) is introduced. This new procedure is based on last square interpolation of nodal displacements resulting from the finite element solution. In fact, it has been observed that the major part of the error affecting stress resultants is introduced when shape functions are derived in order to obtain strains components from displacements. This procedure shows to be ultraconvergent and is extremely cost effective, as it needs in input only nodal displacements directly coming from finite element solution, avoiding any other post-processing in order to obtain stress resultants using the traditional method. Numerical tests in plane stress conditions are than presented showing that the procedure is ultraconvergent and leads to convergent first and second order derivatives of stress resultants. In the end, transverse stress profiles reconstruction using First-order Shear Deformation Theory for laminated plates and three dimensional equilibrium equations is presented. It can be seen that accuracy of this reconstruction depends on accuracy of first and second derivatives of stress resultants, which is not guaranteed by most of available low order plate finite elements. RCP and LSD procedures are than used to compute convergent first and second order derivatives of stress resultants ensuring convergence of reconstructed transverse shear and normal stress profiles respectively. Numerical tests are presented and discussed showing the effectiveness of both procedures.