Gradient elasticity with surface energy: mode-III crack problem

In this paper an anisotropic strain-gradient dependent theory of elasticity is exploited, which contains both volumetric and surface energy gradient dependent terms. The theory is applied to the solution of the mode-III crack problem and is extending previous results by Aifantis and co-workers. The two boundary value problems corresponding to the “unclamped” and “clamped” crack tips, respectively, are solved analytically. It turns out that the first problem is physically questionable for some values of the surface energy parameter, whereas the second boundary value problem is leading to a cusping crack, which is consistent with Barenblatt's theory without the incorporation of artificial assumptions.     

Irreversibility and damage of SAFC-40R steel specimen in uniaxial tension

Macroscopic material damage is detected and assessed for the SAFC-40R steel specimen in uniaxial tension even when the stress responded linearly with strain. As the loading increased monotonically at a rate of 0.2 cm/min, the specimen first absorbed heat from the surrounding and then released heat when the strain is almost five times beyond the so-called “elastic limit”. In other words, the specimen undergoes cooling and heating with reference to the ambient temperature. This phenomenon is predicted theoretically for the first time by application of the energy density theory and the results agreed well with experimental data. Obtained is the H-function that possesses a distinct threshold at time between 21 and 22 seconds after loading. This transition is defined as the onset of disorder at which point the energy dissipation density D increases suddenly by one order of magnitude. The corresponding uniaxial stress and strain are 194.4 MPa and 0.9764·10⁻³ cm/cm, respectively. These values are lower than those normally referred to at the yield point.     

Geomechatronics – Interaction between ground and construction machinery and its application to construction robotics

“Geomechatronics” is a technical field in which “Geotechniques” is fused with “Mechatronics” that is the technical field to promote the automatic control of machines by using the electronics. In the field of “Geomechatronics”, a construction machine, which treats geotechnical materials such as soil and rock, automatically evaluates the properties and conditions of the ground and determines the optimum controlling method of itself for the ground with the base of the machine–ground interaction. Some researches for practical use in the field of geomechatronics are introduced, and then the progressing view of this research and technical filed is explained in this paper.     

L''exactitude du rendement des machines agricoles

Precision of the working height, or depth, of earth-moving equipment and other farm machinery can have a marked effect on the quality of the operation. The main factors affecting a machine's performance in maintaining the correct working height are its geometric parameters and the soil surface upon which it operates. The soil surface profile is defined by the “mean slope”, by the “root-mean-square” (RMS) of the deviations from the mean slope, and by the wavelengths calculated by Fourier analysis. The shortest wavelength with an amplitude above 1 cm was used in this study to find the effect of the soil surface waves on the accuracy of farm machinery operations. The machine's geometric parameters affecting its accuracy are its wheelbase and the lengthwise location of the working point (above or under ground) in relation to the wheels. A computer simulation analysis showed that the minimum deviations from the required height are found at the wheels of the machine. The deviations increased at the middle of the machine's span and rose steeply outside the span. The shorter the wheel-span of the machine relative to the wavelengths, the smaller the deviations will be. Smoothed surfaces have longer wavelengths and smaller deviations and as a result, higher accuracy of the working height, or depth, of the machines working on them.     

Molecular dynamics study of size, temperature and rate dependent thermomechanical properties of copper nanofilms

Molecular dynamics simulations are performed to study the thermomechanical properties of copper nanofilms at different temperatures and extremely-high loading rates. The results show a drastic temperature softening effect on the film strength and modulus. The increase of strain rate could result in a much higher strength while the modulus is relatively less affected. It is shown, based on the stress results, that the observed “smaller is softer” and “smaller is stronger” behaviors of nanofilms might be due to the surface plasticity and the volumetric dislocations, respectively. It is also found that the thinner a nanofilm, the smaller the thermal expansion coefficient. The present work reveals that the quasistatic thermomechanical properties of bulk copper at room temperature might be inadequate for the continuum-based study of thermomechanical response of copper nanofilms due to ultrafast laser heating.     

双向纤维腹板增强复合材料夹层板受弯性能实验研究

