单侧摩擦约束夹杂物的瞬态弹性波散射的时域边界元分析--反平面情况

应用时域边界元法研究了瞬态SH波对单侧摩擦约束夹杂物的散射问题，假设摩擦遵守库仑定量，当入射SH波足够强时界面会出现局部滑移，发现了一种有效的迭代技术以确定边界上的未知区域（滑移和粘着区）。算例计算了无限域中圆柱埋置夹杂物对瞬态SH波的散射问题。     

一般各向异性单侧接触界面上波的反射和折射

研究简谐弹性波在一般各向异性介质单侧接触界面上的反射和折射问题．利用Fouier分析方法将非线性Coulomb摩擦接触边界波动问题化为一组代数方程．给出了确定局部分离、滑移和粘着区的思路和方法及各区域的解；讨论了出现界面局部分离和滑移的条件．对特定材料组合情况进行了详细数值计算，给出了界面力、相对滑移速度、张开位移、高频谐波的反射折射系数等特征参量；考察了平面和反平面波动的耦合及整体滑移等．其中关于高频谐波的结果可对已有实验结果给出很好的定性解释．在大多数情况下，即使对摩擦系数无穷大的粘滞接触界面，分离区端部也总是存在一个很小的滑移区。     

反平面弹性波场中的单侧摩擦界面

在Coulomb摩擦假定下，详细地讨论了反平面波场中单侧摩擦界面问题；用函数方程法求解了这类问题，给出了界面状态的严格解析解。     

垂直入射弹性波在压电体单侧接触界面上的传播特性

弹性波与压电材料接触界面的相互作用问题是工程应用中常见而复杂的问题,入射波足够强会引起界面出现滑移和分离,但滑移和分离的边界未知,边界条件具有非线性特性。通过Fourier分析,将混合边值问题的求解转化为非线性代数方程,利用软件通过迭代修正的方法进行了求解;给出3种状态边界的求解,分析入射波强度、外加应力及电场对界面状态的影响,并对高频谐波的特性进行分析,通过实例对理论推导进行验证,结果显示:入射波强度、外加荷载和电场的大小及摩擦因数均会影响到界面,通过改变这些条件可以控制界面状态,另外检测高频谐波的信号也可以反映界面状态。     

双材料单伸臂叠合梁的界面滑移与自振特性分析

将叠合梁划分为接触区和分离区,接触区界面间的摩擦作用会对叠合梁的滑移、刚度和自振频率产生影响。本文给出了单伸臂叠合梁在均布荷载和集中力作用下,考虑叠合界面摩擦作用的滑移应变和滑移分布的表达式;推导了考虑叠合界面间摩擦力及摩擦力产生的抵抗弯矩共同作用下的截面刚度。假设梁按等波长和等刚度两种形式自由振动,运用传递矩阵法推导出了单伸臂叠合梁自振频率的计算公式。最后,基于一个钢-混凝土组合与叠合试验梁的算例分析了滑移沿梁长的分布及自振频率。结果表明:随界面间摩擦系数的增大,接触区界面的滑移值逐渐减小;摩擦效应对叠合梁自振频率的影响随阶次的升高越来越明显。     

基于边界元法求解三维弹性摩擦接触问题

边界元法求解三维摩擦接触问题,其中一个关键点在于如何确定滑移方向。即当出现相对滑移时,滑移方向如何确定。当前常采用的方法是,粘结点利用切向面力得到滑移方向,滑移点利用切向相对位移得到滑移方向。不过该方法难以保证收敛性。针对这一问题,本文采用滑移方向预测技术得到滑移方向。即以后出现相对滑移时,滑移方向采用预测技术中得到的滑移方向。由于摩擦接触问题和历史加载相关,本文采用增量法求解。不同摩擦系数下的数值结果都证明了本文算法的有效性和收敛性及滑移方向预测技术的有效性。     

垂直入射弹性波在一般各向异性压电材料单侧接触界面上的传播

