基于Huang的晶体塑性有限元-CPFEM-入门-6

同学们，关于老师第6课中的程序修改，“H K L 0 0 1”改为“0 0 1 H K L”，这里是根据模拟者使用的欧拉角定义来的（是从局部转全局？还是全局转局部？），所以大家要根据自己的欧拉角数据来源的定义来判断要不要做出这个改动，而不是机械的死搬老师的操作。

老师，能否讲讲黄教授的umat模型，并如果修改成适合自己的模型，谢谢您。

谁有修改后的UMAT文件呀，我修改以后报错，救救孩子，有偿也行

老师，希望更新一期EBSD的数据导入CPFEM中计算，还希望讲解一下材料本构中参数的拟合与确定，部分模型数据试验很难获取，就用拟合的方法，我看有一部分用Isihgt做的，希望老师讲解一下。

为什么我按照视频里的弄了4个晶粒，拉伸0.5mm，用了97年那版的子程序，内存占用只有16g,abaqus却未反应，也没报错也没警告出现

老师的微信号是多少啊

老师，关于这个提取积分点的滑移方向的程序可以发我一份吗，我这一直报错。谢谢。965339756@qq.com

想请问一下老师。我在for里面按这个 open(10,file='d:\temp.dat',status='old', 2 ACCESS='append')改好后仍然报错 forrtl : severe (104) : incorrect access = specifier value for connected file, unit 10, file d:\temp.dat

老师，请问参数标定利用abaqus与isight联合的话，如何实现晶体整个模型的参数化建模呢？

很好奇 最后算出来有12个滑移系的 晶面方向 跟晶向 那我应该取哪一组的米勒因子来表示晶体最后的取向？还有算出来的应该是全局坐标系下的一组方向吧 就比如说我取第一组｛111｝《0-11》的 最终法线方向跟滑移方向，那是不是应该求｛111｝《0-11》到你取的结果之间的旋转矩阵g然后在把旋转矩阵g转换成欧拉角？