An exact conserving algorithm for nonlinear dynamics with rotational DOFs and general hyperelasticity. Part 2: shells

Following the approach developed for rods in Part 1 of this paper (Pimenta et al. in Comput. Mech. 42:715-732, 2008), this work presents a fully conserving algorithm for the integration of the equations of motion in nonlinear shell dynamics. We begin with a re-parameterization of the rotation field in terms of the so-called Rodrigues rotation vector, allowing for an extremely simple update of the rotational variables within the scheme. The weak form is constructed via non-orthogonal projection, the time-collocation of which ensures exact conservation of momentum and total energy in the absence of external forces. Appealing is the fact that general hyperelastic materials (and not only materials with quadratic potentials) are permitted in a totally consistent way. Spatial discretization is performed using the finite element method and the robust performance of the scheme is demonstrated by means of numerical examples.

An exact conserving algorithm for nonlinear dynamics with rotational DOFs and general hyperelasticity. Part 1: Rods

A fully conserving algorithm is developed in this paper for the integration of the equations of motion in nonlinear rod dynamics. The starting point is a re-parameterization of the rotation field in terms of the so-called Rodrigues rotation vector, which results in an extremely simple update of the rotational variables. The weak form is constructed with a non-orthogonal projection corresponding to the application of the virtual power theorem. Together with an appropriate time-collocation, it ensures exact conservation of momentum and total energy in the absence of external forces. Appealing is the fact that nonlinear hyperelastic materials (and not only materials with quadratic potentials) are permitted without any prejudice on the conservation properties. Spatial discretization is performed via the finite element method and the performance of the scheme is assessed by means of several numerical simulations.

Lanczos algorithm on the grassmann manifold

5th. European Congress on Computational Methods in Applied Sciences and Engineering (ECCOMAS 2008) 8th. World Congress on Computational Mechanics (WCCM8)

Analysis of tongue shape and motion in speech production using statistical modeling

The mechanisms of speech production are complex and have been raising attention from researchers of both medical and computer vision fields. In the speech production mechanism, the articulator’s study is a complex issue, since they have a high level of freedom along this process, namely the tongue, which instigates a problem in its control and observation. In this work it is automatically characterized the tongues shape during the articulation of the oral vowels of Portuguese European by using statistical modeling on MR-images. A point distribution model is built from a set of images collected during artificially sustained articulations of Portuguese European sounds, which can extract the main characteristics of the motion of the tongue. The model built in this work allows under standing more clearly the dynamic speech events involved during sustained articulations. The tongue shape model built can also be useful for speech rehabilitation purposes, specifically to recognize the compensatory movements of the articulators during speech production.

A shock capturing scheme for steady, supercritical, free-surface flow

Abstract A finite difference scheme is presented for the solution of the two-dimensional shallow water equations in steady, supercritical flow. The scheme incorporates numerical characteristic decomposition, is shock capturing by design and incorporates space-marching as a result of the assumption that the flow is wholly supercritical in at least one space dimension. Results are shown for problems involving oblique hydraulic jumps and reflection from a wall.

An algorithm for the shallow water equations with body fitted meshes

An efficient algorithm based on flux difference splitting is presented for the solution of the two-dimensional shallow water equations in a generalised coordinate system. The scheme is based on solving linearised Riemann problems approximately and in more than one dimension incorporates operator splitting. The scheme has good jump capturing properties and the advantage of using body-fitted meshes. Numerical results are shown for flow past a circular obstruction.

A BEM formulation for analysing the coupled stretching-bending problem of plates reinforced by rectangular beams with columns defined in the domain

In this work, a boundary element formulation to analyse plates reinforced by rectangular beams, with columns defined in the domain is proposed. The model is based on Kirchhoff hypothesis and the beams are not required to be displayed over the plate surface, therefore eccentricity effects are taken into account. The presented boundary element method formulation is derived by applying the reciprocity theorem to zoned plates, where beams are treated as thin sub-regions with larger rigidities. The integral representations derived for this complex structural element consider the bending and stretching effects of both structural elements working together. The standard equilibrium and compatibility conditions along interface are naturally imposed, being the bending tractions eliminated along interfaces. The in-plane tractions and the bending and in-plane displacements are approximated along the beam width, reducing the number of degrees of freedom. The columns are introduced into the formulation by considering domain points where tractions can be prescribed. Some examples are then shown to illustrate the accuracy of the formulation, comparing the obtained results with other numerical solutions.

Nitrification processes in groundwaters evidenced by modification on the U contents and U-234/U-238 activity ratios

