汽车内饰材料中挥发性有机物的非等温散发仿真计算方法与流程

1.本发明属于有机物分析领域，尤其是涉及一种汽车内饰材料中挥发性有机物的非等温散发仿真计算方法。


背景技术：

2.座椅、仪表板、地毯、顶棚、密封条等汽车内饰，一般由塑料、皮革、胶粘剂等石化产品制成，内含一定量的苯、甲苯、乙苯、二甲苯、苯乙烯、甲醛、乙醛、丙烯醛等挥发性有机物(voc)成分，将在消费者用车期间持续散发。通过降低装车前内饰所含voc的浓度，可有效改善车内空气质量，保障驾乘人员身体健康。
3.汽车内饰材料内含的voc会在浓度梯度的驱动下向周围空气中散发，该过程受环境温度的影响较大且比较复杂。一方面，环境温度越高，voc在材料中和空气中的扩散系数越大，在材料—空气界面处的分配系数越小，越容易从材料内部扩散到空气中，从而使材料内含的voc浓度越快地减小；另一方面，环境温度越高，材料voc的分子动能越大，从而使材料内含的voc的可散发浓度不降反增。针对上述非等温条件下voc的复杂散发情况，目前尚没有符合气体散发机理的仿真计算方法，严重影响了内饰材料中voc浓度目标设定、内饰材料仓储温度控制方案制定等工作的开展。


技术实现要素：

