Prediction of 42CrMo steel flow stress at high temperature and strain rate

The compressive deformation behavior of 42CrMo steel was investigated at temperatures ranging from 850 to 1150 °C and strain rates from 0.01 to 50 s⁻¹ on Gleeble-1500 thermo-simulation machine. Based on the classical stress–dislocation relation and the kinematics of the dynamic recrystallization, the flow stress constitutive equations of the work hardening-dynamical recovery period and dynamical recrystallization period were established for 42CrMo steel, respectively. The stress–strain curves of 42CrMo steel predicted by the established models are in good agreement with experimental results when the strain rate is relatively low. So, the proposed deformation constitutive equations can be used to establish the hot formation processing parameters for 42CrMo steel.     

Flow stress behavior of porous FVS0812 aluminum alloy during hot-compression

The hot deformation behavior of porous FVS0812 aluminum alloy prepared by spray deposition was studied by means of compression tests on a Gleeble 1500 machine. The samples were hot compressed at temperatures ranging from 573 K to 773 K under various true strain rates of 10⁻⁴–10⁰ s⁻¹. The deformation behaviors are characterized by a significant strain hardening during hot-compression due to the progressive compaction of the pores with increasing compressive strain. A revised formula describing the relationships of the flow stress, strain rate and temperature of the porous alloy at elevated temperatures is proposed by compensation of strain. The theoretical predictions are compared with experimental results, which show good agreement.     

Measurement of deformation gradients in hot rolling of AA3004

In this paper we describe an experimental technique developed to measure the deformation gradients and temperature in a single hot rolling pass of an AA3004 sample that was fitted with an insert. The insert had been previously hand engraved with a 1×1 mm grid pitch, and the analysis of the data digitally captured from the image of the deformed grid enabled the calculation of the components of the deformation gradient tensor. Four steel pins prevented relative motion between the insert and the rest of the sample. No detachment was observed between insert and sample after rolling. The temperature was measured during rolling using two embedded thermocouples, one close to the surface and the other in the centerline. The commercial finite element code ABAQUS was used to create a three-dimensional model of the rolling process. The recorded temperature was compared to the numerical values evaluated after tuning the heat transfer coefficient. The shape of the grid after rolling was checked against the deformed mesh using different fricition coefficients in order to obtain the optimum match. The unusually large length of the insert enabled the rolling process to be stopped halfway so that a picture of the roll-gap area could be obtained. This provided a partially deformed grid that represented the transient state during rolling. The experimentaily determined deformation gradient in this area as well as in the steady-state area agreed well with the finite element oredictions.     

Artificial neural network prediction to the hot compressive deformation behavior of Al–Cu–Mg–Ag heat-resistant aluminum alloy

The behavior of the flow stress of Al-Cu-Mg-Ag heat-resistant aluminum alloys during hot compression deformation was studied by thermal simulation test. The temperature and the strain rate during hot compression were 340-500 °C, 0.001 s⁻¹ to 10 s⁻¹, respectively. Constitutive equations and an artificial neural network (ANN) model were developed for the analysis and simulation of the flow behavior of the Al-Cu-Mg-Ag alloys. The inputs of the model are temperature, strain rate and strain. The output of the model is the flow stress. Comparison between constitutive equations and ANN results shows that ANN model has a better prediction power than the constitutive equations.     

Time dependent flow in equimolar micellar solutions: transient behaviour of the shear stress and first normal stress difference in shear induced structures coupled with flow instabilities

Recently we studied time dependent structural changes that are coupled with flow instabilities (Fischer 1998; Wheeler 1998; Fischer 2000). Within a stability analysis, a classification scheme for the feedback circuit of coupled shear-induced structure and flow instabilities was derived by Schmitt et al. (1995) and applied to our samples. Here, inhomogeneous flow layers of different concentration and viscosity are generated by shear-induced diffusion (spinodal demixing) and, as consequence, one no longer observes a homogeneous solution but a type of shear banding that is seen here for the first time. In this paper we present the behaviour of the first normal stress difference observed in the critical shear-rate regime where transient shear-induced structure is coupled with flow instability. Similar to the oscillations of the shear stresses (strain-controlled rheometer) one observes oscillations in the first normal stress difference. This behaviour indicates that elastic structures are built up and destroyed while the shear-induced structures occur and that the induced phase is more elastic than the initial one. Oscillations of shear stress and first normal stress difference are in phase and indicate that both phenomena are caused by the same mechanism. Received: 30 June 1999/Accepted: 14 December 1999     

