A finite strain contact model for mixed lubricated surfaces in forming processes

For an accurate simulation of forming processes, it is of paramount importance to model the different lubrication regimes that can develop at the contact interface. These might vary from zone to zone of the forming piece, and from one regime to another, resulting in forces of different nature and magnitude. In these cases, the use of the classical Coulomb friction law will be clearly not sufficient to capture, in a suitable manner, the variety of forces applied on the forming piece.     

Finite element simulation for dynamic large elastoplastic deformation in metal powder forming

In this paper, a transient dynamic analysis of the powder compaction process is simulated by a large displacement finite element method based on a total and updated Lagrangian formulation. A combination of the Mohr–Coulomb and elliptical yield cap model, which reflects the stress state and degree of densification, is applied to describe the constitutive model of powder materials. A Coulomb friction law and a plasticity theory of friction in the context of an interface element formulation are employed in the constitutive modelling of the frictional behaviour between the die and powder. Finally, the powder behaviour during the compaction of a plain bush, a rotational flanged and a shaped tip component are analysed numerically. It is shown that the updated Lagrangian formulation, using a combination of the Mohr–Coulomb and elliptical cap model, can be effective in simulating metal powder compaction.     

Elasto-plastic finite element modelling of strip cold rolling using Eulerian fixed mesh technique

An Eulerian fixed mesh finite element technique applicable to metal-forming processes operating under steady-state condition is presented. Different specific features are demonstrated by solving plane-strain rolling problem. The advantage of the Eulerian fixed mesh technique over the updated Lagrangian one in modelling the elastic flattening of rolls is demonstrated. The obtained pressure distribution and the stress field are compared with other numerical and/or experimental results available in the literature with which good agreement is found. It is found that the consideration of the elastic flattening of rolls decreases the difference between the measured and the computed results.     

Numerical simulation of pollutant dispersion in street canyons: Geometric and thermal effects

Numerical investigations on pollutant dispersion in street canyons with emission sources located near the ground level are performed in the present work. Pollutant dispersion problems in urban areas are usually studied considering the street canyon model, which consists of long streets laterally confined by buildings. Significant changes can be observed in wind flow patterns and pollutant concentration fields when thermal and geometric effects are considered. Thus, the objective of this study is to investigate numerically the wind flow and pollutant dispersion for the following cases: (a) a two-dimensional street canyon model considering three different aspect ratios and four different wall heating configurations; (b) a flow domain with two immersed buildings arranged in two distinct configurations; (c) a three-dimensional urban area model composed of a building set and street intersections. Expected flow structures were obtained inside the canyon when different aspect ratios and wall heating configurations were considered. Flow phenomena such as separation/reattachment were observed when two-buildings models were analyzed. Finally, three-dimensional flow structures, with some characteristic that are not observed in two-dimensional models, affecting the pollutant removal, were simulated in the last case, highlighting the relevance of model dimensionality. The wind flow and pollutant dispersion are investigated using a numerical model based on the finite element formulation utilized by some of the authors of this work, which is extended here to deal with problems of heat and mass transport in the urban micro-scale. Turbulence is reproduced using Large Eddy Simulation (LES) and thermal effects on the momentum equations are considered as a buoyancy force, according to Boussinesq approximation.     

A jaw model for the study of the mandibular flexure taking into account the anisotropy of the bone

Knowledge of the complex biomechanical behaviour of the human mandible is of great importance in various clinical situations. The biomechanical and physical behaviour of mandibles have been investigated by different approaches. Some research have been done to evaluate the functional character of mandibles. Methods such as indirect measurement of deformations performed by intraoral appliances and by holographic interferometry have being employed. Other studies evaluated the mechanical properties and material parameters of small cubes of mandibles. One disadvantage of the experiments using strain gauges or holographic interferometry is the inability to determine strains at defined positions within the specimen. Additionally, research in biomechanics by these methods is limited to surface deformations and neither stresses nor dislocations can be measured.In the course of this study, we have investigated the mandibular flexure under mechanical loads using the results of a Finite Element Analysis (FEA). In order to obtain more accurate and realistic results, the bone anisotropy has being taken into account for the mathematical modelling of the jaw.The objective of this study was to establish a non-invasive procedure to predict precisely the complex biomechanical reactions of mandibles under mechanical loading. In order to achieve this aim, a comparison of the numerical data obtained with the experimental values of previous studies was performed. It showed a good correlation between in vitro measurements and mathematical modelling. Then the Finite Elements (FE) model was used to evaluate some mandibular movements (corporal approximation, dorsoventral shear, and corporal rotation in edentulous subjects).It is concluded that the applied procedure of generating the FE model is a valid and accurate non-invasive method to predict different parameters of the complex biomechanical behaviour of human mandibles.     

