作者shane ()
看板Cad_Cae
标题Re: [问题] 脑部组织形变
时间Wed Jan 6 16:22:37 2010
※ 引述《gegemini (遇见小野丽莎莎...)》之铭言：
: ※ 引述《shane ()》之铭言：
: : 先跟各位说声抱歉
: : 由於本人为医学生物背景出生
: : 很多东西都还不是很了解
: : 所以有些问题要请大家给个指导与方向
: : 专案背景:
: : 例:
: : 有位病人脑部有肿瘤，经由医学影像分析後，医师决定要开刀切除，
: : 由术前的影像决定动刀方向与位置，
: : 但当病人进到手术室打开脑壳後，
: : 因种种原因脑部开始形变，
: : 此时肿瘤位置也会与当初预想的有误差，
: : 故须利用外部资讯来on-line tracking脑部的形变。
: : 专案目标:
: : 想利用脑部表面资讯(利用摄影机追踪脑壳打开处头脑外露部分)
: : 或得整体(包含内部)脑部的形变量
: : 问题:
: : 1.3D mesh的建立是否只要输入头脑的轮廓即可?输入的格式(点座标)?
: 现在的科技有比较进步一点，你可以透过专业的绘图软体来绘制 头脑外型图，
: 如果你话的出来的话。你可以透过的应用软体可以有 pro/E SOLIDWORKS UG catia
: : 2.给予一个deform model跟mesh，是否可产出施力後的形变後mesh?
: ^^^^^^^^^^^^
: deform model 指的是可变形体吧，只要只要 mesh後的element具可变形的
: 特徵的话，当然你施力後，element会变形阿
: (有限元素法处理的就是element跟node的资讯)，所以理论上你可以看到形变後mesh
: : 3.有model跟形变前、後的mesh，是否可以解出施力方向、大小?
: 这各应该是可以吧，可能再实际实现上会有些困难
: : 4.符合需求又支援强大的软体有哪些?(ANSYS)?
: 前处理软体：hypermesh 或 patran
: 分析 ：Abaqus 或 marc 或 lsdyna
: 後处理软体：hypermesh 或 patran
: : 5.若之前完全没碰过CAE、CAD软体，自学符合效益吗?
: 不符合效益，不过如果你有这方面的天份当然就有效益拉
: 但是我还建议你找专业的厂商配合吧，可能会比较快有结果(但也只是可能)
: : 谢谢回答了
首先感谢各位前辈回答
先前的一些问题有了初步的解决方法
首先在建solid model方面
目前是藉着自行利用MATLAB开发的程式对3D脑部影像作segmentation
并且利用isosurface找出vertices,faces
将其转存成.txt档并且写成ANSYS的指令
接者将其复制贴入ANSYS的命令列
使ANSYS自动建立keypoint,line,area
最後在ANSYS内建立volume
当然在建立volume时会遇到一些问题:
像triangle的面积为零<=这个问题我把vertex稍作一下移动，已解决
还有facet intersection，想请问有没有适当的解决办法
用Hypermesh?
还是用ANSYS的workingplane分段做mesh然後glue?
再请问一下侦测facet intersection的cretiria是甚麽?
至於在mechanical model的使用上
我们目前将完全按照所找的论文输入参数值
虽然大脑并非简单的linear elastic、hyperelastic、hyperviscoelastic.........
但简单的线性弹性模型目前已足够我们使用
在形变後的各节点位移量上面
ANSYS内有内建的table可以输出
利用其他程式将其内差成一个矩阵形变场应该不成问题
其他比较有困惑的地方如下
1.ANSYS是否内建多种solver，而其会自动选择最适当的solver求解?
2.有办法自己写solver写进ANSYS中来求解吗?
3.有办法自己写一个新的mechanical model并自行定义参数?
4.查网路发现ANSYS有支援parallel processing，这是外插的组件，
还是某版本以上即有支援?
5.不用商业的软体，自己写程式=>model building、discretization、problem solver
的可能性有多大?
--
※ 发信站: 批踢踢实业坊(ptt.cc)
◆ From: 163.25.97.25
※ 编辑: shane 来自: 118.168.45.124 (01/08 00:11)
1F：→ Rexrex:内定就是sparse solver,可以改用PCG,自己写应该没他强吧 01/08 16:41
2F：→ Rexrex:MATLAB也可以求解,你可以参考看看 01/08 16:42
3F：→ Rexrex:软组织很难求解,加油罗~! 01/08 16:45