A comprehensive study on a propagating buckle in externally pressurized pipelines

Buckle propagation is a unique phenomenon occurring in deep-sea pipelines. In previous works, this phenomenon was investigated using a ring technique in which the pipeline was assumed to be in plane strain condition and the energies absorbed in membrane stretching and longitudinal bending were ignored. This paper presents a three-dimensional analysis of the buckle propagation phenomenon with an emphasis to address more complete factors that were not accounted for in the ring analysis. The analyses are based on the available solutions of the transition zone obtained in our previous works. A comprehensive mechanism for buckle propagation phenomenon is described from the point view of plastic stability theory for shells which enables the incorporation of the effects of transverse and longitudinal bending, membrane stretching and material strain hardening. The nondimensionalized buckle propagation pressure is represented in terms of yield coefficient, strain hardening coefficient and membrane stretching factor. It is found that a buckle once initiated in a pipeline may or may not propagate along the pipeline depending on its radius-to-thickness ratio. By comparing with various experimental results the theoretical predictions from this analysis are shown to provide very accurate estimations of the buckle propagation pressure for different materials with diverse geometric parameters and material properties. This paper points to the need for more complete information regarding the effects of transverse bending, membrane stretching and material strain-hardening on the buckle propagation pressure. Upon the requirement of application variations of the yield coefficient, strain hardening coefficient and membrane stretching factor with respect to the radius-to-thickness ratio are sketched out. This eliminates the need for recourse the curves and allows a fast and convenient resolution of buckle propagation pressure for certain pipeline. Most importantly, the present analysis offers the potential for future design of pipelines being at once more rationally and parametrically complete, and yet compact and simple to apply.     

On the effect of axial tension on the propagation pressure of long cylindrical shells

Local imperfections induced to long tubes subjected to high external pressures can lead to local collapse from which a propagating buckle can be initiated. This can result in catastrophic collapse of large sections of the structure. The lowest pressure at which such a buckle will propagate, known as the propagation pressure of the tube, is typically half an order of magnitude lower than the collapse pressure of the intact tube. In a number of modern deep water applications, long tubular structures are subject to high axial tension in addition to external pressure. The paper describes the results of an experimental study in which the propagation pressure of long metal tubes was measured in the presence of a constant, axial tensile force. Tension was found to significantly reduce the propagation pressure. A parametric study of the problem, augmented by a simple model of the phenomenon, has yielded approximate expressions for the propagation pressure in the presence of axial tension.     

Modélisation,identification et simulations éléments finis des phénomènes de délaminage dans les structures composites stratifiées

Modelling, identification and finite element predictions of delamination in laminated composite structures. In order to forecast the delamination initiation and propagation in a finite element context, a previously defined damage meso-modelling of composite laminates is used. At the meso-level, the laminate is described as an assembly of damageable layers and interlaminar interfaces. The present work concerns the modelling, the identification and the finite element predictions of delamination phenomena in composite structures. Finite element predictions are conducted with non-linear geometric and material hypothesis. The interface modelling is implemented in the finite element code Castem 2000 developped by CEA. Classic edge delamination tension and propagation tests are conducted in order to improve the interface damage model approach in the finite elements context.     

On the propagation pressure of long cylindrical shells under external pressure

Local imperfections induced in long tubes subjected to high external pressures can lead to local collapse, from which a propagating buckle can be initiated. This can result in catastrophic collapse of large sections of the structure. The propagation pressure is the lowest pressure at which such a buckle will propagate. For common structural metal tubes with diameter-to-thickness ratios of less than 100, the propagation pressure is typically half an order of magnitude lower than the collapse pressure of the intact tubes. As a result, the design of several tubular structures with external pressure loading requires that the collapse and propagation pressures be accurately known. This paper deals with the experimental and analytical challenges of establishing the propagation pressure. A special purpose three-dimensional analysis, in combination with experimental observations and results, is used to demonstrate a mechanism of initiation of propagating buckles in long tubes, to study the parametric dependence of the propagation pressure and to illustrate the effect of axial tension on the propagation pressure. Propagation pressures predicted with this analysis are used to evaluate the strengths and weaknesses of simpler models developed in the past.     

基于导波多尺度能量熵的钢绞线张拉力识别

为了有效评估在役钢绞线中预应力的损失程度,提出一种识别钢绞线张拉力的导波无损检测方法。在不同张拉力钢绞线上进行了超声导波传播数值模拟及实验,构建以导波多尺度能量熵为特征向量的识别指标进行钢绞线张拉力识别,并考虑了导波传播距离及采集方式的影响。结果表明:不同拉力作用下,导波多尺度能量熵差异显著,识别指标与钢绞线张拉力存在明显的线性关系;识别指标敏感性随导波传播距离的增加而增强,但其提升幅度随传播距离的增加而减小。与有限元模拟结果相比,实测识别指标传播距离增加957.69%,敏感性系数提高了20.3%。采用中心钢丝激励、中心钢丝接收导波的采集方式进行张拉力识别,敏感性更优;识别指标与张拉力之间的线性变化规律受导波采集方式的影响较小。     

