Longitudinal force on an embedded filament

Summary A thin elastic filament embedded in an elastic medium is subjected to a concentrated longitudinal load. For two-dimensional geometry and a relatively stiff filament the application of the load gives rise to a system of wedge-like and cylindrical waves. The dynamic shear stresses at the interface of the filament and the matrix, and in the region of the cylindrical waves are determined by means of Fourier transform techniques and Cagniard's method [2]. At the wave fronts of the wedge-like waves the jumps in the shear stresses are computed. Along the filament, the magnitudes of propagating discontinuities decrease exponentially. Along rays normal to the wave fronts of the wedge-like waves, the magnitudes of propagating discontinuities remain unchanged.     

On the stress field singularity in an incompressible half-space weakened by two near-surface wedge-like cracks

We consider the equilibrium problem for an elastic incompressible half-space weakened by two near-surface wedge-like cracks, whose lie in the same plane perpendicular to the half-space surface and have a common vertex. We use the Papkovich-Neuber representation to reduce the problem to finding two harmonic functions satisfying the mixed boundary conditions. These functions are constructed in spherical coordinates by using a Mehler-Fock type integral representation in Legendre functions. The analytic solution thus obtained permits finding the character of the stress distribution near the common tip of the cracks.     

SIMULATIONS OF MECHANICAL BEHAVIOR OF POLYCRYSTALLINE COPPER WITH NANO-TWINS

Mechanical behavior and microstructure evolution of polycrystalline copper with nano-twins were investigated in the present work by finite element simulations. The fracture of grain boundaries are described by a cohesive interface constitutive model based on the strain gradient plasticity theory. A systematic study of the strength and ductility for different grain sizes and twin lamellae distributions is performed. The results show that the material strength and ductility strongly depend on the grain size and the distribution of twin lamellae microstructures in the polycrystalline copper.     

Nucleation and Growth Behavior of Twin Region Around Yield Point of Polycrystalline Pure Ti

The aim of the present study is to investigate the nucleation and growth behavior of twin region around yield point of polycrystalline pure Ti under deformation. Firstly, we prepare commercial polycrystalline pure Ti plate, and investigate the microstructure and pole figures using an Electron Backscatter Diffraction Patterns device. Secondly, tensile specimens are cut out from 0°, 30°, 45° and 90° relative to plate rolling direction. Then, we measure the macroscopic stress–strain curve, local strain distribution and nucleation and growth of twin region arising in specimens under uniaxial tensile loading. Results show the anisotropic characteristics in those behaviors. Those could be related to c axis in hcp lattice. However, detailed anisotropic mechanism may have something to do with several interactions between slips and twins arising in its body. It is also understood that the avalanche behavior of twin region nucleation occurs as a result of larger twin region formation, with inhomogeneous small twin region nucleation in transient process. Finally, we could suppose the bridge mechanism of deformation behaviors from macroscale to microscale for polycrystalline pure Ti under deformation.     

A crystal plasticity theory for latent hardening by glide twinning through dislocation transmutation and twin accommodation effects

Imagine a residual glide twin interface advancing in a grain under the action of a monotonic stress. Close to the grain boundary, the shape change caused by the twin is partly accommodated by kinks and partly by slip emissions in the parent; the process is known as accommodation effects. When reached by the twin interface, slip dislocations in the parent undergo twinning shear. The twinning shear extracts from the parent dislocation a twinning disconnection, and thereby releases a transmuted dislocation in the twin. Transmutation populates the twin with dislocations of diverse modes. If the twin deforms by double twinning, double-transmutation occurs even if the twin retwins by the same mode or detwins by a stress reversal. If the twin deforms only by slip, transmutation is single. Whether single or double, dislocation transmutation is irreversible. The multiplicity of dislocation modes increases upon strain, since the twin finds more dislocations to transmute upon further slip of the parent and further growth of the twin. Thus, the process induces an increasing latent hardening rate in the twin. Under profuse twinning conditions, typical of double-lattice structures, this rate-increasing latent hardening combined with crystal rotation to hard orientations by twinning is consistent with a regime of increasing hardening rate, known as Regime II or Regime B. In this paper, we formulate governing equation of the above transmutation and accommodation effects in a crystal plasticity framework. We use the dislocation density based model originally proposed by Beyerlein and Tomé (2008) to derive the effect of latent hardening in a transmuting twin. The theory is expected to contribute to surmounting the difficulty that current models have to simultaneously predict under profuse twinning, the stress-strain curves, intermediate deformation textures, and intermediate twin volume fractions.     

