桩基完整性测试的非线性最小二乘优化方法

针对桩基检测中低应变反射波法存在的应力波被桩土吸收和散射而衰减问题,把包含桩土衰减系数和介质不连续界面反射系数的反射波方程转化为非线性最小二乘优化问题,并采用信赖域法进行求解。由于采用信赖域法并采用了Gauss-Newton方向定义信赖域二维子空间,具有较好的收敛速度。通过模型桩算例表明该方法是确实可行的。     

正交各向异性压磁电弹性矩形板的稳定问题

从压磁电弹性力学的三维平衡方程和本构关系出发,推导得到了正交各向异性压磁电弹性矩形板稳定问题的状态方程。求解状态方程并结合边界条件得到了稳定问题的屈曲方程。给出了临界应力的算例。     

平面弹性力学问题的离散元法

根据离散元的基本原理,基于变形体的理论提出了适用于平面弹性力学问题的界面位移、应变和应力模式,建立了求解平面弹性力学问题的离散元方程和相应的迭代求解方法.通过界面位移可以简洁地将位移和力的边界条件引入离散系统的控制方程,也可以方便地求解节点位移.数值算例表明,与具有相同网格的有限元结果相比,离散元能同时给出精度相对较高的应力解和精度相当的位移解.     

A Variational Approach to the Fredholm Alternative for the p-Laplacian near the First Eigenvalue

We are concerned with the existence of a weak solution to the degenerate quasi-linear Dirichlet boundary value problem

Nonlinear interactions and chaotic dynamics of suspended cables with three-to-one internal resonances

Nonlinear planar oscillations of suspended cables subjected to external excitations with three-to-one internal resonances are investigated. At first, the Galerkin method is used to discretize the governing nonlinear integral–partial-differential equation. Then, the method of multiple scales is applied to obtain the modulation equations in the case of primary resonance. The equilibrium solutions, the periodic solutions and chaotic solutions of the modulation equations are also investigated. The Newton–Raphson method and the pseudo-arclength path-following algorithm are used to obtain the frequency/force–response curves. The supercritical Hopf bifurcations are found in these curves. Choosing these bifurcations as the initial points and applying the shooting method and the pseudo-arclength path-following algorithm, the periodic solution branches are obtained. At the same time, the Floquet theory is used to determine the stability of the periodic solutions. Numerical simulations are used to illustrate the cascades of period-doubling bifurcations leading to chaos. At last, the nonlinear responses of the two-degree-of-freedom model are investigated.     

Buckling delamination in compressed multilayers on curved substrates with accompanying ridge cracks

Residually compressed films and coatings are susceptible to buckle delamination. The buckles often have linear or telephone cord morphology. When the films are brittle, such buckles are susceptible to the formation of ridge cracks that extend along their length, terminating close to the propagating front. The ridge-cracked buckles are invariably straight-sided (not telephone cord) and differ in width. Buckle delaminations of this type occur on flat and curved substrates: having greatest technological relevance in the latter. They occur not only in single layer films but also in multilayers, such as thermal barrier systems. Establishing the mechanics of ridge-cracked buckle delaminations for multilayers on curved substrates serves two purposes. (a) It allows the prediction of buckle delamination and spalling for technologically important systems. (b) It provides a test protocol for measuring properties such as the delamination toughness of the interface and the stresses in the layers. Both objectives are addressed in the article: the latter by devising an inverse algorithm. Implementation of the algorithm is demonstrated for diamond-like carbon films on planar glass substrates and a thermal barrier multilayer on a curved superalloy substrate.     

含孔von Karman板中非线性波散射与边值问题

基于ｖｏｎ Ｋａｒｍａｎ板大挠度弯曲理论,利用小参数摄动法,分析研究了含孔ｖｏｎＫａｒｍａｎ板的非线性波散射与动应力集中问题,其中一类可看成是薄板弯曲波动问题的控制方程。当有单频波入射时,由于弯曲应力与膜应力状态的非线性耦合,孔洞会产生高次谐波散射现象。建立了求解本问题的边界积分方程法,利用积分方程法交替求求这两类问题,最终可获得问题的近似分析解。     

矩形区域不等圆形随机堆放的数值模拟

将若干待堆放的圆形堆放到一个给定矩形区域,是一类特殊的Packing问题。针对这一具有NP难度的困难问题,提出一种快速的数值模拟方法以及模拟力学检验的方法。该方法将圆形堆放的力学平衡规律化为几何形体相互接触的几何条件,把大量的圆形按级配及分布律等约束条件,随机堆放形成空间区域。按此方法得到的圆形堆体,在不加外力的情形下即可保持自身的稳定。使用该方法,在边长为15cm的正方形区域内,取圆形的直径范围为0.5cm~4cm,模拟了孔隙率小于25%(面积比)的圆形堆体。最后,将所模拟的圆形堆体试件进行了有限元网格划分及力学加载试算,结果表明该堆体是稳定的。为进一步在细观层次研究圆形堆体提供了一个快速的数值模拟方法。     

An efficient method for solving the mixed direct-inverse problem of the transonic rotational flow in plane casades

The method in [1] has been extended to the case of rotational flow in this paper. A new method for dealing with the shock wave is presented. This method has the advantages of both the shock-fitting and the shock capturing methods. The direct problem and the mixed direct-inverse problem of the rotational flow in a transonic plane cascade at both design and off design conditions are solved, and the results show that the present method has rapid convergence rate and high accuracy even for the flow with moderately strong shocks. The calculations have been carried out on the DPS-8 computer, and for the direct problem, only 50–80 iterations are needed, and 50–80 seconds of CPU time are required.