找回密码
 注册
Simdroid-非首页
查看: 630|回复: 20

[二次开发] 提问讨论:晶体塑性(多晶取向,黄永刚1991年单晶子程序)

[复制链接]
发表于 2018-8-29 15:07:37 | 显示全部楼层 |阅读模式 来自 宁波大学
本帖最后由 燃烧的向日葵 于 2018-8-29 15:24 编辑

提问交流贴:感觉关于晶体塑性的一些精华帖都是以前的老帖,近些年好的帖子非常少,今年4月份开始接触晶体塑性的课题,遇到一些难题总感觉无计可施,这里提出一些问题希望同行的仿友能多指导交流,不喜轻喷!:
  1,关于多晶取向的问题以前有版主讨论过,结合一些毕业论文感觉主要有:
     a:Matlab中建立晶粒取向数据,或者是利用EBSD获得数据,储存在文本文档中,再利用Python脚本读取赋予晶粒
(但是有的仿友提出赋与取向运行后,和没有设置取向得到的结果一样,感觉应该是晶粒取向和子程序UMAT没有结合起来,另外做晶体塑性涉及到织构演变的问题,后处理要提取取向演变后数据做极图或ODF?)
    b:Shawn版主之前回答的提问(提问截图见下图),感觉比较靠谱,但是具体实施还有待进一步研究UMAT子程序。


  2,关于将黄永刚单晶子程序修改为多晶的问题
   看了黄老师的单晶子程序和对应的论文将近四五个月,将迭代公式推导了一遍,才大概整体看懂了点,但是理论和编程基础薄弱,对于将子程序修改为适应于多晶二维Voronoi(目前做二维Voronoi晶体塑性),还是感觉困难很大。看有些文章里面说单晶过渡到多晶,要用到用到自洽模型,但是关于自洽模型的理论公式如何在子程序中体现?感觉目前需要做的是将理论书再好好过一遍(连续介质力学、弹塑性力学等),子程序编程先从简单的开始写,熟练了再修改黄老师的子程序。
希望做过晶体塑性的仿友能不吝赐教!后面学习的心得也会在这个帖子里跟新。

本帖子中包含更多资源

您需要 登录 才可以下载或查看,没有账号?注册

×
发表于 2018-8-29 16:12:45 | 显示全部楼层 来自 河南郑州
Simdroid开发平台
感觉现在很少有人做这个了,关于晶体取向在ABAQUS中如何实现
回复 不支持

使用道具 举报

 楼主| 发表于 2018-9-21 16:19:00 | 显示全部楼层 来自 宁波大学
关于多晶取向的问题,还是通过写脚本来批量赋予比较容易实现,赋予取向后新的问题又出现了,"如何将取向和黄永刚老师的umat联系起来",这和以前帖子的版主提出的问题一样
回复 不支持

使用道具 举报

 楼主| 发表于 2018-9-21 16:20:58 | 显示全部楼层 来自 宁波大学
燃烧的向日葵 发表于 2018-9-21 16:19
关于多晶取向的问题,还是通过写脚本来批量赋予比较容易实现,赋予取向后新的问题又出现了,"如何将取向和 ...

还是历史遗留问题
回复 不支持

使用道具 举报

发表于 2018-9-25 14:20:44 | 显示全部楼层 来自 北京
弄过这个 我是修改的UMAT里那几个参数 不是orientation 很麻烦

点评

感谢回复,对于您说的修改Umat参数,能多讲一些这方面的工作吗?  发表于 2018-9-25 21:40
回复 不支持

使用道具 举报

发表于 2018-10-13 13:20:06 | 显示全部楼层 来自 江苏常州
关于如何给晶粒赋初始取向,黄永刚已经说的很清楚了,通过两个坐标系,一个旋转矩阵
回复 不支持

使用道具 举报

 楼主| 发表于 2018-10-25 18:41:10 | 显示全部楼层 来自 宁波大学
好久没更新了,现在多晶取向通过Python脚本不难实现,前处理不算大问题。主要问题还是在UMAT上,修改UMAT进行模拟,对比实验结果。(关键在能说明什么问题)
回复 不支持