采用真空导入成型工艺,制备出双向纤维腹板增强复合材料夹层板。在四边简支条件下,对无腹板增强、不同腹板高度和间距的试件进行了集中加载实验,研究其受弯承载力、破坏形态等。结果表明:用双向纤维腹板来增强复合材料夹层板,可显著提高试件承载能力,还能有效减缓面层的剥离破坏;与此同时,腹板高度和间距的减小能减缓面板剥离现象;随着腹板高度增加和腹板间距减小,抗弯能力得到显著提高;增加腹板高度还能有效减小夹层板跨中挠度,但腹板间距对跨中挠度的影响并不明显。     

Effective simulation of flexible lateral boundaries in two- and three-dimensional DEM simulations

Discrete element method （DEM） models to simulate laboratory element tests play an important role in advancing our understanding of the mechanics of granular material response, including bonded or cemented, particulate materials. Comparisons of the macro-scale response observed in a real physical test and a ＂virtual＂ DEM-simulated test can calibrate or validate DEM models. The detailed, particle scale information provided in the DEM simulation can then be used to develop our understanding of the material behaviour. It is important to accurately model the physical test boundary conditions in these DEM simulations. This paper specifically considers triaxial tests as these tests are commonly used in soil mechanics. In a triaxial test, the test specimen of granular material is enclosed within a flexible latex membrane that allows the material to deform freely during testing, while maintaining a specified stress condition. Triaxial tests can only be realistically simulated in 3D DEM codes, however analogue, 2D, biaxial DEM simulations are also often considered as it is easier to visualize particle interactions in two dimensions. This paper describes algorithms to simulate the lateral boundary conditions imposed by the latex membrane used in physical triaxial tests in both 2D and 3D DEM simulations. The importance of carefully considering the lateral boundary conditions in DEM simulations is illustrated by considering a 2D biaxial test on a specimen of frictional unbonded disks and a 3D triaxial test on a bonded （cemented） specimen of spheres. The comparisons indicate that the lateral boundary conditions have a more significant influence on the local, particle-scale response in comparison with the overall macro-scale observations.     

Scattering of pulsed Rayleigh surface waves by a cylindrical cavity

A pulsed Rayleigh surface wave of prescribed shape is incident on a cylindrical cavity which is parallel to both the plane free surface and the plane wave front. Multiple reflections at the cylindrical and plane free surface are considered and the resulting displacements and stress components are calculated in the surrounding of the cavity by approximately summing infinite double sums. Use is made of the stationary loading case simulated by a periodic train of wave pulses and its time Fourier series representation and of expansions of all incident and reflected waves in terms of cylindrical wave functions. For reflection, the free surface of the half-space is approximated by a fictitious convex (or concave) cylindrical surface of “large” radius. The wave pattern due to a single pulse loading is constructed from the stationary solution by enforcing homogeneous initial conditions in the half-space ahead of the single loading pulse and by prescribing a wide spacing in the periodically set-forth train of pulses. The numerical results for stresses and dynamic stress magnification factors are especially useful for the interpretation of recent measurements in dynamic photoelasticity.     

Bounds for the effective properties of heterogeneous plates

This paper presents new bounds for heterogeneous plates which are similar to the well-known Hashin–Shtrikman bounds, but take into account plate boundary conditions. The Hashin–Shtrikman variational principle is used with a self-adjoint Green-operator with traction-free boundary conditions proposed by the authors. This variational formulation enables to derive lower and upper bounds for the effective in-plane and out-of-plane elastic properties of the plate. Two applications of the general theory are considered: first, in-plane invariant polarization fields are used to recover the “first-order” bounds proposed by Kolpakov [Kolpakov, A.G., 1999. Variational principles for stiffnesses of a non-homogeneous plate. J. Meth. Phys. Solids 47, 2075–2092] for general heterogeneous plates; next, “second-order bounds” for n-phase plates whose constituents are statistically homogeneous in the in-plane directions are obtained. The results related to a two-phase material made of elastic isotropic materials are shown. The “second-order” bounds for the plate elastic properties are compared with the plate properties of homogeneous plates made of materials having an elasticity tensor computed from “second-order” Hashin–Shtrikman bounds in an infinite domain.     

