Representative elementary volume analysis of polydisperse granular packings using discrete element method

The representative elementary volume （REV） for three-dimensional polydisperse granular packings was determined using discrete element method simulations. Granular mixtures of various sizes and particle size distributions were poured into a cuboid chamber and subjected to uniaxial compression, Findings showed that the minimum REV for porosity was larger compared with the REV for parameters such as coordination number, effective elastic modulus, and pressure ratio. The minimum REV for porosity and other parameters was found to equal 15,10, and 5 times the average grain diameter, respectively. A study of the influence of sample size on energy dissipation in random packing of spheres has also confirmed that the REV size is about 15 times the average grain diameter. The heterogeneity of systems was found to have no effect on the REV for the parameters of interest for the narrow range of coefficient of uniformity analyzed in this paper. As the REV approach is commonly applied in both experimental and numerical studies, determining minimum REV size for polydisperse granular packings remains a crucial issue.

Determining a representative element volume for DEM simulations of samples with non-circular particles

Numerical studies on the number of particles or system size required to attain a representative element volume(REV)for discrete element method(DEM)simulations of granular materials have almost always considered samples with spherical or circular particles.This study considers how many particles are needed to attain a REV for 2D samples of 2-disc cluster particles where the particle aspect ratio(AR)was systematically varied.Dense and loose assemblies of particles were simulated.The minimum REV was assessed both by considering the repeatability of static packing characteristics and the shearing behaviour in biaxial compression tests,and by investigating the effect of sample size on the measured characteris-tics and observed shearing behaviour.The repeatability of the data considered generally improved with increasing sample size.The packing characteristics of the dense samples were more repeatable suggest-ing that the minimum REV reduces with increasing packing density.The minimum REV was observed to be sensitive to the characteristic measured.Although the overall responses of the samples during shear deformation were similar irrespective of the sample sizes,the smaller the sample size,the higher the fluctuations observed in the responses.Analysis of the coefficient of variation of the fluctuations around the critical state stress ratio can provide insight as to whether a REV is attained.The particle AR influences the effect of sample size on shearing characteristics and thus the minimum number of particles required to attain a REV;this can be explained by the influence of AR on the number of contacts within the samples.

Effects of fracture distribution and length scale on the equivalent continuum elastic compliance of fractured rock masses

Fracture systems have strong influence on the overall mechanical behavior of fractured rock masses dueto their relatively lower stiffness and shear strength than those of the rock matrix. Understanding theeffects of fracture geometrical distribution, such as length, spacing, persistence and orientation, isimportant for quantifying the mechanical behavior of fractured rock masses. The relation betweenfracture geometry and the mechanical characteristics of the fractured rock mass is complicated due tothe fact that the fracture geometry and mechanical behaviors of fractured rock mass are stronglydependent on the length scale. In this paper, a comprehensive study was conducted to determine theeffects of fracture distribution on the equivalent continuum elastic compliance of fractured rock massesover a wide range of fracture lengths. To account for the stochastic nature of fracture distributions, threedifferent simulation techniques involving Oda＇s elastic compliance tensor, Monte Carlo simulation （MCS）,and suitable probability density functions （PDFs） were employed to represent the elastic compliance offractured rock masses. To yield geologically realistic results, parameters for defining fracture distributionswere obtained from different geological fields. The influence of the key fracture parameters andtheir relations to the overall elastic behavior of the fractured rock mass were studied and discussed. Adetailed study was also carried out to investigate the validity of the use of a representative elementvolume （REV） in the equivalent continuum representation of fractured rock masses. A criterion was alsoproposed to determine the appropriate REV given the fracture distribution of the rock mass.

The Revival of Traditional Cultural Brand Genes with Media Convergence

Under the background of media convergence,a new path has been established for the development of traditional cultural brands.By exploring the balance point of representative elements of traditional cultural brand genes and media convergence in the form of digital media,this paper analyzes the founding value and advantages of traditional cultural brand genes,transforms the representative elements of traditional cultural brand genes into visual traces,and extends its information characteristics with digital technology.The research gives new vitality to the genes of traditional cultural brands,enhances the social space and cultural value of traditional cultural brands,and drives the development of digital cultural and creative industries.

