基于结构动力修改的桩基检测方法

采用结构动力修改技术研究桩身损伤检测，将其转化为在一定约束条件下拟合桩顶频响函数的优化问题。用遗传算法加以求解。文中给出的数值模拟试验和工程应用实例，显示该方法能在较大范围内可靠地反演出桩身横截面积、桩长、波速等参数，达到检测桩身缺陷的目的，实现测试和分析一体化，使桩基检测易于进行并且结果实用可靠。     

Deteriorated Geometrical Stiffness for Higher Order Finite Elements with Application to Panel Flutter

Higher order elements were first design for linear problems where, in certain situations, they present advantages over the lower order elements. A method to efficiently extend their use to geometrical nonlinear problems as panel flutter and postbuckling behavior is presented. The chaotic and limit-cycle oscillations of an isotropic plate are obtained based on direct integration of the discretized equation of motion. The plate is modeled using the von Karman theory and the geometrical nonlinearities are separated in a nonlinear term of the first kind which manifests especially in the prebuckling and buckling regimes, and a nonlinear term of the second kind which is responsible for the postbuckling behavior. A fifth order, fully compatible element has been used to model thin plates while the inplane loads where introduced through a membrane element. The aerodynamics was modeled using the first order 'piston theory. The method introduces the concept of a deteriorated form of the second geometric matrix which is equivalent to neglecting higher order terms in the strain energy of the plate. This allows for a drastic reduction in the computational effort with no observable loss of accuracy. Well established results in the literature are used to validate the method.     

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

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

新型各向异性核工业石墨断裂力学行为的实验研究和数值分析

本文分别采用激光和白光DSCM(数字散斑相关测量)方法对一种新型各向异性核工业石墨的裂纹尖端位移、应变场进行了实验研究，两种方法都取得了较好的结果。考虑核工业石墨制备过程形成的各向异性特点，本文构建了三维各向异性有限元模型，采用奇异单元，计算模拟石墨裂纹尖端的变形和应力场，通过对实验结果和有限元计算结果的比较可以发现两者具有相近的趋势。     

Numerical study of sliding wear caused by a loaded pin on a rotating disc

A computational approach is proposed to predict the sliding wear caused by a loaded spherical pin contacting a rotating disc, a condition typical of the so-called pin-on-disc test widely used in tribological studies. The proposed framework relies on the understanding that, when the pin contacts and slides on the disc, a predominantly plane strain region exists at the centre of the disc wear track. The wear rate in this plane strain region can therefore be determined from a two dimensional idealisation of the contact problem, reducing the need for computationally expensive three dimensional contact analyses. Periodic unit cell techniques are used in conjunction with a ratchetting-based failure criterion to predict the wear rate in the central plane strain region. The overall three dimensional wear rate of the disc is then determined by scaling the plane strain wear rate with a conversion factor related to the predicted shape of the wear track. The approach is used to predict pin-on-disc test data from an Al-Si coating using a tungsten carbide pin. The predicted results are found to be consistent with measured data.     

Finite element simulations of martensitic phase transitions and microstructures based on a strain softening model

A new approach for modeling multivariant martensitic phase transitions (PT) and martensitic microstructure (MM) in elastic materials is proposed. It is based on a thermomechanical model for PT that includes strain softening and the corresponding strain localization during PT. Mesh sensitivity in numerical simulations is avoided by using rate-dependent constitutive equations in the model. Due to strain softening, a microstructure comprised of pure martensitic and austenitic domains separated by narrow transition zones is obtained as the solution of the corresponding boundary value problem. In contrast to Landau-Ginzburg models, which are limited in practice to nanoscale specimens, this new phase field model is valid for scales greater than 100 nm and without upper bound. A finite element algorithm for the solution of elastic problems with multivariant martensitic PT is developed and implemented into the software ABAQUS. Simulated microstructures in elastic single crystals and polycrystals under uniaxial loading are in qualitative agreement with those observed experimentally.     

桥梁空间分析中预应力效应分析方法研究

