水力射孔对套管强度的影响研究

水力射孔技术是一种新型的完井方式，深穿透水力射孔技术辅助定向水力压裂可以实现油气层改造和油气井增产。但是，水力射孔套管会在孔眼处造成应力集中，使得局部应力过大而削弱了套管和水泥环的强度，可能会影响套管的使用寿命。因此，文章利用有限元软件建立了射孔套管的三维有限元模型，分析了水力射孔后套管管体的应力分布，重点研究了水力射孔参数-孔密、孔径、孔深、射孔方位角等对套管强度的影响规律。分析结果表明，沿着最大水平地应力方向深穿透水力射孔，选择合适的孔密、孔径可将射孔对套管的损害降至最小。研究结果为提高水力射孔完井设计水平提供了一定的参数依据。     

燃气管网动态仿真的研究及应用

研究和了解燃气在管道中的流动，是进行燃气管网系统优化和提高管网输运系统安全性的前提和基础。以流体力学三大守恒方程为基础，建立了等温和非等温条件下的燃气管网稳态和动态仿真理论模型。以有限差分法为基础，得到了上述模型的求解方法，并加入了耗散项，提高了求解方法的稳定性。通过增加初始点和延长出口点的方法，简化了模型的求解。通过在北京六环天然气管网中的初步应用，分别得到了与实际值综合相对误差为1.21％和2.62％的压力和流量仿真值，从而验证了模型及求解方法的可靠性。     

波动方程数值模拟的隐式差分法

在有限差分波动方程数值模拟中,通常采用高阶差分方法来提高空间导数的数值逼近精度,以实现降低数值频散,提高数值模拟精度的目的。首先对差分频散进行了理论分析;然后讨论了估计一阶空间导数的隐式差分格式,并与通常采用的高阶精度显式差分格式进行了对比分析,结果表明,隐式差分格式能够在更宽的波数范围使差分频散控制在可接受的水平,如8阶精度的显式差分格式所适应的波数带宽约为O．55kmax,而隐式差分格式所适应的波数带宽约为0．7kmax;最后通过模型试算,对隐式差分格式的有效性进行了验证。模拟结果表明,用隐式差分格式在一定程度上降低了差分频散,提高了模拟精度。     

基于PC机群的大地电磁Occam反演并行计算研究

Occam反演以其稳定收敛和不依赖于初始模型的特性被广泛应用于大地电磁数据的处理，但偏导数矩阵的计算和拉格朗日乘子的求取导致大量的模型正演，使得反演速度较低。为此，研究了用基于PC机群的并行计算来解决这一问题的方法。首先对Occam反演方法进行了阐述；然后，分析了反演方法中各计算耗费的时间，提出对约占计算量90％的偏导数和拉格朗日乘子进行并行计算的思想，即偏导数计算采用频点计算一级的大粒度并行，拉格朗日乘子扫描和一维搜索分别采用弘值计算一级的大粒度并行和频点计算一级的小粒度并行；给出了并行计算的实现方法，即在PC机群上，利用主一从编程模式实现Occam反演的整体并行计算。在计算中，采用任务组合方式，减少了通信量，较好地实现了负载均衡。在4节点PC机群上，对应于拉格朗日乘子的扫描和一维搜索方式，整体加速比分别达到3．36和2．83。     

丙烯腈装置回收塔内聚合物生成原因的分析

针对丙烯腈生产过程中回收塔内聚合物堵塞降液管造成回收塔运行周期短的问题,通过对丙烯腈回收塔内各组分分布的模拟计算,直观显示了回收塔内各组分的浓度分布,发现了回收塔的聚合段与丙烯醛的富集段相吻合的现象,初步得出了回收塔内局部聚合的原因;通过对塔内聚合物的特性分析和评判,进一步证实了“回收塔内局部聚合严重是由于在聚合段丙烯醛聚集引起”这一现象;根据模拟计算结果和实际运行经验,提出了预防聚合物堵塞降液管的措施,并进行了部分生产验证,为丙烯腈生产企业延长装置运行周期提供了理论指导。     

