The application of B-P constitutive equations in finite element analysis of high velocity impact

We present in this paper the application of B-P constitutive equations in finite element analysis of high velocity impact. The impact process carries out in so quick time that the heat-conducting can be neglected and meanwhile, the functions of temperature in equations need to be replaced by functions of plastic work. The material constants in the revised equations can be determined by comparison of the one-dimensional calculations with the experiments of Hopkinson bar. It can be seen from the comparison of the calculation with the experiment of a tungsten alloy projectile impacting a three-layer plate that the B-P constitutive equations in that the functions of temperature were replaced by the functions of plastic work can be used to analysis of high velocity impact.

Dynamic constitutive equation of GFRP obtained by Lagrange experiment

The note presents a method of constructing dynamic constitutive equations of material by means of Lagrange experiment and analysis. Tests were carried out by a light gas gun and the stress history profiles were recorded on multiple Lagrange positions. The dynamic constitutive equations were deduced from the regression of a series of data which was obtained by Lagrange analysis based upon recorded multiple stress histories. Here constitutive equations of glass fibre reinforced phenolic resin composite(GFRP) in uniaxil strain state under dynamic loading are given. The proposed equations of the material agree well with experimental results.

Energy Functions And Basic Equations For Ferroelectrics

Energy functions (or characteristic functions) and basic equations for ferroelectrics in use today are given by those for ordinary dielectrics in the physical and mechanical communications. Based on these basic equations and energy functions, the finite element computation of the nonlinear behavior of the ferroelectrics has been carried out by several research groups. However, it is difficult to process the finite element computation further after domain switching, and the computation results are remarkably deviating from the experimental results. For the crack problem, the iterative solution of the finite element calculation could not converge and the solutions for fields near the crack tip oscillate. In order to finish the calculation smoothly, the finite element formulation should be modified to neglect the equivalent nodal load produced by spontaneous polarization gradient. Meanwhile, certain energy functions for ferroelectrics in use today are not compatible with the constitutive equations of ferroelectrics and need to be modified. This paper proposes a set of new formulae of the energy functions for ferroelectrics. With regard to the new formulae of the energy functions, the new basic equations for ferroelectrics are derived and can reasonably explain the question in the current finite element analysis for ferroelectrics.

Component Assembling Model and Its Application to Quasi-Brittle Damage

Potential energy can be approximated by ‘‘pair-functional’’ potentials which is composed of pair potentials and embedding energy. Pair potentials are grouped according to discrete directions of atomic bonds such that each group is represented by an orientational component. Meanwhile, another kind of component, the volumetric one is derived from embedding energy. Damage and fracture are the changing and breaking of atomic bonds at the most fundamental level and have been reflected by the changing of these components’ properties. Therefore, material is treated as a component assembly, and its constitutive equations are formed by means of assembling these two kinds of components’ response functions. This material model is referred to as the component assembling model. Theoretical analysis and numerical computing indicate that the proposed model has the capacity of reproducing some results satisfactorily, with the advantages of physical explicitness and intrinsic induced anisotropy, etc.

Study Of The Damage-Induced Anisotropy Of Quasi-Brittle Materials Using The Component Assembling Model

Damage-induced anisotropy of quasi-brittle materials is investigated using component assembling model in this study. Damage-induced anisotropy is one significant character of quasi-brittle materials coupled with nonlinearity and strain softening. Formulation of such complicated phenomena is a difficult problem till now. The present model is based on the component assembling concept, where constitutive equations of materials are formed by means of assembling two kinds of components' response functions. These two kinds of components, orientational and volumetric ones, are abstracted based on pair-functional potentials and the Cauchy - Born rule. Moreover, macroscopic damage of quasi-brittle materials can be reflected by stiffness changing of orientational components, which represent grouped atomic bonds along discrete directions. Simultaneously, anisotropic characters are captured by the naturally directional property of the orientational component. Initial damage surface in the axial-shear stress space is calculated and analyzed. Furthermore, the anisotropic quasi-brittle damage behaviors of concrete under uniaxial, proportional, and nonproportional combined loading are analyzed to elucidate the utility and limitations of the present damage model. The numerical results show good agreement with the experimental data and predicted results of the classical anisotropic damage models.

Energy Principle And Nonlinear Electric-Mechanical Behavior Of Ferroelectric Ceramics

