A theory of plasticity for carbon nanotube reinforced composites

Carbon nanotubes (CNTs) possess exceptional mechanical properties, and when introduced into a metal matrix, it could significantly improve the elastic stiffness and plastic strength of the nanocomposite. But current processing techniques often lead to an agglomerated state for the CNTs, and the pristine CNT surface may not be able to fully transfer the load at the interface. These two conditions could have a significant impact on its strengthening capability. In this article we develop a two-scale micromechanical model to analyze the effect of CNT agglomeration and interface condition on the plastic strength of CNT/metal composites. The large scale involves the CNT-free matrix and the clustered CNT/matrix inclusions, and the small scale addresses the property of these clustered inclusions, each containing the randomly oriented, transversely isotropic CNTs and the matrix. In this development the concept of secant moduli and a field fluctuation technique have been adopted. The outcome is an explicit set of formulae that allows one to calculate the overall stress-strain relations of the CNT nanocomposite. It is shown that CNTs are indeed a very effective strengthening agent, but CNT agglomeration and imperfect interface condition can seriously reduce the effective stiffness and elastoplastic strength. The developed theory has also been applied to examine the size (diameter) effect of CNTs on the elastic and elastoplastic response of the composites, and it was found that, with a perfect interface contact, decreasing the CNT radius would enhance the overall stiffness and plastic strength, but with an imperfect interface the size effect is reversed. A comparison of the theory with some experiments on the CNT/Cu nanocomposite serves to verify the applicability of the theory, and it also points to the urgent need of eliminating all CNT agglomeration and improving the interface condition if the full potential of CNT reinforcement is to be realized.     

An Analysis of Professor J.F. Bell's Research on the Finite Twist and Extension of Thin-walled Polycrystalline Cylindrical Tubes

Professor J.F. Bell's empirical result regarding the rotation factor in the polar decomposition of the deformation gradient for the finite twist–extension of a thin-walled polycrystalline cylindrical metal tube is examined. The correct expression for the rotation is derived and used to show how Bell's result should be interpreted. Some implications for his incremental plasticity equations are also discussed. In particular, they are shown to satisfy appropriate invariance requirements when cast in terms of the variables actually measured by Bell in his experiments. Further consequences of his equations consistent with his data are also derived. Finally, it is shown that his theory furnishes a consistent constitutive statement about the response of isotropic solids provided that the Cauchy stress is constrained to be symmetric.     

Size effects and idealized dislocation microstructure at small scales: Predictions of a Phenomenological model of Mesoscopic Field Dislocation Mechanics: Part II

In Part I of this set of two papers, a model of mesoscopic plasticity is developed for studying initial-boundary value problems of small scale plasticity. Here we make qualitative, finite element method-based computational predictions of the theory. We demonstrate size effects and the development of strong inhomogeneity in simple shearing of plastically constrained grains. Non-locality in elastic straining leading to a strong Bauschinger effect is analyzed. Low shear strain boundary layers in constrained simple shearing of infinite layers of polycrystalline materials are not predicted by the model, and we justify the result based on an examination of the no-dislocation-flow boundary condition. The time-dependent, spatially homogeneous, simple shearing solution of PMFDM is studied numerically. The computational results and an analysis of continuous dependence with respect to initial data of solutions for a model linear problem point to the need for a nonlinear study of a stability transition of the homogeneous solution with decreasing grain size and increasing applied deformation. The continuous-dependence analysis also points to a possible mechanism for the development of spatial inhomogeneity in the initial stages of deformation in lower-order gradient plasticity theory. Results from thermal cycling of small scale beams/films with different degrees of constraint to plastic flow are presented showing size effects and reciprocal-film-thickness scaling of dislocation density boundary layer width. Qualitative similarities with results from discrete dislocation analyses are noted where possible.We discuss the convergence of approximate solutions with mesh refinement and its implications for the prediction of dislocation microstructure development, motivated by the notion of measure-valued solutions to conservation laws.     

Micromechanics, macromechanics and constitutive modeling of the elasto-viscoplastic deformation of rubber-toughened glassy polymers

Under certain conditions, such as sufficiently low temperatures, high loading rates and/or highly triaxial stress states, glassy polymers display an unfavorable characteristic—brittleness. A technique used for reducing the brittleness (increasing the fracture toughness) of these materials is rubber toughening. While there is significant qualitative understanding of the mechanical behavior of rubber-toughened polymers, quantitative modeling tools for the large-strain deformation of rubber-toughened glassy polymers are largely lacking.In this paper, we develop a suite of numerical tools to investigate the mechanical behavior of rubber-toughened glassy polymers, with emphasis on rubber-toughened polycarbonate. The rubber particles are modeled as voids in view of their deformation-induced cavitation early during deformation. A three-dimensional micromechanical model of the heterogeneous microstructure is developed to study the effects of initial rubber particle (void) volume fraction on the underlying elasto-viscoplastic deformation mechanisms in the material, and how these mechanisms influence the macroscopic response of the material. A continuum-level constitutive model is developed for the large-strain elasto-viscoplastic deformation of porous glassy polymers, and it is calibrated against micromechanical modeling results for porous polycarbonate. The constitutive model can be used to study various boundary value problems involving rubber-toughened (porous) glassy polymers. As an example, the case of an axisymmetric notched bar is simulated for the case of polycarbonate with varying levels of initial porosity. The quality of the constitutive model calibration is assessed using a multi-scale modeling approach.     

