Numerical analysis and elastic–plastic deformation behavior of anisotropically damaged solids

The present paper is concerned with the numerical modelling of the large elastic–plastic deformation behavior and localization prediction of ductile metals which are sensitive to hydrostatic stress and anisotropically damaged. The model is based on a generalized macroscopic theory within the framework of nonlinear continuum damage mechanics. The formulation relies on a multiplicative decomposition of the metric transformation tensor into elastic and damaged-plastic parts. Furthermore, undamaged configurations are introduced which are related to the damaged configurations via associated metric transformations which allow for the interpretation as damage tensors. Strain rates are shown to be additively decomposed into elastic, plastic and damage strain rate tensors. Moreover, based on the standard dissipative material approach the constitutive framework is completed by different stress tensors, a yield criterion and a separate damage condition as well as corresponding potential functions. The evolution laws for plastic and damage strain rates are discussed in some detail. Estimates of the stress and strain histories are obtained via an explicit integration procedure which employs an inelastic (damage-plastic) predictor followed by an elastic corrector step. Numerical simulations of the elastic–plastic deformation behavior of damaged solids demonstrate the efficiency of the formulation. A variety of large strain elastic–plastic-damage problems including severe localization is presented, and the influence of different model parameters on the deformation and localization prediction of ductile metals is discussed.     

Experiments and theory in strain gradient elasticity

Conventional strain-based mechanics theory does not account for contributions from strain gradients. Failure to include strain gradient contributions can lead to underestimates of stresses and size-dependent behaviors in small-scale structures. In this paper, a new set of higher-order metrics is developed to characterize strain gradient behaviors. This set enables the application of the higher-order equilibrium conditions to strain gradient elasticity theory and reduces the number of independent elastic length scale parameters from five to three. On the basis of this new strain gradient theory, a strain gradient elastic bending theory for plane-strain beams is developed. Solutions for cantilever bending with a moment and line force applied at the free end are constructed based on the new higher-order bending theory. In classical bending theory, the normalized bending rigidity is independent of the length and thickness of the beam. In the solutions developed from the higher-order bending theory, the normalized higher-order bending rigidity has a new dependence on the thickness of the beam and on a higher-order bending parameter, b_h. To determine the significance of the size dependence, we fabricated micron-sized beams and conducted bending tests using a nanoindenter. We found that the normalized beam rigidity exhibited an inverse squared dependence on the beam's thickness as predicted by the strain gradient elastic bending theory, and that the higher-order bending parameter, b_h, is on the micron-scale. Potential errors from the experiments, model and fabrication were estimated and determined to be small relative to the observed increase in beam's bending rigidity. The present results indicate that the elastic strain gradient effect is significant in elastic deformation of small-scale structures.     

Size effects on strength, toughness and fatigue crack growth of gradient elastic solids

The present study investigates the microstructural size effect on the strength of a bar under axial loading, and on the toughness and crack growth of a beam under three-point bending within the framework of strain gradient elasticity. The gradient responses have been found considerably tougher as compared to the classical theory predictions and the observed deviation increases with increasing values of the non-dimensional parameter g/L (microstructural length over structural length). Based on the analytical solution of the strain energy release rate for the three-point bending case, a new, simple and universal, strain gradient elasticity, brittle fracture criterion and a new, size adjusted fatigue crack growth law have been established. Finally, the analytical predictions of the current modeling compare well with previous experimental data, based on three-point bending tests on single-edge notched concrete beams.     

ANCF索梁单元应变耦合问题与模型解耦

索粱结构在土木工程、航空航天等领域有着广泛的应用.在各类索梁动力学建模方法中,由于绝对节点坐标方法(absolute nodal coordinate formulation,ANCF)能够描述柔性体的大变形和大转动问题,因此非常适合大变形索梁结构的动力学建模.对绝对节点坐标索梁单元的应变进行分析可知,弯曲变形会引起单元内部轴向应变的不均匀分布,即单元轴向应变与弯曲应变相互耦合.这种应变耦合效应使单元产生伪应变能,导致单元刚度增大,造成单元失真.分析不同弯曲角下的单元应变及应变能可知,弯曲变形越大,单元失真越严重.通过构造等效一维杆单元重新描述轴向应变,实现了轴向应变与弯曲应变解耦.在此基础上推导广义弹性力,得到了绝对节点坐标索梁单元的应变解耦模型.对解耦前后的两种梁模型进行静力学和动力学仿真,结果表明;解耦模型消除了单元伪应变,相比原模型表现出更好的收敛性和曲率连续性,在相同单元数目下具有更高的精度.同时由于解耦模型降低了单元刚度,因此相比原模型,速度曲线中不再有高频振动.     

