发明内容
本发明的目的是利用有限元分析软件-COMSOL Multiphysics为基础,对微流道加速度计动、静态特性进行分析,寻找传感器最优性能感测场合的微流道加速度计感应腔设计参数的优化方法。
本发明的步骤是:
(1)建立包含能斯特-普朗克方程和层流场的三维动态模型,该过程包括以下步骤:
(1-1)在COMSOL Multiphysics 4.3a建立一个包含能斯特-普朗克方程和层流场的三维动态模型;
(1-2)在模型中,定义3种变量,分别用 、 、表示;
(1-3)在模型中,将预置求解定义为瞬态模式;
(2)根据微流道感应腔设计尺寸,在步骤(1)所建模型内画出所设计的微流道加速度计微流道感应腔三维结构,该过程包括以下步骤:
(2-1)在模型中,将长度单元设置为μm,角度单元设置为Degrees;
(2-2)在模型中,以坐标原点(0,0)为起点,沿x轴正方向画直径为 μm,高度依次为 μm, μm,μm, μm, μm, μm, μm, μm, μm的圆柱体,分别命名为:a,b,c,d,e,f,g,h,i;
(3)输入所设计微流道加速度计电解液导电率、相对介电常数、电解液密度、粘滞度这四个参数值;
(4)在能斯特-普朗克方程中定义电极、电势及边界条件和约束方程:
电气绝缘:将圆柱体a, c, e, g, i定义为绝缘体,约束方程为:,无通量的约束方程为:,其中,
电场迁移的约束方程为:,设定圆柱体b, h为阳极,电势值为U1,圆柱体d,f为阴极,电势值为U2;
(5)层流场中定义入口、出口、流速及边界条件和约束方程:
入口方程为:,边界条件定义为速度,且定义为方向流入速度,为设一初始值;
出口方程为:,,其中,是温度,是压力,是动力粘度,是速度场,是圆柱体长度,边界条件定义为压力,无粘滞应力,为设一初始值;
边界条件定义为无滑移;
(6)对求解域进行网格剖分,该过程包括以下步骤:
(6-1)对所建的模型进行粗化网格剖分,并对x向、y向、z向的缩放比例进行设置,对穿过删除的控制实体平滑下的迭代数、最大单元深度进行设置;
整个模型建立完毕;
(7)计算:借助有限元仿真软件仿真求解;
(8)计算完毕后在结果中查看并分析计算结果;
(9)根据计算结果优化设计参数。
本发明所提出的方法,具有以下优点:
本发明方法可以在没有微流道加速度计样品情况下,通过计算机仿真得到某些设计参数,进而了解微流道加速度计反应腔内部情况,通过仿真分析设计参数的情况,加快研发速度,降低研发成本。
本发明方法基于有限元分析软件的建立,对多物理场进行耦合,可根据实际设计制定需仿真的结构,再根据设计要求制定需要达到的技术参数,最后给出优化的设计参数。
可在已知微流道加速度计各项设计参数情况下,通过计算机仿真得到某设计参数下的传感器运动情况,并对其动、静态特性进行分析,寻找传感器最优性能的感测场合,使得设计的微流道加速度计内部各种参数达到最优水平,大大提高了研发效率,节约研发成本。
具体实施方式
本发明的步骤是:
(1)建立包含能斯特-普朗克方程和层流场的三维动态模型,该过程包括以下步骤:
(1-1)在COMSOL Multiphysics 4.3a建立一个包含能斯特-普朗克方程和层流场的三维动态模型;
(1-2)在模型中,定义3种变量,分别用 、 、表示;
(1-3)在模型中,将预置求解定义为瞬态模式;
(2)根据微流道感应腔设计尺寸,在步骤(1)所建模型内画出所设计的微流道加速度计微流道感应腔三维结构,该过程包括以下步骤:
(2-1)在模型中,将长度单元设置为μm,角度单元设置为Degrees;
(2-2)在模型中,以坐标原点(0,0)为起点,沿x轴正方向画直径为 μm,高度依次为 μm, μm,μm, μm, μm, μm, μm, μm, μm的圆柱体,分别命名为:a,b,c,d,e,f,g,h,i;
(3)输入所设计微流道加速度计电解液导电率、相对介电常数、电解液密度、粘滞度这四个参数值;
(4)在能斯特-普朗克方程中定义电极、电势及边界条件和约束方程:
电气绝缘:将圆柱体a, c, e, g, i定义为绝缘体,约束方程为:,无通量的约束方程为:,其中,
电场迁移的约束方程为:,设定圆柱体b, h为阳极,电势值为U1,圆柱体d,f为阴极,电势值为U2;
(5)层流场中定义入口、出口、流速及边界条件和约束方程:
入口方程为:,边界条件定义为速度,且定义为方向流入速度,为设一初始值;
出口方程为:,,其中,是温度,是压力,是动力粘度,是速度场,是圆柱体长度,边界条件定义为压力,无粘滞应力,为设一初始值;
边界条件定义为无滑移;
(6)对求解域进行网格剖分,该过程包括以下步骤:
(6-1)对所建的模型进行粗化网格剖分,并对x向、y向、z向的缩放比例进行设置,对穿过删除的控制实体平滑下的迭代数、最大单元深度进行设置;
整个模型建立完毕;
(7)计算:借助有限元仿真软件仿真求解;
(8)计算完毕后在结果中查看并分析计算结果;
(9)根据计算结果优化设计参数。
微流道感应腔是微流道加速度计的核心部件,其内部结构如图1所示,由电解液密封腔、电解液、电极(阳极和阴极)、绝缘隔层以及导线组成:电解液密封腔内填充一定浓度的电解液(通常是水基盐溶液),阳极和阴极安装在电解密封液腔内并被绝缘隔层间隔开,电极两端施加一定电势,并与电解液密封腔外的信号采集处理电路相连。其中电极与绝缘隔层通常被封装在一起,被称为“敏感元件”。敏感元件是微流道感应腔的核心部件,也微流道加速度计的核心敏感部件。
微流道加速度计的工作原理基于外界运动引起的微流道感应腔内电解液的离子浓度差的变化:在密闭腔体内封装电解液(以水基盐溶液为例),并将一对电极放置于电解液中。当电极两侧无电势差时,电解液内的某种有效离子的浓度保存恒定;当电极两侧有电势差时,会在阴极和阳极产生电极反应,从而导致在电极两侧形成上述的某种有效离子的浓度差。该浓度差会随密闭腔体的运动或姿态而发生改变,进而在电极直接产生电流变化。通过测量该电流变化即可获得密闭腔体此时的运动或姿态变化情况。
以下结合附图对本发明做进一步详细的描述:
本发明提出的一种利用有限元分析软件-COMSOL Multiphysics为微流道加速度计建立原理模型的方法,其流程图如图2所示,包括以下各步骤:
(1)建立包含能斯特-普朗克方程和层流场的三维动态模型,该过程包括以下步骤:
(1-1)在COMSOL Multiphysics 4.3a建立一个包含能斯特-普朗克方程和层流场的三维动态模型。
其中,能斯特-普朗克方程为:
,其中,是浓度,是时间,是扩散系数,是电荷数,是迁移率, 是电势,F是Faraday常量,是速度,是变量生产条件。
,其中,F是Faraday常量,是电荷数,是变量生产条件。
,其中,是电荷数,是浓度。
,其中,是通量,是扩散系数,是浓度, 是电荷数,是迁移率, 是电势,F是Faraday常量,是速度向量。
,其中,F是Faraday常量,是电荷数,是扩散系数,是浓度, 是迁移率, 是电势。
层流场方程为:
,其中,是密度,是速度向量,是时间,是圆柱体长度, 是温度,F是Faraday常量。
,其中,是密度,是速度向量。
(1-2)在模型中,定义3种变量,分别用 、 、表示。
(1-3)在模型中,将预置求解定义为瞬态模式。
(2)根据微流道感应腔设计尺寸,在步骤(1)所建模型内画出所设计的微流道加速度计微流道感应腔三维结构,该过程包括以下步骤:
(2-1)在模型中,将长度单元设置为μm,角度单元设置为Degrees。
(2-2)在模型中,以坐标原点(0,0)为起点,沿x轴正方向画直径为 μm,高度依次为 μm, μm,μm, μm, μm, μm, μm, μm, μm的圆柱体,分别命名为:a,b,c,d,e,f,g,h,i。
(3)输入所设计微流道加速度计电解液导电率、相对介电常数、电解液密度、粘滞度这四个参数值。
(4)在能斯特-普朗克方程中定义电极、电势及边界条件和约束方程。
电气绝缘:将圆柱体a, c, e, g, i定义为绝缘体,约束方程为:
无通量的约束方程为:,其中,
电场迁移的约束方程为:
设定圆柱体b, h为阳极,电势值为U1,圆柱体d,f为阴极,电势值为U2。
(5)层流场中定义入口、出口、流速及边界条件和约束方程。
入口方程为:,边界条件定义为速度,且定义为方向流入速度,为设一初始值。
出口方程为:,,其中,是温度,是压力,是动力粘度,是速度场,是圆柱体长度,边界条件定义为压力,无粘滞应力,为设一初始值。
边界条件定义为无滑移。
(6)对求解域进行网格剖分,该过程包括以下步骤:
(6-1)对所建的模型进行粗化网格剖分,并对x向、y向、z向的缩放比例进行设置,对穿过删除的控制实体平滑下的迭代数、最大单元深度进行设置。整个模型建立完毕。
(7)计算。
由于以上方程都是多元偏微分方程,很难直接求解,因此借助有限元仿真软件仿真求解。
(8)计算完毕后在结果中查看并分析计算结果。
(9)根据计算结果优化设计参数。