髋关节发育不良髋臼假体位置与稳定的三维有限元分析

目的利用三维有限元仿真的方法，研究全髋关节置换术治疗成人髋关节发育不良时3 种不同髋臼假体位置与髋臼 假体初始稳定的相关性，为临床手术操作提供参考依据。方法采用专门的生物力学有限元网格划分器，从髋关节发育不良的CT 扫描数据中，建立高度仿真的个性化髋骨三维有限元模型，并模拟在真臼位置安放髋臼假体的解剖位重建、将髋臼内壁打磨穿透内移的中心化重建、和在真性髋臼上方假臼高位重建3种常见的临床手术方案，模拟不同臼杯假体位置对全髋关节置换术后髋臼假体稳定的影响。结果髋臼假体在真臼位置安放的解剖重建，出现应力集中和大剪切应力的可能性最小，穿透内移安放的模拟应力分布略高于真臼重建，而在假体上移高位安放的模拟则出现预测应力大幅增加。结论成人髋关节发育不良全髋关节置换术中应尽可能在真臼位置重建和安放臼杯假体，以达到最优的应力分布和稳定性。     

Comparison of electromigration simulation in test structure and actual circuit

With the rapid increase in circuit complexity, accurate reliability simulator is necessary as reliability testing for fabricated circuit is progressively more time and resource consuming. Simple 2D electromigration (EM) simulator has many limitations and most of the 3D EM simulations in the literature are using simple line-via structure which is only part of the real circuit structure. In this paper, a thorough comparison between simple 3D line-via structure model and 3D circuit model under different test conditions is performed using class-AB amplifier as an example, and the results showed that simple line-via structure simulation may not always correctly represent the real circuit failure site.     

Nonlinear analysis of residual stresses in a rail manufacturing process by FEM

The aim of this paper is to study the residual stresses in an UIC-60 rail and their reduction by means of roller straightening. Both experimental and numerical investigations have been carried out in the past to reveal the formation of dominant longitudinal residual stresses. However, the agreement between both investigations was not particularly good. The finite element method (FEM) has also been used to simulate one, two and three-dimensional analyses of a rail during roller straightening processes. The present model considers the longitudinal movement of a rail through the straightening machine, contact conditions between rail and rollers and kinematic hardening so as to take into account the plastic behaviour of the rail material (steel). These results were compared with the experimental investigations and good agreement was observed. In this respect, this paper presents a novel, more realistic numerical simulation by FEM for the roller straightening process. Finally, an improvement of the straightening process in order to obtain smaller residual stress in the rail section is proposed.     

Modeling and analysis of nonlinear rotordynamics due to higher order deformations in bending

A mathematical model incorporating the higher order deformations in bending is developed and analyzed to investigate the nonlinear dynamics of rotors. The rotor system considered for the present work consists of a flexible shaft and a rigid disk. The shaft is modeled as a beam with a circular cross section and the Euler Bernoulli beam theory is applied with added effects such as rotary inertia, gyroscopic effect, higher order large deformations, rotor mass unbalance and dynamic axial force. The kinetic and strain (deformation) energies of the rotor system are derived and the Rayleigh–Ritz method is used to discretize these energy expressions. Hamilton’s principle is then applied to obtain the mathematical model consisting of second order coupled nonlinear differential equations of motion. In order to solve these equations and hence obtain the nonlinear dynamic response of the rotor system, the method of multiple scales is applied. Furthermore, this response is examined for different possible resonant conditions and resonant curves are plotted and discussed. It is concluded that nonlinearity due to higher order deformations significantly affects the dynamic behavior of the rotor system leading to resonant hard spring type curves. It is also observed that variations in the values of different parameters like mass unbalance and shaft diameter greatly influence dynamic response. These influences are also presented graphically and discussed.     

FE simulation of human trachea swallowing movement before and after the implantation of an endoprothesis

Background and aims

Nowadays interventions associated to the implantation of tracheal stents in patients with airway pathologies, are a very common surgery that in the post-operating period can bear many problems such as migration of the stent, development of granulation tissue at the edges of the stent with overgrowth of the tracheal lumen, or accumulation of secretions inside the prosthesis. Among the movements that trachea carries out, swallowing seems to drive harmful consequences for the tracheal tissues surrounding the prosthesis. In this work a finite element model of a human trachea has been developed and used to analyze its behavior during swallowing.

