A three-dimensional semi-implicit unstructured grid finite volume ocean model

A new three-dimensional semi-implicit finite-volume ocean model has been developed for simulating the coastal ocean circulation, which is based on the staggered C-unstructured non-orthogonal grid in the hor- izontal direction and z-level grid in the vertical direction. The three-dimensional model is discretized by the semi-implicit finite-volume method, in that the free-surface and the vertical diffusion are semi-implicit, thereby removing stability limitations associated with the surface gravity wave and vertical diffusion terms. The remaining terms in the momentum equations are discretized explicitly by an integral method. The partial cell method is used for resolving topography, which enables the model to better represent irregular topography. The model has been tested against analytical cases for wind and tidal oscillation circulation, and is applied to simulating the tidal flow in the Bohal Sea. The results are in good agreement both with the analytical solutions and measurement results.

Three-dimensional forward modeling for magnetotelluric sounding by finite element method

A finite element algorithm combined with divergence condition was presented for computing three-dimensional(3D) magnetotelluric forward modeling. The finite element equation of three-dimensional magnetotelluric forward modeling was derived from Maxwell's equations using general variation principle. The divergence condition was added forcedly to the electric field boundary value problem, which made the solution correct. The system of equation of the finite element algorithm was a large sparse, banded, symmetric, ill-conditioned, non-Hermitian complex matrix equation, which can be solved using the Bi-CGSTAB method. In order to prove correctness of the three-dimensional magnetotelluric forward algorithm, the computed results and analytic results of one-dimensional geo-electrical model were compared. In addition, the three-dimensional magnetotelluric forward algorithm is given a further evaluation by computing COMMEMI model. The forward modeling results show that the algorithm is very efficient, and it has a lot of advantages, such as the high precision, the canonical process of solving problem, meeting the internal boundary condition automatically and adapting to all kinds of distribution of multi-substances.

Three-Dimensional Simulations of RESET Operation in Phase-Change Random Access Memory with Blade-Type Like Phase Change Layer by Finite Element Modeling

An optimized device structure for reducing the RESET current of phase-change random access memory （PCRAM） with blade-type like （BTL） phase change layer is proposed. The electrical thermal analysis of the BTL cell and the blade heater contactor structure by three-dimensional finite element modeling are compared with each other during RESET operation. The simulation results show that the programming region of the phase change layer in the BTL cell is much smaller, and thermal electrical distributions of the BTL cell are more concentrated on the TiN/GST interface. The results indicate that the BTL cell has the superiorities of increasing the heating efficiency, decreasing the power consumption and reducing the RESET current from 0.67mA to 0.32mA. Therefore, the BTL cell will be appropriate for high performance PCRAM device with lower power consumption and lower RESET current.

Autologous nerve graft repair of different degrees of sciatic nerve defect:stress and displacement at the anastomosis in a three-dimensional finite element simulation model

In the repair of peripheral nerve injury using autologous or synthetic nerve grafting, the mag- nitude of tensile forces at the anastomosis affects its response to physiological stress and the ultimate success of the treatment. One-dimensional stretching is commonly used to measure changes in tensile stress and strain; however, the accuracy of this simple method is limited. There- fore, in the present study, we established three-dimensional finite element models of sciatic nerve defects repaired by autologous nerve grafts. Using PRO E 5.0 finite element simulation software, we calculated the maximum stress and displacement of an anastomosis under a 5 N load in 10-, 20-, 30-, 40-mm long autologous nerve grafts. We found that maximum displacement increased with graft length, consistent with specimen force. These findings indicate that three-dimensional finite element simulation is a feasible method for analyzing stress and displacement at the anas- tomosis after autologous nerve grafting.

Evaluation value of three-dimensional finite element model analysis for bone mineral density and bone metabolism activity in patients with osteoporosis

Objective: To study the evaluation value of three-dimensional finite element model analysis for bone mineral density (BMD) and bone metabolism activity in patients with osteoporosis. Methods: A total of 218 patients who were diagnosed with osteoporosis in the hospital between February 2014 and January 2017 were collected as observation group, and 100 healthy volunteers who received physical examination in the hospital during the same period were selected as normal control group. The femoral head of the two groups was analyzed by three-dimensional finite element model, and the femoral head BMD levels and serum bone metabolism index contents were measured. Pearson test was used to evaluate the evaluation value of femoral head three-dimensional finite element model for osteoporosis. Results: The cancellous bone and cortical bone Von Mises stress value of observation group were lower than those of normal control group, and femoral neck BMD value of observation group was lower than that of normal control group;serum bone metabolism index BGP content was lower than that of normal control group while NBAP, TRACP-5b and CTX-1 contents were higher than those of normal control group. Pearson test showed that the cancellous bone and cortical bone Von Mises stress value of patients with osteoporosis were directly correlated with BMD value and bone metabolism index contents. Conclusion: The three-dimensional finite element model analysis resultsof patients with osteoporosis can objectively reflect the femoral headBMD value and bone metabolism activity, and is a reliable way to evaluate the risk of long-term fractures.

