Borehole stability in naturally fractured rocks with drilling mud intrusion and associated fracture strength weakening:A coupled DFN-DEM approach

Borehole instability in naturally fractured rocks poses significant challenges to drilling.Drilling mud invades the surrounding formations through natural fractures under the difference between the wellbore pressure(P w)and pore pressure(P p)during drilling,which may cause wellbore instability.However,the weakening of fracture strength due to mud intrusion is not considered in most existing borehole stability analyses,which may yield significant errors and misleading predictions.In addition,only limited factors were analyzed,and the fracture distribution was oversimplified.In this paper,the impacts of mud intrusion and associated fracture strength weakening on borehole stability in fractured rocks under both isotropic and anisotropic stress states are investigated using a coupled DEM(distinct element method)and DFN(discrete fracture network)method.It provides estimates of the effect of fracture strength weakening,wellbore pressure,in situ stresses,and sealing efficiency on borehole stability.The results show that mud intrusion and weakening of fracture strength can damage the borehole.This is demonstrated by the large displacement around the borehole,shear displacement on natural fractures,and the generation of fracture at shear limit.Mud intrusion reduces the shear strength of the fracture surface and leads to shear failure,which explains that the increase in mud weight may worsen borehole stability during overbalanced drilling in fractured formations.A higher in situ stress anisotropy exerts a significant influence on the mechanism of shear failure distribution around the wellbore.Moreover,the effect of sealing natural fractures on maintaining borehole stability is verified in this study,and the increase in sealing efficiency reduces the radial invasion distance of drilling mud.This study provides a directly quantitative prediction method of borehole instability in naturally fractured formations,which can consider the discrete fracture network,mud intrusion,and associated weakening of fracture strength.The information provided by the numerical approach(e.g.displacement around the borehole,shear displacement on fracture,and fracture at shear limit)is helpful for managing wellbore stability and designing wellbore-strengthening operations.

C^1 natural element method for strain gradient linear elasticity and its application to microstructures

C^1 natural element method （C^1 NEM） is applied to strain gradient linear elasticity, and size effects on mi crostructures are analyzed. The shape functions in C^1 NEM are built upon the natural neighbor interpolation （NNI）, with interpolation realized to nodal function and nodal gradient values, so that the essential boundary conditions （EBCs） can be imposed directly in a Galerkin scheme for partial differential equations （PDEs）. In the present paper, C^1 NEM for strain gradient linear elasticity is constructed, and sev- eral typical examples which have analytical solutions are presented to illustrate the effectiveness of the constructed method. In its application to microstructures, the size effects of bending stiffness and stress concentration factor （SCF） are studied for microspeciem and microgripper, respectively. It is observed that the size effects become rather strong when the width of spring for microgripper, the radius of circular perforation and the long axis of elliptical perforation for microspeciem come close to the material characteristic length scales. For the U-shaped notch, the size effects decline obviously with increasing notch radius, and decline mildly with increasing length of notch.

Hybrid natural element method for viscoelasticity problems

A hybrid natural element method（HNEM） for two-dimensional viscoelasticity problems under the creep condition is proposed. The natural neighbor interpolation is used as the test function, and the discrete equation system of the HNEM for viscoelasticity problems is obtained using the Hellinger–Reissner variational principle. In contrast to the natural element method（NEM）, the HNEM can directly obtain the nodal stresses, which have higher precisions than those obtained using the moving least-square（MLS） approximation. Some numerical examples are given to demonstrate the validity and superiority of this HNEM.

Hybrid natural element method for large deformation elastoplasticity problems

We present the hybrid natural element method（HNEM） for two-dimensional elastoplastic large deformation problems. Sibson interpolation is adopted to construct the shape functions of nodal incremental displacements and incremental stresses. The incremental form of Hellinger–Reissner variational principle for elastoplastic large deformation problems is deduced to obtain the equation system. The total Lagrangian formulation is used to describe the discrete equation system.Compared with the natural element method（NEM）, the HNEM has higher computational precision and efficiency in solving elastoplastic large deformation problems. Some numerical examples are selected to demonstrate the advantage of the HNEM for large deformation elastoplasticity problems.

Numerical Study of Natural Convection in Square Cavity with Inner Bodies Using Finite Element Method

