Numerical simulation on the seismic performance of retrofitted masonry walls based on the combined finite-discrete element method

Due to the long construction life,improper design methods,brittle material properties and poor construction techniques,most existing masonry structures do not perform well during earthquakes.The retrofitting method using an external steel-meshed mortar layer is widely used to retrofit existing masonry buildings.Assessing the seismic performance of masonry walls reinforced by an external steel-meshed mortar layer reasonably and effectively is a difficult subject in the research field of masonry structures.Based on the combined finite-discrete elements method,the numerical models of retrofitted brick walls with four different masonry mortar strengths by an external mortar layer are established.The shear strength of mortar and the contact between the retrofitted mortar layer and the brick blocks are discussed in detail.The failure patterns and load-displacement curves of the retrofitted brick walls were obtained by applying low cycle reciprocating loads to the numerical model,and the bearing capacity and the failure mechanism of the retrofitted walls were obtained by comparing the failure patterns,ultimate bearing capacity,deformability and other aspects with the tests.This study provides a basis for improving the seismic strengthening design method of masonry structures and helps to better assess the seismic performance of masonry structures after retrofitting.

Numerical simulation of hydraulic fracturing and associated microseismicity using finite-discrete element method

Hydraulic fracturing （HF） technique has been extensively used for the exploitation of unconventional oiland gas reservoirs. HF enhances the connectivity of less permeable oil and gas-bearing rock formationsby fluid injection, which creates an interconnected fracture network and increases the hydrocarbonproduction. Meanwhile, microseismic （MS） monitoring is one of the most effective approaches to evaluatesuch stimulation process. In this paper, the combined finite-discrete element method （FDEM） isadopted to numerically simulate HF and associated MS. Several post-processing tools, includingfrequency-magnitude distribution （b-value）, fractal dimension （D-value）, and seismic events clustering,are utilized to interpret numerical results. A non-parametric clustering algorithm designed specificallyfor FDEM is used to reduce the mesh dependency and extract more realistic seismic information.Simulation results indicated that at the local scale, the HF process tends to propagate following the rockmass discontinuities; while at the reservoir scale, it tends to develop in the direction parallel to themaximum in-situ stress. 2014 Institute of Rock and Soil Mechanics, Chinese Academy of Sciences. Production and hosting byElsevier B.V. All rights reserved.

GPGPU-parallelised hybrid finite-discrete element modelling of rock chipping and fragmentation process in mechanical cutting

Mechanical cutting provides one of the most flexible and environmentally friendly excavation methods.It has attracted numerous efforts to model the rock chipping and fragmentation process,especially using the explicit finite element method(FEM) and bonded particle model(BPM),in order to improve cutting efficiency.This study investigates the application of a general-purpose graphic-processing-unit parallelised hybrid finite-discrete element method(FDEM) which enjoys the advantages of both explicit FEM and BPM,in modelling the rock chipping and fragmentation process in the rock scratch test of mechanical rock cutting.The input parameters of FDEM are determined through a calibration procedure of modelling conventional Brazilian tensile and uniaxial compressive tests of limestone,A series of scratch tests with various cutting velocities,cutter rake angles and cutting depths is then modelled using FDEM with calibrated input parameters.A few cycles of cutter/rock interactions,including their engagement and detachment process,are modelled for each case,which is conducted for the first time to the best knowledge of the authors,thanks to the general purpose graphic processing units(GPGPU) parallelisation.The failure mechanism,cutting force,chipping morphology and effect of various factors on them are discussed on the basis of the modelled results.Finally,it is concluded that GPGPU-parallelised FDEM provides a powerful tool to further study rock cutting and improve cutting efficiencies since it can explicitly capture different fracture mechanisms contributing to the rock chipping as well as chip formation and the separation process in mechanical cutting.Moreover,it is concluded that chipping is mostly owed to the mix-mode Ⅰ-Ⅱ fracture in all cases although mode Ⅱ cracks and mode Ⅰ cracks are the dominant failures in rock cutting with shallow and deep cutting depths,respectively.The chip morphology is found to be a function of cutter velocdty,cutting depth and cutter rake angle.

Research on Simulation Methods of Electric Field Intensity on Surface of 10 kV Cable Joint

The electric field intensity (EFI) is important characteristic quantity for evaluating the internal insulation state of cable joints. Based on finite element method, this paper proposes two EFI research methods, field-circuit coupling method and equivalent circuit method. The average EFI of the inner surface of the outer semi-conducting shield can be calculated from the current in the measuring circuit. The relative error between these two methods is about 15%, which roughly proves the consistency of the two methods. Further practical application research enables online monitoring of cable joints.

Numerical Simulation of Projection Welding Processes for Door Hinge of Automobile Based on Coupled Fields Analysis

Projection welding is a variation of electric resistance welding with the dynamic changes of the flow paths for heat and electrical properties with changing temperature caused by the large plastic deformation collapse of projection. As the joint type between the auto door hinge and the inner plate, projection welding may bring welding distortions and would affect the assembly quality of auto body. A comprehensive electric-thermal-mechanical numerical simulation was performed to quantitatively simulate the processes of projection welding by using a coupled finite element method. The mechanism of projection collapse and the formation process of nugget were discussed and good conclusions have been achieved comparing with the test results.

Stability and reliability of pit slopes in surface mining combined with underground mining in Tonglushan mine

Slope stability is of critical importance in the process of surface-underground mining combination. The influence of underground mining on pit slope stability was mainly discussed, and the self-stabilization of underground stopes was also studied. The random finite element method was used to analyze the probability of the rock mass stability degree of both pit slopes and underground stopes. Meanwhile, 3D elasto-plastic finite element method was used to research into the stress, strain and rock mass failure resulting from mining. The results of numerical simulation indicate that the mining of the underground test stope has certain influence on the stability of the pit slope, but the influence is not great. The safety factor of pit slope is decreased by 0.06, and the failure probability of the pit slope is increased by 1.84%. In addition, the strata yielding zone exists around the underground test stope. The results basically conform to the information coming from the field monitoring.

Numerical Simulation of Interaction Between Laminar Flow and Elastic Sheet

A numerical simulation of the interaction between laminar flow with low Reynolds number and a highly flexible elastic sheet is presented. The mathematical model for the simulation includes a three-dimensional finitevolume based fluid solver for incompressible viscous flow and a combined finite-discrete element method for the three-dimensional deformation of solid. An immersed boundary method is used to couple the simulation of fluid and solid. It is implemented through a set of immersed boundary points scattered on the solid surface. These points provide a deformable solid wall boundary for the fluid by adding body force to Navier-Stokes equations. The force from the fluid is also obtained for each point and then applied on the boundary nodes of the solid. The vortex-induced vibration of the highly flexible elastic sheet is simulated with the established mathematical model. The simulated results for both swing pattern and oscillation frequency of the elastic sheet in low Reynolds number flow agree well with experimental data.

Probabilistic analysis of embankment slope stability in frozen ground regions based on random finite element method

Prediction on the coupled thermal-hydraulic fields of embankment and cutting slopes is essential to the assessment on evolution of melting zone and natural permafrost table, which is usually a key factor for permafrost embankment design in frozen ground regions. The prediction may be further complicated due to the inherent uncertainties of material properties. Hence, stochastic analyses should be conducted. Firstly, Karhunen-Loeve expansion is applied to attain the random fields for hydraulic and thermal conductions. Next, the mixed-form modified Richards equation for mass transfer （i.e., mass equation） and the heat transport equation for heat transient flow in a variably saturated frozen soil are combined into one equation with temperature unknown. Furthermore, the finite element formulation for the coupled thermal-hydraulic fields is derived. Based on the random fields, the stochastic finite element analyses on stability of embankment are carried out. Numerical results show that stochastic analyses of embankment stability may provide a more rational picture for the distribution of factors of safety （FOS）, which is definitely useful for embankment design in frozen ground regions.

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

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

基于微观流固耦合的超深层致密砂岩气藏应力敏感性分析

