CN113221420B - 一种基于ansys软件的传感器电容误差评估方法 - Google Patents
一种基于ansys软件的传感器电容误差评估方法 Download PDFInfo
- Publication number
- CN113221420B CN113221420B CN202110533188.6A CN202110533188A CN113221420B CN 113221420 B CN113221420 B CN 113221420B CN 202110533188 A CN202110533188 A CN 202110533188A CN 113221420 B CN113221420 B CN 113221420B
- Authority
- CN
- China
- Prior art keywords
- capacitance
- sensor
- error
- simulation
- value
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/04—Constraint-based CAD
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Measuring Fluid Pressure (AREA)
Abstract
本发明属于传感器和计算机技术领域,具体的说是涉及一种基于ANSYS软件的传感器电容误差评估方法。本发明首先建立传感器仿真模型,划分网格,施加压力载荷,进行静态结构力学仿真,得到中心挠度仿真值W,计算中心挠度理论值W0;判断中心挠度误差是否满足要求,如果不满足则继续优化网格质量;采取微元电容离散法来计算电容仿真值C,计算电容理论值C0;比较电容理论值C0和仿真值C,判断电容误差是否满足要求,如果满足,则误差评估结束;如果不满足,则继续判断仿真电容误差能否满足要求;如果仿真电容误差满足要求,继续优化网格质量,如果不满足,继续优化电容计算方法。本发明克服了目前传感器误差评估中仅有的应力和形变的误差评估的不足之处。
Description
技术领域
本发明属于传感器和计算机技术领域,具体的说是涉及一种基于ANSYS软件的传感器电容误差评估方法。
背景技术
电容式压力传感器具有灵敏度高、温度漂移小、结构更健壮、功耗低、对环境影响的灵敏度低等优点,主要应用于汽车系统、工业控制、环境监测、航天航空等领域,工作环境较为复杂,容易受到高温、盐雾、振动环境载荷影响,很难在直接监测传感器在复杂情况下工作能力,对传感器进行有限元仿真和误差评估,有利于直接从仿真分析中准确了解传感器的工作能力和工作状态。
根据传感器的工作环境,在确保传感器模型的边界约束条件合理的情况下,施加静态压力载荷,对传感器进行静态结构力学仿真,此时传感器的误差主要来源于网格单元的离散误差和电容计算方法的误差。目前,传感器的误差分析主要包括有限元仿真方法误差分析、奇异点单元附件网格离散误差分析和“比较分析”方法,重点研究应力和形变仿真结果误差。由于ANSYS无法直接得到电容仿真结果,需要通过传感器的形变,采用高精度的电容计算方法计算出传感器的电容值,传感器的电容误差分析不仅包括有限元仿真方法的误差和网格离散化误差,还应包括电容计算方法的误差,基于ANSYS仿真软件的传感器电容误差评估方法还未出现。在ANSYS有限元分析软件中,传感器电容计算方法一般通过仿真得到中心挠度和接触半径,再采用挠度理论来计算传感器的仿真电容,由于在接触工作状态下,传感器敏感膜和衬底出现较小的接触穿透,很难找准接触半径,挠度理论计算出的电容准确度低。
发明内容
针对上述问题,本发明提出一种基于ANSYS软件的传感器电容误差评估方法,结合中心挠度误差评估和电容误差评估,判断电容误差是否满足要求和误差的主要来源,通过减小网格离散误差和计算方法误差,从而获得较为精确电容值。
本发明的技术方案是:
一种基于ANSYS软件的传感器电容误差评估方法,包括以下步骤:
S1、根据传感器电容结构参数和工作环境,在ANSYS软件上建立传感器仿真模型,并确定边界约束条件;
S2、对传感器仿真模型进行网格划分;
S3、施加压力载荷,进行静态结构力学仿真,得到中心挠度仿真值W,并通过挠度理论计算中心挠度理论值W0;
S4、比较中心挠度理论值W0和仿真值W,判断中心挠度误差是否满足要求,若是,则进入步骤S5,否则回到步骤S2重新进行网格划分;
S5、提取传感器发生形变后网格点坐标值,采取微元电容离散法来计算电容仿真值C,并通过挠度理论计算电容理论值C0;
S6、比较电容理论值C0和仿真值C,判断电容误差是否满足要求,若是,则误差评估结束,得到满足误差要求的电容值,否则,进入步骤S7;
S7、通过挠度理论计算中心挠度为W时的传感器电容C1,比较微元电容离散法计算的电容仿真值C和挠度理论计算出的电容仿真值C1,判断仿真电容误差能否满足要求,若是,则回到步骤S2重新进行网格划分,否则,进入步骤S5,采用MATLAB中meshgrid和griddata函数将电容器离散成更多的元电容,继续计算电容仿真值C。
进一步的,步骤S5中所述微元电容离散法为映射法和加密法,使用加密法计算传感器空腔的电容,使用映射法求解传感器环带电容;加密法和映射法的具体方式可参考申请号为202010573493.3的中国专利申请,这种计算方法是利用微元思想将电容器离散成有限个电容之和,(简称“微元电容离散法”),该方法能比较精确计算传感器的电容仿真值。
本发明的有益效果是:本发明首先评估中心挠度误差,判断网格离散误差是否满足要求,再对传感器电容误差进行评估,判断仿真电容求解方法误差是否满足要求;本发明补充了目前传感器误差评估中仅有的应力和形变误差评估方法,详细描述了传感器电容误差评估的具体流程;通过挠度理论计算法和微元离散电容法求解仿真电容,判断电容误差的主要来源,为提升仿真结果准确度提供思路;本发明创新性提出了电容式压力传感器的电容误差评估方法,为ANSYS软件电容式压力传感器电容误差提供指导。
附图说明
图1为传感器电容误差评估方法的总体流程示意图;
图2为电容式压力传感器电容误差评估的具体流程示意图;
图3为电容式压力传感器仿真模型示意图;
图4为电容式压力传感器中心挠度值和中心挠度误差结果图;
图5为电容式压力传感器电容值和电容误差结果图。
具体实施方式
下面结合附图和实施例,详细描述本发明的技术方案。
图1是传感器电容误差评估方法的总体流程示意图,首先建立传感器仿真模型并确定边界约束条件;网格划分和优化网格质量;施加压力载荷,进行静态结构力学仿真,得到中心挠度仿真值W,并计算中心挠度理论值W0;比较在中心挠度理论值W0和仿真值W,判断中心挠度误差是否满足要求,如果不满足要求则继续优化网格质量;提取传感器发生形变后的网格点坐标值,采取微元电容离散法来计算电容仿真值C,并通过挠度理论计算电容理论值C0;比较电容理论值C0和仿真值C,判断电容误差是否满足要求,如果满足,则误差评估结束;如果电容误差不满足要求,通过挠度理论计算中心挠度为W时的仿真电容C1,比较微元电容离散法计算的电容仿真值C和挠度理论计算的电容仿真值C1,并判断仿真电容误差能否满足要求;如果仿真电容误差满足要求,误差主要来源网格离散误差,则继续优化网格质量,如果不满足,误差来源于电容计算方法,则继续优化电容计算方法。如图2所示,本发明的具体实施步骤如下:
步骤1、建立传感器仿真模型并确定边界约束条件,电容式压力传感器仿真模型如图3所示,主要由SiC衬底、SiC敏感膜和SiO2绝缘层、SiO2钝化层、空腔和欧姆电极组成,空腔高度为6.3μm,敏感膜半径为1400μm,敏感膜厚度为26μm,传感器的边界约束条件主要由实际工作环境确定。
步骤2、采用高精度的网格划分方法,对电容式压力传感器仿真模型进行网格划分,并优化网格质量,使得网格离散化误差尽可能小。
步骤3、施加载荷,进行静态力学仿真,得到中心挠度仿真值W,同时通过挠度理论计算出中心挠度理论值W0,如图4所示,得到在未接触状态中心挠度曲线。
步骤4、比较中心挠度仿真值W和理论值W0,图4所示的中心挠度误差小于10%,判断中心挠度误差是否满足要求,如果满足要求,则继续优化网格质量。
步骤5、利用挠度理论计算出传感器电容理论值C0,同时,提取出发生形变后传感器敏感膜和绝缘层的坐标值,采用加密法和映射法,通过MATLAB软件计算出传感器的仿真值C,如图5所示,得到传感器的电容曲线。
步骤6、比较电容式压力传感器的电容理论值和仿真值,图5中电容误差小于1%,判断误差是否满足要求,如果满足要求,则传感器的电容误差评估结束。
步骤7、当电容误差不满足要求时,通过挠度理论计算中心挠度为W时,传感器的仿真电容C1。
步骤8、比较用挠度理论计算的C1和用微元离散电容法求解的C,判断仿真电容误差能否满足要求。如果满足要求,则说明仿真电容误差主要来源网格离散误差;如果不满足要求,则说明仿真电容误差主要来源计算方法误差,继续优化计算方法,将传感器离散成更多的电容之和,再次求解传感器电容仿真值。
本次发明克服了ANSYS误差评估中仅有的应力和形变误差评估方法的不足之处,结合中心挠度误差评估和电容误差评估,判断误差的主要来源,得到符合误差要求的电容值,为电容式压力传感器电容误差评估提供指导。
Claims (1)
1.一种基于ANSYS软件的传感器电容误差评估方法,其特征在于,包括以下步骤:
S1、根据传感器电容结构参数和工作环境,在ANSYS软件上建立传感器仿真模型,并确定边界约束条件;所述传感器仿真模型为电容式压力传感器仿真模型,由SiC衬底、SiC敏感膜和SiO2绝缘层、SiO2钝化层、空腔和欧姆电极组成,其中空腔高度为6.3μm,敏感膜半径为1400μm,敏感膜厚度为26μm,传感器的边界约束条件由实际工作环境确定;
S2、对传感器仿真模型进行网格划分;
S3、施加压力载荷,进行静态结构力学仿真,得到中心挠度仿真值W,并通过挠度理论计算中心挠度理论值W0;
S4、比较中心挠度理论值W0和仿真值W,判断中心挠度误差是否满足要求,若是,则进入步骤S5,否则回到步骤S2重新进行网格划分;
S5、提取传感器发生形变后网格点坐标值,采取微元电容离散法来计算电容仿真值C,并通过挠度理论计算电容理论值C0;所述微元电容离散法为映射法和加密法,使用加密法计算传感器空腔的电容,使用映射法求解传感器环带电容;
S6、比较电容理论值C0和仿真值C,判断电容误差是否满足要求,若是,则误差评估结束,得到满足误差要求的电容值,否则,进入步骤S7;
S7、通过挠度理论计算中心挠度为W时的传感器电容C1,比较微元电容离散法计算的电容仿真值C和挠度理论计算出的电容仿真值C1,判断仿真电容误差能否满足要求,若是,则回到步骤S2重新进行网格划分,否则,进入步骤S5,采用MATLAB中meshgrid和griddata函数将电容器离散成更多的元电容,继续计算电容仿真值C。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110533188.6A CN113221420B (zh) | 2021-05-17 | 2021-05-17 | 一种基于ansys软件的传感器电容误差评估方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110533188.6A CN113221420B (zh) | 2021-05-17 | 2021-05-17 | 一种基于ansys软件的传感器电容误差评估方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113221420A CN113221420A (zh) | 2021-08-06 |
CN113221420B true CN113221420B (zh) | 2023-03-03 |
Family
ID=77092279
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110533188.6A Active CN113221420B (zh) | 2021-05-17 | 2021-05-17 | 一种基于ansys软件的传感器电容误差评估方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113221420B (zh) |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103964370A (zh) * | 2013-01-29 | 2014-08-06 | 北京大学 | 一种电容式压力传感器的制备方法 |
CN107515060A (zh) * | 2017-08-29 | 2017-12-26 | 南京工业大学 | 一种电容式压力传感器、线性补偿方法及制备方法 |
CN111737900A (zh) * | 2020-06-22 | 2020-10-02 | 电子科技大学 | 基于有限元仿真后处理的传感器电容值计算方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103162894A (zh) * | 2013-01-14 | 2013-06-19 | 楚雄师范学院 | 一种电容压力传感器 |
-
2021
- 2021-05-17 CN CN202110533188.6A patent/CN113221420B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103964370A (zh) * | 2013-01-29 | 2014-08-06 | 北京大学 | 一种电容式压力传感器的制备方法 |
CN107515060A (zh) * | 2017-08-29 | 2017-12-26 | 南京工业大学 | 一种电容式压力传感器、线性补偿方法及制备方法 |
CN111737900A (zh) * | 2020-06-22 | 2020-10-02 | 电子科技大学 | 基于有限元仿真后处理的传感器电容值计算方法 |
Non-Patent Citations (5)
Title |
---|
SiC高温压力传感器综合应力仿真及可靠性分析;杨荣森;《中国优秀硕士学位论文全文数据库(信息科技辑)》;20200715;I135-163 * |
Stress investigation of the AlGaN/GaN micromachined circular diaphragms of a pressure sensor;J Dzuba1等;《Journal of Micromechanics and Microengineering》;20140905;第25卷(第1期);1-7 * |
厚膜陶瓷电容式压力传感器设计与制备;李鹏等;《传感器世界》;20061219(第12期);12-15 * |
基于ANSYS软件均布载荷作用下方板挠度分析;谢中敏等;《数字技术与应用》;20180905(第09期);50-51 * |
新型厚膜陶瓷电容式压力传感器感压元件研究;唐力强等;《传感器与微系统》;20060620(第6期);39-42 * |
Also Published As
Publication number | Publication date |
---|---|
CN113221420A (zh) | 2021-08-06 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103076131B (zh) | 用于测量大型机械臂大力与小力矩的六维力与力矩传感器 | |
CN107092759B (zh) | 基于重力坝坝基参数反演的坝体位移监测点优化布置方法 | |
GB2593851A (en) | Bearing reaction influence line curvature-based continuous beam damage identification method | |
CN101592552A (zh) | 预测汽车后悬架多轴疲劳的方法 | |
CN106125548A (zh) | 工业机器人动力学模型参数辨识方法 | |
CN103267507A (zh) | 基于有限元分析提取机械结构平面的平面度误差的方法 | |
Chwalowski et al. | Preliminary computational analysis of the (hirenasd) configuration in preparation for the aeroelastic prediction workshop | |
CN113221420B (zh) | 一种基于ansys软件的传感器电容误差评估方法 | |
CN103063333A (zh) | 一种飞机框类零件装配应力检测与控制方法 | |
CN106021689A (zh) | 一种计算结构变形后接触力的方法 | |
CN105550383B (zh) | 一种非定常气动力测量试验系统的设计方法 | |
CN109933883A (zh) | 一种基于多元线性回归的阻性电流增长率计算方法 | |
WO2023226566A1 (zh) | 臂架的控制方法、控制系统及工程机械 | |
CN102270249B (zh) | 识别零部件特征频率的方法 | |
CN110146290A (zh) | 一种基于定位普遍法则的滚动轴承外圈缺陷定位诊断方法 | |
CN104776947A (zh) | 一种机器人腕力传感器温度漂移的补偿方法 | |
CN106053984B (zh) | 确定防滑刹车控制装置高温故障分布的方法 | |
Sah et al. | Process monitoring in stamping operations through tooling integrated sensing | |
CN112765863B (zh) | 基于位姿依赖特性和交叉耦合项的机器人刀尖频响预测方法和系统 | |
CN112179550B (zh) | 一种四支点压电测力仪的优化布置方法 | |
CN112711809A (zh) | 一种舵面载荷筛选方法 | |
CN113358311B (zh) | 基于有限测点和振动模态的板/梁结构横向位移确定方法 | |
CN107168573B (zh) | 一种多压电边界排列三维触压感知方法 | |
CN104462790A (zh) | 疲劳耐久分析的自由表面方法 | |
CN105424365A (zh) | 一种滚动轴承故障迁移轨迹可视化的方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |