- 积分
- 1
- 注册时间
- 2009-8-15
- 仿真币
-
- 最后登录
- 1970-1-1
|
本帖最后由 shero 于 2015-1-23 20:07 编辑
目的:用于对心血管进行有限元模型分析
关键词:ANSYS,二次开发,超弹性(Hyperelastic)
该二次开发的环境:
Win8.1 64位操作系统,Visual_Studio_Pro_2008_ChnSimp,Intel.Visual.Fortran.Compiler.Professional.v11.1.038, ANSYS14.5,安装过程参考师访老师的经验贴,在此特别感谢师访老师:http://forum.simwe.com/forum.php ... 1%E5%BC%80%E5%8F%91
Demiray本构形式:
采用ansys中子程序UserHyper.F进行二次开发,Fortran代码见附件
验证:对厚壁圆筒内壁面施加血压,对比血压与内径的变化曲线,如下图:
另外,给出一个略微复杂一点儿的各向异性本构模型:
该本构模型是Fung指数型,该模型不能够采用子程序UserHyper.F进行二次开发,需要用usermat.F才可以,但是由于涉及到复杂的理论,二次开发的代码编写就比较困难,我还没有编写出来,只能给一个参考资料,其中包含一个Neo-Hookean的usermat.F代码,并给了相应的理论讲解,具体推导我都已经做过,没有问题,现在就剩下使用代码计算不收敛的问题。
()(),
我想知道附件中关于Neo-Hookean的usermat.F二次开发的代码是否正确,我尝试做计算,不能收敛,如果这个代码最终可以做出来没问题,那么Fung的那个超弹性本构做ansys的usermat.F二次开发也是没问题的。希望大家看看“ANSYS_NeoHook_usermat”这个文档里面具体是哪里出了问题,还请大家多多指教,万分感谢!
Objective stress rates的解释(http://en.wikipedia.org/wiki/Objective_stress_rates)
Demiray和Fung本构模型以及验证数据出处:
Holzapfel, G. A., Gasser, T. C., & Ogden, R. W. (2000). A New Constitutive Framework for Arterial Wall Mechanics and a Comparative Study of Material Models. Journal of Elasticity and the Physical Science of Solids, 61(1-3), 1–48. doi:10.1023/A:1010835316564
希望小伙伴们可以一起交流探讨,可以通过以下信息联系我:
邮箱:shao_xiong_@126.com
qq:
|
本帖子中包含更多资源
您需要 登录 才可以下载或查看,没有账号?注册
×
评分
-
1
查看全部评分
-
|