Characterizing the influence of stress-induced microcracks on the laboratory strength and fracture development in brittle rocks using a finite-discrete element method-micro discrete fracture network FDEM-μDFN approach

Heterogeneity is an inherent component of rock and may be present in different forms including mineralheterogeneity, geometrical heterogeneity, weak grain boundaries and micro-defects. Microcracks areusually observed in crystalline rocks in two forms： natural and stress-induced; the amount of stressinducedmicrocracking increases with depth and in-situ stress. Laboratory results indicate that thephysical properties of rocks such as strength, deformability, P-wave velocity and permeability areinfluenced by increase in microcrack intensity. In this study, the finite-discrete element method （FDEM）is used to model microcrack heterogeneity by introducing into a model sample sets of microcracks usingthe proposed micro discrete fracture network （mDFN） approach. The characteristics of the microcracksrequired to create mDFN models are obtained through image analyses of thin sections of Lac du Bonnetgranite adopted from published literature. A suite of two-dimensional laboratory tests including uniaxial,triaxial compression and Brazilian tests is simulated and the results are compared with laboratory data.The FDEM-mDFN models indicate that micro-heterogeneity has a profound influence on both the mechanicalbehavior and resultant fracture pattern. An increase in the microcrack intensity leads to areduction in the strength of the sample and changes the character of the rock strength envelope. Spallingand axial splitting dominate the failure mode at low confinement while shear failure is the dominantfailure mode at high confinement. Numerical results from simulated compression tests show thatmicrocracking reduces the cohesive component of strength alone, and the frictional strength componentremains unaffected. Results from simulated Brazilian tests show that the tensile strength is influenced bythe presence of microcracks, with a reduction in tensile strength as microcrack intensity increases. Theimportance of microcrack heterogeneity in reproducing a bi-linear or S-shape failure envelope and itseffects on the mechanisms leading to spalling damage near an underground opening are also discussed.

Construction of n-sided polygonal spline element using area coordinates and B-net method

In general, triangular and quadrilateral elements are commonly applied in two-dimensional finite element methods. If they are used to compute polycrystalline materials, the cost of computation can be quite significant. Polygonal elements can do well in simulation of the materials behavior and provide greater flexibility for the meshing of complex geometries. Hence, the study on the polygonal element is a very useful and necessary part in the finite element method. In this paper, an n-sided polygonal element based on quadratic spline interpolant, denoted by PS2 element, is presented using the triangular area coordinates and the B-net method. The PS2 element is conforming and can exactly model the quadratic field. It is valid for both convex and non-convex polygonal element, and insensitive to mesh distortions. In addition, no mapping or coordinate transformation is required and thus no Jacobian matrix and its inverse are evaluated. Some appropriate examples are employed to evaluate the performance of the proposed element.

Finite element analysis of the effect of micro-pore defect on linear friction welding of medium carbon steel

Micro-pore is a very common material defect. In the present paper, the temperature fields of medium carbon steel joints with and without micro-pore defect during linear friction welding （LFW） were investigated by using finite element method. The effect of micro-pore defect on the axial shortening of joints during LFW was examined. The x- and y-direction displacements of micro-pore during the LFW process were also studied. In addition, the shape of micro-pore after LFW was observed. The heat conducted from the weld inteace to the specimen interior. The fluctuation range of the temperature curves for the joint with micro-pore is larger than that without micro-pore. Position of micro-pore changes with the change of the friction time. The circular shape of micro-pore becomes oval after welding.

Two 8-node quadrilateral spline elements by B-net method

Isoparametric quadrilateral elements are widely used in the finite element method, but the accuracy of the isoparametric quadrilateral elements will drop obviously deteriorate due to mesh distortions. Spline functions have some properties of simplicity and conformality. Two 8-node quadrilateral elements have been developed using the trian- gular area coordinates and the B-net method, which can ex- actly model the quadratic field for both convex and concave quadrangles. Some appropriate examples are employed to evaluate the performance of the proposed elements. The nu- merical results show that the two spline elements can obtain solutions which are highly accurate and insensitive to mesh distortions.

