Modeling of textural anisotropy in granular materials with stochastic micro-polar hypoplasticity

The paper deals with a numerical analysis of the effect of textural anisotropy on the behaviour of cohesionless granular materials with consideration of shear localization. For a simulation of the mechanical behaviour of a granular material during a monotonic deformation path, an isotropic micro-polar hypoplastic constitutive model was used. To describe textural effects, spatially correlated random fields of the initial void ratio were subject to rotation against the horizontal axis. The 2D random fields were generated using a conditional rejection method. The results were compared with those obtained with an anisotropic micro-polar constitutive model for a uniform distribution of the initial void ratio. The calculations were carried out with an initially dense granular specimen during plane strain compression under constant lateral pressure.     

Dual variational formulation for Trefftz finite element method of elastic materials

A dual variational principle is presented for Trefftz finite element analysis. The proof of the stationary conditions of the variational functional and the theorem on the existence of extremum are provided in this paper. They are boundary displacement condition, surface traction condition and interelement continuity condition. Based on the assumed intraelement and frame fields, element stiffness matrix equation is obtained which can easily be implemented into computer programs for numerical analysis with Trefftz finite element method. Two numerical examples are considered to illustrate the effectiveness and applicability of the proposed element model.     

Boundary effects on behaviour of granular material during plane strain compression

The paper investigates the boundary effect on the behaviour of granular materials during plane strain compression using finite element method. A micro-polar hypoplastic constitutive model was used. The numerical calculations were carried out with different initial densities and boundary conditions. The behaviour of initially dense, medium dense and loose sand specimen with very smooth or very rough horizontal boundary was investigated. The formation of shear zones gave rise to different global and local stress and strain. Comparisons of the mobilized internal friction, dilatancy and non-coaxiality between global and local quantities were made.     

Experimental and numerical study of the rotation and the erosion of fillers suspended in viscoelastic fluids under simple shear flow

When a porous agglomerate immersed in a fluid is submitted to a shear flow, hydrodynamic stresses acting on its surface may cause a size reduction if they exceed the cohesive stress of the agglomerate. The aggregates forming the agglomerate are slowly removed from the agglomerate surface. Such a behaviour is known when the suspending fluid is Newtonian but unknown if the fluid is viscoelastic. By using rheo-optical tools, model fluids, carbon black agglomerates and particles of various shapes, we found that the particles had a rotational motion around the vorticity axis with a period which is independent on shape (flat particles not considered), but which is exponentially increasing with the elasticity of the medium expressed by the Weissenberg number (We). Spherical particles are always rotating for We up to 2.6 (largest investigated We in this study) but elongated particles stop rotating for We>0.9 while orienting along the flow direction. Erosion is strongly reduced by elasticity. Since finite element numerical simulation shows that elasticity increases the local stress around a particle, the origin of the erosion reduction is interpreted as an increase of cohesiveness of the porous agglomerate due to the infiltration of a viscoelastic fluid.     

Interaction of microcracks with a macrocrack yielded in a narrow strip

Plastic zone growth of collinear cracks has had a longstanding interest in ductile fracture. This work further considers yield zone growth in an isotropic, homogeneous elastic–perfectly plastic infinite plate containing a macrocrack with several neighboring microcracks. Normal loading is considered at distances far away from the cracks. The strip yield is adopted where the plastic zone is assumed to be confined to two narrow strips extending from the ends of a finite length crack while the microcracks are assumed to be elastic. The plastic zone length and crack opening displacement are found from asymptotic solution and compared with finite element solution.     

A two-phase model and related numerical tool for the design of soil structures reinforced by stiff linear inclusions

A two-phase macroscopic model is presented for ascertaining the overall linear elastic behaviour of a soil reinforced by stiff long inclusions, with a particular emphasis on the shear and flexural behaviour of the reinforcements. Based on a minimum principle for the potential energy of any two-phase system, expressed as a function of the kinematically admissible displacement and rotation fields, a finite element formulation is established, resulting in the elaboration of a numerical computer code devoted to the simulation of any reinforced soil structure under plane strain conditions. This finite element code is applied to the evaluation of the settlements experienced by a rigid raft foundation placed on top of a soil reinforced by a group of vertical piles, and subject to combined loading conditions. One of the important conclusions which may be drawn from such a quantitative analysis is that the shear and flexural behaviour of the reinforcing piles, as well as the way these piles are connected to the foundation, play a decisive role in the case of lateral loading.     

