Compatibility conditions for a left Cauchy-Green strain field

This paper deals chiefly with various issues pertaining to the existence and uniqueness of a finite deformation that gives rise to a prescribed right or left Cauchy-Green strain-tensor field.Following a review and discussion of available existence and uniqueness theorems appropriate to a pre-assigned right strain field, the extent of uniqueness of a generating deformation is established under minimal smoothness and invertibility assumptions.The remainder of the paper is devoted to the more involved corresponding existence and uniqueness questions for a given left strain-tensor field. These questions are first discussed in a three-dimensional setting and are then resolved for the special class of plane deformations. The results thus obtained stand in marked contrast to their counterparts for a given right strain field.     

Integral Solution of a Problem with Non-Standard Boundary Conditions in an Enhanced Theory of Bending of Elastic Plates

A modified boundary integral equation method is used to solve a specific type of mixed boundary value problem in an enhanced theory of bending of elastic plates in which the effects of transverse shear deformation and transverse normal strain are taken into account. The problem considered is characterized by the fact that a combination of transverse displacement and bending and twisting moments is prescribed on the curve which bounds the middle surface of the plate. Both interior and exterior problems are formulated and the corresponding existence and uniqueness results derived.     

Mode I and mode II crack tip asymptotic fields with strain gradient effects

The strain gradient effect becomes significant when the size of fracture process zone around a crack tip is comparable to the intrinsic material lengthl, typically of the order of microns. Using the new strain gradient deformation theory given by Chen and Wang, the asymptotic fields near a crack tip in an elastic-plastic material with strain gradient effects are investigated. It is established that the dominant strain field is irrotational. For mode I plane stress crack tip asymptotic field, the stress asymptotic field and the couple stress asymptotic field can not exist simultaneously. In the stress dominated asymptotic field, the angular distributions of stresses are consistent with the classical plane stress HRR field; In the couple stress dominated asymptotic field, the angular distributions of couple stresses are consistent with that obtained by Huang et al. For mode II plane stress and plane strain crack tip asymptotic fields, only the stress-dominated asymptotic fields exist. The couple stress asymptotic field is less singular than the stress asymptotic fields. The stress asymptotic fields are the same as mode II plane stress and plane strain HRR fields, respectively. The increase in stresses is not observed in strain gradient plasticity for mode I and mode II, because the present theory is based only on the rotational gradient of deformation and the crack tip asymptotic fields are irrotational and dominated by the stretching gradient. The project supported by the National Natural Science Foundation of China (19704100), National Natural Science Foundation of Chinese Academy of Sciences (KJ951-1-20), CAS K.C. Wong Post-doctoral Research Award Fund and Post-doctoral Science Fund of China     

Multi-Scale Asymptotic Expansion for a Small Inclusion in Elastic Media

The aim of this paper is to present an asymptotic expansion of the influence of a small inclusion of different stiffness in an elastic media. The applicative interest of this study is to provide tools which take into account this influence and correct the deformation without inclusion by additive terms that can be precalculated and which depend only on the shape of the inclusion. We treat two problems: an anti-plane linearized elasticity problem and a plane strain problem. On every expansion order we provide corrective terms modeling the influence of the inclusion using techniques of scaling and multi-scale asymptotic expansions. The resulting expansion is validated by comparing it to a test case obtained by solving the Poisson transmission problem in the case of an inclusion of circular shape using the separation of variables method. Proofs of existence and uniqueness of our fields on unbounded domains are also adapted to the bidimensional Poisson problem and the linear elasticity problem.     

韧性金属裂纹尖端损伤区内应变场测量——数字散斑相关法应用

本文用数字散斑相关方法测量了五种不同幂硬化指数韧性金属材料(铝和铜),双边裂纹尖端细观区域内应变场.对所得结果用韧性损伤模型进行了分析.在此法中以金属自然表面结构为散斑场,不同加载状态的散斑场进行比较,得到相对变形与应变.图象之间相关性 C 是变形参数或是位移及其导数的泛函.使其相关性 C 取最大值的试凑变形即为其真实变形场.这一方法在细观测量中应用得到满意的结果.     

Wedge indentation into elastic-plastic single crystals, 1: Asymptotic fields for nearly-flat wedge

Asymptotic stress and deformation fields under the contact point singularities of a nearly-flat wedge indenter and of a flat punch are derived for elastic ideally-plastic single crystals with three effective in-plane slip systems that admit a plane strain deformation state. Face-centered cubic (FCC), body-centered cubic (BCC), and hexagonal-close packed (HCP) crystals are considered. The asymptotic fields for the flat punch are analogous to those at the tip of a stationary crack, so a potential solution is that the deformation field consists entirely of angular constant stress plastic sectors separated by rays of plastic deformation across which stresses change discontinuously. The asymptotic fields for a nearly-flat wedge indenter are analogous to those of a quasistatically growing crack tip fields in that stress discontinuities can not exist across sector boundaries. Hence, the asymptotic fields under the contact point singularities of a nearly-flat wedge indenter are significantly different than those under a flat punch. A family of solutions is derived that consists entirely of elastically deforming angular sectors separated by rays of plastic deformation across which the stress state is continuous. Such a solution can be found for FCC and BCC crystals, but it is shown that the asymptotic fields for HCP crystals must include at least one angular constant stress plastic sector. The structure of such fields is important because they play a significant role in the establishment of the overall fields under a wedge indenter in a single crystal. Numerical simulations—discussed in detail in a companion paper—of the stress and deformation fields under the contact point singularity of a wedge indenter for a FCC crystal possess the salient features of the analytical solution.     

Wedge indentation into elastic-plastic single crystals. 2: Simulations for face-centered cubic crystals

The asymptotic stress and deformation fields associated with the contact point singularity of a nearly-flat wedge indenter impinging on a specially-oriented single face-centered cubic crystal are derived analytically in a companion paper. In order to investigate the extent of the asymptotic fields, the indentation process is simulated numerically using single crystal plasticity. The quasistatically translating asymptotic fields consist of four angular elastic sectors separated by plastically deforming sector boundaries, as predicted in the companion paper. The asymptotic stress distributions are in accord with the analytical predictions. In addition, simulations are performed for a wedge indenter with a 90° included angle in order to investigate the consequences of finite deformation and finite lattice rotation. Several salient features of the deformation field for the nearly-flat indenter persist in the deformation field for the 90° wedge indenter. The existence of the salient features is validated by comparison to experimental measurements of the lower bound on geometrically necessary dislocation (GND) densities.     

Effect of material non-homogeneity on the inhomogeneous shearing deformation of a Gent slab subjected to a temperature gradient

It is well known that most rubber-like materials are non-homogeneous due to either imperfect manufacturing conditions or the action of severe thermo-oxidative environments in many practical applications. In this study, within the context of finite thermoelasticity, we theoretically analyze the inhomogeneous shearing deformation of a non-homogeneous rubber-like slab subjected to a thermal gradient across its thickness. The major objective of this study is to investigate the effect of the material non-homogeneity, which is the material-coordinate dependence of the material response functions, on the stress-strain fields for a given temperature gradient. First, we show the existence of a simple shearing deformation from which the generalized shear modulus and the generalized thermal conductivity of the slab could be obtained. Based on this information, the Gent material model is generalized to take the material non-homogeneity and the temperature dependence of the stress into account. To analyze the inhomogeneous shearing deformation of the non-homogeneous slab, deformation and temperature fields are postulated; then the decoupled temperature field is obtained analytically by solving the local energy balance equation. Finally, the static equilibrium equations are solved considering the linear temperature field. Our results show that the spatial pattern and the degree of the material non-homogeneity have profound effects on the stress-strain fields. The shear strain becomes nearly homogeneous and the stresses are relatively small for a certain spatial variation of the material non-homogeneity. This result suggests the possibility of designing a novel class of materials: functionally graded rubber-elastic materials (FGREMs).     

????????????????????????????????

在重轨万能轧制过程中,首先建立了简化的三维理论模型. 然后分别给出了轨腰、轨头及轨 底的运动学许可速度场以及相应的应变速度场和剪应变速度强度,并求出了相应变形区的塑 性变形功率、速度间断面上消耗的功率以及由于摩擦产生的摩擦功率. 最后根据上限原理分 别求解了水平辊和两个立辊的轧制力上限解. 通过比较可知,二辊轧制理论公式误差很大因 此不能用于万能轧制过程,而上限法求得的轧制力近似解大于轧制现场数据但最大误差不超 过13\%, 因此根据上限原理进行轧制工艺参数设定及优化是比较可靠的.     

Non-linear constitutive relations for magnetostrictive materials

In this paper, non-linear deformation behavior of magnetostrictive materials is studied and three magnetoelastic coupling constitutive models are developed. The standard square (SS) constitutive model is developed by means of truncating the polynomial expansion of the Gibbs free energy. The hyperbolic tangent (HT) constitutive equations, which involve a hyperbolic tangent magnetic-field dependence, are proposed to model the magnetic-field-induced strain saturation of magnetostrictive materials in the region of intense magnetic fields. A new model based on density of domain switching (DDS) is established in terms of the basic truth that magnetic domain switching underlies magnetostrictive deformation. In this model, it is assumed that the relation between density of domain switching, defined by the quantity of magnetic domains switched by per unit magnetic field and magnetic field can be described by a density function with normal distribution. The moduli in these constitutive models can be determined by a material function that is proposed to describe the dependence of the peak piezomagnetic coefficient on the compressive pre-stress for one-dimensional cases based on the experimental results published. The accuracy of the non-linear constitutive relations is evaluated by comparing the theoretical values with experimental results of a Terfenol-D rod operated under both compressive pre-stress and bias magnetic field. Results indicate that the SS constitutive equations can accurately predict the experimental results under a low or moderate magnetic field while the HT model can, to some extent, reflect the trend of saturation of magnetostrictive strain under a high magnetic field. The model based on DDS, which is more effective in simulating the experimental curves, can capture the main characteristics of the mechanism of magnetoelastic coupling deformation of a Terfenol-D rod, such as the notable dependence of magnetoelastic response on external stress and the saturation of magnetostrictive strain under intense magnetic fields. In addition, the SS constitutive relation for a general three-dimensional problem is discussed and an approach to characterize the modulus tensors is proposed.     

BOUNDARY ELEMENT METHOD FOR BUCKLING EIGENVALUE PROBLEM AND ITS CONVERGENCE ANALYSIS

ＩｎｔｒｏｄｕｃｔｉｏｎＴｈｅｂｕｃｋｌｉｎｇｅｉｇｅｎｖａｌｕｅｐｒｏｂｌｅｍｈａｓｉｍｐｏｒｔａｎｔｓｉｇｎｉｆｉｃａｎｃｅｉｎｔｈｅｓｔａｂｉｌｉｔｙａｎａｌｙｓｉｓｏｆｅｎｇｉｎｅｅｒｉｎｇｓｔｒｕｃｔｕｒｅ.Ｈｅｎｃｅｔｈｅｎｕｍｅｒｉｃａｌｃａｌｃｕｌａｔｉｏｎｆｏｒｔｈｅｓｅｐｒｏｂｌｅｍｓｉｓｅｘｔｒｅｍｅｌｙｍｅａｎｉｎｇｆｕｌｉｎｃｏｍｐｕｔａｔｉｏｎａｌｍｅｃｈａｎｉｃｓ.ＴｈｅｐｒｅｓｅｎｔｃｏｍｐｕｔａｔｉｏｎａｌｍｅｔｈｏｄｓｆｏｃｕｓｏｎＦＥＭ ,ｄｉｆｆｅｒｅｎｃｅｍ…     

