1. Iteration 子程序里面的两个导数是什么。
产生newton-rhapson迭代方法的矩阵。
DHESLF
DHLATN
2. 欧拉角
位移变形率及节点力在总体坐标系中计算,弹性本构是在局部系中计算的。
从局部系转到整体系中通过欧拉角。
3. sets NSLPTL
set<110>{111} to set<110>{100}
set NSLIP
ISPNOR a typical normal to slip plane in this set of slip systems (integer) (INPUT) (1.0,1.0,.0.0)
ISPDIR -- a typical slip direction in this set of slip systems (integer) (INPUT) [1.0,1.0,1.0]
不能超过 three sets of slip systems
4. Denote s as a slip direction and m as normal to a slip plane.
C In a cubic crystal, (s,-m), (-s,m) and (-s,-m) are NOT considered
C independent of (s,m). 啥意思?
5. 问题:
(1) [012] nsdir为啥等于12呢?
(2) [111] 为啥是4个可能的方向呢?
(3) nslip(1)=12ge NSPTL=12还是36个?
Nset到底等于1还是3?
(4) m 和 s 如何后处理成极图?
6. X初始位置,x时间t后的位置。
F=Fe*Fp :变形梯度分解为弹性和塑性部分
L=dF/dt*inv(F) : 速度梯度等于偏v/偏x
L=dFe/dt*
dF/dt=Lp*Fp
Lp=gamma*m×乘n
LP=DP+WP=dF/dt*inv(Fp)
中间构形经历了剪切滑移和晶格旋转。
单晶塑性的流动法则:LP=剪切应变速率*m叉乘n
7. 结果处理
(1) 晶体取向的极图
(2)晶粒的取向变化
(3) 模拟织构演化
(4)多晶体塑性变形的各向异性
(2) 单晶体的实验
8. 参考文献
(1) introduction to dislocation
light Alloys:form traditional alloys to nanocrystals