A Simple Method of Handling Thermomechanical Coupling for Temperature Computation in a Rolling Tire

2004 ◽  
Vol 32 (2) ◽  
pp. 56-68 ◽  
Author(s):  
S. Futamura ◽  
A. Goldstein

Abstract The thermomechanical analysis of a pneumatic tire is a highly complex process due to the effects of temperature on both the mechanical state and the viscoelastic energy dissipation in the tire. This coupled thermomechanical behavior typically requires that rolling tire temperatures be determined iteratively. As a result, a steady-state analysis involves updating the temperature dependent elastic and viscoelastic properties as the solution proceeds. The process is further complicated in a non-steady-state analysis where material properties need to be updated at multiple intervals in time. A simplified method is proposed. First, the sensitivity of the tire elastic response of the tire to changes in material stiffness is characterized using the “deformation index”. Then, using a commercial finite element program and an appropriate user subroutine, heat generation is expressed as a function of the local temperature using a simple algebraic expression involving the temperature dependent material properties and the deformation indices. Temperatures are computed using the finite element program with the coupling information contained in the user subroutine. The result is a simplified method for a fully coupled thermomechanical analysis of a tire for steady-state and transient thermal analysis. The accuracy and the simplicity of the method are demonstrated using a small “tire-like” model. The simplified method is compared to the fully coupled iterative method for a steady-state thermal solution.

2013 ◽  
Vol 457-458 ◽  
pp. 354-357
Author(s):  
Yu Jie Sun ◽  
Qing Chun Cui ◽  
Suo Huai Zhang ◽  
Li Jun Yan

The objective of this paper provides a numerical implementation procedure of thermo-metallurgical-mechanical constitute equation based on additively decomposition of strain rate. Together with phase transformation kinetics, the macro material properties are determined by assigning temperature dependent material properties to each phase and by applying mixture rule to combine. Then the constitute equation is implemented into general purpose implicit finite element program via user material subroutine. The effectiveness of developed computational method is confirmed by a Satoh test simulation. Simulation of Satoh test demonstrates that transformation induce plasticity has significant effect of the evolution of residual stress and can not be neglected for alloy steel during hot working process.


2013 ◽  
Vol 798-799 ◽  
pp. 235-238
Author(s):  
Xiao Hui Chen ◽  
Xu Chen

The current paper reports the results of a numerical simulation and experiment of ratcheting behavior of pressurized straight pipe under reversed bending. A nonlinear isotropic/ kinematic (combined) hardening model is implemented into finite element program ANSYS by writing own user subroutine in FORTRAN language. The results of the numerical simulation is compared with experimental data. A reasonable agreement is noticed between the experimental and the numerical results for the ratcheting behavior of the pressurized straight pipe subjected to reversed bending.


2013 ◽  
Vol 785-786 ◽  
pp. 16-19 ◽  
Author(s):  
Xiao Hui Chen ◽  
Xu Chen

The paper compares numerical simulation with experimental results of pressurized elbow piping subjected to reversed in-plane bending in elastoplastic domain. The modified AbdelKarim-Ohno model is implemented into finite element program ANSYS by writing own user subroutine in FORTRAN language. The modified AbdelKarim-Ohno model may improve the prediction accuracy of ratcheting behavior of pressurized elbow under cyclic loading.


2007 ◽  
Vol 353-358 ◽  
pp. 1013-1016
Author(s):  
Jian Bing Sang ◽  
Su Fang Xing ◽  
Xiao Lei Li ◽  
Jie Zhang

It has been well known that rubber-like material can undergo large deformation and exhibit large nonlinear elastic behavior. Because of the geometrically nonlinear of rubber like material, it is more difficult to analyze it with finite element near the notch tip. What is more, because there are varieties of the strain energy functions, implementation of these models in a general finite element program to meet the need of industry applications can be time consuming. In order to make use of the constitutive equation of Y.C. Gao in 1997 and analyze the notch tip of rubber-like material, a framework to implement the rubber-like material model is established within the general-purpose finite element program MSC.Marc. It will be very convenient to implement this isotropic hyperelastic model into the program with a user subroutine. This paper starts with the theoretical analysis based on the strain energy function given by Y.C. Gao in 1997. A user subroutine is programmed to implement this strain energy function into the program of MSC.Marc, which offer a convenient method to analyze the stress and strain of rubber-like material with the strain energy function that is needed. Though analysis with MSC.Marc, it is found that the result with finite element is consistent with the analytical result that given by Y.C. Gao in 1997, which testify that analyzing rubber like material with this method is reasonable and convenient.


1993 ◽  
Vol 115 (2) ◽  
pp. 147-152 ◽  
Author(s):  
Y.-H. Pao ◽  
S. Badgley ◽  
E. Jih ◽  
R. Govila ◽  
J. Browning