A numerical study of heat transfer problem by natural convection of a fluid inside a square cavity with two inner bodies is presented. This subject is of great interest in the engineering area, mainly in applications involving development of heat exchangers and cooling or heating systems of bodies by natural convection mechanism. Two cases have been studied. The inner bodies are square in case 1 and circular in case 2. In both cases, the bodies are solid and thermally conductive, the cavity lower and upper horizontal surfaces are isothermal with high temperature Th and low temperature Tc, respectively. Both vertical surfaces are adiabatic. A FORTRAN code using Finite Element Method (FEM) is developed to simulate the problem and solve the governing equations. The distributions of stream function, ψ, dimensionless temperature, θ, and vorticity, ω, are determined. Heat transfer is evaluated by analyzing the behavior of the average Nusselt number. The Grashof number and thermal diffusivity ratio are considered in range from 2 × 104 to 105 and from 0.1 to 100, respectively. The fluid is air with Prandtl number fixed in 0.733.

Complex fracture propagation model and plugging timing optimization for temporary plugging fracturing in naturally fractured shale

In this paper,a viscoelasticity-plastic damage constitutive equation for naturally fractured shale is deduced,coupling nonlinear tensile-shear mixed fracture mode.Dynamic perforation-erosion on fluid re-distribution among multi-clusters are considered as well.DFN-FEM(discrete fracture network combined with finite element method)was developed to simulate the multi-cluster complex fractures propagation within temporary plugging fracturing(TPF).Numerical results are matched with field injection and micro-seismic monitoring data.Based on geomechanical characteristics of Weiyuan deep shale gas reservoir in Sichuan Basin,SW China,a multi-cluster complex fractures propagation model is built for TPF.To study complex fractures propagation and the permeability-enhanced region evolution,intersecting and competition mechanisms between the fractures before and after TPF treatment are revealed.Simulation results show that:fracture from middle cluster is restricted by the fractures from side-clusters,and side-clusters plugging is benefit for multi fractures propagation in uniformity;optimized TPF timing should be delayed within a higher density or strike of natural fractures;Within a reservoir-featured natural fractures distribution,optimized TPF timing for most clustered method is 2/3 of total fluid injection time as the optimal plugging time under different clustering modes.

Simulation of Natural Convection Flow with Magneto-Hydrodynamics in a Wavy Top Enclosure with a Semi-Circular Heater

Natural convection flow in enclosure has different applications such as room ventilation, heat exchangers, the cooling system of a building etc. The Finite-Element method based on the Galerkin weighted residual approach is used to solve two-dimensional governing mass, momentum and energy-equations for natural convection flow in the presence of a magnetic field on a roof top with semi-circular heater. In the enclosure the horizontal lower wall was heated, the vertical two walls were adiabatic, inside the semi-circular heater, the wavy top wall cooled. The parameters Rayleigh number, Hartmann number and Prandtl number are considered. The effects of the Hartmann number and Rayleigh number on the streamlines, isotherms, velocity profiles and average Nusselt number are examined graphically. The local Nusselt number and the average Nusselt number of the heated portion of the enclosure with the semi-circular heater are presented in this paper. Finally, for the validation of the existing work, the current results are compared with published results and the auspicious agreement is achieved.

Finite Element Analysis of Magnetohydrodynamic Natural Convection within Semi-Circular Top Enclosure with Triangular Obstacles

The phenomena of magneto-hydrodynamic natural convection in a two-dimensional semicircular top enclosure with triangular obstacle in the rectangular cavity were studied numerically. The governing differential equations are solved by using the most important method which is finite element method (weighted-residual method). The top wall is placed at cold Tc and bottom wall is heated Th. Here the sidewalls of the cavity assumed adiabatic. Also all the wall are occupied to be no-slip condition. A heated triangular obstacle is located at the center of the cavity. The study accomplished for Prandtl number Pr = 0.71;the Rayleigh number Ra = 103, 105, 5 × 105, 106 and for Hartmann number Ha = 0, 20, 50, 100. The results represent the streamlines, isotherms, velocity and temperature fields as well as local Nusselt number.

Non-pneumatic mechanical elastic wheel natural dynamic characteristics and influencing factors

