Multiscale modeling of heterogeneous propellants from particle packing to grain failure using a surface-based cohesive approach

In the present work,a computational framework is established for multiscale modeling and analysis ofsolid propellants.A packing algorithm,considering the ammonium perchlorate(AP) and aluminum(Al) particles asspheres or discs is developed to match the size distributionand volume fraction of solid propellants.A homogenizationtheory is employed to compute the mean stress and strainof a representative volume element(RVE).Using the meanresults,a suitable size of RVE is decided.Without considering the interfaces between particles and matrix,several numerical simulations of the relaxation of propellants are performed.The relaxation effect and the nonlinear mechanicalbehavior of propellants which are dependent on the appliedloads are discussed.A new technology named surface-basedcohesive behavior is proposed to describe the phenomenonof particle dewetting consisting of two ingredients:a damageinitiation criterion and a damage evolution law.Several examples considering contact damage behavior are computedand also nonlinear behavior caused by damaged interfaces isdiscussed in this paper.Furthermore the effects of the critical contact stress,initial contact stiffness and contact failuredistance on the damaged interface model have been studied.     

Nonlinear constitutive force model selection,update and uncertainty quantification for periodically sequential impact applications

Nonlinear Dynamics - Contact–impact events frequently occur between solid contact interfaces in complex dynamic mechanical systems and engineering applications. Impact force models are used...     

Superquadric DEM-SPH coupling method for interaction between non-spherical granular materials and fluids

The research on the coupling method of non-spherical granular materials and fluids aims to predict the particle–fluid interaction in this study. A coupling method based on superquadric elements is developed to describe the interaction between non-spherical solid particles and fluids. The discrete element method (DEM) and the smoothed particle hydrodynamics (SPH) are adopted to simulate granular materials and fluids. The repulsive force model is adopted to calculate the coupling force and then a contact detection method is established for the interaction between the superquadric element and the fluid particle. The contact detection method captures the shape of superquadric element and calculates the distance from the fluid particle to the surface of superquadric element. Simulation cases focusing on the coupling force model, energy transfer, and large-scale calculations have been implemented to verify the validity of the proposed coupling method. The coupling force model accurately represents the water entry process of a spherical solid particle, and reasonably reflects the difference of solid particles with different shapes. In the water entry process of multiple solid particles, the total energy of the water entry process of multiple solid particles tends to be stable. The collapse process of the partially submerged granular column is simulated and analyzed under different parameters. Therefore, this coupling method is suitable to simulate fluid–particle systems containing solid particles with multiple shapes.     

Unified Model of Drainage and Imbibition in 3D Fractionally Wet Porous Media

We develop a grain-based model for capillarity controlled displacement within 3D fractionally wet porous media. The model is based on a novel local calculation of the position of stable fluid–fluid interfaces in contact with multiple spherical grains of arbitrary contact angles. The interface is assumed to be locally spherical between bulk phases; the interface is assumed to be toroidal between pairs of grains (surfaces of pendular rings). Because the calculation of interface position is entirely local and grain-based, it provides a single, generalized, geometric basis for computing pore-filling events during drainage as well as imbibition with both Melrose events (merging of two interfaces) and Haines events (geometric instability). The model is validated against a series of drainage/imbibition experiments (oil/water) on fractionally wet porous media prepared by mixing oil-wet grains with water-wet grains.     

Two-dimensional numerical analyses of double conforming contacts with effect of curvature

Pinned connections and journal bearings with cylindrical conforming contacts are widely used as fundamental building blocks in machines and civil structures of civilization. It is critical to determine their stress fields and contact performances under desired configurations in order to successfully design them. This paper presents a numerical model to deal with such contacts, particularly with the general configuration–double interfaces. Necessary formulae are developed and numerical procedures are described in this paper. Efficient numerical methods, i.e., the discrete convolution-fast Fourier transformation (DC-FFT) method and the conjugate gradient method (CGM), are implemented in the algorithm. Validations are conducted against the Persson’s results and finite element results and demonstrate excellent agreement. This contact analysis can be useful for design engineers to evaluate stress and contact pressure distribution. Furthermore, this efficient method of determining radial displacement can be applied in an elasto-hydrodynamic lubrication analysis.     

结构爆炸毁伤的浸没多介质有限体积物质点法

