岩石类材料直接拉伸破坏的FDEM数值模拟

有限元-离散元耦合算法(FDEM)在岩石裂纹扩展模拟中广泛应用,尤其在直接拉伸模拟中用于参数校核。然而,模拟结果通常受到拉伸速率和网格尺寸的影响。本文首先阐述FDEM基本原理;然后,建立直接拉伸数值计算模型,探究拉伸速率、网格尺寸及高径比对破坏模式和力学参数的影响;最后,给出直接拉伸模拟中建议采用的拉伸速率和网格尺寸。研究结果表明:随着拉伸速率增加,试样抗拉强度、破坏程度和峰值应变率增加,而平均弹性模量减小;在较低拉伸速率下(v≤0.01 m/s),试样仅产生1条主拉伸裂缝;当拉伸速率较高时,有分支裂缝出现在主拉伸裂缝上侧。网格尺寸对抗拉强度、峰值应变率和平均弹性模量的影响较小,但对破坏模式和破坏程度的影响较显著。随着高径比增加,试样的抗拉强度、破坏程度和峰值应变率增大,但平均弹性模量呈相反趋势,并且高径比对破坏程度的影响较为显著。建议在直接拉伸模拟中使用0.01 m/s的拉伸速率,并将网格尺寸限制在1.5 mm以内。

基于Weibull分布的非均质隧洞围岩破裂碎胀FDEM数值模拟研究

岩石材料具有非均质特征,矿物颗粒和胶结材料强度在统计学上服从Weibull分布。基于有限元-离散元耦合数值模拟方法(FDEM)提出了三角形单元弹性模量和四边形节理单元强度服从Weibull分布的参数赋值流程,随后研究了非均质岩样在单轴压缩、三轴压缩和巴西劈裂下的力学特征和破坏行为,最后对非均质隧洞围岩在高应力下的破裂碎胀大变形机制进行了模拟研究。结果表明:提出的随机参数赋值方法可建立不同弹性模量三角形单元数目服从Weibull分布、且空间分布具有随机性的FDEM数值模型;非均质岩样单轴抗压强度、三轴抗压强度、巴西劈裂抗拉强度、等效粘聚力和等效内摩擦角随着非均质度m的增大,均呈指数增大,并向均质材料逼近;高应力非均质隧洞围岩随着m的增大,裂隙网络形态变化不大,均以共轭剪切破裂为主,并伴随少量拉伸裂隙,但围岩破坏程度和裂隙扩展范围呈指数降低。

基于精确缩尺和粗粒化的颗粒材料FDEM模拟方法

颗粒材料具有非连续、离散性等特征,在进行数值模拟时面临着较大的计算压力。通过将精确缩尺准则和粗粒化方法引入到连续-离散耦合(combined finite-discrete element method,FDEM)方法中,旨在为加速基于FDEM的颗粒材料数值模拟提供一种解决方案。基于精确缩尺和粗粒化等理论,推导了FDEM中应遵循的精确缩尺准则,在此基础上分别进行了等粒径颗粒体系及二元颗粒体系的三轴剪切数值试验。试验结果表明,在未引入精确缩尺准则时,粗粒化模型表现的力学响应特征会发生改变,结果出现失真,因此必须对粗粒化模型参数进行修正。引入精确缩尺准则后,粗粒化模型的力学响应特征会得到补正。试验结果论证了FDEM引入精确缩尺准则和粗粒化方法的有效性,即能在近似原始颗粒体系的条件下大幅度提升采用FDEM进行颗粒材料数值模拟的计算效率。基于数值试验结果进行了宏细观力学分析,宏观应力变形和细观接触力相互映证,揭示了精确缩尺和粗粒化方法的细观力学机理。

A coupled cryogenic thermo-hydro-mechanical model for frozen medium:Theory and implementation in FDEM

