为研究大豆在机械通风过程温度变化，以华北地区直径为30 m，装粮高度为20 m的钢筋混凝土浅圆仓为研究对象，基于多孔介质传热传质理论，利用数值模拟软件COMSOL建立大豆储藏过程中粮堆内部温度变化的数值模型，并利用此模型预测了不同初始粮温和装粮高度条件下仓内粮堆温度变化过程。结果表明：粮层温度随着与风道距离的增加，变化幅度由剧烈逐渐变的缓慢。粮堆热阻的存在，导致热量在整个粮堆的传递过程表现出一定的迟滞性，从而使整个粮堆温度分布不均；仓内粮温随通风时间的变化而变化，在机械通风86 h后仓内粮堆温度分别降至4.9 ℃、6.2 ℃、7.3 ℃、8.4 ℃和9.6 ℃，此时与初始粮温相比已分别下降了59%、56%、54%、53%和52%，在相同边界条件下，五种初始粮温中，初始温度较低的粮堆，粮温降幅较大；机械通风过程，不同装粮高度（18 m、19 m和20 m）的粮堆温度变化差异呈现先增大后减小趋势，在机械通风19 h和110 h后，装粮高度相差1 m时，粮堆之间温差为0.1 ℃和0.6 ℃；在机械通风结束后，仓内粮堆温度分别为0.3 ℃、0.5 ℃和0.9 ℃，此时装粮高度相差1 m时，粮堆之间的温差为0.2～0.4 ℃，不同初始边界条件下粮堆温度变化趋势相同。综上，所建立的模型可对对多种初始边界条件的粮堆温度场进行预测。