scholarly journals A three-pronged approach to predict the effect of plastic orthotropy on the formability of thin sheets subjected to dynamic biaxial stretching

2020 ◽  
Author(s):  
Jose Rodriguez-Martinez ◽  
Komi Espoir N'souglo ◽  
Nicolas Jacques

In this paper, we have investigated the effect of material orthotropy on the formability of metallic sheets subjectedto dynamic biaxial stretching. For that purpose, we have devised an original three-pronged methodology which includes a linear stability analysis, a nonlinear two-zone model and ?finite element calculations. We have studied 5 different materials whose mechanical behavior is described with an elastic isotropic, plastic anisotropic constitutive model with yielding based on Hill (1948) criterion. The linear stability analysis and the nonlinear two-zone model are extensions of the formulations developed by Zaera et al. (2015) and Jacques (2020), respectively, to consider Hill (1948) plasticity. The ?finite element calculations are performed with ABAQUS/Explicit (2016) using the unit-cell model developed by Rodriguez-Martinez et al. (2017), which includes a sinusoidal spatial imperfection to favor necking localization. The predictions of the stability analysis and the two-zone model are systematically compared against the ?finite element results --which are considered as the reference approach to validate the theoretical models-- for loading paths ranging from plane strain stretching to equibiaxial stretching, and for different strain rates ranging from 100 s-1 to 50000 s-1. The stability analysis and the two-zone model yield the same overall trends obtained with the ?finite element simulations for the 5 materials investigated, and for most of the strain rates and loading paths the agreement for the necking strains is also quantitative. Notably, the differences between the ?finite element results and the two-zone model rarely go beyond 5%. Altogether, the results presented in this work provide new insights into the mechanisms which control dynamic formability of anisotropic metallic sheets.

2001 ◽  
Author(s):  
Davide Valtorta ◽  
Khaled E. Zaazaa ◽  
Ahmed A. Shabana ◽  
Jalil R. Sany

Abstract The lateral stability of railroad vehicles travelling on tangent tracks is one of the important problems that has been the subject of extensive research since the nineteenth century. Early detailed studies of this problem in the twentieth century are the work of Carter and Rocard on the stability of locomotives. The linear theory for the lateral stability analysis has been extensively used in the past and can give good results under certain operating conditions. In this paper, the results obtained using a linear stability analysis are compared with the results obtained using a general nonlinear multibody methodology. In the linear stability analysis, the sources of the instability are investigated using Liapunov’s linear theory and the eigenvalue analysis for a simple wheelset model on a tangent track. The effects of the stiffness of the primary and secondary suspensions on the stability results are investigated. The results obtained for the simple model using the linear approach are compared with the results obtained using a new nonlinear multibody based constrained wheel/rail contact formulation. This comparative numerical study can be used to validate the use of the constrained wheel/rail contact formulation in the study of lateral stability. Similar studies can be used in the future to define the limitations of the linear theory under general operating conditions.


Author(s):  
K. E. N’souglo ◽  
A. Srivastava ◽  
S. Osovski ◽  
J. A. Rodríguez-Martínez

At high strain rates, the fragmentation of expanding structures of ductile materials, in general, starts by the localization of plastic deformation in multiple necks. Two distinct mechanisms have been proposed to explain multiple necking and fragmentation process in ductile materials. One view is that the necking pattern is related to the distribution of material properties and defects. The second view is that it is due to the activation of specific instability modes of the structure. Following this, we investigate the emergence of necking patterns in porous ductile bars subjected to dynamic stretching at strain rates varying from 10 3   s −1 to 0.5×10 5   s −1 using finite-element calculations and linear stability analysis. In the calculations, the initial porosity (representative of the material defects) varies randomly along the bar. The computations revealed that, while the random distribution of initial porosity triggers the necking pattern, it barely affects the average neck spacing, especially, at higher strain rates. The average neck spacings obtained from the calculations are in close agreement with the predictions of the linear stability analysis. Our results also reveal that the necking pattern does not begin when the Considère condition is reached but is significantly delayed due to the stabilizing effect of inertia.


2017 ◽  
Vol 822 ◽  
pp. 813-847 ◽  
Author(s):  
Azan M. Sapardi ◽  
Wisam K. Hussam ◽  
Alban Pothérat ◽  
Gregory J. Sheard

