The discrete ordinates(S N)method requires numerous angular unknowns to achieve the desired accu-racy for shielding calculations involving strong anisotropy.Our objective is to develop an angular adaptive algorithm in...The discrete ordinates(S N)method requires numerous angular unknowns to achieve the desired accu-racy for shielding calculations involving strong anisotropy.Our objective is to develop an angular adaptive algorithm in the S N method to automatically optimize the angular distribution and minimize angular discretization errors with lower expenses.The proposed method enables linear dis-continuous finite element quadrature sets over an icosahe-dron to vary their quadrature orders in a one-twentieth sphere so that fine resolutions can be applied to the angular domains that are important.An error estimation that operates in conjunction with the spherical harmonics method is developed to determine the locations where more refinement is required.The adaptive quadrature sets are applied to three duct problems,including the Kobayashi benchmarks and the IRI-TUB research reactor,which emphasize the ability of this method to resolve neutron streaming through ducts with voids.The results indicate that the performance of the adaptive method is more effi-cient than that of uniform quadrature sets for duct transport problems.Our adaptive method offers an appropriate placement of angular unknowns to accurately integrate angular fluxes while reducing the computational costs in terms of unknowns and run times.展开更多
Determination of ballistic performance of an armor solution is a complicated task and evolved significantly with the application of finite element methods(FEM) in this research field.The traditional armor design studi...Determination of ballistic performance of an armor solution is a complicated task and evolved significantly with the application of finite element methods(FEM) in this research field.The traditional armor design studies performed with FEM requires sophisticated procedures and intensive computational effort,therefore simpler and accurate numerical approaches are always worthwhile to decrease armor development time.This study aims to apply a hybrid method using FEM simulation and artificial neural network(ANN) analysis to approximate ballistic limit thickness for armor steels.To achieve this objective,a predictive model based on the artificial neural networks is developed to determine ballistic resistance of high hardness armor steels against 7.62 mm armor piercing ammunition.In this methodology,the FEM simulations are used to create training cases for Multilayer Perceptron(MLP) three layer networks.In order to validate FE simulation methodology,ballistic shot tests on 20 mm thickness target were performed according to standard Stanag 4569.Afterwards,the successfully trained ANN(s) is used to predict the ballistic limit thickness of 500 HB high hardness steel armor.Results show that even with limited number of data,FEM-ANN approach can be used to predict ballistic penetration depth with adequate accuracy.展开更多
We review recent advances in the finite element method (FEM) simulations of interactions between waves and structures. Our focus is on the potential theory with the fully nonlinear or second-order boundary condition. ...We review recent advances in the finite element method (FEM) simulations of interactions between waves and structures. Our focus is on the potential theory with the fully nonlinear or second-order boundary condition. The present paper has six sections. A review of previous work on interactions between waves and ocean structures is presented in Section one. Section two gives the mathematical formulation. In Section three, the finite element discretization, mesh generation and the finite element linear system solution methods are described. Section four presents numerical methods including time marching schemes, computation of velocity, remeshing and smoothing techniques and numerical radiation conditions. The application of the FEM to the wave-structure interactions are presented in Section five followed by the concluding remarks in Section six.展开更多
The complex structure and strong heterogeneity of advanced nuclear reactor systems pose challenges for high-fidelity neutron-shielding calculations. Unstructured meshes exhibit strong geometric adaptability and can ov...The complex structure and strong heterogeneity of advanced nuclear reactor systems pose challenges for high-fidelity neutron-shielding calculations. Unstructured meshes exhibit strong geometric adaptability and can overcome the deficiencies of conventionally structured meshes in complex geometry modeling. A multithreaded parallel upwind sweep algorithm for S_(N) transport was proposed to achieve a more accurate geometric description and improve the computational efficiency. The spatial variables were discretized using the standard discontinuous Galerkin finite-element method. The angular flux transmission between neighboring meshes was handled using an upwind scheme. In addition, a combination of a mesh transport sweep and angular iterations was realized using a multithreaded parallel technique. The algorithm was implemented in the 2D/3D S_(N) transport code ThorSNIPE, and numerical evaluations were conducted using three typical benchmark problems:IAEA, Kobayashi-3i, and VENUS-3. These numerical results indicate that the multithreaded parallel upwind sweep algorithm can achieve high computational efficiency. ThorSNIPE, with a multithreaded parallel upwind sweep algorithm, has good reliability, stability, and high efficiency, making it suitable for complex shielding calculations.展开更多
Base d on fluid velocity potential, an ALE finite element formulation for the analysi s of nonlinear sloshing problems has been developed. The ALE kinemat ical description is introduced to move the computational mesh...Base d on fluid velocity potential, an ALE finite element formulation for the analysi s of nonlinear sloshing problems has been developed. The ALE kinemat ical description is introduced to move the computational mesh independently of f luid motion, and the container fixed noninertial coordinate system is employed to establish the governing equations so that the mesh is needed to be updated in this coordinate system only. This leads to a very simple mesh moving algorithm which makes it easy to trace the motion of the moving boundaries and the free su rface without producing undesirable distortion of the computational mesh. The fi nite element method and finite difference method are used spacewise and timewise , respectively. A numerical example involving either forced horizontal oscillati on or forced pitching oscillation of the fluid filled container is presented to illustrate the effectiveness and the robustness of the method. In additi on, this work can be extended for the fluid structure interaction problems.展开更多
A contact bolt model is proposed as a new modeling technique to investigate the complex structure with bolted joints for modal analysis and compared with the coupled bolt model, and the test results are given. Among t...A contact bolt model is proposed as a new modeling technique to investigate the complex structure with bolted joints for modal analysis and compared with the coupled bolt model, and the test results are given. Among these models, the coupled bolt model provides the best accurate responses compared with the experimental results. The contact bolt model shows the best effectiveness and usefulness in view of operational time. The bolt models proposed in this study are adopted for a dynamic characteristic analysis of a large diesel engine consisting of several parts which are connected by many bolts. The dynamic behavior of the entire engine structure was investigated by experiment. The coupled bolt model and the contact bolt model were applied to model the assembly of engine with high preload. The experimental results are in good agreement with the finite element method （FEM） results. Compared with the other models, the contact bolt model presented in this paper is more effective and useful in view of operational time and experience of analysts.展开更多
With the application of Hammer integral formulas of a continuous function on a triangular element, the numerical integral formulas of some discrete functions on the element are derived by means of decomposition and re...With the application of Hammer integral formulas of a continuous function on a triangular element, the numerical integral formulas of some discrete functions on the element are derived by means of decomposition and recombination of base functions. Hammer integral formulas are the special examples of those of the paper.展开更多
结合有限元离散元方法(finite-discrete element method,FDEM),对计算流体力学(computational fluid dynamics,CFD)软件FLOW-3D进行二次开发,建立了基于CFD-DEM的流固耦合模型,模拟了多块石入水、沉降以及触底的动力过程,分析了不同块...结合有限元离散元方法(finite-discrete element method,FDEM),对计算流体力学(computational fluid dynamics,CFD)软件FLOW-3D进行二次开发,建立了基于CFD-DEM的流固耦合模型,模拟了多块石入水、沉降以及触底的动力过程,分析了不同块石等效直径、形状和入水速度对触底速度和反力的影响。研究发现,块石入水后速度迅速减小,并逐渐趋于定值,随后做动态平衡沉降运动,直至与底面发生碰撞。块石抛填的触底速度随等效直径的增大而增大,不同等效直径下球形块石触底速度绝对值最大,其次是纺锤形块石,最小为圆盘形块石。最大触底反力也随等效直径的增大而增大,成非线性关系,通过拟合得到了触底反力的经验公式。显著性分析结果表明,块石等效直径对触底反力影响最大,其次是块石形状,最小为入水速度。展开更多
Impact dynamics of flexible solids is important in engineering practice. Obtaining its dynamic response is a challenging task and usually achieved by numerical methods. The objectives of the study are twofold. Firstly...Impact dynamics of flexible solids is important in engineering practice. Obtaining its dynamic response is a challenging task and usually achieved by numerical methods. The objectives of the study are twofold. Firstly, the discrete singular convolution （DSC） is used for the first time to analyze the impact dynamics. Secondly, the efficiency of various numerical methods for dynamic analysis is explored via an example of a flexible rod hit by a rigid ball. Three numerical methods, including the conventional finite element （FE） method, the DSC algorithm, and the spectral finite element （SFE） method, and one proposed modeling strategy, the improved spectral finite element （ISFE） method, are involved. Numerical results are compared with the known analytical solutions to show their efficiency. It is demonstrated that the proposed ISFE modeling strategy with a proper length of con- ventional FE yields the most accurate contact stress among the four investigated models. It is also found that the DSC algorithm is an alternative method for collision problems.展开更多
Residual stresses can have a strong effect on the usability of machined parts,and the X-ray diffraction(XRD)measuring equipment,which is commonly used to measure residual stresses,is very expensive.This paper presents...Residual stresses can have a strong effect on the usability of machined parts,and the X-ray diffraction(XRD)measuring equipment,which is commonly used to measure residual stresses,is very expensive.This paper presents a method of measuring the residual stresses induced by boring in the internal surface of a tube with much cheaper equipment.The method,called the strain-based method is mainly based on the strains measured on the external surface of the tube.It is proposed on the basis of the very long tube assumption.The finite element method(FEM)analysis is thus used to validate the length of the tube.Guided by the FEM results,an appropriate length of the tube is chosen,and the residual stresses are obtained from both the strain-based method and the XRD method.Stress profiles obtained from both two methods are compared.The comparison result indicates that the profiles of the two methods agree well with each other.Therefore,it can be concluded that the accuracy of the strain-based method is high enough,and it can be applied to residual stress measurement in practice.展开更多
Fatigue strength assessment of a horizontal axis wind turbine(HAWT)composite blade is considered.Fatigue load cases are identified,and loads are calculated by the GH Bladed software which is specified at the IEC61400 ...Fatigue strength assessment of a horizontal axis wind turbine(HAWT)composite blade is considered.Fatigue load cases are identified,and loads are calculated by the GH Bladed software which is specified at the IEC61400 international specification and GL(Germanisher Lloyd)regulations for the wind energy conversion system.Stress analysis is performed with a 3-D finite element method(FEM).Considering Saint-Venant′s principle,a uniform cross section FEM model is built at each critical zone.Stress transformation matrixes(STM)are set up by applied six unit load components on the FEM model separately.STM can be used to convert the external load into stresses in the linear elastic range.The main material of composite wind turbine blade is fiber reinforced plastics(FRP).In order to evaluate the degree of fatigue damage of FRP,the stresses of fiber direction are extracted and the well-known strength criterion-Puck theory is used.The total fatigue damage of each laminate on the critical point is counted by the rain-flow counting method and Miner′s damage law based on general S-N curves.Several sections of a 45.3mblade of a 2 MW wind turbine are studied using the fatigue evaluation method.The performance of this method is compared with far more costly business software FOCUS.The results show that the fatigue damage of multi-axis FRP can be assessed conveniently by the FEM-STM method.And the proposed method gives a reliable and efficient method to analyze the fatigue damage of slender composite structure with variable cross-sections.展开更多
基金supported by the National Natural Science Foundation of China(No.11975097)the Fundamental Research Funds for the Central Universities(No.2019MS038).
文摘The discrete ordinates(S N)method requires numerous angular unknowns to achieve the desired accu-racy for shielding calculations involving strong anisotropy.Our objective is to develop an angular adaptive algorithm in the S N method to automatically optimize the angular distribution and minimize angular discretization errors with lower expenses.The proposed method enables linear dis-continuous finite element quadrature sets over an icosahe-dron to vary their quadrature orders in a one-twentieth sphere so that fine resolutions can be applied to the angular domains that are important.An error estimation that operates in conjunction with the spherical harmonics method is developed to determine the locations where more refinement is required.The adaptive quadrature sets are applied to three duct problems,including the Kobayashi benchmarks and the IRI-TUB research reactor,which emphasize the ability of this method to resolve neutron streaming through ducts with voids.The results indicate that the performance of the adaptive method is more effi-cient than that of uniform quadrature sets for duct transport problems.Our adaptive method offers an appropriate placement of angular unknowns to accurately integrate angular fluxes while reducing the computational costs in terms of unknowns and run times.
基金Otokar Otomotiv ve Savunma Sanayi A.S. for the financial support
文摘Determination of ballistic performance of an armor solution is a complicated task and evolved significantly with the application of finite element methods(FEM) in this research field.The traditional armor design studies performed with FEM requires sophisticated procedures and intensive computational effort,therefore simpler and accurate numerical approaches are always worthwhile to decrease armor development time.This study aims to apply a hybrid method using FEM simulation and artificial neural network(ANN) analysis to approximate ballistic limit thickness for armor steels.To achieve this objective,a predictive model based on the artificial neural networks is developed to determine ballistic resistance of high hardness armor steels against 7.62 mm armor piercing ammunition.In this methodology,the FEM simulations are used to create training cases for Multilayer Perceptron(MLP) three layer networks.In order to validate FE simulation methodology,ballistic shot tests on 20 mm thickness target were performed according to standard Stanag 4569.Afterwards,the successfully trained ANN(s) is used to predict the ballistic limit thickness of 500 HB high hardness steel armor.Results show that even with limited number of data,FEM-ANN approach can be used to predict ballistic penetration depth with adequate accuracy.
文摘We review recent advances in the finite element method (FEM) simulations of interactions between waves and structures. Our focus is on the potential theory with the fully nonlinear or second-order boundary condition. The present paper has six sections. A review of previous work on interactions between waves and ocean structures is presented in Section one. Section two gives the mathematical formulation. In Section three, the finite element discretization, mesh generation and the finite element linear system solution methods are described. Section four presents numerical methods including time marching schemes, computation of velocity, remeshing and smoothing techniques and numerical radiation conditions. The application of the FEM to the wave-structure interactions are presented in Section five followed by the concluding remarks in Section six.
文摘The complex structure and strong heterogeneity of advanced nuclear reactor systems pose challenges for high-fidelity neutron-shielding calculations. Unstructured meshes exhibit strong geometric adaptability and can overcome the deficiencies of conventionally structured meshes in complex geometry modeling. A multithreaded parallel upwind sweep algorithm for S_(N) transport was proposed to achieve a more accurate geometric description and improve the computational efficiency. The spatial variables were discretized using the standard discontinuous Galerkin finite-element method. The angular flux transmission between neighboring meshes was handled using an upwind scheme. In addition, a combination of a mesh transport sweep and angular iterations was realized using a multithreaded parallel technique. The algorithm was implemented in the 2D/3D S_(N) transport code ThorSNIPE, and numerical evaluations were conducted using three typical benchmark problems:IAEA, Kobayashi-3i, and VENUS-3. These numerical results indicate that the multithreaded parallel upwind sweep algorithm can achieve high computational efficiency. ThorSNIPE, with a multithreaded parallel upwind sweep algorithm, has good reliability, stability, and high efficiency, making it suitable for complex shielding calculations.
文摘Base d on fluid velocity potential, an ALE finite element formulation for the analysi s of nonlinear sloshing problems has been developed. The ALE kinemat ical description is introduced to move the computational mesh independently of f luid motion, and the container fixed noninertial coordinate system is employed to establish the governing equations so that the mesh is needed to be updated in this coordinate system only. This leads to a very simple mesh moving algorithm which makes it easy to trace the motion of the moving boundaries and the free su rface without producing undesirable distortion of the computational mesh. The fi nite element method and finite difference method are used spacewise and timewise , respectively. A numerical example involving either forced horizontal oscillati on or forced pitching oscillation of the fluid filled container is presented to illustrate the effectiveness and the robustness of the method. In additi on, this work can be extended for the fluid structure interaction problems.
基金Sponsored by the Ministerial Level Foundation(40402020105)
文摘A contact bolt model is proposed as a new modeling technique to investigate the complex structure with bolted joints for modal analysis and compared with the coupled bolt model, and the test results are given. Among these models, the coupled bolt model provides the best accurate responses compared with the experimental results. The contact bolt model shows the best effectiveness and usefulness in view of operational time. The bolt models proposed in this study are adopted for a dynamic characteristic analysis of a large diesel engine consisting of several parts which are connected by many bolts. The dynamic behavior of the entire engine structure was investigated by experiment. The coupled bolt model and the contact bolt model were applied to model the assembly of engine with high preload. The experimental results are in good agreement with the finite element method （FEM） results. Compared with the other models, the contact bolt model presented in this paper is more effective and useful in view of operational time and experience of analysts.
文摘With the application of Hammer integral formulas of a continuous function on a triangular element, the numerical integral formulas of some discrete functions on the element are derived by means of decomposition and recombination of base functions. Hammer integral formulas are the special examples of those of the paper.
文摘结合有限元离散元方法(finite-discrete element method,FDEM),对计算流体力学(computational fluid dynamics,CFD)软件FLOW-3D进行二次开发,建立了基于CFD-DEM的流固耦合模型,模拟了多块石入水、沉降以及触底的动力过程,分析了不同块石等效直径、形状和入水速度对触底速度和反力的影响。研究发现,块石入水后速度迅速减小,并逐渐趋于定值,随后做动态平衡沉降运动,直至与底面发生碰撞。块石抛填的触底速度随等效直径的增大而增大,不同等效直径下球形块石触底速度绝对值最大,其次是纺锤形块石,最小为圆盘形块石。最大触底反力也随等效直径的增大而增大,成非线性关系,通过拟合得到了触底反力的经验公式。显著性分析结果表明,块石等效直径对触底反力影响最大,其次是块石形状,最小为入水速度。
基金Supported by the National Natural Science Foundation of China(50830201)the Priority Academic Program Development of Jiangsu Higher Education Institutions~~
文摘Impact dynamics of flexible solids is important in engineering practice. Obtaining its dynamic response is a challenging task and usually achieved by numerical methods. The objectives of the study are twofold. Firstly, the discrete singular convolution （DSC） is used for the first time to analyze the impact dynamics. Secondly, the efficiency of various numerical methods for dynamic analysis is explored via an example of a flexible rod hit by a rigid ball. Three numerical methods, including the conventional finite element （FE） method, the DSC algorithm, and the spectral finite element （SFE） method, and one proposed modeling strategy, the improved spectral finite element （ISFE） method, are involved. Numerical results are compared with the known analytical solutions to show their efficiency. It is demonstrated that the proposed ISFE modeling strategy with a proper length of con- ventional FE yields the most accurate contact stress among the four investigated models. It is also found that the DSC algorithm is an alternative method for collision problems.
基金Supported by the National Defense Program of China(C152012C002)the Specialized Research Fund for the Doctoral Program of Higher Education of China(20123218120025)
文摘Residual stresses can have a strong effect on the usability of machined parts,and the X-ray diffraction(XRD)measuring equipment,which is commonly used to measure residual stresses,is very expensive.This paper presents a method of measuring the residual stresses induced by boring in the internal surface of a tube with much cheaper equipment.The method,called the strain-based method is mainly based on the strains measured on the external surface of the tube.It is proposed on the basis of the very long tube assumption.The finite element method(FEM)analysis is thus used to validate the length of the tube.Guided by the FEM results,an appropriate length of the tube is chosen,and the residual stresses are obtained from both the strain-based method and the XRD method.Stress profiles obtained from both two methods are compared.The comparison result indicates that the profiles of the two methods agree well with each other.Therefore,it can be concluded that the accuracy of the strain-based method is high enough,and it can be applied to residual stress measurement in practice.
基金supported jointly by the National Basic Research Program of China(″973″Program)(No.2014CB046200)the Natural Science Foundation of Jiangsu Province(No.BK2014059)+1 种基金the Priority Academic Program Development of Jiangsu Higher Education Institutionsthe National Natural Science Foundation of China(No.11172135)
文摘Fatigue strength assessment of a horizontal axis wind turbine(HAWT)composite blade is considered.Fatigue load cases are identified,and loads are calculated by the GH Bladed software which is specified at the IEC61400 international specification and GL(Germanisher Lloyd)regulations for the wind energy conversion system.Stress analysis is performed with a 3-D finite element method(FEM).Considering Saint-Venant′s principle,a uniform cross section FEM model is built at each critical zone.Stress transformation matrixes(STM)are set up by applied six unit load components on the FEM model separately.STM can be used to convert the external load into stresses in the linear elastic range.The main material of composite wind turbine blade is fiber reinforced plastics(FRP).In order to evaluate the degree of fatigue damage of FRP,the stresses of fiber direction are extracted and the well-known strength criterion-Puck theory is used.The total fatigue damage of each laminate on the critical point is counted by the rain-flow counting method and Miner′s damage law based on general S-N curves.Several sections of a 45.3mblade of a 2 MW wind turbine are studied using the fatigue evaluation method.The performance of this method is compared with far more costly business software FOCUS.The results show that the fatigue damage of multi-axis FRP can be assessed conveniently by the FEM-STM method.And the proposed method gives a reliable and efficient method to analyze the fatigue damage of slender composite structure with variable cross-sections.