Multiscale Mechanics and Optimization of Gastropod Shells

A vast majority of mollusks grow a hard shell for protection. The structure of these shells comprises several levels of hierarchy that increase their strength and their resistance to natural threats. This article focuses on nacreous shells, which are composed of two distinct layers. The outer layer is made of calcite, which is a hard but brittle material, and the inner layer is made of nacre, a tough and ductile material. The inner and outer layers are therefore made of materials with distinct structures and properties. In this article, we demonstrate that this system is optimum to defeat attacks from predators. A two-scale mod- eling and optimization approach was used. At the macroscale, a two-layer finite element model of a seashell was developed to capture shell geometry. At the microscale, a representative volume element of the microstructure of nacre was used to model the elastic modulus of nacre as well as a multiaxial failure criterion, both expressed as function of microstructural parameters. Experiments were also performed on actual shells of red abalone to validate the results obtained from simulations and gain insight into the way the shell fails under sharp perforation. Both optimization and experimental results revealed that the shell displays optimum performance when two modes of failure coincide within the structure. Finally, guidelines for designing two-layer shells were proposed to improve the performance of engineered protective systems undergoing similar structural and loading conditions.

Behaviour of granular matter under gravity-induced stress gradient:A two-dimensional numerical investigation

Gravity is the most important load source in mining and geotechnical engineering,which causes both the stress level and stress gradient inside geomaterials.Different from the stress level,the influence of gravity-induced stress gradient on the behaviour of the material is still unknown.An in-deep study on it will help to promote the understanding of material behaviour,especially for those cases related to unconventional gravity such as terrestrial ng physical modelling and extraterrestrial resource exploitation(g is the terrestrial gravitational acceleration).In this study,a high-order homogenization for granular materials is proposed at first,in which the stress gradient is drawn into the constitutive representation by adopting a representative volume element(RVE).The consolidation and shear strength behaviour of RVE are then investigated by performing numerical biaxial tests.The results show that all the compressibility,shear strength,shear stiffness,volumetric deformation,and critical state behaviour show a stress gradient dependence.A coupling between stress gradient,stress level,and material properties is also observed.These observations suggest that,besides the stress level,extra attention needs to be paid to material responses related to stress gradient during engineering practices.

THEORETICAL PREDICTION OF TOOL-CHIP CONTACT LENGTH IN ORTHOGONAL METAL MACHINING BY COMPUTER SIMULATION

A method for determination of tool-chip contact length is theoreticallypresented in orthogonal metal machining. By using computer simulation and based on the analyses ofthe elastro-plastic deformation with lagrangian finite element method in the deformation zone, theaccumulated representative length of the low layer, the tool-chip contact length of the chipcontacting the tool rake are calculated, experimental studies are also carried out with 0.2 percentcarbon steel. It is shown that the tool-chip contact lengths obtained from computer simulation havea good agreement with those of measured values.

Multiscale Homogenization Analysis of Alkali–Silica Reaction (ASR) Effect in Concrete

The alkali silica reaction (ASR) is one of the major long-term deterioration mechanisms occurring in con- crete structures subjected to high humidity levels, such as bridges and dams. ASR is a chemical reaction between the silica existing inside the aggregate pieces and the alkali ions from the cement paste. This chemical reaction produces ASR gel, which imbibes additional water, leading to gel swelling. Damage and cracking are subsequently generated in concrete, resulting in degradation of its mechanical proper- ties. In this study, ASR damage in concrete is considered within the lattice discrete particle model (LDPM), a mesoscale mechanical model that simulates concrete at the scale of the coarse aggregate pieces. The authors have already modeled successfully ASR within the LDPM framework and they have calibrated and validated the resulting model, entitled ASR-LDPM, against several experimental data sets. In the pre- sent work, a recently developed multiscale homogenization framework is employed to simulate the macroscale effects of ASR, while ASR-LDPM is utilized as the mesoscale model. First, the homogenized behavior of the representative volume element (RVE) of concrete simulated by ASR-LDPM is studied under both tension and compression, and the degradation of effective mechanical properties due to ASR over time is investigated. Next, the developed homogenization framework is utilized to reproduce experimental data reported on the free volumetric expansion of concrete prisms. Finally, the strength degradation of prisms in compression and four-point bending beams is evaluated by both the mesoscale model and the proposed multiscale approach in order to analyze the accuracy and computational ef - ciency of the latter. In all the numerical analyses, different RVE sizes with different inner particle realiza- tions are considered in order to explore their effects on the homogenized response.