Groundwaters from the Botucatu-Piramboia aquifer system located at the Parana sedimentary basin in Brazil and sampled at Ribeirao Preto city, São Paulo State, were chemically and isotopically analysed, with the aim of evaluating how uranium isotopes U-234 and U-238 can be utilized to monitor the progress of nitrate, the main indicator of agrichemical pollution into an aquifer. Based on the isotope dilution analysis, the U-234/U-238 activity ratio and U content data in the groundwaters were utilized to deduce their proportions in a mixture, suggesting that the water not affected significantly by nitrification still dominates the exploitation system involving the studied wells.

Uranium in groundwaters from Botucatu-Piramboia aquifer, Brazil

Groundwater from 60 pumped tubular wells of the Botucatu-Piramboia aquifer system located at the Parana sedimentary basin in Brazil were chemically analysed with the aim of evaluating if the mechanisms related to the migration of uranium can generate concentrations greater than the maximum permissible limit in drinking water, as defined by the Brazilian national standards.

Health risk related to radon in a thermal spa utilizing waters from Guarani aquifer

This investigation was carried out within the Parana sedimentary basin, Brazil, involved the sampling of groundwater and air, and was realized with the purpose of evaluating the radioactivity, due to radon gas, in a thermal spa utilizing the waters from Guarani (Botucatu-Piramboia) aquifer. The results reported here provide additional information relative to that of previous studies focusing on the presence of radionuclides in the aquifer, which have mainly characterized those belonging to uranium and thorium series decay, such as the uranium isotopes (U-238 and U-234), radium isotopes (Ra-226 and Ra-228), radon daughters (Bi-214 and Pb-214) and radon (Rn-222) itself the results obtained were compared with the maximum permissible concentration limits in air and drinking water defined by international standards, such as the guidelines for drinking water quality established by the World Health Organization. The possible processes responsible for the presence of radon in the aquifer were also considered in order to evaluate the data obtained.

An implicit BEM formulation to model strong discontinuities in solids

A boundary element method (BEM) formulation to predict the behavior of solids exhibiting displacement (strong) discontinuity is presented. In this formulation, the effects of the displacement jump of a discontinuity interface embedded in an internal cell are reproduced by an equivalent strain field over the cell. To compute the stresses, this equivalent strain field is assumed as the inelastic part of the total strain. As a consequence, the non-linear BEM integral equations that result from the proposed approach are similar to those of the implicit BEM based on initial strains. Since discontinuity interfaces can be introduced inside the cell independently on the cell boundaries, the proposed BEM formulation, combined with a tracking scheme to trace the discontinuity path during the analysis, allows for arbitrary discontinuity propagation using a fixed mesh. A simple technique to track the crack path is outlined. This technique is based on the construction of a polygonal line formed by segments inside the cells, in which the assumed failure criterion is reached. Two experimental concrete fracture tests were analyzed to assess the performance of the proposed formulation.

Linear analysis of building floor structures by a BEM formulation based on Reissner's theory

In this work, the plate bending formulation of the boundary element method (BEM) based on the Reissner's hypothesis is extended to the analysis of zoned plates in order to model a building floor structure. In the proposed formulation each sub-region defines a beam or a slab and depending on the way the sub-regions are represented, one can have two different types of analysis. In the simple bending problem all sub-regions are defined by their middle surface. on the other hand, for the coupled stretching-bending problem all sub-regions are referred to a chosen reference surface, therefore eccentricity effects are taken into account. Equilibrium and compatibility conditions are automatically imposed by the integral equations, which treat this composed structure as a single body. The bending and stretching values defined on the interfaces are approximated along the beam width, reducing therefore the number of degrees of freedom. Then, in the proposed model the set of equations is written in terms of the problem values on the beam axis and on the external boundary without beams. Finally some numerical examples are presented to show the accuracy of the proposed model.

Environmental education and the health risk question: how do we educate the educators working with this theme?

Our society is presently seeing a trend of continuous increase of risk with severe consequences, especially those linked to environmental aspects. It has obligated us to reflect on this production and its effects, either positive or negative, and, after a long reflection, to think about the adopted model of civilization and the possibilities of changing such a model. The educational institution is among the social departments responsible for this revision process, since it represents one possibility to promote a more critical vision from society, providing it with tools needed for a more responsible action. However, schools are dominated by a traditional work, which does not allow for enough and proper attention to these relevant and current questions, and which involve a critical analysis of the conceptions and values established by our civilization, resulting in the present patterns of science-technology-society relations and their influence on the environment and health. Several research works have indicated a diversity of barriers that obstruct this necessary change. Since formal education plays an important role in the education of society for this theme, in this paper the possibilities and difficulty dealing with this problem are discussed, focusing on the question of teacher education, based on investigations related to teacher education programs for environmental education. The data collected show the urgency for introducing this theme in teacher training programs, taking into account the fact that educators have their own conceptions, values and attitudes that should be considered in these programs.