Fracture criteria identification using an inverse technique method and blanking experiment

In order to optimize the blanking processes, it is important to identify the conditions within the deforming workpiece which may lead to fracture initiation and propagation. Within this framework, numerical simulations are widely used in industries to optimize sheet metal forming processes. However, in order to have a confidence in the results of such simulations, an accurate material model is required. The accuracy of a material model is affected by the constitutive equations and the values of the material parameters. In order to reduce the danger of fracture of metal parts during manufacturing processes, advanced optimal design requires knowledge of critical values of some fracture criteria of the material used. Experimental identification of fracture criteria are currently obtained by performing several complicated tests and long duration of experiments.This study presents a computation methodology allowing for the identification of critical values of fracture criteria in order to simulate crack initiation and propagation generated by shearing mechanisms, which are needed for metal blanking processes simulation. The approach is based on inverse technique using circular blanking experiments and finite element calibration model. The critical values of fracture criteria are obtained in such a way that the finite element force–penetration predicted curve fit the experimental plot deduced from blanking tests. The numerical results obtained by the simulation were compared with experimental ones to verify the validity of the proposed technique for fracture criteria identification.     

Integral buckle arrestors for offshore pipelines: Enhanced design criteria

Integral buckle arrestors are relatively thick wall rings periodically welded in an offshore pipeline at intervals of several hundred meters in order to safeguard the line in case a propagating buckle initiates. They provide additional circumferential rigidity and thus impede downstream propagation of collapse, limiting the damage to the length of pipe separating the two arrestors. The effectiveness of such devices was studied parametrically through experiments and numerical simulations in Park and Kyriakides [On the design of integral buckle arrestors for offshore pipelines. International Journal of Mechanical Sciences 1997;39(6):643–69]. The experiments involved quasi-static propagation of collapse towards an arrestor, engagement of the arrestor, temporary arrest, and the eventual crossing of collapse to the downstream pipe at a higher pressure. The same processes were simulated with finite element models that included finite deformation plasticity and contact. The experimental crossover pressures enriched with numerically generated values were used to develop an empirical design formula for the arresting efficiency of such devices. A recent experimental extension of this work revealed that for some combinations of arrestor and pipe yield stresses, the design formula was overly conservative. Motivated by this finding, a new broader parametric study of the problem was undertaken, which demonstrated that the difference between the pipe and the arrestor yield stress affects significantly the arrestor performance. The original arrestor design formula was then modified to include the new experimental and numerical results producing an expression with a much wider applicability.     

Dynamic crack growth in TDCB specimens

Dynamic crack propagation in tapered double cantilever beam (TDCB) specimens is analysed via beam theory and the finite element method. Steady state and transient solutions of the energy release rate G are given for various load conditions. Finite element analysis is performed to obtain the dynamic G at given crack speed or the crack history for a given fracture toughness. The stress wave effects on the dynamic G are discussed. The beam solutions are compared with the finite element results and some experimental phenomena are explained.     

Fatigue growth simulation for cracks in notched and unnotched round bars

Fatigue crack growth for various cracks in both unnotched and notched round bars is directly modelled by using an automated numerical technique, which calculates the stress intensity factors at a set of points on the crack front through the three-dimensional finite element method and then applies an appropriate fatigue crack growth law to this set of points to obtain a new crack front. This technique also has the capability of automatic remeshing so that the crack propagation can conveniently be followed. The crack geometries modelled in the present paper include a small internal crack near the free boundary and several initially part-elliptical or irregular surface cracks in a smooth round bar under tension, a surface crack in different semi-circularly notched bars under both tension and bending, a surface crack initiated from the root of a V-notched bar and an initially twin crack configuration within a smooth tension bar. Some fatigue growth characteristics relevant to each type of cracks are also revealed. It is demonstrated that the fatigue growth analyses of various cracks commonly occurring in bars can reliably be made by using the automated finite element technique proposed.     

Ultrasonic based structural damage detection using combined finite element and model Lamb wave propagation parameters in composite materials

