Simulation of ratcheting strain to a high number of cycles under biaxial loading

The Ohno–Wang kinematic hardening rule is modified to incorporate the Burlet–Cailletaud radial evanescence term for an improved simulation of the ratcheting behavior. The Delobelle parameter δ^′ is implemented in the modified model to compromise shakedown of the Burlet–Cailletaud hardening rule and over-prediction of the Ohno–Wang model. An evolution equation is proposed for δ^′ to simulate the ratcheting strain over an extended domain of cycles. Ratcheting tests were conducted on S45C steel under four types of nonproportional axial–torsional loading. The new model is found to yield reasonably accurate predictions of ratcheting strain to a much higher number of cycles compared with other studies.     

Non-linear dynamics of a drill-string with uncertain model of the bit–rock interaction

The drill-string dynamics is difficult to predict due to the non-linearities and uncertainties involved in the problem. In this paper a stochastic computational model is proposed to model uncertainties in the bit–rock interaction model. To do so, a new strategy that uses the non-parametric probabilistic approach is developed to take into account model uncertainties in the bit–rock non-linear interaction model. The mean model considers the main forces applied to the column such as the bit–rock interaction, the fluid–structure interaction and the impact forces. The non-linear Timoshenko beam theory is used and the non-linear dynamical equations are discretized by means of the finite element method.     

Drag models for simulating gas–solid flow in the turbulent fluidization of FCC particles

This paper examines the suitability of various drag models for predicting the hydrodynamics of the turbulent fluidization of FCC particles on the Fluent V6.2 platform. The drag models included those of Syamlal–O’Brien, Gidaspow, modified Syamlal–O’Brien, and McKeen. Comparison between experimental data and simulated results showed that the Syamlal–O’Brien, Gidaspow, and modified Syamlal–O’Brien drag models highly overestimated gas–solid momentum exchange and could not predict the formation of dense phase in the fluidized bed, while the McKeen drag model could not capture the dilute characteristics due to underestimation of drag force. The standard Gidaspow drag model was then modified by adopting the effective particle cluster diameter to account for particle clusters, which was, however, proved inapplicable for FCC particle turbulent fluidization. A four-zone drag model (dense phase, sub-dense phase, sub-dilute phase and dilute phase) was finally proposed to calculate the gas–solid exchange coefficient in the turbulent fluidization of FCC particles, and was validated by satisfactory agreement between prediction and experiment.     

Effect of streamline curvature on flow field of a turbulent plane jet in cross-flow

The effects of the streamline curvature and finite edge velocity on the flow field of a turbulent plane jet in cross-flow are studied numerically by incorporating the curvature effect in the k–ε turbulence model. The improvement in the predictions by the streamline curvature model is assessed by comparing its prediction with those by the standard k–ε model. The predictions by both the models are compared with available experimental data. It has been observed that the performance of the k–ε model with streamline curvature modification is superior to the standard k–ε model.     

Effective elastic moduli of two dimensional solids with distributed cohesive microcracks

Effective elastic properties of a defected solid with distributed cohesive micro-cracks are estimated based on homogenization of the Dugdale–Bilby–Cottrell–Swinden (Dugdale–BCS) type micro-cracks in a two dimensional elastic representative volume element (RVE).Since the cohesive micro-crack model mimics various realistic bond forces at micro-scale, a statistical average of cohesive defects can effectively represent the overall properties of the material due to bond breaking or crack surface separation in small scale. The newly proposed model is distinctive in the fact that the resulting effective moduli are found to be pressure sensitive.     

A thermoviscoelastic model for amorphous shape memory polymers: Incorporating structural and stress relaxation

A thermoviscoelastic constitutive model is developed for amorphous shape memory polymers (SMP) based on the hypothesis that structural and stress relaxation are the primary molecular mechanisms of the shape memory effect and its time-dependence. This work represents a new and fundamentally different approach to modeling amorphous SMPs. A principal feature of the constitutive model is the incorporation of the nonlinear Adam–Gibbs model of structural relaxation and a modified Eyring model of viscous flow into a continuum finite–deformation thermoviscoelastic framework. Comparisons with experiments show that the model can reproduce the strain–temperature response, the temperature and strain-rate dependent stress–strain response, and important features of the temperature dependence of the shape memory response. Because the model includes structural relaxation, the shape memory response also exhibits a dependence on the cooling and heating rates.     

Nonequilibrium stretching dynamics of dilute and entangled linear polymers in extensional flow

