大石峡水电站压力隧洞钢筋混凝土衬砌限裂分析

大石峡引水压力隧洞上平段为钢筋混凝土衬砌,内水压力标准值为0.99 MPa,需要设置一定数量钢筋控制衬砌裂缝分布及其宽度。通过采用基于修正摩尔-库伦准则的弹塑性损伤混凝土有限元法和预设裂缝法,建立非线性数值模拟模型,其中用植入式杆单元模拟钢筋力学效应,用具有抗拉和滑移力学特性的节理单元来模拟衬砌上的预设裂缝。结果表明:仅衬砌内层和内外层均布置钢筋对混凝土衬砌塑性损伤区即开裂出现时机和位置基本没有影响,只对塑性区分布的范围稍有影响;布置钢筋方式对钢筋应力影响较小,钢筋混凝土塑性模型计算的两种情况的钢筋应力相差1.6%,预设裂缝模型计算的相差22.1%。两种数值模型计算出的钢筋混凝土衬砌裂缝宽度均小于规范允许值,满足结构限裂要求。

古运河枢纽明挖段暗涵的防裂限裂温度控制研究

对于与隧洞工程同类型的水工暗涵大体积混凝土超静定结构,要实现完全防裂,要求实施的温度控制相当严格,因此成为工程施工中尚待解决的难题。南水北调中线古运河枢纽明挖段暗涵施工,由于采取了优化混凝土配合比、优化工程原材料、优化施工温控及运用大型混凝土工程一体化连续施工装置,加快混凝土的生产、输送和入仓布料浇筑速度等综合措施,使前述难题获得较好解决。工程完建两年,目前尚未发现裂缝。文章就该暗涵工程的温控状态进行分析和研究,得到令人满意的结果。可供同类工程参考。

An adaptive finite element procedure for crack propagation analysis

This paper presents the adaptive mesh finite element estimation method for analyzing 2D linear elastic fracture problems. The mesh is generated by the advancing front method and the norm stress error is taken as a posteriori error estimator for the h-type adaptive refinement. The stress intensity factors are estimated by a displacement extrapolation technique. The near crack tip displacements used are obtained from specific nodes of natural six-noded quarter-point elements which are generated around the crack tip defined by the user. The crack growth and its direction are determined by the calculated stress intensity factors. The maximum circumference theory is used for the latter. In evaluating the accuracy of the estimated stress intensity factors, four cases are tested consisting of compact tension specimen, three-point bending specimen, central cracked plate and double edge notched plate. These were carried out and compared to the results from other studies. The crack trajectories of these specimen tests are also illustrated.

Prediction of central bursting defects in rod extrusion process with upper bound analysis method

The prediction of central bursting defects in the rod extrusion process through conical dies using the upper bound analysisis investigated. A kinematically admissible velocity field, including the radial and angular velocity components, is proposed. A newcriterion is presented to predict the occurrence of the central bursting defects. Parameter bobt, which represents the risk probability ofcracking, is proposed. It is calculated using the shape of the boundary at the entrance by minimizing the total power dissipationduring the extrusion process. When bobt is equal to or greater than bcr, central bursting occurs. Furthermore, the quantitativerelationships between central bursting defects and process parameters （semi die angle, reduction in area and frictional factor） arestudied. The results show that the central bursting defects are affected primarily by the reduction in area and the friction factor. Thepresented criterion is verified by comparing with the FEM simulation data and the results of the published paper.

Distribution of stress and strain between adjacent particles in particulate reinforced metal matrix composites

The distribution of stress and strain between adjacent particles in particulate reinforced metal matrix composites wasinvestigated using cohesive zone models. It is found that the strain of the composite is concentrated in the matrix, and there is aregion with higher strain along the loading path, which can promote the formation of a void near the particles pole. The stress andstrain in matrix near the particles gradually decrease with the increase of the distance between particles. And it is calculated that thereis a critical distance within which the stress and strain fields of the neighboring particles can influence with each other. This criticaldistance increases with the increase of particle size. It is also found that the angle between the tensile direction and the center line ofparticles plays an important role in the stress and strain distribution. The model with the angle of 0° has the greatest influence on thedistribution of stress and strain in the matrix, while the model with the angle of 45° has the least influence on the distribution of stressand strain in the matrix.