Outflow boundary conditions for one-dimensional finite element modeling of blood flow and pressure waves in arteries

Flow and pressure waves, originating due to the contraction of the heart, propagate along the deformable vessels and reflect due to tapering, branching, and other discontinuities. The size and complexity of the cardiovascular system necessitate a “multiscale” approach, with “upstream” regions of interest (large arteries) coupled to reduced-order models of “downstream” vessels. Previous efforts to couple upstream and downstream domains have included specifying resistance and impedance outflow boundary conditions for the nonlinear one-dimensional wave propagation equations and iterative coupling between three-dimensional and one-dimensional numerical methods. We have developed a new approach to solve the one-dimensional nonlinear equations of blood flow in elastic vessels utilizing a space-time finite element method with GLS-stabilization for the upstream domain, and a boundary term to couple to the downstream domain. The outflow boundary conditions are derived following an approach analogous to the Dirichlet-to-Neumann (DtN) method. In the downstream domain, we solve simplified zero/one-dimensional equations to derive relationships between pressure and flow accommodating periodic and transient phenomena with a consistent formulation for different boundary condition types. In this paper, we also present a new boundary condition that accommodates transient phenomena based on a Green’s function solution of the linear, damped wave equation in the downstream domain.     

Use of Scanning Electron Microscope and the Non-local Isotropic Damage Model to Investigate Fracture Process Zone in Notched Concrete Beams

In this study, a Scanning Electron Microscope (SEM) is used to understand the micro level aspect of the Fracture Process Zone (FPZ) in a concrete beam. It is mainly based on the preparation and analyzing samples which are considered as being a very important part of SEM (poor preparation techniques can lead to erroneous diagnosis of the concrete study). Numerically, the fracture of concrete requires the consideration of progressive damage, which is usually modeled by a constitutive law. This latest relies on numerical methods to obtain adequate solutions. It is shown herein that by using the Object Oriented Finite Element Method (OOFEM), obtained results agreed more or less with those of others researchers. On the other side, experimental results compromise those obtained by the use of the non-local isotropic damage model. It is finally proven throughout this study that the FPZ is defined by two parameters: the length and the width.     

Atomistic and continuum studies of stress and strain fields near a rapidly propagating crack in a harmonic lattice

Large-scale atomistic simulations of a mode I crack propagating in a harmonic lattice are presented. The objective of this work is to study the stress and strain fields near a rapidly propagating mode I crack. The asymptotic continuum mechanics solutions of the elastic fields are compared quantitatively with molecular-dynamics simulation results for different crack velocities. It is observed that both atomistic stress and atomistic strain can be successfully related to the corresponding continuum quantities. The study reveals that the atomistic simulation results agree well with the continuum theory predictions, which suggests that the continuum theory can be applied for nano-scale dynamic problems.     

Computational fluid dynamics simulations of respiratory airflow in human nasal cavity and its characteristic dimension study

To study the airflow distribution in human nasal cavity during respiration and the characteristic parameters of nasal structure, three-dimensional, anatomically accurate representations of 30 adult nasal cavity models were recons- tructed based on processed tomography images collected from normal people. The airflow fields in nasal cavities were simulated by fluid dynamics with finite element software ANSYS. The results showed that the difference of human nasal cavity structure led to different airflow distribution in the nasal cavities and variation of the main airstream passing through the common nasal meatus. The nasal resistance in the regions of nasal valve and nasal vestibule accounted for more than half of the overall resistance. The characteristic model of nasal cavity was extracted on the basis of characteristic points and dimensions deduced from the original models. It showed that either the geometric structure or the airflow field of the two kinds of models was similar. The characteristic dimensions were the characteristic parameters of nasal cavity that could properly represent the original model in model studies on nasal cavity.     

The yield surface of viscoelastic and plastic fluids in a vane viscometer

