Development of a Three-Dimensional Semi-Analytical Elastic-Plastic Contact Code

2002 ◽  
Vol 124 (4) ◽  
pp. 653-667 ◽  
Author(s):  
C. Jacq ◽  
D. Ne´lias ◽  
G. Lormand ◽  
D. Girodin

A three-dimensional elastic-plastic contact code based on semi-analytical method is presented and validated. The contact is solved within a Hertz framework. The reciprocal theorem with initial strains is then introduced, to express the surface geometry as a function of contact pressure and plastic strains. The irreversible nature of plasticity leads to an incremental formulation of the elastic-plastic contact problem, and an algorithm to solve this problem is set up. Closed form expression, which give residual stresses and surface displacements from plastic strains, are obtained by integration of the reciprocal theorem. The resolution of the elastic-plastic contact using the finite element (FE) method is discussed, and the semi-analytical code presented in this paper is validated by comparing results with experimental data from the nano-indentation test. Finally, the resolution of the rolling elastic-plastic contact is presented for smooth and dented surfaces and for a vertical or rolling loading. The main advantage of this code over classical FE codes is that the calculation time makes the transient analysis of three-dimensional contact problems affordable, including when a fine mesh is required.

Author(s):  
Eduard Antaluca ◽  
Daniel Ne´lias ◽  
Spiridon Cretu

A three-dimensional numerical model based on a semi-analytical method in the framework of small strains and small displacements with respect of Hertz’s hypotheses is presented for solving an elastic-plastic dented contact with friction. The calculation of surface deformations and pressure distribution, which is the most time consuming step during the elastic-plastic algorithm, is obtained using a method based on a variational principle with a Fast Fourier Transform (FFT) and a Conjugate Gradient Method (CGM). The method is fast enough to allow investigating the effect of a small size surface defect, here a debris denting, on the subsurface elastic-plastic stress state, requiring a fine mesh with around 106 surface grid points. Further, the FFT approach is also involved in the calculation of internal stress state. The plasticity model is based on an incremental load and Von Mises yield criterion. The effects of the contact pressure distribution and residual strain on the geometry of the contacting surfaces yield from the Betti’s reciprocal theorem with initial strain. The code is used to compute a few smooth and dented contacts, with several types of contact interfaces conditions, including frictionless and Coulomb friction. The effects of surface dents and friction on the contact pressure and subsurface stress field are presented and discussed.


2007 ◽  
Vol 129 (4) ◽  
pp. 761-771 ◽  
Author(s):  
Daniel Nélias ◽  
Eduard Antaluca ◽  
Vincent Boucly ◽  
Spiridon Cretu

A three-dimensional numerical model based on a semianalytical method in the framework of small strains and small displacements is presented for solving an elastic-plastic contact with surface traction. A Coulomb’s law is assumed for the friction, as commonly used for sliding contacts. The effects of the contact pressure distribution and residual strain on the geometry of the contacting surfaces are derived from Betti’s reciprocal theorem with initial strain. The main advantage of this approach over the classical finite element method (FEM) is the computing time, which is reduced by several orders of magnitude. The contact problem, which is one of the most time-consuming procedures in the elastic-plastic algorithm, is obtained using a method based on the variational principle and accelerated by means of the discrete convolution fast Fourier transform (FFT) and conjugate gradient methods. The FFT technique is also involved in the calculation of internal strains and stresses. A return-mapping algorithm with an elastic predictor∕plastic corrector scheme and a von Mises criterion is used in the plasticity loop. The model is first validated by comparison with results obtained by the FEM. The effect of the friction coefficient on the contact pressure distribution, subsurface stress field, and residual strains is also presented and discussed.


2007 ◽  
Vol 51 (02) ◽  
pp. 128-136
Author(s):  
Gonghyun Jung

A new numerical model, Q-Weld (Edison Welding Institute, Columbus, OH), which is a shell-element-based elastic analysis, is proposed for the prediction of the distortion induced in ship panels. Based on the results of the three-dimensional thermal-elastic-plastic analyses, it was found that the shell element-based model excluding the geometry of fillet welds and including only transverse and longitudinal plastic strains is valid without significant loss of accuracy. The developed Q-Weld predicts well-agreed distortions with the three-dimensional thermal-elastic-plastic analysis and demonstrates its potential in welding-induced distortion analysis, including buckling analysis.