This paper presents the development of a coupled modeling approach to simulate cryogenic thermo-hydro-mechanical(THM)processes associated with a freezing medium,which is then implemented in the combined finite-discrete element method code(FDEM)for multi-physics simulation.The governing equations are deduced based on energy and mass conservation,and static equilibrium equations,considering water/ice phase change,where the strong couplings between multi-fields are supplemented by critical coupling parameters(e.g.unfrozen water content,permeability,and thermal conductivity).The proposed model is validated against laboratory and field experiments.Results show that the cryogenic THM model can well predict the evolution of strongly coupled processes observed in frozen media(e.g.heat transfer,water migration,and frost heave deformation),while also capturing,as emergent properties of the model,important phenomena(e.g.latent heat,cryogenic suction,ice expansion and distinct three-zone distribution)caused by water/ice phase change at laboratory and field scales,which are difficult to be all revealed by existing THM models.The novel modeling framework presents a gateway to further understanding and predicting the multi-physical coupling behavior of frozen media in cold regions.

A new thermomechanical coupled FDEM model for geomaterials considering continuum-discontinuum transitions

A new thermomechanical(TM)coupled finite-discrete element method(FDEM)model,incorporating heat conduction,thermal cracking,and contact heat transfer,has been proposed for both continuous and discontinuous geomaterials.This model incorporates a heat conduction model that can accurately calculate the thermal field in continuousediscontinuous transition processes within a finite element framework.A modified contact heat transfer model is also included,which accounts for the entire contact area of discrete bodies.To align with the finite strain theory utilized in the FDEM mechanics module,the TM coupling module in the model is based on the multiplicative decomposition of the deformation gradient.The proposed model has been applied to various scenarios,including heat conduction in both continuous and discontinuous media during transient states,thermal-induced strain and stress,and thermal cracking conditions.The thermal field calculation model and the TM coupling model have been validated by comparing the numerical results with experiment findings and analytical solutions.These numerical cases demonstrate the reliability of the proposed model convincingly,making it suitable for use across a wide range of continuous and discontinuous media.

A Review on Sources,Extractions and Analysis Methods of a Sustainable Biomaterial:Tannins

Condensed and hydrolysable tannins are non-toxic natural polyphenols that are a commercial commodity industrialized for tanning hides to obtain leather and for a growing number of other industrial applications mainly to substitute petroleum-based products.They are a definite class of sustainable materials of the forestry industry.They have been in operation for hundreds of years to manufacture leather and now for a growing number of applications in a variety of other industries,such as wood adhesives,metal coating,pharmaceutical/medical applications and several others.This review presents the main sources,either already or potentially commercial of this forestry by-materials,their industrial and laboratory extraction systems,their systems of analysis with their advantages and drawbacks,be these methods so simple to even appear primitive but nonetheless of proven effectiveness,or very modern and instrumental.It constitutes a basic but essential summary of what is necessary to know of these sustainable materials.In doing so,the review highlights some of the main challenges that remain to be addressed to deliver the quality and economics of tannin supply necessary to fulfill the industrial production requirements for some materials-based uses.

河下浅埋厚煤层采动覆岩裂隙-涌水量时空演化的FDEM-CFD耦合分析

河下浅埋厚煤层采动覆岩裂隙分布与涌水量是工作面安全生产的决定性因素之一。数值模拟是二者重要预测方法,其合理性的关键在于建立岩体破坏-裂隙流体耦合理论及相应的模拟方法。以路家村矿15404工作面为研究背景,构建拉/剪应力下非贯通裂隙开裂、贯通裂隙的法向与切向本构关系,并基于二相流质量守恒、动量守恒与状态方程,结合增强的浸没边界算法识别流-固界面,通过流体体积法实现裂隙内流体自由面的追踪和重构。在此基础上形成FDEM-CFD河下采煤覆岩裂隙与涌水量预测数值模型耦合程序,通过相邻工作面地面钻孔冲洗液消耗量观测法验证导水裂隙带发育高度,以及采用大井法理论对涌水量结果进行对比。结果表明:采动岩体破坏-裂隙流体耦合理论及相应的FDEM-CFD程序,可数值实现河下浅埋厚煤层采动覆岩裂隙形成,以及裂隙内流体运移过程。当工作面推进至80~120 m时覆岩内形成贯通地表的导水裂隙。招山河水的主要下泄路径为工作面后方8~20 m位置处的导水裂隙,其斜向采空区、倾角65°~72°。模拟所得采空区涌水量为18.78 m^(3)/h,与大井法计算结果接近。上述成果在路家村矿区得到初步应用,为进一步开展河下浅埋厚煤层防治水工程提供理论支撑。

Drilling-based measuring method for the c-φ parameter of rock and its field application

The technology of drilling tests makes it possible to obtain the strength parameter of rock accurately in situ. In this paper, a new rock cutting analysis model that considers the influence of the rock crushing zone(RCZ) is built. The formula for an ultimate cutting force is established based on the limit equilibrium principle. The relationship between digital drilling parameters(DDP) and the c-φ parameter(DDP-cφ formula, where c refers to the cohesion and φ refers to the internal friction angle) is derived, and the response of drilling parameters and cutting ratio to the strength parameters is analyzed. The drillingbased measuring method for the c-φ parameter of rock is constructed. The laboratory verification test is then completed, and the difference in results between the drilling test and the compression test is less than 6%. On this basis, in-situ rock drilling tests in a traffic tunnel and a coal mine roadway are carried out, and the strength parameters of the surrounding rock are effectively tested. The average difference ratio of the results is less than 11%, which verifies the effectiveness of the proposed method for obtaining the strength parameters based on digital drilling. This study provides methodological support for field testing of rock strength parameters.

A comparison study on structure-function relationship of polysaccharides obtained from sea buckthorn berries using different methods:antioxidant and bile acid-binding capacity

In this study,the structural characters,antioxidant activities and bile acid-binding ability of sea buckthorn polysaccharides(HRPs)obtained by the commonly used hot water(HRP-W),pressurized hot water(HRP-H),ultrasonic(HRP-U),acid(HRP-C)and alkali(HRP-A)assisted extraction methods were investigated.The results demonstrated that extraction methods had significant effects on extraction yield,monosaccharide composition,molecular weight,particle size,triple-helical structure,and surface morphology of HRPs except for the major linkage bands.Thermogravimetric analysis showed that HRP-U with filamentous reticular microstructure exhibited better thermal stability.The HRP-A with the lowest molecular weight and highest arabinose content possessed the best antioxidant activities.Moreover,the rheological analysis indicated that HRPs with higher galacturonic acid content and molecular weight showed higher viscosity and stronger crosslinking network(HRP-C,HRP-W and HRP-U),which exhibited stronger bile acid binding capacity.The present findings provide scientific evidence in the preparation technology of sea buckthorn polysaccharides with good antioxidant and bile acid binding capacity which are related to the structure affected by the extraction methods.

Material point method simulation of hydro-mechanical behaviour in twophase porous geomaterials: A state-of-the-art review

The material point method(MPM)has been gaining increasing popularity as an appropriate approach to the solution of coupled hydro-mechanical problems involving large deformation.In this paper,we survey the current state-of-the-art in the MPM simulation of hydro-mechanical behaviour in two-phase porous geomaterials.The review covers the recent advances and developments in the MPM and their extensions to capture the coupled hydro-mechanical problems involving large deformations.The focus of this review is aiming at providing a clear picture of what has or has not been developed or implemented for simulating two-phase coupled large deformation problems,which will provide some direct reference for both practitioners and researchers.

Volumetric lattice Boltzmann method for pore-scale mass diffusionadvection process in geopolymer porous structures

Porous materials present significant advantages for absorbing radioactive isotopes in nuclear waste streams.To improve absorption efficiency in nuclear waste treatment,a thorough understanding of the diffusion-advection process within porous structures is essential for material design.In this study,we present advancements in the volumetric lattice Boltzmann method(VLBM)for modeling and simulating pore-scale diffusion-advection of radioactive isotopes within geopolymer porous structures.These structures are created using the phase field method(PFM)to precisely control pore architectures.In our VLBM approach,we introduce a concentration field of an isotope seamlessly coupled with the velocity field and solve it by the time evolution of its particle population function.To address the computational intensity inherent in the coupled lattice Boltzmann equations for velocity and concentration fields,we implement graphics processing unit(GPU)parallelization.Validation of the developed model involves examining the flow and diffusion fields in porous structures.Remarkably,good agreement is observed for both the velocity field from VLBM and multiphysics object-oriented simulation environment(MOOSE),and the concentration field from VLBM and the finite difference method(FDM).Furthermore,we investigate the effects of background flow,species diffusivity,and porosity on the diffusion-advection behavior by varying the background flow velocity,diffusion coefficient,and pore volume fraction,respectively.Notably,all three parameters exert an influence on the diffusion-advection process.Increased background flow and diffusivity markedly accelerate the process due to increased advection intensity and enhanced diffusion capability,respectively.Conversely,increasing the porosity has a less significant effect,causing a slight slowdown of the diffusion-advection process due to the expanded pore volume.This comprehensive parametric study provides valuable insights into the kinetics of isotope uptake in porous structures,facilitating the development of porous materials for nuclear waste treatment applications.

Sparse Modal Decomposition Method Addressing Underdetermined Vortex-Induced Vibration Reconstruction Problem for Marine Risers

When investigating the vortex-induced vibration(VIV)of marine risers,extrapolating the dynamic response on the entire length based on limited sensor measurements is a crucial step in both laboratory experiments and fatigue monitoring of real risers.The problem is conventionally solved using the modal decomposition method,based on the principle that the response can be approximated by a weighted sum of limited vibration modes.However,the method is not valid when the problem is underdetermined,i.e.,the number of unknown mode weights is more than the number of known measurements.This study proposed a sparse modal decomposition method based on the compressed sensing theory and the Compressive Sampling Matching Pursuit(Co Sa MP)algorithm,exploiting the sparsity of VIV in the modal space.In the validation study based on high-order VIV experiment data,the proposed method successfully reconstructed the response using only seven acceleration measurements when the conventional methods failed.A primary advantage of the proposed method is that it offers a completely data-driven approach for the underdetermined VIV reconstruction problem,which is more favorable than existing model-dependent solutions for many practical applications such as riser structural health monitoring.

Uncertainties of landslide susceptibility prediction: Influences of random errors in landslide conditioning factors and errors reduction by low pass filter method

In the existing landslide susceptibility prediction(LSP)models,the influences of random errors in landslide conditioning factors on LSP are not considered,instead the original conditioning factors are directly taken as the model inputs,which brings uncertainties to LSP results.This study aims to reveal the influence rules of the different proportional random errors in conditioning factors on the LSP un-certainties,and further explore a method which can effectively reduce the random errors in conditioning factors.The original conditioning factors are firstly used to construct original factors-based LSP models,and then different random errors of 5%,10%,15% and 20%are added to these original factors for con-structing relevant errors-based LSP models.Secondly,low-pass filter-based LSP models are constructed by eliminating the random errors using low-pass filter method.Thirdly,the Ruijin County of China with 370 landslides and 16 conditioning factors are used as study case.Three typical machine learning models,i.e.multilayer perceptron(MLP),support vector machine(SVM)and random forest(RF),are selected as LSP models.Finally,the LSP uncertainties are discussed and results show that:(1)The low-pass filter can effectively reduce the random errors in conditioning factors to decrease the LSP uncertainties.(2)With the proportions of random errors increasing from 5%to 20%,the LSP uncertainty increases continuously.(3)The original factors-based models are feasible for LSP in the absence of more accurate conditioning factors.(4)The influence degrees of two uncertainty issues,machine learning models and different proportions of random errors,on the LSP modeling are large and basically the same.(5)The Shapley values effectively explain the internal mechanism of machine learning model predicting landslide sus-ceptibility.In conclusion,greater proportion of random errors in conditioning factors results in higher LSP uncertainty,and low-pass filter can effectively reduce these random errors.

Numerical manifold method for thermo-mechanical coupling simulation of fractured rock mass

As a calculation method based on the Galerkin variation,the numerical manifold method(NMM)adopts a double covering system,which can easily deal with discontinuous deformation problems and has a high calculation accuracy.Aiming at the thermo-mechanical(TM)coupling problem of fractured rock masses,this study uses the NMM to simulate the processes of crack initiation and propagation in a rock mass under the influence of temperature field,deduces related system equations,and proposes a penalty function method to deal with boundary conditions.Numerical examples are employed to confirm the effectiveness and high accuracy of this method.By the thermal stress analysis of a thick-walled cylinder(TWC),the simulation of cracking in the TWC under heating and cooling conditions,and the simulation of thermal cracking of the SwedishÄspöPillar Stability Experiment(APSE)rock column,the thermal stress,and TM coupling are obtained.The numerical simulation results are in good agreement with the test data and other numerical results,thus verifying the effectiveness of the NMM in dealing with thermal stress and crack propagation problems of fractured rock masses.

A stable implicit nodal integration-based particle finite element method(N-PFEM)for modelling saturated soil dynamics

In this study,we present a novel nodal integration-based particle finite element method(N-PFEM)designed for the dynamic analysis of saturated soils.Our approach incorporates the nodal integration technique into a generalised Hellinger-Reissner(HR)variational principle,creating an implicit PFEM formulation.To mitigate the volumetric locking issue in low-order elements,we employ a node-based strain smoothing technique.By discretising field variables at the centre of smoothing cells,we achieve nodal integration over cells,eliminating the need for sophisticated mapping operations after re-meshing in the PFEM.We express the discretised governing equations as a min-max optimisation problem,which is further reformulated as a standard second-order cone programming(SOCP)problem.Stresses,pore water pressure,and displacements are simultaneously determined using the advanced primal-dual interior point method.Consequently,our numerical model offers improved accuracy for stresses and pore water pressure compared to the displacement-based PFEM formulation.Numerical experiments demonstrate that the N-PFEM efficiently captures both transient and long-term hydro-mechanical behaviour of saturated soils with high accuracy,obviating the need for stabilisation or regularisation techniques commonly employed in other nodal integration-based PFEM approaches.This work holds significant implications for the development of robust and accurate numerical tools for studying saturated soil dynamics.

An inverse analysis of fluid flow through granular media using differentiable lattice Boltzmann method

This study presents a method for the inverse analysis of fluid flow problems.The focus is put on accurately determining boundary conditions and characterizing the physical properties of granular media,such as permeability,and fluid components,like viscosity.The primary aim is to deduce either constant pressure head or pressure profiles,given the known velocity field at a steady-state flow through a conduit containing obstacles,including walls,spheres,and grains.The lattice Boltzmann method(LBM)combined with automatic differentiation(AD)(AD-LBM)is employed,with the help of the GPU-capable Taichi programming language.A lightweight tape is used to generate gradients for the entire LBM simulation,enabling end-to-end backpropagation.Our AD-LBM approach accurately estimates the boundary conditions for complex flow paths in porous media,leading to observed steady-state velocity fields and deriving macro-scale permeability and fluid viscosity.The method demonstrates significant advantages in terms of prediction accuracy and computational efficiency,making it a powerful tool for solving inverse fluid flow problems in various applications.

Robustness Study and Superior Method Development and Validation for Analytical Assay Method of Atropine Sulfate in Pharmaceutical Ophthalmic Solution

Background: The robustness is a measurement of an analytical chemical method and its ability to contain unaffected by little with deliberate variation of analytical chemical method parameters. The analytical chemical method variation parameters are based on pH variability of buffer solution of mobile phase, organic ratio composition changes, stationary phase (column) manufacture, brand name and lot number variation;flow rate variation and temperature variation of chromatographic system. The analytical chemical method for assay of Atropine Sulfate conducted for robustness evaluation. The typical variation considered for mobile phase organic ratio change, change of pH, change of temperature, change of flow rate, change of column etc. Purpose: The aim of this study is to develop a cost effective, short run time and robust analytical chemical method for the assay quantification of Atropine in Pharmaceutical Ophthalmic Solution. This will help to make analytical decisions quickly for research and development scientists as well as will help with quality control product release for patient consumption. This analytical method will help to meet the market demand through quick quality control test of Atropine Ophthalmic Solution and it is very easy for maintaining (GDP) good documentation practices within the shortest period of time. Method: HPLC method has been selected for developing superior method to Compendial method. Both the compendial HPLC method and developed HPLC method was run into the same HPLC system to prove the superiority of developed method. Sensitivity, precision, reproducibility, accuracy parameters were considered for superiority of method. Mobile phase ratio change, pH of buffer solution, change of stationary phase temperature, change of flow rate and change of column were taken into consideration for robustness study of the developed method. Results: The limit of quantitation (LOQ) of developed method was much low than the compendial method. The % RSD for the six sample assay of developed method was 0.4% where the % RSD of the compendial method was 1.2%. The reproducibility between two analysts was 100.4% for developed method on the contrary the compendial method was 98.4%.