Three Phase Composite Cylinder Assemblage Model for Analyzing the Elastic Behavior of MWCNT-Reinforced Polymers

Evolution of computational modeling and simulation has given more emphasis on the research activities related to carbon nanotube(CNT)reinforced polymer composites recently.This paper presents the composite cylinder assemblage(CCA)approach based on continuum mechanics for investigating the elastic properties of a polymer resin reinforced by multi-walled carbon nanotubes(MWCNTs).A three-phase cylindrical representative volume element(RVE)model is employed based on CCA technique to elucidate the effects of inter layers,chirality,interspacing,volume fraction of MWCNT,interphase properties and temperature conditions on the elastic modulus of the composite.The interface region between CNT and polymer matrix is modeled as the third phase with varying material properties.The constitutive relations for each material system have been derived based on solid mechanics and proper interfacial traction continuity conditions are imposed.The predicted results from the CCA approach are in well agreement with RVE-based finite element model.The outcomes reveal that temperature softening effect becomes more pronounced at higher volume fractions of CNTs.

Elastoplastic analysis of AA7075-O aluminum sheet by hybrid micro-scale representative volume element modeling with really-distributed particles and in-situ SEM experimental testing

This paper addresses the challenge of reconstructing randomly distributed second-phase particlestrengthened microstructure of AA7075-O aluminum sheet material for computational analysis.The particle characteristics in 3D space were obtained from focused ion beam and scanning electron microscopy(FIB-SEM)and SEM-based Electron Backscatter Diffraction/Energy Dispersive X-ray Spectrometry(EBSD/EDS)techniques.A theoretical framework for analysis of elastic-plastic deformation of such3D microstructures is developed.Slip-induced shear band formation,void initiation,growth and linkage at large plastic strains during uniaxial tensile loading were investigated based on reconstructed 3D representative volume element(RVE)models with real-distribution of particles and the results compared with experimental observations.In-situ SEM interrupted tension tests along transverse direction(TD)and rolling direction(RD),employing microscopic-digital image correlation(μ-DIC)technique,were carried out to investigate slip bands,micro-voids formation and obtain microstructural strain maps.The resulting local strain maps were analyzed in relation to the experimentally observed plastic flow localization,failure modes and local stress maps from simulations of RVE models.The influences of particle size,shape,orientation,volume fraction as well as matrix-particle interface properties on local plastic deformation,global stress-strain/strain-hardening curves and interfacial failure mechanisms were studied based on 3D RVE models.When possible,the model results were compared with in-situ tensile test data.In general,good agreement was observed,indicating that the real 3D microstructure-based RVE models can accurately predict the plastic deformation and interfacial failure evolution in AA7075-O aluminum sheet.

Stress Analysis of Wire Strands by Mesoscale Mechanics

Steel wire ropes have wide application in a variety of engineering fields such as ocean engineering and civil engineering.The stress calculation for steel wire ropes is of crucial importance when conducting strength and fatigue analyses.In this study,we performed a finite element analysis of single-strand steel wire ropes.For the geometric modeling,we used an analytic geometry of space method.We established helical line equations and used the coordinates of the contact points.The finite-element model was simplified using the periodic law.Periodic boundary conditions were used to simulate a wire strand of infinite length under tensile strain,for which we calculated the cross-sectional stresses and inner forces.The results showed that bending and torsion moments emerged when the wire strand was under tensile load.In some cases,the bending stress reached 18%of the tensile stress,and the torsion stress reached 29%of the tensile stress,which means that the total stress was higher than the nominal stress.Whereas in ear-lier studies,a conservative prediction of nominal stress was not possible,the results of our strength and fatigue analyses were more conservative.

Computational model for short-fiber composites with eigenstrain formulation of boundary integral equations