3-D simulation of spatial stress distribution in an AZ31 Mg alloy sheet under in-plane compression

A complete 3-D crystal plasticity finite element method (CPFEM) that considered both crystallographic slip and deformation twinning was applied to simulate the spatial distribution of the relative amount of slip and twin activities in a polycrystalline AZ31 Mg alloy during in-plane compression. A microstructure mapping technique that considered the grain size distribution and microtexture measured by electron backscatter diffraction (EBSD) technique was used to create a statistically representative 3-D microstructure for the initial configuration. Using a 3-D Monte Carlo method, a 3-D digital microstructure that matched the experimentally measured grain size distribution was constructed. Crystallographic orientations obtained from the EBSD data were assigned on the 3-D digital microstructure to match the experimentally measured misorientation distribution. CPFEM captured the heterogeneity of the stress concentration as well as the slip and twin activities of a polycrystalline AZ31 Mg alloy during in-plane compression.     

Size effects on the martensitic phase transformation of NiTi nanograins

The analysis of nanocrystalline NiTi by transmission electron microscopy (TEM) shows that the martensitic transformation proceeds by the formation of atomic-scale twins. Grains of a size less than about 50 nm do not transform to martensite even upon large undercooling. A systematic investigation of these phenomena was carried out elucidating the influence of the grain size on the energy barrier of the transformation. Based on the experiment, nanograins were modeled as spherical inclusions containing (0 0 1) compound twinned martensite. Decomposition of the transformation strains of the inclusions into a shear eigenstrain and a normal eigenstrain facilitates the analytical calculation of shear and normal strain energies in dependence of grain size, twin layer width and elastic properties. Stresses were computed analytically for special cases, otherwise numerically. The shear stresses that alternate from twin layer to twin layer are concentrated at the grain boundaries causing a contribution to the strain energy scaling with the surface area of the inclusion, whereas the strain energy induced by the normal components of the transformation strain and the temperature dependent chemical free energy scale with the volume of the inclusion. In the nanograins these different energy contributions were calculated which allow to predict a critical grain size below which the martensitic transformation becomes unlikely. Finally, the experimental result of the atomic-scale twinning can be explained by analytical calculations that account for the transformation-opposing contributions of the shear strain and the twin boundary energy of the twin-banded morphology of martensitic nanograins.     

A micromechanical continuum model for the tensile behavior of shape memory metal nanowires

We have previously discovered a novel shape memory effect and pseudoelastic behavior in single-crystalline face-centered-cubic metal (Cu, Ni, and Au) nanowires. Under tensile loading and unloading, these wires can undergo recoverable elongations of up to 50%, well beyond the recoverable strains of 5-8% typical for most bulk shape memory alloys. This phenomenon only exists at the nanoscale and is associated with a reversible lattice reorientation driven by the high surface-stress-induced internal stresses. We present here a micromechanical continuum model for the unique tensile behavior of these nanowires. Based on the first law of thermodynamics, this model decomposes the lattice reorientation process into two parts: a reversible, smooth transition between a series of phase-equilibrium states and a superimposed irreversible, dissipative twin boundary propagation process. The reversible part is modeled within the framework of strain energy functions with multiple local minima. The irreversible, dissipative nature of the twin boundary propagation is due to the ruggedness of strain energy curves associated with dislocation nucleation, glide, and annihilation. The model captures the major characteristics of the unique behavior due to lattice reorientation and accounts for the size and temperature effects, yielding results that are in excellent agreement with the results of molecular dynamics simulations.     

A continuum thermodynamics formulation for micro-magneto-mechanics with applications to ferromagnetic shape memory alloys