超深层致密砂岩气藏具有强应力敏感性,目前常用的研究方法包括压力脉冲实验法和实时在线CT扫描法两种,压力脉冲实验法不能揭示产生应力敏感性的微观机理,而实时在线CT扫描法也无法模拟深部地层高压、高应力的条件。为解决应力敏感性实验研究的不足,基于离散单元法与管道网络模型建立了微观流固耦合算法,编制了模拟器,并对模拟器力学计算和流固耦合模块的正确性进行了验证,分析了应力大小、加载方向和孔隙压力对岩心渗透率的影响,最后从微观上揭示了超深层致密砂岩气藏的应力敏感性机理。研究结果表明:①应力通过增加与之垂直方向上喉道两侧的法向压力,减小喉道的水力半径,进而降低储层的渗透率;②较高的孔隙压力能够阻碍岩石颗粒在应力作用下的移动,从而减缓了孔隙和喉道的变形,使模型保持较高的渗透率;③致密砂岩气藏的渗透率受到应力和地层压力的共同控制,并且具有各向异性,在垂直于最小主应力方向上形成渗透率较大的优势通道;④异常高压阻碍了地应力的压实作用,有利于保护储层孔隙,使地层有较好的储集性能和较高的渗透率。结论认为,根据离散元法结合孔隙网络模型建立的流固耦合方法可为理解超深层致密砂岩应力敏感性提供理论参考,并为超深层致密砂岩气藏的科学高效开发提供指导。

Anisotropic shearing mechanism of Kangding slate:Experimental investigation and numerical analysis

The shear mechanical behavior is regarded as an essential factor affecting the stability of the surrounding rocks in underground engineering.The shear strength and failure mechanisms of layered rock are significantly affected by the foliation angles.Direct shear tests were conducted on cubic slate samples with foliation angles of 0°,30°,45°,60°,and 90°.The effect of foliation angles on failure patterns,acoustic emission(AE)characteristics,and shear strength parameters was analyzed.Based on AE characteristics,the slate failure process could be divided into four stages:quiet period,step-like increasing period,dramatic increasing period,and remission period.A new empirical expression of cohesion for layered rock was proposed,which was compared with linear and sinusoidal cohesion expressions based on the results made by this paper and previous experiments.The comparative analysis demonstrated that the new expression has better prediction ability than other expressions.The proposed empirical equation was used for direct shear simulations with the combined finite-discrete element method(FDEM),and it was found to align well with the experimental results.Considering both computational efficiency and accuracy,it was recommended to use a shear rate of 0.01 m/s for FDEM to carry out direct shear simulations.To balance the relationship between the number of elements and the simulation results in the direct shear simulations,the recommended element size is 1 mm.

孔洞形状对层状岩石力学特性影响的FDEM数值模拟研究

为揭示孔洞形状对层状岩石力学性质和破坏模式的影响,采用有限元-离散元耦合数值模拟方法(FDEM)开展系列数值模拟研究。首先,阐述FDEM模拟层状岩石的基本原理,然后,开展不同层理倾角完整试样的单轴压缩模拟,并与试验结果进行对比验证。最后,对含有圆形、椭圆形、三角形、矩形和方形孔洞的层状岩石试样(β为0°、30°、45°、60°和90°)进行单轴压缩模拟研究。研究结果表明:对于不同孔洞形状的试样,抗压强度和试样破坏程度随层理倾角的增大呈V字型变化趋势,且均在层理倾角为0°时抗压强度取得最大值,在90°时破坏程度最为严重;孔洞的存在严重削弱了试样的力学性能,且削弱幅度与孔洞形状密切相关,其中方形和圆形孔洞对抗压强度的削弱能力最弱。将含不同形状孔洞试样的破坏模式随层理倾角的变化主要分为穿层理面拉剪混合破坏(β=0°)、沿层理面与穿层理面混合拉剪破坏(β=30°)、沿层理面剪切破坏(β=45°和60°)和沿层理面拉剪劈裂破坏(β=90°)。

二尖瓣运动的流固耦合仿真