A computational model is proposed for short-fiber reinforced materials with the eigenstrain formulation of the boundary integral equations （BIE） and solved with the newly developed boundary point method （BPM）. The model is closely derived from the concept of the equivalent inclusion Of Eshelby tensors. Eigenstrains are iteratively determined for each short-fiber embedded in the matrix with various properties via the Eshelby tensors, which can be readily obtained beforehand either through analytical or numerical means. As unknown variables appear only on the boundary of the solution domain, the solution scale of the inhomogeneity problem with the model is greatly reduced. This feature is considered significant because such a traditionally time-consuming problem with inhomogeneity can be solved most cost-effectively compared with existing numerical models of the FEM or the BEM. The numerical examples are presented to compute the overall elastic properties for various short-fiber reinforced composites over a representative volume element （RVE）, showing the validity and the effectiveness of the proposed computational modal and the solution procedure.

MICRO/MACROMECHANICAL PLASTIC LIMIT ANALYSES OF COMPOSITE MATERIALS AND STRUCTURES

The load-bearing capacities f ductile composite materials andstructures are studied by means of a combined micro/macromechanicsapproach. Firstly, on the microscopic scale, the aim is to get themacroscopic strength domains by means of the homogenization theory ofmicromechanics. A representative volume element (RVE) is selected toreflect the microstructures of the composite materials. Byintroducing the homogenization theory into the kinematic limittheorem of plastic limit analysis, an optimization format to directlycalculate the limit loads of the RVE is obtained. And the macroscopicyield criterion can be deter- mined according to the relation betweenmacroscopic and microscopic fields.

Multi-scale finite element simulation on large deformation behavior of wood under axial and transverse compression conditions

Multi-scale finite element method is adopted to simulate wood compression behavior under axial and transverse loading. Representative volume elements (RVE) of wood microfibril and cell are proposed to analyze orthotropic mechanical behavior. Lignin, hemicellulose and crystalline-amorphous cellulose core of spruce are concerned in spruce nanoscale model. The equivalent elastic modulus and yield strength of the microfibril are gained by the RVE simulation. The anisotropism of the crystalline-amorphous cellulose core brings the microfibril buckling deformation during compression loading. The failure mechanism of the cell-wall under axial compression is related to the distribution of amorphous cellulose and crystalline cellulose. According to the spruce cell observation by scanning electron microscope, numerical model of spruce cell is established using simplified circular hole and regular hexagon arrangement respectively. Axial and transverse compression loadings are taken into account in the numerical simulations. It indicates that the compression stress-strain curves of the numerical simulation are consistent with the experimental results. The wood microstructure arrangement has an important effect on the stress plateau during compression process. Cell-wall buckling in axial compression induces the stress value drops rapidly. The wide stress plateau duration means wood is with large energy dissipation under a low stress level. The numerical results show that loading velocity affects greatly wood microstructure failure modes in axial loading. For low velocity axial compression, shear sliding is the main failure mode. For high velocity axial compression, wood occur fold and collapse. In transverse compression, wood deformation is gradual and uniform, which brings stable stress plateau.

STRENGTH PROPERTIES OF JOINTED ROCK MASSES BASED ON THE HOMOGENIZATION METHOD

This paper aims to determine the strength properties of jointed rock masses by means of the homogenization method. To reflect the microstructure of jointed rock masses, a representative element volume （REV） is selected. Assuming being rigid and perfectly plastic and utilizing the Mohr-Coulomb yield criterion for rock and joints, an upper bound limit analysis method is proposed based on the homogenization theory. By using the finite element method and Goodman joint element, a nonlinear mathematical programming with equality constraints is formulated, which can be solved by a direct iterative algorithm. Numerical results show the strength of jointed rock masses behaves anisotropy with different joint directions. This method presents an effective tool for the strength analysis of jointed rock masses.

Multi-scale thermodynamic analysis method for 2D SiC/SiC composite turbine guide vanes

