Simulación de fluidos inmiscibles con el Particle Finite Element Method (PFEM)

Proyecto de investigación realizado a partir de una estancia en el Centro Internacional de Métodos Computacionales en Ingeniería (CIMEC), Argentina, entre febrero y abril del 2007. La simulación numérica de problemas de mezclas mediante el Particle Finite Element Method (PFEM) es el marco de estudio de una futura tesis doctoral. Éste es un método desarrollado conjuntamente por el CIMEC y el Centre Internacional de Mètodos Numèrics en l'Enginyeria (CIMNE-UPC), basado en la resolución de las ecuaciones de Navier-Stokes en formulación Lagrangiana. El mallador ha sido implementado y desarrollado por Dr. Nestor Calvo, investigador del CIMEC. El desarrollo del módulo de cálculo corresponde al trabajo de tesis de la beneficiaria. La correcta interacción entre ambas partes es fundamental para obtener resultados válidos. En esta memoria se explican los principales aspectos del mallador que fueron modificados (criterios de refinamiento geométrico) y los cambios introducidos en el módulo de cálculo (librería PETSc, algoritmo predictor-corrector) durante la estancia en el CIMEC. Por último, se muestran los resultados obtenidos en un problema de dos fluidos inmiscibles con transferencia de calor.

Simulated joint and muscle forces in reversed and anatomic shoulder prostheses.

Reversed shoulder prostheses are increasingly being used for the treatment of glenohumeral arthropathy associated with a deficient rotator cuff. These non-anatomical implants attempt to balance the joint forces by means of a semi-constrained articular surface and a medialised centre of rotation. A finite element model was used to compare a reversed prosthesis with an anatomical implant. Active abduction was simulated from 0 degrees to 150 degrees of elevation. With the anatomical prosthesis, the joint force almost reached the equivalence of body weight. The joint force was half this for the reversed prosthesis. The direction of force was much more vertically aligned for the reverse prosthesis, in the first 90 degrees of abduction. With the reversed prosthesis, abduction was possible without rotator cuff muscles and required 20% less deltoid force to achieve it. This force analysis confirms the potential mechanical advantage of reversed prostheses when rotator cuff muscles are deficient.

Early stage disc degeneration does not have an appreciable affect on stiffness and load transfer following vertebroplasty and kyphoplasty.

Vertebroplasty and kyphoplasty have been reported to alter the mechanical behavior of the treated and adjacent-level segments, and have been suggested to increase the risk for adjacent-level fractures. The intervertebral disc (IVD) plays an important role in the mechanical behavior of vertebral motion segments. Comparisons between normal and degenerative IVD motion segments following cement augmentation have yet to be reported. A microstructural finite element model of a degenerative IVD motion segment was constructed from micro-CT images. Microdamage within the vertebral body trabecular structure was used to simulate a slightly (I = 83.5% of intact stiffness), moderately (II = 57.8% of intact stiffness), and severely (III = 16.0% of intact stiffness) damaged motion segment. Six variable geometry single-segment cement repair strategies (models A-F) were studied at each damage level (I-III). IVD and bone stresses, and motion segment stiffness, were compared with the intact and baseline damage models (untreated), as well as, previous findings using normal IVD models with the same repair strategies. Overall, small differences were observed in motion segment stiffness and average stresses between the degenerative and normal disc repair models. We did however observe a reduction in endplate bulge and a redistribution in the microstructural tissue level stresses across both endplates and in the treated segment following early stage IVD degeneration. The cement augmentation strategy placing bone cement along the periphery of the vertebra (model E) proved to be the most advantageous in treating the degenerative IVD models by showing larger reductions in the average bone stresses (vertebral and endplate) as compared to the normal IVD models. Furthermore, only this repair strategy, and the complete cement fill strategy (model F), were able to restore the slightly damaged (I) motion segment stiffness above pre-damaged (intact) levels. Early stage IVD degeneration does not have an appreciable effect in motion segment stiffness and average stresses in the treated and adjacent-level segments following vertebroplasty and kyphoplasty. Placing bone cement in the periphery of the damaged vertebra in a degenerative IVD motion segment, minimizes load transfer, and may reduce the likelihood of adjacent-level fractures.