This paper presents a combined finite element and model Lamb waves propagation parameters method as a tool for structural health monitoring in composite materials. Modal analysis allows identifying the mode conversions induced by the defects. A simulation combining a lossless finite element approach and Lamb wave propagation parameter for finding natural frequencies and mode shapes of the structures in undamaged and damaged condition is proposed. This analysis is performed on two carbon-fiber-reinforced plastic bars in both undamaged and damaged state, where the two damaged states are (1) having a cut partway through the bar, perpendicular to the long axis of the bar and (2) having a circular hole. The lamb wave propagation parameters are calibrated using the ultrasonic pulse generator test setup. The natural frequencies for the theoretical, finite element and experimental results are compared and close agreement is found between the frequencies obtained experimentally and computationally.     

A coupled stream function-finite element analysis for wire drawing processes

A numerical approach has been developed based on the stream function technique and the finite element analysis to predict required power and temperature rise in wire drawing processes. An admissible velocity field is first proposed using a stream function and then power consumption in the wire drawing is optimized to achieve sensible and unique deformation geometry. In addition, the finite element method together with axi-symmetric Petrov?CGalerkin algorithm is coupled with the deformation model to assess the temperature distribution in both the deforming wire and the die during the process. The work hardening effects are also considered in the model both in the deformation zone and on the velocity discontinuity surfaces. The model can estimate the effects of various process parameters such as drawing velocity and die geometry. In order to evaluate the results of the model, the predictions are compared with the established models including force equilibrium as a lower-bound approach and an upper-bound solution based on the spherical velocity field.     

复合材料低模量蒙皮结构的损伤特性探讨

低模量蒙皮设计概念是改进复合材料结构损伤容限特性的一种有效的设计方法。它不仅可用于改善结构在拉伸、剪切载荷下的损伤容限特性 ,同时也可提高蒙皮的抗冲击损伤特性。本文主要对用低模量蒙皮设计概念设计的复合材料壁板在拉伸和剪切载荷作用下的损伤扩展特性、剩余强度及其设计方法进行探讨。     

Development of high speed computation algorithm for transient analysis of rectangular plates

A new transient analysis method for a rectangular plate structure comprised of a large number of plate elements was developed in order to significantly reduce computational time and memory. This algorithm was derived from the combination of the transfer technique of the transfer mass coefficient method, the modeling technique of the finite element method, and the numerical integration technique of Newmark’s method. In this paper, the algorithm for the transient analysis of a rectangular plate structure is formulated by the proposed method. In order to verify the computational accuracy and efficiency of the proposed method, the results obtained by the proposed method were compared with those obtained by the finite element method and the finite element-transfer matrix method. The proposed method, the finite element-transfer mass coefficient method, could considerably reduce the computation time without the loss of accuracy, in spite of using small computation memory, by using the transfer rules successively.     

Wrinkling initiation and growth in modified Yoshida buckling test: Finite element analysis and experimental comparison

Wrinkling is one of the major defects in sheet metal products and may also play a significant role in the wear of the tool. The initiation and growth of wrinkles are influenced by many factors such as stress ratios, mechanical properties of the sheet material, geometry of the workpiece, contact condition, etc. It is difficult to analyze the wrinkling initiation and growth considering all the factors because the effects of the factors are very complex and the wrinkling behavior may show a wide scatter of data even for small deviations of factors. In this study, the bifurcation theory is introduced for the finite element analysis of wrinkling initiation and growth. All the above-mentioned factors are conveniently considered by the finite element method. The wrinkling initiation is found by checking the determinant of the stiffness matrix at each iteration and the wrinkling behavior is analyzed by successive iteration with the perturbed guess along the eigenvector. The effect of magnitude of perturbation on the wrinkling behavior can be avoided by the Newton-type iteration method. The finite element formulation is based on the incremental deformation theory and elastic-plastic material modeling. The finite element analysis is carried out using the continuum-based resultant shell elements considering the anisotropy of the sheet metal. For the verification of the analysis, the postbuckling of columns and circular plates are analyzed by finite element analysis using the bifurcation algorithm introduced in the study, and the results are compared with the exact solutions. In order to investigate the effects of geometry and stress ratio on the wrinkling initiation and growth, a modified Yoshida buckling test is proposed as an improved effective buckling test. In the modified Yoshida buckling test, the dimensions of the sheet specimen are varied to change the stress ratio and the degree of constraint. The finite element analysis is carried out for the modified Yoshida buckling test and compared with the experimental results.     

高频谐振载荷作用下Ⅰ型疲劳裂纹尖端力学参数变化规律