Delamination fracture analysis in the GII–GIII plane using prestressed transparent composite beams

In this paper, the mixed-mode II/III version of the prestressed end-notched flexure fracture specimen is developed, which combines the well-known end-notched flexure and the modified split-cantilever beam specimens using a special rig. The new beam-like specimen is able to provide any combination of the mode-II and mode-III strain energy release rates. The mode-III part of the strain energy release rate is fixed by using the special rig, which loads the specimen in the plane of the delamination. The mode-II part of the strain energy release rate is provided by the external load using a three-point bending fixture. A simple closed-form solution using beam theory is developed for the strain energy release rates of the new configuration. The applicability and the limitations of the novel fracture mechanical test are demonstrated using unidirectional E-glass/polyester composite specimens. If only crack propagation onset is involved then the mixed-mode II/III prestressed end-notched flexure specimen can be used to obtain the fracture criterion of transparent composite materials in a relatively simple way.     

Large deflections of a beam subject to three-point bending

In this paper, a solution for the equilibrium configuration of an elastic beam subject to three-point bending is given in terms of Jacobi elliptical functions. General equations are derived, and the domain of the solution is established. Several examples that illustrate a use of the solution are discussed. The obtained numerical results are compared with the results of other authors. An approximation formula by which the beam load is given as a polynomial function of beam deflection is also derived. The range of applicability of the approximation is illustrated by numerical examples.     

考虑尺度效应的微梁静态弯曲数值分析

基于修正偶应力理论及表面弹性理论,本文提出了一种新的双曲线剪切变形梁模型,用于均匀微尺度梁的静态弯曲分析。该理论可以直接利用本构关系获得横向剪切应力,满足梁顶部和底部的无应力边界条件,避免了引入剪切修正因子。根据广义Young-Laplace方程建立了梁的内部与表面层的应力连续性条件,单一的变量场可以描述梁的位移模式。通过在位移场中考虑表面层厚度以及表面层的应力连续条件,可以使新模型能够更准确地预测微尺寸和表面能相关的尺度效应。通过Hamilton原理推导出了梁的控制方程和边界条件。应变能除了考虑经典弹性理论,还要考虑微结构效应和表面能。Navier-type的解析解适用于简支边界条件,而基于拉格朗日插值的微分求积法(DQEM)可以研究在不同边界条件下的力学响应。把该数值解与Navier方法得出的解析解作了对比,得出:微尺度梁在考虑表面能或微尺寸效应、不同载荷和梁高变化下的响应一致;当不考虑微结构相关性和表面能效应时,该模型退化为经典的欧拉梁模型。     

A framework for numerical integration of crystal elasto-plastic constitutive equations compatible with explicit finite element codes

In this paper we summarize the elements of a numerical integration scheme for elasto-plastic response of single crystals. This is intended to be compatible with large-scale explicit finite element codes and therefore can be used for problems involving multiple crystals and also overall behavior of polycrystalline materials. The steps described here are general for anisotropic elastic and plastic response of crystals. The crystallographic axes of the lattice are explicitly stored and updated at each time step. A plastic predictor–elastic corrector scheme is used to calculate the plastic strain rates on all active slip systems based on a rate-dependent physics-based constitutive model without the need of further auxiliary assumptions. Finally we present the results of numerous calculations using a physics-based rate- and temperature-dependent model of copper and the effect of elastic unloading, elastic crystal anisotropy, and deformation-induced lattice rotation are emphasized.     

Development of a hybrid model based on mesh and meshfree methods and its application to fluid–elastic structure interaction for free surface waves

In this paper, a hybrid scheme, Fluid–Fluid–Elastic Structure (FFES) model was developed in the time domain to address the wave breaking impact on the structure. The model is developed based on the partitioned approach with different governing equations that describe various regions of the model domain. The fluid–fluid model denotes that two different fluid models were used to describe fluid in the actual physical domain. The method is a physics-based approximation to reduce the computational time, i.e. in the far-field inviscid fluid (fully nonlinear potential flow theory model), and near to the structure, viscous fluid (Navier Stokes model) is used. The coupled model then interacts with the elastic structure (based on Euler–Bernoulli beam theory). The system of equations is strongly coupled both in space and time. The Fluid–Fluid coupling uses an implicit predictor–corrector scheme, and the fluid–structure coupling works based on an iterative scheme. This approach makes the method more robust and for future extension. Three different possibilities for introducing the coupling was identified and implemented. The model was validated against results from the analytical solution and literature. The method proposed is a reliable, robust, and efficient alternative for simulating fluid–structure interaction problems.     