We propose an extension of the FENE-CR model for dilute polymer solutions [M.D. Chilcott, J.M. Rallison, Creeping flow of dilute polymer solutions past cylinders and spheres, J. Non-Newtonian Fluid Mech. 29 (1988) 382–432] and the Rouse-CCR tube model for linear entangled polymers [A.E. Likhtman, R.S. Graham, Simple constitutive equation for linear polymer melts derived from molecular theory: Rolie–Poly equation, J. Non-Newtonian Fluid Mech. 114 (2003) 1–12], to describe the nonequilibrium stretching dynamics of polymer chains in strong extensional flows. The resulting models, designed to capture the progressive changes in the average internal structure (kinked state) of the polymer chain, include an ‘effective’ maximum contour length that depends on local flow dynamics. The rheological behavior of the modified models is compared with various results already published in the literature for entangled polystyrene solutions, and for the Kramers chain model (dilute polymer solutions). It is shown that the FENE-CR model with an ‘effective’ maximum contour length is able to describe correctly the hysteretic behavior in stress versus birefringence in start-up of uniaxial extensional flow and subsequent relaxation also observed and computed by Doyle et al. [P.S. Doyle, E.S.G. Shaqfeh, G.H. McKinley, S.H. Spiegelberg, Relaxation of dilute polymer solutions following extensional flow, J. Non-Newtonian Fluid Mech. 76 (1998) 79–110] and Li and Larson [L. Li, R.G. Larson, Excluded volume effects on the birefringence and stress of dilute polymer solutions in extensional flow, Rheol. Acta 39 (2000) 419–427] using Brownian dynamics simulations of bead–spring model. The Rolie–Poly model with an ‘effective’ maximum contour length exhibits a less pronounced hysteretic behavior in stress versus birefringence in start-up of uniaxial extensional flow and subsequent relaxation.     

Nonlinear analysis of beam–column joints using softened truss model

The aim of this study is to expand the application of the nonlinear softened truss model for membrane elements on beam–column joints. The softened truss model employs three equations for equilibrium, three for compatibility and four equations for the constitutive laws of materials. The constitutive equations for both the concrete and steel are based on the actually observed stress–strain relationships. The model has three important attributes. The first is the nonlinear association of stress and strain. The second, and conceivably more noteworthy, is the softening of concrete in compression due to tensile strains in the perpendicular direction. The third is that the influence of the concrete tensile stresses between cracks on the average stress–strain relationship for reinforcing steel and the influence of orthogonal tensile stresses on the compression stress–strain relationship for concrete can be considered in the model. For beam–column joints, one of the most important factors influencing the behaviour is certainly the bond conditions of the beam bars. In this study, the softened truss model is expanded to take into account the influence of this important factor into account. In the revised version of the model, full strain compatibility does not exist between the steel reinforcement and the surrounding concrete and thus the factors influencing the bond-slip between concrete and reinforcement is adequately considered. The improved softened truss model is applied on 51 exterior beam–column joint tests. It is apparent from the results that the revised model gives very accurate predictions of the shear strength of joints and is an improvement on the existing version of the model proposed by Hsu.     

Formation of Nonmetallic Inclusions of Variable Composition During Unidirectional Solidification of Liquid Steel

A two–dimensional unsteady mathematical model is presented for numerical study of the formation of nitrogen–containing compounds of variable composition during unidirectional solidification of liquid steel. The model allows calculation of the distributions of temperature and concentrations of soluble admixtures, the shape of solidification boundaries, the coordinate of the beginning of the chemical reaction, and the composition and mass of inclusions formed in the uncrystallized zone of the ingot and in the two–phase region.     

Strain-induced passivity breakdown in corrosion crack initiation

Corrosion crack nucleation and growth are modelled as a moving boundary problem. The model incorporates three physical processes––dissolution, passivation and straining––into a continuum mechanical framework. The dissolution triggers surface advance; the passivation restrains the access of the environment to bare metal; the deformation causes for passivity breakdown. Plane cracks nucleating from surface pits in an elastic–plastic material body under fatigue load are considered. The problem is solved using a FEM program and a moving boundary tracking procedure. The model simulates how cracks form and grow in a single continuous course. The geometry of the developed cracks is found independent of the initial pit size. Plasticity is found to influence the curvature at the tip of the nucleated corrosion cracks. The most important evolution length parameter, the width of the corrosion crack, is found to depend on the size constraints of the tracking procedure. It is concluded that the model is deficient for determining all length scales observed in reality. Physical processes to be considered in an advanced model are proposed and discussed.     

Modeling stress state dependent damage evolution in a cast Al–Si–Mg aluminum alloy