Material and methods

In the present work, a complete human trachea finite element model based on experimental study was developed. The real swallowing movement of two patients before and after the implantation of Dumon prosthesis was used to simulate and then analyze the effect that the tracheal implant has on the stress response of the trachea and on the physiological capacity to swallow.

Results

In both studied cases with an implanted Dumon prosthesis, patients showed a decrease of their ability to swallow; one lost 26.4% and the other one 18.9% of their tracheal ascending movements. Besides, the prosthesis implantation caused an increase of the stresses located in the superior contact border between the tracheal wall and the prosthesis. It could be seen that the resulting force equivalent to the elevating tracheal muscle forces for degluting, was around F = 13.5 N for the two patients both before and after the stent implantation.

Conclusion

The implantation of a Dumon prosthesis modifies the mechanical response of the trachea altering its stress distribution and its ascending movement.     

Modelling and numerical simulation of high strength fibre reinforced concrete corbels

In the present study, new constitutive models for high strength steel fibre reinforced concrete (HSSFRC) have been formulated by means of a regression analysis of many experimental data (from literature) by using SPSS-statistical program. This proposed constitutive models have been employed for formulating the material finite element models to study the behaviour of HSSFRC corbels.     

Object-oriented modelling and simulation of heat exchangers with finite element methods

The paper describes the derivation of finite-element models of one-dimensional fluid flows with heat transfer in pipes, using the Galerkin/least-squares approach. The models are first derived for one-phase flows, and then extended to homogeneous two-phase flows. The resulting equations have then been embedded in the context of object-oriented system modelling; this allows one to combine the fluid flow model with a model for other phenomena such as heat transfer, as well as with models of other discrete components such as pumps or valves, to obtain complex models of heat exchangers. The models are then validated by simulating a typical heat exchanger plant.     

The material parameter design and finite element simulation of the quadrilateral thermal cloak device

In this paper, we first design a coordinate transformation and derive the anisotropic material parameters of the quadrilateral thermal cloak according to the transformation thermodynamics principle. Then, since the derived parameters are inherently anisotropic, we eliminate its anisotropy by considering the effective medium theory and use a layered structure of metamaterials composed of only two isotropic materials to design the cloak device. Finally, we simulate the performance of a perfect and layered thermal cloak by the finite element method. To the best of our knowledge, this is the first work to design and simulate the performance of this quadrilateral thermal cloak by the finite element method(FEM).     

Numerical analysis and modeling of large deformation and necking behavior of tensile specimens

The present paper is concerned with numerical simulations of large deformation and necking behavior of axisymmetric tensile specimens. In particular, an efficient framework for the numerical analysis of finite deformation behavior of elastic-rate-independent plastic problems is summarized which is based on a plastic predictor method. Furthermore, numerical modeling of conventional tensile tests as well as their finite deformation and localization phenomena are discussed in some detail, and the results will be compared to those obtained by simplified numerical simulations.     

Mathematical modeling and numerical simulation of CO2 transport through hollow-fiber membranes

Chemical absorption of carbon dioxide was studied theoretically using hollow-fiber membrane contactors in this work. A 2D mathematical model was developed to study CO₂ transport through hollow-fiber membrane contactors. The model considers axial and radial diffusion in the membrane contactor. It also considers convection in the tube and shell side with chemical reaction. The finite element method (FEM) was used to solve the model equations. Modeling predictions were validated with the experimental data obtained from literature for CO₂ absorption in amine aqueous solutions as solvent. The modeling predictions were in good agreement with the experimental data for different values of gas and liquid velocities. The liquid solvents considered for this study include aqueous solutions of monoethanolamine (MEA), diethanolamine (DEA), N-methyldiethanolamine (MDEA), 2-amino-2-methyl-1-propanol (AMP) and potassium carbonate (K₂CO₃). The simulation results indicated that amine aqueous solutions were better than K₂CO₃ aqueous solution for CO₂ absorption. Also simulation results revealed that the removal of CO₂ with aqueous solution of MEA was the highest among the amines solvents. The hollow-fiber membrane contactors showed a great potential in the area of CO₂ absorption.     

Numerical investigation of the deformation characteristics and heat generation in pneumatic aircraft tires : Part I. Mechanical modeling

