各向异性势问题充要的随机边界积分方程

本文针对各向异性势问题提出了一类充分必要的随机边界积分方程。数值计算结果表明在退化尺度附近，充要的随机边界积分方程较习用的随机边界积分方程有较大的优越性。     

势问题的重构核粒子边界无单元法

将重构核粒子法和势问题的边界积分方程方法结合，提出了势问题的重构核粒子边界无单元 法. 推导了势问题的重构核粒子边界无单元法的公式，研究其数值积分方案，建立了重构核 粒子边界无单元法的离散化边界积分方程，并推导了重构核粒子边界无单元法的内点位势的 积分公式. 重构核粒子法形成的形函数具有重构核函数的光滑性，且能再现多项式在插值点 的精确值，所以该方法具有更高的精度. 最后给出了数值算例，验证了所提方法的有效性 和正确性. }     

各向异性势问题充要的随机边界积分方程

本文针对各向异性势问题提出了一类充分必要的随机边界积分方程，数值计算结果表明在退化尺度附近，充要的随机边界积分程较习用的随机边界积分方程有较大的优越性。     

水平均流中细管排放气泡的三维数值模拟

在液体为无粘不可压，流动有势和气体遵循完全气体绝热关系的假定下，本文应用边界积分方程方法数值模拟了水平均流中垂直细管排放气泡的三维动力学问题，边界采用高阶有限元表达。文中介绍了有关泡面法向矢量、切向速度、曲率和接触线等的计算技术。与已知解的比较，表明了这一数值方法的高精度和优越性。算例显示了水平均流对于气泡形状和体积的影响     

《无奇异边界元法》内容简介

孙焕纯等著《无奇异边界元法》一书共有上下两篇 ,上篇阐述虚边界元法的理论、方法与应用。虚边界法有三种 :一般配点法 ;最小二乘配点法 (超额配点法 ) ;最小二乘二重积分法。*分别对弹性空间、弹性平面、薄板、薄壳问题给出了一个从弹性空间方程出发的统一的数值解法 ,抛弃了板、壳理论关于变形和应力的一切假设 ,又对位势问题、弹性平面问题等给出了边界积分方程离散化求解的系数阵元素的解析计算式。下篇针对传统边界元直接法与间接法的边界积分方程的充要性问题进行了论述 ,并对位势、弹性平面和薄板等问题建立了充要积分方程 ;其次是…     

三维常数势边界元中的精确积分

对三维问题边界元方法中应用最广泛的常数边界元的积分提出一种精确积分方法。借助于一个假想的闭合曲面,将特定的势场应用于边界积分方程,发现对于三维问题,常数势项的积分可以化作球面三角型的面积计算,而导数项的积分则可在平面域用极坐标进行。本文方法结果精确,公式简单,同一计算公式可以用来计算非奇异、几乎奇异和奇异积分,统一了积分算法。     

非连续边界元——有限元耦合方法分析

对边界元－有限元耦合方法进行了分析，采用非连续元离散边界积分方程，解决了耦合分析中自由度约束问题，给出了非连续边界元同有限元耦合的具体实施步骤，通过对二维弹性力学和流＝固耦合问题分析，表明了该文方法的有效性。     

缓变主流中三维气泡的非线性振动

空化现象和水下噪声机制与液体中气泡的动力学行为密切相关．在无粘势流的假定下，采用多参数摄动分析，研究了缓变主流中三维气泡的非线性体积模态振动．推导了关于缓变泡形展开的各阶扰动方程，获得了一阶振动的演化方程和一些特殊情况下的解析解；并采用高阶有限元离散的边界积分方程方法，对平面固壁和自由面附近三维气泡的固有频率进行了数值计算     

平面弹性力学充要的随机边界元

将平面弹性力学确定性的充分必要的边界积分方程推广到含材料常数随机的不确定问题中去，给出了位移的均值以及偏差的充分必要的边界积分方程。数值计算结果表明，和确定性的积分方程一样，习用的随机边界积分方程在退化尺度附近，无论是均值还是偏差都存在巨大的误差，而充要的随机边界积分方程则始终保持良好的精度     

弹性理论几类导数边界积分方程之间的变换关系