Non-pneumatic tire appears to have advantages over traditional pneumatic tire in terms of flat proof and maintenance free.A mechanical elastic wheel(MEW) with a non-pneumatic elastic outer ring which functions as air of pneumatic tire was presented.The structure of MEW was non-inflatable integrated configuration and the effect of hinges was accounted for only in tension. To establish finite element model of MEW, various nonlinear factors, such as geometrical nonlinearity, material nonlinearity and contact nonlinearity, were considered. Load characteristic test was conducted by tyre dynamic test-bed to obtain force-deflection curve. And the finite element model was validated through load characteristic test. Natural dynamic characteristics of the MEW and its influencing factors were investigated based on the finite element model. Simulation results show that the finite element model closely matched experimental wheel. The results also show that natural frequency is related to ground constraints, material properties, loads and torques. Influencing factors as above obviously affect the amplitude of mode of vibration, but have little effect on mode of vibration shape. The results can provide guidance for experiment research, structural optimization of MEW.

Investigation of Microstructure, Natural Frequencies and Vibration Modes of Dragonfly Wing

In the present work, a thorough investigation on the microstructural and morphological aspects of dragonfly wings was carried out using scanning electron microscope. Then, based on this study and the previous reports, a precise three-dimensional numerical model was developed and natural frequencies and vibration modes of dragonfly forewing were determined by finite element method. The results shown that dragonfly wings are made of a series of adaptive materials, which form a very complex composite structure. This bio-composite fabrication has some unique features and potential benefits. Furthermore, the numerical results show that the first natural frequency of dragonfly wings is about 168 Hz and bending is the predominant deformation mode in this stage. The accuracy of the present analysis is verified by comparison of calculated results with experimental data. This paper may be helpful for micro aerial vehicle design concerning dynamic response.

Optimization of Bearing Locations for Maximizing First Mode Natural Frequency of Motorized Spindle-Bearing Systems Using a Genetic Algorithm

This paper has developed a genetic algorithm (GA) optimization approach to search for the optimal locations to install bearings on the motorized spindle shaft to maximize its first-mode natural frequency (FMNF). First, a finite element method (FEM) dynamic model of the spindle-bearing system is formulated, and by solving the eigenvalue problem derived from the equations of motion, the natural frequencies of the spindle system can be acquired. Next, the mathematical model is built, which includes the objective function to maximize FMNF and the constraints to limit the locations of the bearings with respect to the geometrical boundaries of the segments they located and the spacings between adjacent bearings. Then, the Sequential Decoding Process (SDP) GA is designed to accommodate the dependent characteristics of the constraints in the mathematical model. To verify the proposed SDP-GA optimization approach, a four-bearing installation optimazation problem of an illustrative spindle system is investigated. The results show that the SDP-GA provides well convergence for the optimization searching process. By applying design of experiments and analysis of variance, the optimal values of GA parameters are determined under a certain number restriction in executing the eigenvalue calculation subroutine. A linear regression equation is derived also to estimate necessary calculation efforts with respect to the specific quality of the optimization solution. From the results of this illustrative example, we can conclude that the proposed SDP-GA optimization approach is effective and efficient.

Dynamic analysis of beam-cable coupled systems using Chebyshev spectral element method

The dynamic characteristics of a beam-cable coupled system are investigated using an improved Chebyshev spectral element method in order to observe the effects of adding cables on the beam. The system is modeled as a double Timoshenko beam system interconnected by discrete springs. Utilizing Chebyshev series expansion and meshing the system according to the locations of its connections, numerical results of the natural frequencies and mode shapes are obtained using only a few elements, and the results are validated by comparing them with the results of a finite-element method. Then the effects of the cable parameters and layout of connections on the natural frequencies and mode shapes of a fixed-pinned beam are studied. The results show that the modes of a beam-cable coupled system can be classified into two types, beam mode and cable mode, according to the dominant deformation. To avoid undesirable vibrations of the cable, its parameters should be controlled in a reasonable range, or the layout of the connections should be optimized.

Impacts of Heat Flux Distribution, Sloping Magnetic Field and Magnetic Nanoparticles on the Natural Convective Flow Contained in a Square Cavity