3D dynamic stress intensity factors at three rectangular cracks in an infinite elastic medium subjected to a time-harmonic stress wave

Summary Dynamic stresses around three coplanar cracks in an infinite elastic medium are determined in the paper. Two of the cracks are equal, rectangular and symmetrically situated on either side of the centrally located rectangular crack. Time-harmonic normal traction acts on each surface of the three cracks. To solve the problem, two kind of solutions are superposed: one is a solution for a rectangular crack in an infinite elastic medium, and the other one is that for two rectangular cracks in an infinite elastic medium. The unknown coefficients in the combined solution are determined by applying the boundary conditions at the surfaces of the cracks. Finally, stress intensity factors are calculated numerically for several crack configurations. Received 14 July 1998; accepted for publication 2 December 1998     

DEFORMATION RIGIDITY OF ASSUMED STRESS MODES IN HYBRID ELEMENTS

The new methods to determine the zero-energy deformation modes in the hybrid elements and the zero-energy stress modes in their assumed stress fields are presented by the natural deformation modes of the elements. And the formula of the additional element deformation rigidity due to additional mode into the assumed stress field is derived. Based on, it is concluded in theory that the zero-energy stress mode cannot suppress the zero-energy deformation modes but increase the extra rigidity to the nonzero-energy deformation modes of the element instead. So they should not be employed to assume the stress field. In addition, the parasitic stress modes will produce the spurious parasitic energy and result the element behaving over rigidity. Thus, they should not be used into the assumed stress field even though they can suppress the zero-energy deformation modes of the element. The numerical examples show the performance of the elements including the zero-energy stress modes or the parasitic stress modes.     

Prediction of residual flow and thermoviscoelastic stresses in injection molding

A general model for predicting the total residual stresses generated during filling and cooling stages of injection-molded parts has been developed. The model takes into account the phenomena associated with non-isothermal stress relaxation. The main hypothesis in our approach is to use the kinematics of a generalized Newtonian fluid at the end of the filling stage as the initial state for the calculation of residual flow stresses. These stresses are calculated using a single integral rheological model (Wagner model). The calculation of stresses developed during the cooling stage is based on a thermoviscoelastic model with structural relaxation. Illustrative results emphasizing the effect of both the melt temperature and the flow rate during the filling stage are presented.     

Squeeze flow of Bingham fluids under slip with friction boundary condition

This paper develops a theoretical analysis of a Bingham fluid in slipping squeeze flow. The flow field decomposition consists in combining a central extensional flow zone in the plane of symmetry and shear flow zones near the plates. It is also considered that the slipping zone is located around a central sticking zone as previously shown from experiments. It is assumed that the shear stress at the plates is constant in the slipping zone and equals a fixed friction yield value. The squeeze force required to compress a Bingham fluid under the slipping behaviour as well as the radial evolution of the transition point between both sticking and slipping zones are finally determined.     

Experimental study of the induced flow birefringence of a suspension of rigid particles at the transition from Couette flow to Taylor vortex flow

The flow birefringence induced in solutions of rigid particles is studied experimentally in the region of the axisymmetrical Taylor vortex flow which arises once the velocity gradient G in the annular gap of a conventional Couette cell reaches a critical value G _c .The measurements are performed for several values of G > G _c and for 10 radial observation points in the annular gap. Solutions of two types of rigid particles are investigated: the first is a suspension of flattened clay particles like bentonite, while the second contains rod-like particles of tobacco mosaic virus (TMV). The variations of the birefringence intensity n and of the extinction angle measured in the domain of the axisymmetrical flow show a different behavior according to the shape of the particle in solution. This fact is confirmed theoretically with a good agreement for the measurements performed with solutions of flat particles.     

Prediction of incompressible flow separation with the approximate factorization technique