导数场边界积分方程通常难以应用，因为存在着超奇异主值积分的计算障碍。弹性理论中有几类不同的位移导数边界积分方程，本文采用算子δij和∈ij(排列张量)作用于这些导数边界积分方程，做一系列变换，原有的超奇异积分被正则化为强奇异积分获解。从而建立了这些位移导数边界积分方程之间的转换关系，它们均可以归结为自然边界积分方程。自然边界积分方程仅存在容易计算的Cauchy主值积分。自然边界积分方程分析可直接获得边界应力和位移导数。     

A TIME DOMAIN BOUNDARY ELEMENT METHOD FOR WATER-SOLID IMPACT ANALYSIS

A reciprocal theorem of dynamics for potential flow problems is first derived by meansof the Laplace transform in which the compressibility of water is taken into account.Based on this the-orem,the corresponding time-space boundary integral equation is obtained.Then,a set of time do-main boundary element equations with recurrence form is immediately formulated through discretiza-tion in both time and boundary.After having carried out the numerical calculation two solutions arefound in which a rigid semicircular cylinder and a rigid wedge with infinite length suffer normal impacton the surface of a half-space fluid.The results show that the present method is more efficient than theprevious ones.     

Time optimal motions of mechanical system with a prescribed trajectory

The problem of time minimization of a holonomic scleronomic mechanical system on a prescribed trajectory between two specified positions in configuration space is solved. The generalized force with restricted coordinates is taken as the controlling force. The application of the Green theorem (the well-known Miele method in flight mechanics) has shown that at every instant at least one control is at its boundary and possesses controlling functions with interruptions. It is assumed that at least one generalized coordinate exists that is monotonous during the interval of movement. An algorithm for numerical computation is presented for assessing the boundary of the admissible domain in the state space, thus, solving the problem of finding the optimal control as a function of time. Numerical integration is, therefore, carried out forward from the start point and backward from the end point by the use of the Runge-Kutta method. The mentioned procedure is illustrated in the example of time minimization for a manipulator which has its tip moving in a straight line. The application of the presented method simplifies solving of this type of problem compared to other methods, for instance, dynamic programming.     

Stokes flow before plane boundary with mixed stick-slip boundary conditions

A general theorem for the Stokes flow over a plane boundary with mixed stick-slip boundary conditions is established. This is done by using a representation for the velocity and pressure fields in the three-dimensional Stokes flow in terms of a biharmonic function and a harmonic function. The earlier theorem for the Stokes flow due to fundamental singularities before a no-slip plane boundary is shown to be a special case of the present theorem. Furthermore, in terms of the Stokes stream function, a corollary of the theorem is also derived, providing a solution to the problem of the axisymmetric Stokes flow along a rigid plane with stick-slip boundary conditions. The formulae for the drag and torque exerted by the fluid on the boundary are established. An illustrative example is given.     

A theorem for a shear stress-free plane boundary in Stokes flow

A theorem for non-axisymmetrical Stokes flow about a shear stress-free plane boundary is established by using a representation for the velocity and pressure fields for the same flow in terms of biharmonic and harmonic functions. A corollary of the theorem is derived which gives the axisymmetrical Stokes flow in terms of the Stokes function about the same boundary. The formulae for drag and torque on the boundary are also given. A few illustrative examples are presented.     

Uniform persistence and flows near a closed positively invariant set

In this paper, the behavior of a continuous flow in the vicinity of a closed positively invariant subset in a metric space is investigated. The main theorem in this part in some sense generalizes previous results concerning classification of the flow near a compact invariant set in a locally compact metric space which was described by Ura-Kimura (1960) and Bhatia (1969). By applying the obtained main theorem, we are able to prove two persistence theorems. In the first one, several equivalent statements are established, which unify and generalize earlier results based on Liapunov-like functions and those about the equivalence of weak uniform persistence and uniform persistence. The second theorem generalizes the classical uniform persistence theorems based on analysis of the flow on the boundary by relaxing point dissipativity and invariance of the boundary. Several examples are given which show that our theorems will apply to a wider varity of ecological models.     

Boundary conditions for axisymmetric circular cylinder in one-dimensional hexagonal quasicrystals

