Analysis of ground vibrations due to underground trains by 2.5D finite/infinite element approach

The 2.5D finite/infinite element approach is adopted to study wave propagation problems caused by underground moving trains. The irregularities of the near field, including the tunnel structure and parts of the soil, are modeled by the finite elements, and the wave propagation properties of the far field extending to infinity are modeled by the infinite elements. One particular feature of the 2.5D approach is that it enables the computation of the three-dimensional response of the half-space, taking into account the load-moving effect, using only a two-dimensional profile. Although the 2.5D finite/infinite element approach shows a great advantage in studying the wave propagation caused by moving trains, attention should be given to the calculation aspects, such as the rules for mesh establishment, in order to avoid producing inaccurate or erroneous results. In this paper, some essential points for consideration in analysis are highlighted, along with techniques to enhance the speed of the calculations. All these observations should prove useful in making the 2.5D finite/infinite element approach an effective one.

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.

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.

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.

2.5-D modeling of cross-hole electromagnetic measurement by nite element method

A finite element method is developed for simulating frequency domain electromagnetic responses due to a dipole source in the 2-D conductive structures. Computing costs are considerably minimized by reducing the full three-dimensional problem to a series of two-dimensional problems. This is accomplished by transforming the problem into y-wave number （Ky） domain using Fourier transform and the y-axis is parallel to the structural strike. In the Ky domain, two coupled partial differential equations for magnetic field Hy and electric field Ey are derived. For a specific value of Ky, the coupled equations are solved by the finite element method with isoparametric elements in the x-z plane. Application of the inverse Fourier transform to the Ky, domain provides the electric and magnetic fields in real space. The equations derived can be applied to general complex two-dimensional structures containing either electric or magnetic dipole source in any direction. In the modeling of the electromagnetic measurement, we adopted a pseudo-delta function to distribute the dipole source current and circumvent the problem of singularity at the source point. Moreover, the suggested method used isoparametric finite elements to accommodate the complex subsurface formation. For the large scale linear system derived from the discretization of the Maxwell＇s equations, several iterative solvers were used and compared to select the optimal one. A quantitative test of accuracy was presented which compared the finite element results with analytic solutions for a dipole source in homogeneous space for different ranges and different wave numbers Ky. to validate the addressed the effects of the distribution range τ of the homogeneous medium. code and check its effectiveness. In addition, we pseudo-delta function on the numerical results in

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.

带式输送机不停机托辊更换车结构设计及特性分析

为提高煤矿带式输送机在更换损坏托辊时的维护效率,减小人工更换产生的安全风险,文中设计了一种可在带式输送机不停机状态下进行托辊更换的更换车。针对不停机更换托辊的工作要求,确定更换车的更换方案,并建立三维模型。采用D-H参数法建立更换车机械手的连杆坐标系,通过MATLAB软件机器人工具箱建立机械手模型,并用蒙特卡洛法确定机械手的工作空间,结果表明,其运动范围满足槽型托辊组3个托辊的更换。采用ADAMS软件对机械手拆卸和安装托辊过程进行仿真,得到损坏托辊和备用托辊的位移曲线,结合夹取结构的运动曲线和运行轨迹,结果表明,机械手的工作过程平稳流畅。采用ANSYS软件对更换车工作中关键受力部件支撑臂和支撑托辊架进行有限元分析,结果表明:两部件强度满足工作要求,保证更换车工作的可靠性。

降雨作用下抗滑桩与三维土工网垫植草护坡坡体的稳定分析

针对抗滑桩与三维土工网垫植草护坡在降雨作用下的稳定性问题,采用有限元强度折减法,利用ABAQUS建立抗滑桩与三维土工网垫植草护坡坡体的三维有限元模型,以坡底水平位移突变情况作为边坡失稳的判断标准,分别比较了不同工况下护坡安全系数随降雨时长的变化规律,以及降雨作用下抗滑桩桩位、桩长对抗滑桩与三维土工网垫植草护坡坡体的稳定性影响,并与未考虑降雨时的影响规律进行对比。结果表明:降雨时长对三维土工网垫植草工况的边坡稳定性影响较为明显;抗滑桩桩长对边坡稳定性影响较大,桩长宜取8 m;无降雨作用下将桩位置于坡中时边坡整体获得最大安全系数,降雨作用下将抗滑桩置于相对位置(L 0.4=4.8 m)时能够起到更好的抗滑效果;抗滑桩与三维土工网垫植草在降雨作用时加固护坡的效果比无降雨时更好,且优于原边坡。