This paper describes one application of the approximate factorization technique to the solution of incompressible steady viscous flow problems in two dimensions. The velocity-pressure formulation of the Navier-Stokes equations written in curvilinear non-orthogonal co-ordinates is adopted. The continuity equation is replaced with one equation for the pressure by means of the artificial compressibility concept to obtain a system parabolic in time. The resulting equations are discretized in space with centred finite differences, and the steady state solution obtained by a time-marching ADI method requiring to solve 3 x 3 block tridiagonal linear systems. An optimized fourth-order artificial dissipation is introduced to damp the numerical instabilities of the artificial compressibility equation and ensure convergence. The resulting solver is applied to the prediction of a wide variety of internal flows, including both streamlined boundaries and sharp corners, and fast convergence and good results obtained for all the configurations investigated.     

Investigation of scale effect for the computation of turbulent flow around a circular cylinder

In order to investigate the scale effect of turbulent flow around a circular cylinder, two similarity numbers （criteria） based on turbulent kinetic and dissipation rates associ- ated with the fluctuation characteristics of turbulence wake are deduced by analyzing the Reynolds averaged NavierStokes equations （RANS）. The RNG k-s models and finite volume method are used to solve the governing equations and the second-order implicit time and upwind space discretization algorithms are used to discrete the governing equations. A numerical computation of flow parameters around a two-dimensional circular cylinder with Reynolds numbers ranging from 102 to l07 is accomplished and the result indicates that the fluctuation of turbulence flow along the center line in the wake of circular cylinder can never be changed with increasing Reynolds numbers when Re ≥ 3 × 10^6. This conclusion is useful for controlling the scale of numerical calculations and for applying model test data to engineering practice.     

Orthogonal basic deformation mode method for zero-energy mode suppression of hybrid stress elements

A set of basic deformation modes for hybrid stress finite elements are directly derived from the element displacement field. Subsequently, by employing the so-called united orthogonal conditions, a new orthogonalization method is proposed. The resulting orthogonal basic deformation modes exhibit simple and clear physical meanings. In addition, they do not involve any material parameters, and thus can be efficiently used to examine the element performance and serve as a unified tool to assess different hybrid elements. Thereafter, a convenient approach for the identification of spurious zero-energy modes is presented using the positive definiteness property of a flexibility matrix. Moreover, based on the orthogonality relationship between the given initial stress modes and the orthogonal basic deformation modes, an alternative method of assumed stress modes to formulate a hybrid element free of spurious modes is discussed. It is found that the orthogonality of the basic deformation modes is the sufficient and necessary condition for the suppression of spurious zero-energy modes. Numerical examples of 2D 4-node quadrilateral elements and 3D 8-node hexahedral elements are illustrated in detail to demonstrate the efficiency of the proposed orthogonal basic deformation mode method.     

Vortices and wall shear stresses in asymmetric and symmetric channels with sinusoidal wavy walls for pulsatile flow at low Reynolds numbers

Numerical simulation and flow visualization were performed to study the dynamical behavior of vortices generated in channels with two different geometries, i.e., a periodically converging–diverging channel and serpentine channel, both having sinusoidal wavy walls. This system for pulsatile flow is used to enhance heat and mass transfer in very viscous liquids. The numerical results predict well the dynamical behavior of vortices and agree with the flow visualizations. For both channels, the vortex expands in each furrow of the channel walls during the deceleration phase and shrinks during the acceleration phase, which leads to fluid exchange between the vortex and the mainstream. The time-averaged vortex strength and wall shear stresses increase, as the frequency of fluid oscillation increases under a fixed oscillatory fraction of the flow rate. However, above a certain value of the frequency, they reversely decrease due to viscous effects. This frequency for the serpentine channel is smaller than that for the converging–diverging channel. The channel geometries are found to have an important effect of the flow characteristics.     

Prediction of refrigerant void fraction in horizontal tubes using probabilistic flow regime maps

A state of the art review of two-phase void fraction models in smooth horizontal tubes is provided and a probabilistic two-phase flow regime map void fraction model is developed for refrigerants under condensation, adiabatic, and evaporation conditions in smooth, horizontal tubes. Time fraction information from a generalized probabilistic two-phase flow map is used to provide a physically based weighting of void fraction models for different flow regimes. The present model and void fraction models in the literature are compared to data from multiple sources including R11, R12, R134a, R22, R410A refrigerants, 4.26–9.58 mm diameter tubes, mass fluxes from 70 to 900 kg/m² s, and a full quality range. The present model has a mean absolute deviation of 3.5% when compared to the collected database.     

Nonlinear deformation of a particulate material in complex stress state

