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.

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.

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.

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.

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.

Benchmarking of two three-dimensional numerical models in time/space domain to predict railway-induced ground vibrations

In the last 30 years,the scientific community has developed and proposed different models and numerical approaches for the study of vibrations induced by railway traffic.Most of them are formulated in the frequency/wave number domain and with a 2.5D approach.Three-dimensional numerical models formulated in the time/space domain are less frequently used,mainly due to their high computational cost.Notwithstanding,these models present very attractive characteristics,such as the possibility of considering nonlinear behaviors or the modelling of excess pore pressure and non-homogeneous and non-periodic geometries in the longitudinal direction of the track.In this study,two 3D numerical approaches formulated in the time/space domain are compared and experimentally validated.The first one consists of a finite element approach and the second one of a finite difference approach.The experimental validation in an actual case situated in Carregado(Portugal)shows an acceptable fitting between the numerical results and the actual measurements for both models.However,there are some differences among them.This study therefore includes some recommendations for their use in practical soil dynamics and geotechnical engineering.

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.

Numerical Investigation of Green Water Loading on Flexible Structures Using Three-Step CFD-BEM-FEM Approach

In this paper, the effect of green water impact on a flexible structure is studied based on three-step computational fluid dynamics(CFD)–boundary element method(BEM)–finite element method(FEM) approach. The impact due to shipping of water on the deck of the vessel is computed using commercial CFD software and used as an external force in coupled BEM-FEM solver. Other hydrodynamic forces such as radiation, diffraction, and Froude-Krylov forces acting on the structure are evaluated using 3 D time domain panel method. To capture the structural responses such as bending moment and shear force, 1 D finite element method is developed. Moreover, a direct integration scheme based on the Newmark–Beta method is employed to get the structural velocity,displacement, etc., at each time step. To check the effect of the green water impact on the structure, a rectangular barge without forward speed is taken for the analysis. The influence is studied in terms of bending moment, shear force, etc. Results show that the effect of green water impact on the bow region can be severe in extreme seas and lead to various structural damages. Similarly,it is also verified that vessel motion affects green water loading significantly and therefore one must consider its effect while designing a vessel.

基于FEM-MEI和正余弦算法的二维电磁成像方法

针对二维介质目标的电磁成像问题,将正余弦算法(Sine Cosine Algorithm,SCA)与有限元方法(Finite Element Method,FEM)和不变性测试方程(Measured Equation of Invariance,MEI)进行结合提出一种新的成像方法。将FEM与MEI进行结合求解二维介质目标的电磁散射正问题,即求解Helmholtz方程。其中,MEI保证边界截断的精度,FEM适用于复杂介质目标的准确模拟。对于电磁散射逆问题,引入SCA并加以改进提出一种新的重构方法。该方法采用等效原理与格林函数的渐近式求得远区散射场,以测量的散射场和计算的散射场最大偏差为目标函数,采用改进的SCA优化介质参数,使目标函数达到最小值,以此重构散射体。为提高计算效率,采用MPI算法进行并行计算。文中采用基准函数展示了改进的SCA算法的快速收敛性,并采用非规则的均匀介质柱目标验证了成像方法的正确性。

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 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.

A Three-Dimensional Numerical Model for Predicting the Weld Bead Geometry Characteristics in Laser Overlap Welding of Low Carbon Galvanized Steel

Laser welding (LW) becomes one of the most economical high quality joining processes. LW offers the advantage of very controlled heat input resulting in low distortion and the ability to weld heat sensitive components. To exploit efficiently the benefits presented by LW, it is necessary to develop an integrated approach to identify and control the welding process variables in order to produce the desired weld characteristics without being forced to use the traditional and fastidious trial and error procedures. The paper presents a study of weld bead geometry characteristics prediction for laser overlap welding of low carbon galvanized steel using 3D numerical modelling and experimental validation. The temperature dependent material properties, metallurgical transformations and enthalpy method constitute the foundation of the proposed modelling approach. An adaptive 3D heat source is adopted to simulate both keyhole and conduction mode of the LW process. The simulations are performed using 3D finite element model on commercial software. The model is used to estimate the weld bead geometry characteristics for various LW parameters, such as laser power, welding speed and laser beam diameter. The calibration and validation of the 3D numerical model are based on experimental data achieved using a 3 kW Nd:Yag laser system, a structured experimental design and confirmed statistical analysis tools. The results reveal that the modelling approach can provide not only a consistent and accurate prediction of the weld characteristics under variable welding parameters and conditions but also a comprehensive and quantitative analysis of process parameters effects on the weld quality. The results show great concordance between predicted and measured values for weld bead geometry characteristics, such as depth of penetration, bead width at the top surface and bead width at the interface between sheets, with an average accuracy greater than 95%.

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.