In the present paper,the effect of the heat flux distribution on the natural convective flow inside a square cavity in the presence of a sloping magnetic field and magnetic nanoparticles is explored numerically.The nondimensional governing equations are solved in the framework of a finite element method implemented using the Galerkin approach.The role played by numerous model parameters in influencing the emerging thermal and concentration fields is examined;among them are:the location of the heat source and its lengthH
,the magnitude of the thermal Rayleigh number,the nanoparticles shape and volume fraction,and the Hartmann number.It is found that the nanofluid velocity becomes higher when the thermal source length,the nanoparticles volume fraction and/or the thermal Rayleigh number are increased,while it decreases as the Hartmann number Ha grows and the position of the heat source moves toward the center of the lower wall of the cavity.Moreover,the temperature of the nanofluid grows with the extension of the thermal source and decreases slowly when the heat flux position moves toward the center of the lower wall.The outcomes of the research also indicate that the average Nusselt number becomes smaller on increasing Hartmann number Ha and heat source length H^(*).The addition of Fe_(3)O_(4) to engine oil leads to a higher rate of heat transfer with respect to the addition of SiO_(2) particles.Blade-shaped nanoparticles generate the highest value of the Nusselt number compared to all the other considered shapes.

Conjugate natural convection of non-Newtonian hybrid nanofluid in wavy-shaped enclosure

The present study concerns the modelization and numerical simulation for the heat and flow exchange characteristics in a novel configuration saturated with a nonNewtonian Ag-MgO hybrid nanofluid.The wavy shaped enclosure is equipped with onequarter of a conducting solid cylinder.The system of equations resulting from the mathematical modeling of the physical problem in its dimensionless form is discretized via the higher-order Galerkin-based finite element method(GFEM).The dependency of various factors and their interrelationships affecting the hydro-thermal behavior and heat exchange rate are delineated.The numerical experiments reveal that the best heat transfer rate is achieved for the pseudo-plastic hybrid nanoliquid with high Rayleigh number and thermal conductivity ratio and low Hartmann number.Besides,the power-law index has a major effect in deteriorating the heat convection at high Rayleigh number.

Natural frequency of tapered high pier considering pile-soil interaction

In seismic design of tapered high pier, the analysis of natural vibration frequency is of great importance. According to the engineering features of tapered high pier in mountainous area, a vibration calculation model was set up considering the tapered pier characteristics and pile-soil interaction. Based on Southwell frequency composition theory, it consists of elastic deformation of bridge pier and the rigid deformation of group piles, which are respectively solved by the finite-element method and energy method, and then the natural frequency is derived. The comparison between the measured and calculated results shows that the calculation errors with and without considering pile-soil interaction are 4.9% and 14.7%, respectively. Additionally, the main parameters (pier height, section variation coefficient and lateral foundation horizontal proportional coefficient) affecting natural frequency were investigated. The result shows that natural frequency ascends with the increase of the lateral foundation horizontal proportional coefficient; and it is quite necessary to consider the pile-soil interaction in natural frequency calculation of tapered high pier.

掺氢环境下含腐蚀缺陷X52天然气管道的结构强度失效机理研究

[目的]旨在探究掺氢环境下含腐蚀缺陷管道的结构强度劣化规律。[方法]基于试验数据分析X52钢在不同掺氢比例下的承载性能,采用非线性有限元方法研究X52掺氢天然气管道在不同腐蚀缺陷下的结构强度失效机理。[结果]研究结果表明:掺氢比例增大到50.0%时,管道在组合载荷下达到失效的位移最小为2.58 mm,结构韧性劣化程度较大;掺氢比例为25.0%时,缺陷深度越大,结构失效程度越严重,通过对正方形和圆形缺陷的对比,发现正方形缺陷的最大应力高于圆形缺陷。[结论]基于材料性能构建的含缺陷管道模型能够实现对不同掺氢环境下天然气管道服役状态结构承载力劣化程度的评估,所做研究可为新能源管道设计及健康监测提供工程化建议。

NEM与IEM耦合方法在地下工程中的应用研究