Experimental investigation on surface deformation of soft half plane indented by rigid wedge

In this work, the mechanical behavior of a block of soft material subject to large deformation from a series of wedge-shaped indenters is evaluated. Data fields acquired from digital image correlation (DIC) are compared with the existing theoretical models. The slope angles of the wedges vary from 5° to 73.5°, and the minimum measurement uncertainties of the DIC system are established in advance to define the accuracy. It is concluded that the assumptions underpinning the analytical theory make it difficult to characterize large deformation of soft materials during contact. The strain fields are also obtained from the measured displacement field and verify the previously postulated existence of two deformation sectors, namely, a so-called shrinkage sector symmetric to the loading axis and an expansion sector, which become smaller with the increasing load and decreasing wedge angle.     

Effects of the strain-hardening exponent on two-parameter characterizations of surface-cracks under large-scale yielding

The influence of strain hardening exponent on two-parameter J-Q near tip opening stress field characterization with modified boundary layer formulation and the corresponding validity limits are explored in detail. Finite element simulations of surface cracked plates under uniaxial tension are implemented for loads exceeding net-section yield. The results from this study provide numerical methodology for limit analysis and demonstrate the strong material dependencies of fracture parameterization under large scale yielding. Sufficient strain hardening is shown to be necessary to maintain J-Q predicted fields when plastic flow progresses through the remaining ligament. Lower strain hardening amplifies constraint loss due to stress redistribution in the plastic zone and increases the ratio of tip deformation to J. The onset of plastic collapse is marked by shape change and/or rapid relaxation of tip fields compared to those predicted by MBL solutions and thus defining the limits of J-Q dominance. A radially independent Q-parameter cannot be evaluated for the low strain hardening material at larger deformations within a range where both cleavage and ductile fracture mechanisms are present. The geometric deformation limit of near tip stress field characterization is shown to be directly proportional to the level of stress the material is capable of carrying within the plastic zone. Accounting for the strain hardening of a material provides a more adjusted and less conservative limit methodology compared to those generalized by the yield strength alone. Results from this study are of relevance to establishing testing standards for surface cracked tensile geometries.     

Generating and Adding Flows on Locally Complete Metric Spaces

As a generalization of a vector field on a manifold, the notion of an arc field on a locally complete metric space was introduced in Bleecker and Calcaterra (J Math Anal Appl, 248: 645–677, 2000). In that paper, the authors proved an analogue of the Cauchy–Lipschitz Theorem, i.e they showed the existence and uniqueness of solution curves for a time independent arc field. In this paper, we extend the result to the time dependent case, namely we show the existence and uniqueness of solution curves for a time dependent arc field. We also introduce the notion of the sum of two time dependent arc fields and show existence and uniqueness of solution curves for this sum.     

Boundary-layer corrections for stress and strain fields in randomly heterogeneous materials

Boundary-layer effects on the effective response of fibre-reinforced media are analysed. The distribution of the fibres is assumed random. A methodology is presented for obtaining non-local effective constitutive operators in the vicinity of a boundary. These relate ensemble averaged stress to ensemble averaged strain. Operators are also developed which re-construct the local fields from their ensemble averages. These require information on the local configuration of the medium. Complete information is likely not to be available, but averages of these operators conditional upon any given local information generate corresponding conditional averages of the fields. Explicit implementation is performed within the framework of an approximation of Hashin-Shtrikman type. Two types of geometry are considered in examples: a half-space and a crack in an infinite heterogeneous medium. These are representative, asymptotically, of the field in the vicinity of any smooth boundary, and in the vicinity of a crack tip, respectively. Results have been obtained for the case of anti-plane deformation, realized by the imposition of either Dirichlet or Neumann conditions on the boundary; those for the Neumann condition are presented and discussed explicitly. The stresses in both fibre and matrix adjacent to a crack tip are shown to differ substantially from the values that would be predicted by ordinary homogenization.     