A hybrid vertex-centered finite volume/element method for viscous incompressible flows on non-staggered unstructured meshes

This paper proposes a hybrid vertex-centered fi- nite volume/finite element method for solution of the two di- mensional （2D） incompressible Navier-Stokes equations on unstructured grids. An incremental pressure fractional step method is adopted to handle the velocity-pressure coupling. The velocity and the pressure are collocated at the node of the vertex-centered control volume which is formed by join- ing the centroid of cells sharing the common vertex. For the temporal integration of the momentum equations, an im- plicit second-order scheme is utilized to enhance the com- putational stability and eliminate the time step limit due to the diffusion term. The momentum equations are discretized by the vertex-centered finite volume method （FVM） and the pressure Poisson equation is solved by the Galerkin finite el- ement method （FEM）. The momentum interpolation is used to damp out the spurious pressure wiggles. The test case with analytical solutions demonstrates second-order accuracy of the current hybrid scheme in time and space for both veloc- ity and pressure. The classic test cases, the lid-driven cavity flow, the skew cavity flow and the backward-facing step flow, show that numerical results are in good agreement with the published benchmark solutions.

Numerical simulation and response analysis of microspherical focused logging in inclined micro-fractured formation

Tight and unconventional reservoirs have become the focus with the progress of petroleum exploration and development.Micro-fractures in these reservoirs can effectively improve reservoir permeability,and well-developed micro-fractures can serve to directly improve productivity.Compared with the centered electrical well logging method,the Micro Spherical Focused Logging(MSFL)is more suitable for microfracture identification due to its high resolution and near borehole wall measuring method.In this study,an anisotropic model is used to depict micro-fractured formation.First,a forward model with microfractured formation,borehole,logging instrument and surrounding rock is established.Subsequently,MSFL responses under different micro-fracture porosity,resistivity,dip angle and borehole radius,are calculated based on the finite element method(FEM).Finally,the MSFL responses under different microfracture parameters are analyzed with the response laws clarified.

Investigation of Micro-mechanical Response of Asphalt Mixtures by a Three-dimensional Discrete Element Model

The micro-mechanical response of asphalt mixtures was studied using the discrete element method. The discrete element sample of stone mastic asphalt was generated first and the vehicle load was applied to the sample. A user-written program was coded with the FISH language in PFC3 D to extract the contact forces within the sample and the displacements of the particles. Then, the contact forces within the whole sample, in asphalt mastic, in coarse aggregates and between asphalt mastic and coarse aggregates were investigated. Finally, the movement of the particles in the sample was analyzed. The sample was divided into 15 areas and a figure was drawn to show how the balls move in each area according to the displacements of the balls in each area. The displacements of asphalt mastic balls and coarse aggregates were also analyzed. The experimental results explain how the asphalt mixture bears vehicle load and the potential reasons why the rutting forms from a micro-mechanical view.

铺层方式对VIHPS制备的GO-CF混杂增强复合材料弯曲性能的影响

本文首先采用真空浸渗热压工艺(VIHPS)制备了不同铺层方式的氧化石墨烯-碳纤维(GO-CF)混杂增强复合材料,并对该材料进行三点弯曲试验,再通过扫描电子显微镜(SEM)表征其微观断口形貌。然后使用Digmat软件的FE模块建立微-纳跨尺度的代表性体积单元(RVE)模型,将该模型导入ANSYS Workbench软件中进行三点弯曲仿真模拟。复合材料三点弯曲试验结果表明,随着沿试样宽度方向(Y向)纤维铺层的增加,复合材料的弯曲强度从433.2 MPa降低至12.7 MPa,弯曲模量从13.8 GPa降低至0.2 GPa,它们分别降低了97.07%和98.55%。最后通过对比试验与仿真结果发现,两种结果的变化趋势基本一致,且仿真结果与试验结果接近。复合材料在弯曲变形过程中,受到拉伸和压缩的共同作用,其中导致复合材料失效的方式主要包括基体开裂、界面脱粘和纤维断裂。

