血管流固耦合分析实例.docx

血管流固耦合分析实例.docx

《血管流固耦合分析实例.docx》由会员分享，可在线阅读，更多相关《血管流固耦合分析实例.docx（17页珍藏版）》请在冰点文库上搜索。

血管流固耦合分析实例

Ansys14workbench血管流固耦合实例

根据收集的一些资料，进行学习后，试着做了这个ansys14workbench的血管流

固耦合模拟，感觉能够耦合上，仅是熟悉流固耦合分析过程，不一定正确，仅供参考，希望大家多讨论。

谢谢！

1、先在proe5中建立血管与血液流体区的模型（两者装配起来），或者直接在workbench中建模。

1

图1模型图

2、新建工程。

在workbench中toolbox中选customsystem，双击FSI:

FluidFlow（fluent）->staticstructure.

图2计算工程

3、修改engineeringdata，因为系统缺省材料是钢，需要构建血管材料，如图

3所示。

先复制steel，而后修改密度1150kg/m3,杨氏模量4.5e8Pa，泊松比

0.3，重新命名，最后在主菜单中点击“updateproject”保存.

图3修改工程材料

4、模型导入，进入gemetry模块‘import外部模型文件。

图4模型导入图

5、进入FLUENT网格:

划分。

在workbench工程视图中的Mesh上点击右键，选择Edit…，如图5所示，进入网格划分meshing界面，如图6所示。

我们这里需要去掉血管部分，只保留血液几何。

詁■岂：

•*■*VepFrehi-t*丁□屮存片*・》，4r■>*

-

WZ21'■'1■■-■

■--■■X*/*X*X*

叫p厂ir■时⑥电m■#*s-jn*.sja

H1*V：

l.4iqpr"~4

Xi

生_T

一一UfiHH

■*：

+p^MnaijIm■'«

■hk*

旳J

™—1

A

图5进入网格划分

BlPE|ect

EJ

It

H^del【A3〕

■等IGedreirz

■-v*丑om回

-d®

7£LC-DFCSlatESrSte?

*.•礁】CcnnezlKria

fljlMrgll

>VSft

lUpddte

Phr-eview

PHid*Bad）gKdpAJO（h?

f

SuppirL号方出1Ol卜bBodies

ijpprrs^Boe

UpdateSelectedPans

图6禁用血管模型

6设置网格方法。

默认是采用ICEMCFD进行网格划分，设置万式如图7所示，截面圆弧边分为12份，纵截面的边均分为10份，网格结果如图8所示。

另外在这个界面中要设置边界的几何面，如inlet、outlet、symmetry

fij*>***€（-!

宜Mod*l（U）

l-lGmnclr,

用aijaop435t

x¥3LOCCftS5S£L[*：

J

ElC«fMnM*5rf»fnt

l+ll.：

»；■*

创Model（A3）

二_昭|Geanezif

!

ftOOMLW

'x.SLDODV^爼㈣E»’^4-匸gSHti呂祥54王…*^&匚DTILZltior1!

-赢Mesh

Ortifari£d

5<«pt

ScopingMethods

GeometiySetecticir

■Geometry

IEdge

Orftnrttsn

bupfXt-a&ed

NO

T^_r

Nm-benof

Nijnibr<-al05前5业

帖

Behavior

W輛

甘翻Ope

NoBim

Detaihof"FilprSizingJ

1-Siring

MopingMHhrd

C^Qrnrtiry7代和沖「

Gcgmrtiry

6Fdqr^

-Oefiriitrth

NC

丙pr

NurrberoilDwisiors

NunraWTOfDivi^iOre?

10

ScHSwior

Htrd

Bi&iType

ItaBhi

图7设置网格划分方式

图8最终出网格

图9边界几何

7、进入fluent

（Li）&e6rriStry

图10进入fluent

关闭mesh,回到fluent工程窗口，右键点击setup,选择edit…，进入fluent。

这里设置瞬态计算，流体为血液（密度1060,动力粘度0.004pas），入口压力

波动（用profile输入），出口压力0Pa,采用k-e湍流模型。

General

瞬态

血液材质

入口压力波动profile。

（（inletpresstransient70）

（time00.0040.0050.0060.010.0150.03）（pressure0150017001500000）

入口边界

图11Fluent中设置

计算时记录入口边界压力的变化，如图12所示，与profile设置近似

rijrtnoiM

Wl

rum

Int^gri

nron

JBJ

ADD

HDCiSO曲

i&ratm

鹿Ltl

ACDj

HUJ

!

*•ir■Is'i-.1Kiri"u■:

i

KNETEFLUErIT14□M.

Fe皿SQTJbkaIransidFif'

图12入口压力变化

面上压力分布（可以分步看）中间平面的压力与速度如图14所示。

流线如图15所示。

图13fluidinterface面上压力分布

dEO-L-HUrh.«-

O.SKns-asLB—^St8^tL»J.tofl-<«K£U-itv匕二$二

匕Hf»i-MDfs人AF^X

MPK3LE?

右2.a$

rlB-nnurEQii冒！

-富Fr・・£!

□□*■=

-■nl*hwu_TVX7

awin-Ei-s«87smg

■dT-fi

*豐I.-O1

*-NlagT-fi

」MOI

^Qss33盂：

c$盖snmui?

亠m,电liiw®aCTn*.=^

函」4甘「凹卡3（MWH）steHArrs;岡

Tucnai

AN.2YBFLUENT」4bEnu雷•^Us-fl4-

£3«-_rs5r0（X?

-

m!

£J$5

二二1=

>'Mu£3

『Z4tn*l

•S9t且

■nlM■»**&!

"3*»g

n3X^.

N;:

亠I

O8<*8

fffiK-?

K■--0rnm・・A』DO?

召3>

函15SB

is.72Miz食fnHFeNT:

Ds-ss-

8、关闭fluent,进入StaticStructure的model,如图16所示。

注意geometry树形节点上，禁用血液计算区域。

Outline

:

:

!

:

is

j……（9BLOOD1[39]

卑BLOODVESSEL[40]p鼻CcordnateSystems彳轴Connections”哪Mesh

.百StaticStructural

yGSettjnfls

白ImportedLoad（Sclutiori）电IrnptxtedP?

es^ure

自血|^clutwn（B6）

'_£|SotutwnInfermatKn

Detailsof'ImportedPressure**

Scope

^ccipingMethod

GeometrySrelection

Geometry

3Faces

-1

-1

Type

ImportedPrfisurf

弓uppr«M

No

TransferDefiniition

CFDSurface

inface二|

CFOData

CFDResultsFile

C:

\U5erj\.HP\DocumenU\blaodll#s\dpO\S¥S\MECH\S...

图16禁用流体区域

8、导入流体载荷

如图17所示，插入流体载荷，血管几何选择内圆面，CFDsurface选择预先设置的inface面，就是将inface的压力数据传递给血管的内圆面。

导入流体载荷后的模型如图17所示。

图17插入流体载荷

9、设置固体约束

约束血管一端为fixed，如图18所示。

进行计算。

计算结果如图18所示

图18血管一端约束

变形

10、总结

在ansys14workbench中进行流固耦合（FSI）计算相当简单，但这仅是单向耦合,如果实现双向耦合的话还是比较复杂的。

计算过程仅供参考，经不起专家的推敲啊，有建议发邮件给我啊，谢谢!