使用道具 举报

发表于 2018-10-25 20:17:11 | 显示全部楼层 来自 河南郑州
兄dei,多晶取向你是怎么赋予的,能具体说一下嘛?
回复 不支持

使用道具 举报

 楼主| 发表于 2018-11-6 15:14:06 | 显示全部楼层 来自 宁波大学
赋予一个晶粒取向在rpy里面有对应的脚本,多晶取向则是将这个脚本循环多次,赋予随机取向的话要用到random函数
回复 不支持

使用道具 举报

发表于 2018-11-7 15:59:29 | 显示全部楼层 来自 山西晋中
cjm1224 发表于 2018-10-13 13:20
关于如何给晶粒赋初始取向,黄永刚已经说的很清楚了,通过两个坐标系,一个旋转矩阵 ...

你好,请教一下,两个坐标系和旋转矩阵都是随机的吗?还是先确定一个坐标系,根据某个旋转矩阵,再确定另外一个是吗?
回复 不支持

使用道具 举报

发表于 2018-11-14 10:19:37 | 显示全部楼层 来自 河南郑州
请问一下,初始硬化模量h0,阶段I的应力τs,初始屈服应力τ0这三个参数,你是怎么确定的
回复 不支持

使用道具 举报

发表于 2018-12-15 21:37:26 | 显示全部楼层 来自 哈尔滨工业大学
关于晶体塑性,目前有一些文献很好,可惜我的数学与力学功底不够。最有代表性的一篇是Overview of constitutive laws, kinematics, homogenization and multiscale methods in crystal plasticity finite-element modeling:Theory, experiments, applications。此外,https://icme.hpc.msstate.edu/med ... /Code:_ABAQUS_CPFEM中有一段程序,另有linux上的软件包DAMASK是开源的。
我目前也在朝这个方向努力,还没有进展。
回复 不支持

使用道具 举报

发表于 2019-3-22 15:53:48 | 显示全部楼层 来自 上海
弋阳 发表于 2018-12-15 21:37
关于晶体塑性,目前有一些文献很好,可惜我的数学与力学功底不够。最有代表性的一篇是Overview of constitu ...

我也是这个方向,可以加QQ交流下么,292502387
回复 不支持

使用道具 举报

发表于 2019-3-27 10:49:59 | 显示全部楼层 来自 江苏常州
晶界处的强化没有解决啊,晶界的产生就是因为取向的不同,故解决这两个就是多晶
回复 不支持

使用道具 举报

发表于 2020-1-17 11:20:34 | 显示全部楼层 来自 辽宁沈阳
太好的帖子
回复 不支持

使用道具 举报

发表于 2020-2-28 10:19:22 | 显示全部楼层 来自 日本
请问楼主有了解多晶的自洽模型吗?
回复 不支持

使用道具 举报

发表于 2020-7-22 21:03:52 | 显示全部楼层 来自 江苏南京
燃烧的向日葵 发表于 2018-11-6 15:14
赋予一个晶粒取向在rpy里面有对应的脚本,多晶取向则是将这个脚本循环多次,赋予随机取向的话要用到random ...

楼主可以指点一下吗
回复 不支持

使用道具 举报

发表于 2020-12-8 10:00:52 | 显示全部楼层 来自 黑龙江
小白弱弱的问一句,直接将黄永刚老师的inp,不加任何修改,输入2016版ABAQUS,可以运行成功吗
回复 不支持

使用道具 举报

发表于 2023-11-23 22:58:24 | 显示全部楼层 来自 辽宁沈阳
楼主你好,请问现在多晶的问题解决了吗?
回复 不支持

使用道具 举报

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

本版积分规则

Archiver|小黑屋|联系我们|仿真互动网 ( 京ICP备15048925号-7 )

GMT+8, 2024-4-18 12:22 , Processed in 0.046989 second(s), 16 queries , Gzip On, MemCache On.

Powered by Discuz! X3.5 Licensed

© 2001-2024 Discuz! Team.

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