The thermal cyclic shear stress/strain hysteresis response and associated steady-state creep parameters of 97Sn-3Cu solder joints have been determined using a beam specimen previously developed by Pao et al. (1992a). The solder joint was subjected to a 40-minute thermal cycling from 40°C to 140°C. A constitutive equation based on elastic and steady-state creep deformation for the solder has been formulated and implemented in a finite element program, ABAQUS, to model the experiment. The results show that the constitutive equation based on one single creep mechanism cannot fully account for the deformation during cooling, as opposed to the case of 90Pb-10Sn where the entire cyclic deformation can be well modeled by a similar constitutive equation (Pao et al., 1992c). This suggests that another creep mechanism is dominant for lower stresses and higher temperature. The thermal fatigue results show that the failure mechanism of 97Sn-3Cu joints is similar to that of 90Pb-10Sn joints, but the number of cycles to failure of 97Sn-Cu solder joints is at least 5 times longer than 90Pb-10Sn solder joints. This indicates the potential application of 97Sn-3Cu in place of 90Pb-10Sn solder.


2004 ◽  
Vol 120 ◽  
pp. 727-735
Author(s):  
F. Frerichs ◽  
Th. Lübben ◽  
U. Fritsching ◽  
H. Lohner ◽  
A. Rocha ◽  
...  

The prediction of mechanical behaviour of specimen during heat treatment by means of numerical simulation requires numerous modules e.g. for heat transfer and mechanical behaviour. The quality of predictions depend on the quality of the applied models within the modules. In this paper the strain hardening model used in the mechanical module will be investigated. For simulation of mechanical behaviour during gas quenching it is first of all necessary to calculate the interaction between gas and specimen. Using simulated flow field and temperature distribution within the gas, the heat transfer coefficient is calculated from computational fluid dynamics. The cooling and further the mechanical behaviour e.g. residual stresses and distortion of the specimen are simulated by a commercial Finite Element program. To investigate strain hardening it is helpful to choose in a first step a material that will not show phase transformations due to heat treatment. Therefore simulation of mechanical behaviour of austenitic cylinders (SAE30300) is investigated. The required thermo-physical properties such as thermal conductivity, density, and specific heat are taken from literature. With the exception of Poisson’s ratio the mechanical properties are measured and calculated by own investigations. For description of the temperature dependent stress strain curves the Ramberg-Osgood model is used. The simulated results are compared with experimental data in order to decide which model better describes the mechanical response, whether the kinematic or isotropic strain hardening.


2004 ◽  
Vol 467-470 ◽  
pp. 733-738 ◽  
Author(s):  
Jason Gruber ◽  
Denise C. George ◽  
Andrew P. Kuprat ◽  
Gregory S. Rohrer ◽  
Anthony D. Rollett

Through simulations with the moving finite element program GRAIN3D, we have studied the effect of anisotropic grain boundary energy on the distribution of boundary types in a polycrystal during normal grain growth. An energy function similar to that hypothesized for magnesia was used, and the simulated grain boundary distributions were found to agree well with measured distributions. The simulated results suggest that initially random microstructures develop nearly steady state grain boundary distributions that have local maxima and minima corresponding to local minima and maxima, respectively, of the energy function.


Author(s):  
Erik Tamsen ◽  
Daniel Balzani

AbstractIn this paper we present a fully-coupled, two-scale homogenization method for dynamic loading in the spirit of FE$$^2$$ 2 methods. The framework considers the balance of linear momentum including inertia at the microscale to capture possible dynamic effects arising from micro heterogeneities. A finite-strain formulation is adapted to account for geometrical nonlinearities enabling the study of e.g. plasticity or fiber pullout, which may be associated with large deformations. A consistent kinematic scale link is established as displacement constraint on the whole representative volume element. The consistent macroscopic material tangent moduli are derived including micro inertia in closed form. These can easily be calculated with a loop over all microscopic finite elements, only applying existing assembly and solving procedures. Thus, making it suitable for standard finite element program architectures. Numerical examples of a layered periodic material are presented and compared to direct numerical simulations to demonstrate the capability of the proposed framework. In addition, a simulation of a split Hopkinson tension test showcases the applicability of the framework to engineering problems.


1991 ◽  
Vol 18 (4) ◽  
pp. 537-555 ◽  
Author(s):  
Diana E. Chernenko ◽  
D. J. Laurie Kennedy

The Canadian Standards Association Standard CAN-CSA-S16.1-M89 “Limit states design of steel structures” assigns welded wide flange (WWF) columns to the column curve for rolled H-shaped sections. This is conservative because of differences in the production of WWF and rolled sections. The WWFs are stipulated to have flame-cut edges. The residual stress pattern, with favourable tensile stresses, results in a delayed loss of stiffness as weak axis inelastic buckling occurs. Thus the weak axis and strong axis buckling curves lie closer together for WWF shapes than they do for rolled H-shapes. Automatic cutting and welding processes result in close tolerances on out-of-straightness. As well, the statistical variations in the geometric properties are favourable. A detailed analysis of data collected from mill records and on-site measurements was made to obtain statistical parameters of relevant geometric and material properties. A finite element program modelling inelastic behaviour, residual strain patterns, out-of-straightness, and material properties was used with the test results of others to establish test-to-predicted ratios of column strengths. Parametric studies provided an assessment of the effect of varying residual strain patterns and column out-of-straightness. This formed the basis for determining the factored compressive resistance of WWF sections for three different slenderness ratios. Key words: columns, inelastic analysis, out-of-straightness, residual stresses, welded wide flange.


Sign in / Sign up

Export Citation Format

Share Document