Vane viscometers are often used to investigate the low shear rate properties of plastic fluids. The shear stress is determined by assuming that the material is held in the space between the vane blades so that it behaves like a rigid cylinder. Experimental evidence supports this assumption and the aim of the present study is to model numerically the yield process in a vane rheometer using viscoelastic and plastic fluids. The finite element method has been used to model the behavior of Herschel-Bulkley (Bingham), Casson and viscoelastic (Maxwell type) fluids. The penalty function approach for the pressure approximation and a rotating reference frame are used together with fine meshes containing more than 1300 elements. The results show that for Herschel-Bulkley (Bingham), and Casson fluids a rotating rigid cylinder of fluid is trapped inside the periphery of the vane, the shear stress is uniformly distributed over the surface of the cylinder. Finally a modified second order fluid is used to simulate the viscoelastic behaviour, anticipated to be an intermediate between the elastic deformation and the plastic flow, to provide a more realistic simulation of the yield process about a vane. In this case, contrast with the concentration of the elastic strain rate at the blade tips, a nearly uniform distribution of the plastic shear rate is still found. This implies that the plastic shear always distributes uniformly during the entire yielding process. Evidently the assumption of uniform shear on a rotating cylinder of material occluded in the blades of a vane is a valid and useful model for many types of fluid possessing a yield stress.     

Numerical study on the dynamics of droplet passing through a cylinder obstruction in confined microchannel flow

The droplet dynamics passing through a cylinder obstruction was investigated with direct numerical simulations with FE-FTM (Finite Element-Front Tracking Method). The effect of droplet size and capillary number (Ca) was studied for both Newtonian and viscoelastic fluids. In the case of Newtonian droplet immersed in Newtonian medium, the droplet breakup induced by the geometric hindrance depends on the droplet size. As Ca increases, the short droplets (1.3 times longer than the channel width) break up while passing through the obstruction. However, the breakup does not occur for longer droplets (1.8 times longer than the channel width). When the viscoelastic fluid characterized by the Oldroyd-B model is considered, the Newtonian droplet immersed in viscoelastic medium breaks up into two smaller droplets while passing through the cylinder obstruction with increasing De_m (Deborah number of the medium). We also show that the normal stress difference plays a key role on the droplet breakup and the droplet extension. The normal stress difference is enhanced in the negative wake region due to the droplet flow, which also promotes droplet extension in that region. This numerical study provides information not only on underlying physics of the droplet flows passing through a cylinder obstruction but also on the useful guidelines for microfluidic applications.     

A theory and a simulation capability for the growth of a solid electrolyte interphase layer at an anode particle in a Li-ion battery

A major mechanism for electrochemical aging of Li-ion batteries is the growth of a solid electrolyte interphase (SEI) layer on the surface of anode particles, which leads to capacity fade and also results in a rise in cell resistance. We have formulated a continuum theory for the growth of an SEI layer—a theory which accounts for the generation of the attendant growth stresses. The theory has been numerically implemented in a finite-element program. This simulation capability for SEI growth is coupled with our previously published chemo-mechanical simulation capability for intercalation of Li-ions in electrode particles. Using this new combined capability we have simulated the formation and growth of an SEI layer during cyclic lithiation and delithiation of an anode particle, and predicted the evolution of the growth stresses in the SEI layer. The evolution of the stress state within the SEI layer and at the SEI/anode-particle interface for spherical- and spheroidal-shaped graphite particles is studied. This knowledge of the local interfacial stresses provides a good estimate for the propensity of potential delamination of an SEI layer from an anode particle.     

Incorporation of twinning into a crystal plasticity finite element model: Evolution of lattice strains and texture in Zircaloy-2

A crystal plasticity finite element code is developed to model lattice strains and texture evolution of HCP crystals. The code is implemented to model elastic and plastic deformation considering slip and twinning based plastic deformation. The model accounts for twinning reorientation and growth. Twinning, as well as slip, is considered to follow a rate dependent formulation. The results of the simulations are compared to previously published in situ neutron diffraction data. Experimental results of the evolution of the texture and lattice strains under uniaxial tension/compression loading along the rolling, transverse, and normal direction of a piece of rolled Zircaloy-2 are compared with model predictions. The rate dependent formulation introduced is capable of correctly capturing the influence of slip and twinning deformation on lattice strains as well as texture evolution.     

Modeling of Heat and Mass Transfer in a Fractured Porous Medium