2011 ◽  
Vol 44 (6) ◽  
pp. 1211-1216 ◽  
Author(s):  
Masoud Allahkarami ◽  
Jay C. Hanan

An efficient method for X-ray diffraction data collection mapping on a given curved surface has been developed. The method uses a laser–video autozalignment system to collect a map of heights on a fine mesh grid. It also reconstructs the surface geometry and determines surface normals on a three-dimensional surface fit. An algorithm was used to calculate the required rotation and tilt angles to coincide the sample normal with the diffraction center before each exposure. A set of diffraction frames collected using this technique was analyzed in order to superimpose a phase transformation map onto a typical zirconia ceramic sample. With this method, mapping of phase and stresses on a complex surface has been demonstrated. This approach is broadly applicable to many important areas of study.


Author(s):  
Ch. Hirsch ◽  
S. Kang ◽  
G. Pointel

The three-dimensional flow in centrifugal impellers is investigated on the basis of a detailed analysis of the results of numerical simulations. In order to gain confidence in this process, an in-depth validation is performed, based on computations of Krain’s centrifugal compressor and of a radial pump impeller, both with vaneless diffusers. Detailed comparisons with available experimental data provide high confidence in the numerical tools and results. The appearance of a high loss ‘wake’ region results from the transport of boundary layer material from the blade surfaces to the shroud region and its location depends on the balance between secondary and tip leakage flows and is not necessarily connected to 3D boundary layer separation. Although the low momentum spots near the shroud can interfere with 3D separated regions, the main outcome of the present analysis is that these are two distinct phenomena. Part I of this paper focuses on the validation base of the numerical approach, based on fine mesh simulations, while Part II presents an analysis of the different contributions to the secondary flows and attempts to estimate their effect on the overall flow pattern.


Author(s):  
Diego Micheli ◽  
Valentino Pediroda ◽  
Stefano Pieri

An automatic approach for the multi-objective shape optimization of microgas turbine heat exchangers is presented. According to the concept of multidisciplinary optimization, the methodology integrates a CAD parametric model of the heat transfer surfaces, a three-dimensional meshing tool, and a CFD solver, all managed by a design optimization platform. The repetitive pattern of the surface geometry has been exploited to reduce the computational domain size, and the constant flux boundary conditions have been imposed to better suit the real operative conditions. A new approach that couples cold and warm fluids in a periodic unitary cell is introduced. The effectiveness of the numerical procedure was verified comparing the numerical results with available literature data. The optimization objectives are maximizing the heat transfer rate and minimizing both friction factor and heat transfer surface. The paper presents the results of the optimization of a 50kWMGT recuperator. The design procedure can be effectively extended and applied to any industrial heat exchanger application.


2017 ◽  
Vol 14 (130) ◽  
pp. 20170202 ◽  
Author(s):  
Joseph Libby ◽  
Arsalan Marghoub ◽  
David Johnson ◽  
Roman H. Khonsari ◽  
Michael J. Fagan ◽  
...  

During the first year of life, the brain grows rapidly and the neurocranium increases to about 65% of its adult size. Our understanding of the relationship between the biomechanical forces, especially from the growing brain, the craniofacial soft tissue structures and the individual bone plates of the skull vault is still limited. This basic knowledge could help in the future planning of craniofacial surgical operations. The aim of this study was to develop a validated computational model of skull growth, based on the finite-element (FE) method, to help understand the biomechanics of skull growth. To do this, a two-step validation study was carried out. First, an in vitro physical three-dimensional printed model and an in silico FE model were created from the same micro-CT scan of an infant skull and loaded with forces from the growing brain from zero to two months of age. The results from the in vitro model validated the FE model before it was further developed to expand from 0 to 12 months of age. This second FE model was compared directly with in vivo clinical CT scans of infants without craniofacial conditions ( n = 56). The various models were compared in terms of predicted skull width, length and circumference, while the overall shape was quantified using three-dimensional distance plots. Statistical analysis yielded no significant differences between the male skull models. All size measurements from the FE model versus the in vitro physical model were within 5%, with one exception showing a 7.6% difference. The FE model and in vivo data also correlated well, with the largest percentage difference in size being 8.3%. Overall, the FE model results matched well with both the in vitro and in vivo data. With further development and model refinement, this modelling method could be used to assist in preoperative planning of craniofacial surgery procedures and could help to reduce reoperation rates.


Sign in / Sign up

Export Citation Format

Share Document