Pipe cooling is an effective method of mass concrete temperature control, but its accurate and convenient numerical simulation is still a cumbersome problem. An improved embedded model, considering the water temperatu...Pipe cooling is an effective method of mass concrete temperature control, but its accurate and convenient numerical simulation is still a cumbersome problem. An improved embedded model, considering the water temperature variation along the pipe, was proposed for simulating the temperature field of early-age concrete structures containing cooling pipes. The improved model was verified with an engineering example. Then, the p-version self-adaption algorithm for the improved embedded model was deduced, and the initial values and boundary conditions were examined. Comparison of some numerical samples shows that the proposed model can provide satisfying precision and a higher efficiency. The analysis efficiency can be doubled at the same precision, even for a large-scale element. The p-version algorithm can fit grids of different sizes for the temperature field simulation. The convenience of the proposed algorithm lies in the possibility of locating more pipe segments in one element without the need of so regular a shape as in the explicit model.展开更多
We propose a high-order enriched partition of unity finite element method for linear and nonlinear time-dependent diffusion problems.The solution of this class of problems often exhibits non-smooth features such as st...We propose a high-order enriched partition of unity finite element method for linear and nonlinear time-dependent diffusion problems.The solution of this class of problems often exhibits non-smooth features such as steep gradients and boundary layers which can be very challenging to recover using the conventional low-order finite element methods.A class of steady-state exponential functions has been widely used for enrichment and its performance to numerically solve these challenges has been demonstrated.However,these enrichment functions have been used only in context of the standard h-version refinement or the so-called q-version refinement.In this paper we demonstrate that the p-version refinement can also be a very attractive option in terms of the efficiency and the accuracy in the enriched partition of unity finite element method.First,the transient diffusion problem is integrated in time using a semi-implicit scheme and the semi-discrete problem is then integrated in space using the p-version enriched finite elements.Numerical results are presented for three test examples of timedependent diffusion problems in both homogeneous and heterogeneous media.The computed results show the significant improvement when using the p-version refined enriched approximations in the finite element analysis.In addition,these results support our expectations for a robust and high-order accurate enriched partition of unity finite element method.展开更多
Nanoplates have been widely used as elementary components for ultrasensitive and ultrafine resolution applications in the field of nano-electro-mechanical systems because of their potentially remarkable mechanical pro...Nanoplates have been widely used as elementary components for ultrasensitive and ultrafine resolution applications in the field of nano-electro-mechanical systems because of their potentially remarkable mechanical properties.The accurate analysis of their mechanical behavior is currently of particular interest in the function design and reliability analysis of nano-scaled devices.To examine the size-dependent bending and vibration behavior of nanoplates with curvilinear and irregular shapes,a new p-version curved C^(1)finite element is formulated in the framework of the nonlocal Kirchhoff plate model.This newly developed element not only enables an accurate geometry representation and easy mesh generation of curvilinear domains but also overcomes the difficulty of imposing C^(1)conformity required by the nonlocal Kirchhoff plate model,particularly on the curvilinear inter-element boundaries.Numerical examples show that this element can produce an exponential rate of convergence even when curved elements are used in the domain discretization.Vast numerical results are presented for nanoplates with various geometric shapes,including rectangular,circular,elliptic,annular,and sectorial.The high accuracy of the present element is verified by comparing the obtained results with analytical and numerical results in the literature.Additionally,a comprehensive parametric analysis is conducted to investigate the influences of nonlocal parameters,plate dimensions,and boundary conditions on the nonlocal behavior of nanoplates.The present element can be envisaged to allow large-scale mechanical simulations of nanoplates,with a guarantee of accuracy and efficiency.展开更多
A novel two level spline method is proposed for semi-linear elliptic equations, where the two level iteration is implemented between a pair of hierarchical spline spaces with different orders. The new two level method...A novel two level spline method is proposed for semi-linear elliptic equations, where the two level iteration is implemented between a pair of hierarchical spline spaces with different orders. The new two level method is implementation in a manner of p-adaptivity. A coarse solution is obtained from solving the model problem in the low order spline space, and the solution with higher accuracy are generated subsequently, via one step Newton or monidifed Newton iteration in the high order spline space. We also derive the optimal error estimations for the proposed two level schemes. At last, the illustrated numerical results confirm our error estimations and further research topics are commented.展开更多
The stress intensity factors(SIFs)for two-dimensional cracks are extracted using the p-version finite element method(P-FEM)and the contour integral method.Several numerical experiments,e.g.,crack initiating from the e...The stress intensity factors(SIFs)for two-dimensional cracks are extracted using the p-version finite element method(P-FEM)and the contour integral method.Several numerical experiments,e.g.,crack initiating from the edge of a circular hole under an unidirectional uniform tension and two equal-length,unequal-length hole-edge cracks,respectively,at a rectangular plate,an inclined centered crack under uniaxial tension at a square plate and a pipeline crack model,are used to demonstrate the accuracy and effectiveness of the approaches.SIFs are presented for the effects of various crack lengths and length-width ratio.Numerical results are analyzed and compared with reference solutions and results obtained by the Voronoi cell finite element method,boundary element method,high-order extended finite element method(high-order XFEM)and commercial finite element software ABAQUS in the available literature.Numerical results are in good agreement with the benchmark problems and show faster convergence rate,higher accuracy and better numerical stability.展开更多
An elasto-plastie impact model based on the p-version finite element method is presented for the collision protection of ocean and offshore structures. The impact force and responses of the impactor-absorber-structure...An elasto-plastie impact model based on the p-version finite element method is presented for the collision protection of ocean and offshore structures. The impact force and responses of the impactor-absorber-structure system can be predicted efficiently and automatically. A cost-effective Cellular Reinforced Concrete Absorber （CRCA） is designed to smooth the impact force and absorb the impact energy. Quasi-static tests show that the concrete absorber has an excellent energy absorbing characteristic. The impact experiment of a scaled offshore oil-piping frame with the proposed concrete absorber is carried out. The simulation results of the elasto-plastie model and the p-version finite element method are in good agreement with the experimental ones. Owing to the plastic deformation of the absorber, the impact force during the impact and responses of the structure are considerably reduced. Further, the proposed impact model and the concrete absorber are applied to the design of collision protection of the sheet-pile groin on the Qiantang River used to weaken the famous Qiantang bore.展开更多
基金supported by the National Natural Science Foundation of China(Grant No.51109071)
文摘Pipe cooling is an effective method of mass concrete temperature control, but its accurate and convenient numerical simulation is still a cumbersome problem. An improved embedded model, considering the water temperature variation along the pipe, was proposed for simulating the temperature field of early-age concrete structures containing cooling pipes. The improved model was verified with an engineering example. Then, the p-version self-adaption algorithm for the improved embedded model was deduced, and the initial values and boundary conditions were examined. Comparison of some numerical samples shows that the proposed model can provide satisfying precision and a higher efficiency. The analysis efficiency can be doubled at the same precision, even for a large-scale element. The p-version algorithm can fit grids of different sizes for the temperature field simulation. The convenience of the proposed algorithm lies in the possibility of locating more pipe segments in one element without the need of so regular a shape as in the explicit model.
文摘We propose a high-order enriched partition of unity finite element method for linear and nonlinear time-dependent diffusion problems.The solution of this class of problems often exhibits non-smooth features such as steep gradients and boundary layers which can be very challenging to recover using the conventional low-order finite element methods.A class of steady-state exponential functions has been widely used for enrichment and its performance to numerically solve these challenges has been demonstrated.However,these enrichment functions have been used only in context of the standard h-version refinement or the so-called q-version refinement.In this paper we demonstrate that the p-version refinement can also be a very attractive option in terms of the efficiency and the accuracy in the enriched partition of unity finite element method.First,the transient diffusion problem is integrated in time using a semi-implicit scheme and the semi-discrete problem is then integrated in space using the p-version enriched finite elements.Numerical results are presented for three test examples of timedependent diffusion problems in both homogeneous and heterogeneous media.The computed results show the significant improvement when using the p-version refined enriched approximations in the finite element analysis.In addition,these results support our expectations for a robust and high-order accurate enriched partition of unity finite element method.
基金the National Major Science and Technology Projects of China(Grant No.J2019-VI-0001-0114)the National Natural Science Foundation of China(Grant Nos.11972004,11772031,11402015)。
文摘Nanoplates have been widely used as elementary components for ultrasensitive and ultrafine resolution applications in the field of nano-electro-mechanical systems because of their potentially remarkable mechanical properties.The accurate analysis of their mechanical behavior is currently of particular interest in the function design and reliability analysis of nano-scaled devices.To examine the size-dependent bending and vibration behavior of nanoplates with curvilinear and irregular shapes,a new p-version curved C^(1)finite element is formulated in the framework of the nonlocal Kirchhoff plate model.This newly developed element not only enables an accurate geometry representation and easy mesh generation of curvilinear domains but also overcomes the difficulty of imposing C^(1)conformity required by the nonlocal Kirchhoff plate model,particularly on the curvilinear inter-element boundaries.Numerical examples show that this element can produce an exponential rate of convergence even when curved elements are used in the domain discretization.Vast numerical results are presented for nanoplates with various geometric shapes,including rectangular,circular,elliptic,annular,and sectorial.The high accuracy of the present element is verified by comparing the obtained results with analytical and numerical results in the literature.Additionally,a comprehensive parametric analysis is conducted to investigate the influences of nonlocal parameters,plate dimensions,and boundary conditions on the nonlocal behavior of nanoplates.The present element can be envisaged to allow large-scale mechanical simulations of nanoplates,with a guarantee of accuracy and efficiency.
文摘A novel two level spline method is proposed for semi-linear elliptic equations, where the two level iteration is implemented between a pair of hierarchical spline spaces with different orders. The new two level method is implementation in a manner of p-adaptivity. A coarse solution is obtained from solving the model problem in the low order spline space, and the solution with higher accuracy are generated subsequently, via one step Newton or monidifed Newton iteration in the high order spline space. We also derive the optimal error estimations for the proposed two level schemes. At last, the illustrated numerical results confirm our error estimations and further research topics are commented.
基金the firancinal support of the National Natural Science Foundation of China(Grant No:51769011)for this work,and the authors are also deeply grateful to the editors and revewerse for tbeir rigorous work and valuable comments.
文摘The stress intensity factors(SIFs)for two-dimensional cracks are extracted using the p-version finite element method(P-FEM)and the contour integral method.Several numerical experiments,e.g.,crack initiating from the edge of a circular hole under an unidirectional uniform tension and two equal-length,unequal-length hole-edge cracks,respectively,at a rectangular plate,an inclined centered crack under uniaxial tension at a square plate and a pipeline crack model,are used to demonstrate the accuracy and effectiveness of the approaches.SIFs are presented for the effects of various crack lengths and length-width ratio.Numerical results are analyzed and compared with reference solutions and results obtained by the Voronoi cell finite element method,boundary element method,high-order extended finite element method(high-order XFEM)and commercial finite element software ABAQUS in the available literature.Numerical results are in good agreement with the benchmark problems and show faster convergence rate,higher accuracy and better numerical stability.
文摘An elasto-plastie impact model based on the p-version finite element method is presented for the collision protection of ocean and offshore structures. The impact force and responses of the impactor-absorber-structure system can be predicted efficiently and automatically. A cost-effective Cellular Reinforced Concrete Absorber （CRCA） is designed to smooth the impact force and absorb the impact energy. Quasi-static tests show that the concrete absorber has an excellent energy absorbing characteristic. The impact experiment of a scaled offshore oil-piping frame with the proposed concrete absorber is carried out. The simulation results of the elasto-plastie model and the p-version finite element method are in good agreement with the experimental ones. Owing to the plastic deformation of the absorber, the impact force during the impact and responses of the structure are considerably reduced. Further, the proposed impact model and the concrete absorber are applied to the design of collision protection of the sheet-pile groin on the Qiantang River used to weaken the famous Qiantang bore.