A thermoviscoelastic model for amorphous shape memory polymers: Incorporating structural and stress relaxation

A thermoviscoelastic constitutive model is developed for amorphous shape memory polymers (SMP) based on the hypothesis that structural and stress relaxation are the primary molecular mechanisms of the shape memory effect and its time-dependence. This work represents a new and fundamentally different approach to modeling amorphous SMPs. A principal feature of the constitutive model is the incorporation of the nonlinear Adam–Gibbs model of structural relaxation and a modified Eyring model of viscous flow into a continuum finite–deformation thermoviscoelastic framework. Comparisons with experiments show that the model can reproduce the strain–temperature response, the temperature and strain-rate dependent stress–strain response, and important features of the temperature dependence of the shape memory response. Because the model includes structural relaxation, the shape memory response also exhibits a dependence on the cooling and heating rates.     

Modeling size-dependent thermo-mechanical behaviors of shape memory polymer Bernoulli-Euler microbeam

The objective of this paper is to model the size-dependent thermo-mechanical behaviors of a shape memory polymer(SMP) microbeam. Size-dependent constitutive equations, which can capture the size effect of the SMP, are proposed based on the modified couple stress theory(MCST). The deformation energy expression of the SMP microbeam is obtained by employing the proposed size-dependent constitutive equation and Bernoulli-Euler beam theory. An SMP microbeam model, which includes the formulations of deflection, strain, curvature, stress and couple stress, is developed by using the principle of minimum potential energy and the separation of variables together. The sizedependent thermo-mechanical and shape memory behaviors of the SMP microbeam and the influence of the Poisson ratio are numerically investigated according to the developed SMP microbeam model. Results show that the size effects of the SMP microbeam are significant when the dimensionless height is small enough. However, they are too slight to be necessarily considered when the dimensionless height is large enough. The bending flexibility and stress level of the SMP microbeam rise with the increasing dimensionless height, while the couple stress level declines with the increasing dimensionless height.The larger the dimensionless height is, the more obvious the viscous property and shape memory effect of the SMP microbeam are. The Poisson ratio has obvious influence on the size-dependent behaviors of the SMP microbeam. The paper provides a theoretical basis and a quantitatively analyzing tool for the design and analysis of SMP micro-structures in the field of biological medicine, microelectronic devices and micro-electro-mechanical system(MEMS) self-assembling.     

Vibrations and Dissipative Heating of a Viscoelastic Beam under a Moving Load

The vibrations and dissipative heating of an infinite viscoelastic beam under a moving load are studied on the basis of Timoshenko beam theory. The influence of transverse-shear strain and rotary inertia on the critical velocities of the moving load, the amplitude of bending vibration, and the temperature of dissipative heating is analyzed__________Translated from Prikladnaya Mekhanika, Vol. 41, No. 1, pp. 69–76, January 2005.     

An experimental study of bending impact waves in beams

An experimental study of the propagation of flexural waves in an elastic beam of circular cross section subjected to an approximate-step-function bending moment is given. The test beam was a low-carbon-steel bar 2 in. diam × 30-ft long and was suspended in a vertical position by a pin located near the upper end. The step moment was applied at the upper pinned end of the bar by an arrangement of two high-pressure, nitrogen-operated cylinders. The strains were measured with strain gages located at eight stations along the bar. The experimental results are compared with results obtained from a solution of the “Timoshenko” beam by W. Flügge and E. E. Zajac for a semi-infinite pinned-end beam subjected to a suddenly applied bending moment at the pinned end. The experimental results are correlated with the predictions of the theoretical solution.     

Crack kinking in sandwich structures under three-point bending

Sandwich beams under three-point bending containing cracks in the core material very close to the upper skin interface are investigated. The cracks considered parallel or with an imperceptible inclination to the longitudinal beam axis and at different distances from the upper skin interface, are analyzed with static non-linear elastic two and three dimensional finite element analyses. From the proposed analyses stress intensity factors are calculated using linear elastic fracture mechanics. It is demonstrated that the crack propagation on the compression side of the core is mainly subjected in shear. The strain energy density criterion is used in order to determine the angle of kinking of the crack into the core.     

Coupled thermoviscoplasticity of glassy polymers in the logarithmic strain space based on the free volume theory

