非连续变形接触问题的分区有限元-块体界面元动力混合算法

本文针对动力情况下不连续变形问题,提出一种新的接触面动力迭代求解方法,即分区有限元与块体界面元动力混合算法。在接触问题基本概念的基础上描述了该混合算法的力学模型。由接触力与接触点对总位移的非线性关系,推导了包含接触力、块体形心刚体加速度和接触节点总位移在内的非线性方程;建立了以界面接触力和块体形心刚体加速度为混合变量与界面结点总位移进行迭代求解的动力界面迭代方程。由此,能够精确求解至界面上接触点对全部失效的状态,从而实现连续—非连续过程的模拟。该混合解法仅将非线性迭代收缩到接触界面上进行,从而大幅提高局部位移非连续问题的计算效率。通过典型算例验证了该混合算法的正确性和有效性。

IBE-PFE方法在某水电站泄洪洞出口边坡稳定分析中的应用

鉴于传统方法在含复杂结构面岩质边坡稳定分析中的局限性,引进分区有限元与块体界面元混合算法,验算边坡的稳定性。针对某水电站泄洪洞出口边坡地形地质条件和结构面分布特点,建立数值模型。采用分区有限元与块体界面元混合算法,并结合强度折减法,逐步降低不连续面抗剪强度,计算得出块体刚体位移和接触点对接触状态。根据块体刚体位移发生突变并结合滑动通道失效点对完全连通,确定了边坡的安全系数和破坏方式。结果表明:边坡稳定安全系数为1.80;边坡最可能以坡顶预设缝、坡脚预设缝以及f_(555)结构面为滑动通道发生失稳破坏。该方法具有准确确定边坡危险滑动面和安全系数的优点。

Parallel hydrodynamic finite element model with an N-Best refining partition scheme

We enhance a robust parallel finite element model for coasts and estuaries cases with the use of N-Best refinement algorithms,in multilevel partitioning scheme.Graph partitioning is an important step to construct the parallel model,in which computation speed is a big concern.The partitioning strategy includes the division of the research domain into several semi-equal-sized sub-domains,minimizing the sum weight of edges between different sub-domains.Multilevel schemes for graph partitioning are divided into three phases:coarsening,partitioning,and uncoarsening.In the uncoarsening phase,many refinement algorithms have been proposed previously,such as KL,Greedy,and Boundary refinements.In this study,we propose an N-Best refinement algorithm and show its advantages in our case study of Xiamen Bay.Compared with original partitioning algorithm in previous models,the N-Best algorithm can speed up the computation by 1.9 times,and the simulation results are in a good match with the in-situ data.

Finite element analysis of temperature and stress fields during selective laser melting process of Al−Mg−Sc−Zr alloy

A 3D finite element model was established to investigate the temperature and stress fields during the selective laser melting process of Al−Mg−Sc−Zr alloy.By considering the powder−solid transformation,temperaturedependent thermal properties,latent heat of phase transformations and molten pool convection,the effects of laser power,point distance and hatch spacing on the temperature distribution,molten pool dimensions and residual stress distribution were investigated.Then,the effects of laser power,point distance and hatch spacing on the microstructure,density and hardness of the alloy were studied by the experimental method.The results show that the molten pool size gradually increases as the laser power increases and the point distance and hatch spacing decrease.The residual stress mainly concentrates in the middle of the first scanning track and the beginning and end of each scanning track.Experimental results demonstrate the accuracy of the model.The density of the samples tends to increase and then decrease with increasing laser power and decreasing point distance and hatch spacing.The optimum process parameters are laser power of 325−375 W,point distance of 80−100μm and hatch spacing of 80μm.

Structural behaviors of steel roof truss exposed to pool fire

Experimental research was conducted to study the structural behaviors of a steel roof truss model without fire-proof coating under pool fire conditions. The data of temperature distribution and displacements of typical members were obtained. It is found that the temperature distribution of environment inside the structure, which is found to be in accordance with the multi-zone model with height, has a decisive effect on the tempera^tre evolution of steel members. Besides, it can also be observed that due to the restriction and coordination among the truss members in the localized fire, the maximum relative deflection, which occurs at the mid-span of the top chord, is relatively slight and has not exceeded 1 mm under experimental conditions. On the other hand, the column experiences a notable thermal expansion during the test. Then, a finite element model is presented and validated by the test results.