Three-dimensional finite element analysis on effects of tunnel construction on nearby pile foundation

A three-dimensional finite element simulation was carried out to investigate the effects of tunnel construction on nearby pile foundation.The displacement controlled model (DCM) was used to simulate the tunneling-induced volume loss effects.The numerical model was verified based on the results of a centrifuge test and a set of parametric studies was implemented based on this model.There is good agreement between the trend of the results of the centrifuge test and the present model.The results of parametric studies show that the tunnelling-induced pile internal force and deformation depend mainly on the pile?tunnel distance,the pile length to tunnel depth ratio and the volume loss.Two different zones are separated by a 45° line projected from the tunnel springline.Within the zone of influence,the pile is subjected to tensile force and large settlement;whereas outside the zone of influence,dragload and small settlement are induced.It is also established that the impact of tunnelling on a pile group is substantially smaller as compared with a single pile in the same location with the rear pile in a group,demonstrating a positive pile group effect.

Three-Dimensional Thermal-Stress Analysis of Semi-infinite Transversely Isotropic Composites

By making use of the direct integration method,an exact analysis of the general three-dimensional thermoelasticity problem is performed for the case of a transversely isotropic homogeneous half-space subject to local thermal and force loadings.The material plane of isotropy is assumed to be parallel to the limiting surface of the halfspace.By reducing the original thermoelasticity equations to the governing ones for individual stress-tensor components,the effect of material anisotropy in the stress field is analyzed with regard to the feasibility requirement,i.e.,the finiteness of the stress field at a distance from the disturbed area.As a result,the solution is constructed in the form of explicit analytical dependencies on the force and thermal loadings for various kinds of transversely isotropic materials and agrees with the basic principles of the continua mechanics.The solution can be efficiently used as a benchmark one for the direct computation of temperature and thermal stresses in transversely isotropic semi-infinite domains,as well as for the verification of solutions constructed by different means.

Simulation-Based Construction of Three-Dimensional Process Model for Punching Cartridge Cases

A method of constructing three-dimensional process model for the punching cartridge cases is presented based on DEFORM simulation analysis. Using DEFORM software,the finite element simulation models for the punching and forming process of cartridge cases are established,and the corresponding simulation result model of each intermediate procedure is obtained by continuously performing the forming process simulation. The simulation model cannot annotate size and process information due to poor interface between DEFORM software and CAD software. Thus,a 3D annotation module is developed with secondary development technology of UG NX software. Consequently,the final process model with dimension and process information is obtained. Then,with the current 3D process management system,the 3D punching and forming process design of cartridge cases can be completed further. An example is also provided to illustrate that the relative error between the simulation process model and the physical model is less than 2%,which proves the validity and reliability of the proposed method in this study.

Implementation of a particle-in-cell method for the energy solver in 3D spherical geodynamic modeling

The thermal evolution of the Earth’s interior and its dynamic effects are the focus of Earth sciences.However,the commonly adopted grid-based temperature solver is usually prone to numerical oscillations,especially in the presence of sharp thermal gradients,such as when modeling subducting slabs and rising plumes.This phenomenon prohibits the correct representation of thermal evolution and may cause incorrect implications of geodynamic processes.After examining several approaches for removing these numerical oscillations,we show that the Lagrangian method provides an ideal way to solve this problem.In this study,we propose a particle-in-cell method as a strategy for improving the solution to the energy equation and demonstrate its effectiveness in both one-dimensional and three-dimensional thermal problems,as well as in a global spherical simulation with data assimilation.We have implemented this method in the open-source finite-element code CitcomS,which features a spherical coordinate system,distributed memory parallel computing,and data assimilation algorithms.

Microstructure-based three-dimensional characterization of chip formation and surface generation in the machining of particulate-reinforced metal matrix composites

Particulate-reinforced metal matrix composites(PRMMCs)are difficult to machine due to the inclusion of hard,brittle reinforcing particles.Existing experimental investigations rarely reveal the complex material removal mechanisms(MRMs)involved in the machining of PRMMCs.This paper develops a three-dimensional(3D)microstructure-based model for investigating the MRM and surface integrity of machined PRMMCs.To accurately mimic the actual microstructure of a PRMMC,polyhedrons were randomly distributed inside the matrix to represent irregular SiC particles.Particle fracture and matrix deformation and failure were taken into account.For the model’s capability comparison,a two-dimensional(2D)analysis was also conducted.Relevant cutting experiments showed that the established 3D model accurately predicted the material removal,chip morphology,machined surface finish,and cutting forces.It was found that the matrix-particle-tool interactions led to particle fractures,mainly in the primary shear and secondary deformation zones along the cutting path and beneath the machined surface.Particle fracture and dilodegment greatly influences the quality of a machined surface.It was also found that although a 2D model can reflect certain material removal features,its ability to predict microstructural variation is limited.

A Three-Dimensional Prediction Method for Thermal Diffusion

A three-dimensional, first order turbulence closure, thermal diffusion model is described in this paper. The governing equations consist of an equation of continuity, three components of momentum, conservation equations for salt, temperature and subgridscale energy, and an equation of state. In the model, according to the hypothesis of Kolmogorov and Prandtl, the viscosity coefficient of turbulent flow of homogeneous fluid is related to the local turbulent energy, and the horizontal and vertical exchange coefficients of mass, heat and momentum are computed with the introduction of subgridscale turbulence energy. The governing equations are solved by finite difference techniques. This model is applied to the Jiaozhou bay to predict thermal pollution by the Huangdao power plant. An instantaneous tidal current field is computed, then the distribution of temperature increment is predicted, and finally the effect of wind stress on thermal discharge is discussed.

3D finite element analysis of a two-surface wear model in fretting tests

This article aims at developing a computationally efficient framework to simulate the erosion of two contact surfaces in three-dimensional(3D),depending on the body resistance.The framework involves finite element(FE)resolution of a fretting problem,wear computation via a non-local criterion including a wear distribution parameter(WDP),as well as updating of the geometry and automatic remeshing.Its originality is based on the capability to capture the damage on each surface and obtain local and global results for a quantitative and qualitative analysis.Numerical simulations are carried out for two 3D contact specimens with different values of WDP.The results highlight the importance of correctly modelling wear:One-surface wear model is sufficient from a global point of view(wear volume),or whenever the wear resistance for a body is much higher than that of another one,whereas a 3D two-surface wear model is essential to capturing local effects(contact pressure,wear footprint,etc.)related to the difference in wear resistance of the bodies.

FINITE ELEMENT MODELLING OF COMPLEX 3D STATIC AND DYNAMIC CRACK PROPAGATION BY EMBEDDING COHESIVE ELEMENTS IN ABAQUS

This study proposes an algorithm of embedding cohesive elements in Abaqus and develops the computer code to model 3D complex cragk propagation in quasi-brittle materials in a relatively easy and efficient manner. The cohesive elements with softening traction-separation relations and damage initiation and evolution laws are embedded between solid elements in regions of interest in the initial mesh to model potential cracks. The initial mesh can consist of tetrahedrons, wedges, bricks or a mixture of these elements. Neither remeshing nor objective crack propagation criteria are needed. Four examples of concrete specimens, including a wedgesplitting test, a notched beam under torsion, a pull-out test of an anchored cylinder and a notched beam under impact, were modelled and analysed. The simulated crack propagation processes and load-displacement curves agreed well with test results or other numerical simulations for all the examples using initial meshes with reasonable densities. Making use of Abaqus＇s rich pre/post- processing functionalities and powerful standard/explicit solvers, the developed method offers a practical tool for engineering analysts to model complex 3D fracture problems.

Experimental and numerical evaluation of multidirectional compressive and flexure behavior of three-dimensional printed concrete

Three-dimensional concrete printing(3DCP)can proliferate the industrialization of the construction sector,which is notoriously conservative and indolent toward changes.However,the mechanical behavior of 3DCP should be characterized and modeled considering the interfaces when its performance is thoroughly compared to that of the existing concrete construction methods.This study presents an experimental and numerical investigation of uniaxial compression and three-point bending(TPB)tests on extruded 3DCP beams in different loading directions.The orientation of translational and depositional interfaces with respect to the direction of loading influenced the strength.Both the elastic and post-damage behavior of the 3DCP specimens were compared with those of the conventionally cast specimen under quasi-static loading conditions.Despite the higher compressive strength of the casted specimen,the flexural strength of the 3DCP specimens was higher.This study employed the finite element and cohesive zone models of the appropriate calibrated traction-separation law to model fracture in the notched TPB specimens.Furthermore,the real-time acoustic emission test revealed the nature of failure phenomenon of three-dimensional-printed specimens under flexion,and accordingly,the cohesive law was chosen.The predicted load-displacement responses are in good agreement with the experimental results.Finally,the effects of cohesive thickness and notch shape on the performance under bending were explored through parametric studies.