While fractured formations are possibly the most important contributors to the production of oil worldwide, modeling fractured formations with rigorous treatments has eluded reservoir engineers in the past. To date, one of the most commonly used fractured reservoir models remains the one that was suggested by Warren and Root nearly four decades ago. In this paper, a new model for fractures embedded in a porous medium is proposed. The model considers the Navier-Stokes equation in the fracture (channel flow) while using the Brinkman equation for the porous medium. Unlike the previous approach, the proposed model does not require the assumption of orthogonality of the fractures (sugar cube assumption) nor does it impose incorrect boundary conditions for the interface between the fracture and the porous medium. Also, the transfer coefficient between the fracture and matrix interface does not need to be specified, unlike the cases for which Darcy's law is used. In order to demonstrate the usefulness of the approach, a two-dimensional model of a fractured formation is developed and numerical simulation runs conducted.

The proposed model is derived through a series of finite element modeling runs for various cases using the Navier-Stokes equation in the channel while maintaining the Brinkman equation in the porous medium. Various cases studied include different fracture orientations, fracture frequencies, and thermal and solutal constraints. The usefulness of the proposed model in modeling complex formations is discussed. Finally, a series of numerical runs also provided validity of the proposed model for the cases in which thermal and solutal effects are important. Such a study of double diffusive phenomena, coupled with forced convection, in the context of fractured formations has not been reported before.     

基于内聚力理论的二维二次界面单元在ABAQUS中的UEL程序实现

依据有限元理论,结合内聚力模型法则,推导出二维二次粘结界面单元在大位移情况下的数值格式,得到用形函数表示的单元位移模式、载荷向量和刚度矩阵,并进行了离散化。基于ABAQUS软件的自定义扩展模块,编制了相应的用户单元子程序UEL,通过数值算例验证了该程序的准确性和有效性。这一成果能为在ABAQUS软件中开展相关数值研究,以及开发其他类型的内聚力界面有限单元提供思路和参考。     

不同轮轨接触模型在车桥耦合振动中的比较

以工程实例为研究对象,建立了整车-整桥系统耦合振动数值分析模型。考虑车轮的跳轨和挤密情况,建立了单边弹簧-阻尼系统弹性轮轨接触模型。采用基于多体系统动力学和有限元法结合的联合仿真技术,计算了两种轮轨接触时动车组列车以不同车速通过大跨度连续桥梁的耦合振动响应。数值计算结果表明:两种轮轨接触模型的桥梁动力响应比较接近;列车的横向轮轨力、轮重减载率和脱轨系数相差较大,当速度为350km/h时,横向轮轨力增大了46.5%,轮重减载率增大了130.8%,脱轨系数增大了24.66%;用单边-弹簧阻尼系统弹性轮轨接触模型更符合实际。     

Evaluation of factors and current approaches related to computerized design of tillage tools: a review

The objectives of this paper are to evaluate the factors that are involved in the tillage process, and to explore the potential approaches for the computer-aided design of tillage tools. An overview related to the dynamic effect on the performance of tillage operations has been conducted. Compared with the analytical methods, the finite element method (FEM) has some advantages for the computerized design of tillage tools. The artificial neural networks (ANN) may be useful for the integrated evaluation of tillage performance with multi-objectives. ANN can be employed for simulation of a dynamic constitutive model and identification of soil conditions for agricultural soils. The integral approach of ANN analysis with FEM is found to be promising for optimizing design of tillage tools.     

Euler–Lagrange/DEM simulation of wood gasification in a bubbling fluidized bed reactor

We present an Euler–Lagrange method for the simulation of wood gasification in a bubbling fluidized bed. The gas phase is modeled as a continuum using the 2D Navier–Stokes equations and the solid phase is modeled by a Discrete Element Method (DEM) using a soft-sphere approach for the particle collision dynamic. Turbulence is included via a Large-Eddy approach using the Smagorinsky sub-grid model. The model takes into account detailed gas phase chemistry, zero-dimensional modeling of the pyrolysis and gasification of each individual particle, particle shrinkage, and heat and mass transfer between the gas phase and the particulate phase. We investigate the influence of wood feeding rate and compare exhaust gas compositions and temperature results obtained with the model against experimental data of a laboratory scale bubbling fluidized bed reactor.