Numerical Investigation of Laser Surface Hardening of AISI 4340 Using 3D FEM Model for Thermal Analysis of Different Laser Scanning Patterns

Laser surface hardening is becoming one of the most successful heat treatment processes for improving wear and fatigue properties of steel parts. In this process, the heating system parameters and the material properties have important effects on the achieved hardened surface characteristics. The control of these variables using predictive modeling strategies leads to the desired surface properties without following the fastidious trial and error method. However, when the dimensions of the surface to be treated are larger than the cross section of the laser beam, various laser scanning patterns can be used. Due to their effects on the hardened surface properties, the attributes of the selected scanning patterns become significant variables in the process. This paper presents numerical and experimental investigations of four scanning patterns for laser surface hardening of AISI 4340 steel. The investigations are based on exhaustive modelling and simulation efforts carried out using a 3D finite element thermal analysis and structured experimental study according to Taguchi method. The temperature distribution and the hardness profile attributes are used to evaluate the effects of heating parameters and patterns design parameters on the hardened surface characteristics. This is very useful for integrating the scanning patterns’ features in an efficient predictive modeling approach. A structured experimental design combined to improved statistical analysis tools is used to assess the 3D model performance. The experiments are performed on a 3 kW Nd:Yag laser system. The modeling results exhibit a great agreement between the predicted and measured values for the hardened surface characteristics. The model evaluation reveals also its ability to provide not only accurate and robust predictions of the temperature distribution and the hardness profile as well an in-depth analysis of the effects of the process parameters.

Developing a 3D-FEM Model for Electromagnetic Analysis of an Axial Flux Permanent Magnet Machine

Recently, many optimal designs for axial flux permanent magnet (AFPM) motors were performed based on finite- element (FE) analysis. Most of the models are based on reduction of 3D problem to 2D problem which is not accurate for design aspects. This paper describes an accurate electromagnetic analysis of a surface mounted, 28 pole AFPM with concentrated stator winding. The AFPM is modeled with three-dimensional finite-element method. This model in-cludes all geometrical and physical characteristics of the machine components. Using this accurate modeling makes it possible to obtain demanded signals for a very high precision analysis. Magnetic flux density, back-EMF, magnetic axial force and cogging torque of the motor are simulated using FLUX-3D V10.3.2. Meanwhile, the model is paramet-ric and can be used for design process and sensitivity analysis.

熔融堆积3D打印成形金属件层间结合研究

金属件熔融堆积3D打印过程中,制件的层间结合性能主要取决于热作用过程,因此对成形过程温度变化进行研究显得尤为重要。该研究基于熔融堆积3D打印成形特点,建立了成形过程有限元分析数学及物理模型,并使用ANSYS有限元软件对不同熔融金属温度、基板温度及堆积速度条件下成形过程温度场变化进行模拟研究。结合相同参数条件下的对比工艺试验,研究了这些参数对最终成形金属实体层间结合性能的影响。研究结果表明:随着熔融金属温度、基板温度的升高,以及堆积速度的增加,成形实体温度上升速度加快,高温热影响区增大,温度梯度减小,实体层间结合及拉伸性能提高,并在熔融金属温度160℃,基板温度90℃,堆积速度16mm/s参数条件下打印出了层间结合良好的铋锡合金实体。

Contrast between 2D inversion and 3D inversion based on 2D high-density resistivity data