预应力技术在桥梁中的广泛使用,使得桥梁结构受力状况复杂,对预应力混凝土桥梁进行空间分析是必要的;但由于预应力筋数量庞大且形状复杂,对预应力空间效应的模拟是预应力混凝土桥梁空间分析的难点。从预应力筋的实际作用机理出发,提出在预应力张拉阶段用等效节点荷载来模拟预应力对结构的作用;而在预应力管道灌浆后,将预应力筋作为结构的一部分,用实体退化单元进行结构分析,并用两算例对本文的方法进行了验证。结果表明,该方法能真正模拟预应力的空间效应,可用于预应力混凝土桥的空间分析。     

Vortex-induced oscillations at low Reynolds numbers: Hysteresis and vortex-shedding modes

Results are presented for the numerical simulation of vortex-induced vibrations (VIVs) of a cylinder at low Reynolds numbers (Re). A stabilized space–time finite-element formulation is utilized to solve the incompressible flow equations in primitive variables. The cylinder, of low nondimensional mass (m^*=10), is free to vibrate in, both, the transverse and in-line directions. To investigate the effect of Re and reduced natural frequency, F_n, two sets of computations are carried out. In the first set of computations the Reynolds number is fixed (=100) and the reduced velocity (U^*=1/F_n) is varied. Hysteresis, in the response of the cylinder, is observed at the low- as well as high-end of the range of reduced velocity for synchronization/lock-in. In the second set of computations, the effect of Reynolds number (50Re500) is investigated for a fixed reduced velocity (U^*=4.92). The effect of the Reynolds number is found to be very significant for VIVs. While the vortex-shedding mode at low Re is 2S (two single vortices shed per cycle), at Re300 and larger, the P+S mode of vortex shedding (a single vortex and one pair of counter-rotating vortices are released in each cycle of shedding) is observed. This is the first time that the P+S mode has been observed for a cylinder undergoing free vibrations. This change of vortex-shedding mode is hysteretic in nature and results in a very large increase in the amplitude of in-line oscillations. Since the flow ceases to remain two-dimensional beyond Re200, it remains to be seen whether the P+S mode of shedding can actually be observed in reality for free vibrations.     

Constitutive equations for martensitic reorientation and detwinning in shape-memory alloys

A crystal-inelasticity-based constitutive model for martensitic reorientation and detwinning in shape-memory alloys (SMAs) has been developed from basic thermodynamics principles. The model has been implemented in a finite-element program by writing a user-material subroutine. We perform two sets of finite-element simulations to model the behavior of polycrystalline SMAs: (1) The full finite-element model where each finite element represents a collection of martensitic microstructures which originated from within an austenite single crystal, chosen from a set of crystal orientations that approximates the initial austentic crystallographic texture. The macroscopic stress-strain responses are calculated as volume averages over the entire aggregate: (2) The Taylor model (J. Inst. Metals 62 (1938) 32) where an integration point in a finite element represents a material point which consist of sets of martensitic microstructures which originated from within respective austenite single-crystals. Here the macroscopic stress-strain responses are calculated through a homogenization scheme.Experiments in tension and compression were conducted on textured polycrystalline Ti-Ni rod initially in the martensitic phase by Xie et al (Acta Mater. 46 (1998) 1989). The material parameters for the constitutive model were calibrated by fitting the tensile stress-strain response from a full finite-element calculation of a polycrystalline aggregate to the simple tension experiment. With the material parameters calibrated the predicted stress-strain curve for simple compression is in very good accord with the corresponding experiment. By comparing the simulated stress-strain response in simple tension and simple compression it is shown that the constitutive model is able to predict the observed tension-compression asymmetry exhibited by polycrystalline Ti-Ni to good accuracy. Furthermore, our calculations also show that the macroscopic stress-strain response depends strongly on the initial martensitic microstructure and crystallographic texture of the material.We also show that the Taylor model predicts the macroscopic stress-strain curves in simple tension and simple compression reasonably well. Therefore, it may be used as a relatively inexpensive computational tool for the design of components made from shape-memory materials.