基于有限断裂法和比例边界有限元法的裂纹扩展模拟

基于有限断裂法和比例边界有限元法提出了一种裂缝开裂过程模拟的数值模型。采用基于有限断裂法的混合断裂准则作为起裂及扩展的判断标准,当最大环向应力和能量释放率同时达到其临界值时,裂缝扩展。结合多边形比例边界有限元法,可以半解析地求解裂尖区域附近的应力场和位移场,在裂尖附近无需富集即可获得高精度的解。计算能量释放率时,只需将裂尖多边形内的裂尖位置局部调整,无需改变整体网格的分布,网格重剖分的工作量降至最少。裂缝扩展步长通过混合断裂准则确定,避免了人为假设的随意性,并可以实现裂缝变步长扩展的模拟,更符合实际情况。通过对四点剪切梁的复合型裂缝扩展过程的模拟,对本文模型进行了验证,并应用于重力坝模型的裂缝扩展模拟,计算结果表明,本文提出的模型简单易行且精度较高。

Numerical simulation and optimization of process parameters in fineblanking process

Fineblanking process is a typical large localized plastic deformation process. Based on its forming characteristics, a numerical model is established and an elasto-plastic simulation is performed using the finite element method （FEM）. The re-meshing method is used when the severe element distortion occurs to facilitate further computation and avoid divergence. The McClintock fracture criterion is adopted to predict and determine the time and site of crack initiation and propagation. Based on this numerical model, the distribution and developing trend of the stress and strain in the shearing zone are predicted. Furthermore, the influence of several process parameters, such as punch-die clearance, edge radius of punch and die, V-ring force, counter force, etc., on the blanked quality is analyzed. The discipline is in accordance with the actual manufacture situation, which can be a guidance to optimization of process parameters.

Fatigue fracture analysis of gear teeth using XFEM

Gear teeth in gear transmission systems suffer seriously from fatigue failure during service. In this work, a 2 D double-tooth model was constructed with periodic boundary conditions. The fatigue fracture behavior of gear teeth was analyzed using the extended finite element method(XFEM), with emphases on the impacts of initial crack geometries and cyclic load factors. The results suggested that the shortest fatigue life is expected for 0° orientation cracks initiating at the maximum principal stress. Cracks that initiate closer to the bottom land of gear tooth are relatively safe. Moreover, to evaluate the fatigue load conditions, load ratio, load range, and mean load should be all taken into considerations. Further XFEM simulation for material selection was performed to guide the gear design. Among various material parameters, the material constant C and tensile strength are the most significant ones in determining the fatigue life.

Stress based forming limit diagram for formability characterization of 6061 aluminum

Two numerical criteria of forming limit diagram（FLD） criterion and ductile fracture criterion（DFC） are presented for FLD prediction of 6061 aluminum. The numerical results are compared with the experimental FLD and also punch＇s load-displacement curve of experimental samples. Experimental FLD of this study is calculated using hemispherical punch test of Hecker. Experimental FLD is converted to FLSD and imported to the Abaqus software to predict necking of samples. Numerical results for FLSD prediction were compared with experimental FLSD. Results show that ductile fracture criterion has higher accuracy for FLD and FLSD prediction of 6061 aluminum. Comparison of numerical and experimental results for force-displacement curve of punch shows that numerical results have a good agreement with experiment.

Fracture properties of epoxy asphalt mixture based on extended finite element method

Crack is found to be a major distress that affects the performance of the epoxy asphalt pavement.An extended finite element method was proposed for investigating the fracture properties of the epoxy asphalt mixture.Firstly,the single-edge notched beam test was used to analyze the temperature effect and calculate the material parameters.Then,the mechanical responses were studied using numerical analysis.It is concluded that 5℃ can be selected as the critical temperature that affects the fracture properties,and numerical simulations indicate that crack propagation is found to significantly affect the stress state of the epoxy asphalt mixture.The maximum principal stress at the crack surface exhibits different trends at various temperatures.Numerical solution of stress intensity factor can well meet the theoretical solution,especially when the temperature is lower than 5℃.

A New 3D Meso-mechanical Modeling Method of Coral Aggregate Concrete Considering Interface Characteristics

On the basis of the three-dimensional(3D)random aggregate&mortar two-phase mesoscale finite element model,C++programming was used to identify the node position information of the interface between the aggregate and mortar elements.The nodes were discretized at this position and the zero-thickness cohesive elements were inserted.After that,the crack energy release rate fracture criterion based on the fracture mechanics theory was assigned to the failure criterion of the interface transition zone(ITZ)elements.Finally,the three-phase mesomechanical model based on the combined finite discrete element method(FDEM)was constructed.Based on this model,the meso-crack extension and macro-mechanical behaviour of coral aggregate concrete(CAC)under uniaxial compression were successfully simulated.The results demonstrated that the meso-mechanical model based on FDEM has excellent applicability to simulate the compressive properties of CAC.

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.

An Improved Ductile Fracture Criterion for Fine-blanking Process

In order to accurately simulate the fine-blanking process, a suitable ductile fracture is significant.So an evaluation strategy based on experimental and corresponding simulation results of tensile, compression, torsion and fine-blanking test is designed to evaluate five typical ductile fracture criteria, which are widely-used in metal forming process.The stress triaxiality and ductile damage of each test specimen are analyzed.The results show that none of these five criteria is sufficient for all tests.Furthermore, an improved fracture criterion based on Rice and Tracey model, taking the influence of both volume change and shape change of voids into account, is proposed.The characterization of this model for fine-blanking process is easily done by the tensile test and the prediction result shows good.

An efficient hybrid model for fractured reservoirs

In this study, an efficient hybrid model is proposed to simulate the fluid flow in the reservoirs with multi-scale fractures, which cannot be easily modeled by neither the continuum models nor the discrete fracture models. In the proposed method, the small fractures are modeled by using an improved Multiple Sub-Region method, which can capture the strongly anisotropy of fracture elements and the effects of border region on the transmissibility and provide more accurate results, on the other hand, the large fractures are modeled explicitly as major fluid conduits by the Embedded Discrete Fracture Model. Then, an efficient numerical algorithm based on the Mimetic Finite Difference method is developed to solve the hybrid method. At the end, several numerical examples are carried out to verify the accuracy and applicability of the proposed numerical model.