A new approach based on resonance technique and modified boundary ele-ment method is presented to calculate the impedance parameter matrix of a microwaveN-port network of waveguide structure.A two port network is take...A new approach based on resonance technique and modified boundary ele-ment method is presented to calculate the impedance parameter matrix of a microwaveN-port network of waveguide structure.A two port network is taken as a numerical ex-ample and the results show that the approach occupys the advantages of high accuracyand less computation effort.展开更多
In this work, a conceptual numerical solution of the two-dimensional wave partial differential equation (PDE) is developed by coupling the Complex Variable Boundary Element Method (CVBEM) and a generalized Fourier ser...In this work, a conceptual numerical solution of the two-dimensional wave partial differential equation (PDE) is developed by coupling the Complex Variable Boundary Element Method (CVBEM) and a generalized Fourier series. The technique described in this work is suitable for modeling initial-boundary value problems governed by the wave equation on a rectangular domain with Dirichlet boundary conditions and an initial condition that is equal on the boundary to the boundary conditions. The new numerical scheme is based on the standard approach of decomposing the global initial-boundary value problem into a steady-state component and a time-dependent component. The steady-state component is governed by the Laplace PDE and is modeled with the CVBEM. The time-dependent component is governed by the wave PDE and is modeled using a generalized Fourier series. The approximate global solution is the sum of the CVBEM and generalized Fourier series approximations. The boundary conditions of the steady-state component are specified as the boundary conditions from the global BVP. The boundary conditions of the time-dependent component are specified to be identically zero. The initial condition of the time-dependent component is calculated as the difference between the global initial condition and the CVBEM approximation of the steady-state solution. Additionally, the generalized Fourier series approximation of the time-dependent component is fitted so as to approximately satisfy the derivative of the initial condition. It is shown that the strong formulation of the wave PDE is satisfied by the superposed approximate solutions of the time-dependent and steady-state components.展开更多
The Complex Variable Boundary Element Method (CVBEM) procedure is extended to modeling applications of the two-dimensional linear diffusion partial differential equation (PDE) on a rectangular domain. The methodology ...The Complex Variable Boundary Element Method (CVBEM) procedure is extended to modeling applications of the two-dimensional linear diffusion partial differential equation (PDE) on a rectangular domain. The methodology in this work is suitable for modeling diffusion problems with Dirichlet boundary conditions and an initial condition that is equal on the boundary to the boundary conditions. The underpinning of the modeling approach is to decompose the global initial-boundary value problem into a steady-state component and a transient component. The steady-state component is governed by the Laplace PDE and is modeled using the Complex Variable Boundary Element Method. The transient component is governed by the linear diffusion PDE and is modeled by a linear combination of basis functions that are the products of a two-dimensional Fourier sine series and an exponential function. The global approximation function is the sum of the approximate solutions from the two components. The boundary conditions of the steady-state problem are specified to match the boundary conditions from the global problem so that the CVBEM approximation function satisfies the global boundary conditions. Consequently, the boundary conditions of the transient problem are specified to be continuously zero. The initial condition of the transient component is specified as the difference between the initial condition of the global initial-boundary value problem and the CVBEM approximation of the steady-state solution. Therefore, when the approximate solutions from the two components are summed, the resulting global approximation function approximately satisfies the global initial condition. In this work, it will be demonstrated that the coupled global approximation function satisfies the governing diffusion PDE. Lastly, a procedure for developing streamlines at arbitrary model time is discussed.展开更多
Mining operation, especially underground coal mining, always has the remarkable risks of ground control. Passive seismic velocity tomography based on simultaneous iterative reconstructive technique (SIRT) inversion is...Mining operation, especially underground coal mining, always has the remarkable risks of ground control. Passive seismic velocity tomography based on simultaneous iterative reconstructive technique (SIRT) inversion is used to deduce the stress redistribution around the longwall mining panel. The mining-induced microseismic events were recorded by mounting an array of receivers on the surface, above the active panel. After processing and filtering the seismic data, the three-dimensional tomography images of the p-wave velocity variations by SIRT passive seismic velocity tomography were provided. To display the velocity changes on coal seam level and subsequently to infer the stress redistribution, these three-dimensional tomograms into the coal seam level were sliced. In addition, the boundary element method (BEM) was used to simulate the stress redistribution. The results show that the inferred stresses from the passive seismic tomograms are conformed to numerical models and theoretical concept of the stress redistribution around the longwall panel. In velocity tomograms, the main zones of the stress redistribution around the panel, including front and side abutment pressures, and gob stress are obvious and also the movement of stress zones along the face advancement is evident. Moreover, the effect of the advance rate of the face on the stress redistribution is demonstrated in tomography images. The research result proves that the SIRT passive seismic velocity tomography has an ultimate potential for monitoring the changes of stress redistribution around the longwall mining panel continuously and subsequently to improve safety of mining operations.展开更多
A new algorithm combining nonlinear Galerkin method and coupling method of finite element and boundary element is introduced to solve the exterior nonstationary Navier-Stokes equations. The regularity of the coupling ...A new algorithm combining nonlinear Galerkin method and coupling method of finite element and boundary element is introduced to solve the exterior nonstationary Navier-Stokes equations. The regularity of the coupling variational formulation and the convergence of the approximate solution corresponding to the algorithm are proved. If the fine mesh h is chosen as coarse mesh H-square, the nonlinear Galerkin method, nonlinearity is only treated on the coarse grid and linearity is treated on the fine grid. Hence, the new algorithm can save a large amount of computational time.展开更多
The application of BEM to calculate acoustic radiation from closed surface in an infinite acoustic medium has the advantages of less memory space, higher accuracy and faster calculating speed. However, the singular in...The application of BEM to calculate acoustic radiation from closed surface in an infinite acoustic medium has the advantages of less memory space, higher accuracy and faster calculating speed. However, the singular integrals and the nonuniqueness problem of surface helmholtz equation at characteristic frequencies should be dealt with. It is suggested in this paper that combination of the interior Helmholtz equation and its derivatives with respect to the coordinate components of an interior point forms an extra equation, in order to solve the problems of acoustic radiation at various frequencies with surface Helmholtz equation; on the area which includes singular point, the singular integrals are converted to ordinary ones by polar coordinate transformation. As an example, the acoustic radiation field of closed axisymmetric surfaces with prescribed surface velocity distribution is calculated in this paper.展开更多
在采用离散元法分析机械部件与颗粒材料接触作用时,需要建立机械部件(边界)的离散元法分析模型。分析可知,机械部件中与颗粒材料接触作用的零件表面,存在不能用初等解析函数表达的非规则曲面。为此,采用推进波前法(AFT:Advancing Front ...在采用离散元法分析机械部件与颗粒材料接触作用时,需要建立机械部件(边界)的离散元法分析模型。分析可知,机械部件中与颗粒材料接触作用的零件表面,存在不能用初等解析函数表达的非规则曲面。为此,采用推进波前法(AFT:Advancing Front Technique)进行非规则曲面网格划分,把非规则曲面离散成三角形平面片的组合,同时添加运动属性和材料特性参数,由此建立非规则曲面边界的离散元法分析模型。在对PRO/E软件进行二次开发的基础上,研制了非规则曲面边界建模软件。通过实例验证,初步证明了基于AFT边界建模方法和软件的可行性,为复杂结构机械部件工作过程的仿真分析奠定了基础。展开更多
在势流假设下,利用摄动原理,求解 Stokes 波中铅垂圆柱的二阶绕射势,得到完整的二阶波浪压力和波浪力。根据铅垂圆柱的几何特性,采用两种求解方法。一是直接利用外场解析表达式求解;二是设置控制面,外场用解析表达式,内场用简单 Green ...在势流假设下,利用摄动原理,求解 Stokes 波中铅垂圆柱的二阶绕射势,得到完整的二阶波浪压力和波浪力。根据铅垂圆柱的几何特性,采用两种求解方法。一是直接利用外场解析表达式求解;二是设置控制面,外场用解析表达式,内场用简单 Green 函数方法求解。两种方法得到的结果吻合良好,验证了计算方法的有效性,可为求解任意三维物体的二阶绕射问题提供参考。展开更多
文摘A new approach based on resonance technique and modified boundary ele-ment method is presented to calculate the impedance parameter matrix of a microwaveN-port network of waveguide structure.A two port network is taken as a numerical ex-ample and the results show that the approach occupys the advantages of high accuracyand less computation effort.
文摘In this work, a conceptual numerical solution of the two-dimensional wave partial differential equation (PDE) is developed by coupling the Complex Variable Boundary Element Method (CVBEM) and a generalized Fourier series. The technique described in this work is suitable for modeling initial-boundary value problems governed by the wave equation on a rectangular domain with Dirichlet boundary conditions and an initial condition that is equal on the boundary to the boundary conditions. The new numerical scheme is based on the standard approach of decomposing the global initial-boundary value problem into a steady-state component and a time-dependent component. The steady-state component is governed by the Laplace PDE and is modeled with the CVBEM. The time-dependent component is governed by the wave PDE and is modeled using a generalized Fourier series. The approximate global solution is the sum of the CVBEM and generalized Fourier series approximations. The boundary conditions of the steady-state component are specified as the boundary conditions from the global BVP. The boundary conditions of the time-dependent component are specified to be identically zero. The initial condition of the time-dependent component is calculated as the difference between the global initial condition and the CVBEM approximation of the steady-state solution. Additionally, the generalized Fourier series approximation of the time-dependent component is fitted so as to approximately satisfy the derivative of the initial condition. It is shown that the strong formulation of the wave PDE is satisfied by the superposed approximate solutions of the time-dependent and steady-state components.
文摘The Complex Variable Boundary Element Method (CVBEM) procedure is extended to modeling applications of the two-dimensional linear diffusion partial differential equation (PDE) on a rectangular domain. The methodology in this work is suitable for modeling diffusion problems with Dirichlet boundary conditions and an initial condition that is equal on the boundary to the boundary conditions. The underpinning of the modeling approach is to decompose the global initial-boundary value problem into a steady-state component and a transient component. The steady-state component is governed by the Laplace PDE and is modeled using the Complex Variable Boundary Element Method. The transient component is governed by the linear diffusion PDE and is modeled by a linear combination of basis functions that are the products of a two-dimensional Fourier sine series and an exponential function. The global approximation function is the sum of the approximate solutions from the two components. The boundary conditions of the steady-state problem are specified to match the boundary conditions from the global problem so that the CVBEM approximation function satisfies the global boundary conditions. Consequently, the boundary conditions of the transient problem are specified to be continuously zero. The initial condition of the transient component is specified as the difference between the initial condition of the global initial-boundary value problem and the CVBEM approximation of the steady-state solution. Therefore, when the approximate solutions from the two components are summed, the resulting global approximation function approximately satisfies the global initial condition. In this work, it will be demonstrated that the coupled global approximation function satisfies the governing diffusion PDE. Lastly, a procedure for developing streamlines at arbitrary model time is discussed.
文摘Mining operation, especially underground coal mining, always has the remarkable risks of ground control. Passive seismic velocity tomography based on simultaneous iterative reconstructive technique (SIRT) inversion is used to deduce the stress redistribution around the longwall mining panel. The mining-induced microseismic events were recorded by mounting an array of receivers on the surface, above the active panel. After processing and filtering the seismic data, the three-dimensional tomography images of the p-wave velocity variations by SIRT passive seismic velocity tomography were provided. To display the velocity changes on coal seam level and subsequently to infer the stress redistribution, these three-dimensional tomograms into the coal seam level were sliced. In addition, the boundary element method (BEM) was used to simulate the stress redistribution. The results show that the inferred stresses from the passive seismic tomograms are conformed to numerical models and theoretical concept of the stress redistribution around the longwall panel. In velocity tomograms, the main zones of the stress redistribution around the panel, including front and side abutment pressures, and gob stress are obvious and also the movement of stress zones along the face advancement is evident. Moreover, the effect of the advance rate of the face on the stress redistribution is demonstrated in tomography images. The research result proves that the SIRT passive seismic velocity tomography has an ultimate potential for monitoring the changes of stress redistribution around the longwall mining panel continuously and subsequently to improve safety of mining operations.
文摘A new algorithm combining nonlinear Galerkin method and coupling method of finite element and boundary element is introduced to solve the exterior nonstationary Navier-Stokes equations. The regularity of the coupling variational formulation and the convergence of the approximate solution corresponding to the algorithm are proved. If the fine mesh h is chosen as coarse mesh H-square, the nonlinear Galerkin method, nonlinearity is only treated on the coarse grid and linearity is treated on the fine grid. Hence, the new algorithm can save a large amount of computational time.
基金The project supported by National Natural Science Foundation of China
文摘The application of BEM to calculate acoustic radiation from closed surface in an infinite acoustic medium has the advantages of less memory space, higher accuracy and faster calculating speed. However, the singular integrals and the nonuniqueness problem of surface helmholtz equation at characteristic frequencies should be dealt with. It is suggested in this paper that combination of the interior Helmholtz equation and its derivatives with respect to the coordinate components of an interior point forms an extra equation, in order to solve the problems of acoustic radiation at various frequencies with surface Helmholtz equation; on the area which includes singular point, the singular integrals are converted to ordinary ones by polar coordinate transformation. As an example, the acoustic radiation field of closed axisymmetric surfaces with prescribed surface velocity distribution is calculated in this paper.
文摘在采用离散元法分析机械部件与颗粒材料接触作用时,需要建立机械部件(边界)的离散元法分析模型。分析可知,机械部件中与颗粒材料接触作用的零件表面,存在不能用初等解析函数表达的非规则曲面。为此,采用推进波前法(AFT:Advancing Front Technique)进行非规则曲面网格划分,把非规则曲面离散成三角形平面片的组合,同时添加运动属性和材料特性参数,由此建立非规则曲面边界的离散元法分析模型。在对PRO/E软件进行二次开发的基础上,研制了非规则曲面边界建模软件。通过实例验证,初步证明了基于AFT边界建模方法和软件的可行性,为复杂结构机械部件工作过程的仿真分析奠定了基础。
文摘在势流假设下,利用摄动原理,求解 Stokes 波中铅垂圆柱的二阶绕射势,得到完整的二阶波浪压力和波浪力。根据铅垂圆柱的几何特性,采用两种求解方法。一是直接利用外场解析表达式求解;二是设置控制面,外场用解析表达式,内场用简单 Green 函数方法求解。两种方法得到的结果吻合良好,验证了计算方法的有效性,可为求解任意三维物体的二阶绕射问题提供参考。