Predicting Finite Element Submodel Boundary Conditions for Contact Models Using Richardson Extrapolation

Author(s):  
Michael W. Sracic ◽  
William J. Elke

Abstract This paper considers an efficient way to apply submodeling methods to finite element models using Richardson Extrapolation. A problem is considered where a rigid cylindrical indenter contacts an elastic half plane (RCEHP). A submodeling method is introduced where the errors of the displacements on the boundaries of the submodel are controlled by employing a best-fit Richardson Extrapolation curve. Specifically, the curve is fit to the convergence relationship of various estimates of submodel boundary displacements. The method is tested on the RCEHP problem, and the results of the model predictions for maximum contact pressure are compared to an analytical and converged global model result. The submodeling method predicted the maximum contact pressure of the RCEHP contact interface to be about 7% higher than the analytical prediction and 5% higher than the converged global model prediction. The error is likely due to the selection of the global and submodel domains, the numerical algorithm used to estimate the Richardson Extrapolation Curve Fits, and the mesh refinements used for the various models. The proposed method solved in about 42.6 minutes while the converged global model solved in 11.19 hours. Future work will aim to provide best practices to reduce error and maximize computational time savings when using the method.

Author(s):  
William J. Elke ◽  
Michael W. Sracic

Abstract Submodeling is a method used in Finite Element Modeling in order to evaluate features of interest, such as fillets or contact interfaces at a reduced cost. Submodeling can be done by creating a full-structure, coarse-mesh, “global” model and solving it. Once the solutions for this model are solved, a model just of the feature of interest, or “submodel”, can be solved using boundary conditions estimated based on global model results. The location of the submodel boundary has large effects on the accuracy of the solution and has been examined by the authors previously, but not for models with contact interfaces. This work uses the submodeling procedure on a model with two cantilever beams sandwiched together with a bolt close to the free end (dubbed “sandwich beam”). Numerous models are produced with different submodel boundary locations in order to better understand how those locations affect the solutions of a model with contact interfaces. The maximum contact pressure was the main metric used to examine the effects and the feature of interest was the bolt. In previous works, it was determined that the global model mesh size and submodel boundary locations are the main sources of error in submodeling [15]. From this work, it was concluded that the inaccuracies of the global model mesh size are magnified when the submodel mesh is too close to the feature of interest. Furthermore, contact pressure tends to be overestimated as the submodel is refined, but it tends to converge as the global model is converged. This work also demonstrated that errors in global model solution (due to meshing) are mitigated when the entire feature of interest is included in the submodel and the submodel boundary locations are far enough away from that feature.


2012 ◽  
Vol 268-270 ◽  
pp. 737-740
Author(s):  
Yang Yu ◽  
Yi Hua Dou ◽  
Fu Xiang Zhang ◽  
Xiang Tong Yang

It is necessary to know the connecting and sealing ability of premium connection for appropriate choices of different working conditions. By finite element method, the finite element model of premium connection is established and the stresses of seal section, shoulder zone and thread surface of tubing by axial tensile loads are analyzed. The results show that shoulder zone is subject to most axial stresses at made-up state, which will make distribution of stresses on thread reasonable. With the increase of axial tensile loads, stresses of thread on both ends increase and on seal section and shoulder zone slightly change. The maximum stress on some thread exceed the yield limit of material when axial tensile loads exceed 400KN. Limited axial tensile loads sharply influence the contact pressures on shoulder zone while slightly on seal section. Although the maximum contact pressure on shoulder zone drop to 0 when the axial tensile load is 600KN, the maximum contact pressure on seal section will keep on a high level.


Author(s):  
Chris Alexander ◽  
Wade Armer ◽  
Stuart Harbert

KOCH Heat Transfer Company contracted Stress Engineering Services, Inc. to perform a design/parameter study of a return bonnet used in hairpin heat exchangers that employs an elliptical flange design. The return bonnet is an important component of the heat exchanger as it can be removed to permit inspection of the heat exchanger tubes. The return bonnet is bolted to the hairpin leg flange. To maintain sealing integrity a gasket is placed between the return bonnet flange and the hairpin leg flange. The sealing efficiency of two return bonnet sizes (24-inch and 30-inch) was investigated in this study using finite element analysis. The sealing efficiency is an indication of how the contact pressure changes circumferentially around the gasket and is calculated by dividing the local contact pressure by the maximum contact pressure calculated in the gasket for each respective design. The study assessed the effects of geometric changes to the mating flanges. Using an iterative design process using finite element analysis, the elliptical flanges were optimized to maximize sealing efficiency. Upon completion of the study, the manufacturer successfully employed the modifications as evidenced with multiple successful hydrotests.


Author(s):  
Louay S. Yousuf ◽  
Yaakob K. H. Dabool

Abstract The bending deflection of the disc cam profile and the dynamic response of the follower were discussed and analyzed for three paths of contact. The objective of this paper was to study the influence of maximum contact pressure on the bending deflection of the cam profile. Numerical simulation was carried out using SolidWorks Software to simulate the follower displacement, velocity and acceleration. Finite element analysis was used taking into account the use of ANSYS package to calculate the bending deflection. The experiment setup had been done through an infrared camera device. The bending deflection of point (18) is bigger than the bending deflection of point (4) because of the bigness of radius of curvature of nose (2).


