CadCaeCam.com是目前CAD/CAE/CAM/PLM类专业网站中,用户最多,技术含量最高的网站之一, 涵盖目前所有常用的C3P类软件技术讨论。

C3P专业门户

 找回密码
 注册

QQ登录

只需一步,快速开始

扫一扫,访问微社区

查看: 3187|回复: 0

生物力学有限元分析软件MedFea 1.0beta版

[复制链接]
发表于 2011-3-31 12:56:29 | 显示全部楼层 |阅读模式
生物力学有限元软件MedFeaV1.0beta_img_73.jpg % h9 W( h/ d5 ~
业务联系 HuangSteve@sohu.com$ F: X. w2 _) B% Q3 M5 }
生物力学有限元软件MedFeaV1.0beta_img_2.jpg   r0 M- I0 a' _! o
1         基于生物力学模型的有限元分析概述
! H/ l6 H5 I$ C! ?) v/ R  B随着有限元技术的发展,有限元在医学研究中的应用也越来越深入。由于生物组织材料具有不可压缩性、非线性、粘弹性、塑性与各向异性等力学性质,精确构造其力学模型是相当困难的,且计算复杂,难于满足实时性,因此针对生物组织的三维非线性有限元分析将是前瞻的与充满挑战性课题。
; U/ L' Z2 N8 _4 F
* W( ]. o' J2 n7 y  {0 f  p3 ?- U1 ~# B尽管有限元分析技术在工程中得到广泛应用,但在生物医学中应用尚处于初始阶段。目前几乎没有专用于生物组织的有限元力学分析开发包或软件,也没有成熟的有限元分析商业软件(如Ansys,abaqus等)提供专用于生物力学的有限元模型分析功能,因此研究生物力学模型与开发针对生物组织的有限元力学软件成为当务之急。, O. ?2 ~& {/ Q) x: w