A “quasi-elastic” affine formulation for the homogenised behaviour of nonlinear viscoelastic polycrystals and composites

The derivation of the overall behaviour of nonlinear viscoelastic (or rate-dependent elastoplastic) heterogeneous materials requires a linearisation of the constitutive equations around uniform per phase stress (or strain) histories. The resulting Linear Comparison Material (LCM) has to be linear thermoviscoelastic to fully retain the viscoelastic nature of phase interactions. Instead of the exact treatment of this LCM (i.e., correspondence principle and inverse Laplace transforms) as proposed by the “classical” affine formulation, an approximate treatment is proposed here. First considering Maxwellian behaviour, comparisons for a single phase as well as for two-phase materials (with “parallel” and disordered morphologies) show that the “direct inversion method” of Laplace transforms, initially proposed by Schapery (1962), has to be adapted to fit correctly exact responses to creep loading while a more general method is proposed for other loading paths. When applied to nonlinear viscoelastic heterogeneous materials, this approximate inversion method gives rise to a new formulation which is consistent with the classical affine one for the steady-state regimes. In the transient regime, it leads to a significantly more efficient numerical resolution, the LCM associated to the step by step procedure being no more thermoviscoelastic but thermoelastic. Various comparisons for nonlinear viscoelastic polycrystals responses to creep as well as relaxation loadings show that this “quasi-elastic” formulation yields results very close to classical affine ones, even for high contrasts.     

Low-density layer formation and “lifting force” effect at micro- and meso-scale levels

The processes at various scale levels in the contact area of interacting objects under high-energy action will be examined from the viewpoint of mesomechanics. Modeling of contact area at atomic- and meso-scale levels was carried out on the base of discrete computational approach (method of particles). Molecular dynamic method was used at the micro-scale level; movable cellular automata method—at the meso-scale level. The gradient of velocity in areas near the surface leads to formation of low density and fragmented areas. This effect is accompanied by the failure of crystal lattice stability and intensive mixing process at the atomic level. The mechanisms of mass transfer in contact area were discussed. The results allow us to explain a host of experimental data of mechanochemistry such as phase formation at friction surface, alloy formation due to contact interaction under “pressure + shear” loading conditions.     

An elastic-plastic analysis of an asymmetric cracked specimen subjected to longitudinal shear

Some recent elastic-plastic analyses of cracked specimens subjected to symmetric mode III loading are extended to include asymmetric loading and geometry. Solutions are given for arbitrary work hardening behaviour in any specimen that is amenable to a linear elastic analysis. It is shown that asymmetry has a major influence on the shape of the plastic zone, but does not affect the J-integral unil the loading is well into the large scale yielding range. In particular the “plastic zone corrected” estimate of J, obtained by elastically solving a problem for a crack longer than the actual one, is shown to remain a valid two-term asymptotic expansion in the presence of asymmetry. The general results are applied to a crack at an angle to a uniform stress field in a power law hardening material. The growth of the plastic zone is displayed graphically for various hardening exponents and crack orientations. No other asymmetric solution is available, but values of J are compared with those obtained from a fully plastic analysis in the symmetric case.     

Energy approach to crack growth in a fiber-reinforced brittle material modeled by an inclusion

Brittle materials randomly reinforced with a low volume fraction of strong, stiff and ductile fibers are considered, with specific reference to fiber-reinforced cements and concrete. Visible cracks in such materials are accompanied by a surrounding damage zone – together these constitute a very complex “crack system”. Enormous effort has been put into trying to understand the micromechanics of such systems. Almost all of these efforts do not deal with the “crack system” propagation behavior as a whole. The propagation process of such a “crack system” includes propagation of the visible crack and the growth of the damage zone. Propagation may take place by lengthening of the visible crack together with the concomitant lengthening of the surrounding damage zone, or simply by broadening of the damage zone while the visible crack length remains unchanged – or simultaneously by growth of both types. A phenomenological completely theoretical model (for an ideal material) is here proposed which can serve to examine the propagation process by means of energy principles, without recourse to the microscopic details of the process. An application of this theoretical approach is presented for the case of a damage zone evolving with a rectangular shape. This shape is chosen because it is expected that it will illustrate the nature of damage evolution and because the computational procedure necessary to follow the growth is the most straightforward.     

