断裂力学的半解析相似边界元法

本文首先对弹性力学的相似边界元法进行了研究，推导了相应的计算公式。与传统的边界元法相比，相似边界元法由于只需在少数单元上进行数值积分，大大减少了计算量。在此基础上，对断裂力学问题，利用裂纹尖端位移场的解析表达式将裂纹尖端节点未知量转化为几个待定常数，提出了半解析相似边界元法，可大大减少最终形成的线性代数方程组的系数矩阵的阶数，进一步减小计算量。最后给出了算例，说明了本文方法的有效性。     

关于边界元方法中柯西主值积分的探讨

本文对边界元方法中的各类积分根据其奇异性作分类，并对主值积分的收敛条件、变量替换等进行了讨论，又给出了变替换附加项显式。文中提供的主值积分配项消奇术在边界元方法中是有普遍意义的。     

扁球壳的边界元几何非线性分析

本文从壳体位移的三个微分方程出发，采用付立叶积分变换的基本解，利用加权残值法推导了几何非线性边界积分方程。这种基本解的壳体边界元法类似于板的非线性边界元法，各种变量物理意义明确，能方便地处理各种复杂边界条件及有开口情况。文末算例说明本文方法的可行性、收敛性和精确性，并与二变量边界单元法或有限元结果相比较，吻合较好。     

中厚板统计分析的随机边界元法

本文建立了分析含随机材料参数并具厚度不均匀性的中厚板问题的随机边界元法,基于Taylor级数展开技术,分析和到广义位移的均值和一阶偏差的积分方程,其中将材料参数的随机性和厚度的不均匀性作为等效荷载处理,从而得到广义边界位移或面力的均值和协方差,并进一步求出部点广义位移和内力的均值和协方差,最后用本文方法计算了两个数例,并对所得结果进行了分析,探讨。     

弹性薄板弯曲问题的边界轮廓法

导出了弹性薄板弯曲问题边界积分方程的另一种形式,基于这种方程,提出了平板弯曲问题的边界轮廓法,讨论了三次边界单元边界轮廓法的计算列式,并给出了计算内力的边界轮廓法方程。该法无需进行数值积分计算,完全避免了角点问题和奇异积分计算。给出的算例,与解析解相比较,证实该方法的有效性。     

关于时域边界元法的若干研究

本文从三维弹性动力学方程的基本奇异解着手，导出了适于计算机计算的求解三弹性动力学问题的边界积分方程（ＢＩＥ），并在理论上提出了ＩＢＥ前缘系数矩阵（５）具有（１）准对角特性，（２）其各元素不随时间而变化。据此，本文给出了用时域边界元法求解弹性动力学问题的新方法。最后，数值算例验证了本文方法的正确性。     

动力分析中的边界元法(Ⅱ)

<正> 3 其他边界方法 第1,2两节较详细介绍了频域中和时域中的直接边界元法。这种边界元法称为直接的,是因为它讨论了具有重要物理意义的量,例如位移和拉力等,从而成为不仅在弹性动力学中而且在这种方法的所有其他力学应用中最广泛采用的边界元法。反之,以没有物理意义的某些虚拟量或源密度表述问题的间接边界元法是不太流行的,尽管事实上它具有比直接法更长的历史,而且与积分方程法方面的经典工作密切相关。     

概率断裂分析的随机边界元法

应用随机边界元法分析材料弹性常数的随机性和裂纹面随机性对应力强度因子的影响。文中首先简介了随机边界元法，给出了具有随机材料或几何参数的弹性体的边界位移或面力的协方差，进而给出了材料参数和裂纹面随机时应力强度因子均值和方差的计算公式。算例中详细讨论了杨氏模量、泊松比及裂纹面的随机性对应力强度因子的影响。     

薄板结构的边界元优化设计

完成了一种用边界元法对弹性薄板结构进行优化设计的计算机分析及彩色绘图系统.本文论述了平板弯曲问题的边界元法及其优化设计.系统采用本文作者提出的三方程协调方案,并精确计算分布荷载的域内积分.本文完成了多个平板弯曲问题的边界元分析算例及板结构边界元优化设计的工程算例.结果表明,本系统方法先进、结果精确可靠,具有十分明显的工程实用价值.     

边界元法应用的若干近期研究

边界元法是在有限元法等其它数值方法的推动和竞争下发展的,作为一种数值方法,发展其实际应用至关重要,本文扼要介绍作者们近期近作若干工作,其中包括：二维弹性体移动接触和滚动接触,网络要群环境下边界元并行计算,以及二维三维边界元实用软件及应用等。     

HC轧机辊间接触与辊系变形的Taylor级数多极边界元数学规划解析

分析HC轧机辊间接触分布和辊系弹性变形对于改善辊间压力分布状态,减少轧故褂檬倜案纳瓢逍畏浅Ｖ匾?醯捎诩扑懔亢艽?使用传统数值方法(有限元法或边界元法)分析辊间接触和辊系变形是非常困难的.本文描述了一种基于点-面接触模型的三维弹性接触Taylor级数多极边界元法,给出了数学规划解析方法,适合大规模弹性接触问题的求解....     

二维弹性随机边界元与可靠性分析

本文利用对随机变量求偏导数的方法，推导了二维弹性随机边界积分方程及其相应的公式，考虑了随机边界条件、材料性能参数随机量以及几何尺寸随机因素等。本文还将所发展的方法用于结构静强度的可靠性分析。算例分析表明，本文数值结果与ＭｏｎｔｅＣａｒｌｏ模拟和理论解相比是十分满意的。     

随机杆系结构几何非线性分析的递推求解方法