Ceramic Matrix Composite （CMC） turbine guide vanes possess multi-scale stress and strain with inhomogeneity at the microscopic scale. Given that the macroscopic distribution cannot reflect the microscopic stress fluctuation, the macroscopic method fails to meet the requirements of stress and strain analysis of CMC turbine guide vanes. Furthermore, the complete thermodynamic properties of 2D woven SiC/SiC-CMC cannot be obtained through experimentation, Accordingly, a method to calculate the thermodynamic properties of CMC and analyze multi-scale stress and strain of the turbine guide vanes should be established. In this study, the multi-scale thermodynamic analysis is investigated. The thermodynamic properties of Chemical Vapor Infiltration （CVI） pro- cessed SiC/SiC-CMC are predicted by a Representative Volume Element （RVE） model with porosity, leading to the result that the relative error between the calculated in-plane tensile modulus and the experimental value is 4.2%. The macroscopic response of a guide vane under given conditions is predicted. The relative error between the predicted strain on the trailing edge and the experimental value is 9.7%. The calculation of the stress distribution of micro-scale RVE shows that the maximum value of microscopic stress, which is located in the interlayer matrix, is more than 1.5 times that of macroscopic stress in the same direction and the microscopic stress distribution of the interlayer matrix is related to the pore distribution of the composite.

Simulations of deformation and damage processes of SiCp/Al composites during tension

The deformation, damage and failure behaviors of 17 vol.% SiCp/2009AI composite were studied by micro- scopic finite element （FE） models based on a representative volume element （RVE） and a unit cell. The RVE having a 3D realistic microstructure was constructed via computational modeling technique, in which an interface phase with an average thickness of 50 nm was generated for assessing the effects of interracial properties. Modeling results showed that the RVE based FE model was more accurate than the unit cell based one. Based on the RVE, the predicted stress-strain curve and the fracture morphology agreed well with the experimental results. Furthermore, lower interface strength resulted in lower flow stress and ductile damage of interface phase, thereby leading to decreased elongation. It was revealed that the stress concentration factor of SiC was -2.0： the average stress in SiC particles reached -1200 MPa, while that of the composite reached -600 MPa.

A fast numerical method of introducing the strengthening effect of residual stress and strain to tensile behavior of metal matrix composites

Thermal residual stress and strain(TRSS)in particle reinforced metal matrix composites(PRMMCs)are believed to cause strengthening effects,according to previous studies.Here,the representative volume element(RVE)based computational homogenization technique was used to study the tensile deformation of PRMMCs with different particle aspect ratios(AR).The influence of TRSS was assessed quantitatively via comparing simulations with or without the cooling process.It was found that the strengthening effect of TRSS was affected by the particle AR.With the average strengthening effect of TRSS,a fast method of introducing the strengthening effect of TRSS to the tensile behavior of PRMMCs was developed.The new method has reduced the computational cost by a factor 2.The effect of TRSS on continuous fiber-reinforced metal matrix composite was found to have a softening-effect during the entire tensile deformation process because of the pre-yield effect caused by the cooling process.

Influence of Lamellar Direction in Pearlitic Steel Wire on Mechanical Properties and Microstructure Evolution

During cold drawing of pearlitic steel wire, the lamellar structure becomes gradually aligned with the draw ing axis, which contributes to the ultra high strength. A direct simulation about the mechanical behaviors and microstructural evolution of pearlitic lamellae was presented. A representative volume element （RVE） containing one pearlitic colony was established based on the real transmission electron microscope （TEM） observation. The deformation of pearlitic colony during tension, shear and wire drawing were successfully simulated. The numerical results show that this metallographic texture leads to a strong anisotropy. The colony has higher yielding stress when the la mellar direction is parallel and perpendicular to the tensile direction. The lamellar evolution is strongly dependent on the initial direction and deformation mode. The formation of typical period shear bands is analyzed. In the wire draw ing, the pearlitic colony at the sub surface experiences a complex strain path： rotation, stretching along the die sur face, and rotation back.

Two-scale modeling of granular materials:A FEM-FEM approach

In the present paper,a homogenization-based two-scale FEM-FEM model is developed to simulate compactions of visco-plastic granular assemblies.The granular structure consisting of two-dimensional grains is modeled by the microscopic finite element method at the small-scale level,and the homogenized viscous assembly is analyzed by the macroscopic finite element method at large-scale level.The link between scales is made using a computational homogenization method.The two-scale FEM-FEM model is developed in which each particle is treated individually with the appropriate constitutive relations obtained from a representative volume element,kinematic conditions,contact constraints,and elimination of overlap satisfied for every particle.The method could be used in a variety of problems that can be represented using granular media.