Finite element approximation of field dislocation mechanics

A tool for studying links between continuum plasticity and dislocation theory within a field framework is presented. A finite element implementation of the geometrically linear version of a recently proposed theory of field dislocation mechanics (J. Mech. Phys. Solids 49 (2001) 761; Proc. Roy. Soc. 459 (2003) 1343; J. Mech. Phys. Solids 52 (2004) 301) represents the main idea behind the tool. The constitutive ingredients of the theory under consideration are simply elasticity and a specification of dislocation velocity and nucleation. The set of equations to be approximated are non-standard in the context of solid mechanics applications. It comprises the standard second-order equilibrium equations, a first-order div-curl system for the elastic incompatibility, and a first-order, wave-propagative system for the evolution of dislocation density. The latter two sets of equations require special treatment as the standard Galerkin method is not adequate, and are solved utilizing a least-squares finite element strategy. The implementation is validated against analytical results of the classical elastic theory of dislocations and analytical results of the theory itself. Elastic stress fields of dislocation distributions in generally anisotropic media of finite extent, deviation from elastic response, yield-drop, and back-stress are shown to be natural consequences of the model. The development of inhomogeneity, from homogeneous initial conditions and boundary conditions corresponding to homogeneous deformation in conventional plasticity, is also demonstrated. To our knowledge, this work represents the first computational implementation of a theory of dislocation mechanics where no analytical results, singular solutions in particular, are required to formulate the implementation. In particular, a part of the work is the first finite element implementation of Kröner's linear elastic theory of continuously distributed dislocations in its full generality.     

基于土质边坡塑性极限分析条分法的可靠度计算方法研究

从土塑性极限分析上限定理出发,以土质边坡塑性极限分析条分法为基础,在引入可靠性分析方法(JC法)后,建立了土质边坡塑性极限分析条分法可靠度极限状态方程,并推导得到了塑性极限分析条分法可靠度计算公式。这种方法实现了定值计算向不确定性计算的转变,同时可靠度指标β的计算结果具有概率的意义。     

Nonlinear Dynamics and Stability of a Two D.O.F. Elastic/Elasto-Plastic Model System

The local and global nonlinear dynamics of a two-degree-of-freedom model system is studied. The undeflected model consists of an inverted T formed by three rigid bars, with the tips of the two horizontal bars supported on springs. The springs exhibit an elasto-plastic response, including the Bauschinger effect. The vertical rigid bar is subjected to a conservative (dead) or non-conservative (follower) force having static and periodic components. First, the method of multiple scales is used for the analysis of the local dynamics of the system with elastic springs. The attention is focused at modal interaction phenomena in weak excitation at primary resonance and in hard sub-harmonic excitation. Three different asymptotic expansions are utilised to get a structural response for typical ranges of excitation parameters. Numerical integration of the governing equations is then performed to validate results of asymptotic analysis in each case. A full global nonlinear dynamics analysis of the elasto-plastic system is performed to reveal the role of plastic deformations in the stability of this system. Static 'force-displacement' curves are plotted and the role of plastic deformations in the destabilisation of the system is discussed. Large-amplitude non-linear oscillations of the elasto-plastic system are studied, including the influence of material hardening and of static and sinusoidal components of the applied force. A practical method is proposed for the study of a non-conservative elasto-plastic system as a non-conservative elastic system with an 'equivalent' viscous damping. This revised version was published online in July 2006 with corrections to the Cover Date.     

Effect of material nonhomogeneities on the HRR dominance

This work studies the asymptotic stress and displacement fields near the tip of a stationary crack in an elastic–plastic nonhomogeneous material with the emphasis on the effect of material nonhomogeneities on the dominance of the crack tip field. While the HRR singular field still prevails near the crack tip if the material properties are continuous and piecewise continuously differentiable, a simple asymptotic analysis shows that the size of the HRR dominance zone decreases with increasing magnitude of material property gradients. The HRR field dominates at points that satisfy |α⁻¹ ∂α/∂x_δ|1/r, |α⁻¹ ∂²α/(∂x_δ ∂x_γ)|1/r², |n⁻¹ ∂n/∂x_δ|1/[r|ln(r/A)|] and |n⁻¹ ∂²n/(∂x_δ ∂x_γ)|1/[r²|ln(r/A)|], in addition to other general requirements for asymptotic solutions, where α is a material property in the Ramberg–Osgood model, n is the strain hardening exponent, r is the distance from the crack tip, x_δ are Cartesian coordinates, and A is a length parameter. For linear hardening materials, the crack tip field dominates at points that satisfy |E_tan⁻¹ ∂E_tan/∂x_δ|1/r, |E_tan⁻¹ ∂²E_tan/(∂x_δ ∂x_γ)|1/r², |E⁻¹ ∂E/∂x_δ|1/r, and |E⁻¹ ∂²E/(∂x_δ ∂x_γ)|1/r², where E_tan is the tangent modulus and E is Young’s modulus.