首页 | 本学科首页   官方微博 | 高级检索  
     

接触问题的三角形载荷离散FFT加速算法
引用本文:陈楠,朱凯,蒋志桢,龚诗雨,李璞,金晓清. 接触问题的三角形载荷离散FFT加速算法[J]. 重庆大学学报(自然科学版), 2024, 47(2): 95-105
作者姓名:陈楠  朱凯  蒋志桢  龚诗雨  李璞  金晓清
作者单位:1.重庆大学,航空航天学院,重庆 400044;2.重庆大学,机械传动国家重点实验室,重庆 400044;3.哈尔滨工业大学 理学院, 广东 深圳 518055
基金项目:国家自然科学基金资助项目(52205192,51875059,11932004);重庆市科技计划资助项目(cstc2020jcyj-msxmX0850)。
摘    要:接触问题控制方程的有效求解,往往涉及到复杂的数学理论知识,而在实际工程应用中,接触应力分布又具有高度随机性。为高效快速求解任意载荷分布下固体的接触响应,基于三角形载荷离散单元,嵌入离散卷积快速傅里叶变换(DC-FFT)算法,提供了一种高精度、高可靠度的计算方法。相比于通常采用的分段均布载荷离散方法,三角形单元的解析求解略显复杂,但能更好地模拟接触载荷任意分布的特性,对于接触边缘处载荷由零递增或递减为零的情况,也可以予以充分考虑。为优化三角形载荷离散单元的求解方法,根据接触影响系数矩阵的“激励-响应”特性,推导了三角形载荷单元和均布载荷单元作用下的应力分量解析解。通过构造包含影响系数矩阵的离散卷积形式应力解,将某一目标节点在所有载荷单元作用下,重复度极高的矩阵运算叠加过程,采用DC-FFT算法来简化加速计算。通过程序编程计算,分析验证了所提出算法的精确度和高效性。

关 键 词:三角形单元  接触应力  DC-FFT  数值解  应力场
收稿时间:2023-02-16

FFT acceleration algorithm for contact problems based on triangular element discretization
CHEN Nan,ZHU Kai,JIANG Zhizhen,GONG Shiyu,LI Pu,JIN Xiaoqing. FFT acceleration algorithm for contact problems based on triangular element discretization[J]. Journal of Chongqing University(Natural Science Edition), 2024, 47(2): 95-105
Authors:CHEN Nan  ZHU Kai  JIANG Zhizhen  GONG Shiyu  LI Pu  JIN Xiaoqing
Affiliation:1. College of Aerospace Engineering; Chongqing University, Chongqing 400044, P. R. China;2.State Key Laboratory of Mechanical Transmissions, Chongqing University, Chongqing 400044, P. R. China;3.School of Science, Harbin Institute of Technology, Shenzhen 518055, Guangdong, P. R. China
Abstract:Effectively solving the governing equations for contact problems often involves complex mathematical theory, while the distribution of contact stress is highly random in practical engineering applications. This study proposes a novel algorithm based on the triangular load discrete element and the discrete convolution fast Fourier transform (DC-FFT) algorithm. This algorithm provides a high-precision and reliable method for efficiently solving the contact response of a solid under any load distribution. Compared to the commonly used uniform load element discrete method, the analytical solution of the triangular element is more complex. However, it better simulates the characteristics of contact load distribution, accounting for situations where the load at the contact edge increases from zero or decreases to zero. The stress component under the action of the triangular and uniform load elements is derived based on the “excitation-response” characteristics of the contact influence coefficient matrix. This information is used to optimize the solution method of the triangular load discrete element. By constructing the stress solution in the form of a discrete convolution, including the influence coefficient matrix, the stress superposition effect of a target node under the action of all elements can be further simplified and accelerated by using the DC-FFT algorithm for highly repetitive matrix calculations. Programming and calculation analysis show that the proposed algorithm based on the triangular load element is accurate and efficient.
Keywords:triangular element  contact stress  DC-FFT  numerical solution  stress field
点击此处可从《重庆大学学报(自然科学版)》浏览原始摘要信息
点击此处可从《重庆大学学报(自然科学版)》下载免费的PDF全文
设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号