Many experimental observations have shown that a single domain in a ferroelectric material switches by progressive movement of domain walls, driven by a combination of electric field and stress. The mechanism of the domain switch involves the following steps: initially, the domain has a uniform spontaneous polarization; new domains with the reverse polarization direction nucleate, mainly at the surface, and grow though the crystal thickness; the new domain expands sideways as a new domain continues to form; finally, the domain switch coalesces to complete the polarization reversal. According to this mechanism, the volume fraction of the domain switching is introduced in the constitutive law of the ferroelectric material and used to study the nonlinear constitutive behavior of a ferroelectric body in this paper. The principle of stationary total potential 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. The mechanical field equation and a new domain switching criterion are obtained from the principle of stationary total potential energy. The domain switching criterion proposed in this paper is an expansion and development of the energy criterion established by Hwang et al. [ 1]. Based on the domain switching criterion, a set of linear algebraic equations for determining 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. If the volume fraction rho(I) of domain switching for each domain is prescribed, the unknown displacement and electric potential can be obtained based on the conventional finite element procedure. It is assumed that a domain switches if the reduction in potential energy exceeds a critical energy barrier. According to the experimental results, the energy barrier will strengthen when the volume fraction of the domain switching increases. The external mechanical and electric loads are increased step by step. The volume fraction rho(I) of domain switching for each element obtained from the last loading step is used as input to the constitutive equations. Then the strain and electric fields are calculated based on the conventional finite element procedure. The finite element analysis is carried out on the specimens subjected to uniaxial coupling stress and electric field. Numerical results and available experimental data are compared and discussed. The present theoretic prediction agrees reasonably with the experimental results.

国内非牛顿流体力学进展

非牛顿流体;;本构方程;;流变测量;;流动稳定性;;收缩流动;;聚合物加工;;石油工业;;生物流体力学

Plastic stress-strain fields and blunting effects during large deformations

Plastic stress-strain fields of two types of steel specimens loaded to large deformations are studied. Computational results demonstrate that, owing to the fact that the hardening exponent of the material varies as strain enlarges and the blunting of the crack tip, the well known HRR stress field in the plane strain model can only hold for the stage of a small plastic strain. Plastic dilatancy is shown to have substantial effects on strain distributions and blunting. To justify the constitutive equations used for analysis and to check the precision of computations, the load-deflection of a three-point bend beam and the load-elongation of an axisymmetric bar notched by a V-shaped cut were tested and recorded. The computed curves are in good accordance with experimental data.

Study of the damage-induced anisotropy of quasi-brittle materials using the component assembling model

Damage-induced anisotropy of quasi-brittle materials is investigated using component assembling model in this study. Damage-induced anisotropy is one significant character of quasi-brittle materials coupled with nonlinearity and strain softening. Formulation of such complicated phenomena is a difficult problem till now. The present model is based on the component assembling concept, where constitutive equations of materials are formed by means of assembling two kinds of components' response functions. These two kinds of components, orientational and volumetric ones, are abstracted based on pair-functional potentials and the Cauchy - Born rule. Moreover, macroscopic damage of quasi-brittle materials can be reflected by stiffness changing of orientational components, which represent grouped atomic bonds along discrete directions. Simultaneously, anisotropic characters are captured by the naturally directional property of the orientational component. Initial damage surface in the axial-shear stress space is calculated and analyzed. Furthermore, the anisotropic quasi-brittle damage behaviors of concrete under uniaxial, proportional, and nonproportional combined loading are analyzed to elucidate the utility and limitations of the present damage model. The numerical results show good agreement with the experimental data and predicted results of the classical anisotropic damage models.

Adiabatic shear localization: occurrence, theories, and applications

The elasto-damage theory of the components assembling model

The potential energy in materials is well approximated by pair functional which is composed of pair potentials and embedding energy. During calculating material potential energy, the orientational component and the volumetric component are derived respectively from pair potentials and embedding energy. The sum of energy of all these two kinds of components is the material potential. No matter how microstructures change, damage or fracture, at the most level, they are all the changing and breaking atomic bonds. As an abstract of atomic bonds, these components change their stiffness during damaging. Material constitutive equations have been formulated by means of assembling all components' response functions. This material model is called the component assembling model. Theoretical analysis and numerical computing indicate that the proposed model has the capacity of reproducing some results satisfactorily, with the advantages of great conceptual simplicity, physical explicitness, and intrinsic induced anisotropy, etc.

Numerical simulation of a single fluid particle rising in non-Newtonian fluids

In this work, a level set method is developed for simulating the motion of a fluid particle rising in non-Newtonian fluids described by generalized Newtonian as well as viscoelastic model fluids. As the shear-thinning model we use a Carreau-Yasuda model, and the viscoelastic effect can be modeled with Oldroyd-B constitutive equations. The control volume formulation with the SIMPLEC algorithm incorporated is used to solve the governing equations on a staggered Eulerian grid. The level set method is implemented to compute the motion of a bubble in a Newtonian fluid as one of typical examples for validation, and the computational results are in good agreement with the reported experimental data．The level set method is also applied for simulating a Newtonian drop rising in Carreau-Yasuda and Oldroyd-B fluids．Numerical results including noticeably negative wake behind the drop and viscosity field are obtained, and compare satisfactorily with the known literature data.

Effect of annealing on the viscoelastic and viscoplastic responses of low-density polyethylene

