P_1-nonconforming triangular finite element method for elliptic and parabolic interface problems

The lowest order Pl-nonconforming triangular finite element method （FEM） for elliptic and parabolic interface problems is investigated. Under some reasonable regularity assumptions on the exact solutions, the optimal order error estimates are obtained in the broken energy norm. Finally, some numerical results are provided to verify the theoretical analysis.

3D Finite Element Simulation of Tunnel Boring Machine Construction Processes in Deep Water Conveyance Tunnel

Applying stiffness migration method,a 3D finite element mechanical model is established to simulate the excavation and advance processes.By using 3D nonlinear finite element method,the tunnel boring machine(TBM) excavation process is dynamically simulated to analyze the stress and strain field status of surrounding rock and segment.The maximum tensile stress of segment ring caused by tunnel construction mainly lies in arch bottom and presents zonal distribution.The stress increases slightly and limitedly in the course of excavation.The maximum and minimum displacements of segment,manifesting as zonal distribution,distribute in arch bottom and vault respectively.The displacements slightly increase with the advance of TBM and gradually tend to stability.

Propagations of Rayleigh and Love waves in ZnO films/glass substrates analyzed by three-dimensional finite element method

Propagation characteristics of surface acoustic waves（SAWs） in ZnO films/glass substrates are theoretically investigated by the three-dimensional（3D） finite element method. At first, for（11ˉ20） ZnO films/glass substrates, the simulation results confirm that the Rayleigh waves along the [0001] direction and Love waves along the [1ˉ100] direction are successfully excited in the multilayered structures. Next, the crystal orientations of the ZnO films are rotated, and the influences of ZnO films with different crystal orientations on SAW characterizations, including the phase velocity, electromechanical coupling coefficient, and temperature coefficient of frequency, are investigated. The results show that at appropriate h/λ, Rayleigh wave has a maximum k^2 of 2.4% in（90°, 56.5°, 0°） ZnO film/glass substrate structure; Love wave has a maximum k^2 of 3.81% in（56°, 90°, 0°） ZnO film/glass substrate structure. Meantime, for Rayleigh wave and Love wave devices, zero temperature coefficient of frequency（TCF） can be achieved at appropriate ratio of film thickness to SAW wavelength. These results show that SAW devices with higher k^2 or lower TCF can be fabricated by flexibly selecting the crystal orientations of ZnO films on glass substrates.

THE h-p VERSION OF THE INFINITE ELEMENT METHOD AND ITS APPLICATION

The paper investigates the h-p version of the infinite element method. An exponential conver gence can be abtained by a small amount of computing work.

Maximum Norm Estimates for Finite Volume Element Method for Non-selfadjoint and Indefinite Elliptic Problems

In this paper, we establish the maximum norm estimates of the solutions of the finite volume element method （FVE） based on the P1 conforming element for the non-selfadjoint and indefinite elliptic problems.

Uniform Convergence for Finite Volume Element Method for Non-selfadjoint and Indefinite Elliptic Problems

In this paper, we prove the existence, uniqueness and uniform convergence of the solution of finite volume element method based on the P1 conforming element for non-selfadjoint and indefinite elliptic problems under minimal elliptic regularity assumption.

Grain refinement of Mg-Al alloys by optimization of process parameters based on three-dimensional finite element modeling of roll casting

To study the influence of roll casting process parameters on temperature and thermal-stress fields for the AZ31 magnesium alloy sheets,three-dimensional geometric and 3D finite element models for roll casting were established based on the symmetry of roll casting by ANSYS software.Meshing method and smart-sizing algorithm were used to divide finite element mesh in ANSYS software.A series of researches on the temperature and stress distributions during solidification process with different process parameters were done by 3D finite element method.The temperatures of both the liquid-solid two-phase zone and liquid phase zone were elevated with increasing pouring temperature.With the heat transfer coefficient increasing,the two-phase region for liquid-solid becomes smaller.With the pouring temperature increasing and the increase of casting speed,the length of two-phase zone rises.The optimized of process parameters（casting speed 2 m/min,pouring temperature 640 ℃ and heat transfer coefficient 15 kW/（m2·℃） with the water pouring at roller exit was used to produce magnesium alloy AZ31 sheet,and equiaxed grains with the average grain size of 50 μm were achieved after roll casting.The simulation results give better understanding of the temperature variation in phase transformation zone and the formation mechanism of hot cracks in plates during roll casting and help to design the optimized process parameters of roll casting for Mg alloy.

Finite Element Simulation of Flexible Roll Forming with Supplemented Material Data and the Experimental Verification