This research focuses on wear of disc brake by method of analytical and finite element. An SUV car of DD6470C disk brake selected. Factors which considered as input parameters, such as the maximum allowable speed of the car as a constant throughout the work and by varying applied brake pressure. Analytical of the distribution of contact temperature along in the radial direction of disc brake caused by applied heat flux solved by using the separation of the variable method. Finite element simulation for contact pressure and von misses stress using; ANSYS workbench for the case of the structural analysis of disc brake was done by applying brake pressure and angular velocity. But, the thermal-structural, the maximum -contact temperature value of disc considered in addition. The contact pressure and von misses stresses were calculated analytically and ANSYS workbench results were presented in contour plot and numerically. The result shows that contact pressure, Von Misses stress and wear increased by increasing brake pressure in the case of the structure. And also parameters increase as increasing of both brake pressure and contact temperature of the disc in case of thermal-structural in all aspects like contact pressure, Von Misses and wear.


Author(s):  
Weiguo Ma ◽  
Baolong Qu ◽  
Feng Guan

A systematic study of the packer rubber contact pressure under a fixed-displacement load is conducted to gain further insight into the packer seal mechanism. A Y221-114 double rubber packer is investigated using the finite element software ANSYS, where a design of experiments method is utilized to study the effects of the friction coefficient. The results show that the friction coefficient of the packer and the tubing had the greatest effect on contact pressure than other factors. Decreasing the rubber friction coefficient is conducive to forming the double rubber seal and increasing the maximum contact pressure working range. However, there is additionally a slight decrease in the value of maximum contact pressure. The results of the study provide valuable insight into the importance of packer design optimization.


2018 ◽  
Vol 165 ◽  
pp. 22029 ◽  
Author(s):  
Dong-Hyung Lee ◽  
Ha-Young Choi ◽  
Seok-Jin Kwon ◽  
Jeong-Won Seo

In the shrink or press-fitted shafts such as railway axles, the rotor of a steam turbine or coupling, a high-stress concentration takes place in the close of contact edge due to relative slip between shaft and boss in a press-fitted shaft and this is a major cause of fatigue failure of the shaft. The object of this paper is to build a finite element analysis model for analysing press-fitted and bending load condition in a pressfitted assembly and is to propose a hub shape optimization method to minimize a contact pressure in the close of shaft contact edge. Numerical asymmetric-axisymmetric finite element model was developed to predict the contact stress state of the press-fitted shaft. Global optimization method, genetic algorithm, and local optimization method, sequential quadratic programming, was applied to the press-fitted assembly to optimize the hub contact edge geometry. The results showed that the maximum contact pressure with optimized hub shape decreased more than 60% compared to conventional hub shape, the maximum contact stress affecting fatigue life reduced about 47%. In addition, hub shape optimization design could be a useful tool, able to increase the load capabilities of press fits concerning wear and fatigue behaviour.


Author(s):  
Ehsan Askari ◽  
Michael S Andersen

Using Hertz contact law results in inaccurate outcomes when applied to the soft conformal hip implants. The finite element method also involves huge computational time and power. In addition, the sliding distance computed using the Euler rotation method does not incorporate tribology of bearing surfaces, contact mechanics and inertia forces. This study, therefore, aimed to develop a nonlinear dynamic model based on the multibody dynamic methodology to predict contact pressure and sliding distance of metal-on-polyethylene hip prosthesis, simultaneously, under normal walking condition. A closed-form formulation of the contact stresses distributed over the articulating surfaces was derived based upon the elastic foundation model, which reduced computational time and cost significantly. Three-dimensional physiological loading and motions, inertia forces due to hip motion and energy loss during contact were incorporated to obtain contact properties and sliding distance. Comparing the outcomes with that available in the literature and a finite element analysis allowed for the validation of our approach. Contours of contact stresses and accumulated sliding distances at different instants of the walking gait cycle were investigated and discussed. It was shown that the contact point at each instant was located within the zone with the corresponding highest accumulated sliding distance. In addition, the maximum contact pressure and area took place at the stance phase with a single support. The stress distribution onto the cup surface also conformed to the contact point trajectory and the physiological loading.


2005 ◽  
Vol 128 (1) ◽  
pp. 115-123 ◽  
Author(s):  
Barbara Zielinska ◽  
Tammy L. Haut Donahue

The goal of this study is to quantify changes in knee joint contact behavior following varying degrees of the medial partial meniscectomy. A previously validated 3D finite element model was used to simulate 11 different meniscectomies. The accompanying changes in the contact pressure on the superior surface of the menisci and tibial plateau were quantified as was the axial strain in the menisci and articular cartilage. The percentage of medial meniscus removed was linearly correlated with maximum contact pressure, mean contact pressure, and contact area. The lateral hemi-joint was minimally affected by the simulated medial meniscectomies. The location of maximum strain and location of maximum contact pressure did not change with varying degrees of partial medial meniscectomy. When 60% of the medial meniscus was removed, contact pressures increased 65% on the remaining medial meniscus and 55% on the medial tibial plateau. These data will be helpful for assessing potential complications with the surgical treatment of meniscal tears. Additionally, these data provide insight into the role of mechanical loading in the etiology of post-meniscectomy osteoarthritis.


Sign in / Sign up

Export Citation Format

Share Document