Investigation on mechanical properties and size effect of nanocrystalline twin copper

The stress-strain relations of nanocrystalline twin copper with variously sized grains and twins are studied by using FEM simulations based on the conventional theory of mechanism-based strain gradient plasticity (CMSG). A model of twin lamellae strengthening zone is proposed and a cohesive interface model is used to simulate grain-boundary sliding and separation. Effects of material parameters on stress-strain curves of polycrystalline twin copper are studied in detail. Furthermore, the effects of both twin lamellar spacing and twin lamellar distribution on the stress-strain relations are investigated under tension loading. The numerical simulations show that both the strain gradient effect and the material hardening increase with decreasing the grain size and twin lamellar spacing. The distribution of twin lamellae has a significant influence on the overall mechanical properties, and the effect is reduced as both the grain size and twin lamellar spacing decrease. Finally, the FEM prediction results are compared with the experimental data.

Energy principle of ferroelectric ceramics and single domain mechanical model

Many physical experiments have shown that the domain switching in a ferroelectric material is a complicated evolution process of the domain wall with the variation of stress and electric field. According to this mechanism, the volume fraction of the domain switching is introduced in the constitutive law of ferroelectric ceramic and used to study the nonlinear constitutive behavior of ferroelectric body in this paper. The principle of stationary total energy is put forward in which the basic unknown quantities are the displacement u (i) , electric displacement D (i) and volume fraction rho (I) of the domain switching for the variant I. Mechanical field equation and a new domain switching criterion are obtained from the principle of stationary total energy. The domain switching criterion proposed in this paper is an expansion and development of the energy criterion. On the basis of the domain switching criterion, a set of linear algebraic equations for the volume fraction rho (I) of domain switching is obtained, in which the coefficients of the linear algebraic equations only contain the unknown strain and electric fields. Then a single domain mechanical model is proposed in this paper. The poled ferroelectric specimen is considered as a transversely isotropic single domain. By using the partial experimental results, the hardening relation between the driving force of domain switching and the volume fraction of domain switching can be calibrated. Then the electromechanical response can be calculated on the basis of the calibrated hardening relation. The results involve the electric butterfly shaped curves of axial strain versus axial electric field, the hysteresis loops of electric displacement versus electric filed and the evolution process of the domain switching in the ferroelectric specimens under uniaxial coupled stress and electric field loading. The present theoretic prediction agrees reasonably with the experimental results given by Lynch.

A coupled model for multiphase fluid flow and sedimentation deformation in oil reservoir and its numerical simulation

A mathematical model for coupled multiphase fluid flow and sedimentation deformation is developed based on fluid-solid interaction mechanism. A finite difference-finite element numerical approach is presented. The results of an example show that the fluid-solid coupled effect has great influence on multiphase fluid flow and reservoir recovery performances, and the coupled model has practical significance for oilfield development.

Three-dimensional crack problem analysis using boundary element method with finite-part integrals

A set of hypersingular integral equations of a three-dimensional finite elastic solid with an embedded planar crack subjected to arbitrary loads is derived. Then a new numerical method for these equations is proposed by using the boundary element method combined with the finite-part integral method. According to the analytical theory of the hypersingular integral equations of planar crack problems, the square root models of the displacement discontinuities in elements near the crack front are applied, and thus the stress intensity factors can be directly calculated from these. Finally, the stress intensity factor solutions to several typical planar crack problems in a finite body are evaluated.

Study of three elastic-plastic constitutive models by non-proportional finite deformations of OFHC copper

Experimental stress-strain data of OFHC copper first under torsion to 13% and then under torsion-tension to about 10% are used to study the characteristics of three elastic-plastic constitutive models: Chaboche's super-positional nonlinear model, Dafalias and Popov's two surface model and Watanabe and Atluri's version of the endochronic model. The three models, originally oriented for infinitesimal deformation, have been extended for finite deformation. The results show (a) the Mises-type yield surface used in the three models brings about significant departure of the predictions from the experimental data; (b) Chaboche's and Dafalias' models are easier than Watanabe and Atluri's model in determining the material parameters in them, and (c) Chaboche's and Watanabe & Atluri's models produce almost the same prediction to the data, while Dafalias' model cannot accurately predict the plastic deformations when a loading path changes in its direction. Copyright (C) 1996 Elsevier Science Ltd