Mechanics of granular materials: The discrete and the continuum descriptions juxtaposed

In recent years a discussion could be followed where the pros and cons of the applicability of the Cosserat continuum model to granular materials were debated [Bardet, J.P., Vardoulakis, I., 2001. The asymmetry of stress in granular media. Int. J. Solids Struct. 38, 353–367; Kruyt, N.P., 2003. Static and kinematics of discrete Cosserat-type granular materials. Int. J. Solids Struct. 40, 511–534; Bagi, K., 2003. Discussion on “The asymmetry of stress in granular media”. Int. J. Solids Struct. 40, 1329–1331; Bardet, J.P., Vardoulakis, I. 2003a. Reply to discussion by Dr. Katalin Bagi. Int. J. Solids Struct. 40, 1035; Kuhn, M., 2003. Discussion on “The asymmetry of stress in granular media”. Int. J. Solids Struct. 40, 1805–1807; Bardet, J.P., Vardoulakis, I., 2003b. Reply to Dr. Kuhn’s discussion. Int. J. Solids Struct. 40, 1809; Ehlers, W., Ramm, E., Diebels, S., D’Addetta, G.A., 2003. From particle ensembles to Cosserat continua: homogenization of contact forces towards stresses and couple stresses. Int. J. Solids Struct. 40, 6681–6702; Chang, C.S., Kuhn, M.R., 2005. On virtual work and stress in granular media. Int. J. Solids Struct. 42, 3773–3793]. The authors follow closely this debate and try, with this paper, to provide a platform where the various viewpoints could find their position. We consider an ensemble of rigid, arbitrarily shaped grains as a set with structure. We establish a basic mathematical framework which allows to express the balance laws and the action–reaction laws for the discrete system in a “global” form, through the concepts of “part”, “granular surface”, “separately additive function” and “flux”. The independent variable in the balance laws is then the arbitrary part of the assembly rather than the single grain. A parallel framework is constructed for Cosserat continua, by applying the axiomatics established by [Noll, W., 1959. The foundation of classical mechanics in the light of recent advances in continuum mechanics. In: The axiomatic method, with special reference to Geometry and Physics, North-Holland Publishing Co., Amsterdam pp. 266–281, Gurtin, M.E., Williams, W.O., 1967. An axiomatic foundation of continuum thermodynamics. Arch. Rat. Mech. Anal. 26, 83–117, Gurtin, M.E., Martins, L.C., 1976. Cauchy’s theorem in classical physics. Arch. Rat. Mech. Anal. 60, 305–324]. The comparison between the two realisations suggests the microscopic interpretation for some features of Cosserat Mechanics, among which the asymmetry of the Cauchy-stress tensor and the couple-stress.     

Simply approximate estimation of local failure position of a free-free slender shell

Dynamic plastic failure characteristics of a space free-free slender shell subjected to intense dynamic loading of suddenly applied pressure unsymmetrical triangle distributed along its span was studied. Both rigid perfectly plastic (r-p-p) analytical method and finite element method based elastic perfectly plastic (e-p-p) material idealization and shell element model were adopted to predict the local failure position in the structure. It was shown that both r-p-p and e-p-p model could estimate a plastic “kink” taking place in the slender shell, which reflects the strain localization of deformation. The comparison for the position of “kink” predicted by using r-p-p and e-p-p methods is found to be reasonable good.     

Stretch and twist of rotating prismatic beam

A turbine blade is modelled as a uniform isotropic prismatic beam of general cross-section and “setting angle” rotating about one end, and is analysed according to the linear theory of elasticity. A semi-inverse solution is presented which reduces the three-dimensional problem to one of two dimensions, and explicit stress and strain components given for the mathematically amenable elliptic cross-section. As expected, the planar stresses σ_x,σ_y, and τ_xy arising from the two-dimensional problem are found to be small. For the general section, the theory predicts small curvature of the blade centre line, and a twisting moment which tends to reduce the “angle of set”.     