结构在爆炸载荷作用下的毁伤现象涉及强非线性激波、固体结构极端变形和破坏破碎、强流固耦合, 给数值计算方法带来了极大的困难与挑战. 针对结构爆炸毁伤问题, 建立了浸没多介质有限体积物质点法(iMMFV-MPM), 采用基于黎曼求解器的多介质有限体积法(MMFVM)模拟爆炸产物和空气的多介质流体, 采用物质点法(MPM)模拟固体结构, 并将提出的基于拉格朗日乘子的连续力浸没边界法(lg-CFIBM)扩展到多介质流体中以处理流固耦合边界条件. 该算法在每个时间步严格满足流固耦合界面处的速度边界条件及动量守恒方程, 不需要重构流固耦合界面, 能够有效地模拟近场爆炸下爆炸产物与结构的相互作用、激波与结构的相互作用和演化以及结构的动态断裂和拓扑变化. 利用iMMFV-MPM对近场爆炸下方形钢筋混凝土靶板的失效模式、外爆载荷下建筑物的毁伤现象以及多腔室内爆炸试验进行了模拟, 模拟结果与相关实验数据吻合良好, 验证了所建立的流固耦合算法的有效性及精度.      

An elastoplastic multi-level damage model for ductile matrix composites considering evolutionary weakened interface

An elastoplastic multi-level damage model considering evolutionary weakened interface is developed in this work to predict the effective elastoplastic behavior and multi-level damage evolution in particle reinforced ductile matrix composites (PRDMCs). The elastoplastic multi-level damage model is micromechanically derived on the basis of the ensemble-volume averaging procedure and the first-order effects of eigenstrains. The Eshelby’s tensor for an ellipsoidal inclusion with slightly weakened interface [Qu, J., 1993a. Eshelby tensor for an elastic inclusion with slightly weakened interfaces. Journal of Applied Mechanics 60 (4), 1048–1050; Qu, J., 1993b. The effect of slightly weakened interfaces on the overall elastic properties of composite materials. Mechanics of Materials, 14, 269–281] is adopted to model particles having mildly or severely weakened interface, and a multi-level damage model [Lee, H.K., Pyo, S.H., in press. Multi-level modeling of effective elastic behavior and progressive weakened interface in particulate composites. Composites Science and Technology] in accordance with the Weibull’s probabilistic function is employed to describe the sequential, progressive weakened interface in the composites. Numerical examples corresponding to uniaxial, biaxial and triaxial tension loadings are solved to illustrate the potential of the proposed micromechanical framework. A series of parametric analysis are carried out to investigate the influence of model parameters on the progression of weakened interface in the composites. Furthermore, the present prediction is compared with available experimental data in the literature to verify the proposed elastoplastic multi-level damage model.     

Strong discontinuity analysis of a class of anisotropic continuum damage constitutive models – Part II: Concrete material application

On the basis of the strong discontinuity analysis, a discrete model expressed in terms of traction vector-displacement jump has been constructed from a continuous model expressed in terms of stress–strain law. In the first part of the paper, this approach has been extended to a class of anisotropic continuum damage constitutive models [1]. In this second part of the paper, the proposed class of discrete anisotropic damage constitutive models is particularized. More precisely, a micromechanical-based anisotropic damage constitutive model is derived. This model accounts in a natural manner for particular crack families orientation. The aims of this paper are (i) to illustrate the capabilities of the proposed discrete enhanced model in reproducing the induced anisotropy appearing in quasi-brittle materials when cracking and (ii) to assess the numerical robustness of the time integration scheme. For this purpose, two numerical examples at the material point level are exposed after a brief presentation of the time integration scheme. The correspondence between the continuous and the discrete model as well as the induced anisotropy features are emphasized.     

Dynamic coupling of fluid and structural mechanics for simulating particle motion and interaction in high speed compressible gas particle laden flow

In this work, structural finite element analyses of particles moving and interacting within high speed compressible flow are directly coupled to computational fluid dynamics and heat transfer analyses to provide more detailed and improved simulations of particle laden flow under these operating conditions. For a given solid material model, stresses and displacements throughout the solid body are determined with the particle–particle contact following an element to element local spring force model and local fluid induced forces directly calculated from the finite volume flow solution. Plasticity and particle deformation common in such a flow regime can be incorporated in a more rigorous manner than typical discrete element models where structural conditions are not directly modeled. Using the developed techniques, simulations of normal collisions between two 1 mm radius particles with initial particle velocities of 50–150 m/s are conducted with different levels of pressure driven gas flow moving normal to the initial particle motion for elastic and elastic–plastic with strain hardening based solid material models. In this manner, the relationships between the collision velocity, the material behavior models, and the fluid flow and the particle motion and deformation can be investigated. The elastic–plastic material behavior results in post collision velocities 16–50% of their pre-collision values while the elastic-based particle collisions nearly regained their initial velocity upon rebound. The elastic–plastic material models produce contact forces less than half of those for elastic collisions, longer contact times, and greater particle deformation. Fluid flow forces affect the particle motion even at high collision speeds regardless of the solid material behavior model. With the elastic models, the collision force varied little with the strength of the gas flow driver. For the elastic–plastic models, the larger particle deformation and the resulting increasingly asymmetric loading lead to growing differences in the collision force magnitudes and directions as the gas flow strength increased. The coupled finite volume flow and finite element structural analyses provide a capability to capture the interdependencies between the interaction of the particles, the particle deformation, the fluid flow and the particle motion.     

