华丰矿动压低强度介质围岩巷道稳定性模拟

采用FEM与BEM耦合技术对华丰矿受采动影响的底极及两帮为煤层的锚网支护顺槽巷道围岩稳定性进行数值模拟。结果表明，巷道的两帮及底板破坏范围大，变形严重，特别是底板的鼓起变形，是巷道围岩控制的重点。数值模拟结果与实际情况基本吻合，为采取行之有效的支护方式对底板进行加固提供了依据。     

BEM analysis of elastic foundation beams on multilayered isotropic soils

This paper presents a boundary element method to analyze the elastic foundation finite beams on 2D plane-strain and 3D multilayered isotropic soils. Starting with the basic equations of elasticity and based on the Fourier transform, the transfer matrix of a single soil layer is derived. According to the boundary conditions and continuity conditions between two adjacent soil layers, the solution of multilayered elastic soils is obtained to be a kernel function of BEM analysis. The elastic foundation beam is modeled as a Bernoulli–Euler beam using the finite difference method. With the displacement and stress condition of coordination between beam and soil, the solution is acquired for beams resting on multilayered soils. Comparing the solution with the published data shows that the solution is in good agreement, and some numerical examples show that the beam behavior is affected vitally by soil–beam stiffness ratio and the stratification of soils.     

Boundary element analysis of rotationally symmetric systems under general boundary conditions

Engineering systems often exhibit cyclic symmetry (or rotational periodicity) of geometrical and physical properties. This paper deals with the use of this feature in the context of the direct boundary element method, and shows how to decompose an elastostatic and steady state b.v. problem with general (not cyclically symmetric) boundary conditions into a number of much smaller problems over a single segment.

There is no need to discretize the interfaces between adjacent segments. The decomposition is achieved by rotational transformations of physical into cyclic components     

补偿阳极位置对舰船腐蚀电场防护效果的影响

为了减少舰船腐蚀电场对船上电磁通信设备的干扰,需要对舰船腐蚀电场进行抑制或消除. 采用三维边界元法建立潜艇腐蚀静电场模型,研究补偿阳极位置对腐蚀电场防护效果的影响,结合电偶极子模型理论解释了静电场防护原理,基于静电场防护原理,分析不同补偿阳极位置产生舰船腐蚀电场差异的原因. 结果表明:补偿阳极的周向位置对电场防护效果影响较小,而纵向位置对电场防护效果有明显的影响;补偿阳极与螺旋桨的纵向距离分别为16.6、10.8、5、1 m时,电场模量的防护效果分别为3.41%、24.88%、49.27%和34.15%.     

高反照率屋顶对城市热岛及空调能耗的影响

为了研究高反照率屋顶这一城市降温技术对夏季城市热岛效应的缓解作用及对空调制冷负荷的影响,以夏季高温频发的杭州市为例,选取2017年7月22日—28日作为典型高温时段,采用中尺度天气预报模式（WRF）耦合考虑建筑物能量交换的多层城市冠层参数化方案（BEP+BEM）进行数值模拟. 结果表明,反照率为0.85的屋顶使城市区域日均降温0.37 °C,城市热岛强度减小0.21 °C,城市热岛效应得到了一定程度的缓解,空调制冷负荷降低约5.3%;建筑物密集的商业区的降温和节能效果均优于低密度居住区;屋顶反照率与城区2 m气温及空调能耗均线性负相关.     

开口薄壁结构声固耦合问题的计算方法

对于开口薄壁结构的声辐射问题,很难用解析法直接计算,可采用多种数值方法进行仿真分析。文中研究了有限元法、边界元法与无限元法几种数值方法及其耦合的理论基础,研究了求解不封闭声腔声辐射问题的几种数值方法。在考虑声固耦合的前提下,分析了间接边界元法、直接边界元内外声场耦合法、有限元与边界元耦合法、有限元与无限元耦合法4种方法求解开口结构声场问题的系统方程及求解过程,以具体开口薄壁结构为实例采用4种方法进行了模拟计算,并对结果进行了分析比较,得出相关结论。     

Multiple pole residue approach for 3D BEM analysis of mathematical degenerate and non‐degenerate materials

In this paper we develop an alternative boundary element method (BEM) formulation for the analysis of anisotropic three‐dimensional (3D) elastic solids. Our implementation is based on the derivation of explicit expressions for the fundamental solution displacements and tractions, of general validity for any class of anisotropic materials, by means of Stroh formalism and Cauchy's residue theory. The resulting fundamental solution remains valid for mathematical degenerate cases when Stroh's eigenvalues are coincident, meanwhile it does not exhibit numerical instabilities for quasi‐degenerate cases when Stroh's eigenvalues are nearly equal. A multiple pole residue approach is followed, leading to general explicit expressions to evaluate the traction fundamental solution for poles of m‐multiplicity. Despite the existence of general displacement solutions in the literature, and for the sake of completeness, the same approach as for the traction solution is considered to derive the displacement fundamental solution as well. Based on these solutions, an explicit BEM approach for the numerical solution of 3D linear elastic problems for solids with general anisotropic behavior is presented. The analysis of cracked anisotropic solids is also considered. Details on the numerical implementation and its validation for degenerate cases are discussed. Copyright © 2010 John Wiley & Sons, Ltd.     

各向异性外问题的非均匀网格的自然边界元法及其耦合法

本文对于无界区域上各向异性外问题提出了在椭圆边界非均匀网格上的自然边界元法及其与有限元法的耦合法,证明相应的收敛定理和误差估计式,并且在这两种方法中引入基于等分布原理的移动网格技巧.最后,通过数值结果表明了误差收敛理论的正确性以及所提方法和技巧的有效性.     

On subsurface box-shaped lined tunnel under incident SH-wave propagation

In this paper, a half-plane time-domain boundary element method is applied to obtain the seismic ground response, including a subsurface box-shaped lined tunnel deployed in a linear homogenous elastic medium exposed to obliquely incident SH-waves. Only the boundary around the tunnel is required to be discretized. To prepare an appropriate model by quadratic elements, a double-node procedure is used to receive dual boundary fields at corners as well as change the direction of the normal vector. After encoding the method in a previously confirmed computer program, a numerical study is carried out to sensitize some effective parameters, including frequency content and incident wave angle for obtaining a surface response. The depth and impedance ratio of the lining are assumed to be unvaried. The responses are illustrated in the time and frequency domains as two/three-dimensional graphs. The results showed that subsurface openings with sharp corners distorted the propagation path of the anti-plane waves to achieve the critical states on the ground surface. The present approach can be proposed to civil engineers for preparing simple underground box-shaped models with angular boundaries.     

Dynamic crack propagation in matrix involving inclusions by a time-domain BEM

The dynamic crack propagation in composites with inclusions is analyzed by a time-domain boundary element method (BEM) in conjunction with the sub-region technique. The crack-growth direction and the crack-tip instantaneous velocity are determined by the maximum circumferential stress criterion. The instantaneous velocity is well smoothed by a bisection technique. New crack-tip elements of inconstant length are added to the active crack-tip to simulate the fast crack growth. Square root shape functions are adopted as to describe the proper asymptotic behavior in the vicinity of the crack-tips. The computation time for the dynamic fracture problems in composites with multiple inclusions is reduced by a numerical method. The influences of the inclusions on the shielding ratio, the crack growing path and the crack-tip instantaneous speed are well investigated.