Through generalizing the method of a decay analysis technique determining the interior solution developed by Gregory and Wan, a set of necessary conditions on the end-data of axisymmetric circular cylinder in one-dimensional (1D) hexagonal quasicrystals (QCs) for the existence of a rapidly decaying solution is established. By accurate solutions for auxiliary regular state, and using the reciprocal theorem and general solution, these necessary conditions for the end-data to induce only a decaying elastostatic state (boundary layer solution) will be translated into appropriate boundary conditions for the circular cylinder with axisymmetric deformations in 1D hexagonal QCs. The results of the present paper enable us to establish a set of correct boundary conditions, most of which are obtained for the first time. Furthermore, corresponding necessary conditions for the case of axisymmetric circular cylinder with transversely isotropic are derived subsequently, and their isotropic elastic counterparts are also obtained. The accuracy of mixed boundary conditions of isotropic axisymmertic circular cylinder is proved at last.     

A Cauchy integral approach to Hele-Shaw problems with a free boundary: The case of zero surface tension

In this paper, we study a nonlinear and nonlocal free-boundary dynamics — the Hele-Shaw problem without surface tension when the fluid domain is either bounded or unbounded. The key idea is to use a global quantity, the Cauchy integral of the free boundary, to capture the motion of the boundary. This Cauchy integral is shown to be linear in time. The free boundary at a fixed time is then recovered from its Cauchy integral at that time. The main tool in our analysis isCherednichenko's theorem concerning the inverse properties of the Cauchy integrals.As products of our approach, we establish the short-time existence and uniqueness of classical solutions for analytic initial boundaries. We also show the non-existence of classical solutions for all smooth but non-analytic initial boundaries when there is a sink at either a finite point or at infinity. When the fluid domain is bounded, all solutions except the circular one break down before all the fluid is sucked out from the sink. Regularity results are also obtained when there is a source at a finite point or at infinity.     

ON THE MODIFIED CASTIGLIANO’S THEOREM

This paper gives the modified Castigliano’s theorem, which is more convenient and more extensive for applications than the classical Castigliano’s theorem.The modifications to the classical Castigliano ’s theorem are in following two respects, the first respect is that the expression of the partial derivative with respect to the concentrated load p of the complementary energy density in the clasical Castigliano ’s theorem is replaced by the expressions with teh products of the external loads and influence functions, the modification brings us greatest simplicity and greatest convenience for calculations under various loads; the second respect is that the expressions with the products of the inhomogeneous boundary displacements and the influence functions are introduced into the classical Castigliano’s theorem, the modification provides the theoretical fundamental for solving the problems of various boundary conditions.We show also the method of how to apply the modified Castigliano’s theorem to solve the problems of the surface structural mechanics.Finally, as a calculated example of the applications of the modified Castigliano’s theorem we solved the equation of the deflection surface of the rectangular plate with two adjacent edges built-in and other two adjacent edges free.     

Trefftz有限元法的研究进展

Trefftz有限元法(Trefftz finite element method,TFEM)是一种高效的数值计算方法,兼有传统有限元法和边界元法的诸多优点.基于双独立插值模式,结合杂交泛函和高斯散度定理,推得仅含边界积分的有限元格式.简述了过去10年间(2007—2016)Trefftz有限元法在单元域内插值函数、源项处理、特殊功能单元以及非各向同性材料等方面的研究进展,并对未来的发展趋势给出了几点展望.     

Velocity-based reciprocal theorems in elastodynamics and BIEM implementation issues

Reciprocal theorems in elastodynamics are introduced as extensions of respective theorems from elastostatics. Inasmuch as the latter is a subset of the former, the aim here is to present an elastodynamic reciprocal theorem that also includes elastostatics as a special case when the time variable becomes irrelevant. This is accomplished by introducing a velocity-based reciprocal theorem, whose basic properties are subsequently explored. The next step is to use this theorem and formulate a numerical approach based on boundary integral equation statements and compare them with existing formulations based on conventional reciprocity relations. The applications presented here involve the standard mechanical oscillator and the unidimensional axial element as two simple, yet important problems of structural dynamics. Along with the numerical results, a thorough stability analysis of the corresponding time-stepping algorithms is formulated. In both cases, the superior performance of the methodologies built on velocity-based reciprocal theorems is clearly demonstrated.