Flexible roll forming is a promising manufacturing method for the production of variable cross section products. Considering the large plastic strain in this forming process which is much larger than that of uniform deformation phase of uniaxial tensile test, the widely adopted method of simulating the forming processes with non-supplemented material data from uniaxial tensile test will certainly lead to large error. To reduce this error, the material data is supplemented based on three constitutive models. Then a finite element model of a six passes flexible roll forming process is established based on the supplemented material data and the original material data from the uniaxial tensile test. The flexible roll forming experiment of a B pillar reinforcing plate is carried out to verify the proposed method. Final cross section shapes of the experimental and the simulated results are compared. It is shown that the simulation calculated with supplemented material data based on Swift model agrees well with the experimental results, while the simulation based on original material data could not predict the actual deformation accurately. The results indicate that this material supplement method is reliable and indispensible, and the simulation model can well reflect the real metal forming process. Detailed analysis of the distribution and history of plastic strain at different positions are performed. A new material data supplement method is proposed to tackle the problem which is ignored in other roll forming simulations, and thus the forming process simulation accuracy can be greatly improved.

OPTIMIZATION OF SHEET METAL FORMING PROCESSES USING FINITE ELEMENT SIMULATIONS

The paper focuses on the combination of the Finite Element simulation and optimization to improve process or product quality. Three different examples to illustrate the developed genetic approach are given. In all three examples is-DYNA3D is used to simulate the process and a general aptimiza- tion sensitivity based strategy is utilized to improve the design. The included examples are: 1) stretch bending of tubes, 2) bulging of tubes, and finally 3) hydromechanical deep drawing. these examples clearly illustrate the potential of systematic optimization in the area of metal processing.

STUDY ON THREE-DIMENSIONAL FINITE BODIES CONTAINING CRACKS USING THE FINITE ELEMENT METHOD OF LINES

The three-dimensional finite element method of lines is presented, and the basic processing description of 3D FEMOL in cracking questions is given in detail. Applications to 3D bodies with cracks indicate that good accuracy can be obtained with relatively coarse girds. In particular, application to the tension specimen shows very good agreement with the evaluation of stress intensity factors, which is better than the results of other methods. This implies a considerable potential for using this method in the 3D analysis of finite geometry solids and suggests a possible extension of this technique to nonlinear material behavior.

Numerical Study of the Biomechanical Behavior of a 3D Printed Polymer Esophageal Stent in the Esophagus by BP Neural Network Algorithm

Esophageal disease is a common disorder of the digestive system that can severely affect the quality of life andprognosis of patients. Esophageal stenting is an effective treatment that has been widely used in clinical practice.However, esophageal stents of different types and parameters have varying adaptability and effectiveness forpatients, and they need to be individually selected according to the patient’s specific situation. The purposeof this study was to provide a reference for clinical doctors to choose suitable esophageal stents. We used 3Dprinting technology to fabricate esophageal stents with different ratios of thermoplastic polyurethane (TPU)/(Poly-ε-caprolactone) PCL polymer, and established an artificial neural network model that could predict the radial forceof esophageal stents based on the content of TPU, PCL and print parameter. We selected three optimal ratios formechanical performance tests and evaluated the biomechanical effects of different ratios of stents on esophagealimplantation, swallowing, and stent migration processes through finite element numerical simulation and in vitrosimulation tests. The results showed that different ratios of polymer stents had different mechanical properties,affecting the effectiveness of stent expansion treatment and the possibility of postoperative complications of stentimplantation.

PBX-3造型粉颗粒压制成型过程仿真分析

为了探究3号高聚物粘结炸药(PBX-3)压制成型过程中的颗粒演变行为和密度演化规律,采用多颗粒有限元法(MPFEM)结合Johnson-Cook(J-C)模型对PBX-3粉末在压制过程中的致密化行为进行了模拟仿真;将仿真结果与Heckel方程拟合,分析了压制过程中粉体的温度与压紧力大小对粉体颗粒致密化质量的影响,并讨论了粉末颗粒中的应力和应变分布。结果表明:在PBX-3致密化过程的滑移阶段、变形阶段以及压实阶段中,升高压制温度能有效减小压制过程中所需的压紧力,并提高药柱成型密度。

An explicit method for numerical simulation of wave equations: 3D wave motion

In this paper, an explicit method is generalized from 1D and 2D models to a 3D model for numerical simulation of wave motion, and the corresponding recursion formulas are developed for 3D irregular grids. For uniform cubic grids, the approach used to establish stable formulas with 2M-order accuracy is discussed in detail, with M being a positive integer, and is illustrated by establishing second order （M=1） recursion formulas. The theoretical results presented in this paper are demonstrated through numerical testing.

A NEW METHOD FOR SOLUTION OF 3D ELASTIC-PLASTIC FRICTIONAL CONTACT PROBLEMS

