Based on local algorithms,some parallel finite element(FE)iterative methods for stationary incompressible magnetohydrodynamics(MHD)are presented.These approaches are on account of two-grid skill include two major phas...Based on local algorithms,some parallel finite element(FE)iterative methods for stationary incompressible magnetohydrodynamics(MHD)are presented.These approaches are on account of two-grid skill include two major phases:find the FE solution by solving the nonlinear system on a globally coarse mesh to seize the low frequency component of the solution,and then locally solve linearized residual subproblems by one of three iterations(Stokes-type,Newton,and Oseen-type)on subdomains with fine grid in parallel to approximate the high frequency component.Optimal error estimates with regard to two mesh sizes and iterative steps of the proposed algorithms are given.Some numerical examples are implemented to verify the algorithm.展开更多
In this paper, we study the mixed element method for Sobolev equations. A time-discretization procedure is presented and analysed and the optimal order error estimates are derived.For convenience in practical computat...In this paper, we study the mixed element method for Sobolev equations. A time-discretization procedure is presented and analysed and the optimal order error estimates are derived.For convenience in practical computation, an alternating-direction iterative scheme of the mixed fi-nite element method is formulated and its stability and converbence are proved for the linear prob-lem. A numerical example is provided at the end of this paper.展开更多
This study develops an optimized finite difference iterative （OFDI） scheme for the two-dimensional （2D） viscoelastic wave equation. The OFDI scheme is obtained using a proper orthogonal decomposition （POD） metho...This study develops an optimized finite difference iterative （OFDI） scheme for the two-dimensional （2D） viscoelastic wave equation. The OFDI scheme is obtained using a proper orthogonal decomposition （POD） method. It has sufficiently high accuracy with very few unknowns for the 2D viscoelastic wave equation. Existence, stability, and convergence of the OFDI solutions are analyzed. Numerical simulations verify efficiency and feasibility of the proposed scheme.展开更多
In this paper, a parallel algorithm with iterative form for solving finite element equation is presented. Based on the iterative solution of linear algebra equations, the parallel computational steps are introduced in...In this paper, a parallel algorithm with iterative form for solving finite element equation is presented. Based on the iterative solution of linear algebra equations, the parallel computational steps are introduced in this method. Also by using the weighted residual method and choosing the appropriate weighting functions, the finite element basic form of parallel algorithm is deduced. The program of this algorithm has been realized on the ELXSI-6400 parallel computer of Xi'an Jiaotong University. The computational results show the operational speed will be raised and the CPU time will be cut down effectively. So this method is one kind of effective parallel algorithm for solving the finite element equations of large-scale structures.展开更多
The purpose of this article is to study the weak and strong convergence of implicit iteration process with errors to a common fixed point for a finite family of asymptotically nonexpansive mappings and nonexpansive ma...The purpose of this article is to study the weak and strong convergence of implicit iteration process with errors to a common fixed point for a finite family of asymptotically nonexpansive mappings and nonexpansive mappings in Banach spaces. The results presented in this article extend and improve the corresponding results of [1, 2, 4-9, 11-15].展开更多
There is a common difficulty in elastic-plastic impact codes such as EPIC[2,3] NONSAP[4], etc.. Most of these codes use the simple linear functions usually taken from static problem to represent the displacement compo...There is a common difficulty in elastic-plastic impact codes such as EPIC[2,3] NONSAP[4], etc.. Most of these codes use the simple linear functions usually taken from static problem to represent the displacement components. In such finite element formulation, the stress components are constant in each element and they are discontinuous in any two neighboring elements. Therefore, the bases of using the virtual work principle in such elements are unreliable. In this paper, we introduce a new method, namely, the compatible stress iterative method, to eliminate the above-said difficulty. The calculated examples show that the calculation using the new method in dynamic finite element analysis of high velocity impact is valid and stable, and the element stiffness can be somewhat reduced.展开更多
A combined method for the fast 3-D finite element modeling of defect responses in nondestructive testing of electromagnetics is presented. The method consists of three numerical techniques: zoom-in technique, differen...A combined method for the fast 3-D finite element modeling of defect responses in nondestructive testing of electromagnetics is presented. The method consists of three numerical techniques: zoom-in technique, difference field technique and iterative solution technique. Utilizing the zoom-in technique, the computational zone focuses on a relatively small domain around the defect. Employing the difference field technique, the axisymmetrical field solution corresponding to the case with no defect can be used to simplify the mesh generation and obtain the modeling results quickly. Using the iterative solution technique, the matrix equation system in the 3-D finite element modeling of nondestructive probe signals can easily be solved. The sample calculation shows that the presented method is highly effective and can consequently save significant computer resources.展开更多
This study has focused on developing numerical procedures for the static and dynamic nonlinear analysis of mooring lines. A geometrically nonlinear finite element method using isoparametric cable element with two node...This study has focused on developing numerical procedures for the static and dynamic nonlinear analysis of mooring lines. A geometrically nonlinear finite element method using isoparametric cable element with two nodes is briefly presented on the basis of the total Lagrangian formulation. The static and dynamic equilibrium equations of mooring lines are established. An incremental-iterative method is used to determine the initial static equilibrium state of cable systems under the action of self weights, buoyancy and current. Also the Newmark method is used for dynamic nonlinear analysis of ocean cables. Numerical examples are presented to validate the present numerical method, and examine the effect of various parameters.展开更多
In this work, the two-dimensional convective Brinkman-Forchheimer equa- tions are considered. The well-posedness for the variational problem and its mixed finite element approximation is established, and the error est...In this work, the two-dimensional convective Brinkman-Forchheimer equa- tions are considered. The well-posedness for the variational problem and its mixed finite element approximation is established, and the error estimates based on the conforming approximation are obtained. For the computation, a one-step Newton (or semi-Newton) iteration algorithm initialized using a fixed-point iteration is proposed. Finally, numerical experiments using a Taylor-Hood mixed element built on a structured or unstructured triangular mesh are implemented. The numerical results obtained using the algorithm are compared with the analytic data, and are shown to be in very good agreement. Moreover, the lid-driven problem at Reynolds numbers of 100 and 400 is considered and analyzed.展开更多
The dynamic calculations of slender marine risers, such as Finite Element Method （FEM） or Modal Expansion Solution Method （MESM）, are mainly for the slender structures with their both ends hinged to the surface an...The dynamic calculations of slender marine risers, such as Finite Element Method （FEM） or Modal Expansion Solution Method （MESM）, are mainly for the slender structures with their both ends hinged to the surface and bottom. However, for the re-entry operation, risers held by vessels are in vertical free hanging state, so the displacement and velocity of lower joint would not be zero. For the model of free hanging flexible marine risers, the paper proposed a Finite Difference Approximation （FDA） method for its dynamic calculation. The riser is divided into a reasonable number of rigid discrete segments. And the dynamic model is established based on simple Euler-Bemoulli Beam Theory concerning tension, shear forces and bending moments at each node along the cylindrical structures, which is extendible for different boundary conditions. The governing equations with specific boundary conditions for riser＇s free hanging state are simplified by Keller-box method and solved with Newton iteration algorithm for a stable dynamic solution. The calculation starts when the riser is vertical and still in calm water, and its behavior is obtained along time responding to the lateral forward motion at the top. The dynamic behavior in response to the lateral parametric excitation at the top is also proposed and discussed in this paper.展开更多
Two basic hypothesises of Taylor-Galerkin Finite Element Method are studied in this paper. One of them which is unreasonable is redefined. The only hypothesis becomes the standpoint of Generalized Finite Element. We u...Two basic hypothesises of Taylor-Galerkin Finite Element Method are studied in this paper. One of them which is unreasonable is redefined. The only hypothesis becomes the standpoint of Generalized Finite Element. We use this idea to analysis stream function-vorticity equations with Modified Taylor-Galerkin Finite Element Method, and give the two-step solving method, which makes the solving process more reasonable than ever before. Several computational examples reveal that the results of this new method are satisfied.展开更多
A three-dimensional （3D） predictor-corrector finite difference method for standing wave is developed. It is applied to solve the 3D nonlinear potential flow equa- tions with a free surface. The 3D irregular tank is ...A three-dimensional （3D） predictor-corrector finite difference method for standing wave is developed. It is applied to solve the 3D nonlinear potential flow equa- tions with a free surface. The 3D irregular tank is mapped onto a fixed cubic tank through the proper coordinate transform schemes. The cubic tank is distributed by the staggered meshgrid, and the staggered meshgrid is used to denote the variables of the flow field. The predictor-corrector finite difference method is given to develop the difference equa- tions of the dynamic boundary equation and kinematic boundary equation. Experimental results show that, using the finite difference method of the predictor-corrector scheme, the numerical solutions agree well with the published results. The wave profiles of the standing wave with different amplitudes and wave lengths are studied. The numerical solutions are also analyzed and presented graphically.展开更多
The method of nonlinear finite element reliability analysis (FERA) of slope stability using the technique of slip surface stress analysis (SSA) is studied. The limit state function that can consider the direction of s...The method of nonlinear finite element reliability analysis (FERA) of slope stability using the technique of slip surface stress analysis (SSA) is studied. The limit state function that can consider the direction of slip surface is given, and the formula-tions of FERA based on incremental tangent stiffness method and modified Aitken accelerating algorithm are developed. The limited step length iteration method (LSLIM) is adopted to calculate the reliability index. The nonlinear FERA code using the SSA technique is developed and the main flow chart is illustrated. Numerical examples are used to demonstrate the efficiency and robustness of this method. It is found that the accelerating convergence algorithm proposed in this study proves to be very efficient for it can reduce the iteration number greatly, and LSLIM is also efficient for it can assure the convergence of the iteration of the reliability index.展开更多
A mixed finite element formulation for viscoelastic flows is derived in this paper, in which the FIC （finite incremental calculus） pressure stabilization process and the DEVSS （discrete elastic viscous stress split...A mixed finite element formulation for viscoelastic flows is derived in this paper, in which the FIC （finite incremental calculus） pressure stabilization process and the DEVSS （discrete elastic viscous stress splitting） method using the Crank-Nicolson-based split are introduced within a general framework of the iterative version of the fractional step algorithm. The SU （streamline-upwind） method is particularly chosen to tackle the convective terms in constitutive equations of viscoelastic flows. Thanks to the proposed scheme the finite elements with equal low-order interpolation approximations for stress-velocity-pressure variables can be successfully used even for viscoelastic flows with high Weissenberg numbers. The XPP （extended Pom-Pom） constitutive model for describing viscoelastic behaviors is particularly integrated into the proposed scheme. The numerical results for the 4：1 sudden contraction flow problem demonstrate prominent stability, accuracy and convergence rate of the proposed scheme in both pressure and stress distributions over the flow domain within a wide range of the Weissenberg number, particularly the capability in reproducing the results, which can be used to explain the ＂die swell＂ phenomenon observed in the polymer injection molding process.展开更多
A finite difference method is presented to simulate transverse vibrations of an axially moving string. By discretizing the governing equation and the equation of stress strain relation at different frictional knots, t...A finite difference method is presented to simulate transverse vibrations of an axially moving string. By discretizing the governing equation and the equation of stress strain relation at different frictional knots, two linear sparse finite difference equation systems are obtained. The two explicit difference schemes can be calculated alternatively, which make the computation much more efficient. The numerical method makes the nonlinear model easier to deal with and of truncation errors, O（△t^2 ＋ △x^2）. It also shows quite good stability for small initial values. Numerical examples are presented to demonstrate the efficiency and the stability of the algorithm, and dynamic analysis of a viscoelastic string is given by using the numerical results.展开更多
In this study, for the purpose of improving the efficiency and accuracy of numerical simulation of massive concrete, the symmetric successive over relaxation-preconditioned conjugate gradient method （SSOR-PCGM） with...In this study, for the purpose of improving the efficiency and accuracy of numerical simulation of massive concrete, the symmetric successive over relaxation-preconditioned conjugate gradient method （SSOR-PCGM） with an improved iteration format was derived and applied to solution of large sparse symmetric positive definite linear equations in the computational process of the finite element analysis. A three-dimensional simulation program for massive concrete was developed based on SSOR-PCGM with an improved iteration format. Then, the programs based on the direct method and SSOR-PCGM with an improved iteration format were used for computation of the Guandi roller compacted concrete （RCC） gravity dam and an elastic cube under free expansion. The comparison and analysis of the computational results show that SSOR-PCGM with the improved iteration format occupies much less physical memory and can solve larger-scale problems with much less computing time and flexible control of accuracy.展开更多
The coupling iteration (CI) of the finite element method(FEM) is used to simulate the magnetic and mechanical characteristics for a GMM actuator. The convergent ability under different prestress and different load typ...The coupling iteration (CI) of the finite element method(FEM) is used to simulate the magnetic and mechanical characteristics for a GMM actuator. The convergent ability under different prestress and different load types is investigated. Then the calculated deformations are compared with the experimental values. The results convince that the CI of FEM is suitable for the simulation of energy coupling and transformation mechanism of the GMM. At last, the output deformation properties are studied under different input currents, showing that there is a good compromise between good linearity and large strain under the prestress 6 MPa.展开更多
In this paper, a new elimination of finite differential equations has been discussed. It applies the numerical direct iteration to obtain the residual equations, in which the number of unknowns has been reduced greatl...In this paper, a new elimination of finite differential equations has been discussed. It applies the numerical direct iteration to obtain the residual equations, in which the number of unknowns has been reduced greatly. The solution process is simple and efficient, and the solution is exact展开更多
In this paper, the modified iteration method is further generalized to the study of axisymmetrical postbuckling of thin circular plates and hereby a new approximate analytic solution of the problem is obtained. Furthe...In this paper, the modified iteration method is further generalized to the study of axisymmetrical postbuckling of thin circular plates and hereby a new approximate analytic solution of the problem is obtained. Further utilizations of this method to postbuckling analyses of plates of more complicated structure are expected.展开更多
The influence of oil film inertia on piston skirt lubrication in a high speedengine is investigated by an iteration method that alternately solves the Navier-Stocks equationsand the Reynolds equa-tion by finite elemen...The influence of oil film inertia on piston skirt lubrication in a high speedengine is investigated by an iteration method that alternately solves the Navier-Stocks equationsand the Reynolds equa-tion by finite element method and difference method. The Reynolds lubricationequation including oil film inertia is developed, in which the inertia coefficient is introduced toinvestigate the effect of oil film inertia. The iteration procedure and finite formulation ofsolving the new Reynolds lubrication equation are given to analyze the effect of oil film on pistonskirt in this kind of engine. The calculation results show that the oil film inertia has someeffects on the friction force, pressure force and load capacity of oil film and its effect isobvious for the last. The Reynolds lubrication equation proposed can be also used to analyze thelubrication performance of the piston skirt in low or medium speed engine and some other lubricationproblems generally excluding oil film inertia with the inertia coefficient being set at zero.展开更多
基金Project supported by the National Natural Science Foundation of China(Nos.11971410 and12071404)the Natural Science Foundation of Hunan Province of China(No.2019JJ40279)+2 种基金the Excellent Youth Program of Scientific Research Project of Hunan Provincial Department of Education(Nos.18B064 and 20B564)the China Postdoctoral Science Foundation(Nos.2018T110073 and 2018M631402)the International Scientific and Technological Innovation Cooperation Base of Hunan Province for Computational Science(No.2018WK4006)。
文摘Based on local algorithms,some parallel finite element(FE)iterative methods for stationary incompressible magnetohydrodynamics(MHD)are presented.These approaches are on account of two-grid skill include two major phases:find the FE solution by solving the nonlinear system on a globally coarse mesh to seize the low frequency component of the solution,and then locally solve linearized residual subproblems by one of three iterations(Stokes-type,Newton,and Oseen-type)on subdomains with fine grid in parallel to approximate the high frequency component.Optimal error estimates with regard to two mesh sizes and iterative steps of the proposed algorithms are given.Some numerical examples are implemented to verify the algorithm.
基金the National Natural Science Foundation of China and China State Key Project for Basic Researches
文摘In this paper, we study the mixed element method for Sobolev equations. A time-discretization procedure is presented and analysed and the optimal order error estimates are derived.For convenience in practical computation, an alternating-direction iterative scheme of the mixed fi-nite element method is formulated and its stability and converbence are proved for the linear prob-lem. A numerical example is provided at the end of this paper.
基金Project supported by the National Natural Science Foundation of China(No.11671106)the Fundamental Research Funds for the Central Universities(No.2016MS33)
文摘This study develops an optimized finite difference iterative （OFDI） scheme for the two-dimensional （2D） viscoelastic wave equation. The OFDI scheme is obtained using a proper orthogonal decomposition （POD） method. It has sufficiently high accuracy with very few unknowns for the 2D viscoelastic wave equation. Existence, stability, and convergence of the OFDI solutions are analyzed. Numerical simulations verify efficiency and feasibility of the proposed scheme.
基金This work has been carried out as of a research project which has been supported by the National Structural Strength & Vibration Laboratory of Xi'an Jiaotong University with National Fund
文摘In this paper, a parallel algorithm with iterative form for solving finite element equation is presented. Based on the iterative solution of linear algebra equations, the parallel computational steps are introduced in this method. Also by using the weighted residual method and choosing the appropriate weighting functions, the finite element basic form of parallel algorithm is deduced. The program of this algorithm has been realized on the ELXSI-6400 parallel computer of Xi'an Jiaotong University. The computational results show the operational speed will be raised and the CPU time will be cut down effectively. So this method is one kind of effective parallel algorithm for solving the finite element equations of large-scale structures.
基金The present studies were supported by the Natural Science Foundation of Zhe-jiang Province (Y605191)the Natural Science Foundation of Heilongjiang Province (A0211)the Key Teacher Creating Capacity Fund of Heilongjiang General College (1053G015)the Scientific Research Foundation from Zhejiang Province Education Committee (20051897)the Starting Foundation of Scientific Research from Hangzhou Teacher's College.
文摘The purpose of this article is to study the weak and strong convergence of implicit iteration process with errors to a common fixed point for a finite family of asymptotically nonexpansive mappings and nonexpansive mappings in Banach spaces. The results presented in this article extend and improve the corresponding results of [1, 2, 4-9, 11-15].
文摘There is a common difficulty in elastic-plastic impact codes such as EPIC[2,3] NONSAP[4], etc.. Most of these codes use the simple linear functions usually taken from static problem to represent the displacement components. In such finite element formulation, the stress components are constant in each element and they are discontinuous in any two neighboring elements. Therefore, the bases of using the virtual work principle in such elements are unreliable. In this paper, we introduce a new method, namely, the compatible stress iterative method, to eliminate the above-said difficulty. The calculated examples show that the calculation using the new method in dynamic finite element analysis of high velocity impact is valid and stable, and the element stiffness can be somewhat reduced.
文摘A combined method for the fast 3-D finite element modeling of defect responses in nondestructive testing of electromagnetics is presented. The method consists of three numerical techniques: zoom-in technique, difference field technique and iterative solution technique. Utilizing the zoom-in technique, the computational zone focuses on a relatively small domain around the defect. Employing the difference field technique, the axisymmetrical field solution corresponding to the case with no defect can be used to simplify the mesh generation and obtain the modeling results quickly. Using the iterative solution technique, the matrix equation system in the 3-D finite element modeling of nondestructive probe signals can easily be solved. The sample calculation shows that the presented method is highly effective and can consequently save significant computer resources.
基金supported by the National Natural Science Foundation of China (Grant No.11072052)the National High Technology Research and Development Program of China (863 Program,Grant No.2006AA09A109-3)
文摘This study has focused on developing numerical procedures for the static and dynamic nonlinear analysis of mooring lines. A geometrically nonlinear finite element method using isoparametric cable element with two nodes is briefly presented on the basis of the total Lagrangian formulation. The static and dynamic equilibrium equations of mooring lines are established. An incremental-iterative method is used to determine the initial static equilibrium state of cable systems under the action of self weights, buoyancy and current. Also the Newmark method is used for dynamic nonlinear analysis of ocean cables. Numerical examples are presented to validate the present numerical method, and examine the effect of various parameters.
基金supported by the National Natural Science Foundation of China(Nos.11461068,11362021,and 11401511)the Doctoral Foundation of Xinjiang Uygur Autonomous Region of China(No.BS110101)
文摘In this work, the two-dimensional convective Brinkman-Forchheimer equa- tions are considered. The well-posedness for the variational problem and its mixed finite element approximation is established, and the error estimates based on the conforming approximation are obtained. For the computation, a one-step Newton (or semi-Newton) iteration algorithm initialized using a fixed-point iteration is proposed. Finally, numerical experiments using a Taylor-Hood mixed element built on a structured or unstructured triangular mesh are implemented. The numerical results obtained using the algorithm are compared with the analytic data, and are shown to be in very good agreement. Moreover, the lid-driven problem at Reynolds numbers of 100 and 400 is considered and analyzed.
基金supported and sponsored jointly by the National Natural Science Foundation of China(Grand Nos.51009092 and 50909061)Doctoral Foundation of the Ministry of Education of China(Grand No.20090073120013)the National High Technology Research and Development Program of China(863Program,Grand No.2008AA092301-1)
文摘The dynamic calculations of slender marine risers, such as Finite Element Method （FEM） or Modal Expansion Solution Method （MESM）, are mainly for the slender structures with their both ends hinged to the surface and bottom. However, for the re-entry operation, risers held by vessels are in vertical free hanging state, so the displacement and velocity of lower joint would not be zero. For the model of free hanging flexible marine risers, the paper proposed a Finite Difference Approximation （FDA） method for its dynamic calculation. The riser is divided into a reasonable number of rigid discrete segments. And the dynamic model is established based on simple Euler-Bemoulli Beam Theory concerning tension, shear forces and bending moments at each node along the cylindrical structures, which is extendible for different boundary conditions. The governing equations with specific boundary conditions for riser＇s free hanging state are simplified by Keller-box method and solved with Newton iteration algorithm for a stable dynamic solution. The calculation starts when the riser is vertical and still in calm water, and its behavior is obtained along time responding to the lateral forward motion at the top. The dynamic behavior in response to the lateral parametric excitation at the top is also proposed and discussed in this paper.
文摘Two basic hypothesises of Taylor-Galerkin Finite Element Method are studied in this paper. One of them which is unreasonable is redefined. The only hypothesis becomes the standpoint of Generalized Finite Element. We use this idea to analysis stream function-vorticity equations with Modified Taylor-Galerkin Finite Element Method, and give the two-step solving method, which makes the solving process more reasonable than ever before. Several computational examples reveal that the results of this new method are satisfied.
基金supported by the Yunnan Provincial Applied Basic Research Program of China(No. KKSY201207019)
文摘A three-dimensional （3D） predictor-corrector finite difference method for standing wave is developed. It is applied to solve the 3D nonlinear potential flow equa- tions with a free surface. The 3D irregular tank is mapped onto a fixed cubic tank through the proper coordinate transform schemes. The cubic tank is distributed by the staggered meshgrid, and the staggered meshgrid is used to denote the variables of the flow field. The predictor-corrector finite difference method is given to develop the difference equa- tions of the dynamic boundary equation and kinematic boundary equation. Experimental results show that, using the finite difference method of the predictor-corrector scheme, the numerical solutions agree well with the published results. The wave profiles of the standing wave with different amplitudes and wave lengths are studied. The numerical solutions are also analyzed and presented graphically.
基金supported by the National Natural Science Foundation of China (No. 50748033)the Specific Foundation for PhD of Hefei University of Technology (No. 2007GDBJ044), China
文摘The method of nonlinear finite element reliability analysis (FERA) of slope stability using the technique of slip surface stress analysis (SSA) is studied. The limit state function that can consider the direction of slip surface is given, and the formula-tions of FERA based on incremental tangent stiffness method and modified Aitken accelerating algorithm are developed. The limited step length iteration method (LSLIM) is adopted to calculate the reliability index. The nonlinear FERA code using the SSA technique is developed and the main flow chart is illustrated. Numerical examples are used to demonstrate the efficiency and robustness of this method. It is found that the accelerating convergence algorithm proposed in this study proves to be very efficient for it can reduce the iteration number greatly, and LSLIM is also efficient for it can assure the convergence of the iteration of the reliability index.
基金the National Natural Science Foundation of China (10672033,10590354,90715011 and 10272027)the National Key Basic Research and Development Program (2002CB412709)
文摘A mixed finite element formulation for viscoelastic flows is derived in this paper, in which the FIC （finite incremental calculus） pressure stabilization process and the DEVSS （discrete elastic viscous stress splitting） method using the Crank-Nicolson-based split are introduced within a general framework of the iterative version of the fractional step algorithm. The SU （streamline-upwind） method is particularly chosen to tackle the convective terms in constitutive equations of viscoelastic flows. Thanks to the proposed scheme the finite elements with equal low-order interpolation approximations for stress-velocity-pressure variables can be successfully used even for viscoelastic flows with high Weissenberg numbers. The XPP （extended Pom-Pom） constitutive model for describing viscoelastic behaviors is particularly integrated into the proposed scheme. The numerical results for the 4：1 sudden contraction flow problem demonstrate prominent stability, accuracy and convergence rate of the proposed scheme in both pressure and stress distributions over the flow domain within a wide range of the Weissenberg number, particularly the capability in reproducing the results, which can be used to explain the ＂die swell＂ phenomenon observed in the polymer injection molding process.
基金Project supported by the National Natural Science Foundation of China(No.10472060)the Natural Science Foundation of Shanghai Municipality(No.04ZR14058)the Shanghai Leading Academic Discipline Project(No.Y0103)
文摘A finite difference method is presented to simulate transverse vibrations of an axially moving string. By discretizing the governing equation and the equation of stress strain relation at different frictional knots, two linear sparse finite difference equation systems are obtained. The two explicit difference schemes can be calculated alternatively, which make the computation much more efficient. The numerical method makes the nonlinear model easier to deal with and of truncation errors, O（△t^2 ＋ △x^2）. It also shows quite good stability for small initial values. Numerical examples are presented to demonstrate the efficiency and the stability of the algorithm, and dynamic analysis of a viscoelastic string is given by using the numerical results.
基金supported by the National Natural Science Foundation of China (Grant No.50808066)
文摘In this study, for the purpose of improving the efficiency and accuracy of numerical simulation of massive concrete, the symmetric successive over relaxation-preconditioned conjugate gradient method （SSOR-PCGM） with an improved iteration format was derived and applied to solution of large sparse symmetric positive definite linear equations in the computational process of the finite element analysis. A three-dimensional simulation program for massive concrete was developed based on SSOR-PCGM with an improved iteration format. Then, the programs based on the direct method and SSOR-PCGM with an improved iteration format were used for computation of the Guandi roller compacted concrete （RCC） gravity dam and an elastic cube under free expansion. The comparison and analysis of the computational results show that SSOR-PCGM with the improved iteration format occupies much less physical memory and can solve larger-scale problems with much less computing time and flexible control of accuracy.
基金This project is supported by National Natural Science Foundation of China (No.50077019).
文摘The coupling iteration (CI) of the finite element method(FEM) is used to simulate the magnetic and mechanical characteristics for a GMM actuator. The convergent ability under different prestress and different load types is investigated. Then the calculated deformations are compared with the experimental values. The results convince that the CI of FEM is suitable for the simulation of energy coupling and transformation mechanism of the GMM. At last, the output deformation properties are studied under different input currents, showing that there is a good compromise between good linearity and large strain under the prestress 6 MPa.
文摘In this paper, a new elimination of finite differential equations has been discussed. It applies the numerical direct iteration to obtain the residual equations, in which the number of unknowns has been reduced greatly. The solution process is simple and efficient, and the solution is exact
文摘In this paper, the modified iteration method is further generalized to the study of axisymmetrical postbuckling of thin circular plates and hereby a new approximate analytic solution of the problem is obtained. Further utilizations of this method to postbuckling analyses of plates of more complicated structure are expected.
文摘The influence of oil film inertia on piston skirt lubrication in a high speedengine is investigated by an iteration method that alternately solves the Navier-Stocks equationsand the Reynolds equa-tion by finite element method and difference method. The Reynolds lubricationequation including oil film inertia is developed, in which the inertia coefficient is introduced toinvestigate the effect of oil film inertia. The iteration procedure and finite formulation ofsolving the new Reynolds lubrication equation are given to analyze the effect of oil film on pistonskirt in this kind of engine. The calculation results show that the oil film inertia has someeffects on the friction force, pressure force and load capacity of oil film and its effect isobvious for the last. The Reynolds lubrication equation proposed can be also used to analyze thelubrication performance of the piston skirt in low or medium speed engine and some other lubricationproblems generally excluding oil film inertia with the inertia coefficient being set at zero.