Two series of tensile tests with constant crosshead speeds (ranging from 5 to 200 mm/min) and tensile relaxation tests (at strains from 0.03 to 0.09) were performed on low-density polyethylene in the subyield region of deformations at room temperature. Mechanical tests were carried out on nonannealed specimens and on samples annealed for 24 h at the temperatures T = 50, 60, 70, 80, and 100 degreesC. Constitutive equations were derived for the time-dependent response of semicrystalline polymers at isothermal deformations with small strains. A polymer is treated as an equivalent heterogeneous network of chains bridged by temporary junctions (entanglements, physical crosslinks, and lamellar blocks). The network is thought of as an ensemble of mesoregions linked with each other. The viscoelastic behavior of a polymer is modeled as a thermally induced rearrangement of strands (separation of active strands from temporary junctions and merging of dangling strands with the network). The viscoplastic response reflects sliding of junctions in the network with respect to their reference positions driven by macrostrains. Stress-strain relations involve five material constants that were found by fitting the observations.

The viscoelastic and viscoplastic behavior of low-density polyethylene

Three series of tensile tests with constant cross-head speeds (ranging from 5 to 200 mm/min), tensile relaxation tests (at strains from 0.03 to 0.09) and tensile creep tests (at stresses from 2.0 to 6.0 MPa) are performed on low-density polyethylene at room temperature. Constitutive equations are derived for the time-dependent response of semicrystalline polymers at isothermal deformation with small strains. A polymer is treated as an equivalent heterogeneous network of chains bridged by temporary junctions (entanglements, physical cross-links and lamellar blocks). The network is thought of as an ensemble of meso-regions linked with each other. The viscoelastic behavior of a polymer is modelled as thermally-induced rearrangement of strands (separation of active strands from temporary junctions and merging of dangling strands with the network). The viscoplastic response reflects mutual displacement of meso-domains driven by macro-strains. Stress-strain relations for uniaxial deformation are developed by using the laws of thermodynamics. The governing equations involve five material constants that are found by fitting the observations. Fair agreement is demonstrated between the experimental data and the results of numerical simulation.

喜马拉雅西构造结区岩石圈演化三维有限元动力学数值模拟

The past three decades have seen numerous attempts to numerically model stress and strain patterns in the lithosphere of the Earth on both global and regional scales. This efforts have been indispensable in identifying the features we need to include in our endeavour to develop better models of our planet’s lithosphere and they have also raised our awareness for the many unresolved issue in the deep geodynamical issues that need to be addressed in the future. Nonetheless, in most models, the lithosphere is treated as a single layer with depth-averaged properties, and as the same distribution in the stress and strain fields, and as deforming under plane strain. All these above make a great hander for its reality and degree of recognition. As the beginning in this paper, some principal numerical models and results on the evolution of Tibetan plateau are reviewed and analyzed. Then, the geological and geophysical expedition on the Western Himalayan Syntaxis is briefly reviewed. Furthermore, we analysis the feature in deep geophysical field studies in this area and adjacent regions. Because, for most continents, stress models driven by plate boundary forces have successfully reproduced the main characteristics of the stress and strain field, we present a set of three-dimensional models of lithosphere system for a simplified geometry of the Western Himalayan Syntaxis area and its adjacent regions, where we try to match the first-order characteristics of the stress and strain fields of lithosphere since 10 Ma, and deformation and geodynamical evolution process in former 2Ma. Of course, the kinematic boundary conditions of the stress models driven by plate boundary forces were applied. The rheology plays a significant role in the lithospheric tectonics, which lead to different rheological parameters were used in different works although the have the same constitutive equations in models. So, in this paper we do not aim to produce all characteristics of the Western Himalayan Syntaxis areas’ stress and strain fields by the choices of various parameters, but rather the dynamic response between various rheological parameters and stress and strain fields. We have chosen to concentrate on the importance of rheology and lateral strength variations for lithospheric stress and strain patterns and use our findings to build a model of the Western Himalayan Syntaxis areas. In doing so, we want to go beyond purely elastic models or purely viscoelastic models. Compared the results of the crust viscosity in the Western Himalayan Syntaxis areas, we believed that, when various viscoelastic models are adopted, the selection of the coefficient of viscosity in the Western Syntaxis area has important influence on the its uplifts and evolutions. A wider uplift ranges and gently elevation was observed at the same time when a lower viscosity was used in our models, and vice versa. Data of stress magnitudes are not available, but it is clear that the stress levels must be at or below the failure threshold of rock under compression. Under these criteria, the calculation results show that the viscosity in the Western Syntaxis area should be smaller than 1023Pa.s When elastic model is adopted in relatively rigid Tarim basin, obvious changes are induced to the stress and strain fields of the whole Western Syntaxis area. We found that rigid block of lithosphere reduced stress levels within its interior and that, at the edges of such regions, stress orientation can change. Furthermore there is no evidence that such rigid regions act as stress barriers in that they shield areas in opposite sides of the structure from the influence of one another. In our models, the upper crustal material of the Western Syntaxis area does not turns to move westward. Whereas, because of the stress and strain fields have been decoupling at the interior of the lithosphere, we can get the results that the deep material must not move westward.