高铁荷载下路堤和路堑段地面振动特性现场测试与数值分析

对巴黎-布鲁塞尔高速铁路地面振动开展现场测试,对比分析了高铁运行引起的路堤和路堑段地面振动特性以及传播衰减规律。基于2.5维有限元法基本原理,推导并建立了高铁荷载下路堤和路堑段2.5维有限元动力计算模型,详细讨论了路堤和路堑设计参数对地面振动特性的影响。结果表明:距轨道较近处的地面振动幅值主要受列车轴重影响,而距轨道较远处则由中间车相邻转向架轮对的叠加效应主导。高铁运行引起的路堤和路堑段地面振动随与轨道距离增加而减小,中高频成分振动的衰减速率显著大于其他频段。路堤和路堑段地面振动1阶主频主要由列车荷载的基频f_1主导,随车速提高逐渐向高频方向移动。路堤和路堑对地面振动的衰减规律有重要影响,路堤段地面振动随与轨道间距的衰减过程可分为两个阶段,而路堑段地面振动衰减过程则可分为三个阶段,且振动在堑顶处有较为明显的局部放大现象。当与路堤段轨道间距≥19.0 m(路堑段≥23.0 m)后,地面振动Z振级VL_Z均小于80.0 dB。同时,高铁荷载下地面振动随路堤高度或路堑深度的增加而减小,当超过某一限值时,继续增加路堤高度或路堑深度难以有效减小地面振动。路堤段地面振动随路堤弹性模量的增大而减小。此外,在确保稳定性的前提下,较陡的路堑边坡更有利于减小高铁运行引起的地面振动。

磁矢量势的全球电磁感应时间域三维有限元正演计算方法

全球电磁感应测深方法能获得地球深部电导率结构分布,目前被广泛应用于地球内部结构与热状态研究.地磁台站与地磁卫星观测的电磁感应数据为时间序列信号,在时间域分析全球电磁感应数据,特别是分析频谱极宽的磁暴脉冲响应,具有天然的优势.然而,当前的全球电磁感应数据解释技术一般在频率域进行,缺少时间域中的研究成果.为了弥补缺少时间域全球电磁感应数据解释方法的问题,本文开发了一种基于磁场矢量势的全球电磁感应时间域三维有限元并行正演求解方法,具备高精度快速计算源于地球外部时变电流源的地球感应电磁场时间序列的能力,特别适合于计算与分析频谱极宽的磁暴脉冲电磁感应时变响应.首先,建立基于磁场矢量势的全球电磁感应时间域微分控制方程,结合磁层外部电流源的物理属性建立边界条件和初始条件,从而构建出全球电磁感应时间域初始边界值问题.然后,利用四面体矢量有限元技术和无条件稳定的隐式后退欧拉公式,分别实现磁场矢量势的空间域和时间域离散,进而获得不同时刻的实系数大型有限元线性方程组,借助于高性能并行直接求解器,快速高精度地求解不同时刻的磁场矢量势与感应磁场.最后,利用理论模型验证本文算法的正确性.利用Dst磁暴环电流指数建立的时间域电流源与真实地球三维电性模型,研究"澳科一号"后续卫星200 km轨道对中国和澳大利亚下方地幔转换带高导体的探测能力,结果表明这些地幔转换带异常体在200 km卫星轨道高度能产生明显异常.综上所述,本文开发的时间域全球电磁感应方法不仅具备精确快速计算全球感应电磁场时间序列响应的能力,还能为反演与解释"澳科一号"等我国地磁卫星观测数据提供技术支持.

基于非结构化网格的2.5D直流电阻率模拟

结合大型通用有限元软件ANSYS所提供的强大网格剖分功能，实现了2D地质模型的非结构化三角形网格剖分，并编写了2．5D直流电阻率有限元法FORTRAN语言计算程序。经计算对比表明，在节点数基本相同时，非结构化三角形网格比传统的矩形网格、矩形-三角形网格具有更高的计算精度，可以很好地拟合复杂地形和地质体边界。在此基础上，采用非结构化三角形网格，计算了起伏地形下2．5D地电模型的视电阻率异常，并利用比较法进行了地形改正。与水平地形时的结果对比表明，比较法可以较好地消除地形影响，突出局部地质体的异常。

