Due to the low permeability of tight reservoirs,throats play a significant role in controlling fluid flow.Although many studies have been conducted to investigate fluid flow in throats in the microscale domain,compara...Due to the low permeability of tight reservoirs,throats play a significant role in controlling fluid flow.Although many studies have been conducted to investigate fluid flow in throats in the microscale domain,comparatively fewer works have been devoted to study the effect of adsorption boundary layer(ABL)in throats based on the digital rock method.By considering an ABL,we investigate its effects on fluid flow.We build digital rock model based on computed tomography technology.Then,microscopic pore structures are extracted with watershed segmentation and pore geometries are meshed through Delaunay triangulation approach.Finally,using the meshed digital simulation model and finite element method,we investigate the effects of viscosity and thickness of ABL on microscale flow.Our results demonstrate that viscosity and thickness of ABL are major factors that significantly hinder fluid flow in throats.展开更多
Determination of ballistic performance of an armor solution is a complicated task and evolved significantly with the application of finite element methods(FEM) in this research field.The traditional armor design studi...Determination of ballistic performance of an armor solution is a complicated task and evolved significantly with the application of finite element methods(FEM) in this research field.The traditional armor design studies performed with FEM requires sophisticated procedures and intensive computational effort,therefore simpler and accurate numerical approaches are always worthwhile to decrease armor development time.This study aims to apply a hybrid method using FEM simulation and artificial neural network(ANN) analysis to approximate ballistic limit thickness for armor steels.To achieve this objective,a predictive model based on the artificial neural networks is developed to determine ballistic resistance of high hardness armor steels against 7.62 mm armor piercing ammunition.In this methodology,the FEM simulations are used to create training cases for Multilayer Perceptron(MLP) three layer networks.In order to validate FE simulation methodology,ballistic shot tests on 20 mm thickness target were performed according to standard Stanag 4569.Afterwards,the successfully trained ANN(s) is used to predict the ballistic limit thickness of 500 HB high hardness steel armor.Results show that even with limited number of data,FEM-ANN approach can be used to predict ballistic penetration depth with adequate accuracy.展开更多
We review recent advances in the finite element method (FEM) simulations of interactions between waves and structures. Our focus is on the potential theory with the fully nonlinear or second-order boundary condition. ...We review recent advances in the finite element method (FEM) simulations of interactions between waves and structures. Our focus is on the potential theory with the fully nonlinear or second-order boundary condition. The present paper has six sections. A review of previous work on interactions between waves and ocean structures is presented in Section one. Section two gives the mathematical formulation. In Section three, the finite element discretization, mesh generation and the finite element linear system solution methods are described. Section four presents numerical methods including time marching schemes, computation of velocity, remeshing and smoothing techniques and numerical radiation conditions. The application of the FEM to the wave-structure interactions are presented in Section five followed by the concluding remarks in Section six.展开更多
Base d on fluid velocity potential, an ALE finite element formulation for the analysi s of nonlinear sloshing problems has been developed. The ALE kinemat ical description is introduced to move the computational mesh...Base d on fluid velocity potential, an ALE finite element formulation for the analysi s of nonlinear sloshing problems has been developed. The ALE kinemat ical description is introduced to move the computational mesh independently of f luid motion, and the container fixed noninertial coordinate system is employed to establish the governing equations so that the mesh is needed to be updated in this coordinate system only. This leads to a very simple mesh moving algorithm which makes it easy to trace the motion of the moving boundaries and the free su rface without producing undesirable distortion of the computational mesh. The fi nite element method and finite difference method are used spacewise and timewise , respectively. A numerical example involving either forced horizontal oscillati on or forced pitching oscillation of the fluid filled container is presented to illustrate the effectiveness and the robustness of the method. In additi on, this work can be extended for the fluid structure interaction problems.展开更多
A contact bolt model is proposed as a new modeling technique to investigate the complex structure with bolted joints for modal analysis and compared with the coupled bolt model, and the test results are given. Among t...A contact bolt model is proposed as a new modeling technique to investigate the complex structure with bolted joints for modal analysis and compared with the coupled bolt model, and the test results are given. Among these models, the coupled bolt model provides the best accurate responses compared with the experimental results. The contact bolt model shows the best effectiveness and usefulness in view of operational time. The bolt models proposed in this study are adopted for a dynamic characteristic analysis of a large diesel engine consisting of several parts which are connected by many bolts. The dynamic behavior of the entire engine structure was investigated by experiment. The coupled bolt model and the contact bolt model were applied to model the assembly of engine with high preload. The experimental results are in good agreement with the finite element method （FEM） results. Compared with the other models, the contact bolt model presented in this paper is more effective and useful in view of operational time and experience of analysts.展开更多
We present a novel method for the solution of the diffusion equation on a composite AMR mesh. This approach is suitable for including diffusion based physics modules to hydrocodes that support ALE and AMR capabilities...We present a novel method for the solution of the diffusion equation on a composite AMR mesh. This approach is suitable for including diffusion based physics modules to hydrocodes that support ALE and AMR capabilities. To illustrate, we proffer our implementations of diffu- sion based radiation transport and heat conduction in a hydrocode called ALE-AMR. Numerical experiments conducted with the diffusion solver and associated physics packages yield 2nd order convergence in the L2 norm.展开更多
Considering the stochastic spatial variation of geotechnical parameters over the slope, a Stochastic Finite Element Method (SFEM) is established based on the combination of the Shear Strength Reduction (SSR) concept a...Considering the stochastic spatial variation of geotechnical parameters over the slope, a Stochastic Finite Element Method (SFEM) is established based on the combination of the Shear Strength Reduction (SSR) concept and quasi-Monte Carlo simulation. The shear strength reduction FEM is superior to the slice method based on the limit equilibrium theory in many ways, so it will be more powerful to assess the reliability of global slope stability when combined with probability theory. To illustrate the performance of the proposed method, it is applied to an example of simple slope. The results of simulation show that the proposed method is effective to perform the reliability analysis of global slope stability without presupposing a potential slip surface.展开更多
In order to study the laws of the extrusion pressure changing with the extrusion parameters in the process of hydrostatic extrusion for the tungsten alloys, the large deformation elasto plastic theory and the sof...In order to study the laws of the extrusion pressure changing with the extrusion parameters in the process of hydrostatic extrusion for the tungsten alloys, the large deformation elasto plastic theory and the software of ANSYS 5 5 are used to carry out the numerical simulation research. The laws of the extrusion pressure changing with the extrusion parameters, such as the die angle, extrusion ratio, and friction coefficient, are obtained. The simulation results are in good agreement with the experimental ones, and the simulated results are believable.展开更多
Numerical simulations are performed on the interface with large deformation induced by the interaction between a moving shock and two consecutive bubbles. The high performance of the level set method for multi-materia...Numerical simulations are performed on the interface with large deformation induced by the interaction between a moving shock and two consecutive bubbles. The high performance of the level set method for multi-material interfaces is demonstrated. Discontinuous Galerkin finite element method is used to solve Euleri- an equations. And the fifth-order weighted essentially non-oscillatory （WENO） scheme is used to solve the level set equation for capturing multi-material interfaces. The ghost fluid method is used to deal with the interfacial boundary condition. Results are obtained for two bubble interacting with a moving shock. The contours of the constant density and the pressure at different time are given. In the computational domain, three different cases are considered, i.e. two helium bubbles, a helium bubble followed by an R22 bubble in the direction of the moving shock, and an R22 bubble followed by a helium bubble. Computational results indicate that multi-mate- rial interfaces can be properly captured by the level set method. Therefore, for problems involving the flow of three different materials with two different interfaces, each interface separating two different materials can be similarly handled.展开更多
Flat metal strips are deformed progressively into profiles with certain geometry sections by a series of successive rotating rolls at the room temperature, and it is called the cold roll-forming. An effective method i...Flat metal strips are deformed progressively into profiles with certain geometry sections by a series of successive rotating rolls at the room temperature, and it is called the cold roll-forming. An effective method is given for roll design to avoid major defects via finite element method （FEM） simulation using a booting model. The simulation gives reasonable fit to the actual product in the two major roll-forming defects, i.e. , edge wave and springback. The redesigned rolls using the approach of multi-stand FEM simulation can effectively control these two defects.展开更多
For studying the dynamic performance of subsea umbilical cable laying system and achieving the goal of cable tension and laying speed control,the rigid finite element method is used to discrete and transform the syste...For studying the dynamic performance of subsea umbilical cable laying system and achieving the goal of cable tension and laying speed control,the rigid finite element method is used to discrete and transform the system into a rigid-flexible coupling multi-body system which consists of rigid elements and spring-damping elements.The mathematical model of subsea umbilical cable laying system kinematic chain is presented with the second order Lagrange equation in the joint coordinate system,and dynamic modeling and simulation is performed with ADAMS.The dynamic analysis is conducted assuming the following three statuses:ideal laying,practical laying under wave disturbance,and practical laying with tension compensation.Results show that motion disturbances of the laying budge under sea waves,especially with heaving and pitching,will cause relatively serious fluctuations in cable tension and laying speed.Tension compensation,i.e.,active back tension torque control can restrict continuous tension increasing or decreasing effectively and rapidly,thus avoiding cable breach or buckling.展开更多
论文建立了基于比例边界有限元法(scaled boundary finite element methods,SBFEM)框架的非局部宏微观损伤模型,考虑材料细观物理参数的空间变异性,探讨了材料参数的空间变异性对结构开裂过程的影响。结果表明:考虑材料参数空间变异性后...论文建立了基于比例边界有限元法(scaled boundary finite element methods,SBFEM)框架的非局部宏微观损伤模型,考虑材料细观物理参数的空间变异性,探讨了材料参数的空间变异性对结构开裂过程的影响。结果表明:考虑材料参数空间变异性后,裂纹扩展路径具有不确定性,建议的模型能够很好地反应材料内在的随机性;随着结构受力情况的复杂化和结构本体缺陷的增多,裂纹开裂模式的变异性也会增大。自相关长度和参数变异系数对结构开裂分析结果有重要影响。展开更多
针对二维介质目标的电磁成像问题,将正余弦算法(Sine Cosine Algorithm,SCA)与有限元方法(Finite Element Method,FEM)和不变性测试方程(Measured Equation of Invariance,MEI)进行结合提出一种新的成像方法。将FEM与MEI进行结合求解二...针对二维介质目标的电磁成像问题,将正余弦算法(Sine Cosine Algorithm,SCA)与有限元方法(Finite Element Method,FEM)和不变性测试方程(Measured Equation of Invariance,MEI)进行结合提出一种新的成像方法。将FEM与MEI进行结合求解二维介质目标的电磁散射正问题,即求解Helmholtz方程。其中,MEI保证边界截断的精度,FEM适用于复杂介质目标的准确模拟。对于电磁散射逆问题,引入SCA并加以改进提出一种新的重构方法。该方法采用等效原理与格林函数的渐近式求得远区散射场,以测量的散射场和计算的散射场最大偏差为目标函数,采用改进的SCA优化介质参数,使目标函数达到最小值,以此重构散射体。为提高计算效率,采用MPI算法进行并行计算。文中采用基准函数展示了改进的SCA算法的快速收敛性,并采用非规则的均匀介质柱目标验证了成像方法的正确性。展开更多
基金National Natural Science Foundation of China(No.51674280,51774308,51704033,51722406,51950410591)Shandong Provincial Natural Science Foundation(ZR2019JQ21,JQ201808)+3 种基金the Fundamental Research Funds for the Central Universities(No.20CX02113A)National Science and Technology Major Project(2016ZX05014-000407)Program for Changjiang Scholars and Innovative Research Team in University(IRT_16R69)PetroChina Innovation Foundation(No.2018D-5007-0210)。
文摘Due to the low permeability of tight reservoirs,throats play a significant role in controlling fluid flow.Although many studies have been conducted to investigate fluid flow in throats in the microscale domain,comparatively fewer works have been devoted to study the effect of adsorption boundary layer(ABL)in throats based on the digital rock method.By considering an ABL,we investigate its effects on fluid flow.We build digital rock model based on computed tomography technology.Then,microscopic pore structures are extracted with watershed segmentation and pore geometries are meshed through Delaunay triangulation approach.Finally,using the meshed digital simulation model and finite element method,we investigate the effects of viscosity and thickness of ABL on microscale flow.Our results demonstrate that viscosity and thickness of ABL are major factors that significantly hinder fluid flow in throats.
基金Otokar Otomotiv ve Savunma Sanayi A.S. for the financial support
文摘Determination of ballistic performance of an armor solution is a complicated task and evolved significantly with the application of finite element methods(FEM) in this research field.The traditional armor design studies performed with FEM requires sophisticated procedures and intensive computational effort,therefore simpler and accurate numerical approaches are always worthwhile to decrease armor development time.This study aims to apply a hybrid method using FEM simulation and artificial neural network(ANN) analysis to approximate ballistic limit thickness for armor steels.To achieve this objective,a predictive model based on the artificial neural networks is developed to determine ballistic resistance of high hardness armor steels against 7.62 mm armor piercing ammunition.In this methodology,the FEM simulations are used to create training cases for Multilayer Perceptron(MLP) three layer networks.In order to validate FE simulation methodology,ballistic shot tests on 20 mm thickness target were performed according to standard Stanag 4569.Afterwards,the successfully trained ANN(s) is used to predict the ballistic limit thickness of 500 HB high hardness steel armor.Results show that even with limited number of data,FEM-ANN approach can be used to predict ballistic penetration depth with adequate accuracy.
文摘We review recent advances in the finite element method (FEM) simulations of interactions between waves and structures. Our focus is on the potential theory with the fully nonlinear or second-order boundary condition. The present paper has six sections. A review of previous work on interactions between waves and ocean structures is presented in Section one. Section two gives the mathematical formulation. In Section three, the finite element discretization, mesh generation and the finite element linear system solution methods are described. Section four presents numerical methods including time marching schemes, computation of velocity, remeshing and smoothing techniques and numerical radiation conditions. The application of the FEM to the wave-structure interactions are presented in Section five followed by the concluding remarks in Section six.
文摘Base d on fluid velocity potential, an ALE finite element formulation for the analysi s of nonlinear sloshing problems has been developed. The ALE kinemat ical description is introduced to move the computational mesh independently of f luid motion, and the container fixed noninertial coordinate system is employed to establish the governing equations so that the mesh is needed to be updated in this coordinate system only. This leads to a very simple mesh moving algorithm which makes it easy to trace the motion of the moving boundaries and the free su rface without producing undesirable distortion of the computational mesh. The fi nite element method and finite difference method are used spacewise and timewise , respectively. A numerical example involving either forced horizontal oscillati on or forced pitching oscillation of the fluid filled container is presented to illustrate the effectiveness and the robustness of the method. In additi on, this work can be extended for the fluid structure interaction problems.
基金Sponsored by the Ministerial Level Foundation(40402020105)
文摘A contact bolt model is proposed as a new modeling technique to investigate the complex structure with bolted joints for modal analysis and compared with the coupled bolt model, and the test results are given. Among these models, the coupled bolt model provides the best accurate responses compared with the experimental results. The contact bolt model shows the best effectiveness and usefulness in view of operational time. The bolt models proposed in this study are adopted for a dynamic characteristic analysis of a large diesel engine consisting of several parts which are connected by many bolts. The dynamic behavior of the entire engine structure was investigated by experiment. The coupled bolt model and the contact bolt model were applied to model the assembly of engine with high preload. The experimental results are in good agreement with the finite element method （FEM） results. Compared with the other models, the contact bolt model presented in this paper is more effective and useful in view of operational time and experience of analysts.
基金the auspices of the U.S.Department of Energy by Lawrence Livermore National Laboratory under Contract DEAC52-07NA27344LBNL under DE-AC0205CH11231 was supported by the Director,Office ofScience of the U.S.Department of Energy and the Petascale Initiative in Computational Science and Engineeringthe National Energy Research Scientific Computing Center,supported by the Office of Science,U.S.Department of Energy under Contract No.DE-AC02-05CH11231.
文摘We present a novel method for the solution of the diffusion equation on a composite AMR mesh. This approach is suitable for including diffusion based physics modules to hydrocodes that support ALE and AMR capabilities. To illustrate, we proffer our implementations of diffu- sion based radiation transport and heat conduction in a hydrocode called ALE-AMR. Numerical experiments conducted with the diffusion solver and associated physics packages yield 2nd order convergence in the L2 norm.
文摘Considering the stochastic spatial variation of geotechnical parameters over the slope, a Stochastic Finite Element Method (SFEM) is established based on the combination of the Shear Strength Reduction (SSR) concept and quasi-Monte Carlo simulation. The shear strength reduction FEM is superior to the slice method based on the limit equilibrium theory in many ways, so it will be more powerful to assess the reliability of global slope stability when combined with probability theory. To illustrate the performance of the proposed method, it is applied to an example of simple slope. The results of simulation show that the proposed method is effective to perform the reliability analysis of global slope stability without presupposing a potential slip surface.
文摘In order to study the laws of the extrusion pressure changing with the extrusion parameters in the process of hydrostatic extrusion for the tungsten alloys, the large deformation elasto plastic theory and the software of ANSYS 5 5 are used to carry out the numerical simulation research. The laws of the extrusion pressure changing with the extrusion parameters, such as the die angle, extrusion ratio, and friction coefficient, are obtained. The simulation results are in good agreement with the experimental ones, and the simulated results are believable.
基金Supported by the National Natural Science Foundation of China(10476011)~~
文摘Numerical simulations are performed on the interface with large deformation induced by the interaction between a moving shock and two consecutive bubbles. The high performance of the level set method for multi-material interfaces is demonstrated. Discontinuous Galerkin finite element method is used to solve Euleri- an equations. And the fifth-order weighted essentially non-oscillatory （WENO） scheme is used to solve the level set equation for capturing multi-material interfaces. The ghost fluid method is used to deal with the interfacial boundary condition. Results are obtained for two bubble interacting with a moving shock. The contours of the constant density and the pressure at different time are given. In the computational domain, three different cases are considered, i.e. two helium bubbles, a helium bubble followed by an R22 bubble in the direction of the moving shock, and an R22 bubble followed by a helium bubble. Computational results indicate that multi-mate- rial interfaces can be properly captured by the level set method. Therefore, for problems involving the flow of three different materials with two different interfaces, each interface separating two different materials can be similarly handled.
文摘Flat metal strips are deformed progressively into profiles with certain geometry sections by a series of successive rotating rolls at the room temperature, and it is called the cold roll-forming. An effective method is given for roll design to avoid major defects via finite element method （FEM） simulation using a booting model. The simulation gives reasonable fit to the actual product in the two major roll-forming defects, i.e. , edge wave and springback. The redesigned rolls using the approach of multi-stand FEM simulation can effectively control these two defects.
基金Supported by the National Science and Technology Major Project:Development of Large Oil-gas Fields&Coal-bed Methane(No.2011ZX05056)
文摘For studying the dynamic performance of subsea umbilical cable laying system and achieving the goal of cable tension and laying speed control,the rigid finite element method is used to discrete and transform the system into a rigid-flexible coupling multi-body system which consists of rigid elements and spring-damping elements.The mathematical model of subsea umbilical cable laying system kinematic chain is presented with the second order Lagrange equation in the joint coordinate system,and dynamic modeling and simulation is performed with ADAMS.The dynamic analysis is conducted assuming the following three statuses:ideal laying,practical laying under wave disturbance,and practical laying with tension compensation.Results show that motion disturbances of the laying budge under sea waves,especially with heaving and pitching,will cause relatively serious fluctuations in cable tension and laying speed.Tension compensation,i.e.,active back tension torque control can restrict continuous tension increasing or decreasing effectively and rapidly,thus avoiding cable breach or buckling.
文摘论文建立了基于比例边界有限元法(scaled boundary finite element methods,SBFEM)框架的非局部宏微观损伤模型,考虑材料细观物理参数的空间变异性,探讨了材料参数的空间变异性对结构开裂过程的影响。结果表明:考虑材料参数空间变异性后,裂纹扩展路径具有不确定性,建议的模型能够很好地反应材料内在的随机性;随着结构受力情况的复杂化和结构本体缺陷的增多,裂纹开裂模式的变异性也会增大。自相关长度和参数变异系数对结构开裂分析结果有重要影响。
文摘针对二维介质目标的电磁成像问题,将正余弦算法(Sine Cosine Algorithm,SCA)与有限元方法(Finite Element Method,FEM)和不变性测试方程(Measured Equation of Invariance,MEI)进行结合提出一种新的成像方法。将FEM与MEI进行结合求解二维介质目标的电磁散射正问题,即求解Helmholtz方程。其中,MEI保证边界截断的精度,FEM适用于复杂介质目标的准确模拟。对于电磁散射逆问题,引入SCA并加以改进提出一种新的重构方法。该方法采用等效原理与格林函数的渐近式求得远区散射场,以测量的散射场和计算的散射场最大偏差为目标函数,采用改进的SCA优化介质参数,使目标函数达到最小值,以此重构散射体。为提高计算效率,采用MPI算法进行并行计算。文中采用基准函数展示了改进的SCA算法的快速收敛性,并采用非规则的均匀介质柱目标验证了成像方法的正确性。
