The electromagnetic(EM)telemetry systems,employed for real-time data transmission from the borehole and the earth surface during drilling,are widely used in measurement-while-drilling(MWD)and logging-while-drilling(LW...The electromagnetic(EM)telemetry systems,employed for real-time data transmission from the borehole and the earth surface during drilling,are widely used in measurement-while-drilling(MWD)and logging-while-drilling(LWD).Several numerical methods,including the method of moments(MoM),the electric field integral equation(EFIE)method,and the finite-element(FE)method have been developed for the simulation of EM telemetry systems.The computational process of these methods is complicated and time-consuming.To solve this problem,we introduce an axisymmetric semi-analytical FE method(SAFEM)in the cylindrical coordinate system with the virtual layering technique for rapid simulation of EM telemetry in a layered earth.The proposed method divides the computational domain into a series of homogeneous layers.For each layer,only its cross-section is discretized,and a high-precision integration method based on Riccati equations is employed for the calculation of longitudinally homogeneous sections.The block-tridiagonal structure of the global coefficient matrix enables the use of the block Thomas algorithm,facilitating the efficient simulation of EM telemetry problems in layered media.After the theoretical development,we validate the accuracy and efficiency of our algorithm through a series of numerical experiments and comparisons with the Multiphysics modeling software COMSOL.We also discussed the impact of system parameters on EM telemetry signal and demonstrated the applicability of our method by testing it on a field dataset acquired from Dezhou,Shandong Province,China.展开更多
In this paper, we propose a hybrid PML （H-PML） combining the normal absorption factor of convolutional PML （C-PML） with tangential absorption factor of Mutiaxial PML （M-PML）. The H-PML boundary conditions can be...In this paper, we propose a hybrid PML （H-PML） combining the normal absorption factor of convolutional PML （C-PML） with tangential absorption factor of Mutiaxial PML （M-PML）. The H-PML boundary conditions can better suppress the numerical instability in some extreme models, and the computational speed of finite-element method and the dynamic range are greatly increased using this HPML. We use the finite-element method with a hybrid PML to model the acoustic reflection of the interface when wireline and well logging while drilling （LWD）, in a formation with a reflector outside the borehole. The simulation results suggests that the PS- and SP- reflected waves arrive at the same time when the inclination between the well and the outer interface is zero, and the difference in arrival times increases with increasing dip angle. When there are fractures outside the well, the reflection signal is clearer in the subsequent reflection waves and may be used to identify the fractured zone. The difference between the dominant wavelength and the model scale shows that LWD reflection logging data are of higher resolution and quality than wireline acoustic reflection logging.展开更多
A least-squares finite-element method （LSFEM） for the non-conservative shallow-water equations is presented. The model is capable of handling complex topography, steady and unsteady flows, subcritical and supercriti...A least-squares finite-element method （LSFEM） for the non-conservative shallow-water equations is presented. The model is capable of handling complex topography, steady and unsteady flows, subcritical and supercritical flows, and flows with smooth and sharp gradient changes. Advantages of the model include： （1） sources terms, such as the bottom slope, surface stresses and bed frictions, can be treated easily without any special treatment; （2） upwind scheme is no needed; （3） a single approximating space can be used for all variables, and its choice of approximating space is not subject to the Ladyzhenskaya-Babuska-Brezzi （LBB） condition; and （4） the resulting system of equations is symmetric and positive-definite （SPD） which can be solved efficiently with the preconditioned conjugate gradient method. The model is verified with flow over a bump, tide induced flow, and dam-break. Computed results are compared with analytic solutions or other numerical results, and show the model is conservative and accurate. The model is then used to simulate flow past a circular cylinder. Important flow charac-teristics, such as variation of water surface around the cylinder and vortex shedding behind the cylinder are investigated. Computed results compare well with experiment data and other numerical results.展开更多
A modeling tool for simulating three-dimensional land frequency-domain controlled-source electromagnetic surveys,based on a finite-element discretization of the Helmholtz equation for the electric fields,has been deve...A modeling tool for simulating three-dimensional land frequency-domain controlled-source electromagnetic surveys,based on a finite-element discretization of the Helmholtz equation for the electric fields,has been developed.The main difference between our modeling method and those previous works is edge finite-element approach applied to solving the three-dimensional land frequency-domain electromagnetic responses generated by horizontal electric dipole source.Firstly,the edge finite-element equation is formulated through the Galerkin method based on Helmholtz equation of the electric fields.Secondly,in order to check the validity of the modeling code,the numerical results are compared with the analytical solutions for a homogeneous half-space model.Finally,other three models are simulated with three-dimensional electromagnetic responses.The results indicate that the method can be applied for solving three-dimensional electromagnetic responses.The algorithm has been demonstrated,which can be effective to modeling the complex geo-electrical structures.This efficient algorithm will help to study the distribution laws of3-D land frequency-domain controlled-source electromagnetic responses and to setup basis for research of three-dimensional inversion.展开更多
Numerical solution of shallow-water equations （SWE） has been a challenging task because of its nonlinear hyperbolic nature, admitting discontinuous solution, and the need to satisfy the C-property. The presence of s...Numerical solution of shallow-water equations （SWE） has been a challenging task because of its nonlinear hyperbolic nature, admitting discontinuous solution, and the need to satisfy the C-property. The presence of source terms in momentum equations, such as the bottom slope and friction of bed, compounds the difficulties further. In this paper, a least-squares finite-element method for the space discretization and θ-method for the time integration is developed for the 2D non-conservative SWE including the source terms. Advantages of the method include： the source terms can be approximated easily with interpolation functions, no upwind scheme is needed, as well as the resulting system equations is symmetric and positive-definite, therefore, can be solved efficiently with the conjugate gradient method. The method is applied to steady and unsteady flows, subcritical and transcritical flow over a bump, 1D and 2D circular dam-break, wave past a circular cylinder, as well as wave past a hump. Computed results show good C-property, conservation property and compare well with exact solutions and other numerical results for flows with weak and mild gradient changes, but lead to inaccurate predictions for flows with strong gradient changes and discontinuities.展开更多
A new method regarding mesomechanics finite-element research is proposed to predict the peak shear strength of mudded intercalation materials on a mesoscopic scale. Based on geometric and mechanical parameters, along ...A new method regarding mesomechanics finite-element research is proposed to predict the peak shear strength of mudded intercalation materials on a mesoscopic scale. Based on geometric and mechanical parameters, along with the strain failure criteria obtained by sample＇s deformation characteristics, uniaxial compression tests on the sample were simulated through a finite-element model, which yielded values consistent with the data from the laboratory uniaxial compression tests, implying that the method is reasonable. Based on this model, a shear test was performed to calculate the peak shear strength of the mudded intercalation, consistent with values reported in the literature, thereby providing a new approach for investigating the mechanical properties of mudded intercalation materials.展开更多
For two-dimension nonlinear convection diffusion equation, a two-grid method of characteristics finite-element solution was constructed. In this method the nonlinear iterations is only to execute on the coarse grid an...For two-dimension nonlinear convection diffusion equation, a two-grid method of characteristics finite-element solution was constructed. In this method the nonlinear iterations is only to execute on the coarse grid and the fine-grid solution can be obtained in a single linear step. For the nonlinear convection-dominated diffusion equation, this method can not only stabilize the numerical oscillation but also accelerate the convergence and improve the computational efficiency. The error analysis demonstrates if the mesh sizes between coarse-grid and fine-grid satisfy the certain relationship, the two-grid solution and the characteristics finite-element solution have the same order of accuracy. The numerical is more efficient than that of characteristics example confirms that the two-grid method finite-element method.展开更多
A new finite element method （FEM） of B-spline wavelet on the interval （BSWI） is proposed. Through analyzing the scaling functions of BSWI in one dimension, the basic formula for 2D FEM of BSWI is deduced. The 2D F...A new finite element method （FEM） of B-spline wavelet on the interval （BSWI） is proposed. Through analyzing the scaling functions of BSWI in one dimension, the basic formula for 2D FEM of BSWI is deduced. The 2D FEM of 7 nodes and 10 nodes are constructed based on the basic formula. Using these proposed elements, the multiscale numerical model for foundation subjected to harmonic periodic load, the foundation model excited by external and internal dynamic load are studied. The results show the pro- posed finite elements have higher precision than the tradi- tional elements with 4 nodes. The proposed finite elements can describe the propagation of stress waves well whenever the foundation model excited by extemal or intemal dynamic load. The proposed finite elements can be also used to con- nect the multi-scale elements. And the proposed finite elements also have high precision to make multi-scale analysis for structure.展开更多
In this paper,we introduce new stable mixed finite elements of any order on polytopal mesh for solving second-order elliptic problem.We establish optimal order error estimates for velocity and super convergence for pr...In this paper,we introduce new stable mixed finite elements of any order on polytopal mesh for solving second-order elliptic problem.We establish optimal order error estimates for velocity and super convergence for pressure.Numerical experiments are conducted for our mixed elements of different orders on 2D and 3D spaces that confirm the theory.展开更多
Nowadays,an increasing number of ships and marine structures are manufactured and inevitably operated in rough sea.As a result,some phenomena related to the violent fluid-elastic structure interactions(e.g.,hydrodynam...Nowadays,an increasing number of ships and marine structures are manufactured and inevitably operated in rough sea.As a result,some phenomena related to the violent fluid-elastic structure interactions(e.g.,hydrodynamic slamming on marine vessels,tsunami impact on onshore structures,and sloshing in liquid containers)have aroused huge challenges to ocean engineering fields.In this paper,the moving particle semi-implicit(MPS)method and finite element method(FEM)coupled method is proposed for use in numerical investigations of the interaction between a regular wave and a horizontal suspended structure.The fluid domain calculated by the MPS method is dispersed into fluid particles,and the structure domain solved by the FEM method is dispersed into beam elements.The generation of the 2D regular wave is firstly conducted,and convergence verification is performed to determine appropriate particle spacing for the simulation.Next,the regular wave interacting with a rigid structure is initially performed and verified through the comparison with the laboratory experiments.By verification,the MPS-FEM coupled method can be applied to fluid-structure interaction(FSI)problems with waves.On this basis,taking the flexibility of structure into consideration,the elastic dynamic response of the structure subjected to the wave slamming is investigated,including the evolutions of the free surface,the variation of the wave impact pressures,the velocity distribution,and the structural deformation response.By comparison with the rigid case,the effects of the structural flexibility on wave-elastic structure interaction can be obtained.展开更多
In this study, the limit state equation for tensile reliability analysis of the foundation surface of a gravity dam was established. The possible crack length was set as the action effect and allowable crack length wa...In this study, the limit state equation for tensile reliability analysis of the foundation surface of a gravity dam was established. The possible crack length was set as the action effect and allowable crack length was set as the resistance in the limit state. The nonlinear FEM was used to obtain the crack length of the foundation surface of the gravity dam, and the linear response surface method based on the orthogonal test design method was used to calculate the reliability, providing a reasonable and simple method for calculating the reliability of the serviceability limit state. The Longtan RCC gravity dam was chosen as an example. An orthogonal test, including eleven factors and two levels, was conducted, and the tensile reliability was calculated. The analysis shows that this method is reasonable.展开更多
Pilot biomechanical design of biomaterials for artificial nucleus prosthesiswas carried out based on the 3D finite-element method. Two 3D models of lumbar intervertebral discrespectively with a real human nucleus and ...Pilot biomechanical design of biomaterials for artificial nucleus prosthesiswas carried out based on the 3D finite-element method. Two 3D models of lumbar intervertebral discrespectively with a real human nucleus and with the nucleus removed were developed and validatedusing published experimental and clinical data. Then the models with a stainless steel nucleusprosthesis implanted and with polymer nucleus prostheses of various properties implanted were usedfor the 3D finite-element biomechanical analysis. All the above simulation and analysis were carriedout for the L4/L5 disc under a human worst--daily compression load of 2000 N. The results show thatthe polymer materials with Young's modulus of elasticity E = 0.1-100 MPa and Poisson's ratio v=0.35-0.5 are suitable to produce artificial nucleus prosthesis in view of biomechanicalconsideration.展开更多
The energy produced by the melting stretching disks surface has a wide range of commercial applications,including semi-conductor material preparation,magma solidification,permafrost melting,and frozen land refreezing,...The energy produced by the melting stretching disks surface has a wide range of commercial applications,including semi-conductor material preparation,magma solidification,permafrost melting,and frozen land refreezing,among others.In view of this,in the current communication we analyzed magnetohydrodynamic flow ofMaxwell nanofluid between two parallel rotating disks.Nanofluids are important due to their astonishing properties in heat conduction flows and in the enhancement of electronic and manufacturing devices.Furthermore,the distinct tinysized particles Al_(2)O_(3)and TiO_(2)in theMaxwell water-based fluid for enhancing the heat transfer rate are analyzed.The heat equation is developed in the occurrence of thermal radiation.The influences of melting impacts are incorporated.The mathematical model is developed in the form of partial differential expressions then converted to ordinary differential equations by employing tool of similarity variables.Finite element method(FEM)is chosen for solving the nonlinear governing ordinary differential equations(ODEs)with necessary conditions.The consequence of flow parameters against the velocity profiles and heat transport field is considered.The noted novelty of this communication is to discuss the thermal transfer of Maxwell nanofluid model through double stretching disks with thermal radiation and melting phenomenon.Further,Al_(2)O_(3)/water and TiO_(2)/water are considered in the modeling.展开更多
The structural de formation of Lu’ an mining area is characterized by a remarkable feature of zoning along E-W direction, in the east limb of Qinshui basin, Shanxi Province, China. The re gional tectonic stress field...The structural de formation of Lu’ an mining area is characterized by a remarkable feature of zoning along E-W direction, in the east limb of Qinshui basin, Shanxi Province, China. The re gional tectonic stress fields and basement tectonics are two fundamental factors to control the cover tectonic framework. This paper uses the finite-element method with a elastic-plastic pIan problem model to simulate the three periods of stress fields resulting from field geological study’ Based on these works, the formation and evolution of tectonic framework of Lu’ an mining area have been discussed.展开更多
Fluid-structure interaction (FSI) problems caused by fluid impact loads are com- monly existent in naval architectures and ocean engineering fields. For instance, the impact loads due to non-linear fluid motion in a l...Fluid-structure interaction (FSI) problems caused by fluid impact loads are com- monly existent in naval architectures and ocean engineering fields. For instance, the impact loads due to non-linear fluid motion in a liquid sloshing tank potentially affect the structural safety of cargo tanks or vessels. The challenges of numerical study on FSI problems involve not only multidisciplinary features, but also accurate description of non-linear free surface. A fully Lagrangian particle-based method , the moving particle semi-implicit and finite element coupled method ( MPS-FEM), is developed to numerically study the FSI problems. Taking into account the advantage of the Lagrangian method for large deformations of both fluid and solid bounda- ties, the MPS method is used to simulate the fluid field while the finite element method(FEM) to calculate the structure field. Besides, the partitioning strategy is employed to couple the MPS and FEM modules. To validate accuracy of the proposed algorithm, a benchmark case is numer- ically investigated. Both the patterns of free surface and the deflections of the elastic structures are in good agreement with the experimental data. Then, the present FSI solver is applied to the comparative study of the mitigating effects of rigid baffles and elastic baffles on the sloshing motions and impact loads.展开更多
The numerical analytic research approach of stress-strain state of anisotropic composite finite element area with different boundary conditions on the surface, is represented below. The problem is solved by using a sp...The numerical analytic research approach of stress-strain state of anisotropic composite finite element area with different boundary conditions on the surface, is represented below. The problem is solved by using a spatial model of the elasticity theory. Differential equation system in partial derivatives reduces to one-dimensional problem using spline collocation method in two coordinate directions. Boundary problem for the system of ordinary higher-order differential equation is solved by using the stable numerical technique of discrete orthogonalization.展开更多
高位滑坡对建筑集群的冲击破坏时常导致严重的人员伤亡,基于光滑粒子流体动力学-离散元法-有限元法(smoothed particle hydrodynamics-discrete element method-finite element method,SPH-DEM-FEM)耦合的数值模型,开展了高位滑坡对框...高位滑坡对建筑集群的冲击破坏时常导致严重的人员伤亡,基于光滑粒子流体动力学-离散元法-有限元法(smoothed particle hydrodynamics-discrete element method-finite element method,SPH-DEM-FEM)耦合的数值模型,开展了高位滑坡对框架结构建筑群的冲击过程、建筑结构破坏机理、冲击力时程与框架柱关键点应力和弯矩等动力机制研究。研究结果表明:SPH-DEM-FEM耦合数值方法能够有效地模拟碎石土滑坡中土(SPH)石(DEM)混合物的抛射弹跳、爬高绕流冲击运动过程。考虑了常规建筑垂直、平行于滑坡流向的三排建筑组合布局,位于滑坡近端的纵向排列建筑表现为连续性倾倒破坏,横向排列的建筑则呈现整体倾倒破坏;因前排建筑群对滑坡冲击能量的耗散及滑坡自身摩擦耗能,位于滑坡后端建筑表现为引流面墙体和前排柱发生局部破坏,结构保持稳定,损毁程度依次为上游无建筑缓冲耗能的建筑>有横向排列的建筑>有纵向排列的建筑;纵向、横向排列的建筑冲击力衰减幅度分别31%、21%。横向框架建筑整体倾倒的损毁机制表现为框架柱的直接剪断或节点塑形铰链失效;纵向框架建筑连续性倾倒的损毁机制表现为前排框架柱的失效引起后排框架柱轴向压力和极限弯矩增加,持续冲击荷载超过其极限弯矩致使后排框架柱发生弯曲破坏,最终结构倾倒。系统能量在动能、内能和摩擦耗能间转化,其中摩擦耗能占65.5%,结构耗能占23.6%,动能快速下降与内能急剧增加是建筑破坏的关键特征。展开更多
基金supported by the Major Research Project on Scientific Instrument Development of the National Natural Science Foundation of China(42327901)National Natural Science Foundation of China(42030806,42074120,41904104,423B2405).
文摘The electromagnetic(EM)telemetry systems,employed for real-time data transmission from the borehole and the earth surface during drilling,are widely used in measurement-while-drilling(MWD)and logging-while-drilling(LWD).Several numerical methods,including the method of moments(MoM),the electric field integral equation(EFIE)method,and the finite-element(FE)method have been developed for the simulation of EM telemetry systems.The computational process of these methods is complicated and time-consuming.To solve this problem,we introduce an axisymmetric semi-analytical FE method(SAFEM)in the cylindrical coordinate system with the virtual layering technique for rapid simulation of EM telemetry in a layered earth.The proposed method divides the computational domain into a series of homogeneous layers.For each layer,only its cross-section is discretized,and a high-precision integration method based on Riccati equations is employed for the calculation of longitudinally homogeneous sections.The block-tridiagonal structure of the global coefficient matrix enables the use of the block Thomas algorithm,facilitating the efficient simulation of EM telemetry problems in layered media.After the theoretical development,we validate the accuracy and efficiency of our algorithm through a series of numerical experiments and comparisons with the Multiphysics modeling software COMSOL.We also discussed the impact of system parameters on EM telemetry signal and demonstrated the applicability of our method by testing it on a field dataset acquired from Dezhou,Shandong Province,China.
基金supported by the National Natural Science Foundation of China(No.41204094)Science Foundation of China University of Petroleum,Beijing(No.2462015YQ0506)
文摘In this paper, we propose a hybrid PML （H-PML） combining the normal absorption factor of convolutional PML （C-PML） with tangential absorption factor of Mutiaxial PML （M-PML）. The H-PML boundary conditions can better suppress the numerical instability in some extreme models, and the computational speed of finite-element method and the dynamic range are greatly increased using this HPML. We use the finite-element method with a hybrid PML to model the acoustic reflection of the interface when wireline and well logging while drilling （LWD）, in a formation with a reflector outside the borehole. The simulation results suggests that the PS- and SP- reflected waves arrive at the same time when the inclination between the well and the outer interface is zero, and the difference in arrival times increases with increasing dip angle. When there are fractures outside the well, the reflection signal is clearer in the subsequent reflection waves and may be used to identify the fractured zone. The difference between the dominant wavelength and the model scale shows that LWD reflection logging data are of higher resolution and quality than wireline acoustic reflection logging.
基金the National Science Council ot Taiwan,China for funding this research(Project no.:NSC 94-2218-E-035-011)
文摘A least-squares finite-element method （LSFEM） for the non-conservative shallow-water equations is presented. The model is capable of handling complex topography, steady and unsteady flows, subcritical and supercritical flows, and flows with smooth and sharp gradient changes. Advantages of the model include： （1） sources terms, such as the bottom slope, surface stresses and bed frictions, can be treated easily without any special treatment; （2） upwind scheme is no needed; （3） a single approximating space can be used for all variables, and its choice of approximating space is not subject to the Ladyzhenskaya-Babuska-Brezzi （LBB） condition; and （4） the resulting system of equations is symmetric and positive-definite （SPD） which can be solved efficiently with the preconditioned conjugate gradient method. The model is verified with flow over a bump, tide induced flow, and dam-break. Computed results are compared with analytic solutions or other numerical results, and show the model is conservative and accurate. The model is then used to simulate flow past a circular cylinder. Important flow charac-teristics, such as variation of water surface around the cylinder and vortex shedding behind the cylinder are investigated. Computed results compare well with experiment data and other numerical results.
基金Projects(41674080,41674079)supported by the National Natural Science Foundation of China
文摘A modeling tool for simulating three-dimensional land frequency-domain controlled-source electromagnetic surveys,based on a finite-element discretization of the Helmholtz equation for the electric fields,has been developed.The main difference between our modeling method and those previous works is edge finite-element approach applied to solving the three-dimensional land frequency-domain electromagnetic responses generated by horizontal electric dipole source.Firstly,the edge finite-element equation is formulated through the Galerkin method based on Helmholtz equation of the electric fields.Secondly,in order to check the validity of the modeling code,the numerical results are compared with the analytical solutions for a homogeneous half-space model.Finally,other three models are simulated with three-dimensional electromagnetic responses.The results indicate that the method can be applied for solving three-dimensional electromagnetic responses.The algorithm has been demonstrated,which can be effective to modeling the complex geo-electrical structures.This efficient algorithm will help to study the distribution laws of3-D land frequency-domain controlled-source electromagnetic responses and to setup basis for research of three-dimensional inversion.
基金the National Science Council of Taiwan for funding this research (NSC 96-2221-E-019-061).
文摘Numerical solution of shallow-water equations （SWE） has been a challenging task because of its nonlinear hyperbolic nature, admitting discontinuous solution, and the need to satisfy the C-property. The presence of source terms in momentum equations, such as the bottom slope and friction of bed, compounds the difficulties further. In this paper, a least-squares finite-element method for the space discretization and θ-method for the time integration is developed for the 2D non-conservative SWE including the source terms. Advantages of the method include： the source terms can be approximated easily with interpolation functions, no upwind scheme is needed, as well as the resulting system equations is symmetric and positive-definite, therefore, can be solved efficiently with the conjugate gradient method. The method is applied to steady and unsteady flows, subcritical and transcritical flow over a bump, 1D and 2D circular dam-break, wave past a circular cylinder, as well as wave past a hump. Computed results show good C-property, conservation property and compare well with exact solutions and other numerical results for flows with weak and mild gradient changes, but lead to inaccurate predictions for flows with strong gradient changes and discontinuities.
基金Funded by the National Natural Science Foundation of China(No.51574201)the State Key Laboratory of Geohazard Prevention and Geoenvironment Protection(Chengdu University of Technology)(KLGP2015K006)the Scientific and Technical Youth Innovation Group(Southwest Petroleum University)(2015CXTD05)
文摘A new method regarding mesomechanics finite-element research is proposed to predict the peak shear strength of mudded intercalation materials on a mesoscopic scale. Based on geometric and mechanical parameters, along with the strain failure criteria obtained by sample＇s deformation characteristics, uniaxial compression tests on the sample were simulated through a finite-element model, which yielded values consistent with the data from the laboratory uniaxial compression tests, implying that the method is reasonable. Based on this model, a shear test was performed to calculate the peak shear strength of the mudded intercalation, consistent with values reported in the literature, thereby providing a new approach for investigating the mechanical properties of mudded intercalation materials.
文摘For two-dimension nonlinear convection diffusion equation, a two-grid method of characteristics finite-element solution was constructed. In this method the nonlinear iterations is only to execute on the coarse grid and the fine-grid solution can be obtained in a single linear step. For the nonlinear convection-dominated diffusion equation, this method can not only stabilize the numerical oscillation but also accelerate the convergence and improve the computational efficiency. The error analysis demonstrates if the mesh sizes between coarse-grid and fine-grid satisfy the certain relationship, the two-grid solution and the characteristics finite-element solution have the same order of accuracy. The numerical is more efficient than that of characteristics example confirms that the two-grid method finite-element method.
基金supported by the National Natural Science Foundation of China (51109029,51178081,51138001,and 51009020)the State Key Development Program for Basic Research of China (2013CB035905)
文摘A new finite element method （FEM） of B-spline wavelet on the interval （BSWI） is proposed. Through analyzing the scaling functions of BSWI in one dimension, the basic formula for 2D FEM of BSWI is deduced. The 2D FEM of 7 nodes and 10 nodes are constructed based on the basic formula. Using these proposed elements, the multiscale numerical model for foundation subjected to harmonic periodic load, the foundation model excited by external and internal dynamic load are studied. The results show the pro- posed finite elements have higher precision than the tradi- tional elements with 4 nodes. The proposed finite elements can describe the propagation of stress waves well whenever the foundation model excited by extemal or intemal dynamic load. The proposed finite elements can be also used to con- nect the multi-scale elements. And the proposed finite elements also have high precision to make multi-scale analysis for structure.
基金supported in part by the National Science Foundation Grant DMS-1620016supported in parts by HKSAR grant Q81Q and JRI of The Hong Kong Polytechnic University.
文摘In this paper,we introduce new stable mixed finite elements of any order on polytopal mesh for solving second-order elliptic problem.We establish optimal order error estimates for velocity and super convergence for pressure.Numerical experiments are conducted for our mixed elements of different orders on 2D and 3D spaces that confirm the theory.
基金supported by the National Natural Science Foundation of China(51879159,51490675,11432009,and 51579145)Chang Jiang Scholars Program(T2014099)+3 种基金Shanghai Excellent Academic Leaders Program(17XD1402300)Program for Professor of Special Appointment(Eastern Scholar)at Shanghai Institutions of Higher Learning(2013022)Innovative Special Project of Numerical Tank of Ministry of Industry and Information Technology of China(2016-23/09)Lloyd’s Register Foundation for doctoral student
文摘Nowadays,an increasing number of ships and marine structures are manufactured and inevitably operated in rough sea.As a result,some phenomena related to the violent fluid-elastic structure interactions(e.g.,hydrodynamic slamming on marine vessels,tsunami impact on onshore structures,and sloshing in liquid containers)have aroused huge challenges to ocean engineering fields.In this paper,the moving particle semi-implicit(MPS)method and finite element method(FEM)coupled method is proposed for use in numerical investigations of the interaction between a regular wave and a horizontal suspended structure.The fluid domain calculated by the MPS method is dispersed into fluid particles,and the structure domain solved by the FEM method is dispersed into beam elements.The generation of the 2D regular wave is firstly conducted,and convergence verification is performed to determine appropriate particle spacing for the simulation.Next,the regular wave interacting with a rigid structure is initially performed and verified through the comparison with the laboratory experiments.By verification,the MPS-FEM coupled method can be applied to fluid-structure interaction(FSI)problems with waves.On this basis,taking the flexibility of structure into consideration,the elastic dynamic response of the structure subjected to the wave slamming is investigated,including the evolutions of the free surface,the variation of the wave impact pressures,the velocity distribution,and the structural deformation response.By comparison with the rigid case,the effects of the structural flexibility on wave-elastic structure interaction can be obtained.
文摘In this study, the limit state equation for tensile reliability analysis of the foundation surface of a gravity dam was established. The possible crack length was set as the action effect and allowable crack length was set as the resistance in the limit state. The nonlinear FEM was used to obtain the crack length of the foundation surface of the gravity dam, and the linear response surface method based on the orthogonal test design method was used to calculate the reliability, providing a reasonable and simple method for calculating the reliability of the serviceability limit state. The Longtan RCC gravity dam was chosen as an example. An orthogonal test, including eleven factors and two levels, was conducted, and the tensile reliability was calculated. The analysis shows that this method is reasonable.
文摘Pilot biomechanical design of biomaterials for artificial nucleus prosthesiswas carried out based on the 3D finite-element method. Two 3D models of lumbar intervertebral discrespectively with a real human nucleus and with the nucleus removed were developed and validatedusing published experimental and clinical data. Then the models with a stainless steel nucleusprosthesis implanted and with polymer nucleus prostheses of various properties implanted were usedfor the 3D finite-element biomechanical analysis. All the above simulation and analysis were carriedout for the L4/L5 disc under a human worst--daily compression load of 2000 N. The results show thatthe polymer materials with Young's modulus of elasticity E = 0.1-100 MPa and Poisson's ratio v=0.35-0.5 are suitable to produce artificial nucleus prosthesis in view of biomechanicalconsideration.
基金This work was sponsored in part by National Natural Science Foundation of China(No.51869031)Science and Technology Research Program of Chongqing Municipal Education Commission(Grant No.KJQN201903801)+1 种基金Huzhou Key Laboratory of Green Building TechnologyThis work is financially supported by the Government College University,Faisalabad and Higher Education Commission,Pakistan.
文摘The energy produced by the melting stretching disks surface has a wide range of commercial applications,including semi-conductor material preparation,magma solidification,permafrost melting,and frozen land refreezing,among others.In view of this,in the current communication we analyzed magnetohydrodynamic flow ofMaxwell nanofluid between two parallel rotating disks.Nanofluids are important due to their astonishing properties in heat conduction flows and in the enhancement of electronic and manufacturing devices.Furthermore,the distinct tinysized particles Al_(2)O_(3)and TiO_(2)in theMaxwell water-based fluid for enhancing the heat transfer rate are analyzed.The heat equation is developed in the occurrence of thermal radiation.The influences of melting impacts are incorporated.The mathematical model is developed in the form of partial differential expressions then converted to ordinary differential equations by employing tool of similarity variables.Finite element method(FEM)is chosen for solving the nonlinear governing ordinary differential equations(ODEs)with necessary conditions.The consequence of flow parameters against the velocity profiles and heat transport field is considered.The noted novelty of this communication is to discuss the thermal transfer of Maxwell nanofluid model through double stretching disks with thermal radiation and melting phenomenon.Further,Al_(2)O_(3)/water and TiO_(2)/water are considered in the modeling.
文摘The structural de formation of Lu’ an mining area is characterized by a remarkable feature of zoning along E-W direction, in the east limb of Qinshui basin, Shanxi Province, China. The re gional tectonic stress fields and basement tectonics are two fundamental factors to control the cover tectonic framework. This paper uses the finite-element method with a elastic-plastic pIan problem model to simulate the three periods of stress fields resulting from field geological study’ Based on these works, the formation and evolution of tectonic framework of Lu’ an mining area have been discussed.
文摘Fluid-structure interaction (FSI) problems caused by fluid impact loads are com- monly existent in naval architectures and ocean engineering fields. For instance, the impact loads due to non-linear fluid motion in a liquid sloshing tank potentially affect the structural safety of cargo tanks or vessels. The challenges of numerical study on FSI problems involve not only multidisciplinary features, but also accurate description of non-linear free surface. A fully Lagrangian particle-based method , the moving particle semi-implicit and finite element coupled method ( MPS-FEM), is developed to numerically study the FSI problems. Taking into account the advantage of the Lagrangian method for large deformations of both fluid and solid bounda- ties, the MPS method is used to simulate the fluid field while the finite element method(FEM) to calculate the structure field. Besides, the partitioning strategy is employed to couple the MPS and FEM modules. To validate accuracy of the proposed algorithm, a benchmark case is numer- ically investigated. Both the patterns of free surface and the deflections of the elastic structures are in good agreement with the experimental data. Then, the present FSI solver is applied to the comparative study of the mitigating effects of rigid baffles and elastic baffles on the sloshing motions and impact loads.
文摘The numerical analytic research approach of stress-strain state of anisotropic composite finite element area with different boundary conditions on the surface, is represented below. The problem is solved by using a spatial model of the elasticity theory. Differential equation system in partial derivatives reduces to one-dimensional problem using spline collocation method in two coordinate directions. Boundary problem for the system of ordinary higher-order differential equation is solved by using the stable numerical technique of discrete orthogonalization.
文摘高位滑坡对建筑集群的冲击破坏时常导致严重的人员伤亡,基于光滑粒子流体动力学-离散元法-有限元法(smoothed particle hydrodynamics-discrete element method-finite element method,SPH-DEM-FEM)耦合的数值模型,开展了高位滑坡对框架结构建筑群的冲击过程、建筑结构破坏机理、冲击力时程与框架柱关键点应力和弯矩等动力机制研究。研究结果表明:SPH-DEM-FEM耦合数值方法能够有效地模拟碎石土滑坡中土(SPH)石(DEM)混合物的抛射弹跳、爬高绕流冲击运动过程。考虑了常规建筑垂直、平行于滑坡流向的三排建筑组合布局,位于滑坡近端的纵向排列建筑表现为连续性倾倒破坏,横向排列的建筑则呈现整体倾倒破坏;因前排建筑群对滑坡冲击能量的耗散及滑坡自身摩擦耗能,位于滑坡后端建筑表现为引流面墙体和前排柱发生局部破坏,结构保持稳定,损毁程度依次为上游无建筑缓冲耗能的建筑>有横向排列的建筑>有纵向排列的建筑;纵向、横向排列的建筑冲击力衰减幅度分别31%、21%。横向框架建筑整体倾倒的损毁机制表现为框架柱的直接剪断或节点塑形铰链失效;纵向框架建筑连续性倾倒的损毁机制表现为前排框架柱的失效引起后排框架柱轴向压力和极限弯矩增加,持续冲击荷载超过其极限弯矩致使后排框架柱发生弯曲破坏,最终结构倾倒。系统能量在动能、内能和摩擦耗能间转化,其中摩擦耗能占65.5%,结构耗能占23.6%,动能快速下降与内能急剧增加是建筑破坏的关键特征。