Advection upwinding splitting method to solve a compressible two‐fluid model

In this study, the advection upwinding splitting method (AUSM) is modified for the resolution of two‐phase mixtures with interfaces. The compressible two‐fluid model proposed by Saurel and Abgrall is chosen as the model equations. Dense and dilute phases are described in terms of the volume fraction and equations of state to represent multi‐phase mixtures. Test cases involving an air–water shock tube, water faucet, and dilute particulate turbulent flows through a 90° bend are used to verify the current work. It is shown that the AUSM based on flux differences (AUSMD) contains the mechanism to correctly capture the contact discontinuity and interfaces between phases. In addition, a successful application to dilute particulate turbulence flows by the AUSMD is demonstrated. Copyright © 2001 John Wiley & Sons, Ltd.     

The solid variational principles of the discrete from—The variational principles of the discrete analysis by the finite element method

This paper suggests a new solid variational principle of discrete form. Basing on the true case of the discrete analysis by the finite element method and considering the variable boundaries of the elements and the unknown functions of piecewise approximation, the unknown functions have various discontinuities at the interfaces between successive element.Thus, we have used mathematical technique of variable boundary with discontinuity of the unknown functions, based on the conditions that the first variation vanishes immediately, to establish the solid variation principles of discrete form. It generalizes the classical and non-classical variational principles. Successive equations that have to be satisfied by the unknown functions are the convergency necessary conditions for the finite elements method (including conforming and non-conforming). They expand that convergency necessary conditions of the compatibility conditions in the internal interfaces.     

Statistical damage theory of 2D lattices: Energetics and physical foundations of damage parameter

The paper presents an in-depth analysis of two-dimensional disordered lattices of statistical damage mechanics for the study of quasi-brittle materials. The strain energy variation in correspondence to damage formation is thoroughly examined and all the different contributions to the net energy changes are identified and analyzed separately. We demonstrate that the introduction of a new defect in the microstructure produces a perturbation of the microscopic random fields according to a principle of maximum energy dissipation. A redistribution parameter η is introduced to measure the load redistribution capability of the microstructure. This parameter can be estimated from simulation data of detailed models. This energetic framework sets the stage for the investigation of the statistical foundations of the damage parameter as well as the damage localization. Logical statistical arguments are developed to derive two analytical models (a maximum field model and a mean field one) for the estimate of the damage parameter via a bottom-up approach that relates the applied load to the microstructural disorder. Simulation data provided input to the statistical models as well as the means of validation. Simulated tensile tests of honeycomb lattices with mechanical disorder demonstrate that long-range interactions amongst sets of microcracks with different orientations play a fundamental role already in damage nucleation as well as in the homogeneous–heterogeneous transition. A functional “hierarchy of sets” of grain boundaries, based on their orientation in relation to the applied load, seems to emerge from this study. Results put in evidence the ability of discrete models of capturing seamlessly the damage anisotropy. The ideas exposed inhere should be useful to develop a full rational model for disordered lattices and, later, to extend the approach to discrete models with solid elements. The findings suggest that statistical damage mechanics might aid in the quest of reliable and physically sound constitutive relations of damage, even in synergy with micromechanics.     

Preventing numerical oscillations in the flux‐split based finite difference method for compressible flows with discontinuities,II

Problems in the characteristic‐wise flux‐split based finite difference method when compressible flows with contact discontinuities or material interfaces are computed were presented and analyzed. The current analysis showed the following: (i) Even with the local characteristic decomposition technique, numerical errors could be caused by point‐wise flux vector splitting (FVS) methods, such as the Steger–Warming FVS or the van Leer FVS. Therefore, the Lax–Friedrichs type FVS method is required. (ii) If the isobars of a material are vertical lines, the combination of using the local characteristic decomposition and the global Lax–Friedrichs FVS can avoid velocity and pressure oscillations of contact discontinuities in this material for weighted essentially non‐oscillatory (WENO) schemes. (iii) For problems with material interfaces, the quasi‐conservative approach can be realized using characteristic‐wise flux‐split based finite difference WENO schemes if nonlinear WENO schemes in genuinely nonlinear characteristic fields can be guaranteed to be the same and the decomposition equation representing material interfaces is discretized properly. Copyright © 2015 John Wiley & Sons, Ltd.     

High-resolution algorithms of sharp interface treatment for compressible two-phase flows

In this paper, a kind of arbitrary high order derivatives (ADER) scheme based on the generalised Riemann problem is proposed to simulate multi-material flows by a coupling ghost fluid method. The states at cell interfaces are reconstructed by interpolating polynomials which are piece-wise smooth functions. The states are treated as the equivalent of the left and right states of the Riemann problem. The contact solvers are extrapolated in the vicinity of contact points to facilitate ghost fluids. The numerical method is applied to compressible flows with sharp discontinuities, such as the collision of two fluids of different physical states and gas–liquid two-phase flows. The numerical results demonstrate that unexpected physical oscillations through the contact discontinuities can be prevented effectively and the sharp interface can be captured efficiently.     