Nonlinear Finite Element Study of Piles in Integral Abutment Bridges, HR-227, Part 2, 1982

The highway departments of the states which use integral abutments in bridge design were contacted in order to study the extent of integral abutment use in skewed bridges and to survey the different guidelines used for analysis and design of integral abutments in skewed bridges. The variation in design assumptions and pile orientations among the various states in their approach to the use of integral abutments on skewed bridges is discussed. The problems associated with the treatment of the approach slab, backfill, and pile cap, and the reason for using different pile orientations are summarized in the report. An algorithm based on a state-of-the-art nonlinear finite element procedure previously developed by the authors was modified and used to study the influence of different factors on behavior of piles in integral abutment bridges. An idealized integral abutment was introduced by assuming that the pile is rigidly cast into the pile cap and that the approach slab offers no resistance to lateral thermal expansion. Passive soil and shear resistance of the cap are neglected in design. A 40-foot H pile (HP 10 X 42) in six typical Iowa soils was analyzed for fully restrained pile head and pinned pile head. According to numerical results, the maximum safe length for fully restrained pile head is one-half the maximum safe length for pinned pile head. If the pile head is partially restrained, the maximum safe length will lie between the two limits. The numerical results from an investigation of the effect of predrilled oversized holes indicate that if the length of the predrilled oversized hole is at least 4 feet below the ground, the vertical load-carrying capacity of the H pile is only reduced by 10 percent for 4 inches of lateral displacement in very stiff clay. With no predrilled oversized hole, the pile failed before the 4-inch lateral displacement was reached. Thus, the maximum safe lengths for integral abutment bridges may be increased by predrilling. Four different typical Iowa layered soils were selected and used in this investigation. In certain situations, compacted soil (> 50 blow count in standard penetration tests) is used as fill on top of natural soil. The numerical results showed that the critical conditions will depend on the length of the compacted soil. If the length of the compacted soil exceeds 4 feet, the failure mechanism for the pile is similar to one in a layer of very stiff clay. That is, the vertical load-carrying capacity of the H pile will be greatly reduced as the specified lateral displacement increases.

Statistical properties of population differentiation estimators under stepwise mutation in a finite island model.

Microsatellite loci mutate at an extremely high rate and are generally thought to evolve through a stepwise mutation model. Several differentiation statistics taking into account the particular mutation scheme of the microsatellite have been proposed. The most commonly used is R(ST) which is independent of the mutation rate under a generalized stepwise mutation model. F(ST) and R(ST) are commonly reported in the literature, but often differ widely. Here we compare their statistical performances using individual-based simulations of a finite island model. The simulations were run under different levels of gene flow, mutation rates, population number and sizes. In addition to the per locus statistical properties, we compare two ways of combining R(ST) over loci. Our simulations show that even under a strict stepwise mutation model, no statistic is best overall. All estimators suffer to different extents from large bias and variance. While R(ST) better reflects population differentiation in populations characterized by very low gene-exchange, F(ST) gives better estimates in cases of high levels of gene flow. The number of loci sampled (12, 24, or 96) has only a minor effect on the relative performance of the estimators under study. For all estimators there is a striking effect of the number of samples, with the differentiation estimates showing very odd distributions for two samples.

Impact of the Femoral Head Position on Moment Arms in Total Hip Arthroplasty: A Parametric Finite Element Study.