H_2S Stress Corrosion Tests of Welded Joint for X65 Pipeline Steel and Finite Element Numerical Analysis of Crack Tip

The microstructure of welded joint is surveyed and the mechanical properties of X65 pipeline steel are studied in this paper, which provides experimental basis of performance effect on stress corrosion. H 2S stress corrosion cracking (SCC) tests on the steel are carried out in the environment based on NACE TM- 01-77 solution. The threshold stress intensity factor and crack propagation velocity for base metal and HAZ are obtained. The susceptibility of welded joint for X65 pipeline steel to H 2S stress corrosion cracking is investigated. The programming package ANSYS of finite element model (FEM) is used to perform the three-dimensional elastic-plastic finite element analysis of WOL specimens. Stress field and concentration of hydrogen distribution property of the crack tip are obtained.

Experimental and finite element analysis for fracture of coating layer of galvannealed steel sheet

Mechanical properties of galvannealed (GA) steel sheet used for automotive exposed panel and predicted failure phenomenon of its coating layer were evaluated using finite element method. V-bending test was performed to understand better the fracture of coating layer of GA steel sheet during plastic deformation. Yield strength of the coating layer was calculated by using a relative difference between hardness of coating layer measured from the nano-indentation test and that of substrate. To measure shearing strength at the interface between substrate and coating layer, shearing test with two specimens attached by an adhesive was carried out. Using the mechanical properties measured, a series of finite element analyses coupled with a failure model was performed. Results reveal that the fracture of coating layer occurs in an irregular manner at the region where compressive deformation is dominant. Meanwhile, a series of vertical cracks perpendicular to material surface are observed at the tensile stressed-region. It is found that 0.26-0.28 of local equivalent plastic strain exists at the coating and substrate at the beginning of failure. The fracture of coating layer depends on ductility of the coating layer considerably as well.

Simulation of Corrosion Crack in Pretensioned Prestressed Concrete Members

Corrosion of PC tendon embedded in concrete deteriorates performance of PC （prestressed concrete） members. Corrosion crack observed in pretensioned PC members is different from that in reinforced concrete and post tensioned PC members. The purpose of this study is to clarify mechanisms of corrosion crack in pretensioned PC members. Contribution of compressive stress induced by prestressing force and arrangement of PC tendons on corrosion crack in pretensioned PC members is investigated through experiment and numerical analysis. Three dimensional finite element analyses were carried out to simulate corrosion crack in PC specimen. Based on the analytical and experimental results, the influencing factor of corrosion crack in pretensioned PC members is discussed. It was concluded that the contribution of compressive stress in concrete by prestressing force on corrosion crack is not significant and that the dominant reason for wider corrosion cracks in pretensioned PC members is the arrangement of PC tendons with small intervals.

An Analysis on the Transmission Shaft Fracture of Screw Drill

The transmission shaft of the underground screw drill fractured when milling-shoe and grinding bridge plug was applied to the coiled tubes in a horizontal well of Sichuan province, but the position of this transmission shaft fracture did not occur at the minor-diameter retracting position and reducing position that were easy to fracture. An analysis of the transverse planes of the fracture found that the cause of the transmission shaft fracture surface was that the defects of the initial surface were propagating to cracks and gave rise to the fracture under torque load. To specifically know well the strength of the transmission shaft under damaging fracture, a statics analysis was conducted on the transmission shaft through ANSYS finite element simulation software, and the finite element models under no cracks, different-depth circumferential cracks, and similar situation of the transverse planes of fracture were established respectively. An analysis of the crack-free finite element model found that the fracture of the transmission shaft was really not caused by the self-structure of the transmission shaft; an analysis of circumferential crack finite element model found that strong stress concentration would appear in the tip of cracks, and the value of the stress would increase along with the increase of the circumferential crack depth, the stress of the entire crack top tended to fluctuate like waves, and also the strength of the transmission shaft was greatly impacted by the presence of cracks; an analysis of the similar crack finite element model found that stress concentration would appear in the tip of cracks, and the initial cracks always started to propagate from the tip of the external surface of the transmission shaft and would propagate inward until the propagating areas of two cracks overlapped, and finally reached the position of transient interruption, and then the transmission shaft fractured completely and the fracture strength was onlv 1/5 of that under no cracks.