Forming condition of transient saturated zone and its distribution in residual slope under rainfall conditions

Rainfall, as one of the most significant factors triggering the residual soil slope failure, leads to not only the reduction of soil shear strength, but also the increase of soil weight and the decrease of matric suction as well. All these modifications in soil properties have important influence on the slope stability. The water infiltration and redistribution inside the slope are the preconditions of the slope stability under rainfall conditions. Based on the numerical simulation via finite element method, the water infiltration process under rainfall conditions was studied in the present work. The emphases are the formation, distribution and dissipation of transient saturated zone. As for the calculation parameters, the SWCC and the saturated permeability have been determined by pressure plate test and variable head test respectively. The entire process(formation, development, dissipation) of the transient saturated zone was studied in detail. The variations of volumetric water content, matric suction and hydraulic gradient inside the slope, and the eventually raise of groundwater table were characterized and discussed, too. The results show that the major cause of the formation of transient saturated zone is ascribed to the fact that the exudation velocity of rainwater on the wetting front is less than the infiltration velocity of rainfall; as a result, the water content of the soil increases. On the other hand, the formation and extension of transient saturated zone have a close relationship with rainfall intensity and duration. The results can help the geotechnical engineers for the deeper understanding of the failure of residual slope under rainfall condition. It is also suggested that the proper drainage system in the slope may be the cost-effective slope failure mitigation method.

Triangular element partition method with consideration of crack tip

In fracture simulation,how to model the pre-existing cracks and simulate their propagation without remeshing is an important topic.The newly developed triangular element partition method(TEPM)provides an efficient approach to this problem.It firstly meshes the cracked body regardless of the geometry integrity of the interesting object with triangular elements.After the meshing procedure is completed,some elements are intersected by cracks.For the element intersected by a crack,the TEPM takes the element partition technique to incorporate the discontinuity into the numerical model without any interpolation enrichment.By this approach,the TEPM can simulate fracture without mesh modification.In the TEPM,all the cracked elements are treated as the usual partitioned elements in which the crack runs through.The virtual node pairs(the intersection points of crack faces and elements)at the opposite faces of the crack move independently.Their displacements are respectively determined by their neighbor real nodes(nodes formatted in the original mesh scheme)at the same side of the crack.However,among these cracked elements,the element containing a crack tip,referred to as the crack tip element thereafter,behaves differently from those cut through by the crack.Its influence on the singular field at the vicinity of the fracture tip becomes increasingly significant with the element size increasing.In the crack tip element,the virtual node pair at the crack tip move consistently before fracture occurs while the virtual node pair separate and each virtual node moves independently after the fracture propagates.Accordingly,the crack tip element is automatically transformed into the usual partitioned element.In the present paper,the crack tip element is introduced into the TEPM to account for the effect of the crack tip.Validation examples indicate that the present method is almost free from the element size effect.It can reach the same precision as the conventional finite element method under the same meshing scheme.But the TEPM is much more efficient and convenient than the conventional finite element method because the TEPM avoids the troubles that the conventional finite element method suffers,e.g.,the meshing problem of cracked body,modification of mesh scheme,etc.Though the extended finite element method can also avoid these troubles,it introduces extra degrees of freedom due to node interpolation enrichment.Due to the simplicity of the present TEPM,it is believed that its perspective should be highly inspiring.

Contact stress non-uniformity of wafer surface for multi-zone chemical mechanical polishing process

A finite element analysis(FEA)model is developed for the chemical-mechanical polishing(CMP)process on the basis of a 12-in five-zone polishing head.The proposed FEA model shows that the contact stress non-uniformity is less dependent on the material property of the membrane and the geometry of the retaining ring.The larger the elastic modulus of the pad,the larger contact stress non-uniformity of the wafer.The applied loads on retaining ring and zone of the polishing head significantly affect the contact stress distribution.The stress adjustment ability of a zone depends on its position.In particular,the inner-side zone has a high stress adjustment ability,whereas the outer-side zone has a low stress adjustment ability.The predicted results by the model are shown to be consistent with the experimental data.Analysis results have revealed some insights regarding the performance of the multi-zone CMP.