A continuum thermodynamics formulation for micromagnetics coupled with mechanics is devised to model the evolution of magnetic domain and martensite twin structures in ferromagnetic shape memory alloys. The theory falls into the class of phase-field or diffuse-interface modeling approaches. In addition to the standard mechanical and magnetic balance laws, two sets of micro-forces and their associated balance laws are postulated; one set for the magnetization order parameter and one set for the martensite order parameter. Next, the second law of thermodynamics is analyzed to identify the appropriate material constitutive relationships. The proposed formulation does not constrain the magnitude of the magnetization to be constant, allowing for spontaneous magnetization changes associated with strain and temperature. The equations governing the evolution of the magnetization are shown to reduce to the commonly accepted Landau-Lifshitz-Gilbert equations for the case where the magnetization magnitude is constant. Furthermore, the analysis demonstrates that under certain limiting conditions, the equations governing the evolution of the martensite-free strain are shown to be equivalent to a hyperelastic strain gradient theory. Finally, numerical solutions are presented to investigate the fundamental interactions between the magnetic domain wall and the martensite twin boundary in ferromagnetic shape memory alloys. These calculations determine under what conditions the magnetic domain wall and the martensite twin boundary can be dissociated, resulting in a limit to the actuating strength of the material.     

On the modeling of equilibrium twin interfaces in a single-crystalline magnetic shape-memory alloy sample—III: Magneto-mechanical behaviors

This is the concluding part III of a series of papers. The aim of the current paper is to simulate and analyze the procedure of variant reorientation in a magnetic shape-memory alloy (MSMA) sample and to predict the response of the sample subject to various loading conditions. The sample to be considered in this paper has a 3D cuboid shape and is subject to typical magneto-mechanical loading conditions. Variant reorientation in the sample is realized through twin interface movements. To investigate the key features of twin interface movements, the properties of configurational forces on the twin interfaces are analyzed. For both the stress-assisted MFIS tests and the field-assisted quasi-elasticity tests, the magneto-mechanical behavior of the MSMA sample during the whole loading procedure is simulated by using the finite element method. The influence of the initial variant distribution in the sample on its global response is discussed. The obtained numerical results are compared with the experimental results. It can be seen that the model predictions can fit the experimental results both at a qualitative as well as at a quantitative level.     

Modelling the behaviour of polycrystalline austenitic steel with twinning-induced plasticity effect

A micromechanical model using the scale transition method in elastoviscoplasticity has been developed to describe the behaviour of those austenitic steels that display a TWIP effect. A physically based constitutive equation at the grain scale is proposed considering two inelastic strain modes: crystallographic slip and twinning. The typical organizations of microtwins observed in electron microscopy are considered, and the twin–slip as well as the twin–twin interactions are accounted for. The parameters for slip are first fitted on the uniaxial tensile response obtained at intermediate temperatures (when twinning is inhibited). Then, the parameters associated with twinning are identified using the stress–strain curve at room temperature. The simulated results in both macro and micro scales are in good agreement with experimentally obtained results.     

Critical thickness for misfit twinning in an epilayer

A Somigliana dislocation dipole model is developed to determine the critical thickness for misfit twin formation in an epilayer with different elastic constants from its substrate. The critical dipole arm length is determined by minimizing the twin formation energy for a given epilayer thickness and lattice mismatch strain, while a zero value of the minimum formation energy determines the critical thickness for misfit twinning. The results obtained by the Somigliana dislocation dipole model are roughly consistent with those by the previous dislocation-based twinning model.     

On the modeling of equilibrium twin interfaces in a single-crystalline magnetic shape memory alloy sample. I: theoretical formulation

In this paper, a variational approach is proposed to study the response of a single-crystalline magnetic shape memory alloy (MSMA) sample subject to external forces and magnetic fields. Especially, some criteria are derived to model the (quasi-static) movements of twin interfaces in the sample. By considering the compatibility condition, twin interfaces between two martensite variants are found to be flat planes with given normal vectors. To adopt the variational method, a total energy functional for the whole magneto-mechanical system is proposed. By calculating the variations of the total energy functional with respect to the independent variables, the equilibrium equations and the evolution laws for the internal variables can be derived. By further considering the variation of the total energy functional with respect to the variant distribution, some criteria for twin interface movements can be derived. The governing system of the current model is then formulated by composing the equilibrium equations, the evolution laws for the internal variables and the twin interface movement criteria. To show the validity of the governing system, some analytical results are constructed under certain simplified conditions, which can be used to simulate the magneto-mechanical response of the MSMA sample.     