Finite element analysis for blasting law of jointed rock mass

On the basis of Mises strength theory,rock models are built including vertical,horizontal and diagonal joints to simulate jointed rock mass under blasting load by using FEM. The dynamic procedures of jointed rock mass under blasting are quantified and the effective stress-time curves of typical elements are compared to analyze the barrier of joints to the stress wave. The blasting law was studied according to the process of computer simulation and the effect of blasting,and some suggestions were given for solving the problems of overbreak and underbreak.

Computational Fracture Analysis of an AFM-Specimen under Mixed Mode Loading Conditions

Fracture processes in ship-building structures are in many cases of a 3-D character. A finite element （FE） model of an all fracture mode （AFM） specimen was built for the study of 3-D mixed mode crack fracture behavior including modes Ⅰ,Ⅱ, and Ⅲ. The stress intensity factors （SIFs） were calculated by the modified virtual crack closure integral （MVCCI） method, and the crack initiation angle assessment was based on a recently developed 3-D fracture criterion--the Richard criterion. It was shown that the FE model of the AFM-specimen is applicable for investigations under general mixed mode loading conditions, and the computational results of crack initiation angles are in agreement with some available experimental findings. Thus, the applicability of the FE model of the AFM-specimen for mixed mode loading conditions and the validity of the Richard criterion can be demonstrated.

Fracture analysis of rock mass based on 3-D nonlinear Finite Element Method

Traditional fracture analysis is based on fracture mechanics and damage mechanics. They focus on the propagation of the fracture. However, their propagation criterions are not easily applied in practice and the current analysis is limited in planar problem. This paper presents a new theory that the occurrence of the unbalanced force (derived from the Deformation Reinforcement Theory) could be the criterion of the initiation of the fracture, and the distribution area and propagation of the unbalanced force could be the indication of the fracture propagation direction. By aggregate analysis with Stress Intensity Factor (SIF) criterion, the unbalanced force actually is the opposite external load that is the SIF difference incurred between the external loads and permitted by the structure. Numerical simulation and physical experiments on pre-fracture cuboid rock specimens proved that the occurrence of the unbalanced force could be the initiation of the fracture. Mesh size dependence was also considered by analysis of different mesh size finite element gravity dam models. Furthermore, the theory was applied to the feasibility analysis of the Baihetan arch dam together with physical experiments in order to evaluate the fracture propagation of dam heel. The results show that it is an effective way to use unbalanced force to analyze the fracture initiation and propagation when performing 3-dimensional nonlinear FEM calculation.

Nonlinear finite element analysis of crack growth at the interface of rubber-like bimaterials

This paper presents the characteristics of the crack growth at the interface of rubber-rubber and rubber-steel bimaterials undertensile deformation using the non-linear finite element method. By using the commercial finite element software ABAQUS,the J integral calculations are carried out for the initial interface crack in the interfaces in-between two Neo-Hookean materials,two Mooney-Rivlin materials, Neo-Hookean and Mooney-Rivlin rubbers, Neo-Hookean and Polynomial, Mooney-Rivlin andPolynomial, and the Mooney-Rivlin and steel bi-materials. The computational results of the maximum J integral directionaround the crack tip illustrate the possible direction of crack growth initiation. Furthermore, it is found that the crack bends tothe softer rubber material at a certain angle with the initial crack direction if the crack depth is relatively small. For the crackwith a larger depth, the crack propagates to grow along the interface in-between the bimaterials.

Stability analysis of the temperature cracks in Xiaowan arch dam

In the Xiaowan arch dam there are massive temperature cracks nearly parallel to the dam axis. Obviously, whether the cracks may spread or not during the water storage process is one of the crucial factors for the safety of a dam. In this paper, a new type of crack element, in which the contact component is implicitly included into the concrete component, is proposed to simulate the effects of the existing cracks. The crack element is proved by numerical example to share the merits of both conventional contact elements and joint elements. With a finite element model of the cracked arch dam together with its rock foundation established, the transient displacement and stress fields of the dam are obtained. The complicated rock foundation, the construction process of the arch dam, the massive cracks, the transient temperature field, as well as the water storage process have been taken into consideration in the simulation. In addition to the global model, several sub-models for typical crack tips are also generated with finer elements placed around the tips. Thus, more accurate displacement and stress distribution are obtained by simultaneous sub-model simulation. Based on the calculation of stress intensity factor for crack tips by extension method, the temperature cracks in the Xiaowan arch dam are finally proved to be stable.

