血流动力学有限元后处理的集成化开发-北京工业大学学报.PDF

第35卷第5期 北京工业大学学报 Vb丑．35No．5 UNIVERSITY 2009年5月 JOURNALOFBEIJING OFTECHNOLOGY 2009 May 血流动力学有限元后处理的集成化开发 乔爱科1，陈为彩2，张宏斌1 (1．北京工业大学生命科学与生物工程学院。北京100124； 2．北京工业大学机械工程与应用电子技术学院。北京100124) 摘要：为了提高血流动力学有限元分析可视化后处理的效率，基于ANSYS、利用VC、调用MATLAB进行了集成 化软件开发．以ANSYS作为设计平台．借助VC设计参数化输入界面，调用MATLAB绘图功能实现了多种图形的 绘制．以ANSYS绘制的矢量图、等值线图。以MATLAB绘制的压力脉搏波时空分布图说明该程序完全可用．集成 化的后处理软件可以大大提高血流动力学可视化的效率． 关键词：应用软件；有限元后处理；血流动力学 中图分类号：Q66 文献标识码：A 文章编号：0254—0037(2009)05一0679一06 血流动力学数值模拟中会产生大量的结果信息．采用科学计算可视化的方法对这些信息进行处理是很 重要的[1-5】．科学计算可视化是将科学计算过程中所产生的数值转换成为计算机屏幕上直观几何图像的 过程，或者说是研究如何把科学数据转换成可视的、便于理解的信息的计算方法弦71．计算可视化主要针 对三维大规模数据场的可视化处理技术及其应用，核心是将三维数据转换为图像．它不但能加快数据处 理的速度，实现人与计算机之间的数据通讯，而且能实现数值模拟过程的直观化以及计算过程的交 互化[7|． ANSYS软件是目前数值模拟比较常用的工具，具有非常优良的后处理功能，可以显示彩色等值线图、 梯度图、矢量图、粒子流迹图、立体切片图、透明及半透明图等图形．但是，在血流动力学有限元分析方面。 ANSYS后处理的功能有限，尤其是对三维流场的壁面剪切应力时空分布、壁面剪切应力梯度分布、压力脉 搏波、速度脉搏波等图形的处理显得力不从心．与此相比，MATLAB显示了灵活的绘图功能[8-10】． APDL(ANSYS parametricdesign 施加参数化载荷，以及参数化后处理结果显示，从而实现参数化分析的全过程【1卜12】，但是使用APDL编程 不直观、不方便．VisualC++6．0是一种面向对象的程序设计语言，具有良好的图形用户界面，应用其开 发程序可以大大提高软件设计能力及开发速度．为了解决ANSYS分析中的问题，提高血流动力学有限元 分析可视化后处理的效率，作者基于ANSYS、利用VisualC++6．0和MATLAB，对血流动力学有限元后 处理进行了集成化开发． 1程序实现方法 ANSYS采用参数化后处理．VC提供输入参数化设计界面，由此界面可以输入或修改后处理的相关 参数。这些赋值后的参数被传递到ANSYS的参数化命令流文件中．运行ANSYS命令流文件得到的后处 理结果包括：速度矢量图、速度等值线图、速度剖面图、速度沿某一指定路径分布的曲线图、二次流、压力分 布、壁面剪切应力等图形，动画显示以及文本数据；在VC中利用这些数据调用MATLAB绘图，得到包括 壁面剪切应力时空分布、壁面剪切应力梯度分布、压力脉搏波和速度脉搏波时空分布．整个程序体系的基 本构成如图1所示． 收稿日期：2007—10．15． 基金项目：国家自然科学基金资助项目；北京市自然科学基金资助项目(3062003。3092004)． 作者简介：乔爱科(1967一)。男．山西平定人。副教授． 万方数据 北京工业大学学报 2009焦 启动程序 输入相关的量 自动形成APDL代码 堡塑型辈