Internal state variable rate equations are cast in a continuum framework to model void nucleation, growth, and coalescence in a cast Al–Si–Mg aluminum alloy. The kinematics and constitutive relations for damage resulting from void nucleation, growth, and coalescence are discussed. Because damage evolution is intimately coupled with the stress state, internal state variable hardening rate equations are developed to distinguish between compression, tension, and torsion straining conditions. The scalar isotropic hardening equation and second rank tensorial kinematic hardening equation from the Bammann–Chiesa–Johnson (BCJ) Plasticity model are modified to account for hardening rate differences under tension, compression, and torsion. A method for determining the material constants for the plasticity and damage equations is presented. Parameter determination for the proposed phenomenological nucleation rate equation, motivated from fracture mechanics and microscale physical observations, involves counting nucleation sites as a function of strain from optical micrographs. Although different void growth models can be included, the McClintock void growth model is used in this study. A coalescence model is also introduced. The damage framework is then evaluated with respect to experimental tensile data of notched Al–Si–Mg cast aluminum alloy specimens. Finite element results employing the damage framework are shown to illustrate its usefulness.     

Effect of Evaporation of Liquid Droplets on the Distribution of Parameters in a Two–Species Laminar Flow

A calculation model was developed, and the heat– and mass–transfer characteristics in a laminar air—vapor—droplet flow moving in a round tube were studied numerically. The distributions of parameters of the two–phase flow over the tube radius were obtained for varied initial concentrations of the gas phase. The calculated heat and mass transfer is compared to experimental data and calculations of other authors. It is shown that evaporation of droplets in a vapor—gas flow leads to a more intense heat release as compared to a one–species vapor—droplet flow and one–phase vapor flow     

A penny-shaped cohesive crack model for material damage

A novel micromechanics based damage model is proposed to address failure mechanism of defected solids with randomly distributed penny-shaped cohesive micro-cracks (Barenblatt–Dugdale type). Energy release contribution to the material damage process is estimated in a representative volume element (RVE) under macro hydrostatic stress state. Macro-constitutive relations of RVE are derived via self-consistent homogenization scheme, and they are characterized by effective nonlinear elastic properties and a class of pressure sensitive plasticity which depends on crack opening volume fraction and Poisson’s ratio. Several distinguished features of the present model are compared with Gurson model and Gurson–Tvergaard–Needleman (GTN) model, showing that the proposed model can better capture material degradation and catastrophic failure due to cohesive micro-crack growth and coalescence.     

Numerical Study of Hydrodynamics and Heat and Mass Transfer of a Ducted Gas–Vapor‐Droplet Flow

A computational model has been developed to predict heat and mass transfer and hydrodynamic characteristics of a turbulent gas–vapor–droplet flow. Turbulent characteristics of the gas phase are computed using the k– model of turbulence. It is shown that, with increasing inlet droplet diameter, the rate of heat transfer between the duct surface and the vapor–gas mixture decreases appreciably, whereas the wall friction increases only insignificantly. The predicted values agree fairly well with available experimental and numerical data     

Identification of the Ten Inertia Parameters of a Rigid Body

Standard experiments for identifying inertia parameters of a rigid body only provide a subset of the inertia parameters of the body [1–10]. In addition, they do not use in the estimation process the complete information included in the equations of motion of the rigid test body. The objective of the work described in this paper is the simultaneous, automatic experimental identification of the ten inertia parameters of a rigid body using the complete information hidden in the nonlinear model equations of the test body. This task has been solved in several steps:– mathematical modelling of the special motions of a rigid body in space. These model equations have been mapped into a form suitable for identification purposes (identification hypothesis)– design of a special measurement robot for performing the identification experiments– laboratory experiments providing test data used for the identification experiments– identification of the inertia parameters and accuracy tests.The accuracy of the identified parameters is satisfactory.     

Non-Darcian Effects on the Onset of Convection in a Porous Layer with Throughflow

The Brinkman extended Darcy model including Lapwood and Forchheimer inertia terms with fluid viscosity being different from effective viscosity is employed to investigate the effect of vertical throughflow on thermal convective instabilities in a porous layer. Three different types of boundary conditions (free–free, rigid–rigid and rigid–free) are considered which are either conducting or insulating to temperature perturbations. The Galerkin method is used to calculate the critical Rayleigh numbers for conducting boundaries, while closed form solutions are achieved for insulating boundaries. The relative importance of inertial resistance on convective instabilities is investigated in detail. In the case of rigid–free boundaries, it is found that throughflow is destabilizing depending on the choice of physical parameters and the model used. Further, it is noted that an increase in viscosity ratio delays the onset of convection. Standard results are also obtained as particular cases from the general model presented here.     

Model of the internal circulation of a liquid drop in a gas stream

An analytical model for calculating the internal flow in a drop is proposed. The model is checked by numerically solving the exact boundary-value problem for the Navier-Stokes equations.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 6, pp. 167–170, November–December, 1989.