目的 建立包含左心和血液的二尖瓣理想模型,用流固耦合仿真研究二尖瓣在血流中的运动特性。方法 基于解剖学参数建立二尖瓣、左心和血液模型,流固耦合仿真采用有限元结合浸没边界法,使用有限元软件LS-DYNA模拟二尖瓣运动,获取形态学、力学和血液动力学参数,并与结构仿真结果进行对比。结果 两种仿真下二尖瓣形态学结果差异较大,流固耦合结果与超声影像吻合。流固耦合仿真和结构仿真的瓣叶应力分布结果一致,最大第一主应力分别为1.48、1.53 MPa,相对误差为3.27%。左心流场有较为复杂的涡旋结构,舒张期二尖瓣最大流速为1.02 m/s,与健康人体生理数据(0.89±0.15) m/s相吻合。结论 二尖瓣流固耦合仿真可以获取更贴近于生理的形态学结果;流固耦合仿真可以提供临床诊断不可或缺的流场参数信息;单研究瓣叶应力分布问题时,结构仿真更高效。

基于4H-SiC的低横向灵敏度压阻式加速度计设计仿真研究

文中基于4H-SiC设计了一种开槽式低横向灵敏度压阻式加速度计,通过分析支撑梁位置及开槽尺寸对加速度计性能的影响确定了相关结构参数,并对该结构采用有限元方法进行了热-力耦合仿真,确定了该结构在600℃温度场中施加不同载荷时的理论电压输出。研究结果表明:该结构相较于已有的双端四梁结构X轴方向的横向灵敏度由3.3%降至2.6%,Y轴方向横向灵敏度由29.9%降至1.5%,传感器芯片在600℃高温环境下输入载荷与输出电压呈线性关系,100 g载荷下的理论输出电压为3.65 mV。

基于ABAQUS的钛合金管板焊接过程模拟仿真

针对蒸汽发生单元钛合金管板焊接过程中焊接变形预测与控制问题,开展了焊接过程有限元仿真研究。基于顺序耦合的热弹塑性有限元法,采用ABAQUS对典型管板结构进行了不同焊接顺序下的焊接仿真分析,利用校核后的组合热源定义热流密度函数,计算得到焊接温度场、应力场及变形。结果显示,采用高斯热源与椭球体组合热源可以较为精确模拟熔池形貌,管板焊缝区域存在应力集中,且应力峰值位于管板被夹持区域,大小为631.2 MPa,超过材料屈服极限。单个环焊缝在不同路径下的应力分布因受到周围焊缝热输入的影响而出现复杂波动。采用从上至下的对称焊接顺序可以有效减小焊接变形量,最大变形量降幅可达24.3%,为后续工艺优化提供了参考依据。

某小口径火炮弹丸挤进过程动力学特性分析

为研究某小口径火炮不同坡膛结构下弹丸挤进过程的动力学特性,建立了考虑挤进阻力变化的内弹道方程组并通过vuamp子程序引入到有限元求解中,实现了挤进过程中内弹道计算与弹丸挤进有限元仿真耦合求解。通过对弹丸全膛运动的仿真分析验证了模型的可行性,随后得到并分析不同坡膛锥度、膛线深度及阴线宽度下挤进过程中弹丸的速度、加速度、挤进阻力和膛压曲线及变化影响规律。结果表明,坡膛锥度对挤进过程中弹丸的动力学参数影响不是单调的,而膛线深度及阴线宽度的改变则单调地影响弹丸动力学参数,对优化坡膛结构及弹带设计具有一定的指导意义。

Collaborative optimization of linear vibrating screen screening efficiency and dynamic response stability based on coupled DEM-MBK simulation

The dynamic response stability of the vibrating screen is an important factor affecting the screening effect and the structural performance of the vibrating screen.In this paper,to improve dynamic response stability and screening efficiency,we optimized the configuration of linear screening process parameters based on the co-optimization method with dual objectives via the virtual experiment.Firstly,a coupled DEM-MBK simulation model was established according to the dynamics of linear screen,and the dy-namic response law of the screen machine under material impact was investigated.Secondly,the quantitative index of dynamic response stability according to the time-domain characteristics of the centroid amplitude was established.The trend and significance of three types of screening process pa-rameters,including excitation,damping and structure,on the screening efficiency and dynamic response stability were analyzed through virtual orthogonal experiments.Finally,a parameter configuration scheme to achieve co-optimization was proposed based on the comprehensive balance method.The virtual experiment results show that the screening efficiency and dynamic response stability of the proposed scheme are improved by 3.28%and 49.07%,respectively,compared with the empirical parameter configuration.Obviously,the co-optimization method can maintain high screening efficiency and dynamic response stability at the same time,which is beneficial to improve the service life of the screen surface and screen body.

