数学建模例题_之_饮酒驾驶模型[1]

合集下载
  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。

饮酒驾驶模型

摘要

本文针对酒后驾车造成交通事故死亡率高,以及根据国家质量检验检疫局发布的饮酒后驾车新标准,建立了饮酒后血液中酒精含量的数学模型。通过了解酒精在体内吸收,分布和排除的动态过程,及这些过程与人体内酒精反应的定量关系建立微分方程,运用药物动力学原理建立单室和双室模型。得出血液中的酒精含量)(t C ,与进入体内总酒量)(t x 、时间t 的函数关系式:

单室模型:()()()()k k e e x k t x t C a t k kt a a --==--0

双室模型:()()n n p n p t p

t p

t v t x v t x AUC AUC n n

∆⎪⎪⎭⎫

⎛++=--1001

本文还运用了 Wagner-Nelson 法(待吸收的百分数对时间作图法),与题中给出的参考数据在计算机运行的结果作对比。

本文还解决了如下问题:

1、从模型分析了大李第二次被判为饮酒驾车是因为二次饮酒,而使血液中酒精含量累积而超标。

2、对喝了低度酒多长时间驾车违反规则作了量化分析;

3、从单室模型得出了一个血液中酒精含量峰值计算公式:

()k k k gk t a a -=303.2max

4、用本文的模型对天天喝酒能否开车作了讨论。 本文最后对模型的优点和不足作了评价。

一、问题提出

据报载,2003年全国道路交通事故死亡人数为10.4372万,其中因饮酒驾车造成的占有相当的比例。

大李在中午12点喝了一瓶啤酒,下午6点检查时符合新的驾车标准,紧接着他在吃晚饭时又喝了一瓶啤酒,为了保险起见他呆到凌晨2点才驾车回家,又一次遭遇检查时却被定为饮酒驾车,这让他既懊恼又困惑,为什么喝同样多的酒,两次检查结果会不一样呢?

请你参考下面给出的数据(或自己收集资料)建立饮酒后血液中酒精含量的数学模型,并讨论以下问题:

1. 对大李碰到的情况做出解释;

2. 在喝了3瓶啤酒或者半斤低度白酒后多长时间内驾车就会违反上述标准,在以下情况下回答:

1)酒是在很短时间内喝的;

2)酒是在较长一段时间(比如2小时)内喝的。 3. 怎样估计血液中的酒精含量在什么时间最高。

4. 根据你的模型论证:如果天天喝酒,是否还能开车?

5.根据你做的模型并结合新的国家标准写一篇短文,给想喝一点酒的司机如何驾车提出忠告。

二、问题假设

1、机体分为中心室(I 室)和周边室(II 室),两个室的容积(即血液体积或药物分布容积)的过程中保持不变[1]。

2、药物从一室向另一室的转移速率,及向体外的排除速率,与该室的血药浓度成正比。

3、酒精含量的变化基本只受消除速度常数支配。

4、假定消除只发生在中心室,两个房室内酒精初始量都为零(即没有喝酒)。

5、酒在体内运动的配置和消除都是药物动力学过程。

6、人都是在精神状态正常情况下喝酒。

7、酒精可在整个机体内以同速度达到平衡。

三、符号定义

v :房室表观分布容积; k :酒精消除速度常数; a k :酒精吸收速度常数;

cp k :酒精转移速度常数(pc k );

)(t f :t 时刻体内吸收酒精的速度; m C :血酒浓度的最高峰值;

)(t C : 时刻的t 血液中酒精含量; )(t x :进入体内的总酒量; 0x :一次喝下的酒量;

()t x a :t 时刻体内吸收的酒精量; ()t x c :t 时刻中心室内的酒精量; ()t x p :t 时刻周边室内的酒精量;

n t :第n 次喝酒的时刻;

m t :血液浓度达到最高峰值的时刻;

I :已经代谢排泄酒物总量; p AUC :一次喝酒后的吸收总量;

四、模型建立

(一)、单室模型

将人的机理作为一个房室处理的模型,人喝酒后,酒精需要一定的吸收过程,可建立模型图(1):

图(1) 依条件及示意图,得到单室模型;

()kx t f dt

dx

-= (1) ()x x =0

()()v t x t C = (2) ()v x C C 000== 酒精逐渐进入血液循环后;

()()t k a a a a e k t x k t f x -==0* (3) 得到:

()()()k k e e x k t x a t k kt a a --=--0 (4) 将(2)式代入(3)得

()()()()k k e e x k t x t C a t k kt a a --==--0 (5)

根据动力学原理的有关计算方法,总结出的血液中酒精含量最大峰值和达到最大峰值时间计算公式[2]

x e C m kt 0max -= (6) ()k k k gk t a a -=303.2max (7)

(二)、双室模型[3]

二室模型假设酒精进入体内后在两个房室内配置,一个中心室,另一个是外周室,酒精在体内的配置和消除都是一级动力学过程,但酒精的吸收可以是任意的,见图(2):

图(2)

按照质量平衡原理,时间t -0范围内吸收进入体内的总酒量)(t x 为

()()()()I t X t X dt t f t x p c t

++==⎰0 (8)

其中

()dt t X k I t

c ⎰=0 (9)

代入式(5) 并在等号两边同时除以表观分布容积V 得到

()()()()v t x dt t c k t c v t x p t ++=⎰0 (10)

其中血液中酒精含量()t c =()

v

t x c 。

根据式(7), 当∞→t 时,计算酒精吸收分数的公式为:

()()()()⎰

⎰∞

+

+=0

0)(dt

t c k v t x dt t c k t c x t x p t

(11) 我们运用Wagner-Nelson 方法求解,对此,我们在算法作如下基本假设:在时间

1-n t 和n t 之间外周室酒精量()t x p 可以用线性插值近似逼近,因此()

t v

t x p -曲线下的面积

n n t t p

AUC 1

-可用梯形法进行运算

()()()()2111

1

---⎪⎪⎭⎫ ⎝⎛+==⎰

--n n n p n p t t p t t p

t t v t x v t x dt v t x AUC n

n n n (12) 则

()()()()()21

10

000

111---⎪⎪⎭

⎫ ⎝⎛++=+==---⎰⎰⎰n n n p n p t p

t t p t p t p t p

t t v t x v t x AUC dt v

t x dt v t x dt v t x AUC n n n n n n (13)

为了叙述方便,令

2

1--=∆n n n t

t t

则有

()()n n p n p t p

t p

t v t x v t x AUC AUC n n ∆⎪⎪⎭

⎫ ⎝⎛++=--10

1 (14)

相关文档
最新文档