Micromechanics-based elastic damage modeling of particulate composites with weakened interfaces

A micromechanical framework is proposed to predict the effective elastic behavior and weakened interface evolution of particulate composites. The Eshelby’s tensor for an ellipsoidal inclusion with slightly weakened interface [Qu, J., 1993a. Eshelby tensor for an elastic inclusion with slightly weakened interfaces. Journal of Applied Mechanics 60 (4), 1048–1050; Qu, J., 1993b. The effect of slightly weakened interfaces on the overall elastic properties of composite materials. Mechanics of Materials 14, 269–281] is adopted to model spherical particles having imperfect interfaces in the composites and is incorporated into the micromechanical framework. Based on the Eshelby’s micromechanics, the effective elastic moduli of three-phase particulate composites are derived. A damage model is subsequently considered in accordance with the Weibull’s probabilistic function to characterize the varying probability of evolution of weakened interface between the inclusion and the matrix. The proposed micromechanical elastic damage model is applied to the uniaxial, biaxial and triaxial tensile loadings to predict the various stress–strain responses. Comparisons between the present predictions with other numerical and analytical predictions and available experimental data are conducted to assess the potential of the present framework.     

固体火箭发动机柔性接头摆动密封可靠性研究

柔性接头由弹性件与增强件交替粘接而成,是固体火箭发动机进行推力矢量控制的重要装置,因而柔性接头的摆动密封性能对固体火箭发动机而言至关重要.为研究固体火箭发动机柔性接头摆动过程中的密封可靠性,以内聚力模型作为粘接界面的本构模型,通过计算柔性接头各界面的损伤情况及界面间的接触应力,并定义界面节点单元间、界面间以及柔性接头的...     

Discrete and continuous models for in plane loaded random elastic brickwork

The aim of this paper is to study non-periodic masonries – typical of historical buildings – by means of a perturbation approach and to evaluate the effect of a random perturbation on the elastic response of a periodic masonry wall. The random masonry is obtained starting from a periodic running bond pattern. A random perturbation on the horizontal positions of the vertical interfaces between the blocks which form the masonry wall is introduced. In this way, the height of the blocks is uniform, while their width in the horizontal direction is random. The perturbation is limited such as each block has still exactly 6 neighboring blocks. In a first discrete model, the blocks are modeled as rigid bodies connected by elastic interfaces (mortar thin joints). In other words, masonry is seen as a “skeleton” in which the interactions between the rigid blocks are represented by forces and moments which depend on their relative displacements and rotations. A second continuous model is based on the homogenization of the discrete model. Explicit upper and lower bounds on the effective elastic moduli of the homogenized continuous model are obtained and compared to the well-known effective elastic moduli of the regular periodic masonry. It is found that the effective moduli are not very sensitive to the random perturbation (less than 10%). At the end, the Monte Carlo simulation method is used to compare the discrete random model and the continuous model at the structural level (a panel undergoing in plane actions). The randomness of the geometry requires the generation of several samples of size L of the discrete masonry. For a sample of size L, the structural discrete problem is solved using the same numerical procedure adopted in [Cecchi, A., Sab, K., 2004. A comparison between a 3D discrete model and two homogenized plate models for periodic elastic brickwork, International Journal of Solids Structures 41 (9–10), 2259–2276] and the average solution over the samples gives an estimation which depends on L. As L increases, an asymptotic limit is reached. One issue is to find the minimum size for L and to compare the asymptotic average solution to the one obtained from the continuous homogenized model.     

Finite element analysis of flow,heat transfer,and free interfaces in an electron-beam vaporization system for metals

A numerical analysis is made of the liquid flow and energy transport in a system to evaporate metals. The energy from an electron-beam heats an axisymmetric metal disk supported by a water-cooled platform. Metal evaporates from the surface of a hot pool of liquid which is surrounded by a shell of its own solid. Flow in the pool is strongly driven by temperature-induced buoyancy and capillary forces, and is located in the transition region between laminar and turbulent flow. The evaporation rate is strongly influenced by the locations of the free boundaries. A modified finite element method is used to calculate the steady state flow and temperature fields coupled with the interface locations. The mesh is structured with spines that stretch and pivot as the interfaces move. The discretized equations are arranged in an ‘arrow’ matrix and are solved using the Newton–Raphson method. The electron-beam power and platform contact resistance are varied for cases involving the evaporation of aluminum. The results reveal the interaction of liquid flow, heat transfer and free interfaces. © 1998 John Wiley & Sons, Ltd.