基于Micro-CT成像的听骨链有限元建模方法

目的利用Micro-CT扫描数据建立中耳听骨链三维有限元模型。方法用Micro-CT扫描成人颞骨标本,将获得的图像数据通过MicroView、Mimics等软件进行三维重建。再将模型转入其有限元分析(finite elementan alysis,FEA)模块的Remesh环境中进行调整、细化及面网格优化,通过SOLIDEWORK软件转换为实体网格。结果初步形成的三维几何模型可较清晰地辨别鼓室腔、听小骨和内耳系统,但部分图像不同程度存在噪点。最终形成了FEA软件可识别的成人中耳听骨链三维有限元网格模型,网格划分后的听骨链有限元模型中,完整听骨链的节点数降低,由原来包括805个不合理节点在内的12498个节点降低到2050个,在此基础上建立的有限元模型共1350个8节点四面体单元。结论结合Micro-CT技术及Mimics软件的三维建模方法可以快速获得较精确的听骨链三维数据,是建立听骨链三维有限元模型的有效途径。

运用Micro-CT技术建立上颌第一磨牙根管系统的有限元模型

目的:采用有限元的方法建立上颌磨牙根管系统的三维有限元模型。方法:采用Micro-CT扫描,应用Mimics8.1、Geomagic studio8、Unigraphics NX及Ansys9.0等逆向工程软件分别得到上颌第一磨牙根管系统三维外形、实体模型和有限元模型。结果:重建了根管系统的三维有限元模型,网格质量好,具有高度的精确度,并由此证实模型的可行性。结论:将Micro-CT扫描技术、数字图像处理技术与有限元方法结合起来,建立出高精度的上颌第一磨牙根管系统的三维有限元模型,为口腔生物力学提供了研究手段。

Micro-CT技术结合逆向工程软件和CAD软件在建立模拟根管治疗术后牙体组织三维有限元模型中的应用

目的:通过Micro-CT技术和逆向工程软件及CAD软件建立牙齿根管治疗术后的三维有限元模型。方法:通过Micro-CT对获得的完整人下颌第一磨牙进行扫描,将获得的图片导入Mimics 10.01软件进行初步拟合,之后导入Geomagic 10软件中精修细化获得NUBS曲面。在Pro/Engineer Wildfire 4.0软件中根据根管器械的形态学参数建立模拟根管预备后的根管形态。最后,将精修细化的曲面文件和在Pro/E中建立的根管系统在ANSYS前处理软件Workbench中进行布尔运算和切割,并且进行属性定义和网格划分。结果:建立了根管锥度为0.02、根尖直径为0.35mm的下颌第一磨牙的精细有限元模型。结论:Micro-CT技术结合逆向工程软件和CAD软件可有效建立出牙齿模拟根管治疗后精细三维有限元模型。

Theoretical analysis of the sound absorption characteristics of periodically stiffened micro-perforated plates

The vibro-acoustic responses and sound absorption characteristics of two kinds of periodically stiffened micro-perforated plates are analyzed theoretically. The connected periodical structures of the stiffened plates can be ribs or block-like structures. Based on fundamental acoustic formulas of the micro-perforated plate of Maa and Takahashi, semi-analytical models of the vibrating stiffened plates are developed in this paper. Approaches like the space harmonic method, Fourier transforms and finite element method （FEM） are adopted to investigate both kinds of the stiffened plates. In the present work, the vibro-acoustic responses of micro-perforated stiffened plates in the wavenumber space are expressed as functions of plate displacement amplitudes. After approximate numerical solutions of the amplitudes, the vibration equations and sound absorption coefficients of the two kinds of stiffened plates in the physical space are then derived by employing the Fourier inverse transform. In numerical examples, the effects of some physical parameters, such as the perforation ratio, incident angles and periodical distances etc., on the sound absorption performance are examined. The proposed approaches are also validated by comparing the present results with solutions of Takahashi and previous studies of stiffened plates. Numerical results indicate that the flexural vibration of the plate has a signif- icant effect on the sound absorption coefficient in the water but has little influence in the air.