BACKGROUND: Although the importance of accurate femoral reconstruction to achieve a good functional outcome is well documented, quantitative data on the effects of a displacement of the femoral center of rotation on moment arms are scarce. The purpose of this study was to calculate moment arms after nonanatomical femoral reconstruction. METHODS: Finite element models of 15 patients including the pelvis, the femur, and the gluteal muscles were developed. Moment arms were calculated within the native anatomy and compared to distinct displacement of the femoral center of rotation (leg lengthening of 10 mm, loss of femoral offset of 20%, anteversion ±10°, and fixed anteversion at 15°). Calculations were performed within the range of motion observed during a normal gait cycle. RESULTS: Although with all evaluated displacements of the femoral center of rotation, the abductor moment arm remained positive, some fibers initially contributing to extension became antagonists (flexors) and vice versa. A loss of 20% of femoral offset led to an average decrease of 15% of abductor moment. Femoral lengthening and changes in femoral anteversion (±10°, fixed at 15°) led to minimal changes in abductor moment arms (maximum change of 5%). Native femoral anteversion correlated with the changes in moment arms induced by the 5 variations of reconstruction. CONCLUSION: Accurate reconstruction of offset is important to maintaining abductor moment arms, while changes of femoral rotation had minimal effects. Patients with larger native femoral anteversion appear to be more susceptible to femoral head displacements.

Modelling syntectonic sedimentation: combining a discrete element model of tectonic. Deformation and process-based sedimentary Model in 3D.

This paper presents a new numerical program able to model syntectonic sedimentation. The new model combines a discrete element model of the tectonic deformation of a sedimentary cover and a process-based model of sedimentation in a single framework. The integration of these two methods allows us to include the simulation of both sedimentation and deformation processes in a single and more effective model. The paper describes briefly the antecedents of the program, Simsafadim-Clastic and a discrete element model, in order to introduce the methodology used to merge both programs to create the new code. To illustrate the operation and application of the program, analysis of the evolution of syntectonic geometries in an extensional environment and also associated with thrust fault propagation is undertaken. Using the new code, much more complex and realistic depositional structures can be simulated together with a more complex analysis of the evolution of the deformation within the sedimentary cover, which is seen to be affected by the presence of the new syntectonic sediments.

Numerical simulation of two-chamber plasma sources using finite element method

Numerical simulation of plasma sources is very important. Such models allows to vary different plasma parameters with high degree of accuracy. Moreover, they allow to conduct measurements not disturbing system balance.Recently, the scientific and practical interest increased in so-called two-chamber plasma sources. In one of them (small or discharge chamber) an external power source is embedded. In that chamber plasma forms. In another (large or diffusion chamber) plasma exists due to the transport of particles and energy through the boundary between chambers.In this particular work two-chamber plasma sources with argon and oxygen as active mediums were onstructed. This models give interesting results in electric field profiles and, as a consequence, in density profiles of charged particles.

Simulation and material calibration of ultra high strength steel (UHSS) S960 welded joint under static tensile test utilizing finite element method (FEM)

The aim of this work was to calibrate the material properties including strength and strain values for different material zones of ultra-high strength steel (UHSS) welded joints under monotonic static loading. The UHSS is heat sensitive and softens by heat due to welding, the affected zone is heat affected zone (HAZ). In this regard, cylindrical specimens were cut out from welded joints of Strenx® 960 MC and Strenx® Tube 960 MH, were examined by tensile test. The hardness values of specimens’ cross section were measured. Using correlations between hardness and strength, initial material properties were obtained. The same size specimen with different zones of material same as real specimen were created and defined in finite element method (FEM) software with commercial brand Abaqus 6.14-1. The loading and boundary conditions were defined considering tensile test values. Using initial material properties made of hardness-strength correlations (true stress-strain values) as Abaqus main input, FEM is utilized to simulate the tensile test process. By comparing FEM Abaqus results with measured results of tensile test, initial material properties will be revised and reused as software input to be fully calibrated in such a way that FEM results and tensile test results deviate minimum. Two type of different S960 were used including 960 MC plates, and structural hollow section 960 MH X-joint. The joint is welded by BöhlerTM X96 filler material. In welded joints, typically the following zones appear: Weld (WEL), Heat affected zone (HAZ) coarse grained (HCG) and fine grained (HFG), annealed zone, and base material (BaM). Results showed that: The HAZ zone is softened due to heat input while welding. For all the specimens, the softened zone’s strength is decreased and makes it a weakest zone where fracture happens while loading. Stress concentration of a notched specimen can represent the properties of notched zone. The load-displacement diagram from FEM modeling matches with the experiments by the calibrated material properties by compromising two correlations of hardness and strength.