高温高湿环境人体热负荷预警仿真模型

裴裴裴之Miraitowa

已于 2022-12-06 15:10:36 修改

阅读量346

点赞数 1

文章标签： matlab

于 2022-12-04 22:14:03 首次发布

本文链接：https://blog.csdn.net/m0_50852718/article/details/128178099

版权

热预警+Tanabe 65节点模型

背景/依据

Tanabe 65节点模型，将人体分为头部、胸部、骨盆、左右肩、左右手臂、左右手、左右大腿、左右小腿、和左右足共16个节段。根据组织特性将各段分为四个同心层：核心、肌肉、脂肪和皮肤。此外，还有一中央血池。
论文详情可以查询论文所在网址

公式以及简图

$\begin{aligned} & \frac{d T(i, 1)}{d t}= \frac{Q(i, 1)-B(i, 1)-D(i, 1)-R E S(i, 1)}{C(i, 1)} \\ & \frac{d T(i, 2)}{d t}= \frac{Q(i, 2)-B(i, 2)+D(i, 1)-D(i, 2)}{C(i, 2)} \\ & \frac{d T(i, 3)}{d t}= \frac{Q(i, 3)-B(i, 3)+D(i, 2)-D(i, 3)}{C(i, 3)} \\ & \frac{d T(i, 4)}{d t}=\frac{Q(i, 4)-B(i, 4)+D(i, 3)-Q_t(i, 4)-E(i, 4)}{C(i, 4)} \\ & \frac{d T(65)}{d t}=\frac{\sum_{i=1}^{16} \sum_{j=1}^4 B(i, j)}{C(65)} \end{aligned}$
模型展示
此外，论文中缺失一些参数，可根据MSP论文进行补充和更正，如：
$\begin{gathered} Q b b=0.778\left(\text { met }=58.2 W / \mathrm{m}^2\right) \\ W(i, j)=58.2(\operatorname{met}-Q b b) * A_{D u}(i) * \operatorname{Met} f(i) \end{gathered}$
因为 met 是全身代谢率, 根据 $M S P$ 模型, met $R_{\text {rest }}+0.0858 H R,(m e t)$ 又因为论文 word 中, $R_{\text {rest }}=92$ 次 $\min , H R$ 在 $40^{\circ} \mathrm{C}$ 达到 120 次 $\mathrm{min}$ , 故, met $= 3.6598$ .

展示

在这里插入图片描述

部分代码

classdef HeatShow < matlab.apps.AppBase

    % Properties that correspond to app components
    properties (Access = public)
        UIFigure                      matlab.ui.Figure
        dangerminutesEditField        matlab.ui.control.NumericEditField
        dangerminutesEditFieldLabel   matlab.ui.control.Label
        DropDown                      matlab.ui.control.DropDown
        Label_2                       matlab.ui.control.Label
        warninglampLamp               matlab.ui.control.Lamp
        warninglampLampLabel          matlab.ui.control.Label
        workingminutesEditField       matlab.ui.control.NumericEditField
        workingminutesEditFieldLabel  matlab.ui.control.Label
        AmbienttemperatureEditField   matlab.ui.control.NumericEditField
        AmbienttemperatureEditField_2Label  matlab.ui.control.Label
        startButton                   matlab.ui.control.Button
        Label                         matlab.ui.control.Label
        UIAxes                        matlab.ui.control.UIAxes
    end

    
    properties (Access = private)
        T_e % 环境温度
        dangerMin % 危险时间
        T % 运算温度
    end
    
  

    % Callbacks that handle component events
    methods (Access = private)

        % Button pushed function: startButton
        function startButtonPushed(app, event)
            app.T_e=app.AmbienttemperatureEditField.Value;