The solution of 3 D elastic-plastic frictional contact problems belongs to the un specified boundary problems where the interaction between two kinds of nonlinearities should occur. Considering the difficulties for the solution of 3 D frictional contact problems, the key part is the determination of the tangential slip states at the contact points, and a great amount of computing work is needed for a high accuracy result. A new method based on a combination of programming and iteration methods, which are respectively known as two main kinds of methods for contact analysis, was put forward to deal with 3 D elastic-plastic contact problems. Numerical results demonstrate the efficiency of the algorithm illustrated here.

基于SP_(3)方法的动力堆中子噪声分析程序研究

本文基于SP_(3)方程和有限元方法研制了中子噪声分析程序CORCA-NOISE。以国际原子能机构发布的三维压水堆基准例题对程序进行了测试,数值结果表明该程序对通量、keff等稳态参数的计算具有很高的精度。在此基础上基于IAEA 2D基准例题对中子噪声分析进行了初步的计算,给出了在给定输入扰动条件下的中子噪声振幅、相位的空间分布。本程序的研发为下一步深入研究动力堆中子噪声问题打下了基础。

应用B-P材料模型分析高温合金GH536试样的热粘塑性特性

将 Bondner- Partom模型与等效应力函数算法相结合 ,对热粘塑性材料的非弹性变形和应力分析的有限元方法进行了研究 ;并对该模型引入大型工程分析软件 MARC的方法进行了研究。应用所编程序对高温合金 GH5 36带缺口试件的热弹塑性行为进行了分析。研究表明 ,应用 B-

新型3-PRC柔性并联微操作平台的研究

采用压电陶瓷作为驱动器设计了一种带有二级杠杆放大机构的3-PRC柔性并联微操作平台,该平台能实现空间三自由度的微小平移运动,放大机构的理论与有限元仿真分析放大倍数分别为8.772和8.245。首先采用矢量法建立了该机构的运动学模型,得到了其运动学正、逆解和雅克比矩阵。然后利用MATLAB软件分析并绘制了平台的工作空间,最后对机构的解耦性进行了理论分析并利用有限元软件进行了验证。结果表明该机构的耦合性误差小,可以实现空间三自由度微纳米级别的运动,具有较大的运动空间和良好的运动解耦性。

A family of 3D H^2-nonconforming tetrahedral finite elements for the biharmonic equation

In this article,a family of H^2-nonconforming finite elements on tetrahedral grids is constructed for solving the biharmonic equation in 3 D.In the family,the Pl polynomial space is enriched by some high order polynomials for all l≥3 and the corresponding finite element solution converges at the order l-1 in H2 norm.Moreover,the result is improved for two low order cases by using P6 and P7 polynomials to enrich P4 and P5 polynomial spaces,respectively.The error estimate is proved.The numerical results are.provided to confirm the theoretical findings.

Evaluating performance of lignite pillars with 2D approximation techniques and 3D numerical analyses

This paper attempts to investigate the use of approximate 2D numerical simulation techniques for the evaluation of lignite pillar geomechanical response, formed via the room and pillar mining method.Performance and applicability of the developing methodology are assessed through benchmarking with a more direct and accurate 3D numerical model. This analysis utilizes an underground lignite mine which is being developed in soft rock environment. Through the decisions made for the optimum room and pillar layout, the design process highlights the strong points and the weaknesses of 2D finite element analysis, and provides useful recommendations for future reference. The interpretations of results demonstrate that 2D approximation techniques come near quite well to the actual 3D problem.However, external load approximation technique seems to fit even better with the respective outcomes from the 3D analyses.

Optimization of buckling load for laminated composite plates using adaptive Kriging-improved PSO:A novel hybrid intelligent method

An effective hybrid optimization method is proposed by integrating an adaptive Kriging(A-Kriging)into an improved partial swarm optimization algorithm(IPSO)to give a so-called A-Kriging-IPSO for maximizing the buckling load of laminated composite plates(LCPs)under uniaxial and biaxial compressions.In this method,a novel iterative adaptive Kriging model,which is structured using two training sample sets as active and adaptive points,is utilized to directly predict the buckling load of the LCPs and to improve the efficiency of the optimization process.The active points are selected from the initial data set while the adaptive points are generated using the radial random-based convex samples.The cell-based smoothed discrete shear gap method(CS-DSG3)is employed to analyze the buckling behavior of the LCPs to provide the response of adaptive and input data sets.The buckling load of the LCPs is maximized by utilizing the IPSO algorithm.To demonstrate the efficiency and accuracy of the proposed methodology,the LCPs with different layers(2,3,4,and 10 layers),boundary conditions,aspect ratios and load patterns(biaxial and uniaxial loads)are investigated.The results obtained by proposed method are in good agreement with the literature results,but with less computational burden.By applying adaptive radial Kriging model,the accurate optimal resultsebased predictions of the buckling load are obtained for the studied LCPs.