Determination of T-stress using finite element analysis

A method that uses finite element analysis to determine the non-singular stress （T-stress） at a crack tip is proposed in this study. T-stress includes two components： the Tx-stress parallel to the tangent of the crack at its tip and the Ty-stress perpendicular to this tangent. The effects of contact and friction on both the Tx- and Ty-stresses on the crack flanks are considered in the method. Because the method uses a single standard elastic finite element analysis derived directly from the equation of the stress fields around the crack tip and does not require any assumptions or simplification, it can be used to determine the T-stress for any given geometry and loading condition. Theoretical results are used to calibrate the results, which exhibited good agreement and to discuss the T-stress computational methodology. Furthermore, the Tx- and Ty-stresses in center-cracked Brazilian disc （CCBD） specimens subjected to diametrical or partially distributed compression were numerically computed, and the effects of contact and friction on the Tx- and Ty-stresses are discussed.

Variation of stress intensity factors for several interface crack problems under arbitrary material combinations

In this paper the variation of interface stress intensity factors is considered under arbitrary material combinations for several problems such as central internal interface cracks, periodic interface cracks, and edge interface cracks. The stress intensity factors for all these problems have been determined in terms of Dundurs＇ parameters a, b by the application of FEM in the authors＇ previous papers. In this paper, the variations of the stress intensity factors for all these problems are discussed under arbitrary a, b, then the maximum and minimum values of dimensionless stress intensity factors FI,max, FII,max, FI,min, FII,max are indicated on the map of a, b. As an example, it is found that FI,max always appears when （a, b）= （0.2,0.3）, and FI,min always appears when （a, b）= （1,0）, for the central internal interface cracks independent of crack length.

Stress intensity factors under combined bending and torsion moments

This paper discusses stress intensity factor （SIF） calculations for surface cracks in round bars subjected to combined torsion and bending loadings. Different crack aspect ratios, a/b, ranging from 0.0 to 1.2 and relative crack depths, a/D, ranging from 0.1 to 0.6 were considered. Since the loading was non-symmetrical for torsion loadings, a whole finite element model was constructed. Then, the individual and combined bending and torsion loadings were remotely applied to the model. The equivalent SIF method, F＊EQ, was then used explicitly to combine the individual SIFs from the bending and torsion loadings. A comparison was then carried out with the combined SIF, F＊FE, obtained using the finite element analysis （FEA） under similar loadings. It was found that the equivalent SIF method successfully predicted the combined SIF for Mode I. However, discrepancies between the results determined from the different approaches occurred when Fm was involved. It was also noted that the predicted F＊ EE using FEA was higher than the F＊EQ predicted through the equivalent SIF method due to the difference in crack face interactions.

A 3D nonlinear finite element method for the dynamic analysis of rotating rotor with a transverse crack

A general and efficient method is presented in this paper for studying the effects of unbalance on the breathing mechanism of crack.Based on 3D finite element models combined with a nonlinear contact approach for crack modeling, the method is free from theassumption of weight-dominance and can be used to gain deep insights into the breathing mechanism of crack. In order to greatlyreduce the computational time, a complex free-interface component mode synthesis (CMS) method is employed to reduce theorder of the model. Based on the proposed method, the effects of unbalance on the breathing mechanism of crack are discussed.Numerical results show that the unbalance can lead to significant changes in the breathing of crack, even when the unbalance force is about an order of magnitude smaller than the self-weight. Moreover, the level and orientation of the unbalance have also remarkable effects on the breathing behaviors of crack. Besides, a new universal non-steady breathing phenomenon of crack is firstly found in this paper, which denotes that the breathing speed of a crack is fluctuated over one revolution when there exists residual unbalance in the cracked rotor.