Traditional rigid body limit equilibrium method （RBLEM） was adopted for the stability evaluation and analysis of rock slope under earthquake scenario. It is not able to provide the real stress distribution of the st...Traditional rigid body limit equilibrium method （RBLEM） was adopted for the stability evaluation and analysis of rock slope under earthquake scenario. It is not able to provide the real stress distribution of the structure, while the strength reduction method relies on the arbitrary decision on the failure criteria. The dynamic limit equilibrium solution was proposed for the stability analysis of sliding block based on 3-D multi-grid method, by incorporating implicit stepping integration FEM. There are two independent meshes created in the analysis： One original 3-D FEM mesh is for the simulation of target structure and provides the stress time-history, while the other surface grid is for the simulation of sliding surface and could be selected and designed freely. As long as the stress time-history of the geotechnical structure under earthquake scenario is obtained based on 3-D nonlinear dynamic FEM analysis, the time-history of the force on sliding surface could be derived by projecting the stress time-history from 3-D FEM mesh to surface grid. After that, the safety factor time-history of the sliding block will be determined through applying limit equilibrium method. With those information in place, the structure＇s aseismatic stability ean be further studied. The above theory and method were also applied to the aseismatic stability analysis of Dagangshan arch dam＇s right bank high slope and compared with the the result generated by Quasi-static method. The comparative analysis reveals that the method not only raises the FEM＇s capability in accurate simulation of complicated geologic structure, but also increases the flexibility and comprehensiveness of limit equilibrium method. This method is reliable and recommended for further application in other real geotechnical engineering.展开更多
A supported framework of a gyroscope＇s rotor is designed and the B-Spline wavelet finite element model of nonlinear supported magnetic field is worked out. A new finite element space is studied in which the scaling f...A supported framework of a gyroscope＇s rotor is designed and the B-Spline wavelet finite element model of nonlinear supported magnetic field is worked out. A new finite element space is studied in which the scaling function of the B-spline wavelet is considered as the shape function of a tetrahedton. The magnetic field is spited by an artificial absorbing body which used the condition of field radiating, so the solution is unique. The resolution is improved via the varying gradient of the B-spline function under the condition of unchanging gridding. So there are some advantages in dealing with the focus flux and a high varying gradient result from a nonlinear magnetic field. The result is more practical. Plots of flux and in the space is studied via simulating the supported system model. The results of the study are useful in the research of the supported magnetic system for the gyroscope rotor.展开更多
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.展开更多
Based on Timoshenko＇s beam theory and Vlasov＇s thin-walled member theory, a new model of spatial thin-walled beam element is developed for analyzing geometrical and physical nonlinearity, which incorporates an inter...Based on Timoshenko＇s beam theory and Vlasov＇s thin-walled member theory, a new model of spatial thin-walled beam element is developed for analyzing geometrical and physical nonlinearity, which incorporates an interior node and independent interpolations of bending angles and warp and takes diversified factors into consideration, such as traverse shear deformation, torsional shear deformation and their coupling, coupling of flexure and torsion, and the second shear stress. The geometrical nonlinear strain is formulated in updated Lagarange （UL） and the corresponding stiffness matrix is derived. The perfectly plastic model is used to account for physical nonlinearity, and the yield rule of von Mises and incremental relationship of Prandtle-Reuss are adopted. Elastoplastic stiffness matrix is obtained by numerical integration based on the finite segment method, and a finite element program is compiled. Numerical examples manifest that the proposed model is accurate and feasible in the analysis of thin-walled structures.展开更多
The element energy projection （EEP） method for computation of super- convergent resulting in a one-dimensional finite element method （FEM） is successfully used to self-adaptive FEM analysis of various linear probl...The element energy projection （EEP） method for computation of super- convergent resulting in a one-dimensional finite element method （FEM） is successfully used to self-adaptive FEM analysis of various linear problems, based on which this paper presents a substantial extension of the whole set of technology to nonlinear problems. The main idea behind the technology transfer from linear analysis to nonlinear analysis is to use Newton＇s method to linearize nonlinear problems into a series of linear problems so that the EEP formulation and the corresponding adaptive strategy can be directly used without the need for specific super-convergence formulation for nonlinear FEM. As a re- sult, a unified and general self-adaptive algorithm for nonlinear FEM analysis is formed. The proposed algorithm is found to be able to produce satisfactory finite element results with accuracy satisfying the user-preset error tolerances by maximum norm anywhere on the mesh. Taking the nonlinear ordinary differential equation （ODE） of second-order as the model problem, this paper describes the related fundamental idea, the imple- mentation strategy, and the computational algorithm. Representative numerical exam- ples are given to show the efficiency, stability, versatility, and reliability of the proposed approach.展开更多
The numerical mode of nonlinear wave transformation based on both the Laplace equation for water field and the Bernoulli equation for water surface is a kind of time-domain boundary problem with initial conditions. An...The numerical mode of nonlinear wave transformation based on both the Laplace equation for water field and the Bernoulli equation for water surface is a kind of time-domain boundary problem with initial conditions. And the basis for establishing the numerical mode of nonlinear wave in time domain is to trace the position of wave free surface and to calculate the instantaneous surface height and surface potential function. This paper firstly utilizes the ‘0-1' combined BEM to separate the boundary by means of discretization of Green's integral equation based on the Laplace equation, then separates the free surface of wave with FEM and derives the FEM equation of wave surface that satisfies the nonlinear boundary conditions. By jointly solving the above BEM and FEM equations, the wave potential and surface height could be obtained with iteration in time domain. Thus a new kind of nonlinear numerical mode is established for calculating wave transformation. The wave test in the numerical wave tank shows that the numerical simulation with this mode is of high accuracy.展开更多
In the paper, a weak coupling numerical model is developed for the study of the nonlinear dynamic interaction between water waves and permeable sandy seabed. The wave field solveris based on the VOF （Volume of Fluid...In the paper, a weak coupling numerical model is developed for the study of the nonlinear dynamic interaction between water waves and permeable sandy seabed. The wave field solveris based on the VOF （Volume of Fluid） method for continuity equation and the two-dimensional Reynolds Averaged Navier Stokes （RANS） equations with a k-ε closure. The free surface of cnoidal wave is traced through the PLIC-VOF （P/ecewise Linear/nterface Construction）. Blot＇s equations have been applied to solve the sandy seabed, and the u-p fmite dement formulations are derived by the application of the Galerkin weighted-residual procedure. The continuity of the pressure on the interface between fluid and porous medium domains is considered. Laboratory tests were performed to verify the proposed numerical model, and it is shown that the pore-water pressures and the wave heights computed by the VOF-FEM models are in good agreement with the experimental results. It is found that the proposed model is effective in predicting the seabed-nonlinear wave interaction and is able to handle the wave-breakwater-seabed interaction problem.展开更多
SUN Da-peng BAO Wei-bin, WU Hao and LI Yu-cheng （ In this paper the 0-1 combined BEM is adopted to subdivide the computational domain boundary, and to discretize the Green＇s integral expression based on Laplace equ...SUN Da-peng BAO Wei-bin, WU Hao and LI Yu-cheng （ In this paper the 0-1 combined BEM is adopted to subdivide the computational domain boundary, and to discretize the Green＇s integral expression based on Laplace equation. The FEM is used to subdivide the wave surface and deduce the surface equation which satisfies the nonlinear boundary conditions on the surface. The equations with potential function and wave surface height as an unknown quantity by application of Taylor expansion approach can be solved by iteration within the time step. In m-time iteration within the computational process of time step （n-1）At to nat, the results of the previous iteration are taken as the initial value of the two-order unknown terms in the present iteration. Thus, an improved tracking mode of nonlinear wave surface is estabIished, and numerical results of wave tank test indicate that this mode is improved obviously and is more precise than the previous numerical model which ignored the two-order unknown terms of wave surface location and velocity potential function in comparison with the theoretical values.展开更多
基金Project(2013-KY-2) supported by the State Key Laboratory of Hydroscience and Engineering of Hydroscience, ChinaProject(50925931)supported by the National Funds for Distinguished Young Scientists, China
文摘Traditional rigid body limit equilibrium method （RBLEM） was adopted for the stability evaluation and analysis of rock slope under earthquake scenario. It is not able to provide the real stress distribution of the structure, while the strength reduction method relies on the arbitrary decision on the failure criteria. The dynamic limit equilibrium solution was proposed for the stability analysis of sliding block based on 3-D multi-grid method, by incorporating implicit stepping integration FEM. There are two independent meshes created in the analysis： One original 3-D FEM mesh is for the simulation of target structure and provides the stress time-history, while the other surface grid is for the simulation of sliding surface and could be selected and designed freely. As long as the stress time-history of the geotechnical structure under earthquake scenario is obtained based on 3-D nonlinear dynamic FEM analysis, the time-history of the force on sliding surface could be derived by projecting the stress time-history from 3-D FEM mesh to surface grid. After that, the safety factor time-history of the sliding block will be determined through applying limit equilibrium method. With those information in place, the structure＇s aseismatic stability ean be further studied. The above theory and method were also applied to the aseismatic stability analysis of Dagangshan arch dam＇s right bank high slope and compared with the the result generated by Quasi-static method. The comparative analysis reveals that the method not only raises the FEM＇s capability in accurate simulation of complicated geologic structure, but also increases the flexibility and comprehensiveness of limit equilibrium method. This method is reliable and recommended for further application in other real geotechnical engineering.
文摘A supported framework of a gyroscope＇s rotor is designed and the B-Spline wavelet finite element model of nonlinear supported magnetic field is worked out. A new finite element space is studied in which the scaling function of the B-spline wavelet is considered as the shape function of a tetrahedton. The magnetic field is spited by an artificial absorbing body which used the condition of field radiating, so the solution is unique. The resolution is improved via the varying gradient of the B-spline function under the condition of unchanging gridding. So there are some advantages in dealing with the focus flux and a high varying gradient result from a nonlinear magnetic field. The result is more practical. Plots of flux and in the space is studied via simulating the supported system model. The results of the study are useful in the research of the supported magnetic system for the gyroscope rotor.
文摘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.
基金supported by the National Natural Science Foundation of China (50725826)Specific Research on Cable-reinforced Membranes with Super Span and Complex Single-shell Structures of Expo Axis (08dz0580303)Shanghai Postdoctoral Fund (10R21416200)
文摘Based on Timoshenko＇s beam theory and Vlasov＇s thin-walled member theory, a new model of spatial thin-walled beam element is developed for analyzing geometrical and physical nonlinearity, which incorporates an interior node and independent interpolations of bending angles and warp and takes diversified factors into consideration, such as traverse shear deformation, torsional shear deformation and their coupling, coupling of flexure and torsion, and the second shear stress. The geometrical nonlinear strain is formulated in updated Lagarange （UL） and the corresponding stiffness matrix is derived. The perfectly plastic model is used to account for physical nonlinearity, and the yield rule of von Mises and incremental relationship of Prandtle-Reuss are adopted. Elastoplastic stiffness matrix is obtained by numerical integration based on the finite segment method, and a finite element program is compiled. Numerical examples manifest that the proposed model is accurate and feasible in the analysis of thin-walled structures.
基金supported by the National Natural Science Foundation of China(Nos.51378293,51078199,50678093,and 50278046)the Program for Changjiang Scholars and the Innovative Research Team in University of China(No.IRT00736)
文摘The element energy projection （EEP） method for computation of super- convergent resulting in a one-dimensional finite element method （FEM） is successfully used to self-adaptive FEM analysis of various linear problems, based on which this paper presents a substantial extension of the whole set of technology to nonlinear problems. The main idea behind the technology transfer from linear analysis to nonlinear analysis is to use Newton＇s method to linearize nonlinear problems into a series of linear problems so that the EEP formulation and the corresponding adaptive strategy can be directly used without the need for specific super-convergence formulation for nonlinear FEM. As a re- sult, a unified and general self-adaptive algorithm for nonlinear FEM analysis is formed. The proposed algorithm is found to be able to produce satisfactory finite element results with accuracy satisfying the user-preset error tolerances by maximum norm anywhere on the mesh. Taking the nonlinear ordinary differential equation （ODE） of second-order as the model problem, this paper describes the related fundamental idea, the imple- mentation strategy, and the computational algorithm. Representative numerical exam- ples are given to show the efficiency, stability, versatility, and reliability of the proposed approach.
文摘The numerical mode of nonlinear wave transformation based on both the Laplace equation for water field and the Bernoulli equation for water surface is a kind of time-domain boundary problem with initial conditions. And the basis for establishing the numerical mode of nonlinear wave in time domain is to trace the position of wave free surface and to calculate the instantaneous surface height and surface potential function. This paper firstly utilizes the ‘0-1' combined BEM to separate the boundary by means of discretization of Green's integral equation based on the Laplace equation, then separates the free surface of wave with FEM and derives the FEM equation of wave surface that satisfies the nonlinear boundary conditions. By jointly solving the above BEM and FEM equations, the wave potential and surface height could be obtained with iteration in time domain. Thus a new kind of nonlinear numerical mode is established for calculating wave transformation. The wave test in the numerical wave tank shows that the numerical simulation with this mode is of high accuracy.
基金The study was financially supported by the National Natural Science Foundation of China(Grant Nos.10202003 and 50479015)Program for New Century Excellent Talents in University(NCET-05-0710)
文摘In the paper, a weak coupling numerical model is developed for the study of the nonlinear dynamic interaction between water waves and permeable sandy seabed. The wave field solveris based on the VOF （Volume of Fluid） method for continuity equation and the two-dimensional Reynolds Averaged Navier Stokes （RANS） equations with a k-ε closure. The free surface of cnoidal wave is traced through the PLIC-VOF （P/ecewise Linear/nterface Construction）. Blot＇s equations have been applied to solve the sandy seabed, and the u-p fmite dement formulations are derived by the application of the Galerkin weighted-residual procedure. The continuity of the pressure on the interface between fluid and porous medium domains is considered. Laboratory tests were performed to verify the proposed numerical model, and it is shown that the pore-water pressures and the wave heights computed by the VOF-FEM models are in good agreement with the experimental results. It is found that the proposed model is effective in predicting the seabed-nonlinear wave interaction and is able to handle the wave-breakwater-seabed interaction problem.
基金supported by the National Natural Science Foundation of China (Grant No. 50921001)
文摘SUN Da-peng BAO Wei-bin, WU Hao and LI Yu-cheng （ In this paper the 0-1 combined BEM is adopted to subdivide the computational domain boundary, and to discretize the Green＇s integral expression based on Laplace equation. The FEM is used to subdivide the wave surface and deduce the surface equation which satisfies the nonlinear boundary conditions on the surface. The equations with potential function and wave surface height as an unknown quantity by application of Taylor expansion approach can be solved by iteration within the time step. In m-time iteration within the computational process of time step （n-1）At to nat, the results of the previous iteration are taken as the initial value of the two-order unknown terms in the present iteration. Thus, an improved tracking mode of nonlinear wave surface is estabIished, and numerical results of wave tank test indicate that this mode is improved obviously and is more precise than the previous numerical model which ignored the two-order unknown terms of wave surface location and velocity potential function in comparison with the theoretical values.