An explicit finite volume element method for solving characteristic level set equation on triangular grids

Level set methods are widely used for predicting evolutions of complex free surface topologies,such as the crystal and crack growth,bubbles and droplets deformation,spilling and breaking waves,and two-phase flow phenomena.This paper presents a characteristic level set equation which is derived from the two-dimensional level set equation by using the characteristic-based scheme.An explicit finite volume element method is developed to discretize the equation on triangular grids.Several examples are presented to demonstrate the performance of the proposed method for calculating interface evolutions in time.The proposed level set method is also coupled with the Navier-Stokes equations for two-phase immiscible incompressible flow analysis with surface tension.The Rayleigh-Taylor instability problem is used to test and evaluate the effectiveness of the proposed scheme.

The effect of agglomerate on micro-structural evolution in solid-state sintering

Discrete element method （DEM） is used in the present paper to simulate the microstructural evolution of a planar layer of copper particles during sintering. Formation of agglomerates and the effect of their rearrangement on densification are mainly focused on. Comparing to the existing experimental observations, we find that agglomerate can form spontaneously in sintering and its rearrangement could accelerate the densification of compacts. Snapshots of numerical simulations agree qualitatively well with experimental observations. The method could be readily extended to investigate the effect of agglomerate on sintering in a three- dimensional model, which should be very useful for understanding the evolution of microstructure of sintering systems.

Pseudo-beam method for compressive buckling characteristics analysis of space inflatable load-carrying structures

This paper extends Le van＇s work to the case of nonlinear problem and the complicated configuration. The wrinkling stress distribution and the pressure effects are also included in our analysis. Pseudo-beam method is presented based on the inflatable beam theory to model the inflatable structures as a set of inflatable beam elements with a prestressed state. In this method, the discretized nonlinear equations are given based upon the virtual work principle with a 3-node Timoshenko＇s beam model. Finite element simulation is performed by using a 3-node BEAM189 element incorporating ANSYS nonlinear program. The pressure effect is equivalent included in our method by modifying beam element cross-section parameters related to pressure. A benchmark example, the bending case of an inflatable cantilever beam, is performed to verify the accuracy of our proposed method. The comparisons reveal that the numerical results obtained with our method are close to open published analytical and membrane finite element results. The method is then used to evaluate the whole buckling and the loadcarrying characteristics of an inflatable support frame subjected to a compression force. The wrinkling stress and region characteristics are also shown in the end. This method gives better convergence characteristics, and requires much less computation time. It is very effective to deal with the whole load-carrying ability analytical problems for large scale inflatable structures with complex configuration.

Finite Element Analysis of Elastic-Plastic Contact Mechanic Considering the Effect of Contact Geometry and Material Properties

Each surface of roughness has different shape of asperity which is modeled with various shapes of analytical models. In this paper, the differences among various models of shape of asperity investigate using the Finite Element Method (FEM) and various analytical models. The contact stresses in rough surfaces are calculated analytically using various asperity shape models. Finite element analysis is also carried out assuming three types of material properties namely, the linear, the elastic-perfect plastic and the elastic-nonlinear hardening. The analytical results are compared with the results obtained by the finite element method. The results illustrate for using a deterministic approach which the numerical models are suitable. In hertz model, the result of force is very big in interface of causing deformation plastic, while Model Zhao has almost same result with FEM nonlinear property model. It is observed that the results obtained from Zhao’s model are generally in a better agreement with the results obtained from various finite element models especially in elastic-plastic and plastic zones, hence it may be concluded that Zhao’s model can be used for analyzing the rough surfaces in contact mechanics.

Parametric study on single shot peening by dimensional analysis method incorporated with finite element method