Cascading multi-segment rupture process of the 2023 Turkish earthquake doublet on a complex fault system revealed by teleseismic P wave back projection method

In this study,the vertical components of broadband teleseismic P wave data recorded by China Earthquake Network are used to image the rupture processes of the February 6th,2023 Turkish earthquake doublet via back projection analysis.Data in two frequency bands(0.5-2 Hz and 1-3 Hz)are used in the imaging processes.The results show that the rupture of the first event extends about 200 km to the northeast and about 150 km to the southwest,lasting~90 s in total.The southwestern rupture is triggered by the northeastern rupture,demonstrating a sequential bidirectional unilateral rupture pattern.The rupture of the second event extends approximately 80 km in both northeast and west directions,lasting~35 s in total and demonstrates a typical bilateral rupture feature.The cascading ruptures on both sides also reflect the occurrence of selective rupture behaviors on bifurcated faults.In addition,we observe super-shear ruptures on certain fault sections with relatively straight fault structures and sparse aftershocks.

Hybrid Strategy of Partitioned and Monolithic Methods for Solving Strongly Coupled Analysis of Inverse and Direct Piezoelectric and Circuit Coupling

The inverse and direct piezoelectric and circuit coupling are widely observed in advanced electro-mechanical systems such as piezoelectric energy harvesters.Existing strongly coupled analysis methods based on direct numerical modeling for this phenomenon can be classified into partitioned or monolithic formulations.Each formulation has its advantages and disadvantages,and the choice depends on the characteristics of each coupled problem.This study proposes a new option:a coupled analysis strategy that combines the best features of the existing formulations,namely,the hybrid partitioned-monolithic method.The analysis of inverse piezoelectricity and the monolithic analysis of direct piezoelectric and circuit interaction are strongly coupled using a partitioned iterative hierarchical algorithm.In a typical benchmark problem of a piezoelectric energy harvester,this research compares the results from the proposed method to those from the conventional strongly coupled partitioned iterative method,discussing the accuracy,stability,and computational cost.The proposed hybrid concept is effective for coupled multi-physics problems,including various coupling conditions.

Gas Chromatographic Method for Identification and Quantification of Commonly Used Residual Solvents in Pharmaceuticals Products

Background: Impurities are not expected in the final pharmaceutical products. All impurities should be regulated in both drug substances and drug products in accordance with pharmacopeias and ICH guidelines. Three different types of impurities are generally available in the pharmaceutical’s product specification: organic impurities, inorganic impurities, and residual solvents. Residual solvents are organic volatile chemicals used or generated during the manufacturing of drug substances or drug products. Purpose: The aim of this study is to develop a cost-effective gas chromatographic method for the identification and quantification of some commonly used solvents—methanol, acetone, isopropyl alcohol (IPA), methylene chloride, ethyl acetate, tetrahydrofuran (THF), benzene, toluene, and pyridine—in pharmaceutical product manufacturing. This method will be able to identify and quantify the multiple solvents within a single gas chromatographic procedure. Method: A gas chromatography (GC) equipped with a headspace sampler and a flame ionization detector, and a column DB 624, 30-meter-long × 0.32-millimeter internal diameter, 1,8 μm-thick, Brand-Agilent was used to develop this method. The initial GC oven temperature was 40°C and held for 5 minutes. It was then increase to 80˚C at a rate of 2˚C per minute, followed by a further increase to 225˚C at a rate of 30˚C per minute, with a final hold at 225˚C for 10 minutes. Nitrogen was used as a carrier gas at a flow rate of 1.20 mL per minute. Dimethyl sulfoxide (DMSO) was selected as sample solvent. Results: The developed method is precise and specific. The percent RSD for the areas of six replicate injections of this gas chromatographic method was within 10.0 and the recovery result found within 80.0% to 120.0%.