A finite deformation analysis of polycrystal by considering the influence of grain boundary

By combining grain boundary (GB) and its influence zone, a micromechanic model for polycrystal is established for considering the influence of GB. By using the crystal plasticity theory and the finite element method for finite deformation, numerical simulation is carried out by the model. Calculated results display the microscopic characteristic of deformation fields of grains and are in qualitative agreement with experimental results.

Micromechanics analysis of particulate-reinforced composites and their failure mechanisms

Stress fields and failure mechanisms have been investigated in composites with particles either surface treated or untreated under uniaxial tension. Previous experimental observation of failure mechanisms in a composite with untreated particles showed that tensile cracks occurred mostly at the polar region of the particle and grew into interfacial debonding. In a composite with surface-treated particles, however, shear yielding and shear cracking proceeded along the interphase-matrix interface at the polar area of the matrix and thus may improve the mechanical behaviour of the material. The finite element calculations showed that octahedral shear stress at the polar and longitudinal areas of the particle treated by coupling agents is much larger than that of materials with untreated particles, and the shear stress distribution around the interface is sensitive to the interphase property. The results suggest that a th ree-phase model can describe the composites with surface-treated fillers.

Three-dimensional finite element analysis of composites with coated spherical inclusions

A three-dimensional finite element analysis has been used to determine the internal stresses in a three-phase composite. The stresses have been determined for a variety of interphase properties, the thicknesses of the interphase and the volume fractions of particles. Young's modulus has been calculated from a knowledge of these stresses and the applied deformation. The calculations show that stress distributions in the matrix and the mechanical properties are sensitive to the interphase property in the three-phase composites. The interfacial stresses in the three-dimensional analysis are in agreement with results obtained by an axisymmetric analysis. The predicted bulk modulus in three-dimensional analysis agrees well with the theoretical solution obtained by Qui and Weng, but it presents a great divergence from that in axisymmetric analyses. An investigation indicates that this divergence may be caused by the difference in the unit cell structure between two models. A comparison of the numerically predicted bulk and shear modulus for two-phase composites with the theoretical results indicates that the three-dimensional analysis gives quite satisfactory results.

Triple-region structure for turbulent-flow in a square duct - a finite-element approach

The Reynolds-averaged Navier-Stokes equations for describing the turbulent flow in a straight square duct are formulated with two different turbulence models. The governing equations are then expanded as a multi-deck structure in a plane perpendicular to the streamwise direction, with each deck characterized by its dominant physical forces as commonly carried out in analytical work using triple-deck expansion. The resulting equations are numerically integrated using higher polynomial (H-P) finite element technique for each cross-sectional plane to be followed by finite difference representation in the streamwise direction until a fully developed state is reached. The computed results using the two different turbulence models show fair agreement with each other, and concur with the vast body of available experimental data. There is also general agreement between our results and the recent numerical works anisotropic k-epsilon turbulence model.

A perturbational h4 exponential finite-difference scheme for the convective diffusion equation

A perturbational h4 compact exponential finite difference scheme with diagonally dominant coefficient matrix and upwind effect is developed for the convective diffusion equation. Perturbations of second order are exerted on the convective coefficients and source term of an h2 exponential finite difference scheme proposed in this paper based on a transformation to eliminate the upwind effect of the convective diffusion equation. Four numerical examples including one- to three-dimensional model equations of fluid flow and a problem of natural convective heat transfer are given to illustrate the excellent behavior of the present exponential schemes, the h4 accuracy of the perturbational scheme is verified using double precision arithmetic.

Domain decomposition hybrid method for numerical simulation of bluff body flows:theoretical model and application

