In this paper, we represent a new numerical method for solving the nonstationary Stokes equations in an unbounded domain. The technique consists in coupling the boundary integral and finite element methods. The variat...In this paper, we represent a new numerical method for solving the nonstationary Stokes equations in an unbounded domain. The technique consists in coupling the boundary integral and finite element methods. The variational formulation and well posedness of the coupling method are obtained. The convergence and optimal estimates for the approximation solution are provided.展开更多
A novel hybrid-stress finite element method is proposed for constructing simple 4-node quadrilateral plane elements, and the new element is denoted as HH4-3fl here. Firstly, the theoretical basis of the traditional hy...A novel hybrid-stress finite element method is proposed for constructing simple 4-node quadrilateral plane elements, and the new element is denoted as HH4-3fl here. Firstly, the theoretical basis of the traditional hybrid-stress elements, i.e., the Hellinger-Reissner variational principle, is replaced by the Hamilton variational principle, in which the number of the stress variables is reduced from 3 to 2. Secondly, three stress parameters and corresponding trial functions are introduced into the system equations. Thirdly, the displacement fields of the conventional bilinear isoparametric element are employed in the new models. Finally, from the stationary condition, the stress parameters can be expressed in terms of the displacement parameters, and thus the new element stiffness matrices can be obtained. Since the required number of stress variables in the Hamilton variational principle is less than that in the Hellinger-Reissner variational principle, and no additional incompatible displacement modes are considered, the new hybrid-stress element is simpler than the traditional ones. Furthermore, in order to improve the accuracy of the stress solutions, two enhanced post-processing schemes are also proposed for element HH4-3β. Numerical examples show that the proposed model exhibits great improvements in both displacement and stress solutions, implying that the proposed technique is an effective way for developing simple finite element models with high performance.展开更多
A hydromechanical interface element is proposed for the consideration of the hydraulic-mechanical coupling effect along the interface.The fully coupled governing equations and the relevant finite element formulations ...A hydromechanical interface element is proposed for the consideration of the hydraulic-mechanical coupling effect along the interface.The fully coupled governing equations and the relevant finite element formulations are derived in detail for the interface element.All the involved matrices are of the same form as those of a solid element,which makes the incorporation of the model into a finite element program straightforward.Three examples are then numerically simulated using the interface element.Reasonable results confirm the correctness of the proposed model and motivate its application in hydromechanical contact problems in the future.展开更多
The quadrilateral discrete Kirchhoff thin plate bending element DKQ is based on the isoparametric element Q8, however, the accuracy of the isoparametric quadrilateral elements will drop significantly due to mesh disto...The quadrilateral discrete Kirchhoff thin plate bending element DKQ is based on the isoparametric element Q8, however, the accuracy of the isoparametric quadrilateral elements will drop significantly due to mesh distortions. In a previous work, we constructed an 8-node quadrilateral spline element L8 using the triangular area coordinates and the B- net method, which can be insensitive to mesh distortions and possess the second order completeness in the Cartesian co- ordinates. In this paper, a thin plate spline element is devel- oped based on the spline element L8 and the refined tech- nique. Numerical examples show that the present element indeed possesses higher accuracy than the DKQ element for distorted meshes.展开更多
A specific computational program SAFEM was developed based on semi-analytical finite element （FE） method for analysis of asphalt pavement structural responses under static loads. The reliability and efficiency of th...A specific computational program SAFEM was developed based on semi-analytical finite element （FE） method for analysis of asphalt pavement structural responses under static loads. The reliability and efficiency of this FE program was proved by comparison with the general commercial FE software ABAQUS. In order to further reduce the computational time without decrease of the accuracy, the infinite element was added to this program. The results of the finite-infinite element coupling analysis were compared with those of finite element analysis derived from the verified FE program, The study shows that finite-infinite element coupling analysis has higher reliability and efficiency.展开更多
This paper deals with Raviart-Thomas element (Q(2,1) x Q(1,2) - Q(1) element). Apart from its global superconvergence property of fourth order, we prove that a postprocessed extrapolation can globally increased the ac...This paper deals with Raviart-Thomas element (Q(2,1) x Q(1,2) - Q(1) element). Apart from its global superconvergence property of fourth order, we prove that a postprocessed extrapolation can globally increased the accuracy by fifth order.展开更多
The paper presents the variational formulation and well posedness of the coupling method offinite elements and boundary elements for radiation problem. The convergence and optimal errorestimate for the approximate sol...The paper presents the variational formulation and well posedness of the coupling method offinite elements and boundary elements for radiation problem. The convergence and optimal errorestimate for the approximate solution and numerical experiment are provided.展开更多
A finite element / boundary element-modified modal decomposition method (FBMMD) is presented for predicting the vibration and sound radiation from submerged shell of revolution. Improvement has been made to accelerate...A finite element / boundary element-modified modal decomposition method (FBMMD) is presented for predicting the vibration and sound radiation from submerged shell of revolution. Improvement has been made to accelerate the convergence to FBMD method by means of introducing the residual modes which take into accaunt the quasi -state contributiort of all neglected modes. As an example, the vibration and sound radiation of a submerged spherical shell excited by axisymmetric force are studied in cases of ka=l,2,3 and 4. From the calculated results we see that the FBMMD method shows a significant improvement to the accuracy of surface sound pressure, normal displacement and directivity patterns of radiating sound, especially to the directivity patterns.展开更多
Natural geological structures in rock(e.g.,joints,weakness planes,defects)play a vital role in the stability of tunnels and underground operations during construction.We investigated the failure characteristics of a d...Natural geological structures in rock(e.g.,joints,weakness planes,defects)play a vital role in the stability of tunnels and underground operations during construction.We investigated the failure characteristics of a deep circular tunnel in a rock mass with multiple weakness planes using a 2D combined finite element method/discrete element method(FEM/DEM).Conventional triaxial compression tests were performed on typical hard rock(marble)specimens under a range of confinement stress conditions to validate the rationale and accuracy of the proposed numerical approach.Parametric analysis was subsequently conducted to investigate the influence of inclination angle,and length on the crack propagation behavior,failure mode,energy evolution,and displacement distribution of the surrounding rock.The results show that the inclination angle strongly affects tunnel stability,and the failure intensity and damage range increase with increasing inclination angle and then decrease.The dynamic disasters are more likely with increasing weak plane length.Shearing and sliding along multiple weak planes are also consistently accompanied by kinetic energy fluctuations and surges after unloading,which implies a potentially violent dynamic response around a deeply-buried tunnel.Interactions between slabbing and shearing near the excavation boundaries are also discussed.The results presented here provide important insight into deep tunnel failure in hard rock influenced by both unloading disturbance and tectonic activation.展开更多
A simulative analysis coupled with experiment on behaviors of a soil bed cut by a model bulldozer blade is carried out using the finite element/distinct element method（FE/DEM） facility built in the ELFEN package. Be...A simulative analysis coupled with experiment on behaviors of a soil bed cut by a model bulldozer blade is carried out using the finite element/distinct element method（FE/DEM） facility built in the ELFEN package. Before simulation, tensile/compression, triaxial compression and the soil specimens are examined through uniaxial direct shear tests to obtain model characteristics and relevant parameters, then soil cutting experiments are carried out via a mini-soil bin system with a soil bed of 60/120 mm in width and 10 mm in depth cut by a 1/9 scale model bulldozer blade moving with the velocity of 10 mm/s. The soil constitutive model includes the tensile elastic model for tensile breakage and the compressive elastoplastic relationship with Mohr-Coulomb criterion. The cutting length in simulation is set as 1/4 of that in the experiment divided into 1 869 triangular elements. The comparison between the simulated results and experimental ones shows that the used model is capable of analyzing soil dynamic behaviors qualitatively, and the predicted fracturing profiles in general conform to the experiment. Hence the feasibility for analyzing soil fracturing behaviors in tillage or other similar processes is validated.展开更多
In this paper, the coupled thermo-mechanical （TM） processes in the AEspoe Pillar Stability Experiment （APSE） carried out by the Swedish Nuclear Fuel and Waste Management Company （SKB） were simulated using both c...In this paper, the coupled thermo-mechanical （TM） processes in the AEspoe Pillar Stability Experiment （APSE） carried out by the Swedish Nuclear Fuel and Waste Management Company （SKB） were simulated using both continuum and discontinuum based numerical methods. Two-dimensional （2D） and three- dimensional （3D） finite element method （FEM） and 2D distinct element method （DEM） with particles were used. The main objective for the large scale in situ experiment is to investigate the yielding strength of crystalline rock and the formation of the excavation disturbed/damaged zone （EDZ） during excavation of two boreholes, pressurizing of one of the boreholes and heating. For the DEM simulations, the heat flow algorithm was newly introduced into the original code. The calculated stress, displacement and temperature distributions were compared with the ones obtained from in situ measurements and FEM simulations. A parametric study for initial microcracks was also performed to reproduce the spalling phenomena observed in the APSE.展开更多
Top-down crack in asphalt pavements has been reported as a widespread mode of failure.A solid understanding of the mechanisms of crack growth is essential to predict pavement performance in the context of thickness de...Top-down crack in asphalt pavements has been reported as a widespread mode of failure.A solid understanding of the mechanisms of crack growth is essential to predict pavement performance in the context of thickness design,as well as in the design and optimization of mixtures.Using the coupled element free Galerkin (EFG) and finite element (FE) method,top-down crack propagation in asphalt pavements is numerically simulated on the basis of fracture mechanics.A parametric study is conducted to isolate the effects of overlay thickness and stiffness,base thickness and stiffness on top-down crack propagation in asphalt pavements.The results show that longitudinal wheel loads are disadvantageous to top-down crack because it increases the compound stress intensity factor (SIF) at the tip of top-down crack and shortens the crack path,and thus the fatigue life descends.The SIF experiences a process "sharply ascending—slowly descending—slowly ascending—sharply ascending again" with the crack propagating.The thicker the overlay or the base,the lower the SIF; the greater the overlay stiffness,the higher the SIF.The crack path is hardly affected by stiffness of the overlay and base.展开更多
The interaction problem among fractures under the action of compressional stress is studied in this paper by using the finite element method and boundary element method respectively.The mechanical criteria which diffe...The interaction problem among fractures under the action of compressional stress is studied in this paper by using the finite element method and boundary element method respectively.The mechanical criteria which differentiate between the independent fractures and fracture systems and their computation methods are presented in this paper.The proportional conditions between length and spacing of fractures that exist interaction for several kinds of fracture groups of different geometric arrangement are given.The effect of interaction among fractures on the displacement field,stress field and strain energy distribution are computed.The relations between the fracture system of conjugate array and conjugate earthquakes are also discussed in this paper.展开更多
Computational models provide additional tools for studying the brain,however,many techniques are currently disconnected from each other.There is a need for new computational approaches that span the range of physics o...Computational models provide additional tools for studying the brain,however,many techniques are currently disconnected from each other.There is a need for new computational approaches that span the range of physics operating in the brain.In this review paper,we offer some new perspectives on how the embedded element method can fill this gap and has the potential to connect a myriad of modeling genre.The embedded element method is a mesh superposition technique used within finite element analysis.This method allows for the incorporation of axonal fiber tracts to be explicitly represented.Here,we explore the use of the approach beyond its original goal of predicting axonal strain in brain injury.We explore the potential application of the embedded element method in areas of electrophysiology,neurodegeneration,neuropharmacology and mechanobiology.We conclude that this method has the potential to provide us with an integrated computational framework that can assist in developing improved diagnostic tools and regeneration technologies.展开更多
The paper presents two methods for the formulation of free vibration analysis of collecting electrodes of precipitators.The first,called the hybrid finite element method, combines the finit element method used for cal...The paper presents two methods for the formulation of free vibration analysis of collecting electrodes of precipitators.The first,called the hybrid finite element method, combines the finit element method used for calculations of spring deformations with the rigid finite element method used to reflect mass and geometrical features,which is called the hybrid finite element method.As a result,a model with a diagonal mass matrix is obtained.Due to a specific geometry of the electrodes,which are long plates of complicated shapes,the second method proposed is the strip method which is a semi-analytical method.The strip method allows us to formulate the equations of motion with a considerably smaller number of generalized coordinates.Results of numerical calculations obtained by both methods are compared with those obtained using commercial software like ANSYS and ABAQUS.Good compatibility of results is achieved.展开更多
In this paper, we extend the reduced basis methods for parameter dependent problems to the parareal in time algorithm introduced by Lions et al. [12] and solve a nonlinear evolutionary parabolic partial differential e...In this paper, we extend the reduced basis methods for parameter dependent problems to the parareal in time algorithm introduced by Lions et al. [12] and solve a nonlinear evolutionary parabolic partial differential equation. The fine solver is based on the finite element method or spectral element method in space and a semi-implicit Runge-Kutta scheme in time. The coarse solver is based on a semi-implicit scheme in time and the reduced basis approximation in space. Of[line-online procedures are developed, and it is proved that the computational complexity of the on-line stage depends only on the dimension of the reduced basis space （typically small）. Parareal in time algorithms based on a multi-grids finite element method and a multi-degrees finite element method are also presented. Some numerical results are reported.展开更多
文摘In this paper, we represent a new numerical method for solving the nonstationary Stokes equations in an unbounded domain. The technique consists in coupling the boundary integral and finite element methods. The variational formulation and well posedness of the coupling method are obtained. The convergence and optimal estimates for the approximation solution are provided.
基金supported by the National Natural Science Foundation of China (10872108,10876100)the Program for New Century Excellent Talents in University (NCET-07-0477)the National Basic Research Program of China (2010CB832701)
文摘A novel hybrid-stress finite element method is proposed for constructing simple 4-node quadrilateral plane elements, and the new element is denoted as HH4-3fl here. Firstly, the theoretical basis of the traditional hybrid-stress elements, i.e., the Hellinger-Reissner variational principle, is replaced by the Hamilton variational principle, in which the number of the stress variables is reduced from 3 to 2. Secondly, three stress parameters and corresponding trial functions are introduced into the system equations. Thirdly, the displacement fields of the conventional bilinear isoparametric element are employed in the new models. Finally, from the stationary condition, the stress parameters can be expressed in terms of the displacement parameters, and thus the new element stiffness matrices can be obtained. Since the required number of stress variables in the Hamilton variational principle is less than that in the Hellinger-Reissner variational principle, and no additional incompatible displacement modes are considered, the new hybrid-stress element is simpler than the traditional ones. Furthermore, in order to improve the accuracy of the stress solutions, two enhanced post-processing schemes are also proposed for element HH4-3β. Numerical examples show that the proposed model exhibits great improvements in both displacement and stress solutions, implying that the proposed technique is an effective way for developing simple finite element models with high performance.
基金supported by the Innovation Plan for Postgraduate Students sponsored by the Education Department of Jiangsu Province,China (CX08B 107Z)
文摘A hydromechanical interface element is proposed for the consideration of the hydraulic-mechanical coupling effect along the interface.The fully coupled governing equations and the relevant finite element formulations are derived in detail for the interface element.All the involved matrices are of the same form as those of a solid element,which makes the incorporation of the model into a finite element program straightforward.Three examples are then numerically simulated using the interface element.Reasonable results confirm the correctness of the proposed model and motivate its application in hydromechanical contact problems in the future.
基金supported by the National Natural Science Foundation of China(11001037,11102037,11290143)the Fundamental Research Funds for the Central Universities(DUT13LK07)
文摘The quadrilateral discrete Kirchhoff thin plate bending element DKQ is based on the isoparametric element Q8, however, the accuracy of the isoparametric quadrilateral elements will drop significantly due to mesh distortions. In a previous work, we constructed an 8-node quadrilateral spline element L8 using the triangular area coordinates and the B- net method, which can be insensitive to mesh distortions and possess the second order completeness in the Cartesian co- ordinates. In this paper, a thin plate spline element is devel- oped based on the spline element L8 and the refined tech- nique. Numerical examples show that the present element indeed possesses higher accuracy than the DKQ element for distorted meshes.
基金represented by German Federal Highway Research Institute (BASt)financed by the Federal Minister of Transport and Digital Infrastructure (BMVI)conducted under FE 04.0259/2012/NGB
文摘A specific computational program SAFEM was developed based on semi-analytical finite element （FE） method for analysis of asphalt pavement structural responses under static loads. The reliability and efficiency of this FE program was proved by comparison with the general commercial FE software ABAQUS. In order to further reduce the computational time without decrease of the accuracy, the infinite element was added to this program. The results of the finite-infinite element coupling analysis were compared with those of finite element analysis derived from the verified FE program, The study shows that finite-infinite element coupling analysis has higher reliability and efficiency.
文摘This paper deals with Raviart-Thomas element (Q(2,1) x Q(1,2) - Q(1) element). Apart from its global superconvergence property of fourth order, we prove that a postprocessed extrapolation can globally increased the accuracy by fifth order.
基金This research was supported in part by the Institute for Mathematics and its applications with funds provided by NSF, USA
文摘The paper presents the variational formulation and well posedness of the coupling method offinite elements and boundary elements for radiation problem. The convergence and optimal errorestimate for the approximate solution and numerical experiment are provided.
文摘A finite element / boundary element-modified modal decomposition method (FBMMD) is presented for predicting the vibration and sound radiation from submerged shell of revolution. Improvement has been made to accelerate the convergence to FBMD method by means of introducing the residual modes which take into accaunt the quasi -state contributiort of all neglected modes. As an example, the vibration and sound radiation of a submerged spherical shell excited by axisymmetric force are studied in cases of ka=l,2,3 and 4. From the calculated results we see that the FBMMD method shows a significant improvement to the accuracy of surface sound pressure, normal displacement and directivity patterns of radiating sound, especially to the directivity patterns.
基金Projects(52004143,51774194)supported by the National Natural Science Foundation of ChinaProject(2020M670781)supported by the China Postdoctoral Science Foundation+2 种基金Project(SKLGDUEK2021)supported by the State Key Laboratory for GeoMechanics and Deep Underground Engineering,ChinaProject(U1806208)supported by the NSFC-Shandong Joint Fund,ChinaProject(2018GSF117023)supported by the Key Research and Development Program of Shandong Province,China。
文摘Natural geological structures in rock(e.g.,joints,weakness planes,defects)play a vital role in the stability of tunnels and underground operations during construction.We investigated the failure characteristics of a deep circular tunnel in a rock mass with multiple weakness planes using a 2D combined finite element method/discrete element method(FEM/DEM).Conventional triaxial compression tests were performed on typical hard rock(marble)specimens under a range of confinement stress conditions to validate the rationale and accuracy of the proposed numerical approach.Parametric analysis was subsequently conducted to investigate the influence of inclination angle,and length on the crack propagation behavior,failure mode,energy evolution,and displacement distribution of the surrounding rock.The results show that the inclination angle strongly affects tunnel stability,and the failure intensity and damage range increase with increasing inclination angle and then decrease.The dynamic disasters are more likely with increasing weak plane length.Shearing and sliding along multiple weak planes are also consistently accompanied by kinetic energy fluctuations and surges after unloading,which implies a potentially violent dynamic response around a deeply-buried tunnel.Interactions between slabbing and shearing near the excavation boundaries are also discussed.The results presented here provide important insight into deep tunnel failure in hard rock influenced by both unloading disturbance and tectonic activation.
基金This project is supported by National Natural Science Foundation of China (No. 10372113)Royal Society-NSFC China-UK Joint Project (No. 16468).
文摘A simulative analysis coupled with experiment on behaviors of a soil bed cut by a model bulldozer blade is carried out using the finite element/distinct element method（FE/DEM） facility built in the ELFEN package. Before simulation, tensile/compression, triaxial compression and the soil specimens are examined through uniaxial direct shear tests to obtain model characteristics and relevant parameters, then soil cutting experiments are carried out via a mini-soil bin system with a soil bed of 60/120 mm in width and 10 mm in depth cut by a 1/9 scale model bulldozer blade moving with the velocity of 10 mm/s. The soil constitutive model includes the tensile elastic model for tensile breakage and the compressive elastoplastic relationship with Mohr-Coulomb criterion. The cutting length in simulation is set as 1/4 of that in the experiment divided into 1 869 triangular elements. The comparison between the simulated results and experimental ones shows that the used model is capable of analyzing soil dynamic behaviors qualitatively, and the predicted fracturing profiles in general conform to the experiment. Hence the feasibility for analyzing soil fracturing behaviors in tillage or other similar processes is validated.
基金conducted within the context of the international DECOVALEX Project (DEvelopment of COupled models and their VALidation against EXperiments)financed by Japan Atomic Energy Agency (JAEA) who was also one of the Funding Organizations of the projectChrister Anders-son from Swedish Nuclear Fuel and Waste Management Co.(SKB),Sweden
文摘In this paper, the coupled thermo-mechanical （TM） processes in the AEspoe Pillar Stability Experiment （APSE） carried out by the Swedish Nuclear Fuel and Waste Management Company （SKB） were simulated using both continuum and discontinuum based numerical methods. Two-dimensional （2D） and three- dimensional （3D） finite element method （FEM） and 2D distinct element method （DEM） with particles were used. The main objective for the large scale in situ experiment is to investigate the yielding strength of crystalline rock and the formation of the excavation disturbed/damaged zone （EDZ） during excavation of two boreholes, pressurizing of one of the boreholes and heating. For the DEM simulations, the heat flow algorithm was newly introduced into the original code. The calculated stress, displacement and temperature distributions were compared with the ones obtained from in situ measurements and FEM simulations. A parametric study for initial microcracks was also performed to reproduce the spalling phenomena observed in the APSE.
基金Project (Nos. 50908093 and 50778077) supported by the National Natural Science Foundation of China
文摘Top-down crack in asphalt pavements has been reported as a widespread mode of failure.A solid understanding of the mechanisms of crack growth is essential to predict pavement performance in the context of thickness design,as well as in the design and optimization of mixtures.Using the coupled element free Galerkin (EFG) and finite element (FE) method,top-down crack propagation in asphalt pavements is numerically simulated on the basis of fracture mechanics.A parametric study is conducted to isolate the effects of overlay thickness and stiffness,base thickness and stiffness on top-down crack propagation in asphalt pavements.The results show that longitudinal wheel loads are disadvantageous to top-down crack because it increases the compound stress intensity factor (SIF) at the tip of top-down crack and shortens the crack path,and thus the fatigue life descends.The SIF experiences a process "sharply ascending—slowly descending—slowly ascending—sharply ascending again" with the crack propagating.The thicker the overlay or the base,the lower the SIF; the greater the overlay stiffness,the higher the SIF.The crack path is hardly affected by stiffness of the overlay and base.
文摘The interaction problem among fractures under the action of compressional stress is studied in this paper by using the finite element method and boundary element method respectively.The mechanical criteria which differentiate between the independent fractures and fracture systems and their computation methods are presented in this paper.The proportional conditions between length and spacing of fractures that exist interaction for several kinds of fracture groups of different geometric arrangement are given.The effect of interaction among fractures on the displacement field,stress field and strain energy distribution are computed.The relations between the fracture system of conjugate array and conjugate earthquakes are also discussed in this paper.
基金support provided by Computational Fluid Dynamics Research Corporation(CFDRC)under a sub-contract funded by the Department of Defense,Department of Health Program through contract W81XWH-14-C-0045
文摘Computational models provide additional tools for studying the brain,however,many techniques are currently disconnected from each other.There is a need for new computational approaches that span the range of physics operating in the brain.In this review paper,we offer some new perspectives on how the embedded element method can fill this gap and has the potential to connect a myriad of modeling genre.The embedded element method is a mesh superposition technique used within finite element analysis.This method allows for the incorporation of axonal fiber tracts to be explicitly represented.Here,we explore the use of the approach beyond its original goal of predicting axonal strain in brain injury.We explore the potential application of the embedded element method in areas of electrophysiology,neurodegeneration,neuropharmacology and mechanobiology.We conclude that this method has the potential to provide us with an integrated computational framework that can assist in developing improved diagnostic tools and regeneration technologies.
基金Research is financed from the project NR03-0036-04/2008
文摘The paper presents two methods for the formulation of free vibration analysis of collecting electrodes of precipitators.The first,called the hybrid finite element method, combines the finit element method used for calculations of spring deformations with the rigid finite element method used to reflect mass and geometrical features,which is called the hybrid finite element method.As a result,a model with a diagonal mass matrix is obtained.Due to a specific geometry of the electrodes,which are long plates of complicated shapes,the second method proposed is the strip method which is a semi-analytical method.The strip method allows us to formulate the equations of motion with a considerably smaller number of generalized coordinates.Results of numerical calculations obtained by both methods are compared with those obtained using commercial software like ANSYS and ABAQUS.Good compatibility of results is achieved.
文摘In this paper, we extend the reduced basis methods for parameter dependent problems to the parareal in time algorithm introduced by Lions et al. [12] and solve a nonlinear evolutionary parabolic partial differential equation. The fine solver is based on the finite element method or spectral element method in space and a semi-implicit Runge-Kutta scheme in time. The coarse solver is based on a semi-implicit scheme in time and the reduced basis approximation in space. Of[line-online procedures are developed, and it is proved that the computational complexity of the on-line stage depends only on the dimension of the reduced basis space （typically small）. Parareal in time algorithms based on a multi-grids finite element method and a multi-degrees finite element method are also presented. Some numerical results are reported.