为研究高频谐振式疲劳裂纹扩展试验中带有Ⅰ型预制裂纹的紧凑拉伸(CT)试件裂纹尖端力学参数的变化规律,利用动态有限元方法,采用ANSYS和MATLAB软件编写程序,计算了CT试件在高频恒幅正弦交变载荷作用下,在一个应力循环及裂纹扩展到不同长度时裂纹尖端区域的位移、应变场及裂纹尖端的应力强度因子,并分析了其变化规律。在计算裂纹尖端应力强度因子时,首先采用静态有限元方法和理论公式验证了有限元建模和计算的正确性,然后采用动态有限元方法研究了裂纹扩展过程中裂纹尖端应力强度因子的变化规律。最后进行了高频谐振式疲劳裂纹扩展试验,采用动态高精度应变仪测量了裂纹扩展到不同阶段时裂纹尖端点的应变,并对有限元计算结果进行了验证。研究结果表明：在稳态裂纹扩展阶段,高频谐振载荷作用下Ⅰ型疲劳裂纹尖端位移、应变及应力强度因子均为与载荷同一形式的交变量;随着裂纹的扩展,Ⅰ型疲劳裂纹尖端的位移、应变及应力强度因子幅不断增大;静态应力强度因子有限元计算值和理论值的误差为2.51%,裂纹尖端点应变有限元计算结果和试验结果最大误差为2.93% 。     

Vibration response and harmonic wave propagation of ultrasonic arc drivers

A piezoelectric curvilinear arc driver designed for an ultrasonic curvilinear motor is evaluated in this study. A design of piezoelectric curvilinear arc driver is proposed and its governing equations, vibration behaviour and wave propagation are investigated. Then, analysis of forced vibration response or driving characteristics to harmonic excitations in the modal domain is conducted. Finite element modelling and analysis of the arc driver are also discussed. Analytical results of free vibration characteristics are compared favourably with the finite element results. Harmonic analyses of the three finite element models reveal changes of dynamic behaviours of three models and also imply operating frequencies with a significant travelling wave component. Parametric study of mathematical and finite element simulation results suggests that stable travelling waves can be generated to drive a rotor on the proposed curvilinear arc motor system.     

基板预热对激光金属沉积成形过程热应力的影响

为降低沉积过程的热应力,抑制成形过程中裂缝的产生,研究基板预热对激光金属沉积成形(Laser metal deposition shaping, LMDS)过程热应力的影响具有非常重要的意义。根据有限元分析中的“单元生死”思想,利用APDL (ANSYS parametric design language)编程建立多道多层激光金属沉积成形过程的数值模拟模型,深入探讨基板未预热和预热到400 ℃时对成形过程热应力的影响。计算结果表明,基板预热到400 ℃可以显著降低成形过程中试样的热应力变化波动性,试样的Von Mises热应力最大值可降低10%左右,其中x方向热应力最大值可降低8.5%左右,z方向热应力最大值可降低8.1%左右。在与模拟过程相同的条件下,利用自行研制的激光金属沉积成形设备进行了成形试验,成形试验的结果与模拟结果基本吻合。     

板中Lamb波激励响应计算方法

采用基于半解析有限元的激励响应计算方法来求解薄板中Lamb波的激励响应结果,替代常规的实验方法,实现对薄板中Lamb波传播特性的仿真分析。与三维有限元仿真方法相比,采用激励响应计算方法仿真导波在波导介质中的传播过程可以节省计算量,提高效率。通过求解薄板中Lamb波的一般均质方程,基于频谱叠加原理,可以计算薄板中Lamb波激励响应结果。分别提取激励响应计算结果和实验测量数据,通过小波分析计算两种数据的Lamb波群速度。经验证,激励响应仿真计算结果与实验数据有很好的一致性,通过群速度频散曲线对比两种方法获取的数据,模态分析结论一致。激励响应计算方法还可以仿真计算任意截面波导介质中导波的传播过程,具有较好的通用性。     

基于单激发端多接收端压电阵列的板内损伤检测

基于超声导波技术研究了板状结构的损伤检测。理论分析了平板中Lamb波的传播特点;利用有限元耦合场模拟方法模拟了平板中Lamb波的传播及裂纹对Lamb波的影响;根据有限元模拟结果,分析了裂纹尺寸对Lamb波传播的影响;为检测损伤的位置,建立了单激发端多接收端的压电阵列;结合概率成像算法,建立了损伤分布场,该分布场的成像结果直观地反映了损伤的位置和严重程度。     

A study on the impact behavior of adhesively-bonded composite materials

In this paper, a unidirectional carbon-fiber composite is both experimentally and numerically investigated to study the nonlinear material behavior of impacted double cantilever beam (DCB) specimens. For the impact analysis, the load and the displacement applied from pin onto end block as well as the crack energy release rate are measured and compared with the finite element analysis results. The energy release rate is a critical measure of the resistance to crack propagation, which can be estimated by the force and displacement at the crack tip. It is found that the energy release rate measured from impact tests on the specimens is well predicted by the finite element model suggested in this study.