针对纯自然单元法(NEM)在处理地下工程无限域或半无限域问题时需要人为确定边界条件而带来计算误差的问题,引入无限元(IEM)模拟无穷远处边界条件,与自然元相结合形成耦合分析方法,并基于Laplace插值编制了相应程序;通过深埋圆形隧洞算例分析表明,与纯自然单元法相比,耦合方法显著提高了计算精度,选取较小的分析区域就可得到令人满意的结果,使得在满足工程要求前提下,采用耦合方法可以节省计算工作量,从而进一步增强了自然单元法处理地下工程问题的能力.

煤层气水平井扇形洞穴完井应力—渗透率演化机理

我国煤层渗透率低且地质条件复杂,采用常规油气储层改造的开发方式难度大、技术适应性差。近年来,基于应力释放的煤层气改造新方法“煤层气水平井水力喷射造穴”很好地解决了这一技术瓶颈问题,但是造穴卸压—增渗的作用机制及其主控地质因素尚不明晰。为此,考虑了煤岩层理和天然裂隙的影响,采用有限元—离散元耦合方法(Finite-Discrete Element Method,FDEM)建立了煤层气水平井扇形洞穴完井数值模型,探究了造穴后岩体的应力演化历程和储层的卸压—增渗机制,并对比分析了不同储层参数(孔隙压缩系数、储层强度、弱面强度和地应力场)对应力释放的影响规律。研究结果表明:(1)围岩演化过程为造穴后岩体收缩,储层发生应力重构,围岩强度逐渐降低,岩体内部发生新生裂隙萌生和原生裂隙扩展,形成开挖损伤区和应力释放区;(2)参数敏感性分析表明孔隙压缩系数是决定造穴完井储层适应性的关键,弱面强度、储层强度和地应力场分布决定了围岩的应力演化模式和裂缝扩展形态;(3)造穴卸压后储层增渗机制为穴周裂缝提升导流能力,储层应力释放提升基质渗透率。结论认为,模型首次综合考虑了地层特点、造穴过程和煤岩裂隙的影响,研究结果揭示了煤层造穴后的应力演化过程及其卸压、增渗作用机制,深化了对煤层气水平井洞穴完井增产机理的认识,对我国煤层储层改造具有重要的工程参考价值。

基于比例边界有限元的复合梁自由振动频率计算

将比例边界有限元方法(SBFEM)拓展用于计算复合梁的自由振动频率.该方法将梁简化为一维模型,并且仅选用x和z方向的弹性线位移作为基本未知量.从弹性力学基本方程出发,通过比例边界坐标、虚功原理和对偶变量技术推导得到了复合梁的一阶常微分比例边界有限元动力控制方程,其通解为解析的矩阵指数函数.利用Padé级数求解矩阵指数函数可得各个梁层的动力刚度矩阵,根据自由度匹配原则组装得到复合梁的整体刚度和质量矩阵.求解特征值方程,最终可得复合梁的自由振动频率.该方法对复合梁的层数和边界条件均无限制,具有广泛的适用性.将该文的解与三层、四层和十层复合梁振动频率的数值参考解以及阶梯型悬臂梁固有频率的实验实测值进行对比,验证了比例边界有限元算法的准确性、高效性和快速收敛性.

高墩大跨度连续刚构桥的动力力学特性分析思考

连续刚构桥是预应力混凝土连续梁桥中一种特殊结构型式,即是一种桥墩主墩与上部结构主梁固结的预应力混凝土连续梁桥。连续刚构桥多用于山区环境中,充分利用其跨径大、桥墩高的优点,完成对河谷、峡沟的跨越。其中,连续刚构桥的桥墩高度较大、桥墩柔性较强,且河谷、峡沟处风力往往较大,对连续刚构桥梁的动力影响较为突出。因此,连续刚构桥的动力研究分析,是连续刚构桥研究中的一个重难点。本论文为了研究高墩大跨连续刚构桥动力特性的影响,以国内某大桥为工程背景,利用有限元软件Midas Civil 2023建立主桥三维模型,从桥梁自震频率、地震反应谱、车桥耦合动力、桥墩高度、动力特性优化等方面分析该桥梁的动力特性。结果表明:○1主桥的自振频率与墩梁刚度比的变化呈正相关关系,变化速率呈先大后小的趋势;○2主桥的1阶纵向频率受主墩高度变化影响也很大,变化幅度大于50%。本论文研究结论可为高墩大跨连续刚构桥的动力分析设计提供参考和指导。