位错密度梯度结构Cu单晶微柱压缩的三维离散位错动力学模拟
位错密度晶体塑性模型及在微成形模拟中的应用
的金属成形模拟中
上海交通大学工学博士论文答辩
1-5
课题意义及背景
1
晶体塑性模型及积分算法 1
位错密度晶体塑性模型
1
非局部位错密度晶体塑性模型 1
塑性微成形实验不模拟
1
结论不展望
1
有限变形晶体塑性理论 率相关晶体塑性积分算法研究现状 准隐式积分算法 单个单元压缩和剪切测试 等误差图分析 轧制模拟及ODF预测 筒形件拉深模拟及制耳预测 计算效率和稳定性 总结
FEM solver
Yes
1560(B) 1688(T)
8 ABAQUS/Standard
Yes
1602(B) 3079(T)
1
ABAQUS/Explicit
Yes
1602(B) 3079(T)
30 ABAQUS/Explicit
Yes
1602(B) 3079(T)
1
ABAQUS/Explicit
CPU time(s) 15380
当前构型Cauchy应力
2-1
有限变形晶体塑性理论
宏细观应力应变
● 位错沿着原子密排面-密排斱向(滑秱系)滑动 ● FCC晶体有12个滑秱系 ● 每个滑秱系由滑秱面的
法线 和滑秱斱向 组成
位错滑移引起的分切应变张量
上海交通大学工学博士论文答辩
Schmid张量
2-2
有限变形晶体塑性理论
现象学晶体塑性流动模型和硬化模型 ● 塑性功等效原理
隐式算法： Huang, 1991; Kalidindi et al., 1992 优点：无条件稳定，求解精度有保证 缺点：难收敛(无法得到准确的Jacobian矩阵)，计算效率低
显式算法： Grujicic and Batchu,2002; Rossiter et al., 2010. 优点：结构简单，本构斱程计算耗时最少 缺点：条件稳定，易发散，降低模型的整体增量步长，需要大量增量步
基于离散位错动力学的铜晶体材料力学特性研究
摘要近年来，随着航空航天、精密机械等产业的迅猛发展，对微机电系统需求量大幅增加，同时也对微机电系统提出了更高的性能要求。
铜晶体材料被广泛应用在微机电系统中作为基本构件，其尺寸在微米及亚微米尺度。
在该尺度下铜晶体表现出与尺度相关的力学特性，无法使用传统宏观力学进行解释。
因此，为了提高微机电系统性能及可靠度，迫切需要对小尺度铜晶体材料的力学性能开展研究，研究其与尺度相关的力学特性，为微机电系统的设计提供理论及实验基础。
针对小尺度下的材料仿真，离散位错动力学方法在拥有很高计算效率的同时，能分析晶体内部位错运动行为，是一种有效的介观尺度研究方法。
然而，针对小尺度材料测试中常用的微拉伸、纳米压痕两种实验方法，在离散位错动力学层面上对其研究并不完备。
本文中使用自主开发的离散位错动力学软件对这两种典型实验条件进行研究，对进一步了解小尺度铜晶体材料塑性力学特性的内在机理具有积极意义。
本文首先基于2.5D离散位错动力学框架，结合ABAQUS软件的接触算法，使用MATLAB软件及PYTHON语言开发了小尺度下铜晶体离散位错动力学-有限元耦合计算程序。
在该程序的基础上，开发相应的位错初始化模块、云图、曲线后处理模块。
在自主开发的铜晶体单向拉伸离散位错动力学-有限元耦合计算程序的基础上，建立了铜晶体单向拉伸有限元模型。
通过耦合计算，本文分析了铜晶体单向拉伸中的应力分布及内部位错运动机理。
针对铜晶体厚度、取向两方面分别展开研究，揭示了影响铜晶体单向拉伸性能的内在机理。
建立了铜晶体纳米压痕仿真模型，得到了纳米压痕过程中位错行为规律及压痕力-位移曲线。
明确了纳米压痕力-位移曲线与位错增殖速度及平均间距两个因素间的联系，并从位错源密度、初始位错密度两个方面研究其对纳米压痕结果的影响方式。
最后，本文对单晶铜进行纳米压痕实验，得到了单晶铜的力-位移数据。
使用扫描电子显微镜分析压痕周围的表面形貌，使用透射电子显微镜分析压痕周围位错线的分布规律，通过高分辨率透射电镜在原子堆垛层面上观察位错的具体状态，验证了仿真结果的正确性。
纳米多孔γ-TiAl单晶力学性能的分子动力学模拟
纳米多孔γ-TiAl单晶力学性能的分子动力学模拟左迪【摘要】利用分子动力学模拟方法,对多孔γ-TiAl单晶沿[001]晶向拉伸时的力学性能进行了研究,发现温度的增加会降低多孔γ-TiAl单晶的力学性能,屈服应力与杨氏模量均随温度的增加呈现降低的趋势.在多孔γ-TiAl单晶发生屈服后,剪切位错环在孔洞表面形核,导致塑性变形的发生.随着应变的增加,形成具有密排六方结构的双原子层厚的堆垛层错.位错抽取算法计算显示多孔γ-TiAl单晶中的位错均为1/6<112> Shockley不全位错,且随着应变的增加,位错密度呈增加的趋势.【期刊名称】《有色金属加工》【年(卷),期】2017(046)004【总页数】5页(P20-24)【关键词】分子动力学模拟;γ-TiAl单晶;屈服应力;位错【作者】左迪【作者单位】天水师范学院土木工程学院,甘肃天水 741001【正文语种】中文【中图分类】TG146.2从 20 世纪 70 年代以来，轻质的高温结构材料就受到高度重视，而金属间化合物由于其原子的长程有序排列和原子间金属键与共价键的共存性，使其能够同时兼顾金属的塑性和陶瓷的高温强度[1-3]。
γ-TiAl金属间化合物就是典型代表，其体积质量不到镍基合金的50%，具有轻质、高比强、高比刚、耐蚀、耐磨、耐高温以及优异的抗氧化性等优点，并具有优异的常温和高温力学性能，使用温度可达到700 ℃～1000 ℃，成为当代航空航天工业、兵器工业以及民用工业等领域优秀的候选高温结构材料之一，有望用于喷气发动机和涡轮等航空航天、汽车工业的耐高温部件以及超高速飞行器的翼、壳体等，具有重要的工程化应用潜力[4]。
但在现场室温环境时的钛铝合金由于材料自身不均匀、微裂纹或者环境腐蚀等因素影响使得其脆性较大，延展性较小，在拉伸试验过程中延伸率往往不足1%即会断裂，远小于理论值，机械加工非常困难，致使其难以在工厂企业中批量生产[5]。
梯度纳米孪晶Cu的力学性能和变形机制
梯度纳米孪晶Cu在冲击韧性测试中表现优异，其冲击功比传统粗晶Cu 高出约20%。
03
结论
梯度纳米孪晶结构在提高Cu的韧性的同时，对其在冲击载荷下的抗断
裂能力也有所增强。
强度测试
1 2
实验方法
采用拉伸试验机进行强度测试，设定试样尺寸为 1mm×10mm×20mm，以0.5mm/min的速率 进行拉伸，记录屈服强度和抗拉强度。
研究不足与展望
目前对梯度纳米孪晶Cu的力学性能和 变形机制的研究仍停留在实验阶段，缺 乏分子动力学模拟和第一性原理计算等
理论研究的支持。
对于梯度纳米孪晶Cu的制备工艺、微 观结构和力学性能之间的相互影响关系
仍需进一步深入研究。
未来可以探索梯度纳米孪晶Cu在不同 温度、应变速率和应力状态下的力学性 能和变形机制，为实际工程应用提供更
弹性变形
在较低应变下，梯度纳米孪晶Cu 表现出弹性变形行为，其弹性模 量取决于孪晶界的分布和相互作 用。
孪晶演化
初始孪晶
在梯度纳米孪晶Cu中，初始孪晶主 要通过堆垛层错和肖克莱不全位错的 形成而产生。
孪晶扩展
在应变作用下，孪晶扩展主要通过肖 克莱不全位错的连续运动来实现，同 时伴随着孪晶界的迁移。
梯度纳米孪晶Cu的力 学性能和变形机制
汇报人： 2023-11-20
目 录
• 引言 • 梯度纳米孪晶Cu的制备与表征 • 力学性能测试与分析 • 变形机制研究 • 力学性能与变形机制的关系 • 研究结论与展望 • 参考文献
01
引言
研究背景和意义
梯度纳米材料是一种具有梯度结 构的论
梯度纳米孪晶Cu具有优异的力学性能，包括高强度、高韧性和良好的抗疲劳性能。
动力学晶格蒙特卡洛方法模拟Cu薄膜生长
动力学晶格蒙特卡洛方法模拟Cu薄膜生长吴子若;程鑫彬;王占山【期刊名称】《光子学报》【年(卷),期】2010(39)1【摘要】利用动力学晶格蒙特卡洛方法模拟了Cu薄膜在Cu(100)面上的三维生长过程。
模型中考虑了四个动力学过程:原子沉积、增原子迁移、双原子迁移和台阶边缘原子迁移,各动力学过程发生的概率由多体势函数确定。
讨论了基底温度、沉积速率及原子覆盖率对Cu原子迁移、成核和表面岛生长等微观生长机制的影响;获得了Cu薄膜的表面形貌图并计算了表面粗糙度。
模拟结果表明,随基底温度升高或沉积速率下降,岛的平均尺寸增大,数目减少,形状更加规则。
低温时,Cu薄膜表现为分形的离散生长,高温时,Cu原子迁移能力增强形成密集的岛。
Cu薄膜表面粗糙度随着基底温度的升高而迅速减小;当基底温度低于某一临界温度时,表面粗糙度随原子覆盖率或沉积速率的增大而增大;当基底温度超过临界温度时,表面粗糙度随原子覆盖率或沉积速率的变化很小,基本趋于稳定。
【总页数】5页(P62-66)【关键词】薄膜生长;原子迁移;KLMC模拟;沉积速率;基底温度;原子覆盖率;表面粗糙度【作者】吴子若;程鑫彬;王占山【作者单位】同济大学精密光学工程技术研究所【正文语种】中文【中图分类】O484.1【相关文献】1.元胞方法与蒙特卡洛方法相结合的薄膜生长过程模拟∗ [J], 阮聪;孙晓民;宋亦旭2.Zr-Cu薄膜生长及力学性能的分子动力学模拟 [J], 谢璐;安豪杰;秦勤;臧勇3.基于动力学蒙特卡洛方法模拟化学气相沉积生长Cu薄膜的条件（英文） [J], 李朋;杨振雄;赵虎4.Cu表面生长Ag薄膜过程的分子动力学模拟 [J], 孙贺明;白照印;宋海洋;张国香5.动力学蒙特卡洛(KMC)模拟薄膜生长 [J], 周雪飞;吴冲;唐朝云;孔垂岗;邱贝贝;杨云;卢贵武因版权原因，仅展示原文概要，查看原文内容请购买。
板条裂纹问题的位错模拟解法
板条裂纹问题的位错模拟解法
王钟羡;马丽娜
【期刊名称】《机械强度》
【年(卷),期】2010()1
【摘要】假设沿分叉裂纹各分支和板条边界有某位错分布,利用半平面内分叉裂纹问题的复势函数,将板条分叉裂纹问题转化为半平面内的多分叉裂纹问题处理。
根据板条边界和裂纹面上的应力边界条件,建立以集中位错强度和分布位错密度为未知函数的Cauchy型奇异积分方程,利用半开型积分法则将该奇异积分方程化为代数方程求解。
最后,由位错密度函数得到各裂纹分支端的应力强度因子值。
文中分别给出集中力和分布力作用情况下内分叉和边缘分叉裂纹的3个算例,其极限情况的计算结果与精确解是一致的。
【总页数】5页(P139-143)
【关键词】分叉裂纹;板条;奇异积分方程;应力强度因子;位错
【作者】王钟羡;马丽娜
【作者单位】江苏大学理学院
【正文语种】中文
【中图分类】O346.1
【相关文献】
1.裂纹尖端位错发射与运动的分子动力学模拟 [J], 徐荣鹏;曾祥国;陈华燕;范镜泓
2.裂纹顶端的异号位错和无位错区的象力理论 [J], 龙期威;熊良钺
3.α-Fe纯Ⅰ型裂纹扩展中位错形成过程的三维分子动力学模拟 [J], 王建伟;尚新春;吕国才
4.裂纹尖端塑性区和无位错区的微观模拟 [J], 钱才富;姜忠军;陈平;段成红;崔文勇
5.位错和位错偶沿单─滑移系从裂纹尖端的发射 [J], 钱才富
因版权原因，仅展示原文概要，查看原文内容请购买。
211274798_FeCoCrNiCu高熵合金纳米压痕力学响应的分子动力学模拟★
加载速度
中图分类号:
O3
4
1
文献标识码:
A
/
DO
I:
1
0.
1
6
3
5
5
c
n
k
i.
i
s
s
n
1
0
0
7
9
4
3
2
t
u
t.
2
0
2
3.
0
3.
0
2
0
j.
y
文章编号:
1
0
0
7
9
4
3
2(
2
0
2
3)
0
3
0
5
6
2
0
8
Nano
i
nd
en
t
a
t
i
on Me
chan
i
c
a
lPr
ope
r
t
i
e
so
fFeCoCrNiCuHi
gh
En
t
r
opyAl
l
oy
s:
A Mo
l
e
c
u
l
a
rDynami
c
sS
t
udy
a
a
a,
b
a
a,
b
SUN Yewe
i ,ZHANG Lumi
ng ,MAShe
ng
i ,XIN Ha
o
guo ,ZHANG Tuanwe
(
a.
In
s
t
i
铜剧烈塑性变形微观缺陷演化机理的分子动力学模拟
铜剧烈塑性变形微观缺陷演化机理的分子动力学模拟
近些年来,随着科技的发展,对材料性能的要求越来越高。
而纳米金属材料在越来越多的领域都得到了广泛的应用,正是因为纳米晶体材料具有的特殊力学性质,如晶粒尺寸在纳米级别的材料,它的硬度,强度和延展性等都远远优于普通材料。
如果能够更加清楚地了解金属材料微观缺陷对其性能的影响,就可以通过改变材料的微结构来提高金属的性能。
因此,了解纳米金属材料的力学性质进而熟悉微结构变形机制,对于设计出满足人们需求的应用材料具有重大的意义。
本文从微观尺度出发,运用分子动力学模拟方法,对单晶Cu的剧烈塑性变形过程进行模拟计算,研究材料塑性变形过程中位错的微结构演化机理,并通过可视化软件对该材料不同阶段的微结构模拟过程进行分析。
研究了纳米面心立方金属的变形机制中位错结构的演化,即位错增殖、位错塞积、晶界的运动等等,从而探究其对宏观性能的影响,最后与单晶Ni的模拟结果进行对比。
模拟结果表明:在塑性变形过程中,位错主要以弗兰克—里德增殖源机制进行增殖;单晶Cu和Ni在拉伸阶段,滑移面与滑移方向趋于拉伸轴线方向,在压缩阶段,滑移面逐渐与压力轴线的方向垂直,同时结合铜和镍位错密度的对比发现,材料中出现的晶界越多,位错密度越大;对应力与密度的关系进行分析,得出材料塑性变形时,材料内部必然存在着位错的开动、增殖与消亡,而这些位错的运动最终达到结构细化的目的。
金属材料中位错动力学行为的数值模拟研究
金属材料中位错动力学行为的数值模拟研究金属材料是各种工业领域中最基础的材料之一。
在材料科学领域中，对金属材料中位错动力学行为的研究已经成为一个热门话题。
位错是晶体中的一个特殊的缺陷，它是晶体中的一条线，其周围的晶体结构发生了扭曲，从而导致原子的错位。
位错在各种材料的力学行为中起着至关重要的作用。
金属材料中位错的动力学行为是影响金属的力学性能的一个重要因素。
因此，不断地深化我们对位错动力学行为的研究，对于提高金属材料的力学性能具有非常重要的意义。
位错的动力学行为在过去几十年中一直是材料科学领域中的一个研究热点。
最初，人们主要是通过用透射电镜等成像技术来研究位错的行为。
而现在，随着计算机科学技术的不断发展，越来越多的科学家们开始使用计算机模拟技术研究位错的动力学行为。
这种方法不仅可以实现对复杂的三维结构进行模拟，而且可以掌握到位错行为的更为微观的细节。
位错的动力学行为中，位错滑移是其中的重要过程之一。
如果位错在晶体中进行滑动，那么材料就会发生变形。
位错的滑移通常受到多个因素的影响，比如温度、应力等。
在材料科学领域，人们通常使用分子动力学法、蒙特卡罗法、有限元法等来掌握位错滑移的动力学行为。
分子动力学法是一种基于牛顿力学原理的计算方法。
它可以非常精确地模拟分子之间的相互作用以及其运动方式。
在分子动力学法中，位错被看作是晶体中一连串不断变化的原子。
使用分子动力学法，科学家们可以观察到位错滑移的详细过程，这对于更好地分析位错的滑移机制具有重要作用。
蒙特卡罗法是一种模拟物理过程的方法，可以用于计算位错的滑移过程。
这种方法基于量子力学理论以及统计力学原理，能够模拟材料的温度、压强、化学的变化。
蒙特卡罗法主要是通过对各个原子位置的概率分布进行随机采样来模拟位错的运动。
有限元法是一种数值计算方法，可以分析任意形状的物体在外部力作用下的变形情况。
在使用该方法进行位错滑移的研究时，将晶体内部的位错看作是一个点，通过计算力学变形进行位错滑移的模拟。
位错密度演化模型
位错密度演化模型位错密度演化模型是研究材料结构演化和力学性质变化的重要工具。
位错是晶体中的一种结构缺陷，它类似于原子的错位，常常影响到材料的强度、韧性和塑性等力学性质。
位错密度演化模型的研究可以揭示位错在材料中的分布和演化过程，为改善材料的力学性能和设计新材料提供理论依据。
位错密度演化模型的基本原理是考虑位错的生成、迁移和消失过程。
位错可以通过多种途径生成，如外力加载、热处理、相变等。
当外力作用到材料上时，位错会被引入材料中，从而影响到材料的局部结构。
通过模拟位错的迁移过程，可以研究位错在材料中的运动规律和能量状态，进而预测材料的力学性质。
位错密度演化模型通常基于一系列微观动力学方程，包括位错的运动方程、能量方程和位错密度方程等。
这些方程可以描述位错的运动、相互作用和消失过程。
通过数值方法和计算机模拟，可以解析位错密度演化模型，得到位错密度随时间和空间的分布。
位错密度演化模型的研究对于了解材料力学性质变化的机理非常重要。
通过研究位错密度随时间和空间的演化，可以揭示位错的生成机制和运动规律，进而优化材料的力学性能。
例如，在材料设计中，可以通过调控位错密度的分布，来提高材料的强度和塑性。
同时，位错密度演化模型也可以用于材料的寿命评估和失效预测，有助于提高材料的可靠性和安全性。
位错密度演化模型的研究还可以与其他材料模型相结合，如晶体塑性模型、断裂模型等，来建立更加全面的材料性能预测和设计方法。
通过对不同材料体系的位错密度演化模型的比较和分析，可以揭示不同材料力学性质的本质差异，为材料的选择和设计提供指导。
总之，位错密度演化模型是研究材料结构演化和力学性质变化的重要工具。
通过模拟位错的生成、迁移和消失过程，可以揭示位错在材料中的分布和演化规律，为改善材料的力学性能和设计新材料提供理论依据。
位错密度演化模型的研究对于提高材料的可靠性和安全性，以及优化材料的力学性能具有重要指导意义。
Al-Cu合金凝固微观组织的三维模拟及优化
Al-Cu合金凝固微观组织的三维模拟及优化卜晓兵;李落星;张立强;朱必武;王水平【摘要】A 3D cellular automaton-finite element (CA-FE) method was used to predict the solidification structure of Al-Cu alloy at different casting radii, and optimized the solution diffusion coefficient in solid phase. The results show that when the casting radius changes from 10 mm to 30 mm, the proportion of columnar grains changes from 65.2% to 21.4%. With increasing heat transfer coefficient from 500 W/(m2·K) to 5 000 W/(m2·K), the proportion of columnar grains also gradually increases. The simulated results are in accord with the experimental ones well, and can accurately reflect the distribution, proportion, size of equiaxed and columnar grains.%采用CA-FE模型,对同一铸件不同尺寸的Al-2％Cu合金凝固微观组织进行三维模拟及预测,并引入空位形成能对固相扩散系数进行了优化.结果表明,当铸件直径分别为20、40和60 mm时,柱状晶占总晶体比率分别为65.2％、52.1％和21.4％,逐渐减少,而等轴晶所占比率逐渐增加；当换热系数由500 W/(m2.K)增大到5000W/(m2.K)时,同一铸件中不同直径的铸件凝固组织中柱状晶组织比例显著增大.模拟结果与实验结果吻合较好,能够较为准确地反映等轴晶和柱状晶的分布位置、比例和大小.【期刊名称】《中国有色金属学报》【年(卷),期】2011(021)009【总页数】7页(P2195-2201)【关键词】Al-Cu合金;3D元胞自动机;固相,扩散系数;凝固组织【作者】卜晓兵;李落星;张立强;朱必武;王水平【作者单位】湖南大学汽车车身先进设计与制造国家重点实验室,长沙410082;湖南大学材料科学与工程学院,长沙410082;湖南大学汽车车身先进设计与制造国家重点实验室,长沙410082;湖南大学材料科学与工程学院,长沙410082;湖南大学材料科学与工程学院,长沙410082;湖南大学材料科学与工程学院,长沙410082;湖南大学材料科学与工程学院,长沙410082【正文语种】中文【中图分类】TG292铸件中的细小晶粒能够提高材料的强度与韧性，因而预测微观组织形成是控制晶粒尺寸、形貌和铸件产品质量的重要前提。
Al2Cu拉伸变形的分子动力学模拟
基金项目：国家自然科学基金资助项目(11362017)；江西省自然科学基金资助项目(2010GQC0803)；江西省教育厅科技落地计划科学前沿项目 (KJLD12073)
收稿日期：2017-07-24；修订日期：2018-03-23 通信作者：刘晓波，教授，博士；电话：0791-83953108；E-mail：liuxb2000@
第 28 卷第 9 期
刘晓波，等：Al2Cu 拉伸变形的分子动力学模拟
1747
及 Al 含量对 Ti-Al 应力−应变关系及其塑性i
Fi (i )
1 2
ij (Rij )
j( ji)
(1)
尺寸下多晶银纳米线的拉伸变形行为，分析了晶粒尺
DOYAMA[11]采用分子动力学模拟单晶 Cu 和 Fe 的拉 伸过程，结果表明凹槽是位错和裂纹的来源。文玉华 等[12]在模拟纳米晶 Cu 的拉伸变形时，发现强度随晶 粒尺寸减小而降低。DIAO 等[13]通过分子动力学模拟 研究单晶金纳米丝的屈服机制，认为金纳米丝的塑性 变形主要是通过不全位错的运动进行。CAO 等[14]用分 子动力学模拟了纳米孪晶铜在受单向均匀拉伸载荷下 的变形机制，发现孪晶对位错滑移具有阻挡作用，使 得孪晶纳米线得到强化。KIM 等[15]通过模拟纳米多晶 Mg 的拉伸过程，认为应力条件的变化会引发滑移、 孪晶和晶界迁移等不同的变形机制。陈明等[16]通过分 子动力方法，分别对无孔洞和有孔洞的纳米单晶铜杆， 运用 EAM(Embedded atom model)势函数研究了其拉 伸特性。CHEN 等[17]研究了单轴拉伸载荷加载下，铜 (001)/镍(001)旋转晶界所形成的不同晶格失配网络对 其力学性能的影响。ZHANG 等[18]对含孪晶的铜纳米 线的分子动力学模拟结果表明，在拉伸载荷作用下， 孪晶界间距越小，纳米线在塑性变形区域内的应力越 低。WANG 等[19]基于 EAM 势，采用分子动力学方法 对超细镍纳米线在(100)晶向的拉伸性能进行研究，并 对其温度相关性和拉伸应变率相关性进行探讨。张晓 泳等[20]利用分子动力学模拟研究 Ti-Al 纳米杆的单向 拉伸变形过程，比较分析不同拉伸速率、拉伸温度以
基于位错密度变化的动态再结晶动力学确定方法_王健
平均位错密度的平方根成正比的关系求得。若 s 可
通过计算得到, 就能求得动态再结晶分数。实际上,
发生了动态再结晶的晶粒在长大过程中, 也要经历
加工硬化和动态回复过程, 当应变达到临界应变时, 会再次发生动态再结晶, 因此式 ( 1) 右侧第 二项
sX dyn 不能直接确定, 本文利用 Y anagim ot o 的增量 理论[ 6] , 经过推导, 并编写了程序进行计算, 使问
得到的流变应力曲线, 计算平均位错密度。得到一种有别于定量金相法的, 快速估算金属材料动态再结晶动力学的
方法, 运用该方法于普碳钢热压缩实验中, 在不同的变形条件下成功获得了其动态再结晶动力学方程, 并将结果和
先前 有关学者所做的研究结果进行比较得出, 本文方 法用于 模拟金 属热变 形过程中 的微观 组织演 变, 其在 实验和
X dyn = 1 - ex p[ - G( - c ) m ]
( 2)
c 动态再结晶临界应变
G 表征动态再结晶速率的材料参数与变
形金属的合金含量、变形温度 和应变
速度等条件有关
m
材料常数, 本文取 m= 2[ 5]
可以通过热压缩实验测得, n 可以利 用应变
小于临界应变时的流变应力曲线, 以及真应力值与
塑性工程学报
第 16 卷
型, 对于含有不同合金元素的钢种也并不通用。这 样, 在更短的时间内, 用更少量的实验来估算金属 材料动态再结晶动力学方程就显得很必要。本文给 出了一种利用流变应力曲线求解动态再结晶动力学 模型的方法, 以增量的形式建立了动态再结晶体积 分数与平均位错密度之间的联系, 并利用热压缩实 验得到的流变应力计算平均位错密度。这样可以直 接从热压缩实验得到的流变应力曲线上, 估算动态 再结晶动力学。研究中以 Avrami 曲线的形 式给出 了再结晶体积分数的表达式, 将其结果与有关学者 所作的研究做了对比分析。
位错动力学
牛顿位错动力学把晶体理想化为简谐近似中的正 则系综， 则系综 ， 每个位错线的应力场可以通过所有位错 段对应力贡献的线性叠加而得到。 段对应力贡献的线性叠加而得到。 位移场看作时间-空间相关的态变量 看作时间 空间相关的态变量， 把位移场看作时间 空间相关的态变量，把Hooke 定律和空位化学势作为状态方程。 作为状态方程 定律和空位化学势作为状态方程。 假设位错段满足局域力学平衡 条件， 假设 位错段满足局域力学平衡条件 ， 把 牛顿运动 位错段 满足局域力学平衡条件 定律和菲克扩散定律作为结构演化方程。 定律和菲克扩散定律作为结构演化方程。 作为结构演化方程
2 Newton位错动力学 位错动力学
2.2 位错分割原理
每一个位错线可以近似为在三维上相互连 接的、一段一段的直位错段的序列。 接的、一段一段的直位错段的序列。 由于位错不可能在晶体中无缺陷完善区域 终止，所以位错必须满足： 终止，所以位错必须满足： 位错段组成一个闭合环或者排列成半无限 位错条； 位错条；
Straight segment base
Lattice base
Spline segment base
ቤተ መጻሕፍቲ ባይዱ
2 Newton位错动力学 位错动力学
2.1 简介
• 考虑到弹性各向异性及晶体对称性，利用 Drown定理以及Asaro和Barneet理论中的积分 Drown定理以及Asaro和Barneet理论中的积分 定理以及Asaro 方法， 方法 ， 一般可以给出这些位错段的位移场和 应力场。 应力场。 • 每个位错线的应力场可以通过 所有位错段对 每个位错线的应力场可以通过所有位错段 对 所有位错段 线性叠加得到 应力贡献的线性叠加得到。 应力贡献的线性叠加得到。
纳米尺度单晶铜材料表面切削特性分子动力学模拟
， 微机械 、 微构件的疲劳强度 、
使用寿命及结构稳定性等都受到了影响 。 随着切削加工 技 术 的 发 展 和 进 步 ， 切削加工 精 度 得 到 了 明 显 的 提 高， 加 工 尺 度 则 不 断 减 小， 采用传统的试验方法来研究超精密切削加工过程 十分困难 ， 并且很多微观现象无法直接获得
０ 引言
由于纳米技术的发展 ， 微机械在航空航天 、 医 疗、 通信等领 域 得 到 了 广 阔 的 应 用 。 然 而 当 微 机 械、 微构件尺度减小到纳米量级 ， 由于尺寸效应和 表面效应的作用
［ １ － ２］
， 模拟方法是一种描述 微 观 现 象 的 有 效 方 ｉ ｃ ｓ ＭＤ） 法， 已成为研究纳米级切削机理的一个重要工具 。
这是由于采用较高的切削速切削速度50m切削速度100m切削速度200m时单晶铜缺陷分布及底部视图度导致单晶铜基体内部的位错形核及运动的时间较短而引起的13切削力反映切屑的去除过程是解释切削现象的重要物理参数下面根据模拟结果研究不同的切削速度对单晶铜材料表面切削特性的影响
— — 李 勇 杨晓京 纳米尺度单晶铜材料表面切削特性分子动力学模拟研究 —
基体原子与刀具原子 （ Ｃ ｕ －Ｃ）之间的作用采 用 Ｍｏ ｒ ｓ ｅ势函数来描述 ：
２ ｒ ｒ ｒ ｒ － － α（ α（ ｉ ｉ ０） ０）］ ｊ－ ｊ－ Ｅｖ（ ｒ ｅ ｅ ＝ ＥＤ ［ －２ ｉ ｊ）
１ 分子动力学建模与计算
１． １ 模型建立 建立单晶铜纳米切削过程三维分子动力学模 型， 如图 １ 所示 。
究了纳米切削过 程 中 ， 切削速度对切削过程的影 响规律 。 单晶铜 材 料 是 制 作 微 机 械 、 微构件常用 的材料 ， 具有良好的信号传输 、 塑性加工和抗疲劳 被广泛应用于国防精密仪器 、 民用通信设 等性能 ， ·７ ２ １·
《晶格动力学讲稿》课件
晶格动力学的经典理论模型
01
02
03
弹性力学模型
将晶格视为弹性连续体， 通过弹性力学理论来描述 晶格的动力学行为。
哈密顿量模型
将晶格中的原子或分子的 运动视为经典力学中的粒 子运动，通过哈密顿量来 描述晶格的动力学行为。
牛顿运动方程
通过牛顿运动方程来描述 晶格中原子或分子的运动 ，可以更精确地描述晶格 的动力学行为。
晶格振动与声子
晶格振动模式
晶格是由原子或分子构成的周期性结构，原子或分子的相对位置发生变化时会产生晶格振 动。晶格振动模式可以分为光学模式和声学模式。
声子与能量传递
声子在晶格中传播时，会与其他粒子发生相互作用，传递能量和动量。声子是晶格中能量 传递的重要载体。
声子的统计性质
在热力学平衡状态下，声子遵循Bose-Einstein统计分布，表现出粒子数的涨落和热辐射 等特性。了解声子的统计性质对于理解晶格动力学的热力学性质和相变机制具有重要意义 。
材料力学性质的预测
弹性模量
晶格动力学理论可以用于预测材料的弹性模量，通过分析晶 格振动对材料宏观力学性能的影响，可以计算出弹性模量的 大小。
强度与韧性
利用晶格动力学的理论模型，可以预测材料的强度和韧性， 了解材料在不同受力条件下的力学行为和破坏机制。
06
未来研究方向与展望
新型计算方法的探索
机器学习算法的应用
基态求解
密度泛函理论的核心是求解多电 子系统的基态，常用的方法有 Kohn-Sham方法和HartreeFock方法等。
交换关联泛函
密度泛函理论中，需要选择合适 的交换关联泛函来描述电子间的 相互作用，不同的泛函会导致不 同的计算结果。
05
激光温喷丸诱发单晶铜位错扩展的分子动力学模拟
第42卷第7期2015年7月Vol.42,No.7July,2015中国激光CHINESE JOURNAL OF LASERS 激光温喷丸诱发单晶铜位错扩展的分子动力学模拟孟宪凯周建忠黄舒盛杰苏纯杨祥伟江苏大学机械工程学院,江苏镇江212013摘要为了研究面心立方金属在热力耦合作用下的微观组织，采用分子动力学方法对激光温喷丸在单晶铜内诱导的位错扩展进行了模拟，并使用中心对称参数研究了变形量对位错扩展的影响规律，以及温度对塑性变形诱导位错扩展的作用，进而从加工硬化的角度分析了其对激光温喷丸强化效果的影响。
结果表明，激光温喷丸过程中，塑性变形以“空位簇-不全位错-堆垛层错”的方式诱导产生不全位错与堆垛层错。
随塑性变形增加，激光温喷丸的硬化效果逐渐增加。
并且，温度可以促进塑性变形过程中位错的形核与扩展，200℃以内，激光温喷丸的硬化效果随温度的增加而增强，温度增至250℃时，硬化效果由于位错湮灭而减弱。
关键词激光光学;激光温喷丸;位错扩展;分子动力学中图分类号TG17文献标识码Adoi:10.3788/CJL201542.0702003Molecular Dynamics Simulation of Dislocation Development inMonocrystalline Copper induced by Warm Laser PeeningMeng Xiankai Zhou Jianzhong Huang Shu Sheng Jie Su Chun Yang XiangweiSchool of Mechanical Engineering,Jiangsu University,Zhenjiang,Jiangsu 212013,ChinaAbstract In order to study the microstructure of face center cubic metals induced by thermo-mechanical impact,molecular dynamics method was used to investigate the dislocation evolution of monocrystalline copper treatedby warm laser peening (WLP).The effects of plastic strain on the dislocation development were researched throughcentrosymmetry parameter.Moreover,the effects of temperature on the dislocation development induced by plasticstrain during WLP were also investigated.Finally,the strengthening effects of WLP were analyzed by workhardening.It is shown that partial dislocation and stacking faults are induced by plastic strain in the way of “vacancygroup-partial dislocation-stacking faults ”.As the plastic strain increases,the amount of stacking faults and theirstacking extent increases,which leads to a better effect of WLP in work hardening.A higher temperature canpromote the nucleation and development of dislocations induced by WLP.The hardening effect of WLP increaseswith the temperature within 200℃,and the hardening effect reduces while it increases to 250℃because of thedislocation annihilation.Key words laser optics;warm laser peening;dislocation development;molecular dynamicsOCIS codes 140.3540;120.6810;160.3900;220.4610收稿日期:2014-12-15;收到修改稿日期:2015-01-26基金项目:国家自然科学基金(51175236,31405204)、高等学校博士学科点专向科研基金(20123227110022)、中国博士后科学基金(2014T70477,2013M540417)、江苏省科技支撑项目(BE2013097)、江苏省自然科学基金(BK20141166)作者简介:孟宪凯(1986—)，男，博士研究生，主要从事激光冲击强化方面的研究。
位错三维表征技术
部畸变会造成 X 射线相位发生变化, 利用迭代相位重建
位错的 三 维 应 变 场。 BCDI 技 术 无 需 使 用 衍 射 光 学 元
等缺陷的形貌 [38] 。 基于同步辐射 X 射线的位错三维表征
一技术可以在大于 180°倾转范围内收集衍射信息, 能够
泛的方法之一。 X 射线穿透晶体材料时, 晶体内部的局
技术特别是基于 TEM 的位错三维表征技术未来的发展趋
判断位错之间是否相交, 采用衍衬分析方法确定的柏氏
势和技术特点进行了展望。
矢量等位错晶体学信息无法与其空间几何信息相耦合,
图 1 不同尺度上位错的表征方法: ( a) 表面法光镜成像 [2] , ( b) 缀饰法光镜成像 [3] , ( c) X 射线成像 [2] , ( d) 扫描电子显微镜-电子
[6]
, 可
( 图 1) 。 然而, 采用上述表征 方 法 所 得 到 的 位 错 图 像,
ion beam-electron channeling contrast imaging, FIB-ECCI) 系
应用实例, 对比讨论上述重构技术在分辨率、 定量表征
或者是位错与样品表面的截面图像, 或者是位错结构在
滑移面和滑移方向等; ③ 能量参量, 如位错位移场、 应
418
第 40 卷
中国材料进展
力 / 应变场、 核芯能和弹性应变能及其空间分布等。 参量
这些局限性在很大程度上限制了人们对于众多位错相关
的多样性、 关联性和复杂性深刻影响着位错的形态和空
科学现象和规律的准确认识。 因此, 开发先进的位错多
间组态特征, 同时深刻影响着外场作用下位错的动态行
(2. Electron Microscopy Center, Chongqing University, Chongqing 400044, China)
梯度位错晶胞结构
梯度位错晶胞结构
梯度位错晶胞结构是一种金属材料的晶体结构，这种结构通常基于全位错的增殖运动、孪生、剪切或相变来协调。
梯度位错晶胞结构在塑性变形过程中激活了不全位错-层错诱导塑性变形机制。
变形初期，亚十纳米细小层错即从位错胞壁萌生、滑移扩展，其密度随拉伸应变增加而增加，逐渐演变成超高密度三维层错（与少量孪晶界）网格，直至布满整个晶粒。
超高密度细小层错/孪晶的形成有效协调其塑性变形、细化初始位错结构、阻碍其他缺陷运动而贡献强度和加工硬化。