Fatigue strength of welded joints based on local, semi-local and nominal approaches

Analysis based on the so-called “local approach” is made to estimate the fatigue strength of welded joints. Numerical analyses or strain gauges are employed for finding the stress and/or strain state in the vicinity of the weld toe. The notch stress intensity factor (NSIF) approach applied to fillet welded joints, as far as the opening angle between the weld and the main plate surface is constant (e.g. 135°, typical for many fillet welds), is able to rationalise the fatigue strength data both for different joint geometries and absolute dimensions. The NSIF approach has been previously developed as an extension of the Linear Elastic Fracture Mechanics (LEFM) to open V-notches and is based on the exponential local stress field around the V-notch tip. Several different “local approaches”, although simpler and more practical than the NSIF, are based on the stress (or strain) values determined beyond the exponential local one. To distinguish such approaches from the NSIF based one, we define the former as semi-local or nominal approaches while the latter is a local approach. The paper underlines that the local approaches, differently from the other ones, are able to unify in a single scatter band the fatigue strength data obtained from welded joints having different geometry and absolute dimensions.     

Micro-plasticity of surface steps under adhesive contact: Part II—Multiple-dislocation mediated contact hardening

The study of micro-plastic behavior of rough surface contacts is the critical link towards a fundamental understanding of contact, friction, adhesion, and surface failures at small length scales. In the companion paper (Yu, H.H., Shrotriya, P., Gao, Y.F., Kim, K.-S., 2007. Micro-plasticity of surface steps under adhesive contact. Part I. Surface yielding controlled by single-dislocation nucleation. J. Mech. Phys. Solids 55, 489–516), we have studied the onset of surface yielding due to single-dislocation nucleation from a stepped surface under adhesive contact. Here we analyze the contact hardening behavior due to multiple dislocations in a two-dimensional dislocation model. Continuum micro-mechanical analyses are used to derive the configurational force on the dislocation, while a modified Rice–Thomson criterion is used to model dislocation nucleation. Dislocations nucleated from the surface step are stabilized and pile up as a result of the balance between the resolved driving force and the non-zero lattice resistance in the solid. The dislocation pileup will exert a strong back stress to prevent further dislocation nucleation and thus lead to the contact hardening behavior, the degree of which depends on the slip-plane orientation. Particularly, we find that dislocation interactions between two slip planes can make the contact loading order-of-magnitude easy to nucleate multiple dislocations, which is thus named “latent softening”. A mechanistic explanation shows that the latent softening is closely related to the stress-concentration mode mixity at the surface step. Dislocation nucleation will modify the geometric characteristics of the surface step, so that the contact-induced stress state near the step, as described by the mode mixity, changes, which influences the subsequent dislocation nucleation. Our calculations show that the dislocation pileup on one slip plane can even cause the spontaneous dislocation nucleation on the other slip plane without further increase of the contact load. Furthermore, it is found that rough surface contacts at small length scale can lead to the dislocation segregation and the formation of a surface tensile sub-layer. The discrete-dislocation model presented here and in the companion paper provides novel insights in bridging the atomistic simulations and continuum plastic flow analysis of surface asperity contact.     

Fatigue crack growth in high and low strength steel under torsional loading

During loading of a crack in mode III the crack surfaces in contact slide against each other giving rise to friction, abrasion and mutual support, thereby reducing the effective stress at the crack tip (“sliding mode crack closure”). This phenomenon was investigated in a high strength steel (AISI 4340) and in a low strength steel (AISI C1018) in circumferentially notched specimens under pure cyclic torsion and combined loading (cyclic torsion plus static axial load). The influence of sliding mode crack closure on fatigue crack propagation is shown and “true” crack growth values (without the sliding mode crack closure influence) are determined on the basis of an extrapolation procedure. Explanations are given for causes of the various fracture modes observed, such as “factory roof” fracture, macroscopically flat mode III fracture and “lamella” fracture. Finally the scientific and technical importance of sliding mode crack closure is demonstrated.