热障涂层热力耦合条件下的应力仿真计算

有限元模拟是研究热障涂层(thermal barrier coatings,TBCs)界面处热生长氧化层(thermally grown oxides,TGO)应力演变的有效手段之一,可为探索TBCs失效机制提供理论支撑。采用先进的热障涂层多因素耦合设备对热障涂层圆管试样模拟发动机工况进行热力耦合循环实验和热循环实验。利用有限元软件ABAQUS对包含真实初始TGO形貌的热障涂层进行有限元建模,分析实验过程中TGO的应力和变形规律。结果表明:在不考虑TGO及界面开裂的情况下,无论是热力耦合模型还是热循环模型,随着循环次数的增加,Mises应力均增加;加热过程中TGO受拉,冷却过程中TGO受压;加热过程的应力均远小于冷却至室温时的应力。经过相同的循环次数,热力耦合模型中的应力值均高于热循环模型中的应力值。经过20,45,70次热循环后,冷却至室温时,TGO应力分别达到2.85,3.65,3.55 GPa,而经过相同次数的热力耦合循环后,冷却至室温时,相同位置的TGO应力分别达到4.01,5.0,4.81 GPa。与热循环相比,在热力耦合条件下,经过相同循环次数冷却至室温时的TGO应力显著增加。

寒旱区湿陷性土地基渠道衬砌结构改造的数值模拟研究

寒旱区输水渠道防渗衬砌结构极易发生冻胀破坏而影响灌区供水安全和稳定运行。为提高湿陷性黄土地基渠道的抗冻胀性能,以宁夏固海扬水灌区衬砌渠道为例,通过分析原型渠道冻胀破坏特征,基于提高防渗体抗冻胀变形能力的思路,提出了三种不同边坡系数、断面结构形式和边坡筑构材料的防渗结构优化方案,并应用ABAQUS软件,分别计算分析了不同方案的衬砌结构受热力耦合作用下的冻胀变形特征与受力状态。结果表明,在湿陷性黄土地区,弧底梯形渠道的抗冻胀性能普遍优于梯形断面渠道;在渠道防渗衬砌体基土之间加设一定厚度的水泥砂浆垫层,能有效减弱衬砌结构在冻胀作用下的应力集中现象;适当减小渠道边坡系数能有效减少渠道内的冻胀应力峰值;在三种衬砌改造结构中,其峰值应力减幅分别为70.49%、63.36%、27.68%,最大冻胀量降低率分别为7.70%、5.05%、-2.82%,选取边坡系数为1∶1.75的弧底梯形断面型式,并在混凝土防渗面板与基土间设置水泥砂浆垫层时,其抗冻胀效果最佳,可在同类地区推广应用。该研究结果为跨越湿陷性黄土地基的灌区渠道结构改造提供参考,也可为同类地区渠道的防冻胀设计提供依据。

考虑土体非线性山-谷耦合复杂场地空间相关多点地震动模拟

采用考虑局部场地非线性地震波散射效应的空间相关多点地震动模拟方法,对山-谷耦合场地的地震动进行定量模拟分析。非线性地震波散射效应基于有限元-间接边界元耦合法和等效线性土体模型进行计算。结果表明,山体的存在会增大沉积内的地震动反应,对短周期的地震波的放大效应更为显著,增幅可达50%。相较于单一沉积谷地地形的沉积谷地外部,山-谷耦合场地中山体地震响应有所减小。相比线性模型,考虑土体非线性时,沉积谷地地震动强度较线性时的结果有所降低,反应谱峰值位置也会发生改变,但山体影响更为显著。不同山体形状对谷地影响不同,对于谷地内部分点位影响的差别可达40%。该研究为山-谷耦合场地工程结构抗震设防提供了部分理论依据和计算方法。