建立了随机静力作用下考虑几何非线性的随机杆系结构的随机非线性平衡方程. 将和 位移耦合的随机割线弹性模量以及随机响应量表示为非正交多项式展开式，运用传统的摄动方法获 得了关于非正交多项式展式的待定系数的确定性的递推方程. 在求解了待定系数后，利用非 正交多项式展开式和正交多项式展开式的关系矩阵，可以很方便地得到未知响应量的二阶统计矩. 两杆结构和平面桁架拱的算例结果表明，当随机量涨落较大时，递推随机有限元方法比基于 二阶泰勒展开的摄动随机有限元方法更逼近蒙特卡洛模拟结果，显示了该方法对几何非线性 随机问题求解的有效性.     

Nonlinear analysis of structural elements under unilateral contact constraints by a Ritz type approach

A nonlinear modal solution methodology capable of solving equilibrium and stability problems of uni-dimensional structural elements (beams, columns and arches) with unilateral contact constraints is presented in this work. The contact constraints are imposed by an elastic foundation of the Winkler type, where special attention is given to the case in which the foundation reacts in compression only, characterizing the contact as unilateral. A Ritz type approach with moveable boundaries, where the coordinates defining the limits of the contact regions are considered as additional variables of the problem, is proposed to solve this class of unilateral contact problems. The methodology is illustrated by particular problems involving beams, beam-columns and arches, and the results are compared with available results obtained by finite element and mathematical programming techniques. It is concluded that the Ritz type approach proposed is particularly suited for the analysis of structural problems where the number, but not the length, of the contact regions between the bodies are known a priori. Therefore, it can substitute in these cases finite element applications and be used as a benchmark for more general and complex formulations as well.     

A 3D FEM–BEM rolling contact formulation for unstructured meshes

This work presents a new formulation for solving 3D steady-state rolling contact problems. The convective terms for computing the tangential slip velocities involved in the rolling problem, are evaluated using a new approximation inspired in numerical fluid dynamics techniques for unstructured meshes. Moreover, the elastic influence coefficients of the surface points in contact are approached by means of the finite element method (FEM) and/or the boundary element method (BEM). The contact problem is based on an Augmented Lagrangian Formulation and the use of projection functions to establish the contact restrictions. Finally, the resulting nonlinear equations set is solved using the generalized Newton method with line search (GNMls), presenting some acceleration strategies as: a new and more simplified projection operator, which makes it possible to obtain a quasi-complementarity of the contact variables, reducing the number of contact problem unknowns, and using iterative solvers. The presented methodology is validated solving some rolling contact problems and analyzed for some unstructured mesh examples.     

Mapping forces in a 3D elastic assembly of grains

Our understanding of the elasticity and rheology of disordered materials, such as granular piles, foams, emulsions or dense suspensions relies on improving experimental tools to characterise their behaviour at the particle scale. While 2D observations are now routinely carried out in laboratories, 3D measurements remain a challenge. In this paper, we use a simple model system, a packing of soft elastic spheres, to illustrate the capability of X-ray microtomography to characterise the internal structure and local behaviour of granular systems. Image analysis techniques can resolve grain positions, shapes and contact areas; this is used to investigate the material's microstructure and its evolution upon strain. In addition to morphological measurements, we develop a technique to quantify contact forces and estimate the internal stress tensor. As will be illustrated in this paper, this opens the door to a broad array of static and dynamical measurements in 3D disordered systems.     

Stochastic representation of the mechanical properties of irregular masonry structures

A procedure for the stochastic characterization of the elastic moduli of plane irregular masonry structures is presented in this paper. It works in the field of the random composite materials by considering the masonry as a mixture of stones (or bricks) and mortars. Once that the elastic properties of each constituent are known (deterministically or stochastically), the definition of the overall masonry elastic properties requires the knowledge of the random field describing the irregular geometry distribution. This last one is obtained by a software, implemented ad hoc, that, starting from a colour digital photo of the masonry and using the instruments of the digital image processing techniques, gives the random features of this field in both the space and frequency domain. The definition of the stochastic properties of masonry structures may be very useful both for the application of the stochastic homogenization techniques and for the direct stochastic analysis of the structures.     

Block pivot methods for solving frictional contact problems

Based on elementary group theory, the block pivot methods for solving two-dimensional elastic frictional contact problems are presented in this paper. It is proved that the algorithms converge within a finite number of steps when the friction coefficient is “relative small”. Unlike most mathematical programming methods for contact problems, the block pivot methods permit multiple exchanges of basic and nonbasic variables. The project supported by the National Natural Science Foundation of China     

摩擦接触问题的比例边界等几何B可微方程组方法

摩擦接触问题是计算力学领域最具挑战性的问题之一,接触系统的泛函具有非线性、非光滑的特点,导致接触算法的收敛性与精确性难以保证.因此将比例边界等几何分析(scaled boundary isogeometric analysis,SBIGA)与B可微方程组(B dierential equation,BDE)相结合,提出了求解二维摩擦接触问题的比例边界等几何B可微方程组方法.在比例边界等几何坐标变换的基础上,通过虚功原理推导了关于边界控制点变量的接触平衡方程,表示成B可微方程组形式的接触条件可被严格满足,求解B可微方程组的算法的收敛性有理论保证.此比例边界等几何B可微方程组方法(SBIGA-BDE)只需在接触体边界进行等几何离散,使问题降低一维,能精确描述接触边界,并可通过节点插入算法进行真实接触区域的识别.此外,由于几何建模和数值分析使用相同的基函数,节约了划分网格的时间.以赫兹接触问题和悬臂梁摩擦接触问题为例,通过与解析解及数值计算软件ANSYS计算结果进行对比,验证了该方法求解二维摩擦接触问题的有效性及高精度等特点.