A numerical methodology to study the temperature rise in aircraft tires under free rolling conditions is presented in this article. In the first part, we study the deformation characteristics of the tire to determine the heat generation due to the inelastic deformation (viscoelasticity), which is assumed to be the main source of heat generation under free rolling conditions. The heat generation is then used as an input to solve the heat transfer problem which is addressed in the second part of this article. A methodology which considers a 2-D formulation with the contribution of out-of-plane forces is presented. This methodology allows for a significant reduction of the computational requirements of 3-D analysis while capturing the main 3-D effects. The deformation pattern as well as the stress and strain fields are presented for a cross-section of the tire at several locations.     

Identification of diffusion parameters in welded joints of titanium and its alloys

Titanium and its alloys may experience a progressive embrittlement leading to catastrophic failure. Such embrittlement is a consequence of hydrogen migration and the accumulation of brittle titanium hydrides in the presence of a sharp alloy composition gradient near a weld fusion line. The governing equations for such migration are presented. The identification problem for parameters responsible for this phenomenon is formulated. The sensitivity equations are given in terms of finite elements to solve the problem of identification. Based on numerical procedures the parameters of welded joints due to hydrogen diffusion are identified.     

Modeling and simulation of transient responses of a flexible beam floating in finite depth water under moving loads

The dynamic response of a free–free flexible beam floating in an unbounded water domain under the effect of moving loads is numerically analyzed. The water is assumed compressible and inviscid. The surface disturbance satisfies a linear free surface wave condition and an undisturbed condition at infinity. In the present work, a finite element procedure was developed directly in time domain and implemented to solve the two-dimensional problem of the transient behavior of an elastic beam floating on the surface of finite deep water under the passage of a moving force with uniform speed. The presented data demonstrates the applicability of the proposed mathematical model and numerical approach. The influences on the dynamic responses of floating beam of some factors were studied.     

Numerical investigation of the deformation characteristics and heat generation in pneumatic aircraft tires Part II. Thermal modeling

A numerical procedure to determine the temperature rise in aircraft tires under free rolling conditions is presented in this article. Energy dissipation from cyclic inelastic deformation is considered the main heat generation source. This modeling considers the deformation process of the tire to be a steady-state problem, where all concurrent cycles are assumed to be the same as the first. The inelastic energy is determined by imposing a phase lag between the strain and the stress fields. The phase lag is assumed to be frequency independent in the range of interest, in keeping with the experimental observations in aircraft tire materials. It is further assumed that the inelastic energy is completely converted into volumetric heat input for a transient thermal conduction analysis. A conduction model is described and results are compared against thermocouple data recorded by Clark and Dodge [1].     

Nonlinear explicit analysis and study of the behaviour of a new ring-type brake energy dissipator by FEM and experimental comparison

The aim of this paper is to comprehensively analyse the performance of a new ring-type brake energy dissipator through the finite element method (FEM) (formulation and finite element approximation of contact in nonlinear mechanics) and experimental comparison. This new structural device is used as a system component in rockfall barriers and fences and it is composed of steel bearing ropes, bent pipes and aluminium compression sleeves. The bearing ropes are guided through pipes bent into double-loops and held by compression sleeves. These elements work as brake rings. In important events the brake rings contract and so dissipate residual energy out of the ring net, without damaging the ropes. The rope’s breaking load is not diminished by activation of the brake. The full understanding of this problem implies the simultaneous study of three nonlinearities: material nonlinearity (plastic behaviour) and failure criteria, large displacements (geometric nonlinearity) and friction-contact phenomena among brake ring components. The explicit dynamic analysis procedure is carried out by means of the implementation of an explicit integration rule together with the use of diagonal element mass matrices. The equations of motion for the brake ring are integrated using the explicit central difference integration rule. The presence of the contact phenomenon implies the existence of inequality constraints. The conditions for normal contact are and gλ=0, where λ is the normal traction component and g is the gap function for the contact surface pair. To include frictional conditions, let us assume that Coulomb’s law of friction holds pointwise on the different contact surfaces, μ being the dynamic coefficient of friction. Next, we define the non-dimensional variable τ by means of the expression τ=t/μλ, where μλ is the frictional resistance and t is the tangential traction component. In order to find the best brake performance, different dynamic friction coefficients corresponding to the pressures of the compression sleeves have been adopted and simulated numerically by FEM and then we have compared them with the results from full-scale experimental tests. Finally, the most important conclusions of this study are given.