Shot peening is a widely used surface treatment method by generating compressive residual stress near the surface of metallic materials to increase fatigue life and re- sistance to corrosion fatigue, cracking, etc. Compressive re- sidual stress and dent profile are important factors to eval- uate the effectiveness of shot peening process. In this pa- per, the influence of dimensionless parameters on maximum compressive residual stress and maximum depth of the dent were investigated. Firstly, dimensionless relations of pro- cessing parameters that affect the maximum compressive residual stress and the maximum depth of the dent were de- duced by dimensional analysis method. Secondly, the in- fluence of each dimensionless parameter on dimensionless variables was investigated by the finite element method. Fur- thermore, related empirical formulas were given for each di- mensionless parameter based on the simulation results. Fi- nally, comparison was made and good agreement was found between the simulation results and the empirical formula, which shows that a useful approach is provided in this pa- per for analyzing the influence of each individual parameter.

3D random Voronoi grain-based models for simulation of brittle rock damage and fabric-guided micro-fracturing

A grain-based distinct element model featuring three-dimensional （3D） Voronoi tessellations （randompoly-crystals） is proposed for simulation of crack damage development in brittle rocks. The grainboundaries in poly-crystal structure produced by Voronoi tessellations can represent flaws in intact rockand allow for numerical replication of crack damage progression through initiation and propagation ofmicro-fractures along grain boundaries. The Voronoi modelling scheme has been used widely in the pastfor brittle fracture simulation of rock materials. However the difficulty of generating 3D Voronoi modelshas limited its application to two-dimensional （2D） codes. The proposed approach is implemented inNeper, an open-source engine for generation of 3D Voronoi grains, to generate block geometry files thatcan be read directly into 3DEC. A series of Unconfined Compressive Strength （UCS） tests are simulated in3DEC to verify the proposed methodology for 3D simulation of brittle fractures and to investigate therelationship between each micro-parameter and the model＇s macro-response. The possibility of numericalreplication of the classical U-shape strength curve for anisotropic rocks is also investigated innumerical UCS tests by using complex-shaped （elongated） grains that are cemented to one another alongtheir adjoining sides. A micro-parameter calibration procedure is established for 3D Voronoi models foraccurate replication of the mechanical behaviour of isotropic and anisotropic （containing a fabric） rocks. 2014 Institute of Rock and Soil Mechanics, Chinese Academy of Sciences. Production and hosting byElsevier B.V. All rights reserved.

Finite Element Simulation of the Polymeric Liquid Flow in Coat-hanger Die

A three - dimenslonal finite element model is developed to deal with the polymeric liquid flow in coat - hanger die. This model is used to predict the flow behavior of the 2％ CMC/watsr solution in the coat - hanger die with linearly taper manifolds and its validity is experimentally verified quantitatively and qualitatively by using Laser Doppler Velocimetry and Particle Image Velocimetry respectively.

基于SBFEM的剪切型裂纹动态开裂模拟

剪切型断裂是岩土工程中常见的破坏模式,了解剪切破坏机理并准确预测剪切型裂纹的萌生、扩展过程对保障工程结构的安全性与稳定性具有重要意义.文章建立了基于比例边界有限元法(scaled boundary finite element methods,SBFEM)和非局部宏-微观损伤模型的剪切型裂纹动态开裂模拟方法,定义了基于偏应变概念的物质点对的正伸长量,可作为预测剪切型裂纹扩展行为的动态开裂准则,一点的损伤定义为该点影响域范围内连接的物质键损伤的加权平均值,而物质键的损伤则与基于偏应变概念的物质点对的正伸长量相关联,并引入能量退化函数建立结构域几何拓扑损伤与能量损失之间的关系,将拓扑损伤与应力应变联系起来,通过能量退化函数修正了SBFEM的刚度系数矩阵,得到了子域在损伤状态下的刚度矩阵,推导了考虑结构损伤的SBFEM动力控制方程,采用Newmark隐式算法对控制方程进行时间离散.最后,通过3个典型算例验证了建议的模型可较好地模拟剪切型断裂问题,能够很好地捕捉剪切型裂纹的扩展路径,并得到较为准确的载荷-位移曲线.