大型弧形闸门静力特性有限元分析

以某水电站溢洪道弧形闸门为例,建立了大型弧形闸门的有限元分析模型,应用ANSYS软件对其在各种工况下的支铰反力、闸门应力、位移挠度以及屈曲失稳等静力特性进行了计算分析,为完善闸门结构的设计提供了依据,对同类闸门结构的设计计算有参考意义。     

Fatigue crack propagation of multiple coplanar cracks with the coupled extended finite element/fast marching method

A numerical technique for modeling fatigue crack propagation of multiple coplanar cracks is presented. The proposed method couples the extended finite element method (X-FEM) [Int. J. Numer. Meth. Engng. 48 (11) (2000) 1549] to the fast marching method (FMM) [Level Set Methods & Fast Marching Methods: Evolving Interfaces in Computational Geometry, Fluid Mechanics, Computer Vision, and Materials Science, Cambridge University Press, Cambridge, UK, 1999]. The entire crack geometry, including one or more cracks, is represented by a single signed distance (level set) function. Merging of distinct cracks is handled naturally by the FMM with no collision detection or mesh reconstruction required. The FMM in conjunction with the Paris crack growth law is used to advance the crack front. In the X-FEM, a discontinuous function and the two-dimensional asymptotic crack-tip displacement fields are added to the finite element approximation to account for the crack using the notion of partition of unity [Comput. Meth. Appl. Mech. Engng. 139 (1996) 289]. This enables the domain to be modeled by a single fixed finite element mesh with no explicit meshing of the crack surfaces. In an earlier study [Engng. Fract. Mech. 70 (1) (2003) 29], the methodology, algorithm, and implementation for three-dimensional crack propagation of single cracks was introduced. In this paper, simulations for multiple planar cracks are presented, with crack merging and fatigue growth carried out without any user-intervention or remeshing.     

Prediction of brittle-to-ductile transitions in polystyrene

In this study it is attempted to predict brittle-to-ductile transitions (BDTs) in polystyrene blends, induced either by an increase in temperature or by a decrease in inter-particle distance. A representative, two-dimensional volume element (RVE) of a polystyrene matrix with 20% circular voids, is deformed in tension. During deformation a hydrostatic-stress based craze-nucleation criterion [1] is evaluated. The simulations demonstrate that crazes initiate at low temperatures while a transition from crazing to shear yielding (BDT) is found around 75 °C. The numerical results correlate well with tensile tests on similar heterogeneous polystyrene. The presence of an absolute length, as experimentally found, is more difficult to explain. Near a free surface a T_g-depression is measured for polystyrene and also the resistance to indentation in polystyrene is lower than expected from bulk properties. Both observations are rationalised by an enhanced segmental mobility of chains near a free surface. As a consequence of these findings, an absolute length-scale could be incorporated in the numerical simulations. For simplicity, the length-scale is modelled by taking a temperature gradient over a thin layer near the internal free surfaces of the RVE. Deformation of the RVE with different absolute length-scales shows that indeed also the experimentally found brittle-to-ductile transition can be predicted if the ligament thickness between the inclusions (‘voids’) in polystyrene is below a critical value of ca. 15 nm.     

电信企业信息化统一模型构建策略

电信企业信息化的实质就是赋予电信企业精细化的企业资源整合能力，本以此认识为出发点，通过定义基本的信息化业务元素，需求参数和基本信息处理机制，构建出带有行业共性的电信企业信息化统一模型，该模型针对目前电信企业信息化中常见的和潜在的问题提供了有效的解决措施。     

面积坐标法构造含转角自由度的四结点膜元

以四边形面积坐标作为工具,构造了两个含转角自由度的广义协调四边形单元AQ4和lAQ4。它们通过强式分片检验,与同类单元相比,具有很高的计算精度,能消除梯形闭锁现象,有很强的抗网格畸变的能力。