9 c) U4 S- j9 Z! J经过近五年(2007~2011)奋战,本团队开发出专用于医学类生物力学分析的有限元软件MedFea(1.0beta版),其提供完整有限元分析功能,同时简化有限元分析步骤,用户能够快速学以致用。同时结合本团队所开发医学图像与生物模型处理软件MedSys的图像处理功能,便可实现一体化处理流程:
# u% U8 k5 _  ~, u0 c
& D! P2 X  d3 J6 R( {3 B图像处理(如分割,配准等)[MedSys] à 网格处理(如重构,简化与划分等)[MedFea] à 有限元前处理(如边界条件,载荷与分析步等)[MedFea] à 有限元分析[MedFea] à有限元后处理(如动画与录制等)[MedFea]。3 j0 s3 ?8 b  D2 t0 f- t
生物力学有限元软件MedFeaV1.0beta_img_9.jpg 1 V6 D. X, [7 m; a% v) k
% j9 G' M; B7 J3 @, L! F
此外,本团队开发的医学图像与生物模型处理软件处于后期开发与调试阶段,预期三个月后发布MedSys1.0beta版本,同样欢迎试用。在MedSys1.0后版本将会实现生理信号分析功能,目前生理信号分析算法也在改进中。% d( N) X6 M/ _) v
生物力学有限元软件MedFeaV1.0beta_img_11.jpg 生物力学有限元软件MedFeaV1.0beta_img_12.jpg
3 ^% W7 [2 U+ D' C2 ]: w. R% z2         生物力学有限元分析软件MedFea1.0功能
- [$ L" N1 N* D- A( ~2.1   MedFea1.0功能
8 l, I& {2 n: \" {, X4 zMedFea包括主页、图像、面网格、框网格、体网格、前处理、分析与后处理等模块,其中
/ z" U, L4 A* J- s( }  ~
6 z  X  p/ \7 r0 E' P: U9 U; F(1)  主页模块用于文件导入、视图设置与面板控制等;- I! d0 c$ W1 ]# L0 v  ~
8 t1 D" o5 A. o! ^3 L* x, Q' e
(2)  图像模块用于查看与管理图像,MedFea并没有图像处理功能(可使用医学图像与生物模型处理软件MedSys),其主要结合图像来进行材料赋值操作;
7 R( Q% R/ o% d8 W  Q1 E0 x- X) x, p+ j% n3 |3 w
(3)  面网格模块用于面网格相关处理,包括面重构与四面体网格划分等;
* a' b  q4 y3 h" w4 \
; V' K3 E# x. ]% @(4)  框网格模块用于六面体网格划分等;0 z. {( z( j/ [) M8 M6 ]6 E3 L
$ ^& e# ?* V" C$ ~! i
(5)  体网格模块用于评估网格质量、选择集合(单元集、结点集、面集与边集)与材料模型创建及其关联赋值等;0 R* P' ~# R' k/ C4 h5 H8 d

: D5 J, [3 P8 k0 ]  Z5 l+ q$ O+ E(6)  前处理模块用于边界条件、载荷、连接接触、刚性约束、初始条件与分析步等设置;
! y2 U4 L6 h' |$ W1 }) K* Z0 I& u: o
(7)  分析模块主要用于有限元分析;! S4 G4 ?3 G! K
. S6 C, n0 w* t& r' g
(8)  后处理模块主要用于查看、动画播放、录制与导出分析结果等;, O2 Y6 w! s4 v

/ {2 S# T) Z( V* P# E# m下面通过两大实例来描述软件处理流程与相关模块功能应用:2 T3 o8 u( Z  j+ n" A

& `+ i* I2 V. N. R9 S2 Q面网格处理(如重构,简化与四面体划分等)/框网格处理(如六面体划分等) à 体网格处理(如质量评估等) à 有限元前处理(如边界条件,载荷与分析步等) à 有限元分析 à 有限元后处理(如动画与录制等);+ Q* k1 @1 _5 c) K3 n4 x
! f  O3 i' m0 m4 o" P: ]) w
2.2   分析实例1
7 C& w; H0 t: i本实例主要分析骨头受力情况,实例中固定骨头一侧,而另一侧指定位置施加作用力。分析前先使用MedSys软件进行图像分割与生成面网格模型,并导出stl格式面网格。- B$ |. p& ?( X2 d6 n# U% E
4 K$ _% P* B6 x$ F7 L
1)       导入面网格; t" z$ y8 Q0 Z5 k
生物力学有限元软件MedFeaV1.0beta_img_24.jpg
( X. P! o3 Y6 M' H) w4 G! `2)       重构面网格。导入面网格后则可进行面网格优化处理。MedFea提供面重构技术可以有效提供面网格质量,其重构后面网格单元近似等边形。重构面网格也可以用于消减网格单元。
' A, a8 v: B! t7 }) z" e# O' ?# x  z  D$ d) b
MedFea面重构算法经过持续更进,其相当成熟与高效。下图为面网格(64772单元数与166607结点数)重构前后结果:# m. c% L: X& ]* u. u, r
生物力学有限元软件MedFeaV1.0beta_img_39.jpg 0 D4 G" \$ n6 k( S
3)       网格划分。对面网格进行网格划分,得到体网格。
' t- C, _6 [6 O# b 生物力学有限元软件MedFeaV1.0beta_img_40.jpg 2 o. B- ~# E1 p; K8 j8 D) i
4)       网格质量评估。网格划分后则可以通过质量评估来评估其体网格质量,高质量网格才能保证后续有限元分析完成。若网格质量不高,则可重新设置参数进行网格划分,或选择其它网格划分算法,直到满意为止。
7 M/ O3 {% p6 i9 p1 d5 T 生物力学有限元软件MedFeaV1.0beta_img_41.jpg 生物力学有限元软件MedFeaV1.0beta_img_42.jpg 生物力学有限元软件MedFeaV1.0beta_img_43.jpg 生物力学有限元软件MedFeaV1.0beta_img_44.jpg 生物力学有限元软件MedFeaV1.0beta_img_45.jpg 生物力学有限元软件MedFeaV1.0beta_img_46.jpg 生物力学有限元软件MedFeaV1.0beta_img_47.jpg 生物力学有限元软件MedFeaV1.0beta_img_48.jpg 生物力学有限元软件MedFeaV1.0beta_img_49.jpg
6 X9 M) j$ @4 h4 J% s, Y: ]5)       选取结点集。选取结点集nodeset1与nodeset2,其分别用于边界条件与载荷。
* o$ R1 Y; @7 i0 v" [" W0 J
; d, k# @6 P) s# D  s6)       材料模型。MedFea自动生成默认材料模型Initial_Materials与默认单元集Initial_Elements(其包括体网格中全部网格单元)。选择Initial_Materials与Initial_Elements,并点击菜单进行关联,即设置体网格中全部网格单元的材料模型为Initial_Materials。# N: G8 A+ D- t  D" ]% N! @: ~
, ?7 H3 G, {4 }( B- z0 p3 C, u
7)       添加边界条件。设置结点集nodeset1的位移x/y/z为固定位移,即设置其值为0,其意义该结点集中结点在分析过程中位移保持不变。
% |; p$ Y& G' }9 C1 m6 C  ]  g6 n) n+ f- o' t& `
8)       添加载荷。设置结点集nodeset2的z方向作用力为-100,其意义该结点集中结点受到作用力(0,0,-100),即作用力大小为100,方向为-z。
- w/ A. w/ V7 T( o& ]8 R) Z3 @) M$ u- ], n2 m8 }5 }- s* w
9)       添加分析步。MedFea默认创建全局分析步,但至少还需要创建一个计算分析步。其中设置线性算法为SuperLU。
% k5 P9 z, z+ ?. H
8 _/ H; ^3 n, [; H2 U  E; ~5 k10)    有限元分析
# h" }) D' r, H$ H
  o, T3 `/ S, U' v11)    后处理  