基于各向异性和极化介质的2.5维CSAMT正演数值模拟

传统的CSAMT在进行正、反演数值模拟时通常假设地下介质为电各向同性且无IP效应的介质,而实际的地下介质往往同时存在电各向异性和IP效应,因此开展各向异性和IP效应共同作用下的CSAMT电磁场响应特征的研究势在必行。文中基于含源的Maxwell方程组,推导了2.5维各向异性介质的电磁场偏微分方程,采用伽辽金有限元法实现了各向异性的2.5维CSAMT数值模拟,进而讨论了欧拉角和各向异性系数对2.5维CSAMT响应的影响。通过引入Cole-Cole模型分别研究了不同极化参数对CSAMT响应的影响,最后讨论了各向异性和IP效应同时存在时的电磁场响应特征。研究结果表明,各向异性和IP效应均会对CSAMT的响应结果产生较为严重的影响。当各向异性和IP效应同时存在时,根据欧拉角和主轴电阻率的取值,二者的效果将相互叠加或抵消。文中的研究结果对于提高CSAMT资料处理与解释精度具有重要的理论意义。

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.

2.5D forward modeling and inversion of frequency-domain airborne electromagnetic data

Frequency-domain airborne electromagnetics is a proven geophysical exploration method.Presently,the interpretation is mainly based on resistivity-depth imaging and onedimensional layered inversion;nevertheless,it is difficult to obtain satisfactory results for two- or three-dimensional complex earth structures using 1D methods.3D forward modeling and inversion can be used but are hampered by computational limitations because of the large number of data.Thus,we developed a 2.5D frequency-domain airborne electromagnetic forward modeling and inversion algorithm.To eliminate the source singularities in the numerical simulations,we split the fields into primary and secondary fields.The primary fields are calculated using homogeneous or layered models with analytical solutions,and the secondary（scattered） fields are solved by the finite-element method.The linear system of equations is solved by using the large-scale sparse matrix parallel direct solver,which greatly improves the computational efficiency.The inversion algorithm was based on damping leastsquares and singular value decomposition and combined the pseudo forward modeling and reciprocity principle to compute the Jacobian matrix.Synthetic and field data were used to test the effectiveness of the proposed method.

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.

Analytical Model of Elastic Modulus and Coefficient of Thermal Expansion for 2.5D C/SiC Composite

To make better use of 2.5D C/SiC composites in industry, it is necessary to understand the mechanical properties. A finite element model＇of 2.5D composites is established, by considering the fiber undulation and the porosity in 2.5D C/SiC composites. The fiber direction of warp is defined by cosine function to simulate the undulation of warp, and based on uniform strain assumption, analytical model of the elastic modulus and coefficient of thermal expansion （CTE） for 2.5D C/SiC composites were established by using dual- scale model. The result is found to correlate reasonably well with the predicted results and experimental results. The parametric study also demonstrates the effects of the fiber volume fraction, distance of warp yarn, and porosity in micro-scale on the mechanical properties and the coefficients of thermal expansion.

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.

Nonlinear Inversion for Complex Resistivity Method Based on QPSO-BP Algorithm

The significant advantage of the complex resistivity method is to reflect the abnormal body through multi-parameters, but its inversion parameters are more than the resistivity tomography method. Therefore, how to effectively invert these spectral parameters has become the focused area of the complex resistivity inversion. An optimized BP neural network (BPNN) approach based on Quantum Particle Swarm Optimization (QPSO) algorithm was presented, which was able to improve global search ability for complex resistivity multi-parameter nonlinear inversion. In the proposed method, the nonlinear weight adjustment strategy and mutation operator were used to enhance the optimization ability of QPSO algorithm. Implementation of proposed QPSO-BPNN was given, the network had 56 hidden neurons in two hidden layers (the first hidden layer has 46 neurons and the second hidden layer has 10 neurons) and it was trained on 48 datasets and tested on another 5 synthetic datasets. The training and test results show that BP neural network optimized by the QPSO algorithm performs better than the BP neural network without initial optimization on the inversion training and test models, and the mean square error distribution is better. At the same time, a double polarized anomalous bodies model was also used to verify the feasibility and effectiveness of the proposed method, the inversion results show that the QPSO-BP algorithm inversion clearly characterizes the anomalous boundaries and is closer to the values of the parameters.