The paper outlines a constitutive model for finite thermo-visco-plastic behavior of amorphous glassy polymers and considers details of its numerical implementation. In contrast to existing kinematical approaches to finite plasticity of glassy polymers, the formulation applies a plastic metric theory based on an additive split of Lagrangian Hencky-type strains into elastic and plastic parts. The analogy between the proposed formulation in the logarithmic strain space and the geometrically linear theory of plasticity, makes this constitutive framework very transparent and attractive with regard to its numerical formulation. The characteristic strain hardening of the model is derived from a polymer network model. We consider the particularly simple eight chain model, but also comment on the recently developed microsphere model. The viscoplastic flow rule in the logarithmic strain space uses structures of the free volume flow theory, which provides a highly predictive modeling capacity at the onset of viscoplastic flow. The integration of this micromechanically motivated approach into a three-dimensional computational model is a key concern of this work. We outline details of the numerical implementation of this model, including elements such as geometric pre- and post-transformations to/from the logarithmic strain space, a thermomechanical operator split algorithm consisting of an isothermal mechanical predictor followed by a heat conduction corrector and finally, the consistent linearization of the local update algorithm for the dissipative variables as well as its relationship to the global tangent operator. The performance of the proposed formulation is demonstrated by means of a spectrum of numerical examples, which we compare with our experimental findings.     

Basic equations for thermoviscoelastic plates with distributed actuators under monoharmonic loading

The basic equations for thin-walled thermoviscoelastic plates with distributed piezoelectric actuators under monoharmonic mechanical and electric loads are derived. The thermomechanical behavior of materials is described using the concept of complex characteristics. Variational methods of solving nonlinear problems of active damping of the bending vibrations of plates are considered. The effect of dissipative heating on the damping of the axisymmetric bending vibrations of a circular plate is examined as an example Translated from Prikladnaya Mekhanika, Vol. 45, No. 2, pp. 107–123, February 2009.     

On a Theory of Thermoviscoelastic Materials with Voids

In this paper we extend the theory of elastic materials with voids to the case when the time derivative of the strain tensor and the time derivative of the gradient of the volume fraction are included in the set of independent constitutive variables. First, the basic equations of the nonlinear theory of thermoviscoelastic materials with voids are established. Then, the linearized version of the theory is derived. We establish a uniqueness result and the continuous dependence of solution upon the initial data and supply terms. A solution of the field equations is also presented.     

饱和多孔弹性Timoshenko梁的大挠度分析

基于微观不可压饱和多孔介质理论和弹性梁的大挠度变形假设,考虑梁剪切变形效应,在梁轴线不可伸长和孔隙流体仅沿轴向扩散的限定下,建立了饱和多孔弹性Timoshenko梁大挠度弯曲变形的非线性数学模型.在此基础上,利用Galerkin截断法,研究了两端可渗透简支饱和多孔Timoshenko梁在突加均布横向载荷作用下的拟静态弯曲,给出了饱和多孔 Timoshenko梁弯曲变形时固相挠度、弯矩和孔隙流体压力等效力偶等随时间的响应.比较了饱和多孔Timoshenko梁非线性大挠度和线性小挠度理论以及饱和多孔 Euler-Bernoulli梁非线性大挠度理论的结果,揭示了他们间的差异,指出当无量纲载荷参数q＞l0时,应采用饱和多孔Timoshenko梁或Euler-Bernoulli梁的大挠度数学模型进行分析,特别的,当梁长细比λ＜30时,应采用饱和多孔Timoshenko梁大挠度数学模型进行分析.     

Basic equations of the theory of thermoviscoelastic plates with distributed sensors and actuators

The basic equations of the theory of thermoviscoelastic thin-walled plates with piezoelectric sensors and actuators under monoharmonic mechanical and electric loading are derived using the Kirchhoff–Love hypotheses. The thermomechanical behavior of passive and piezoactive materials is described using the concept of complex characteristics. Methods of solving nonlinear problems of active damping of thermomechanical vibrations of plates with sensors and actuators are considered. The effect of dissipative heating on the damping of axisymmetric vibrations of a thermoviscoelastic solid circular plate is analyzed as an example     

A note on the problem of pure bending for linear elastic materials with voids

This note concerns the problem of quasi-static pure bending of a beam in the context of the complete theory of linear elastic materials with voids presented in [1]. It is shown here that the solution in the context of the complete theory of [1] is coincident with the pure bending solution of classical elasticity for small time, and that the solution for large time is the bending solution given in [1], a solution which neglected the rate effect in the complete theory of [1]. In between these two limit solutions the rate effect moderates a monotonic transition.