& v6 d9 J2 r3 b6 ^; B# I3 ]  [ 生物力学有限元软件MedFeaV1.0beta_img_50.jpg 生物力学有限元软件MedFeaV1.0beta_img_51.jpg
! o. N' N! T" [8 r& V  ! Y& t! N' a: A% \; x" ?
% `+ Z/ L& Q) ]

3 ?# W7 i, f, j
5 c: G. E( q$ M; B) Z/ a
4 S0 K& T5 z* R$ L7 g; {- t1 q8 V2 G! i0 o6 w
2.3   分析实例28 A/ g  z" O9 ?; Q% a, _5 K% ~
本实例主要模拟两根血管相交移动与形变过程,并分析其形变过程中各应力等。实例中固定一血管两端x/y/z方向,另一血管两端y/z方向,而x方向可自由移动,并设置其预期移动-6;同时设置连接接触面来避免两血管在移动过程出现交错现象,并设置接触面相关摩擦参数等。
) Y% I- B1 I8 h* u1 B. C
3 y8 [; e7 J0 V, n 生物力学有限元软件MedFeaV1.0beta_img_53.jpg 生物力学有限元软件MedFeaV1.0beta_img_54.jpg 生物力学有限元软件MedFeaV1.0beta_img_55.jpg 生物力学有限元软件MedFeaV1.0beta_img_56.jpg 生物力学有限元软件MedFeaV1.0beta_img_57.jpg 生物力学有限元软件MedFeaV1.0beta_img_58.jpg 生物力学有限元软件MedFeaV1.0beta_img_59.jpg 生物力学有限元软件MedFeaV1.0beta_img_60.jpg 生物力学有限元软件MedFeaV1.0beta_img_61.jpg 生物力学有限元软件MedFeaV1.0beta_img_62.jpg $ T% \) ^- S  S' w
6 @# g8 M0 j1 C5 f  N
1)       导入体网格,本实例中直接导入六面体网格,故不需要再进行网格划分等;, i) g" t( Q# e0 m/ ]8 N1 x2 u

4 X* x& L! X9 O7 y. D* k9 B
3 w: Z5 S3 P2 e6 ?  V' e3 }, B
* @# G7 w& e' V( e& W: {: n6 K& T5 |8 c6 }- y5 g5 W8 b/ {5 P
2)       选取结点集。分别选取两根血管的两端结点集cbknodeset1与cbknodeset2,其将用于边界条件中固定约束中;
9 R& P5 o* t, q) v
" o" d  e, e+ v# _! n) h+ {* S8 e6 I  P% G( |

+ X! \' {. H0 s- Q
& I8 o, @7 N4 z9 z% r% R* h3)       选取单元集。分别选取两根血管的单元集cbkelset1与cbkelset2,其将用于材料模型赋值,即每根血管分别使用不同材料模型;
! A/ d! R& X% H! a6 _: N, M0 p
+ y, _- S- f5 c' H
! v$ a* G: p; M- Z4 Y

$ [# Z# m% p3 W3 ]! t5 S, b+ B4)       选取面集。分别选取两根血管相互接触部分的面集cbkfaceset1与cbkfaceset2,其将用于连接接触处理;
; A: i% k  a  S8 A
7 G' }1 P& i; f* }1 E0 v- _7 S5)       创建材料模型。分别创建两类材料模型cbknh1与cbknh2;点击面板中按钮控件,并选择材料模型类型与输入材料模型名称则可创建材料模型对象;其中材料模型ID值不需要用户输入,也不能修改,软件自动分配与维护; 然后设置材料模型参数,直接修改属性窗体控件中控件值即可生效;最后设置单元集的材料模型,即并联单元集与材料模型,选择单元集与其对应材料模型则可赋值;! O, s* K3 j+ t8 V2 z& u/ b
生物力学有限元软件MedFeaV1.0beta_img_63.jpg 7 O, y( U7 T2 r5 R8 {
生物力学有限元软件MedFeaV1.0beta_img_64.jpg
- K* `+ G( a+ N+ _4 U! E: I3 J3 @8 u 生物力学有限元软件MedFeaV1.0beta_img_65.jpg . }  V3 x5 m" x5 U
6)       添加边界条件。在MedFea有限元分析前需要设置相关边界条件、载荷与连接(接触)及分析步等。添加结点集cbknodeset1的x/y/x方向为固定约束,即均设置其值为0; 设置后点击按钮应用; 同时添加结点集cbknodeset2的y/x方向为固定约束,即均设置其值为0;而设置其x为规定约束,即设置-6值,其非0意义是该结点集当且仅当在x方向上移动;+ i2 J! N- U+ X0 e8 c# {

