CN112417777B - 房水引流装置及其结构优化方法、系统、终端及介质 - Google Patents
房水引流装置及其结构优化方法、系统、终端及介质 Download PDFInfo
- Publication number
- CN112417777B CN112417777B CN202011267286.1A CN202011267286A CN112417777B CN 112417777 B CN112417777 B CN 112417777B CN 202011267286 A CN202011267286 A CN 202011267286A CN 112417777 B CN112417777 B CN 112417777B
- Authority
- CN
- China
- Prior art keywords
- aqueous humor
- fluid
- flow
- fluid flow
- outlet
- 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
- 210000001742 aqueous humor Anatomy 0.000 title claims abstract description 131
- 238000005457 optimization Methods 0.000 title claims abstract description 95
- 238000000034 method Methods 0.000 title claims abstract description 88
- 239000012530 fluid Substances 0.000 claims abstract description 160
- 238000004364 calculation method Methods 0.000 claims abstract description 40
- 239000011148 porous material Substances 0.000 claims abstract description 17
- 238000009826 distribution Methods 0.000 claims description 33
- 238000004422 calculation algorithm Methods 0.000 claims description 24
- 230000015654 memory Effects 0.000 claims description 22
- 238000004088 simulation Methods 0.000 claims description 16
- 238000004590 computer program Methods 0.000 claims description 11
- 230000008569 process Effects 0.000 claims description 9
- 239000006185 dispersion Substances 0.000 claims description 5
- 238000012805 post-processing Methods 0.000 claims description 5
- 238000007781 pre-processing Methods 0.000 claims description 5
- 238000013433 optimization analysis Methods 0.000 claims description 4
- 208000010412 Glaucoma Diseases 0.000 abstract description 16
- 238000002513 implantation Methods 0.000 abstract description 11
- 230000000638 stimulation Effects 0.000 abstract description 5
- 238000013461 design Methods 0.000 description 34
- 230000004907 flux Effects 0.000 description 9
- 238000010586 diagram Methods 0.000 description 8
- 238000012937 correction Methods 0.000 description 7
- 210000001519 tissue Anatomy 0.000 description 7
- 238000005516 engineering process Methods 0.000 description 5
- 206010030043 Ocular hypertension Diseases 0.000 description 4
- 210000004027 cell Anatomy 0.000 description 4
- 230000008878 coupling Effects 0.000 description 4
- 238000010168 coupling process Methods 0.000 description 4
- 238000005859 coupling reaction Methods 0.000 description 4
- 238000004458 analytical method Methods 0.000 description 3
- 230000006870 function Effects 0.000 description 3
- 230000004410 intraocular pressure Effects 0.000 description 3
- 230000004048 modification Effects 0.000 description 3
- 238000012986 modification Methods 0.000 description 3
- 206010016654 Fibrosis Diseases 0.000 description 2
- 230000008901 benefit Effects 0.000 description 2
- 210000005252 bulbus oculi Anatomy 0.000 description 2
- 230000007547 defect Effects 0.000 description 2
- 238000009792 diffusion process Methods 0.000 description 2
- 230000004761 fibrosis Effects 0.000 description 2
- 230000006872 improvement Effects 0.000 description 2
- 230000001575 pathological effect Effects 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 230000001360 synchronised effect Effects 0.000 description 2
- 238000009827 uniform distribution Methods 0.000 description 2
- 206010003694 Atrophy Diseases 0.000 description 1
- 201000004569 Blindness Diseases 0.000 description 1
- 208000010415 Low Vision Diseases 0.000 description 1
- 210000002159 anterior chamber Anatomy 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 230000037444 atrophy Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000011438 discrete method Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 210000001508 eye Anatomy 0.000 description 1
- 208000030533 eye disease Diseases 0.000 description 1
- 230000003176 fibrotic effect Effects 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 102000051759 human factor J Human genes 0.000 description 1
- 108700008420 human factor J Proteins 0.000 description 1
- 238000011835 investigation Methods 0.000 description 1
- 210000001365 lymphatic vessel Anatomy 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 210000001328 optic nerve Anatomy 0.000 description 1
- 238000005192 partition Methods 0.000 description 1
- 230000037361 pathway Effects 0.000 description 1
- 230000002085 persistent effect Effects 0.000 description 1
- 230000002035 prolonged effect Effects 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 230000036573 scar formation Effects 0.000 description 1
- 230000037390 scarring Effects 0.000 description 1
- 230000011218 segmentation Effects 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 238000001356 surgical procedure Methods 0.000 description 1
- 230000004382 visual function Effects 0.000 description 1
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/28—Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
- G16H50/00—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics
- G16H50/50—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics for simulation or modelling of medical disorders
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/08—Fluids
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Public Health (AREA)
- Medical Informatics (AREA)
- Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Geometry (AREA)
- Biomedical Technology (AREA)
- Pure & Applied Mathematics (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Mathematical Optimization (AREA)
- General Engineering & Computer Science (AREA)
- Mathematical Analysis (AREA)
- Fluid Mechanics (AREA)
- Mathematical Physics (AREA)
- Data Mining & Analysis (AREA)
- Databases & Information Systems (AREA)
- Computing Systems (AREA)
- Pathology (AREA)
- Algebra (AREA)
- Epidemiology (AREA)
- General Health & Medical Sciences (AREA)
- Primary Health Care (AREA)
- External Artificial Organs (AREA)
Abstract
本发明提供了一种房水引流装置及其结构优化方法和系统，包括：S1，确定目标对象，明确结构优化要求；S2，创建或修正所述目标对象的3D几何模型，建立流体流通域模型；S3，对流体流通域模型进行网格划分以及有限体积数值计算；S4，分析流体流通域模型的流动特性，如满足结构优化要求，完成结构优化；如不满足结构优化要求，重新执行上述S2和S3，直至满足结构优化要求。同时提供了一种相应的终端和存储介质。本发明目的在于实现沟槽流道或孔道各出口流体流量分配基本均匀，减少青光眼房水引流装置植入术对术区微环境的扰动或刺激，提高青光眼房水引流装置植入术手术成功率。
Description
技术领域
本发明涉及医工交叉技术领域，具体地，涉及一种房水引流装置及其结构优化方法、系统、终端及介质。
背景技术
青光眼是全球第二大致盲性眼病。病理性眼压增高是青光眼的主要危险因素，眼内压间断或持续性升高会给眼球各部分组织和视功能带来损害，将导致视神经萎缩、视野缩小、视力减退，甚至失明。在青光眼治疗中，最直接的方法是降低眼压，而其中尤以难治性青光眼的眼压控制最为棘手。针对难治性青光眼，临床一线常采用房水引流器械植入术，即通过房水引流器械的植入，建立房水引流和扩散的人工通路，将房水引流到眼球外部，再通过各种组织管道(如淋巴管或静脉血管等)扩散循环到全身，从而降低眼内压。目前手术失败的主要原因是术区组织纤维化及瘢痕形成，阻碍房水的引流和扩散；据统计，房水引流装置植入术五年手术成功率只有不到50％。
针对上述临床面临的困局，技术人员在前期研发了一种青光眼房水引流装置(CN204319045U)，以期通过沟槽设计减少房水对组织微环境的扰动或刺激，并进一步减少组织纤维化和瘢痕形成。然而，实验证明沟槽出口处的房水流量分配不均，流量较多的沟槽出口房水易于冲击周围组织，刺激附近的组织纤维化，包裹引流装置，降低房水引流装置植入术的手术成功率。因此，沟槽设计的房水引流装置的结构设计合理性和科学性有待提高，有必要从理论计算方面对结构优化给出方向和方案。
目前没有发现同本发明类似技术的说明或报道，也尚未收集到国内外类似的资料。
发明内容
本发明针对现有技术中存在的上述不足，提供了一种房水引流装置及其结构优化方法、系统、终端及介质。
本发明是通过以下技术方案实现的。
一种房水引流装置的结构优化方法，包括：
S1，确定目标对象，明确结构优化要求；
S2，创建或修正所述目标对象的3D几何模型，建立流体流通域模型；
S3，对流体流通域模型进行网格划分以及有限体积数值计算；
S4，分析流体流通域模型的流动特性，如满足结构优化要求，完成结构优化；如不满足结构优化要求，重新执行上述S2和S3，直至满足结构优化要求。
优选地，所述S1中，目标对象包括：引流管以及具有弧度的引流盘，所述引流盘表面具有放射状沟槽流道和/或孔道；
所述结构优化要求为：所述引流盘表面的各沟槽流道和/或孔道出口的流体流量分配基本均匀。
优选地，所述S2中，创建或修正所述目标对象的3D几何模型，包括如下任一项或任意多项：
-引流盘长轴和/或短轴的尺寸；
-引流盘厚度和/或上下表面曲率半径；
-沟槽流道的宽度、深度和/或底曲面曲率半径；
-孔道的直径和/或位置；
-引流管的管壁厚度和/或内径。
优选地，所述S2中，建立流体流通域模型的方法，包括：
基于房水过流区域确定流体流通域的几何结构，采用三维建模方法创建或重建流体流通域模型。
优选地，所述S3中，对流体流通域模型进行网格划分的方法，包括：
采用计算流体动力学的数值模拟软件的前处理模块对所述流体流通域模型进行网格划分，基于所述流体流通域模型的不规则性，采用非结构化网格对所述流体流通域模型进行空间离散。
优选地，所述S3中，对流体流通域模型进行有限体积数值计算的方法，包括：
采用计算流体动力学的数值模拟软件的流体计算模块对流体属性、流体流通域边界、模拟过程中的算法、精度和网格划分后的空间离散格式进行参数设置，并根据设置的参数，采用有限体积数值计算的方法获取房水引流装置的出入口流量数据。
优选地，所述对流体流通域模型进行有限体积数值计算的方法中，还包括如下任意一项或任意多项：
-所述流体属性包括：流体密度和/或粘度；
-所述流体流通域边界包括：压力入口和/或压力出口；
-所述模拟过程中的算法采用SIMPLE算法；
-所述网格划分后的空间离散格式采用second-order upwind格式；
-所述有限体积法的数值算法采用非耦合隐式算法。
优选地，所述S4中，采用计算流体动力学的数值模拟软件的后处理模块分析流体流通域模型的流动特性。
优选地，采用沟槽流道或孔道出口与入口的流量比来表示房水引流装置的流体流量分配情况，选用平均差公式：计算平均差，用于衡量各出口流量比偏离所述平均值的程度，定量表述沟槽流道或孔道出口的流体流量分配均匀度；其中，xi表示为各出口与入口的流量比，/>为相同入口压力下各出口流量比的算术平均数，n为出口数量。
优选地，基于所述沟槽流道或孔道出口的流体流量分配均匀度，综合评价所述结构优化方法的合理性和有效性。
根据本发明的第二个方面，提供了一种房水引流装置的结构优化系统，包括：
目标获取模块，该模块用于确定目标对象，明确结构优化要求；
流体流通域模型模块，该模块用于创建或修正所述目标对象的3D几何模型，建立流体流通域模型；
模型参数设置模块，该模块用于对流体流通域模型进行网格划分以及有限体积数值计算；
结构优化分析模块，该模块用于分析流体流通域模型的流动特性，当流动特性满足结构优化要求时，完成对目标对象的结构优化。
根据本发明的第三个方面，提供了一种房水引流装置，采用上述任一项所述的结构优化方法进行结构优化。
根据本发明的第四个方面，提供了一种终端，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，所述处理器执行所述程序时可用于执行上述任一项所述的方法。
根据本发明的第五个方面，一种计算机可读存储介质，其上存储有计算机程序，该程序被处理器执行时可用于执行上述任一项所述的方法。
由于采用了上述技术方案，本发明与现有技术相比，具有如下有益效果：
本发明提供的房水引流装置及其结构优化方法、系统、终端及介质，初始引流装置创建几何结构时的具体尺寸是一种房水引流装置的初步结构参数，以均匀流场为目标，最终经过一系列优化后得到的优化模型的尺寸为一种结构合理的房水引流装置结构参数，建立并验证了一种行之有效的房水引流装置结构设计优化方法。
本发明提供的房水引流装置及其结构优化方法、系统、终端及介质，基于CFD技术的支持，对房水引流装置几何结构的合理性和有效性进行数值计算，并提出可能的优化改进意见，最终实现沟槽流道或孔道各出口流体流量分配基本均匀，减少青光眼房水引流装置植入术对术区微环境的扰动或刺激，提高青光眼房水引流装置植入术手术成功率。
实施本发明的任一产品并不一定需要同时达到以上所述的所有优点。
附图说明
通过阅读参照以下附图对非限制性实施例所作的详细描述，本发明的其它特征、目的和优点将会变得更明显：
图1为本发明一优选实施例中房水引流装置的结构优化方法的流程图；
图2-1为本发明一优选实施例中房水引流装置的结构优化方法的目标对象；
图2-2为本发明一优选实施例中房水引流装置的剖视图；
图2-3为本发明一优选实施例中房水引流装置中流体计算域的三维网格划分图及其局部放大图；
图2-4为本发明一优选实施例中房水引流装置的流体流道图；
图2-5为本发明一优选实施例中房水引流装置的计算结果图，其中：(a)为各出口截面积；(b)为入口压力15mmHg时各出口流量比例；(c)为各出口流量平均差；(d)为入口压力21mmHg时各出口流量比例；
图3-1为本发明一优选实施例中第一种改进型房水引流装置的流体流道图；
图3-2为本发明一优选实施例中第一种改进型房水引流装置的计算结果图：(a)为各出口截面积；(b)为入口压力15mmHg时各出口流量比例；(c)为各出口流量平均差；(d)为入口压力21mmHg时各出口流量比例；
图4-1为本发明一优选实施例中第二种改进型房水引流装置的流体流道图；
图4-2为本发明一优选实施例中第二种改进型房水引流装置的计算结果图：(a)为各出口截面积；(b)为入口压力15mmHg时各出口流量比例；(c)为各出口流量平均差；(d)为入口压力21mmHg时各出口流量比例；
图5-1为本发明一优选实施例中第三种改进型房水引流装置的流体流道图；
图5-2为本发明一优选实施例中第三种改进型房水引流装置的计算结果图：(a)为各出口截面积；(b)为入口压力15mmHg时各出口流量比例；(c)为各出口流量平均差；(d)为入口压力21mmHg时各出口流量比例；
图6-1为本发明一优选实施例中第四种改进型房水引流装置的几何模型；
图6-2为本发明一优选实施例中第四种改进型房水引流装置的剖视图；
图6-3为本发明一优选实施例中第四种改进型房水引流装置的流体流道图；
图6-4为本发明一优选实施例中第四种改进型房水引流装置的计算结果图：(a)为各出口截面积；(b)为入口压力15mmHg时各出口流量比例；(c)为各出口流量平均差；(d)为入口压力21mmHg时各出口流量比例；
图6-5为本发明一优选实施例中第四种改进型房水引流装置的流线分布图；
图7为本发明一优选实施例中房水引流装置的结构优化的组成模块示意图。
具体实施方式
下面对本发明的实施例作详细说明：本实施例在以本发明技术方案为前提下进行实施，给出了详细的实施方式和具体的操作过程。应当指出的是，对本领域的普通技术人员来说，在不脱离本发明构思的前提下，还可以做出若干变形和改进，这些都属于本发明的保护范围。
本发明一实施例提供了一种房水引流装置的结构优化方法，该方法拟通过CFD(Computational Fluid Dynamics，计算流体动力学)模拟计算与数值分析技术建立一种房水引流装置的结构优化方法，系统化计算和优化房水引流装置的表面结构特征，以期实现沟槽状房水引流装置均匀引流的目的，提高房水引流装置植入术手术成功率。
本实施例旨在克服现有房水引流装置存在的引流不均匀的缺点，提供的房水引流装置的结构优化方法，实现房水引流装置的结构优化，使其沟槽流道或孔道出口的流体流量分配基本均匀。
本实施例提供的种房水引流装置的结构优化方法，包括如下步骤：
步骤(1)确定目标对象，明确结构优化要求；
步骤(2)创建或修正所述目标对象的3D几何模型，建立流体流通域模型；
步骤(3)对流体流通域模型进行网格划分以及有限体积数值计算；
步骤(4)分析流体流通域模型的流动特性，如满足设计要求，完成结构设计优化；若不满足设计要求，重新提出结构设计优化方案，执行上述步骤(2)和步骤(3)，直至满足设计要求。
作为一优选实施例，步骤(1)中，目标对象包括一个引流管及一个有弧度的引流盘，所述引流盘表面有放射状沟槽流道和/或孔道。优化要求为引流装置各沟槽流道或孔道出口的流体流量分配基本均匀。
作为一优选实施例，步骤(2)中，目标对象的3D几何模型的创建或修正包括引流盘长轴和短轴的尺寸、引流盘厚度和上下表面曲率半径，所述沟槽的宽度、槽的深度、槽底曲面曲率半径，所述孔道的直径和位置，引流管的管壁厚度及内径。
作为一优选实施例，步骤(2)中，几何模型被建立或修正后，基于房水过流区域确定流域几何结构，采用三维建模软件创建或重建流体流通域。所述的流体流通域包括引流管和引流槽，即房水过流区域。
作为一优选实施例，步骤(3)中，采用流体力学领域的数值模拟软件的前处理模块对所述的流通域进行网格划分，基于所述几何模型的不规则性，采用非结构化网格对流通域进行空间离散。
作为一优选实施例，步骤(3)中，采用流体力学领域的数值模拟软件的流体分析模块实现所述有限体积数值计算，其中包括：对流体属性、流体流通域边界、模拟过程中的算法、精度和网格分割后的空间离散格式进行参数设置，根据设置的参数，采用有限体积数值计算的方法获取房水引流装置的出入口流量数据。
作为一优选实施例，对流体流通域模型进行有限体积数值计算的方法中，还包括如下任意一项或任意多项：
-流体属性包括：流体密度和/或粘度；
-流体流通于边界包括：压力入口和/或压力出口；
-模拟过程中求解器的设置采用SIMPLE算法；
-网格划分后的空间离散格式采用second-order upwind格式；
-所述有限体积法的数值算法采用非耦合隐式算法。
进一步地，模拟过程中的算法包含但不限于SIMPLE算法。用于数值模拟计算，得出一组出入口流量数据，用于后处理模块分析验证结构设计的合理性。
进一步地，有限体积法包括三种数值算法：非耦合隐式算法、耦合显式算法和耦合隐式算法，本实施例优选采用非耦合隐式算法。
作为一优选实施例，步骤(4)中，采用流体力学领域的数值模拟软件的后处理模块分析所述的计算数据，如满足设计要求，即完成结构优化，若不满足设计要求，提出进一步结构设计优化方案，返回步骤(2)、(3)，重新开展仿真计算，直至满足设计要求。
作为一优选实施例，本实施例所提供的方法，采用沟槽流道或孔道出口与入口的流量比来表示所述房水引流装置的流体流量分配情况，选用平均差公式：(xi表示为各出口与入口的流量比，/>为相同入口压力下各出口流量比的算术平均数，n为出口数量)计算平均差，用于衡量各出口流量比偏离所述平均值的程度，定量表述沟槽流道或孔道出口的流体流量分配均匀度。
作为一优选实施例，基于沟槽流道或孔道出口的流体流量分配均匀度，综合评价该结构优化方法的合理性和有效性。
本实施例提供的房水引流装置的结构优化方法，将房水引流装置的初步结构参数，通过该优化方法的处理，最终得到具有最合理结构的房水引流装置的一系列参数。
在本发明部分实施例中：
CFD数值模拟方法
(1)控制方程
房水在青光眼引流装置中的流动遵循质量方程以及动量方程两个基本控制方程，本实施例中不考虑温度在装置中对房水流动的影响，故未涉及能量方程。在稳态不可压缩流动情况下，质量守恒方程的形式如式(1)所示:
式中，u,v,w是速度矢量在x,y,z方向的分量。该定律可表述为：单位时间内流体微元体中质量的增加，等于单位时间内流入该微元体的净质量。
纳维-斯托克斯方程(Navier-Stokes equation，N-S方程)是描述粘性不可压缩流体动量守恒的运动方程。它的矢量形式为(2)：
式中，ρ是流体密度，p是压力，常数μ是动力黏性系数。N-S方程概括了粘性不可压缩流体流动的普遍规律，因而在流体力学中具有特殊意义。
(2)计算模型及网格划分
图2-1展示了初始计算模型以及流道进出口位置，入口为引流管前端伸入前房的位置，高眼压时房水从引流管入口流入引流装置，房水从各个出口流出扩散到其它组织。鉴于结构的对称性，左右两侧对称出口命名采用同种数字，如1,1’。本具体应用实例中结果陈述时重点讨论单侧出口。图2-2为引流装置的剖面图。
图2-3所示为计算域网格分布，计算域的离散方法可分为两类：非结构化和结构化网格，考虑到模型结构的不规则性，采用非结构网格对计算区域进行划分，为保证网格质量对各个出口进行了局部加密，图2-3中右侧为出口2’和出口7的网格加密局部放大图。整个流场计算域的网格总数为291万左右。
(3)边界条件
本数值模拟的边界条件为压力出口和压力入口，出口压力参照已有的文献将其设置为9mmHg(参考静脉压和滤过阻力)。临床上认为高于21mmHg是病理性的高眼压，15mmHg为正常眼压，因为所设计的引流装置内有阀门，在一定的眼压下(<21mmHg)阀门打开，当房水压力降低到一定程度时(>＝15mmHg)阀门关闭，所以本实施例采用的入口压力取值在15～21mmHg之间。
本实施例的模拟采用Fluent 3D双精度求解器对流动区域的参数进行设置，房水用粘性流体来代替，其密度为990kg/m3，粘度为7.1ⅹ10-4Pa s。压力和速度采用压力速度耦合的SIMPLE算法，此算法使用压力和速度之间的相互校正关系来强制质量守恒并获取压力场。离散格式采用一般用于非结构网格的second-order upwind格式，精度较高。欠松弛因子采用默认值0.3。
所述的SIMPLE算法是一种压力修正法，通过“先猜想后修正”的方法得到压力场，并求解离散化的动量方程。具体计算过程如下：
连续性方程的离散化格式为：
Jf为通过面f的质量通量，ρvn，后面用面通量表示。
为了进一步研究，有必要将面的速度与单元中心存储的速度值相关联。单元中心速度对面的线性插值导致压力的非物理检测。速度的面值不是线性平均的；相反，使用加权因子系数αp使得动量加权平均。那么面通量Jf可以写成：
其中pc0，pc1，vn,c0，vn,c1分别为面两侧两个单元的压力和法向速度，在这些单元中包含了速度的影响，df是/>的函数，/>即面两边单元的动量方程系数ap的平均值。
SIMPLE算法用假设的压力p*来求解动量方程，得到面通量J* f，由式(4)得到：
不满足连续性方程。因此，对面通量J* f进行修正，修正因子为J'f，修正后的面通量为，
Jf＝J* f+J'f (6)
满足连续性方程。其中，
J'f＝df(p'c0-p'c1) (7)
p'为压力单元修正。
将通量修正方程(6)和(7)代入离散连续性方程(3)，得到单元内压力修正p'的离散方程:
其中源项b为进入单元的净流量:
压力校正方程(8)可用代数多重网格方法求解。一旦得到一个解决方案，单元压力和面通量的修正为：
p＝p*+αpp' (10)
Jf＝J* f+df(p'c0-p'c1) (11)
αp为压力欠松弛因子。在每次迭代过程中，修正后的面通量Jf完全满足离散连续性方程。
(4)输出参数
平均差是用来描述总体的变异程度的指标，指样本中各测量值与其算术平均数离差绝对值的算术平均数，用公式表示为：本实施例中xi表示为各出口流量比率，n为出口数量，/>为相同入口压力下各个出口流量比率的算术平均数。本实施例采用平均差来表示各出口流量分配的均匀程度。
下面结合附图，并以“沟槽状青光眼房水引流装置”为具体应用实例，对本发明上述实施例所提供的技术方案进一步详细的说明。
如图1所示，本发明实施例提供了一种沟槽状青光眼房水引流装置结构优化方法包括以下步骤：
步骤(1)确定目标对象，明确结构优化要求
步骤(2)创建或修正所述目标对象的3D几何模型，提取流体流通域；
步骤(3)开展流体流通域模型的网格划分以及有限体积数值计算；
步骤(4)分析流动特性，如满足设计要求，完成结构优化；分析流动特性，若不满足设计要求，重新提出结构设计优化方案，重新执行上述步骤(2)和步骤(3)，直至满足设计要求。
本发明实施例中：
目标对象为所述的一种沟槽状青光眼房水引流装置，如图2-1所示，包含引流管8和引流盘9，引流盘9为椭圆盘状结构，靠近引流管一端较厚，远离引流管的一端较薄；引流盘上表面设有引流沟槽或孔道，从引流管近端向远端编号为1和1’(近沟槽)，2和2’(近孔道1)，3和3’(近孔道2)，4和4’(近孔道3)，5和5’(远沟槽)，6和6’(远孔道)，以及7(直沟槽)。
结构优化要求为各出口的流体流量分配基本均匀。
步骤(2)中，流体流通域包括引流管、沟槽及孔道，即房水过流区域。
步骤(3)中，采用流体力学领域的数值模拟软件的前处理模块对步骤(2)所述的流体流通域几何结构进行网格划分，基于所述几何模型的不规则性，采用非结构化网格对流通域进行空间离散。步骤(3)中，采用流体力学领域的数值模拟软件的流体分析模块实现所述有限体积数值计算，具体包括：对流体属性、流体流通域边界、模拟过程中的算法、精度和离散格式等参数进行设置，基于设置的参数，采用有限体积算法，获取房水引流装置的出入口流量等计算数据。
步骤(4)中，采用流体力学领域的数值模拟软件的后处理模块分析所述的计算数据，如满足设计要求，即完成结构优化，若不满足设计要求，提出进一步结构设计优化方案，返回步骤(2)和步骤(3)，重新开展仿真计算，直至满足设计要求。
步骤(4)中，采用沟槽流道或孔道出口与入口的流量比来表示所述房水引流装置的流体流量分配情况，选用平均差公式：(xi表示为各出口与入口的流量比，为相同入口压力下各出口流量比的算术平均数，n为出口数量)计算平均差，用于衡量各出口流量比偏离所述平均值的程度，定量表述沟槽流道或孔道出口的流体流量分配均匀度。
采用一定的入口压力，通过本发明所述的结构优化方法，经上述流体流通域范围，计算所得的沟槽流道或孔道出口的流体流量比及平均差的详细数据；根据上述数据，利用数值模拟软件计算沟槽流道或孔道出口的流体流量分配均匀度，综合评价结构设计优化方案的合理性和有效性，进一步提出可能的优化方案。
下面结合具体应用实例，对本发明上述实施例所提供的技术方案进一步描述如下。
具体应用实例1
本实例主要以基本款房水引流装置为优化对象，通过本发明上述实施例所提供的计算方法，计算房水引流装置的压力场和沟槽流道或孔道出口的流体流量分配均匀度，综合评价基本款房水引流装置结构的合理性和有效性。
本实例中基本款房水引流装置，如图2-1所示，包含引流管8和引流盘9，引流盘9为椭圆盘状结构，长轴18.0mm，短轴13.0mm；靠近引流管一端比较厚，厚度为1.38mm，用于插入引流管；远离引流管的一端比较薄，厚度为0.43mm；引流盘上表面因各部分厚度不同呈现出一定的弧度形状，下表面为平面。按照设计，沟槽宽1.2mm；最大深度为0.96mm，最大槽深位于沟槽1和1’位置；最小深度为0.03mm，位于引流盘沟槽7位置。底面圆形引流出口2、3、4的半径为0.6mm。具体流体流道如图2-4所示。
各出口面积为：出口1(1.17mm2)、出口2(1.04mm2)、出口3(1.04mm2)、出口4(1.04mm2)、出口5(0.38mm2)、出口6(1.13mm2)和出口7(0.04mm2)，如图2-5a所示。
采用15、16、17、18、19、20和21mmHg的入口压力，通过本发明上述实施例所提供的结构优化方法，经上述流体流通域范围，计算所得的沟槽流道或孔道出口的流体流量比及平均差的详细数据，如图2-5b，2-5c和2-5d所示。
从图2-5结果可见，具有对称结构的1和1’(近沟槽)，2和2’(近孔道1)，3和3’(近孔道2)，4和4’(近孔道3)，5和5’(远沟槽)，6和6’(远孔道)，分别具有基本相同的流体百分数，对于基本款房水引流装置，基于各沟槽或孔道出口流量比计算的平均差比较高，说明流场均匀性很不好，尤其以7号口最小、1号出口较小和6号最大。
具体应用实例2
针对具体应用实例1中出现的问题，在具体实例2中，我们通过调整引流盘曲率来使房水引流装置更贴合眼球曲率，同时调整沟槽或孔道出口大小，进而调节流场均匀度。
本实例以第一种改进型房水引流装置为优化对象，通过本发明上述实施例所提供的计算方法，计算房水引流装置的压力场和沟槽流道或孔道出口的流体流量分配均匀度，综合评价第一种改进型房水引流装置结构的合理性和有效性。
第一种改进型房水引流装置包含引流管8和引流盘9，引流盘9为椭圆盘状结构，长轴18.0mm，短轴13.0mm；靠近引流管一端比较厚，厚度为1.34mm，用于插入引流管；远离引流管的一端比较薄，厚度为0.43mm；引流盘短轴方向下表面弧度为13.0mm，上表面弧度为18.0mm，长轴方向上下表面弧度分别为18.59mm和18.23mm；沟槽深度为引流盘厚度的1/2。按照设计，沟槽宽1.2mm；最大槽深位于沟槽1位置，最大深度为0.51mm；最小沟槽深度位于引流盘沟槽5位置，最小深度为0.16mm。底面圆形引流出口2、3、4的半径为0.6mm。具体流道如图3-1所示。
各出口面积为：出口1(0.59mm2)、出口2(0.8mm2)、出口3(0.76mm2)、出口4(0.82mm2)、出口5(0.21mm2)、出口6(0.79mm2)和出口7(0.26mm2)，如图3-2a所示。
采用15、16、17、18、19、20和21mmHg的入口压力，通过本发明上述实施例所提供的结构优化方法，经上述流体流通域范围，计算所得的沟槽流道或孔道出口的流体流量比及平均差的详细数据，如图3-2b，3-2c和3-2d所示；
优化结果表明，出口1和7的流量比有所提高，流体流量平均差为8.63-8.75％，但是出口2和4却严重减少，可能是因为出口3影响了流体的流向，因此在后续优化中可以考虑去掉出口3。
具体应用实例3
针对具体应用实施2中出现的问题，在本实例中，通过去除引流盘出口3，形成第二种改进型房水引流装置，来使房水引流装置流体分布更均匀，进而调节流场均匀度。
本实例主要以第二种改进型房水引流装置为优化对象，通过本发明上述实施例所提供的计算方法，计算房水引流装置的压力场和沟槽流道或孔道出口的流体流量分配均匀度，综合评价第二种改进型房水引流装置结构的合理性和有效性。
第二种改进型房水引流装置包含引流管8和引流盘9，引流盘9为椭圆盘状结构，长轴18.0mm，短轴13.0mm；靠近引流管一端比较厚，厚度为1.34mm，用于插入引流管；远离引流管的一端比较薄，厚度为0.43mm；引流盘下表面弧度为13mm，上表面弧度为18mm；沟槽深度为引流盘厚度的1/2。按照设计，沟槽宽1.2mm；最大槽深位于沟槽1位置，最大深度为0.51mm；最小沟槽深度位于引流盘沟槽5位置，最小深度为0.16mm。底面圆形引流出口2、3、4的半径为0.6mm。具体流道如图4-1所示。
各出口面积为：出口1(0.59mm2)、出口2(0.8mm2)、出口4(0.82mm2)、出口5(0.21mm2)、出口6(0.79mm2)和出口7(0.26mm2)，见图4-2a。
采用15、16、17、18、19、20和21mmHg的入口压力，通过本发明上述实施例所提供的结构优化方法，经上述流体流通域范围，计算所得的沟槽流道或孔道出口的流体流量比及平均差，部分数据如图4-2b和4-2d所示；流体流量比平均差，如图4-2c所示。
优化结果表明，各出口流体流量均匀度有所提高，流体流量平均差缩小为7.86-8.05％；但是1和1’(近沟槽)的流量比太高，有必要调整沟槽底面，以减少出口截面积。
具体应用实施4
针对具体应用实施3中出现的问题，在本实例中，通过构建统一的沟槽底面曲面利于样品加工，并且调整各出口截面积，形成第三种改进型房水引流装置，来使房水引流装置流体分布更均匀，进而调节流场均匀度。
本实例主要以第三种改进型房水引流装置为优化对象，通过本发明上述实施例所提供的计算方法，计算房水引流装置的压力场和沟槽流道或孔道出口的流体流量分配均匀度，综合评价第二种改进型房水引流装置结构的合理性和有效性。
第三种改进型房水引流装置包含引流管8和引流盘9，引流盘9为椭圆盘状结构，长轴18.0mm，短轴13.0mm；引流盘短轴方向上表面曲率为12.5mm，下表面与凹槽的曲率均为12.92mm,长轴方向上下表面曲率分别为26.04mm和24.92mm。按照设计，沟槽宽1.2mm；最大槽深位于沟槽前端靠近主流道位置(黑色框标记位置)，最大深度为1.05mm，沟槽1位置的深度为0.7mm；最小沟槽深度位于引流盘沟槽7位置，最小深度为0.15mm；槽底距离下表面最小距离是0.26mm；引流通道距离上表面最小距离为0.29mm，距离下表面最小距离为0.23mm，引流通道的内径为0.8mm；靠近引流管一端比较厚，厚度为1.49mm，用于插入引流管；远离引流管的一端比较薄，厚度为0.55mm；。底面圆形引流出口2、3、4的半径为0.6mm。具体流道如图5-1所示。
各出口面积为：出口1(0.87mm2)、出口2(0.52mm2)、出口4(0.53mm2)、出口5(0.8mm2)、出口6(0.55mm2)和出口7(0.16mm2)，见图5-2a。
采用15、16、17、18、19或24mmHg的入口压力，通过本发明上述实施例所提供的结构优化方法，经上述流体流通域范围，计算所得的沟槽流道或孔道出口的流体流量比及平均差，部分数据如图5-2b和5-2d所示；流体流量比平均值，如图5-2c所示。
优化结果表明，各出口流体流量均匀度有所提高，流体流量平均差缩小为2.2-3.4％，此优化模型已经达到良好的流量分配均匀效果，但为了与实际加工相适应，有必要进一步调整引流盘厚度与曲率，以达到流体分配均匀及结构贴合实际情况的目的。
具体应用实施5
针对具体应用实施4中出现的问题，在本实例中，通过优化沟槽底面曲面曲率半径，形成第四种改进型房水引流装置，来使房水引流装置流体分布更均匀，进而调节流场均匀度。
本实例主要以第四种改进型房水引流装置为优化对象，通过本发明上述实施例所提供的计算方法，计算房水引流装置的沟槽流道或孔道出口的流体流量分配比，计算流量平均差，综合评价第四种改进型房水引流装置结构的合理性和有效性。
第四种改进型房水引流装置包含引流管8和引流盘9，引流盘9为椭圆盘状结构，长轴18.0mm，短轴13.0mm；引流盘曲率包括下表面曲率、引流槽底面曲面和引流盘上表面曲率，椭圆短轴方向的引流盘下表面曲率半径为17.0mm、引流槽底面曲面曲率半径为17.0mm、引流盘上表面曲率12.5mm，引流盘椭圆长轴方向的下表面和上表面的曲率半径分别为34.47mm和32.93mm。按照设计，沟槽宽1.2mm；最大槽深位于沟槽1位置，最大深度为1.44mm；最小沟槽深度位于引流盘沟槽7位置，最小深度为0.30mm；槽底距离下表面最小距离是0.25mm；引流通道距离上表面最小距离为0.25mm，距离下表面最小距离为0.45mm，引流通道的内径为0.8mm；底面圆形引流出口2和4的半径为0.6mm。具体流道如图6-3所示。
各出口面积为：出口1(1.79mm2)、出口2(0.8mm2)、出口4(0.81mm2)、出口5(1.12mm2)、出口6(0.83mm2)和出口7(0.36mm2)，见图6-4a。
采用15、16、17、18、19、20或21mmHg的入口压力，通过本发明上述实施例所提供的结构优化方法，经上述流体流通域范围，计算所得的沟槽流道或孔道出口的流体流量比及平均差，部分数据如图6-4b和6-4d所示；流体流量比平均值，如图6-4c所示。
优化结果表明，各出口流体流量均匀度有明显提高，改进型沟槽状青光眼房水引流装置的流线分布图如图6-5所示，流体流量平均差进一步缩小为2.1-3.0％，基本符合实际应用。
根本发明另一实施例提供了一种房水引流装置的结构优化系统，如图7所示，可以包括：目标获取模块、流体流通域模型模块、模型参数设置模块以及结构优化分析模块。
其中：
目标获取模块，该模块用于确定目标对象，明确结构优化要求；
流体流通域模型模块，该模块用于创建或修正所述目标对象的3D几何模型，建立流体流通域模型；
模型参数设置模块，该模块用于对流体流通域模型进行网格划分以及有限体积数值计算；
结构优化分析模块，该模块用于分析流体流通域模型的流动特性，当流动特性满足结构优化要求时，完成对目标对象的结构优化。
本发明第三个实施例提供了一种房水引流装置，该装置采用本发明上述实施例中任一项结构优化方法进行结构优化。
本发明第四个实施例提供了一种终端，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，所述处理器执行所述程序时可用于执行本发明上述实施例中任一项所述的方法。
可选地，存储器，用于存储程序；存储器，可以包括易失性存储器(英文：volatilememory)，例如随机存取存储器(英文：random-access memory，缩写：RAM)，如静态随机存取存储器(英文：static random-access memory，缩写：SRAM)，双倍数据率同步动态随机存取存储器(英文：Double Data Rate Synchronous Dynamic Random Access Memory，缩写：DDR SDRAM)等；存储器也可以包括非易失性存储器(英文：non-volatile memory)，例如快闪存储器(英文：flash memory)。存储器用于存储计算机程序(如实现上述方法的应用程序、功能模块等)、计算机指令等，上述的计算机程序、计算机指令等可以分区存储在一个或多个存储器中。并且上述的计算机程序、计算机指令、数据等可以被处理器调用。
上述的计算机程序、计算机指令等可以分区存储在一个或多个存储器中。并且上述的计算机程序、计算机指令、数据等可以被处理器调用。
处理器，用于执行存储器存储的计算机程序，以实现上述实施例涉及的方法中的各个步骤。具体可以参见前面方法实施例中的相关描述。
处理器和存储器可以是独立结构，也可以是集成在一起的集成结构。当处理器和存储器是独立结构时，存储器、处理器可以通过总线耦合连接。
本发明第五个实施例提供了一种计算机可读存储介质，其上存储有计算机程序，该程序被处理器执行时可用于执行本发明上述实施例中任一项所述的方法。。
本发明上述实施例提供了一种针对房水引流装置结构优化设计的房水引流装置及其结构优化方法、系统、终端及介质，包括确定目标对象，明确优化设计要求，创建3D几何模型，提取流体流通域，开展流通域三维网格划分，进行有限体积数值计算，分析流动特性。如满足设计要求，完成结构设计优化，若不满足，提出结构设计优化方案，修正3D几何模型，重新开展仿真计算，直至满足设计要求。该方法能够实现一种房水引流装置的结构设计优化，使各出口房水流量分配更加均匀，减少房水引流装置对生物机体微环境的扰动或刺激，延长装置使用寿命。
本发明上述实施例提供的房水引流装置及其结构优化方法、系统、终端及介质，针对房水引流装置存在的房水流量分配不均匀的问题，提供了一种房水引流装置及其结构优化方法，基于CFD技术的支持，对房水引流装置几何结构的合理性和有效性进行数值计算，并提出可能的优化改进意见，最终实现沟槽流道或孔道各出口流体流量分配基本均匀，减少青光眼房水引流装置植入术对术区微环境的扰动或刺激，提高青光眼房水引流装置植入术手术成功率。
本发明上述实施例所提供的方法，不限于沟槽状青光眼房水引流装置，也对其它房水引流装置的结构优化提供了一种解决思路。
需要说明的是，本发明提供的方法中的步骤，可以利用系统中对应的模块、装置、单元等予以实现，本领域技术人员可以参照方法的技术方案实现系统的组成，即，方法中的实施例可理解为构建系统的优选例，在此不予赘述。
本领域技术人员知道，除了以纯计算机可读程序代码方式实现本发明提供的系统及其各个装置以外，完全可以通过将方法步骤进行逻辑编程来使得本发明提供的系统及其各个装置以逻辑门、开关、专用集成电路、可编程逻辑控制器以及嵌入式微控制器等的形式来实现相同功能。所以，本发明提供的系统及其各项装置可以被认为是一种硬件部件，而对其内包括的用于实现各种功能的装置也可以视为硬件部件内的结构；也可以将用于实现各种功能的装置视为既可以是实现方法的软件模块又可以是硬件部件内的结构。
以上对本发明的具体实施例进行了描述。需要理解的是，本发明并不局限于上述特定实施方式，本领域技术人员可以在权利要求的范围内做出各种变形或修改，这并不影响本发明的实质内容。
Claims (11)
1.一种房水引流装置的结构优化方法，其特征在于，包括：
S1，确定目标对象，明确结构优化要求；
S2，创建或修正所述目标对象的3D几何模型，建立流体流通域模型；
S3，对流体流通域模型进行网格划分以及有限体积数值计算；
S4，分析流体流通域模型的流动特性，如满足结构优化要求，完成结构优化；如不满足结构优化要求，重新执行上述S2和S3，直至满足结构优化要求；
所述S1中，目标对象包括：引流管以及具有弧度的引流盘，所述引流盘表面具有放射状沟槽流道和/或孔道；
所述结构优化要求为：所述引流盘表面的各沟槽流道和/或孔道出口的流体流量分配基本均匀；
所述S3中，对流体流通域模型进行网格划分的方法，包括：采用计算流体动力学的数值模拟软件的前处理模块对所述流体流通域模型进行网格划分，基于所述流体流通域模型的不规则性，采用非结构化网格对所述流体流通域模型进行空间离散；
所述S3中，对流体流通域模型进行有限体积数值计算的方法，包括：采用计算流体动力学的数值模拟软件的流体计算模块对流体属性、流体流通域边界、模拟过程中的算法、精度和网格划分后的空间离散格式进行参数设置，并根据设置的参数，采用有限体积数值计算的方法获取房水引流装置的出入口流量数据。
2.根据权利要求1所述的一种房水引流装置的结构优化方法，其特征在于，所述S2中，创建或修正所述目标对象的3D几何模型，包括如下任一项或任意多项：
-引流盘长轴和/或短轴的尺寸；
-引流盘厚度和/或上下表面曲率半径；
-沟槽流道的宽度、深度和/或底曲面曲率半径；
-孔道的直径和/或位置；
-引流管的管壁厚度和/或内径。
3.根据权利要求1所述的一种房水引流装置的结构优化方法，其特征在于，所述S2中，建立流体流通域模型的方法，包括：
基于房水过流区域确定流体流通域的几何结构，采用三维建模方法创建或重建流体流通域模型。
4.根据权利要求1所述的一种房水引流装置的结构优化方法，其特征在于，所述对流体流通域模型进行有限体积数值计算的方法中，还包括如下任意一项或任意多项：
-所述流体属性包括：流体密度和/或粘度；
-所述流体流通域边界包括：压力入口和/或压力出口；
-所述模拟过程中采用SIMPLE算法；
-所述网格划分后的空间离散格式采用second-orderupwind格式；
-所述有限体积数值计算采用非耦合隐式算法。
5.根据权利要求1所述的一种房水引流装置的结构优化方法，其特征在于，所述S4中，采用计算流体动力学的数值模拟软件的后处理模块分析流体流通域模型的流动特性。
6.根据权利要求1-5任一项所述的一种房水引流装置的结构优化方法，其特征在于，采用沟槽流道或孔道出口与入口的流量比来表示房水引流装置的流体流量分配情况，选用平均差公式：计算平均差，用于衡量各出口流量比偏离所述平均差的程度，定量表述沟槽流道或孔道出口的流体流量分配均匀度；其中，xi表示为各出口与入口的流量比，为相同入口压力下各出口流量比的算术平均数，n为出口数量。
7.根据权利要求6所述的一种房水引流装置的结构优化方法，其特征在于，基于所述沟槽流道或孔道出口的流体流量分配均匀度，综合评价所述结构优化方法的合理性和有效性。
8.一种房水引流装置的结构优化系统，其特征在于，包括：
目标获取模块，该模块用于确定目标对象，明确结构优化要求；
流体流通域模型模块，该模块用于创建或修正所述目标对象的3D几何模型，建立流体流通域模型；
模型参数设置模块，该模块用于对流体流通域模型进行网格划分以及有限体积数值计算；
结构优化分析模块，该模块用于分析流体流通域模型的流动特性，当流动特性满足结构优化要求时，完成对目标对象的结构优化；
所述目标对象包括：引流管以及具有弧度的引流盘，所述引流盘表面具有放射状沟槽流道和/或孔道；
所述结构优化要求为：所述引流盘表面的各沟槽流道和/或孔道出口的流体流量分配基本均匀；
对流体流通域模型进行网格划分的方法，包括：采用计算流体动力学的数值模拟软件的前处理模块对所述流体流通域模型进行网格划分，基于所述流体流通域模型的不规则性，采用非结构化网格对所述流体流通域模型进行空间离散；
对流体流通域模型进行有限体积数值计算的方法，包括：采用计算流体动力学的数值模拟软件的流体计算模块对流体属性、流体流通域边界、模拟过程中的算法、精度和网格划分后的空间离散格式进行参数设置，并根据设置的参数，采用有限体积数值计算的方法获取房水引流装置的出入口流量数据。
9.一种房水引流装置，其特征在于，采用权利要求1-7中任一项所述的结构优化方法进行结构优化。
10.一种终端，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，其特征在于，所述处理器执行所述程序时可用于执行权利要求1-7中任一项所述的方法。
11.一种计算机可读存储介质，其上存储有计算机程序，其特征在于，该程序被处理器执行时可用于执行权利要求1-7中任一项所述的方法。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011267286.1A CN112417777B (zh) | 2020-11-13 | 2020-11-13 | 房水引流装置及其结构优化方法、系统、终端及介质 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011267286.1A CN112417777B (zh) | 2020-11-13 | 2020-11-13 | 房水引流装置及其结构优化方法、系统、终端及介质 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN112417777A CN112417777A (zh) | 2021-02-26 |
| CN112417777B true CN112417777B (zh) | 2024-01-09 |
Family
ID=74832365
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202011267286.1A Active CN112417777B (zh) | 2020-11-13 | 2020-11-13 | 房水引流装置及其结构优化方法、系统、终端及介质 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112417777B (zh) |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2001022362A1 (en) * | 1999-09-20 | 2001-03-29 | The Board Of Trustees Of The University Of Illinois | Circulation model and applications |
| CA2795835A1 (en) * | 2010-04-30 | 2011-11-03 | Exxonmobil Upstream Research Company | Method and system for finite volume simulation of flow |
| WO2017026995A1 (en) * | 2015-08-07 | 2017-02-16 | Halliburton Energy Services, Inc. | Modelling of fluid elasticity effects on proppant dynamics |
| WO2017055252A1 (de) * | 2015-09-28 | 2017-04-06 | Carl Zeiss Meditec Ag | Vorrichtung und verfahren zur fixierung einer relativen geometrischen lage eines auges zu einem ophthalmologischen diagnose- und/oder therapiesystem |
| CN209967025U (zh) * | 2018-12-29 | 2020-01-21 | 中山大学中山眼科中心 | 一种眼部化学伤结膜囊冲洗器 |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US8157742B2 (en) * | 2010-08-12 | 2012-04-17 | Heartflow, Inc. | Method and system for patient-specific modeling of blood flow |
-
2020
- 2020-11-13 CN CN202011267286.1A patent/CN112417777B/zh active Active
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2001022362A1 (en) * | 1999-09-20 | 2001-03-29 | The Board Of Trustees Of The University Of Illinois | Circulation model and applications |
| CA2795835A1 (en) * | 2010-04-30 | 2011-11-03 | Exxonmobil Upstream Research Company | Method and system for finite volume simulation of flow |
| WO2011136861A1 (en) * | 2010-04-30 | 2011-11-03 | Exxonmobil Upstream Research Company | Method and system for finite volume simulation of flow |
| WO2017026995A1 (en) * | 2015-08-07 | 2017-02-16 | Halliburton Energy Services, Inc. | Modelling of fluid elasticity effects on proppant dynamics |
| WO2017055252A1 (de) * | 2015-09-28 | 2017-04-06 | Carl Zeiss Meditec Ag | Vorrichtung und verfahren zur fixierung einer relativen geometrischen lage eines auges zu einem ophthalmologischen diagnose- und/oder therapiesystem |
| CN209967025U (zh) * | 2018-12-29 | 2020-01-21 | 中山大学中山眼科中心 | 一种眼部化学伤结膜囊冲洗器 |
Non-Patent Citations (1)
| Title |
|---|
| 不同自动化场景对管制员眼动行为影响研究;林思远;胡荣锦;王艳军;;人类工效学(第02期);全文 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN112417777A (zh) | 2021-02-26 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| EP2942006B1 (en) | Method and system for non-invasive computation of hemodynamic indices for coronary artery stenosis | |
| Sánchez et al. | Numerical model of aqueous humor drainage: effects of collector channel position | |
| CN112417777B (zh) | 房水引流装置及其结构优化方法、系统、终端及介质 | |
| US20200013514A1 (en) | A system and method of artificial intelligence and tomography imaging of human cornea | |
| Philip et al. | Hemodynamic simulation of abdominal aortic aneurysm on idealised models: Investigation of stress parameters during disease progression | |
| Sánchez et al. | Numerical simulation of the aqueous humor flow in the eye drainage system; a healthy and pathological condition comparison. | |
| Fang et al. | Stroke risk evaluation for patients with atrial fibrillation: Insights from left atrial appendage with fluid-structure interaction analysis | |
| Tang et al. | Geometric model and numerical study of aqueous humor hydrodynamics in the human eye | |
| He et al. | Finite element analysis on fluid filtration in system of permeable curved capillary and tissue | |
| Zhang et al. | Fluid-structure interaction simulation of aqueous outflow system in response to juxtacanalicular meshwork permeability changes with a two-way coupled method | |
| Sacco et al. | The role of conventional and unconventional adaptive routes in lowering of intraocular pressure: Theoretical model and simulation | |
| Song et al. | Numerical simulation of multi-field coupling in aqueous humor under the condition of dynamic pressure | |
| US20190209380A1 (en) | Apparatus for Changing the Refractive Power of the Cornea | |
| Basson et al. | An aqueous humour fluid dynamic study for normal and glaucomatous eye conditions | |
| Vilalta-Alonso et al. | Correlation between Hemodynamic Stresses and Morphometric Indices as a Predictor Potential of Abdominal Aortic Aneurysm Rupture | |
| Meng et al. | Three-dimensional choroid neovascularization growth prediction from longitudinal retinal OCT images based on a hybrid model | |
| Tomaszewski et al. | Analysis of artery blood flow before and after angioplasty | |
| CN117010299B (zh) | 基于血流动力学耦合模型的脑组织血流情况预测系统 | |
| US11369381B2 (en) | Tailored venous anastomosis for arteriovenous grafts | |
| 石垣秀記 et al. | Experimental study on blood coagulation inhibitory effect of venous air-trap chamber of hemodialysis circuit | |
| Gandhi | Modeling erythrocyte-released ATP in mouse aorta, and development of CAD models of vascular networks from experimental images | |
| Causin et al. | An embedded surface method for valve simulation. application to stenotic aortic valve estimation | |
| Köppl et al. | Numerical modelling of a peripheral arterial stenosis using dimensionally reduced models and machine learning techniques | |
| Nagy et al. | Methods, Validation and Clinical Implementation of a Simulation Method of Cerebral Aneurysms | |
| Huaman et al. | Coupling 1D blood circulation model and substance absorption model to study drug metabolization |
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 | ||
| TA01 | Transfer of patent application right | ||
| TA01 | Transfer of patent application right |
Effective date of registration: 20230823 Address after: 201600 No.5 Dongjing Road, Songjiang Industrial Zone, Songjiang District, Shanghai Applicant after: SHANGHAI HAOHAI BIOLOGICAL TECHNOLOGY Co.,Ltd. Address before: 200031 No. 83, Xuhui District, Shanghai, Fenyang Road Applicant before: EYE, EAR, NOSE, AND THROAT HOSPITAL OF FUDAN University |
|
| GR01 | Patent grant | ||
| GR01 | Patent grant |