弹性波垂直于压电材料摩擦接触界面入射会引起界面局部滑移或分离,导致问题具有边界非线性。应用傅里叶分析技术将问题的求解转化为一组代数方程,通过迭代-修正的方法确定粘着、滑移、分离区的分布位置,给出了各种状态产生的条件及规律。由于边界非线性导致波形畸变并诱发高频谐波,对某些各向异性压电材料,在某些条件下单一入射波遇到界面产生再极化,导致耦合P波及SV波产生,文中给出了高频谐波解及波形极化的产生条件。外加荷载和外加电场通过机电耦合效应对界面产生影响,通过实例分析,给出了外加作用力和电场对一般各向异性材料的影响规律,并分析了因为材料常数不同而引起的影响规律的不同。     

近场波动模拟的一种应力人工边界

采用平面波和远场散射波混合透射，引入无限介质线弹性本构关系建立了一种应力人工 边界条件. 其优点在于边界结点反应与内部有限元结点反应采用相同的积分格式计算，有限 元积分方法稳定时不存在人工边界失稳问题. 数值算例表明：边界精度高于现有黏性边 界、黏弹性人工边界，以及一、二阶透射人工边界.     

弹性波作用下压电体摩擦接触界面滑移特性分析

弹性波与压电体摩擦接触界面相互作用会引起界面滑移或分离,滑移和分离位置的分布与外加压力、剪力、电场及入射波的条件有关。应用Fourier分析及matlab软件给出了滑移范围的解,给出了粘着、滑移或分离的判定条件,通过算例分析了滑移和分离出现的影响因素和条件,为工程应用提供了理论依据。     

板材冲压成形变形局部化与动边界摩擦约束

提出适于具有强烈弯曲效应的大应变，大位移与界面摩擦三重非线性空间曲壳有限元列式。给出一种动边界摩擦约束的元级代数消元方法，就板材在方盒型冲头作用中压胀形过程，引入几种现代非经典塑性本构理论以及空单元技术，研究了成形过程中的变形局部化与断裂问题。     

DYNAMIC INTERACTION OF PLANE WAVES WITH A UNILATERALLY FRICTIONALLY CONSTRAINED INCLUSION-TIME DOMAIN BOUNDARY ELEMENT ANALYSIS

A 2D time domain boundary element method (BEM) is developed to solve the transient scattering of plane waves by a unilaterally frictionally constrained inclusion. Coulomb friction is assumed along the contact interface. The incident wave is assumed strong enough so that localized slip and separation take place along the interface. The present problem is in effect a nonlinear boundary value problem since the mixed boundary conditions involve unknown intervals (slip, separation and stick regions). In order to determine the unknown intervals, an iterative technique is developed. As an example, we consider the scattering of a circular cylinder embeddedin an infinite solid.     

流体流动的边界滑移问题研究进展

最近十几年来,随着现代微/纳米测试以及分子动力学模拟技术的出现和发展,流体流动的边界滑移问题研究获得了突 飞猛进的发展.边界滑移相关研究大体可分为3个方面: 实验、分子动力学模拟和理论数值分析,前两者主要以发现边界滑移现象、探索边界滑移的产生机理以及各因素对边 界滑移的影响规律为主要研究目的,而后者主要研究边界滑移的物理模型、相关问题的计算方法以及边界滑移对流体系 统流体动力学行为的影响.本文首先简要回顾了液体流动的边界滑移及其相关问题的早期研究历史,随后对边界滑移问题 的研究现状进行了综述,最后展望了该领域今后的研究重点及其应用前景.     

Inverse scattering problem from an impedance crack via a composite method

In this paper we consider an inverse scattering problem from an open arc with impedance boundary conditions on both sides of the crack. Our aim is to recover both the impedance function and the unknown crack simultaneously from the far-field pattern with only one incident wave. Making the most out of the direct problem, a straightforward method of iterative nature is developed for the inverse problem. The ill-posedness of this problem is considered by incorporating the Tikhonov regularization. Numerical examples are provided at the end of the paper to show the feasibility of our method.     

Analysis of transient wave scattering from an inclusion with an unilateral smoothly contact interface by BEM

Summary A 2D time-domain Boundary Element Method (BEM) is applied to solve the problem of transient scattering of plane waves by an inclusion with a unilateral smooth contact interface. The incident wave is assumed strong enough so that localized separations take place along the interface. The present problem is indeed a nonlinear boundary value problem since the mixed boundary conditions involve unknown intervals (separation and contact regions). In order to determine the unknown intervals, an iterative technique is developed. As an example, we consider the scattering of plane waves by the cross section of a circular cylinder embedded in an infinite solid. Numerical results for the near field solutions are presented. The distortion of the response waves and the variation of the interface states are discussed. The financial support by the China National Natural Science Foundation under Grant No. 19872001 and No. 59878004 is gratefully acknowledged. The second author is also grateful to the support of the National Science Fund for Distinguished Young Scholars under Grant No. 10025211.     