An algorithm is proposed to determine the effective deformation properties and stress-strain state of particulate composite materials with physically nonlinear components and complex stress state. The laws that govern the deformation of particulate composites are studied. A particulate composite is considered a two-component material of random structure. Its effective properties are determined by conditional averaging. The nonlinear equations that incorporate the physical nonlinearity of the components are solved by the method of successive approximations. The relationship between macrostresses and macrostrains is established. The effective deformation properties of a particulate composite as a function of the volume fractions of the components and stress state are studied __________ Translated from Prikladnaya Mekhanika, Vol. 42, No. 3, pp. 50–60, March 2006.     

Prediction of normal stresses under large amplitude oscillatory shear flow

The dynamic response of viscoelastic fluids under large amplitude oscillatory shear (LAOS) has been a subject of long history. In the LAOS flow, the analysis has been mostly focused on shear stress, possibly due to the lack of accurate measurement of normal stress. However, the normal stress may become larger than shear stress at high-strain amplitudes, and thus it is important that we have a good understanding of the normal stress behavior. Furthermore, with the advancement in the instrumentation, it has become possible to get more reliable data. The purpose of this paper is to develop a research platform to analyze and to understand the normal stress behavior of complex fluids under LAOS flow. In this study, we utilized the Giesekus model as a representative constitutive model, and investigated its diverse responses. We defined the dynamic properties corresponding to normal stress, in a similar way to define dynamic moduli from shear stress, and examine their behavior with various analyzing tools. Experimental data were also compared with model predictions. Despite the fact that it is not yet possible to compare all of the predictions because of instrumental limitation, the prediction has been found to fit well with the experimental data. This study is expected to provide a useful framework for further understanding the nonlinear behavior of complex fluids at large deformation.     

基于辛弹性力学解析本征函数的有限元应力磨平方法

在实际工程结构的结构强度与优化等力学数值分析中,应力计算结果的精度是非常重要的。有限元法是得到最广泛应用的一类数值方法,并形成了众多通用的有限元程序系统。这些程序系统采用的几乎都是基于最小总势能的位移法,虽然其分析给出的有限元位移场具有较高的精度,但所得到的有限元应力场的精度较位移场大大降低。基于极坐标辛对偶体系所提供的平面弹性力学的解析辛本征展开解,并借用有限元程序系统所得到的节点位移,本文提出了一个应力分析的改进方法。数值结果表明,本方法给出的应力分析精度得到大幅提高,并具有良好的数值稳定性,可用于有限元程序系统的后处理,以提高应力尤其是关键区域应力的分析精度。     

碳钢的高应力冲击磨损行为研究

在高应力冲击磨损条件下，金属材料的磨损表面很容易形成硬度很高（ＨＶ≥１０００）的白层．采用新的试验方法，研究了碳钢在高应力冲击条件下的磨损特征，并且从分析白层的结构出发，讨论了白层的形成机制和失效模式．结果表明，材料的高应力冲击磨损质量损失与冲击次数之间具有线性关系，韧性较好的低碳钢的耐磨性比高碳钢的好，原因是硬度越高，白层在磨损过程中越容易发生剥层与剥落．白层是一种高度塑性变形层，其内部精细结构的变形特征证实了白层的塑性变形机制     

某调距导管桨系泊工况水力性能和流激直发声预报

采用经验证的计算流体力学方法,对某艏辅推调距导管桨设计螺距和系泊工况螺距的水动力性能进行了有效预报,并对系泊工况装船桨流激噪声进行了分析。系泊工况下,由于导管桨的抽吸作用在导管外壁近壁面区域存在与导管内部流动方向相反的逆向流动,且导管桨尾流场速度梯度分布不均匀、流动紊乱,此时桨叶与导管的推力之比约为1.2∶1。系泊工况船+桨的瞬态流场脉动信息表明,导管桨各部件噪声源强度均表现出从1倍到4倍叶频依次下降的规律,最强幅值集中在桨叶导边和导管内壁;在远场声源级频谱曲线中轴向测点线谱较高峰值位置体现出导管桨进流流场的流动特性。对比分析该艏辅推整体和各部件宽带声源级指向性,可知旋转部件(桨叶、桨榖)对总噪声级的贡献较大,静止部件是径向测点噪声的主要贡献源。