The 2D data processing adopted by the high-density resistivity method regards the geological structures as two degrees, which makes the results of the 2D data inversion only an approximate interpretation;the accuracy and effect can not meet the precise requirement of the inversion. Two typical models of the geological bodies were designed, and forward calculation was carried out using finite element method. The forward-modeled profiles were obtained. 1% Gaussian random error was added in the forward models and then 2D and 3D inversions using a high-density resistivity method were undertaken to realistically simulate field data and analyze the sensitivity of the 2D and 3D inversion algorithms to noise. Contrast between the 2D and 3D inversion results of least squares inversion shows that two inversion results of high-density resistivity method all can basically reflect the spatial position of an anomalous body. However, the 3D inversion can more effectively eliminate the influence of interference from Gaussian random error and better reflect the distribution of resistivity in the anomalous bodies. Overall, the 3D inversion was better than 2D inversion in terms of embodying anomalous body positions, morphology and resistivity properties.

3D GIS与有限元模拟无缝耦合方法及其在隧道工程中的应用研究

为实现数值模拟获得的各种场数据与GIS空间信息的综合集成,有必要对GIS中融合有限元数值计算的方法开展研究。本文采用GRASS GIS、VTK等开源软件,通过Python语言建立了GIS与有限元模拟的复合系统开发框架,针对隧道开挖编写了3DGIS与有限元数值计算无缝耦合分析程序(3DPyGisTFEM)。给出了由中心线上多条轮廓线生成隧道模型的算法,提出了对各地层DEM与隧道模型进行剖切运算获得隧道开挖地质剖面的方法。引入Delaunay网格剖分算法,对GIS中TIN模型不断改良,获得了适合有限元模拟的地质剖面三角形网格。基于多个剖面计算结果,运用插值手段,形成研究区域场数据,结合GIS空间分析方法,对场数据分布规律进行综合分析,实现了GIS与有限元的无缝耦合。黄家寨隧道工程实例表明,相对坡度、隧道中心线缓冲区等因素,地表距隧道拱顶高差的空间分布特征与地表沉降规律相关性更强。

基于MeshPy的3DGIS与三维有限元数值计算无缝耦合方法

3DGIS与有限元数值计算耦合方法及其系统开发,是目前GIS与岩土工程交叉学科的重要研究内容。通过对MeshPy、GRASSGIS、vtk Python以及Scipy等多个开源库的探索与编程实践,以Python为"黏合剂",构建了GIS与三维有限元数值计算无缝耦合的平台框架。在满足几何形状、Delaunay特性、角度、面积以及体积等约束条件下,通过MeshPy包含的Triangle与TetGen库获得能用于有限元数值计算的三角形与四面体网格,并阐述了限定Delaunay三角化(CDT)节点插入算法与四面体网格质量控制方法。以排土场为例,针对CDT网格建立的地表模型,结合vtkPython封装的模块与算法,提出了根据节点分布位置综合采用线性与三角形重心坐标的插值方法,有效消除了其他插值手段导致的部分地表高程畸变的现象。借助Scipy提供的矩阵运算模块,开发了邓肯-张E-?模型,给出了采用中点增量法进行非线性求解的算法,并通过对室内大型三轴剪切试验的数值模拟验证了程序的正确性。最后,对准东露天煤矿北部排土场堆排过程开展三维有限元数值计算,结合GRASSGIS对排弃物料高度与厚度空间分布规律的分析结果,对排土场边坡位移变化规律进行初步研究,结果表明,排土场垂直位移与排弃料高度及厚度相关,而最大水平位移则出现在各台阶坡体边缘位置,与实际情况一致。

Improving the structural efficiency of textured threedimensional concrete printing wall by architectural design

Three-dimensional concrete printing(3DCP)technology begins to be adopted into construction application worldwide.Recent studies have focused on producing a higher concrete quality and offering a user-friendly construction process.Still,the 3DCP construction cost is unlikely to be lower than that of conventional construction,which is especially important for projects where the cost is sensitive.To broaden the 3DCP construction applications,reduction of the quantity of 3DCP material usage is needed.This work aims to perform structural analysis of several patterns of geometric textured 3DCP shell wall structures.21 different cantilevered textured patterns of 3DCP shell wall structures were architecturally designed and then subjected to structural analysis by a finite element method(FEM).The results indicated that by designing appropriate patterns,the structural performance to weight ratio could be improved up to 300%.The study therefore offers an innovative design process for constructing 3DCP housing and suggests preconstruction analysis methods for 3DCP shell wall structures.