Kukiełka Krzysztof
Department of Production Engineering, Faculty of Mechanical Engineering, Koszalin University of Technology, Racławicka 15-17, 75-620 Koszalin, Poland.
Materials (Basel). 2023 Jun 27;16(13):4647. doi: 10.3390/ma16134647.
The article concerns the application of the FEM method for the prediction of stress and deformation states in a workpiece during the thread rolling process (TR). The analysis covered a new kinematic variant of the TR process in which the basket of the head rotates and is torque-driven, while the workpiece is stationary and the head with the rollers moves axially relative to the workpiece. The TR process was considered as a geometrical and physical non-linear initial and boundary problem with non-linear, moving, and variable in time and space boundary conditions. The boundary conditions in the contact areas of the tool with the workpiece were unknown. An updated Lagrange (UL) description was used to describe the physical phenomena at a typical incremental step. The states of strain and strain rate were described by non-linear relationships without linearization. New discrete systems of motion and deformation equations of the object in the TR were introduced, which take into account the change in the stiffness of the system during the TR process. This equation was solved by the central differences method (explicit). The material parameters were estimated during tensile tests to determine the characteristics of the C45 steel, and a new semi-empirical method was used to determine the relationship yield stress, effective true strain, and effective true strain rate in the thread rolling process. A modified Cowper-Symonds material model was also used to model the displacement process of the wedge on an elastic/visco-plastic body reflecting the TR process. A non-linear dependency of material hardening module depending on strain and strain rate was introduced. To confirm the plane state of deformation and spatial state of stress, an experimental investigation was carried out. The computer models were validated, and a good convergence of the results was obtained. Applications in the ANSYS/LS-Dyna program were developed to simulate the TR process. The developed applications enable a comprehensive time analysis of the states of displacement, strain, and stress occurring in an object consisting of a workpiece (shaft) and a tool (roller) for the case of a plane strain state and a spatial stress state. Exemplary results of numerical analyzes are presented to explain the influence of the friction coefficient on the condition of the thread quality, and the state of deformations and stresses were shown.
本文涉及有限元法在螺纹滚压过程(TR)中预测工件应力和变形状态的应用。分析涵盖了TR过程的一种新运动学变体,其中头部的篮筐旋转并由扭矩驱动,而工件静止,带有滚轮的头部相对于工件轴向移动。TR过程被视为一个几何和物理非线性的初始和边界问题,具有非线性、移动且随时间和空间变化的边界条件。工具与工件接触区域的边界条件未知。采用更新拉格朗日(UL)描述来描述典型增量步的物理现象。应变和应变率状态通过非线性关系描述,无需线性化。引入了TR中物体新的运动和变形离散方程组,该方程组考虑了TR过程中系统刚度的变化。此方程通过中心差分法(显式)求解。在拉伸试验期间估计材料参数以确定C45钢的特性,并采用一种新的半经验方法来确定螺纹滚压过程中的屈服应力、有效真应变和有效真应变率之间的关系。还使用改进的Cowper-Symonds材料模型对反映TR过程的弹性/粘塑性体上楔块的位移过程进行建模。引入了材料硬化模量随应变和应变率的非线性依赖关系。为了确认变形的平面状态和应力的空间状态,进行了实验研究。对计算机模型进行了验证,并获得了良好的结果收敛性。开发了ANSYS/LS-Dyna程序中的应用来模拟TR过程。所开发的应用能够对平面应变状态和空间应力状态下由工件(轴)和工具(滚轮)组成的物体中发生的位移、应变和应力状态进行全面的时间分析。给出了数值分析的示例结果,以解释摩擦系数对螺纹质量状况的影响,并展示了变形和应力状态。