1.一种浴缸温控模型的构建方法，其特征在于：包括如下步骤，步骤1、构建浴缸内水温度关于空间和时间上的方程模型为u′＝Δ3u+f

其中，u是一个三维时空变量函数u(x，y，z，t)，u′为u的导数，Δ3是三维拉普拉斯算符F(x，y，z，t)是源汇项，x，y，z，t分别为三维空间坐标与时间，设上式有分离变量解：u＝f(x，y，t)g(z，t)其中，g(z，t)是切向偏差系数；

步骤2、建立含热源的热传导方程模型；

步骤2.1建立泛定方程

首先建立含热源的热传导方程：

2

u′＝aΔu+f

其中，Δ表示拉普拉斯算子、Δ3表示三维拉普拉斯算子，a为

2 ‑7 ‑1 ‑1

热传导系数，且a≈2.1429×10 W·m ·K ，F为dt时间内，在(x，y，z)处的体积元内产生的热量；

步骤2.2、定解条件的确定

确定上述含热源的热传导方程在已知初始条件和边界条件下的定解问题；

初始条件：表示初始时刻浴缸内温度的分布情况，u(x，y，z，t)|t＝0＝U0其中U0为浴缸内的初始温度；

边界条件：

设浴缸内水的上表面与空气的接触面为Γ1，浴缸内水的其余面与浴缸壁的接触面为Γ2；空气介质e的三维时空间u的值为u1；则从浴缸中的水流到空气介质中的热量和两者的温度差成正比，即：dQ＝γ(u‑u1)dsdt

其中，γ为热交换系数；取正值，s为面积；ds为面积元，即单位面积；则流过接触面Γ1的热量为：即

其中 为温度变化率，表示面外法向方向，k为比例系数；

浴缸中的水的上表面与空气的边界条件为：浴缸中的水与浴缸壁界面的边界条件为：其中σ水‑气和σ水‑缸分别为水的上表面与空气热传导系数和水与浴缸壁界面的热传导系数；

步骤3、建立浴缸平均温度

定义

式中Ω是三维区域[0，X0]×[0，Y0]×[0，Z0]，其中X0，Y0，Z0为浴缸的长宽高，uxyz为浴缸中水的温度， 为没有人时浴缸中水的平均温度，即函数u取值为步骤4、构建浴缸在人体活动及加水情况下的温度模型步骤4.1、构建人体模型，设人的温差温度为恒定的Up，即人与空气的温差为恒定的Up，除去人体部分的浴缸中平均水温为：为含有人时的浴缸的平均温度， 为除去人体外浴缸中的水所在区域；

Vω表示浴缸中水所占的空间区域大小，且有：步骤4.2、建立人体运动时的浴缸温度模型为：表示由于运动，dt时间后浴缸中的水温是运动前水温的平均值；

步骤4.3、建立注入热水时浴缸温度模型；

设热水在落入浴缸中的时候，维持柱体形状不变，设热水流速为Q，设等效热柱半径为r，则在Δt时间内注入热水的量为：2

QΔt＝πrZ0

式中，Δt为时间维度的剖分间隔；

设热水温度为Uh，以每一层中水龙头正下方对应点(x0，y0)为圆心，以r为半径的区域，作为恒定热源区域，设某一层中水温分布为矩阵U(i，j，t)，设水的流速为v，则经过时间Δt后，该层中的某一点(i，j)运动至点(i′，j)，其中i与i′满足(i′‑i)Δx＝vΔt

其中i是对[0，X0]进行等距划分后的第i个格点，j是对[0，Y0]进行等距划分后的第j个格点，则注入热水的温度模型为：式中Δt与Δx分别为时间维度与x方向的剖分间隔；

步骤5、计算水温舒适度；

设人每次打开阀门的时间间隔Tk，开启阀门的次数为Nk，则在一次沐浴过程中，添加的热水质量Mω可以表示为：Mω＝QNkTk

即

2

Mω＝πrhΔtNkTk

设初始水温为加水的目标温度，得到衡量沐浴过程中水温的舒适度函数ζ(Nk，Tk)，式中ε为正常数，T0为沐浴总时间，ζ为ζ(Nk，Tk)；

为没有人时浴缸中水的平均温度随时间t的变化函数，则规划目标函数为：

2.根据权利要求1所述浴缸温控模型的构建方法，其特征在于：步骤2中泛定方程的构建过程为：设dt时间内，沿面积元ds的外法线方向流过的热量 与该面积元两侧的温度变化率 成正比，比例系数为k，则有其中，u为温度函数u(x，y，z，t)的简记；是面外法向方向； 为面积元ds的外法线方向的单位向量；

对于一个封闭的三维区域Ω，在dt时间内其内部的热量的变化dQ，通过对体积元的闭合面积分得到：进一步对时间积分，得到从t1到t2时刻体积元内部的热量变化Q1，由高斯公式对上式作如下变化：

式中：∑表示平面区域[0，X0]×[0，Y0]；

由于物体一段时间内变化的热量等于物体比热容c、物体质量m和物体这段时间内变化温度Δu的乘积，所以体积元内的热量变化表示为Q2Q2＝cmΔu＝∫∫∫Ωcρ[u(x，y，z，t2)‑u(x，y，z，t1)]dxdydz其中ρ为物体的密度，Ω为一个封闭的体积元；

对Q2进行等价变形得到：

由于Q1＝Q2，即：

得到：

展开表达式得：

由上述推导出：若物体内部有热源，dt时间内，在(x，y，z)处的体积元内产生的热量为F(x，y，z，t)，则得到含热源的热传导方程。

3.根据权利要求2所述浴缸温控模型的构建方法，其特征在于：步骤3中浴缸温度的三维模型的构建过程为：设不同深度的水的热传导性质只与其深度有关，与其平面位置无关，依此对于每个空间点的温度信息ui，将其分离为温度‑位置函数f与温度‑深度h的函数，如下：u(xi，yi，zi，t)＝f(xi，yi，t)h(zi，t)+λ(xi，yi，zi，t)，(i＝1，2，...)其中xi，yi，zi为某个空间点的空间位置，λ为扰动误差，h(z，t)为高度为z的一层的变温系数；

依据上式，对三维空间进行分层处理，对每一层水面看作二维热传导问题，即对f(x，y，t)进行求解；而对于纵向的热交换，采用一维热传导问题的解决方法h(z，t)进行求解；

则建立竖直方向一维热传导方程模型为：表示为浴缸在竖直方向的平均温度，uz为浴缸在竖直方向的温度，Z0为浴缸竖直方向尺寸；

平面二维热传导模型为：

表示浴缸在平面内的平均温度，uxy为浴缸在平面内的温度，式中∑表示平面区域[0，X0]×[0，Y0]，X0Y0为浴缸平面方向尺寸；

根据上述一维热传导方程模型、二维热传导模型得浴缸温度的三维模型。

4.根据权利要求1至3之一所述浴缸温控模型的构建方法，其特征在于：步骤4中人体模型整体呈“L”形，由两个长方体组成。