The effect of time-dependent twinning on low temperature (<0.25 ∗ Tm) creep of an alpha-titanium alloy

The low-temperature (less than one-fourth of the melting temperature) creep deformation behavior of hexagonally close-packed (HCP) α-Ti–1.6 wt.% V was investigated. Creep tests were performed at various temperatures between room temperature and 205 °C at 95% of the respective yield stress at the different temperatures. The creep strain rate was found to increase with increasing temperature. Scanning and transmission electron microscopy revealed that slip and unusually slow twin growth, or time-dependent twinning, are active deformation mechanisms for the entire temperature range of this investigation. The activation energy for creep of this alloy was calculated to identify the rate-controlling deformation mechanism, and was found to increase with increasing creep strain. At low strain, the activation energy for creep was found to be close to the previously calculated activation energy for slip. At high strain, the calculated activation energy indicates that both slip and twinning are significant deformation mechanisms. The appearance of twinning at high strains is explained by a model for twin nucleation by dislocation pileups.     

Characterization of Single-Phase Flow Through Carbonate Rocks: Quantitative Comparison of NMR Flow Propagator Measurements with a Realistic Pore Network Model

We report here the quantitative comparisons between the measured NMR flow propagator of a carbonate rock and the flow propagator calculated with a porous network extracted from the micro-CT image of the twin plug. We developed a numerical model based on a particle tracking algorithm in pore space. The particle tracking in throats is described using the first arrival time distribution. As pores have an important volume fraction in the sample considered, we implemented a time-delay mechanism for particle transport in the pores. We consider that the nodes have volume and there is a transport of the tracking particles inside the nodes, which leads to an “apparent” time-delay. Simulations of flow propagator show good agreement with low field NMR experiments performed on the twin plug of the sample used for pore network extraction with a single adjustable parameter (that describes the dynamics in the pores). These results lead us to a better understanding of the connection between pore structure and the behavior of NMR flow propagator in fluid-saturated rocks and are essential in interpreting the experimental data and correlating NMR parameters to petrophysical properties.     

Edge crack opening in an elastic plane with a wedge-like cut in contact with a punch

The equilibrium of an elastic plane with a wedge-like cut and an internal or edge crack on the symmetry axis was studied in [1] in the case of punch indentation in the lateral faces of the cut at a distance from the cut tip. In [1], the systemof singular integral equations of the problemwas solved numerically by the mechanical quadrature method. In this paper, the generalized Wiener-Hopf method [2] is used to obtain the analytic solution of a similar problem in the case of an edge crack under punch pressure on parts of the cut lateral faces adjacent to the cut tip. Some special cases of this problem were considered earlier without a crack [3, 4] or a punch [5, 6].     

孪生诱发塑性对V-5Cr-5Ti合金应变硬化行为的影响

钒合金(V-Cr-Ti)作为潜在重要的聚变反应堆用结构材料, 近年来受到广泛的关注. 为了研究 V-5Cr-5Ti 合金不同应变率压缩下的应变硬化行为, 特别是孪生对塑性变形的影响, 以位错密度和孪晶演化为基础, 建立了该合金的应变硬化模型. 模型中考虑了孪晶中的位错滑移对材料塑性应变的贡献. 模拟结果表明, 由于孪生诱发塑性, 从而使动态压缩时的位错密度小于准静态加载时的, 这使得 V-5Cr-5Ti 合金在动态压缩时的应变硬化率比准静态加载时的小. 当孪晶形成后, 位错滑移引起的塑性应变率随应变增大而增大, 并逐渐接近加载应变率, 而孪生引起的塑性应变率则随应变增大而减小.      

Development of inhomogeneous fields under postcritical deformation of steel specimens in extension

The paper deals with experimental studies of inhomogeneous strain fields in the “neck” region in extension of plane and cylindrical specimens at the postcritical stage, which directly precedes the fracture, by using a video system and a digital image correlation method. We consider the problems of interpretation of extension diagrams obtained for specimens of various length with strain inhomogeneity in the working section of the specimen with a “neck.”We obtain experimental data about the distribution of longitudinal, transverse, and shear strain fields and the displacement fields in the region of the plane specimen localization by using the digital image correlation method.