Part I of the present work dealt with the initial stage of plastic deformation with grain subdivision into two band families. It covered the situation with different average strains in the individual grains (the non-T...
详细信息
Part I of the present work dealt with the initial stage of plastic deformation with grain subdivision into two band families. It covered the situation with different average strains in the individual grains (the non-Taylor case) and the situation with identical average strains in the individual grains (the Taylor case). For the latter situation part I included solutions with 5 plus 3 and with 4 plus 4 active slip systems in the two band families, respectively. In part II we deal with the subsequent stage of plastic deformation (finite strains) for the Taylor case with four active slip systems in each band family. In the subsequent stage the cooperation between the two band families leads to an energetic advantage for grain subdivision, and it leads to a simulated texture with decreased sharpness in accordance with experimental observations. Finally, the physics behind grain subdivision are discussed on the basis of the results in Parts I and II and various general considerations. (C) 2001 Elsevier Science Ltd. All rights reserved.
This paper deals with a thermodynamically consistent numerical formulation for coupled thermoplastic problems including phase-change phenomena and frictional contact. The final goal is to get an accurate, efficient an...
详细信息
This paper deals with a thermodynamically consistent numerical formulation for coupled thermoplastic problems including phase-change phenomena and frictional contact. The final goal is to get an accurate, efficient and robust numerical model, able for the numerical simulation of industrial solidification processes. Some of the current issues addressed in the paper are the following. A fractional step method arising from an operator split of the governing differential equations has been used to solve the nonlinear coupled system of equations, leading to a staggered product formula solution algorithm. Nonlinear stability issues are discussed and isentropic and isothermal operator splits are formulated. Within the isentropic split, a strong operator split design constraint is introduced, by requiring that the elastic and plastic entropy, as well as the phase-change induced elastic entropy due to the latent heat, remain fixed in the mechanical problem. The formulation of the model has been consistently derived within a thermodynamic framework. All the material properties have been considered to be temperature dependent. The constitutive behavior has been defined by a thermoviscous/ elastoplastic free energy function, including a thermal multiphase change contribution. Plastic response has been modeled by a J2 temperature dependent model, including plastic hardening and thermal softening. The constitutive model proposed accounts for a continuous transition between the initial liquid state, the intermediate mushy state and the final solid state taking place in a solidification process. In particular, a pure viscous deviatoric model has been used at the initial fluid-like state. A thermomecanical contact model, including a frictional hardening and temperature dependent coupled potential, is derived within a fully consistent thermodinamical theory. The numerical model has been implemented into the computational finite element code COMET developed by the authors. numerical s
A problem very often arising in applications is presented: finding the minimal root of an equation with the objective function being multiextremal and nondifferentiable. Applications from the field of electronic measu...
详细信息
A problem very often arising in applications is presented: finding the minimal root of an equation with the objective function being multiextremal and nondifferentiable. Applications from the field of electronic measurements are given. Three methods based on global optimization ideas are introduced for solving this problem. The first one uses an a priori estimate of the global Lipschitz constant. The second method adaptively estimates the global Lipschitz constant. The third algorithm adaptively estimates local Lipschitz constants during the search. All the methods either find the minimal root or determine the global minimizers (in the case when the equation under consideration has no roots). Sufficient convergence conditions of the new methods to the desired solution are established. numerical results including wide experiments with test functions, stability study, and a real-life applied problem are also presented.
A detailed survey is presented of the literature on attitude representation dating from the early work of Euler and Hamilton to recent publications in fields such as navigation and control. The scope is limited to the...
详细信息
A detailed survey is presented of the literature on attitude representation dating from the early work of Euler and Hamilton to recent publications in fields such as navigation and control. The scope is limited to the development of the aircraft kinematic transformation equations in terms of four different attitude representations, including the well-known Enter angles, the Euler-axis rotation parameters, the direction cosines, and the Euler-Rodrigues quaternion. The emphasis is directed at the application of the quaternion formulation to aircraft kinematics. Results are presented that reinforce observations that the quaternion formulation, typically implemented to eliminate singularities associated with the Enter angle formulation, is far superior to the other commonly used formulations based on computational efficiency alone. A development of quaternion constraints necessary to independently constrain roll, pitch, yaw, bank angle, elevation angle, and/or azimuth angle is presented. For verification of simulation codes, a general closed-form solution to the quaternion formulation, for the case of constant rotation, is also presented. Additionally, a discussion is provided of numerical integration methods and numerical errors for the quaternion formulation. This discussion is especially important for simulations that may still utilize a common error reduction scheme originally developed for analog computers.
The discrete-time optimal projection equations for systems with white parameters are strengthened. For the class of minimal ms (mean square) stabilizing compensators the strengthened discrete-time optimal projection e...
详细信息
ISBN:
(纸本)085296708X
The discrete-time optimal projection equations for systems with white parameters are strengthened. For the class of minimal ms (mean square) stabilizing compensators the strengthened discrete-time optimal projection equations are proved to be equivalent to first-order necessary optimality conditions for optimal reduced-order dynamic compensation of systems with white parameters. The conventional discrete-time optimal projection equations are proved to be weaker. As a result solutions of the conventional discrete-time optimal projection equations may not correspond to optimal reduced-order compensators. To compute optimal reduced-order compensators two numerical algorithms are proposed. One is a homotopy algorithm and one is based on iteration of the strengthened discrete-time optimal projection equations. The latter algorithm is a generalization of the algorithm that solves the full-order problem, which in turn is a generalization of the algorithm that solves the two Riccati equations of full-order LQG control through iteration. Therefore the efficiency of these three types of algorithms is comparable. It is demonstrated that, despite the strengthening of the optimal projection equations, the optimal reduced-order compensation problem, in general, may posses multiple extrema.
Diffraction problems in analytical treatments and numeric have been constant topics for a very long time and there are a lot of different analytical expressions as well as numerical algorithms. The modern approaches a...
详细信息
ISBN:
(纸本)0819430382
Diffraction problems in analytical treatments and numeric have been constant topics for a very long time and there are a lot of different analytical expressions as well as numerical algorithms. The modern approaches are concentrated in the area of particle sizing, fibber optical sizing and loss evaluation in laser physics and techniques. Some problems between the fibber optics dimensioning and particles of cylindrical shape dimensioning are chosen and based on literature;some ameliorations in numerical approaches are made. The evolved expressions can be used for different laser wavelengths, different dimensions of disturbance objects tin scattering and diffraction phenomena) as well as different indices of refractions where the particles tin a very general way) are stratified. Some existing programs for scattering are discussed and the new computer architecture is presented here.
The design of an algorithm for the numerical solution of a generalized quadratic algebraic matrix Riccati equation is presented. The approach is based on probability-1 homotopy methods. The algorithm is illustrated wi...
详细信息
The design of an algorithm for the numerical solution of a generalized quadratic algebraic matrix Riccati equation is presented. The approach is based on probability-1 homotopy methods. The algorithm is illustrated with numerical examples.
Some numerical solutions of acoustic propagation problems using linearized Euler equations are studied. The two-dimensional Euler equations are linearized around a known stationary mean flow. The computed solution is ...
详细信息
Some numerical solutions of acoustic propagation problems using linearized Euler equations are studied. The two-dimensional Euler equations are linearized around a known stationary mean flow. The computed solution is obtained by using is dispersion-relation-preserving scheme in space, combined with a fourth-order Runge-Kutta algorithm in time. This numerical integration leads to very good results in terms of accuracy, stability, and low storage. The implementation of source terms in these equations is studied very carefully in various configurations, inasmuch as the final goal is to improve and to validate the stochastic noise generation and radiation model. In this approach, the turbulent velocity field is modeled by a sum of random Fourier modes through a source term in the linearized Euler equations to predict the noise from subsonic flows. The radiation of a point source in a subsonic and a supersonic uniform mean flow is investigated. The numerical estimates are shown to be in excellent agreement with the analytical solutions. Then, the emphasis is on the ability of the method to describe correctly the multipolar structure of aeroacoustic sources. The radiation of dipolar and guadrupolar extended sources is, thus, studied. Next, a typical problem in jet noise is considered with the propagation of acoustic waves in a sheared mean flow. The numerical solution compares favorably with ray tracing. Finally, a nonlinear formulation of Euler's equations is solved to limit the growth of instability waves excited by the acoustic source terms.
We consider the changing face of computing and the evolving interface of numerical algorithms from subroutines to: spreadsheets, high-level scripting languages, visual programming, Web browsers. Through examples we co...
详细信息
We consider the changing face of computing and the evolving interface of numerical algorithms from subroutines to: spreadsheets, high-level scripting languages, visual programming, Web browsers. Through examples we consider the movement from the subroutine language interface to sophisticated computational environments and trace how NAG libraries and products have evolved along this path. (C) 2000 IMACS. Published by Elsevier Science B.V. All rights reserved.
作者:
Yang, KSInha Univ
Dept Engn Mech Nam Ku Inchon 402751 South Korea
Instability and transition in an obstructed channel how are investigated using direct numerical simulation. The flow geometry under consideration is a plane channel with two-dimensional thin obstacles (baffles) mounte...
详细信息
Instability and transition in an obstructed channel how are investigated using direct numerical simulation. The flow geometry under consideration is a plane channel with two-dimensional thin obstacles (baffles) mounted symmetrically in the vertical direction and periodically in the streamwise direction. The flow is steady and symmetric at low Reynolds numbers. Above a certain Reynolds number, the flow undergoes a Hopf bifurcation, leading to a solution periodic in time. At high Reynolds numbers the unsteady flow exhibits a space-time symmetry-breaking bifurcation. A secondary instability is also observed at high Reynolds numbers, which is believed to be responsible for a subsequent chaotic breakdown of the flow. To study the secondary instability, we take the periodic solution, which results from the Hopf bifurcation, as the basic flow to be disturbed with small-amplitude random noise. Depending on the Reynolds number the basic Row becomes unstable to three-dimensional disturbances, which results in a chaotic how The numerical results obtained are consistent with experimental findings currently available.
暂无评论