Cross-wedge rolling (CWR) is a metal process of ro ta ry forming. To produce a part, one cylindrical billet should be placed between t wo counterrotating and wedge-shape dies, which move tangentially relative each oth...Cross-wedge rolling (CWR) is a metal process of ro ta ry forming. To produce a part, one cylindrical billet should be placed between t wo counterrotating and wedge-shape dies, which move tangentially relative each other. The billet suffers plastic deformation (essentially, localized compressio n) during its rotation between the rotating dies. Compared to other numerical si mulation methods, the finite element method (FEM) has advantages in solving gene ral problems with complex shapes of the formed parts. In cross-wedge rolling, t here are four stages in the workpiece deformation process, namely knifing, guidi ng, stretching and sizing stage. It is time-consuming and expensive to design t he CWR process by trial and error method. The application of numerical simul ation for the CWR process will help engineers to efficiently improve the process development. Tselikov, Hayama, Jain and Kobayashi, and Higashimo applied the sl ip-line theory in study of CWR process analysis. Zb.pater studied CWR process i ncluding upsetting by upper-bound method. The above numerical simulation were b ased on the two-dimensional plain-strain assumption ignored the metal flow in workpiece axial direction. Therefore, the complex three-dimensional stress and deformation involved in CWR processes were not presented. Compared to other nume rical simulation methods, the finite element method (FEM) has advantages in solv ing general problems with complex shapes of the formed parts. As yet, a few 3-D finite element simulation studies on CWR process have been reported in literatu res. In this paper, the process of cross wedge rolling (CWR) has been simulated and analyzed by 3D rigid-plastic finite element method. Considering the charact eristic of CWR, the static implicit FEM program is selected. The models proposed in this study uses the commercial code DEFORM 3D to simulate the CWR process. T his is an implicit Lagrangian finite element code, which includes many new enhan cements functions. A new method of utilizing multiple processors using the MPI s tandard has been implemented. Automatic switching between the two different defo rmation solvers (Sparse Solver and Conjugate Gradient Solver) has also been impl emented in order to increase the speed of simulations. In this paper, all stages in CWR process are simulated to be able to closely understand and analyze the a ctual CWR process. For simulating all forming stages in CWR process, the dynam ic adaptive remeshing technology for tetrahedral solid elements was applied. T he stress distributions in cross section of forming workpiece are analyzed to in terpret fracture or rarefaction in the center of workpiece. Authors also analyze d the time-torque curve and the laws of load changing.展开更多
The compaction and stress generation on terrain were always investigated based on empirical approaches or testing methods for tire/soil interaction.However,the analysis should be performed for various tires and at dif...The compaction and stress generation on terrain were always investigated based on empirical approaches or testing methods for tire/soil interaction.However,the analysis should be performed for various tires and at different soil strengths.With the increasing capacity of numerical computers and simulation software,finite element modeling of tire/terrain interaction seems a good approach for predicting the effect of change on the parameters.In this work,an elaborated 3D model fully complianning with the geometry of radial tire 115/60R13 was established,using commercial code Solidwork Simulation.The hyper-elastic and incompressible rubber as tire main material was analyzed by Moony-Rivlin model.The Drucker-Prager yield criterion was used to model the soil compaction.Results show that the model realistically predicts the laboratory tests outputs of the modeled tire on the soft soil.展开更多
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.展开更多
高位滑坡对建筑集群的冲击破坏时常导致严重的人员伤亡,基于光滑粒子流体动力学-离散元法-有限元法(smoothed particle hydrodynamics-discrete element method-finite element method,SPH-DEM-FEM)耦合的数值模型,开展了高位滑坡对框...高位滑坡对建筑集群的冲击破坏时常导致严重的人员伤亡,基于光滑粒子流体动力学-离散元法-有限元法(smoothed particle hydrodynamics-discrete element method-finite element method,SPH-DEM-FEM)耦合的数值模型,开展了高位滑坡对框架结构建筑群的冲击过程、建筑结构破坏机理、冲击力时程与框架柱关键点应力和弯矩等动力机制研究。研究结果表明:SPH-DEM-FEM耦合数值方法能够有效地模拟碎石土滑坡中土(SPH)石(DEM)混合物的抛射弹跳、爬高绕流冲击运动过程。考虑了常规建筑垂直、平行于滑坡流向的三排建筑组合布局,位于滑坡近端的纵向排列建筑表现为连续性倾倒破坏,横向排列的建筑则呈现整体倾倒破坏;因前排建筑群对滑坡冲击能量的耗散及滑坡自身摩擦耗能,位于滑坡后端建筑表现为引流面墙体和前排柱发生局部破坏,结构保持稳定,损毁程度依次为上游无建筑缓冲耗能的建筑>有横向排列的建筑>有纵向排列的建筑;纵向、横向排列的建筑冲击力衰减幅度分别31%、21%。横向框架建筑整体倾倒的损毁机制表现为框架柱的直接剪断或节点塑形铰链失效;纵向框架建筑连续性倾倒的损毁机制表现为前排框架柱的失效引起后排框架柱轴向压力和极限弯矩增加,持续冲击荷载超过其极限弯矩致使后排框架柱发生弯曲破坏,最终结构倾倒。系统能量在动能、内能和摩擦耗能间转化,其中摩擦耗能占65.5%,结构耗能占23.6%,动能快速下降与内能急剧增加是建筑破坏的关键特征。展开更多
In order to investigate the process of incremental sheet forming (ISF) through both experimental and numerical approaches, a three-dimensional elasto-plastic finite element model (FEM) was developed to simulate the pr...In order to investigate the process of incremental sheet forming (ISF) through both experimental and numerical approaches, a three-dimensional elasto-plastic finite element model (FEM) was developed to simulate the process and the simulated results were compared with those of experiment. The results of numerical simulations, such as the strain history and distribution, the stress state and distribution, sheet thickness distribution, etc, were discussed in details, and the influences of process parameters on these results were also analyzed. The simulated results of the radial strain and the thickness distribution are in good agreement with experimental results. The simulations reveal that the deformation is localized around the tool and constantly remains close to a plane strain state. With decreasing depth step, increasing tool diameter and wall inclination angle, the axial stress reduces, leading to less thinning and more homogeneous plastic strain and thickness distribution. During ISF, the plastic strain increases stepwise under the action of the tool. Each increase in plastic strain is accompanied by hydrostatic pressure, which explains why obtainable deformation using ISF exceeds the forming limits of conventional sheet forming.展开更多
To explore the influence of double liquid quenching on the cutting performance of the 7A09 aluminum alloy,quasi-static compression and dynamic impact tests were carried out on the 7A09 aluminum alloy after double liqu...To explore the influence of double liquid quenching on the cutting performance of the 7A09 aluminum alloy,quasi-static compression and dynamic impact tests were carried out on the 7A09 aluminum alloy after double liquid quenching using an MTS810.23 universal testing machine and split-Hopkinson pressure bar(SHPB).The experimental data were fitted to obtain the Johnson–Cook constitutive model parameters of the alloy.Simulations of the machining process were carried out using the Deform-3D finite element software.The results showed that the rheological stress increased with the increase in strain rate and the decrease in temperature.The increase in the cutting speed and feed caused the cutting temperature to rise sharply,whereas the influence of the cutting amount on the cutting temperature was weak.Because of the presence of chip nodules,there was extremum in the cutting force vs cutting speed curves.The increase in the feed and cutting depth increased the cutting area Ac,so the cutting force also increased.The simulation results were verified by experiments.The simulation predictions were in good agreement with the test values,and the cutting force and temperature variations with the cutting parameters were the same.Thus,the correctness of the 7A09 aluminum alloy finite element model was verified.展开更多
Being aimed at the inside wall wrinkling and sinking phenomenon of palladium-yttrium alloy thin wall spiral tube used for preparation of high purity hydrogen, extraction of hydrogen isotope, and purification and separ...Being aimed at the inside wall wrinkling and sinking phenomenon of palladium-yttrium alloy thin wall spiral tube used for preparation of high purity hydrogen, extraction of hydrogen isotope, and purification and separation of hydrogen in the winding process, this article analyzed the reasons for above phenomena, established a numerical simulation model of winding process of above tube, using elastic-plastic Finite Element method analyzed the max. tensile stress and max. compression stress and their locations, thereby provides a theory base for the control of working forming course of thin wall spiral tube.展开更多
In this paper, the internal fluid motion of a jet system is described by the Navier Stokes mechanics equations. For the simulation of the motion, the penalty function finite element method is used, and the velocity ve...In this paper, the internal fluid motion of a jet system is described by the Navier Stokes mechanics equations. For the simulation of the motion, the penalty function finite element method is used, and the velocity vectors and stream function curves are obtained. Using the Prandtl theory, this paper derives the free jet velocity and the jet bunch width in a half-space, the latter of which is amended by experiment. The results obtained in this paper are applied to micro-type high pressure water jet cleaner and the ejector of rocket engine.展开更多
针对二维介质目标的电磁成像问题,将正余弦算法(Sine Cosine Algorithm,SCA)与有限元方法(Finite Element Method,FEM)和不变性测试方程(Measured Equation of Invariance,MEI)进行结合提出一种新的成像方法。将FEM与MEI进行结合求解二...针对二维介质目标的电磁成像问题,将正余弦算法(Sine Cosine Algorithm,SCA)与有限元方法(Finite Element Method,FEM)和不变性测试方程(Measured Equation of Invariance,MEI)进行结合提出一种新的成像方法。将FEM与MEI进行结合求解二维介质目标的电磁散射正问题,即求解Helmholtz方程。其中,MEI保证边界截断的精度,FEM适用于复杂介质目标的准确模拟。对于电磁散射逆问题,引入SCA并加以改进提出一种新的重构方法。该方法采用等效原理与格林函数的渐近式求得远区散射场,以测量的散射场和计算的散射场最大偏差为目标函数,采用改进的SCA优化介质参数,使目标函数达到最小值,以此重构散射体。为提高计算效率,采用MPI算法进行并行计算。文中采用基准函数展示了改进的SCA算法的快速收敛性,并采用非规则的均匀介质柱目标验证了成像方法的正确性。展开更多
针对高速电机和飞轮储能系统等对空间利用率要求较高的场合,提出一种新型异极径向混合磁轴承(heteropolar radial hybrid magnetic bearing,HRHMB).首先,建立该磁轴承的等效磁路模型,通过解析磁场得出其电流刚度、位移刚度及电磁力,并...针对高速电机和飞轮储能系统等对空间利用率要求较高的场合,提出一种新型异极径向混合磁轴承(heteropolar radial hybrid magnetic bearing,HRHMB).首先,建立该磁轴承的等效磁路模型,通过解析磁场得出其电流刚度、位移刚度及电磁力,并通过有限元仿真验证其有效性;然后,在相同约束条件下与传统偏置磁轴承进行对比,分析磁轴承的刚度特性和空间利用率;最后,通过有限元仿真研究新型磁轴承径向两自由度间的电磁力耦合,并与传统磁轴承进行对比.研究结果表明:在相同承载力等约束条件下,该新型磁轴承的体积仅为传统磁轴承的0.87倍,其电磁力在控制电流和转子位移影响下的相对误差值为6.5%,而传统磁轴承的电磁力相对误差为13.6%,表明新型磁轴承径向两自由度的电磁力耦合小于传统磁轴承,解耦效果良好.展开更多
文摘Cross-wedge rolling (CWR) is a metal process of ro ta ry forming. To produce a part, one cylindrical billet should be placed between t wo counterrotating and wedge-shape dies, which move tangentially relative each other. The billet suffers plastic deformation (essentially, localized compressio n) during its rotation between the rotating dies. Compared to other numerical si mulation methods, the finite element method (FEM) has advantages in solving gene ral problems with complex shapes of the formed parts. In cross-wedge rolling, t here are four stages in the workpiece deformation process, namely knifing, guidi ng, stretching and sizing stage. It is time-consuming and expensive to design t he CWR process by trial and error method. The application of numerical simul ation for the CWR process will help engineers to efficiently improve the process development. Tselikov, Hayama, Jain and Kobayashi, and Higashimo applied the sl ip-line theory in study of CWR process analysis. Zb.pater studied CWR process i ncluding upsetting by upper-bound method. The above numerical simulation were b ased on the two-dimensional plain-strain assumption ignored the metal flow in workpiece axial direction. Therefore, the complex three-dimensional stress and deformation involved in CWR processes were not presented. Compared to other nume rical simulation methods, the finite element method (FEM) has advantages in solv ing general problems with complex shapes of the formed parts. As yet, a few 3-D finite element simulation studies on CWR process have been reported in literatu res. In this paper, the process of cross wedge rolling (CWR) has been simulated and analyzed by 3D rigid-plastic finite element method. Considering the charact eristic of CWR, the static implicit FEM program is selected. The models proposed in this study uses the commercial code DEFORM 3D to simulate the CWR process. T his is an implicit Lagrangian finite element code, which includes many new enhan cements functions. A new method of utilizing multiple processors using the MPI s tandard has been implemented. Automatic switching between the two different defo rmation solvers (Sparse Solver and Conjugate Gradient Solver) has also been impl emented in order to increase the speed of simulations. In this paper, all stages in CWR process are simulated to be able to closely understand and analyze the a ctual CWR process. For simulating all forming stages in CWR process, the dynam ic adaptive remeshing technology for tetrahedral solid elements was applied. T he stress distributions in cross section of forming workpiece are analyzed to in terpret fracture or rarefaction in the center of workpiece. Authors also analyze d the time-torque curve and the laws of load changing.
文摘The compaction and stress generation on terrain were always investigated based on empirical approaches or testing methods for tire/soil interaction.However,the analysis should be performed for various tires and at different soil strengths.With the increasing capacity of numerical computers and simulation software,finite element modeling of tire/terrain interaction seems a good approach for predicting the effect of change on the parameters.In this work,an elaborated 3D model fully complianning with the geometry of radial tire 115/60R13 was established,using commercial code Solidwork Simulation.The hyper-elastic and incompressible rubber as tire main material was analyzed by Moony-Rivlin model.The Drucker-Prager yield criterion was used to model the soil compaction.Results show that the model realistically predicts the laboratory tests outputs of the modeled tire on the soft soil.
基金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.
文摘高位滑坡对建筑集群的冲击破坏时常导致严重的人员伤亡,基于光滑粒子流体动力学-离散元法-有限元法(smoothed particle hydrodynamics-discrete element method-finite element method,SPH-DEM-FEM)耦合的数值模型,开展了高位滑坡对框架结构建筑群的冲击过程、建筑结构破坏机理、冲击力时程与框架柱关键点应力和弯矩等动力机制研究。研究结果表明:SPH-DEM-FEM耦合数值方法能够有效地模拟碎石土滑坡中土(SPH)石(DEM)混合物的抛射弹跳、爬高绕流冲击运动过程。考虑了常规建筑垂直、平行于滑坡流向的三排建筑组合布局,位于滑坡近端的纵向排列建筑表现为连续性倾倒破坏,横向排列的建筑则呈现整体倾倒破坏;因前排建筑群对滑坡冲击能量的耗散及滑坡自身摩擦耗能,位于滑坡后端建筑表现为引流面墙体和前排柱发生局部破坏,结构保持稳定,损毁程度依次为上游无建筑缓冲耗能的建筑>有横向排列的建筑>有纵向排列的建筑;纵向、横向排列的建筑冲击力衰减幅度分别31%、21%。横向框架建筑整体倾倒的损毁机制表现为框架柱的直接剪断或节点塑形铰链失效;纵向框架建筑连续性倾倒的损毁机制表现为前排框架柱的失效引起后排框架柱轴向压力和极限弯矩增加,持续冲击荷载超过其极限弯矩致使后排框架柱发生弯曲破坏,最终结构倾倒。系统能量在动能、内能和摩擦耗能间转化,其中摩擦耗能占65.5%,结构耗能占23.6%,动能快速下降与内能急剧增加是建筑破坏的关键特征。
基金Project(50175034) supported by the National Natural Science Foundation of China
文摘In order to investigate the process of incremental sheet forming (ISF) through both experimental and numerical approaches, a three-dimensional elasto-plastic finite element model (FEM) was developed to simulate the process and the simulated results were compared with those of experiment. The results of numerical simulations, such as the strain history and distribution, the stress state and distribution, sheet thickness distribution, etc, were discussed in details, and the influences of process parameters on these results were also analyzed. The simulated results of the radial strain and the thickness distribution are in good agreement with experimental results. The simulations reveal that the deformation is localized around the tool and constantly remains close to a plane strain state. With decreasing depth step, increasing tool diameter and wall inclination angle, the axial stress reduces, leading to less thinning and more homogeneous plastic strain and thickness distribution. During ISF, the plastic strain increases stepwise under the action of the tool. Each increase in plastic strain is accompanied by hydrostatic pressure, which explains why obtainable deformation using ISF exceeds the forming limits of conventional sheet forming.
基金Projects(51575289,51705270)supported by the National Natural Science Foundation of ChinaProject(2019GHY112068)supported by the Key Research and Development of Shandong,China
文摘To explore the influence of double liquid quenching on the cutting performance of the 7A09 aluminum alloy,quasi-static compression and dynamic impact tests were carried out on the 7A09 aluminum alloy after double liquid quenching using an MTS810.23 universal testing machine and split-Hopkinson pressure bar(SHPB).The experimental data were fitted to obtain the Johnson–Cook constitutive model parameters of the alloy.Simulations of the machining process were carried out using the Deform-3D finite element software.The results showed that the rheological stress increased with the increase in strain rate and the decrease in temperature.The increase in the cutting speed and feed caused the cutting temperature to rise sharply,whereas the influence of the cutting amount on the cutting temperature was weak.Because of the presence of chip nodules,there was extremum in the cutting force vs cutting speed curves.The increase in the feed and cutting depth increased the cutting area Ac,so the cutting force also increased.The simulation results were verified by experiments.The simulation predictions were in good agreement with the test values,and the cutting force and temperature variations with the cutting parameters were the same.Thus,the correctness of the 7A09 aluminum alloy finite element model was verified.
文摘Being aimed at the inside wall wrinkling and sinking phenomenon of palladium-yttrium alloy thin wall spiral tube used for preparation of high purity hydrogen, extraction of hydrogen isotope, and purification and separation of hydrogen in the winding process, this article analyzed the reasons for above phenomena, established a numerical simulation model of winding process of above tube, using elastic-plastic Finite Element method analyzed the max. tensile stress and max. compression stress and their locations, thereby provides a theory base for the control of working forming course of thin wall spiral tube.
文摘In this paper, the internal fluid motion of a jet system is described by the Navier Stokes mechanics equations. For the simulation of the motion, the penalty function finite element method is used, and the velocity vectors and stream function curves are obtained. Using the Prandtl theory, this paper derives the free jet velocity and the jet bunch width in a half-space, the latter of which is amended by experiment. The results obtained in this paper are applied to micro-type high pressure water jet cleaner and the ejector of rocket engine.
文摘针对二维介质目标的电磁成像问题,将正余弦算法(Sine Cosine Algorithm,SCA)与有限元方法(Finite Element Method,FEM)和不变性测试方程(Measured Equation of Invariance,MEI)进行结合提出一种新的成像方法。将FEM与MEI进行结合求解二维介质目标的电磁散射正问题,即求解Helmholtz方程。其中,MEI保证边界截断的精度,FEM适用于复杂介质目标的准确模拟。对于电磁散射逆问题,引入SCA并加以改进提出一种新的重构方法。该方法采用等效原理与格林函数的渐近式求得远区散射场,以测量的散射场和计算的散射场最大偏差为目标函数,采用改进的SCA优化介质参数,使目标函数达到最小值,以此重构散射体。为提高计算效率,采用MPI算法进行并行计算。文中采用基准函数展示了改进的SCA算法的快速收敛性,并采用非规则的均匀介质柱目标验证了成像方法的正确性。
文摘针对高速电机和飞轮储能系统等对空间利用率要求较高的场合,提出一种新型异极径向混合磁轴承(heteropolar radial hybrid magnetic bearing,HRHMB).首先,建立该磁轴承的等效磁路模型,通过解析磁场得出其电流刚度、位移刚度及电磁力,并通过有限元仿真验证其有效性;然后,在相同约束条件下与传统偏置磁轴承进行对比,分析磁轴承的刚度特性和空间利用率;最后,通过有限元仿真研究新型磁轴承径向两自由度间的电磁力耦合,并与传统磁轴承进行对比.研究结果表明:在相同承载力等约束条件下,该新型磁轴承的体积仅为传统磁轴承的0.87倍,其电磁力在控制电流和转子位移影响下的相对误差值为6.5%,而传统磁轴承的电磁力相对误差为13.6%,表明新型磁轴承径向两自由度的电磁力耦合小于传统磁轴承,解耦效果良好.