! a/ c& h! K/ N2 G6 X1 I' f
3 V5 [9 X9 I% m; G8 F5 G) a 生物力学有限元软件MedFeaV1.0beta_img_66.jpg # q$ v" ?' [; H. r1 p

* C& p1 B  {5 H* x) `/ w4 V! q  d  ]# a$ G; g5 e! ~
7)       添加连接接触。添加连接接触面来避免两血管在移动过程出现交错现象;点击连接面板中按钮则可创建连接接触对象,选择连接接触类型并设置接触相关参数,如摩擦系数与选择接触面等,其中接触面为前面所选取面单元集cbkfaceset1与cbkfaceset2。
+ s2 I; ]7 C' g6 T
4 Y- n9 M8 M6 X+ _6 K/ D' f  X 生物力学有限元软件MedFeaV1.0beta_img_67.jpg 7 W' U! I' I  i2 ]/ g( Y7 ?/ K! G& X

+ V1 J: R  z$ {" ]+ m" v2 Q+ L/ M7 T4 Z" Z0 f0 r
8)       添加分析步。创建与设置分析步,其中全局分析步(初始分析步)默认生成,并且用户至少还需要创建一个分析步。所有作用于全局分析步的边界条件、载荷与连接接触等均全局有效,否则仅指定分析步分析过程有效。
" i% ?. U/ t; ]+ `2 A 生物力学有限元软件MedFeaV1.0beta_img_68.jpg   |7 ^) e+ ?9 ^$ Q# q$ [& S3 _

3 [0 a$ \5 d* `% o" u& D9)       有限元分析。选择模型并点击分析计算即可,当然可以在分析前查看有限元文件,其使用vtk文件格式。+ x, h) r) R/ N  r' s8 B
生物力学有限元软件MedFeaV1.0beta_img_69.jpg
( |7 z" B7 t- C, {3 p10)    有限元后处理。
, a2 Z- h0 z2 y 生物力学有限元软件MedFeaV1.0beta_img_70.jpg / {8 B* S( {: S+ ?2 ]
生物力学有限元软件MedFeaV1.0beta_img_71.jpg ' ?$ ?: s- y% p" i% g, S

$ F! o( s4 C3 [3 ]$ \3 j# ` 生物力学有限元软件MedFeaV1.0beta_img_72.jpg
& d' h. Y2 B. L  ]& h$ A+ _0 n

- d8 W: C8 G# i2 t3 g: g* D
% d/ ~- G4 I3 p* q: E9 X* r- G8 b7 `7 R

/ _7 y0 ]# @) z# S- F% x, a
; Y$ h8 R. L) q& V( @! q  x

评分

参与人数 1C3P币 +50 收起 理由
bio + 50 赞一个!

查看全部评分

您需要登录后才可以回帖 登录 | 注册

本版积分规则

QQ|手机版|小黑屋|CadCaeCamPlm专业门户网站 ( 蜀ICP备05008798号 ) 网站速度测试

GMT+8, 2019-6-19 14:08 , Processed in 0.107021 second(s), 26 queries .

Powered by Discuz! X3.4 Licensed

© 2001-2017 Comsenz Inc.

快速回复 返回顶部 返回列表