The discrete vortex method is not capable of precisely predicting the bluff body flow separation and the fine structure of flow field in the vicinity of the body surface. In order to make a theoretical improvement over the method and to reduce the difficulty in finite-difference solution of N-S equations at high Reynolds number, in the present paper, we suggest a new numerical simulation model and a theoretical method for domain decomposition hybrid combination of finite-difference method and vortex method. Specifically, the full flow. field is decomposed into two domains. In the region of O(R) near the body surface (R is the characteristic dimension of body), we use the finite-difference method to solve the N-S equations and in the exterior domain, we take the Lagrange-Euler vortex method. The connection and coupling conditions for flow in the two domains are established. The specific numerical scheme of this theoretical model is given. As a preliminary application, some numerical simulations for flows at Re=100 and Re-1000 about a circular cylinder are made, and compared with the finite-difference solution of N-S equations for full flow field and experimental results, and the stability of the solution against the change of the interface between the two domains is examined. The results show that the method of the present paper has the advantage of finite-difference solution for N-S equations in precisely predicting the fine structure of flow field, as well as the advantage of vortex method in efficiently computing the global characteristics of the separated flow. It saves computer time and reduces the amount of computation, as compared with pure N-S equation solution. The present method can be used for numerical simulation of bluff body flow at high Reynolds number and would exhibit even greater merit in that case.

A consistent model of isolated force-free magnetic arches

The magnetic flux tube concentrating strong magnetic field is the basic configuration of magneticfield in the solar atmosphere. In the present paper, the equilibrium of isolated magnetic flux tube inthe solar atmosphere is discussed. In the viewpoint of mathematics, the boundary condition is nonlinearand the position of boundary needs to be determined by the physical condition although the equation ofmagnetic potential is linear for the linear force-free field. Analytical solutions to the arches of bothuniform circular cross-section and non-uniform cross section have been obtained. The results show thatthe nonlinear problem may have or not have any solution according to different azimuthal components of the magnetic field; the number of solutions to the nonlinear problem is four at most, and two in some cases. In the present paper, the analytical solutions to the approximations of both fat and slender arches are given in detail, and the general features of magnetic arch structure are shown.

Characteristics of a CW HF chemical-laser calculated from a simplified two-dimensional model

A two-dimensional simplified model of an HF chemical laser is introduced. Using an implicit finite difference scheme, the solution of two adjacent parallel streams with diffusion mixing and chemical reaction is generated. A contour of mixing and reaction boundary is obtained without presupposition. The distribution of the HF(v) concentrations, gas temperature and the optical small signal gain (alpha sub V, J) on the flowing plane (X, Y) are presented. Compared with the solution solved directly from a set of Navier-Stokes equations, the results of these two methods agree with each other qualitatively. The influences of the different velocity, temperature (T sub 0) and composition of the two streams on the small signal gain after the nozzle exit are investigated. It is interesting that for larger J with a fixed v, the peaks of alpha sub v-T sub 0 profiles move towards higher T sub 0. The computing method is simple and only a short computing time is needed.

Starting Nozzle Flow Simulation Using K-G Two-Equation Turbulence Model

The starting process of two-dimensional nozzle flows has been simulated with Euler, laminar and k - g two-equation turbulence Navier-Stokes equations. The flow solver is based on a combination of LUSGS subiteration implicit method and five spatial discretized schemes, which are Roe, HLLE, MHLLE upwind schemes and AUSM+, AUSMPW schemes. In the paper, special attention is for the flow differences of the nozzle starting process obtained from different governing equations and different schemes. Two nozzle flows, previously investigated experimentally and numerically by other researchers, are chosen as our examples. The calculated results indicate the carbuncle phenomenon and unphysical oscillations appear more or less near a wall or behind strong shock wave except using HLLE scheme, and these unphysical phenomena become more seriously with the increase of Mach number. Comparing the turbulence calculation, inviscid solution cannot simulate the wall flow separation and the laminar solution shows some different flow characteristics in the regions of flow separation and near wall.

Measurement of Electron Density Profile Behind Strong Shock Waves with a Langmuir Probe

At the shock velocity range of 7~9km/s, the variations of electron density behind strong normal shock waves are measured in a low-density shock tube by using the Langmuir electrostatic probe technique. The electron temperature, calculated based on Park’s three-temperature model, is used in interpreting the probe current data. The peak electron densities determined in the present experiment are shown to be in a good agreement with those predicted by Lin’s calculation. The experimentally obtained ratios of the characteristic ionization distance to the mean free path of freestream ahead of the shock wave are found to be in a good agreement with the existing experiments and Park’s calculation.