✨ 长期致力于滑动摩擦界面、瞬态温度场构建、红外测温、发射率修正研究工作擅长数据搜集与处理、建模仿真、程序编写、仿真设计。✅ 专业定制毕设、代码✅如需沟通交流点击《获取方式》1基于红外辐射亮温的发射率动态修正方法提出了EDC-FFT算法利用摩擦过程中接触界面辐射亮温的时间序列通过快速傅里叶变换提取周期性分量再结合已知的摩擦副材料发射率-温度关系曲线多项式拟合阶次为3动态反演出真实发射率。在端面滑动摩擦试验中下试样材料为GCr15钢上试样为铜基粉末冶金载荷50N转速200r/min。修正前的辐射温度偏差达±25℃修正后的温度偏差缩小至±4℃。该方法还针对不同表面粗糙度Ra 0.4-1.6μm建立了发射率补偿表最大补偿量达0.12。2接触界面瞬态温度场有限体积法重构建立了三维瞬态热传导模型采用有限体积法离散时间步长0.01s空间网格尺寸0.5mm×0.5mm×0.2mm。边界条件接触界面热流密度根据摩擦功率动态分配分配系数采用Archard磨损模型计算。以红外热像仪测量的下试样侧表面温度作为校正条件通过共轭梯度法反演接触界面的热流分布。在300秒摩擦过程中重构的接触界面最高温度达380℃而侧表面最高温度仅120℃。与埋入式热电偶实测值370℃对比相对误差2.7%。该重构结果揭示了热点区域集中在半径中段的环带宽度约3mm。3滑动接触界面温度与工况参数的线性回归方程及测温系统构建基于80组不同载荷30-200N、转速100-500r/min和摩擦系数0.1-0.45的试验数据建立了多元线性回归方程T_interface 25 1.2*P 0.45*n 350*μ其中P为载荷(N)n为转速(r/min)μ为摩擦系数模型R20.91。利用该方程可以仅通过测量下试样侧表面温度T_side来估算界面温度T_interface 2.8*T_side - 40。基于此设计了简易测温系统采用单个红外探头对准下试样侧表面采样率100Hz信号经放大滤波后由STM32读取。在验证试验中估算的界面温度与反演值的最大偏差为22℃满足工程快速评估需求。import numpy as np from scipy.fft import fft, ifft from scipy.optimize import curve_fit def edc_fft_correction(radiance_time_series, emissivity_poly): # radiance_time_series: 辐射亮温时序, emissivity_poly: 发射率多项式系数 N len(radiance_time_series) fft_vals fft(radiance_time_series) # 提取主要频率分量 magnitude np.abs(fft_vals) dominant_freq np.argmax(magnitude[1:N//2]) 1 fft_filtered fft_vals.copy() fft_filtered[np.abs(fft_filtered) np.max(magnitude)*0.2] 0 filtered_signal np.real(ifft(fft_filtered)) # 假设真实温度与滤波后信号的关系 T_real_guess filtered_signal 10 # 偏移 def emissivity_func(T): return np.polyval(emissivity_poly, T) # 求解发射率 from scipy.optimize import fsolve T_corrected np.zeros_like(radiance_time_series) for i, T_rad in enumerate(radiance_time_series): def eq(T): return T_rad - T * emissivity_func(T) T_corrected[i] fsolve(eq, T_real_guess[i])[0] return T_corrected def fvm_temperature_field(q_interface, material_props): # 有限体积法简化示意一维瞬态导热 dx 0.5e-3 # 网格尺寸 dt 0.01 alpha material_props[alpha] # 热扩散系数 n_nodes 100 T np.ones(n_nodes) * 25 for step in range(1000): T_new T.copy() for i in range(1, n_nodes-1): T_new[i] T[i] alpha*dt/dx**2 * (T[i1]-2*T[i]T[i-1]) T_new[0] T[0] 2*alpha*dt/dx**2 * (q_interface[step] - (T[0]-T[1])/dx) # 边界热流 T T_new return T def regression_equation(P, n, mu): return 25 1.2*P 0.45*n 350*mu )
滑动摩擦副温度场模型应用优化【附仿真】
✨ 长期致力于滑动摩擦界面、瞬态温度场构建、红外测温、发射率修正研究工作擅长数据搜集与处理、建模仿真、程序编写、仿真设计。✅ 专业定制毕设、代码✅如需沟通交流点击《获取方式》1基于红外辐射亮温的发射率动态修正方法提出了EDC-FFT算法利用摩擦过程中接触界面辐射亮温的时间序列通过快速傅里叶变换提取周期性分量再结合已知的摩擦副材料发射率-温度关系曲线多项式拟合阶次为3动态反演出真实发射率。在端面滑动摩擦试验中下试样材料为GCr15钢上试样为铜基粉末冶金载荷50N转速200r/min。修正前的辐射温度偏差达±25℃修正后的温度偏差缩小至±4℃。该方法还针对不同表面粗糙度Ra 0.4-1.6μm建立了发射率补偿表最大补偿量达0.12。2接触界面瞬态温度场有限体积法重构建立了三维瞬态热传导模型采用有限体积法离散时间步长0.01s空间网格尺寸0.5mm×0.5mm×0.2mm。边界条件接触界面热流密度根据摩擦功率动态分配分配系数采用Archard磨损模型计算。以红外热像仪测量的下试样侧表面温度作为校正条件通过共轭梯度法反演接触界面的热流分布。在300秒摩擦过程中重构的接触界面最高温度达380℃而侧表面最高温度仅120℃。与埋入式热电偶实测值370℃对比相对误差2.7%。该重构结果揭示了热点区域集中在半径中段的环带宽度约3mm。3滑动接触界面温度与工况参数的线性回归方程及测温系统构建基于80组不同载荷30-200N、转速100-500r/min和摩擦系数0.1-0.45的试验数据建立了多元线性回归方程T_interface 25 1.2*P 0.45*n 350*μ其中P为载荷(N)n为转速(r/min)μ为摩擦系数模型R20.91。利用该方程可以仅通过测量下试样侧表面温度T_side来估算界面温度T_interface 2.8*T_side - 40。基于此设计了简易测温系统采用单个红外探头对准下试样侧表面采样率100Hz信号经放大滤波后由STM32读取。在验证试验中估算的界面温度与反演值的最大偏差为22℃满足工程快速评估需求。import numpy as np from scipy.fft import fft, ifft from scipy.optimize import curve_fit def edc_fft_correction(radiance_time_series, emissivity_poly): # radiance_time_series: 辐射亮温时序, emissivity_poly: 发射率多项式系数 N len(radiance_time_series) fft_vals fft(radiance_time_series) # 提取主要频率分量 magnitude np.abs(fft_vals) dominant_freq np.argmax(magnitude[1:N//2]) 1 fft_filtered fft_vals.copy() fft_filtered[np.abs(fft_filtered) np.max(magnitude)*0.2] 0 filtered_signal np.real(ifft(fft_filtered)) # 假设真实温度与滤波后信号的关系 T_real_guess filtered_signal 10 # 偏移 def emissivity_func(T): return np.polyval(emissivity_poly, T) # 求解发射率 from scipy.optimize import fsolve T_corrected np.zeros_like(radiance_time_series) for i, T_rad in enumerate(radiance_time_series): def eq(T): return T_rad - T * emissivity_func(T) T_corrected[i] fsolve(eq, T_real_guess[i])[0] return T_corrected def fvm_temperature_field(q_interface, material_props): # 有限体积法简化示意一维瞬态导热 dx 0.5e-3 # 网格尺寸 dt 0.01 alpha material_props[alpha] # 热扩散系数 n_nodes 100 T np.ones(n_nodes) * 25 for step in range(1000): T_new T.copy() for i in range(1, n_nodes-1): T_new[i] T[i] alpha*dt/dx**2 * (T[i1]-2*T[i]T[i-1]) T_new[0] T[0] 2*alpha*dt/dx**2 * (q_interface[step] - (T[0]-T[1])/dx) # 边界热流 T T_new return T def regression_equation(P, n, mu): return 25 1.2*P 0.45*n 350*mu )
