裂缝各向异性介质中P-SV转换波正演模拟

裂缝介质正演是研究裂缝性油藏的各向异性、裂缝检测方法、裂缝AVO反演等的前提。在裂缝方位各向异性(EDA)介质模型的基础上，重点研究了纵波震源激发的地震波所产生P-SV波的有限差分正演模拟方法。从各向异性单斜对称系统的刚性矩阵推导出的波动方程，可以满足EDA介质模型不同方位角入射的二维差分模拟要求。在有限差分模拟中，采用了普通的Kelly差分格式、二阶各向异性吸收边界条件和一般的纵波震源的解析形式，并应用修正的通量传输校(FCT)方法来消除频散。同时，对纵波震源的解析式进行了技术处理，取得了较好的效果。     

基于遗传算法的土钉支护结构优化设计研究

通过建立三维有限元模型,对土钉支护的变形和受力性能进行分析,得出土钉力的合理分布。在此基础上,考虑土钉支护的水平位移,采用遗传算法,对土钉支护的结构优化设计进行研究,建立了土钉支护结构的优化设计数学模型,并编制了相应的计算程序。通过算例分析,并与基于极限平衡分析的优化结果相比较,得出合理的土钉支护结构设计的参考结论。     

横向各向同性介质中的初至波旅行时计算

本文直接采用Thomsen导出的横向各向同性介质中的相速度、群角与相角、群速度与相速度的精确函数关系,提出了一种新的具有任意各向异性强度的横向各向同性介质中的初至波旅行时计算方法,它不同于Faria等人给出的采用群速度近似公式(Byun)的仅适应弱各向异性的初至波旅行时计算方法。该新算法在计算具有任意各向异性强度的横向各向同性介质中的初至波旅行时时,并不显著地增大计算时间。分析对比由新算法得到的初至波旅行时的等值线图和由交错网格有限差分波场模拟得到的波场快照表明,无论是对弱各向异性还是对强弱各向异性,新算法都具有较高的精度;随着各向异性强度的增大,新算法的计算结果与采用近似公式计算结果的差值也随之增大。     

高压射流冲击破碎岩石的有限元计算

高压射流冲击破岩是一个复杂的非线性问题。通过MSC．Marc建立模型，分别使用动态接触作为非线性冲击载荷模拟高压水射流冲击岩石。通过动力计算分析，根据岩石内部应力的变化使用Hoffman失效准则研究了岩石的破碎过程。分析表明，增加射流冲击速度可以提高射流破碎岩石的效率，当射流速度迭到某一临界值时，射流水锤作用使得岩石发生大块破碎。     

Bifurcation points and bifurcated branches by an asymptotic numerical method and Padé approximants

The aim of this work is to develop a reliable and fast algorithm to compute bifurcation points and bifurcated branches. It is based upon the asymptotic numerical method (ANM) and Padé approximants. The bifurcation point is detected by analysing the poles of Padé approximants or by evaluating, along the computed solution branch, a bifurcation indicator well adapted to ANM. Several examples are presented to assess the effectiveness of the proposed method, that emanate from buckling problems of thin elastic shells. Especially problems involving large rotations are discussed. Copyright © 2004 John Wiley & Sons, Ltd.     

A p‐version,first order shear deformation,finite element for geometrically non‐linear vibration of curved beams

A p‐version, hierarchical finite element for curved, moderately thick, elastic and isotropic beams is introduced. The convergence properties of the element are analysed and some results are compared with results published elsewhere or calculated using a commercial finite element package. It is verified that, with the proposed element, shear locking does not affect the computation of the natural frequencies and that low dimensional, accurate models are obtainable. Geometrically non‐linear vibrations due to finite deformations, which occur for harmonic excitations with frequencies close to the first three natural frequencies of vibration, are investigated using Newmark's method. The influence of the thickness, longitudinal inertia and curvature radius on the dynamic behaviour of curved beams are studied. Copyright © 2004 John Wiley & Sons, Ltd.     

Second-order accurate constraint formulation for subdivision finite element simulation of thin shells

We present a new method for enforcing boundary conditions within subdivision finite element simulations of thin shells. The proposed framework is demonstrated to be second-order accurate with respect to increasing refinement in the displacement and energy norm for simply supported, clamped, free and symmetric boundary conditions. Second-order accuracy on the boundary is consistent with the accuracy of subdivision-based approaches for the interior of a body. Our proposed framework is applicable to both triangular and quadrilateral refinement schemes, and does not impose any topological requirements upon the underlying subdivision control mesh. Several examples from an obstacle course of benchmark problems are used to demonstrate the convergence of the scheme. Copyright © 2004 John Wiley & Sons, Ltd.     

A direct stiffness analysis of a composite beam with partial interaction

The use of the conventional semi-analytical stiffness method in finite element analysis, in which interpolation polynomials are used to develop the stiffness relationships, leads to problems of curvature locking when beam-type elements are developed for composite members with partial interaction between the materials of which it is comprised. The curvature locking phenomenon that occurs for composite steel–concrete members is quite well reported, and the general approach to minimizing the undesirable ramifications of curvature locking has been to use higher-order polynomials with increasing numbers of internal nodes. This paper presents an alternate formulation based on a direct stiffness approach rather than starting from pre-defined interpolation polynomials, and which does not possess the undesirable locking characteristics. The formulation is based on a more general approach for a bi-material composite flexural member, whose constituent materials are joined by elastic shear connection so as to provide partial interaction. The stiffness relationships are derived, and these are applied to a simply supported and a continuous steel–concrete composite beam to demonstrate the efficacy of the method, and in particular its ability to model accurately both very flexible and very stiff shear connection that causes difficulties when implemented in competitive semi-analytical algorithms. Copyright © 2004 John Wiley & Sons, Ltd.     

Application of piece‐wise linear weight functions for 2D 8‐node quadrilateral element in contact problems

The present study is a continuation of our previous work with the aim to reduce problems caused by standard higher order elements in contact problems. The difficulties can be attributed to the inherent property of the Galerkin method which gives uneven distributions of nodal forces resulting in oscillating contact pressures. The proposed remedy is use of piece‐wise linear weight functions. The methods to establish stiffness and/or mass matrix for 8‐node quadrilateral element in 2D are presented, i.e. the condensing and direct procedures. The energy and nodal displacement error norms are also checked to establish the convergence ratio. Interpretation of calculated contact pressures is discussed. Two new 2D 8‐node quadrilateral elements, QUAD8C and QUAD8D, are derived and tested in many examples, which show their good performance in contact problems. Copyright © 2004 John Wiley & Sons, Ltd.