Comparison of Test Specimens for Characterizing the Dynamic Properties of Rubber

Dynamic material properties inferred via experiment can be strongly influenced by the choice of test specimen geometry unless care is taken to ensure that mechanical fields (stress, strain, etc.) within the specimen adequately reflect the ideal homogeneous deformation state. In this work, finite element models of simple shear, cylindrical compression, simple tension, and bi-conical shear test specimens were analyzed in order to quantify the relative conformity of each specimen to its corresponding ideal. Three metrics of conformity were evaluated, based respectively on the distributions of stress, strain, and strain energy density. The results show that a simple shear specimen provides superior conformity. Other factors involved in the selection of test specimen geometry are also discussed. Such factors include relative linearity and symmetry of measured stress–strain data, grip slip, and heat build up.     

Controllability,uniqueness and existence of the incremental response: A mathematical criterion for elastoplastic constitutive laws

In this paper the onset of instabilities in elastoplastic materials is theoretically studied and a conceptual basis for understanding the physical implications of a loss of uniqueness and/or existence of the incremental response is provided. For this purpose, the concept of test controllability is reinterpreted and mixed stress–strain loading programmes are accounted for. A set of scalar indices, the moduli of instability, related with the inception of an unstable response is introduced and their dependency on the loading programme is explicitly illustrated. The paper shows that the use of these newly defined scalar measures provides support for an alternative definition for mechanical stability, which is closely related with the mathematical notions of existence and uniqueness of the predicted incremental response. In the final section, some mathematical properties of the moduli of instability are discussed, suggesting a novel reinterpretation of other well established theories and providing additional tools for the future application of the proposed framework.     

A Dedicated DIC Methodology for Characterizing Plastic Deformation in Single Crystals

This paper focuses on the development of an appropriate Digital Image Correlation (DIC) methodology based on Image Registration and dedicated for characterizing the plastic deformation in single crystals. A pure nickel single crystal specimen is plastically deformed in tension and investigated by DIC technique. Based on the measured kinematic fields, the proposed method enables to identify the slip activity on the crystal surface and to locate precisely the slip band interfaces at microscale which behave as kinematic discontinuities. The computed displacement data are projected on a well-defined physical basis containing slip details, then the strain fields can be derived directly from a set of analytical functions. The possible errors in displacement induced by this projection approach are evaluated. Finally, some results of the evaluated strain fields are presented. It demonstrates that the developed DIC methodology allows quantitative characterization of a heterogeneous deformation process and promotes further relationships to be established between slip activity and strain field evolution in single crystals.     

On the Dirichlet problem for incompressible elastic materials

The linearized equations governing the deformations of incompressible elastic bodies are discussed. The Dirichlet problem is formulated for this system of equations using the theory of elliptic systems due to Douglis and Nirenberg. A uniqueness theorem is proved. Necessary and sufficient conditions for uniqueness of solution to the Dirichlet problem are obtained for small deformations of a Mooney-Rivlin material which has been subjected to a finite homogeneous biaxial deformation.     

Micropolar nematic model for polarized liquid crystals

A micromorphic electroelastic model for polarized liquid crystals is proposed on the basis of a representation of electric multipoles in terms of microdeformation. Nematic liquid crystals are modeled as micropolar continua endowed with intrinsic electric dipole and quadrupole. A nonlinear dimensionless problem for a homogeneous nematic layer is formulated and solved numerically. The existence of a threshold electric potential is discussed, and the corresponding linearized system is also obtained to compare results on small values of deformation and electric field. Differently from common results of the classical continuum approach, asymmetric deformations and electric potentials within the layer are obtained due to the occurrence of non-null intrinsic quadrupole.