This study seeks to characterise the breakdown of the steady two-dimensional solution in the flow around a 180-degree sharp bend to infinitesimal three-dimensional disturbances using a linear stability analysis. The stability analysis predicts that three-dimensional transition is via a synchronous instability of the steady flows. A highly accurate global linear stability analysis of the flow was conducted with Reynolds number $\mathit{Re}<1150$ and bend opening ratio (ratio of bend width to inlet height) $0.2\leqslant \unicode[STIX]{x1D6FD}\leqslant 5$. This range of $\mathit{Re}$ and $\unicode[STIX]{x1D6FD}$ captures both steady-state two-dimensional flow solutions and the inception of unsteady two-dimensional flow. For $0.2\leqslant \unicode[STIX]{x1D6FD}\leqslant 1$, the two-dimensional base flow transitions from steady to unsteady at higher Reynolds number as $\unicode[STIX]{x1D6FD}$ increases. The stability analysis shows that at the onset of instability, the base flow becomes three-dimensionally unstable in two different modes, namely a spanwise oscillating mode for $\unicode[STIX]{x1D6FD}=0.2$ and a spanwise synchronous mode for $\unicode[STIX]{x1D6FD}\geqslant 0.3$. The critical Reynolds number and the spanwise wavelength of perturbations increase as $\unicode[STIX]{x1D6FD}$ increases. For $1<\unicode[STIX]{x1D6FD}\leqslant 2$ both the critical Reynolds number for onset of unsteadiness and the spanwise wavelength decrease as $\unicode[STIX]{x1D6FD}$ increases. Finally, for $2<\unicode[STIX]{x1D6FD}\leqslant 5$, the critical Reynolds number and spanwise wavelength remain almost constant. The linear stability analysis also shows that the base flow becomes unstable to different three-dimensional modes depending on the opening ratio. The modes are found to be localised near the reattachment point of the first recirculation bubble.


1997 ◽  
Vol 352 ◽  
pp. 265-281 ◽  
Author(s):  
A. M. H. BROOKER ◽  
J. C. PATTERSON ◽  
S. W. ARMFIELD

A non-parallel linear stability analysis which utilizes the assumptions made in the parabolized stability equations is applied to the buoyancy-driven flow in a differentially heated cavity. Numerical integration of the complete Navier–Stokes and energy equations is used to validate the non-parallel theory by introducing an oscillatory heat input at the upstream end of the boundary layer. In this way the stability properties are obtained by analysing the evolution of the resulting disturbances. The solutions show that the spatial growth rate and wavenumber are highly dependent on the transverse location and the disturbance flow quantity under consideration. The local solution to the parabolized stability equations accurately predicts the wave properties observed in the direct simulation whereas conventional parallel stability analysis overpredicts the spatial amplification and the wavenumber.


2013 ◽  
Vol 736 ◽  
pp. 464-494 ◽  
Author(s):  
P. Pearce ◽  
J. Daou

AbstractWe investigate the Rayleigh–Bénard convection problem within the context of a diffusion flame formed in a horizontal channel where the fuel and oxidizer concentrations are prescribed at the porous walls. This problem seems to have received no attention in the literature. When formulated in the low-Mach-number approximation the model depends on two main non-dimensional parameters, the Rayleigh number and the Damköhler number, which govern gravitational strength and reaction speed respectively. In the steady state the system admits a planar diffusion flame solution; the aim is to find the critical Rayleigh number at which this solution becomes unstable to infinitesimal perturbations. In the Boussinesq approximation, a linear stability analysis reduces the system to a matrix equation with a solution comparable to that of the well-studied non-reactive case of Rayleigh–Bénard convection with a hot lower boundary. The planar Burke–Schumann diffusion flame, which has been previously considered unconditionally stable in studies disregarding gravity, is shown to become unstable when the Rayleigh number exceeds a critical value. A numerical treatment is performed to test the effects of compressibility and finite chemistry on the stability of the system. For weak values of the thermal expansion coefficient $\alpha $, the numerical results show strong agreement with those of the linear stability analysis. It is found that as $\alpha $ increases to a more realistic value the system becomes considerably more stable, and also exhibits hysteresis at the onset of instability. Finally, a reduction in the Damköhler number is found to decrease the stability of the system.


2016 ◽  
Vol 27 (05) ◽  
pp. 1650050 ◽  
Author(s):  
Guanghan Peng

A new lattice model is proposed by taking into account the interruption probability with passing for two-lane freeway. The effect of interruption probability with passing is investigated about the linear stability condition and the mKdV equation through linear stability analysis and nonlinear analysis, respectively. Furthermore, numerical simulation is carried out to study traffic phenomena resulted from the interruption probability with passing in two-lane system. The results show that the interruption probability with passing can improve the stability of traffic flow for low reaction coefficient while the interruption probability with passing can destroy the stability of traffic flow for high reaction coefficient on two-lane highway.


