【数值仿真】基于有限差分法的钻削过程物体热源分布计算

科研GGB

已于 2024-05-13 13:37:57 修改

阅读量196

点赞数

分类专栏：热传导/红外文章标签：热力学 matlab 仿真器

于 2023-05-15 15:39:07 首次发布

本文链接：https://blog.csdn.net/tkl32172/article/details/130685437

版权

热传导/红外专栏收录该内容

2 篇文章 2 订阅

订阅专栏

基于有限差分法的钻孔过程物体热源分布计算

钻削过程[1]

代码预览

clc;clear;close all;

% 需要计算的时间节点个数
f = 3000;
% 方形试块尺寸L(m) 空间时间步长 deltat deltax
L = [0.015 0.015 0.01]; deltat = 0.005 ;deltax = 0.0005;
% 钻头热源尺寸
LHeat = [0.00106/2 0.006/2 0.001426];
% 密度rou k为热导率 Cp为比热容
rou = 1800; k = [6.92,5.613,5.613];  Cp = 1000;
% 钻头给进速度
v = 0.020/60;
% 计算点云坐标MeshPoint
lx=L(1);ly=L(2);lz=L(3);
[xt,yt,zt] = meshgrid(0:1:lx/deltax,0:1:ly/deltax,0:1:lz/deltax);
xa = xt(:); za = zt(:);ya = yt(:);
% MeshPoint0 = [xa,ya,za];
MeshPoint = [xa,ya,za].*deltax - [L(1)/2 L(2)/2 L(3)] ;
% 是否显示网格
DisplayMesh = 0;
% 网格点连接关系MeshAttach
MeshAttach = HexahedralElement([xa,ya,za],DisplayMesh,k); % 点数多时，显示网格会卡崩

计算结果

在这里插入图片描述

下载联系

Q 191927818 打包打开直接跑20￥，简单问题会直接回复，讲解100￥/h，欢迎咨询

参考文献

[1] Liu J , Chen G , Ji C , et al. An investigation of workpiece temperature variation of helical milling for carbon fiber reinforced plastics (CFRP)[J]. International Journal of Machine Tools & Manufacture, 2014, 86:89-103.