Stochastic analysis of excavation-induced wall deflection and box culvert settlement considering spatial variability of soil stiffness

In this study, a three-dimensional (3D) finite element modelling (FEM) analysis is carried out to investigate the effects of soil spatial variability on the response of retaining walls and an adjacent box culvert due to a braced excavation. The spatial variability of soil stiffness is modelled using a variogram and calibrated by high-quality experimental data. Multiple random field samples (RFSs) of soil stiffness are generated using geostatistical analysis and mapped onto a finite element mesh for stochastic analysis of excavation-induced structural responses by Monte Carlo simulation. It is found that the spatial variability of soil stiffness can be described by an exponential variogram, and the associated vertical correlation length is varied from 1.3 m to 1.6 m. It also reveals that the spatial variability of soil stiffness has a significant effect on the variations of retaining wall deflections and box culvert settlements. The ignorance of spatial variability in 3D FEM can result in an underestimation of lateral wall deflections and culvert settlements. Thus, the stochastic structural responses obtained from the 3D analysis could serve as an effective aid for probabilistic design and analysis of excavations.

基于结合面的龙门五面加工中心虚拟建模及静刚度研究

静刚度是机床性能的重要指标,有必要在设计阶段对其进行预估。针对正在设计中的Linear MC6000型龙门五面加工中心机床,在CAD软件UGNX2.0环境中建立了其各部件及整机的三维实体模型。在有限元分析软件ABAQUS6.51环境中,用弹簧单元来模拟滑块与直线导轨的接触刚度,建立了机床的有限元模型,并对其静刚度进行了分析。

基于ANSYS的塔式起重机臂架有限元参数化建模与分析

针对塔式起重机在有限元分析过程中建模工作量大的问题,将参数化思想引入有限元分析过程中,采用有限元参数化分析文本自动生成技术,以通用编程语言为基础编制程序,根据输入参数自动生成有限元建模分析所需的参数化文件,进行结构建模与分析,并以塔机臂架结构为例说明方法的可行性。

高压水射流作用下岩石破碎机理及过程的数值模拟研究

根据连续介质力学和有限元理论,给出了高压水射流破岩系统中流体和岩石的控制方程,并建立了相应的有限元列式.运用连续损伤力学和细观损伤力学理论,建立了适用于水射流破岩全过程分析的岩石损伤模型以及宏细观损伤的耦合模式.数值计算的结果较真实地反映了水射流破岩过程中,岩石的动态响应以及水射流动力学特性的演化过程,普通连续水射流破碎岩石主体所用的时间为毫秒量级,破岩的主要形式是卸载及射流冲击所产生的拉伸破坏,并呈“阶跃式”发展.数值计算与相关试验结果基本吻合,表明该分析方法是可行的,可用来指导高压水射流破岩理论的进一步研究及应用.

受面内载荷的层压复合材料的层间应力分析

研究任意多向铺层层压复合材料层间应力的有限元分析方法。分别考虑两种建模方法,一种是采用8结点等参实体单元逐层进行离散化,每个子层都视为正交异性材料;另一种方法是在子层面内采用4结点板元进行离散化,而在层间采用刚性元和弹簧元进行离散化,称为准三维模型。计算实例为两种不同铺层的T300QY8911碳纤维双马树脂层压板条,其中一个板条连续,另一个含圆孔,均受到面内拉伸载荷作用。采用的计算软件为MSC.Nastran。计算结果表明,在受拉板条的自由边和孔口附近存在不能忽视的层间应力。两种建模方法有各自的特点,计算结果可以互相印证,而后一种方法更加灵活和有更广泛的适应性。

裂纹梁动力学仿真的有限元模型

针对弥散裂缝模型对裂纹影响区域没有定义的不足,基于应变能等效原则,提出了一种确定裂纹单元长度的计算方法,使这种传统、粗糙的建模方法更加实用.得出了裂纹影响区域相对裂纹深度的变化规律,当量纲一裂纹深度为0.4时,裂纹影响区域最大,约为0.443倍梁高.算例表明,采用改进模型的仿真结果和试验结果吻合较好.