On the application of slip boundary condition on curved boundaries

Hydrodynamic simulations of sloshing phenomena often involve the application of slip boundary condition at the wetted surfaces. If these surfaces are curved, the ambiguous nature of the normal vector in the discretized problem can interfere with the application of such a boundary condition. Even the use of consistent normal vectors, preferred from the point of view of conservation, does not assure good approximation of the continuum slip condition in the discrete problem, and non‐physical recirculating flow fields may be observed. As a remedy, we consider the Navier slip condition, and more successfully, the so‐called BC‐free boundary condition. Copyright © 2004 John Wiley & Sons, Ltd.     

Geometrically necessary dislocations in viscoplastic single crystals and bicrystals undergoing small deformations

In this study we develop a gradient theory of small-deformation single-crystal plasticity that accounts for geometrically necessary dislocations (GNDs). The resulting framework is used to discuss grain boundaries. The grains are allowed to slip along the interface, but growth phenomenona and phase transitions are neglected. The bulk theory is based on the introduction of a microforce balance for each slip system and includes a defect energy depending on a suitable measure of GNDs. The microforce balances are shown to be equivalent to nonlocal yield conditions for the individual slip systems, yield conditions that feature backstresses resulting from energy stored in dislocations. When applied to a grain boundary the theory leads to concomitant yield conditions: relative slip of the grains is activated when the shear stress reaches a suitable threshold; plastic slip in bulk at the grain boundary is activated only when the local density of GNDs reaches an assigned threshold. Consequently, in the initial stages of plastic deformation the grain boundary acts as a barrier to plastic slip, while in later stages the interface acts as a source or sink for dislocations. We obtain an exact solution for a simple problem in plane strain involving a semi-infinite compressed specimen that abuts a rigid material. We view this problem as an approximation to a situation involving a grain boundary between a grain with slip systems aligned for easy flow and a grain whose slip system alignment severely inhibits flow. The solution exhibits large slip gradients within a thin layer at the grain boundary.     

Meshless Galerkin analysis of Stokes slip flow with boundary integral equations

This paper presents a novel meshless Galerkin scheme for modeling incompressible slip Stokes flows in 2D. The boundary value problem is reformulated as boundary integral equations of the first kind which is then converted into an equivalent variational problem with constraint. We introduce a Lagrangian multiplier to incorporate the constraint and apply the moving least‐squares approximations to generate trial and test functions. In this boundary‐type meshless method, boundary conditions can be implemented exactly and system matrices are symmetric. Unlike the domain‐type method, this Galerkin scheme requires only a nodal structure on the bounding surface of a body for approximation of boundary unknowns. The convergence and abstract error estimates of this new approach are given. Numerical examples are also presented to show the efficiency of the method. Copyright © 2009 John Wiley & Sons, Ltd.     

The boundary element method applied to elastic contact problems with friction

In the current paper the boundary integral equations (BIE) for elastic contact problems with friction are derived from the incremental virtual work principle. After introducing contact conditions of adhesion and slip into BIE all variants on boundary are made to discretize by quadratic isoparametric boundary element. In the current paper not only an auto-increment loading law is presented but also the iterative calculation laws for open, slip and adhesion condition are given. The results of numerical examples are satisfactory.     

Necessary and sufficient condition for the correct formulation of boundary integral equations of harmonic functions

A necessary and sufficient condition for the correct formulation of boundary integral equations of harmonic functions is established in the present paper. A super-determined problem of harmonic functions is proposed for the first time. Then a necessary and sufficient condition for the existence of solution of the super-determined problem is proved. At the same time, it is a necessary and sufficient condition for the correct formulation of boundary integral equations with direct unknown quantities. A relation between boundary integral equations and variational principles is discovered for the first time. And a one-to-one correspondence between boundary integral equations with direct and indirect unknown quantities is indicated. The concept and route of the present paper can be applied to other boundary value problems possessing variational principles.