4.有鉴于此，本发明旨在提出一种汽车内饰材料中挥发性有机物的非等温散发仿真计算方法，以解决内饰材料中voc浓度目标较难设定、内饰材料仓储温度控制方案较难制定的问题。
5.为达到上述目的，本发明的技术方案是这样实现的：
6.汽车内饰材料中挥发性有机物的非等温散发仿真计算方法，具体方法如下：
7.s1、设定材料voc散发过程中的逐时温度c(t)，确定仿真温度范围，进行温度段的划分，计算各温度段的材料voc初始可散发浓度，并设定各温度段材料voc初始可散发浓度潜能；
8.s2、计算仿真温度范围内，各逐时温度下的材料voc扩散系数dm(c(t))、各逐时温度下的材料voc分配系数k(c(t))；
9.s3、计算仿真温度范围，各逐时温度下的空气中voc扩散系数da(c(t))；
10.s4、计算材料voc仿真散发初始时刻温度c(t0)下，材料voc初始可散发浓度cm(c(t0))；
11.s5、通过有限元分析构建材料voc散发几何仿真模型，设置材料网格、材料—空气界面处网格、空气网格，根据cm(c(t0))、dm(c(t0))、k(c(t0))、da(c(t0))，并结合菲克第二扩散定律及亨利定律，对每种类型的网格进行网格划分并设定网格节点参数，设定仿真散发的总时长t_total和仿真计算的时间步进t_step；
12.s6、仿真计算得到散发了一个时间步进后各网格节点上材料voc可散发浓度cm及
空气voc散发浓度ca；
13.s7、根据散发了一个时间步进后的温度c(t-t_step)及小于等于该温度c(t-t_step)的各温度段的初始可散发浓度潜能，将各网格节点上的材料voc可散发浓度cm按各温度段初始可散发浓度潜能cm0p_temp的相对比例分配，并赋值给温度小于等于c(t-t_step)的各温度段初始可散发浓度cm0v_temp；
14.s8、根据步骤s7获得散发了一个时间步进后的各温度段初始可散发浓度cm0v_temp，计算下一个时步开始时的温度c(t)，计算小于等于该温度c(t)的各温度段的初始可散发浓度的加和，并把加和的结果更新给各网格节点上的材料voc可散发浓度cm；
15.s9、判断当前的仿真时长t是否小于仿真总时长t_total，若是，则更新材料voc散发几何仿真模型网格节点dm(c(t))、k(c(t))、da(c(t))后，重复步骤s6-s8；若否，则仿真计算结束，输出各网格节点上的各温度段材料voc逐时可散发浓度cm0v_temp(t)、材料voc逐时可散发浓度cm(t)及空气voc逐时散发浓度ca(t)，计算材料网格中各网格节点逐时可散发浓度平均值及空气网格中各网格节点逐时散发浓度平均值。
16.进一步的，步骤s1中，确定仿真温度范围，进行温度段的划分，计算各温度段的材料voc初始可散发浓度，并设定各温度段材料voc初始可散发浓度潜能，具体方法如下：
17.s11、开展材料voc散发测试，得到25℃和65℃恒温条件下的材料voc散发参数，包括25℃条件下的初始可散发浓度cm0_25mg/m3、65℃条件下的初始可散发浓度cm0_65mg/m3、25℃条件下的扩散系数dm_25m2/s、65℃条件下的扩散系数dm_65m2/s、25℃条件下的分配系数k_25、65℃条件下的分配系数k_65；
18.s12、执行温度效应方程参数拟合，根据25℃和65℃恒温条件下的材料voc散发参数得到用于计算任意温度下的材料voc初始可散发浓度的cm0_a和cm0_b、材料内voc扩散系数的dm_a和dm_b、材料—空气界面处的voc分配系数的k_a和k_b；
19.其中：
20.cm0_b＝-2518.1
×
(log(cm0_65/cm0_25)+0.5
×
0.1259)
21.cm0_a＝cm0_25/(298.15^(-0.5)
×
exp(cm0_b/298.15))
22.dm_b＝-2518.1
×
(log(dm_65/dm_25)-1.25
×
0.1259)
23.dm_a＝dm_25/(298.15^(1.25)
×
exp(dm_b/298.15))
24.k_b＝-2518.1
×
(log(k_65/k_25)-0.5
×
0.1259)
25.k_a＝k_25/(298.15^(0.5)
×
exp(k_b/298.15))；
26.s13、设定仿真温度范围为(25～65)℃，进行温度段的划分，具体方法为划分(25～65)℃范围内每1℃划分一个温度段，temp＝25,26,
…
,65；
27.s14、根据材料voc初始可散发浓度参数，计算各温度段的材料voc初始可散发浓度，并设定各温度段材料voc初始可散发浓度潜能。
28.进一步的，步骤s14中，根据材料voc初始可散发浓度参数，计算各温度段的材料voc初始可散发浓度，具体计算方法：
29.根据cm0_a和cm0_b，计算得到(25～65)℃范围内各温度段的初始可散发浓度cm0v_25mg/m3，cm0v_26mg/m3，cm0v_27mg/m3，
……
，cm0v_65mg/m3；cm0v_25＝cm0_a
×
(25+273.15)^-0.5
×
exp(cm0_b/(25+273.15))；cm0v_temp’＝cm0_a
×
(temp+273.15)^-0.5
×
exp(cm0_b/(temp+273.15))-cm0_a
×
(temp-1+273.15)^-0.5
×
exp(cm0_b/(temp-1+
273.15))，temp’＝26,27,
…
,65；
30.设定各温度段材料voc初始可散发浓度潜能，具体方法：
31.设定(25～65)℃范围内各温度段材料voc初始可散发浓度潜能cm0p_25mg/m3，cm0p_26mg/m3，cm0p_27mg/m3，
……
，cm0p_65mg/m3，其中：
32.cm0p_temp＝cm0v_temp，temp＝25,26,
…
,65。
33.进一步的，步骤s2中，计算仿真温度范围内，各逐时温度下的材料voc扩散系数dm(c(t))及各逐时温度下的材料voc分配系数k(c(t))，具体计算公式：
34.dm(c(t))＝dm_a
×
(c(t)+273.15)^1.25
×
exp(dm_b/(c(t)+273.15))
35.k(c(t))＝k_a
×
(c(t)+273.15)^0.5
×
exp(k_b/(c(t)+273.15))；
36.步骤s3中，计算各逐时温度下的空气中voc扩散系数da(c(t))，具体方法为，查阅文献获得材料voc在25℃空气中的扩散系数da_25，
37.da(c(t))＝da_25
×
((c(t)+273.15)/298.15)^1.81。
38.进一步的，步骤s4中，计算材料voc仿真散发初始时刻温度c(t0)下，材料voc初始可散发浓度cm(c(t0))，具体计算公式：
39.cm(c(t0))＝cm0_a
×
(c(t0)+273.15)^-0.5
×
exp(cm0_b/(c(t0)+273.15))。
40.进一步的，步骤s5中，通过有限元分析构建材料voc散发几何仿真模型，设置材料网格、材料—空气界面处网格、空气网格，根据cm(c(t0))、dm(c(t0))、k(c(t0))、da(c(t0))，并结合菲克第二扩散定律及亨利定律，对每种类型的网格进行网格划分并设定网格节点参数，设定仿真散发的总时长t_total s和仿真计算的时间步进t_step s，具体方法：
41.s51、在有限元分析软件中构建材料voc散发几何模型，设置材料网格、材料—空气界面处网格、空气网格，对每种类型的网格进行网格划分；
42.s52、在各网格节点上设定相关参数，包括：在材料网格节点上设定仿真散发初始t0时刻的初始可散发浓度cm(c(t0))和扩散系数dm(c(t0))；在材料—空气界面处网格节点上设定分配系数k(c(t0))；在空气网格节点上设定初始浓度ca＝0和扩散系数da(c(t0))；
43.s53、设定仿真散发的总时长t_total s和仿真计算的时间步进t_step s；
44.s54、在材料和空气网格节点上应用菲克第二扩散定律，在材料—空气界面处网格节点上应用亨利定律。
45.进一步的，步骤s7中，根据散发了一个时间步进后的温度c(t-t_step)及小于等于该温度c(t-t_step)的各温度段的初始可散发浓度潜能，将各网格节点上的材料voc可散发浓度cm按各温度段初始可散发浓度潜能cm0p_temp的相对比例分配，并赋值给温度小于等于c(t-t_step)的各温度段初始可散发浓度cm0v_temp，具体公式：
[0046][0047]
cm0v_temp’＝if(c(t-t_step)《＝temp
’‑
1,cm0v_temp’，
[0048][0049]
temp’＝26,27,
…
,65，temp＝temp’:65；
[0050]
temp＝temp’:65的含义是temp从temp’至65以1为步长依次取值；
[0051]
步骤s8中，根据步骤s7获得散发了一个时间步进后的各温度段初始可散发浓度cm0v_temp，计算下一个时步开始时的温度c(t)，计算小于等于该温度c(t)的各温度段的初
始可散发浓度的加和，并把加和的结果更新给各网格节点上的材料voc可散发浓度cm，具体公式为，
[0052][0053]
进一步的，步骤s9中，计算材料网格中各网格节点逐时浓度平均值及空气网格中各网格节点逐时散发浓度平均值，具体方法如下：
[0054]
s91、根据各温度段材料voc逐时可散发浓度cm0v_temp(t)、材料voc逐时可散发浓度cm(t)，计算材料网格中各网格节点逐时浓度平均值，获得材料voc在逐时变化的温度条件c(t)下的逐时可散发浓度和仿真温度范围内各温度段的逐时可散发浓度，通过调节温度条件，获得材料voc可散发浓度的管控目标值；
[0055]
s92、根据空气voc逐时散发浓度ca(t)，计算空气网格中各网格节点逐时散发浓度平均值，获得材料voc在变化的温度条件下经一定时间散发到周围空气中的空气voc逐时浓度值。
[0056]
进一步的，步骤s91中，计算材料网格中各网格节点逐时浓度平均值，获得材料voc在逐时变化的温度条件c(t)下的逐时可散发浓度平均值和仿真温度范围内各温度段的逐时可散发浓度平均值，具体方法：
[0057]
记材料网格节点数为x，各网格节点命名分别为n1、n2、
……
、nx；
[0058]
计算材料各网格节点上的逐时浓度的平均值，设定cm_avg(t)为材料voc在逐时变化的温度条件c(t)下的逐时可散发浓度平均值，具体公式：
[0059][0060]
设定cm0v_25_avg(t)、cm0v_26_avg(t)、
……
、cm0v_65_avg(t)依次表示(25～65)℃范围内各温度段的逐时可散发浓度平均值，具体计算公式为，
[0061][0062]
temp＝25,26,
…
,65，t＝0:t_step:t_total，
[0063]
t＝0:t_step:t_total含义是时间变量t从0至t_total以t_step为步长依次取值。
[0064]
进一步的，步骤s92中，计算空气网格中各网格节点逐时散发浓度平均值，获得材料voc在变化的温度条件下经一定时间散发到周围空气中的空气voc逐时浓度值，具体方法：
[0065]
记空气网格节点数为y，各网格节点分别为a1、a2、
……
、ay；
[0066]
设定空气各网格节点上的逐时浓度的平均值ca_avg(t)，具体计算公式为，
[0067][0068]
ca_j(t)是t时刻在空气网格的j网格节点上读取的空气voc浓度值，j表示为网格节点a1、a2、
……
、ay中的一个网格节点。
[0069]
相对于现有技术，本发明所述的汽车内饰材料中挥发性有机物的非等温散发仿真计算方法，具有以下有益效果：
[0070]
(1)本发明所述的汽车内饰材料中挥发性有机物的非等温散发仿真计算方法，采用有限元分析仿真计算方法符合voc气体在固体材料中的散发机理，因此仿真结果的可靠
性较高，进而使基于仿真结果设定的内饰材料中voc浓度的管控目标值及其配套的仓储温度控制方案的可靠性较高。
[0071]
(2)本发明所述的汽车内饰材料中挥发性有机物的非等温散发仿真计算方法，由于固体材料中的voc一般认为有两种存在状态，一种状态是存在于固体内部微孔中的空气相voc所处于的可散发状态；另一种状态是吸附结合在固体上的吸附相voc所处于的不可散发状态；不可散发状态可以向可散发状态转化，但需要克服一个吸附能量势垒，通过升高温度增加吸附相voc的分子动能，即可克服这个能量势垒，本技术计算各温度段材料voc初始可散发浓度，将各温度段的材料voc初始可散发浓度作为voc容器，体现了这两种状态在非等温散发条件下的互相转化；
[0072]
(3)本发明所述的汽车内饰材料中挥发性有机物的非等温散发仿真计算方法，由于不可散发状态的voc与可散发状态的voc相互转化所涉及到的吸附能量势垒是与固体材料属性及voc属性有关的特征常数，不随voc浓度的变化而改变，该常数型能量势垒会直接影响可散发状态的voc向各温度段的voc容器中分配时的分配比例，本技术中根据仿真时步开始时的温度c(t-t_step)和小于等于该温度c(t-t_step)的各温度段的初始可散发浓度潜能，将材料各网格节点上的浓度cm按各初始可散发浓度潜能的相对比例分配，并对应赋值给温度小于等于c(t-t_step)的各温度段初始可散发浓度的绝对大小，通过将各温度段的初始可散发浓度潜能作为分配比例的确定因子，提高仿真结果的准确性。
附图说明
[0073]
构成本发明的一部分的附图用来提供对本发明的进一步理解，本发明的示意性实施例及其说明用于解释本发明，并不构成对本发明的不当限定。在附图中：
[0074]
图1为本发明实施例所述的汽车内饰材料中挥发性有机物的非等温散发仿真计算方法示意图。
具体实施方式
[0075]
需要说明的是，在不冲突的情况下，本发明中的实施例及实施例中的特征可以相互组合。
[0076]
在本发明的描述中，需要理解的是，术语“中心”、“纵向”、“横向”、“上”、“下”、“前”、“后”、“左”、“右”、“竖直”、“水平”、“顶”、“底”、“内”、“外”等指示的方位或位置关系为基于附图所示的方位或位置关系，仅是为了便于描述本发明和简化描述，而不是指示或暗示所指的装置或元件必须具有特定的方位、以特定的方位构造和操作，因此不能理解为对本发明的限制。此外，术语“第一”、“第二”等仅用于描述目的，而不能理解为指示或暗示相对重要性或者隐含指明所指示的技术特征的数量。由此，限定有“第一”、“第二”等的特征可以明示或者隐含地包括一个或者更多个该特征。在本发明的描述中，除非另有说明，“多个”的含义是两个或两个以上。
[0077]
在本发明的描述中，需要说明的是，除非另有明确的规定和限定，术语“安装”、“相连”、“连接”应做广义理解，例如，可以是固定连接，也可以是可拆卸连接，或一体地连接；可以是机械连接，也可以是电连接；可以是直接相连，也可以通过中间媒介间接相连，可以是两个元件内部的连通。对于本领域的普通技术人员而言，可以通过具体情况理解上述术语
在本发明中的具体含义。
[0078]
下面将参考附图并结合实施例来详细说明本发明。
[0079]
1)开展材料voc散发测试，得到25℃和65℃恒温条件下的材料voc关键散发参数，包括25℃条件下的初始可散发浓度cm0_25mg/m3、65℃条件下的初始可散发浓度cm0_65mg/m3、25℃条件下的扩散系数dm_25m2/s、65℃条件下的扩散系数dm_65m2/s、25℃条件下的分配系数k_25、65℃条件下的分配系数k_65。
[0080]
2)执行温度效应方程参数拟合，得到用于计算任意温度下的材料内voc初始可散发浓度的cm0_a和cm0_b，得到用于计算任意温度下的材料内voc扩散系数的dm_a和dm_b，得到用于计算任意温度下的材料—空气界面处的voc分配系数的k_a和k_b。其中：
[0081]
cm0_b＝-2518.1
×
(log(cm0_65/cm0_25)+0.5
×
0.1259)
ꢀꢀꢀꢀꢀꢀ
(式1)
[0082]
cm0_a＝cm0_25/(298.15^(-0.5)
×
exp(cm0_b/298.15))
ꢀꢀꢀꢀꢀꢀ
(式2)
[0083]
dm_b＝-2518.1
×
(log(dm_65/dm_25)-1.25
×
0.1259)
ꢀꢀꢀꢀꢀꢀ
(式3)
[0084]
dm_a＝dm_25/(298.15^(1.25)
×
exp(dm_b/298.15))
ꢀꢀꢀꢀꢀꢀ
(式4)
[0085]
k_b＝-2518.1
×
(log(k_65/k_25)-0.5
×
0.1259)
ꢀꢀꢀꢀꢀꢀ
(式5)
[0086]
k_a＝k_25/(298.15^(0.5)
×
exp(k_b/298.15))
ꢀꢀꢀꢀꢀꢀ
(式6)
[0087]
3)根据cm0_a和cm0_b，计算得到(25～65)℃范围内各温度段的初始可散发浓度cm0v_25mg/m3，cm0v_26mg/m3，cm0v_27mg/m3，
……
，cm0v_65mg/m3等。后缀符号“v”是变量variable的首字母，这些参数的值将随着材料voc的散发而变化更新，反映了不同温度段的材料voc初始可散发浓度的绝对大小。其中：
[0088]
cm0v_25＝cm0_a
×
(25+273.15)^-0.5
×
exp(cm0_b/(25+273.15))
ꢀꢀꢀꢀ
(式7)
[0089][0090]
4)设定(25～65)℃范围内各温度段材料voc初始可散发浓度潜能cm0p_25mg/m3，cm0p_26mg/m3，cm0p_27mg/m3，
……
，cm0p_65mg/m3等。后缀符号“p”是潜能potential的首字母，这些参数的值是常数，不随材料voc的散发而变化更新，反映了不同温度段的材料voc初始可散发浓度的相对大小。其中：
[0091][0092]
5)设定材料voc散发过程中的逐时温度c(t)℃。
[0093]
6)基于材料voc仿真散发初始t0时刻的温度c(t0)，计算仿真散发的初始可散发浓度cm(c(t0))mg/m3。
[0094]
cm(c(t0))＝cm0_a
×
(c(t0)+273.15)^-0.5
×
exp(cm0_b/(c(t0)+273.15))
ꢀꢀ
(式10)
[0095]
7)计算(25～65)℃范围内各逐时温度下的材料voc关键散发参数，包括扩散系数dm(c(t))m2/s和分配系数k(c(t))。其中：
[0096]
dm(c(t))＝dm_a
×
(c(t)+273.15)^1.25
×
exp(dm_b/(c(t)+273.15))
ꢀꢀꢀꢀ
(式11)
[0097]
k(c(t))＝k_a
×
(c(t)+273.15)^0.5
×
exp(k_b/(c(t)+273.15))
ꢀꢀꢀꢀ
(式12)
[0098]
8)文献检索得到voc在25℃空气中的扩散系数da_25m2/s，并计算各逐时温度下的空气中voc扩散系数da(c(t))m2/s。
[0099]
da(c(t))＝da_25
×
((c(t)+273.15)/298.15)^1.81
ꢀꢀꢀꢀ
(式13)
[0100]
9)在有限元分析软件中构建材料voc散发几何模型并进行网格划分，在各网格节点上设定相关参数，包括：1)在材料网格节点上设定仿真散发初始t0时刻的初始可散发浓度cm(c(t0))mg/m3和扩散系数dm(c(t0))m2/s；2)在材料—空气界面处网格节点上设定分配系数k(c(t0))；3)在空气网格节点上设定初始浓度ca＝0mg/m3和扩散系数da(c(t0))m2/s。
[0101]
10)设定仿真散发的总时长t_total s和仿真计算的时间步进t_step s。
[0102]
11)在材料和空气网格节点上应用菲克第二扩散定律，在材料—空气界面处网格节点上应用亨利定律，
[0103]
12)计算得到散发了一个时步后各网格节点上的浓度cm mg/m3及ca mg/m3。
[0104]
13)根据仿真散发了一个时步后的温度c(t-t_step)和小于等于该温度c(t-t_step)的各温度段材料voc初始可散发浓度潜能cm0p_25mg/m3、cm0p_26mg/m3等，将材料各网格节点上的浓度cm按各温度段材料voc初始可散发浓度潜能的相对比例分配，并赋值给温度小于等于温度c(t-t_step)的各温度段初始可散发浓度cm0v_25mg/m3、cm0v_26mg/m3等。
[0105]
其中：
[0106][0107][0108]
if条件语句之后的temp＝25:65的含义是temp从25至65以1为步长依次取值，直至if条件语句判断为真后停止后续取值。
[0109]
14)根据步骤s7获得散发了一个时间步进后的各温度段初始可散发浓度cm0v_temp，计算下一个时步开始时的温度c(t)，计算小于等于该温度c(t)的各温度段的初始可散发浓度的加和，并把加和的结果更新给各网格节点上的材料voc可散发浓度cm mg/m3，
[0110][0111]
15)判断当前的仿真时长t是否小于仿真总时长t_total。若是，则更新材料网格节点上的扩散系数为dm(c(t))m2/s，更新材料—空气边界处网格节点上的分配系数为k(c(t))，更新空气网格节点上的扩散系数为da(c(t))m2/s，然后再依次执行12～14；若否，则仿真计算结束，输出各网格节点上的逐时浓度cm(t)mg/m3、ca(t)mg/m3、cm0v_25(t)mg/m3、cm0v_26(t)mg/m3、
……
、cm0v_65(t)mg/m3。
[0112]
16)记材料网格节点数为x，各网格节点分别为n1、n2、
……
、nx，计算材料各网格节点上的逐时浓度的平均值cm_avg(t)mg/m3、cm0v_25_avg(t)mg/m3、cm0v_26_avg(t)mg/m3、
……
、cm0v_65_avg(t)mg/m3，即为材料voc在逐时变化的温度条件c(t)下的逐时可散发浓度和(25～65)℃范围内各温度段的逐时可散发浓度。通过改变逐时变化的温度条件c
(t)，即改变内饰材料的仓储温度控制方案，可以仿真计算得到材料的逐时可散发浓度和(25～65)℃范围内各温度段的逐时可散发浓度，基于这些温度—可散发浓度数据即可准确设定内饰材料中voc浓度的管控目标值及其配套的仓储温度控制方案。其中：
[0113][0114][0115]
式(17)的含义是时间变量t从0至t_total以t_step为步长依次取值，将同一时刻下的各网格节点上的可散发浓度cm_h(t)mg/m3相加后，再除以材料网格节点数x。
[0116]
17)记空气网格节点数为y，各网格节点分别为a1、a2、
……
、ay，计算空气各网格节点上的逐时浓度的平均值ca_avg(t)mg/m3，即为材料voc在变化的温度条件c(t)下经一定时间散发到周围空气中导致的空气voc逐时浓度。
[0117][0118]
以上所述仅为本发明的较佳实施例而已，并不用以限制本发明，凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。