Author(s):  
K.B. Tsiberkin

The stability of incompressible fluid plane-parallel flow over a layer of a saturated porous medium is studied. The results of a linear stability analysis are described at different porosity values. The considered system is bounded by solid wall from the porous layer bottom. Top fluid surface is free and rigid. A linear stability analysis of plane-parallel stationary flow is presented. It is realized for parameter area where the neutral stability curves are bimodal. The porosity variation effect on flow stability is considered. It is shown that there is a transition between two main instability modes: long-wave and short-wave. The long-wave instability mechanism is determined by inflection points within the velocity profile. The short-wave instability is due to the large transverse gradient of flow velocity near the interface between liquid and porous medium. Porosity decrease stabilizes the long wave perturbations without significant shift of the critical wavenumber. Simultaneously, the short-wave perturbations destabilize, and their critical wavenumber changes in wide range. When the porosity is less than 0.7, the inertial terms in filtration equation and magnitude of the viscous stress near the interface increase to such an extent that the Kelvin-Helmholtz analogue of instability becomes the dominant mechanism for instability development. The stability band realizes in narrow porosity area. It separates the two branches of the neutral curve.


Author(s):  
A Amamou ◽  
M Chouchane

Floating ring bearings are used to support and guide rotors in several high-speed rotating machinery applications. They are usually credited for lower heat generation and higher vibration suppressing ability. Similar to conventional hydrodynamic bearings, floating ring bearings may exhibit unstable behaviour above a certain stability critical speed. Linear stability analysis is usually applied to predict the stability threshold speed. Non-linear stability analysis, however, is needed to predict the presence and the size of stable limit cycles above the stability threshold speed or unstable limit cycles below the stability critical speed. The prediction of limit cycles is an important step in bearing stability analysis. In this article, a non-linear dynamic model is derived and used to investigate the stability of a perfectly balanced symmetric rigid rotor supported by two identical floating ring bearings near the critical stability boundaries. The fluid film hydrodynamic reactions of the floating ring bearings are modelled by applying the short bearing theory and the half Sommerfeld solution. Hopf bifurcation theory is then utilized to determine the existence and the approximate size of stable and unstable limit cycles in the neighbourhood of the stability critical speed depending on the bearing design parameters. Numerical integration of the non-linear equations of motion is then carried out in order to compare the trajectories obtained by numerical integration to those obtained analytically using Hopf bifurcation analysis. Stability boundary curves for typical bearing design parameters have been decomposed into boundaries with supercritical stable limit cycles and boundaries with subcritical unstable limit cycles. The shape and size of the limit cycles for selected bearing parameters are presented using both analytical and numerical approaches. This article shows that floating ring stability boundaries may exhibit either stable supercritical limit cycles or unstable subcritical limit cycles predictable by Hopf bifurcation.


2021 ◽  
Author(s):  
Jose Rodriguez-Martinez ◽  
Oana Cazacu ◽  
Nitin Chandola ◽  
Komi Espoir N'souglo

In this paper, we have investigated the effect of the third invariant of the stress deviator on the formation of necking instabilities in isotropic metallic plates subjected to plane strain tension. For that purpose, we have performed finite element calculations and linear stability analysis for initial equivalent strain rates ranging from 10^−4 s−1 to 8 · 10^4 s−1. The plastic behavior of the material has been escribed with the isotropic Drucker yield criterion [11], which depends on both the second and third invariant of the stress deviator, and a parameter c which determines the ratio between the yield stresses in uniaxial tension and in pure shear \sigma_T /\tau_Y . For c = 0, Drucker yield criterion [11] reduces to the von Mises yield criterion [32] while for c = 81/66, the Hershey-Hosford (m = 6) yield criterion [19, 22] is recovered. The results obtained with both finite element calculations and linear stability analysis show the same overall trends and there is also quantitative agreement for most of the loading rates considered. In the quasi-static regime, while the specimen elongation when necking occurs is virtually insensitive to the value of the parameter c, both finite element results and analytical calculations using Considère criterion [10] show that the necking strain increases as the parameter c decreases, bringing out the effect of the third invariant of the stress deviator on the formation of quasi-static necks. In contrast, at high initial equivalent strain rates, when the influence of inertia on the necking process becomes important, both finite element simulations and linear stability analysis show that the effect of the third invariant is reversed, notably for long necking wavelengths, with the specimen elongation when necking occurs increasing as the parameter c increases, and the necking strain decreasing as the parameter c decreases.


Sign in / Sign up

Export Citation Format

Share Document