Solving thermal and phase change problems with the eXtended finite element method

 The application of the eXtended finite element method (X-FEM) to thermal problems with moving heat sources and phase boundaries is presented. Of particular interest is the ability of the method to capture the highly localized, transient solution in the vicinity of a heat source or material interface. This is effected through the use of a time-dependent basis formed from the union of traditional shape functions with a set of evolving enrichment functions. The enrichment is constructed through the partition of unity framework, so that the system of equations remains sparse and the resulting approximation is conforming. In this manner, local solutions and arbitrary discontinuities that cannot be represented by the standard shape functions are captured with the enrichment functions. A standard time-projection algorithm is employed to account for the time-dependence of the enrichment, and an iterative strategy is adopted to satisfy local interface conditions. The separation of the approximation into classical shape functions that remain fixed in time and the evolving enrichment leads to a very efficient solution strategy. The robustness and utility of the method is demonstrated with several benchmark problems involving moving heat sources and phase transformations. Received 20 May 2001 / Accepted 19 December 2001     

Material state remapping in computational solid mechanics

A computational procedure for remapping material state information from one finite element mesh to another is described. The procedure is useful in connection with evolving meshes for inelastic problems, as for example occur in the context of fracture simulation and adaptive mesh refinement. The proposed method is based on weak enforcement of equality between corresponding fields on the two meshes, where piecewise‐constant fields on both meshes are generalized from the quadrature‐point values. The essential algorithmic problem is that of calculating the volume partition of an arbitrary convex region with respect to a covering set of disjoint convex regions. Instead of geometrically resolving the associated intersections, the problem is herein approximated by a constrained optimization problem, which may be readily and efficiently solved computationally. This formulation is a main contribution of the paper. Computational examples are given that illustrate the effectiveness of the proposed procedure. Copyright © 2002 John Wiley & Sons, Ltd.     

铝电解槽磁流体稳定性有限元分析

本文采用有限元法求解分层磁流体小扰动方程系统特征值，研究铝电解槽稳定性。研究结果表明，在没有电磁力时，铝液表面波动（内波）是稳定的。存在电磁力时，它会激发铝液表面波动的某些低频长波分量不稳定性，高频短波不受影响。被激发的不稳定波动频率和扰动增长率与电流密度，磁感应强度垂直分量，铝电解槽长宽比，电解质和铝液的密度差，厚度等因素密切相关。减少电流密度，磁感应强度，增加电极距离和铝液厚度可以提高铝电解槽的稳定性。这对铝电解槽优化设计和在线控制有重要意义。     

Lower bound limit analysis using non‐linear programming

This paper describes a new formulation, based on linear finite elements and non‐linear programming, for computing rigorous lower bounds in 1, 2 and 3 dimensions. The resulting optimization problem is typically very large and highly sparse and is solved using a fast quasi‐Newton method whose iteration count is largely independent of the mesh refinement. For two‐dimensional applications, the new formulation is shown to be vastly superior to an equivalent formulation that is based on a linearized yield surface and linear programming. Although it has been developed primarily for geotechnical applications, the method can be used for a wide range of plasticity problems including those with inhomogeneous materials, complex loading, and complicated geometry. Copyright © 2002 John Wiley & Sons, Ltd.     

Signed line intercept count

Elements of surfaces that bound a phase (β) in a two phase mixture (α+β) may be classified as: (a) convex (++) if both principal curvatures are positive; (b) concave (? ?) if both are negative; and (c) saddle (+ ?) if one is positive and the other negative. This classification excludes the limiting cases for which one or both of the principal curvatures is zero. The traces of these surfaces that form the boundaries of the β areas on a representative two dimensional section may also be: (a) convex (+) if the local curvature is positive; or (b) concave (?) if it is negative. Line intercepts may be tabulated separately for intersections with convex (+) and concave (?) segments of boundary. This paper presents a derivation of fundamental stereological formulae that relate these counting measurements to three-dimensional geometric properties of the structure they sample.     

隔河岩水电厂220 kV线路保护配置及运行分析

介绍了隔河岩电厂220kV线路保护的原理设计及运行情况，以及采用了双重微机化继电保护装置后带来的优越性，简要描述了CSL-101A高频保护和RCS-931A光纤保护的应用特点。     

海洋石油平台水下夹桩器本体结构分析

水下夹桩器是一种广泛应用于海洋石油平台建设的液压夹具。水下夹桩器本体焊接在导管架上，支撑夹桩器液压缸。在水下夹桩器本体设计中，采用有限元法完成了对本体强度和刚度的结构分析；建立了有限元分析模型，用ANSYS软件进行结构分析，得到了分析结果，并将该结果与理论计算结果进行了比较，验证了有限元分析技术的有效性。设计结果满足工程要求，为改进和优化结构提供了理论依据，节省了分析计算时间。     

单桩腿有限元分析

为保证桩腿在冬季冰区海域的安全使用,对桩腿结构进行冰载荷作用下的强度分析。应用ANSYS软件建立桩腿、地基、冰相互作用的仿真计算模型,实现桩腿的力学分析。对桩腿进行了强度校核,计算出冰载荷作用下桩腿的变形及其应力状态。计算结果表明,桩腿在此冰载荷作用下强度满足安全要求。     

潜孔锤冲击作用下土体变形的有限元分析

潜孔锤冲击挤密钻进过程中周围土体的变形机理是一个高度非线性问题。针对潜孔锤在土中的面-面接触问题,运用Drucker-Prager屈服准则为判据,将土体视为DP材料,通过钻具与土之间的接触及施加冲击载荷来模拟钻进过程,从而建立了潜孔锤在土中冲击挤密钻进的有限元模型。利用非线性瞬态动力学有限元分析手段,对钻进过程中冲击力对土的作用加以研究,分析了钻具周围土体的应力、应变以及钻具与土体的摩擦情况。结果表明,在冲击载荷作用下,锥形钻头前端土体应力集中,并形成应力泡,冲击功主要用于竖向压缩土体,而水平方向的挤压作用相对较小,摩擦应力在钻头锥面附近出现峰值,是影响钻进速度的主要原因之一。     

两相电流互感器配电网零序电流构造方法

中性点非直接接地配电网大多只安装A，C两相电流互感器（TA），而现有单相接地故障选线方法基本都以零序电流为基础实现。如何在两相TA配电网中准确检测并隔离接地故障线路，成为长期困扰实际运行的技术难题。通过对暂态电流的频域和时域分析，找出B相暂态电流与A，C相暂态电流的对应关系，提出在A相和C相接地时由A相和C相暂态电流构造B相暂态电流，进而得到较精确零序电流的方法。通过大量仿真实验验证了所提出的构造方法的正确性。零序电流的成功构造为实现两相TA配电网单相接地故障选线提供了基础。