本发明涉及一种板材织构的取向,具体讲涉及一种修正板材织构欧拉角的方法。
背景技术:
研究板材的织构往往需要切制样品再用电子显微镜(sem)电子衍射。在扫描电子显微镜下,电子束作用在大角度倾向的样品表层区,在一次背散射电子与点阵面的相互作用中发生衍射,形成菊池花样。分析菊池花样的获知晶粒的晶体结构、取向等信息,其中取向信息通常用欧拉角
精度低,一般都是通过目测是否旋转到合理状态;
对于多个欧拉角同时存在偏差时,需要尝试很多次才能旋转到较理想的位置;
偏差角度太大时,如偏差角度大于45°时,容易旋转到垂直样品的方位,造成数据错误。
技术实现要素:
针对现有技术的不足,本申请人设计了一种修正板材织构欧拉角的方法,所述方法可以对欧拉角存在的偏差进行自动校正,大大减小了误差和复杂性;克服了大角度等难以调整的数据的调整难题;而且对欧拉空间内的所有欧拉角统一调整,大大提高了效率。
本发明的发明目的是通过下述技术方案进行实现的:
本发明人提供了一种修正板材织构欧拉角的方法,所述方法包括下述步骤:
(1)统计欧拉角在欧拉空间内的峰值中心的欧拉角;
(2)将所述步骤(1)中所得的欧拉角进行两两匹配后计算两个所述欧拉角间的位相差矩阵x;
(3)根据对称位相差矩阵x’求解所述旋转矩阵m;
(4)利用所述旋转矩阵m旋转所述欧拉空间内的所有欧拉角。
优选的,所述步骤(1)中的所述峰值为取向分布密度较大的欧拉空间坐标。
优选的,所述峰值的数量大于等于4。
优选的,所述步骤(2)中所述位相差矩阵x的计算方法为:
x=a1-1*a2
其中:a1、a2分别为所述步骤(1)中统计的欧拉角对应的矩阵。
优选的,所述步骤(3)中的对称位相差矩阵x’是从所述步骤(2)所得的所有位相差矩阵x中筛选得到的。
优选的,所述筛选包括:
①检查所述位相差矩阵是否为对称矩阵;
②检查所述位相差矩阵是否与矩阵λ相似
其中:
优选的,所述步骤①和所述步骤②中对矩阵的检查允许误差存在,选择所述误差最小的位相差矩阵x作为对称位相差矩阵x’。
优选的,所述步骤(3)中求解旋转矩阵m的具体方法为:
①设所述旋转矩阵
②根据公式
③根据公式φ=arccos(l)和l求出所述旋转矩阵m对应的欧拉角的第二个角度φ;
④根据公式
⑤将步骤(1)中所有欧拉角按照
⑥从所述步骤⑤中旋转后的欧拉角对应的矩阵中筛选出一组镜面对称的矩阵p和q;
所述
p13=q23、
p23=q13、
p33=q33和
(p32+q32)/(p31+q31)=-(p31-q31)/(p32-q32)=(p21+q22)/(p11+q21)=-(p11-p21)/(p12-q22)=(p22+q12)/(p21+q11)=-(p22-q12)/(p21-q11);
⑦根据公式
求出所述旋转矩阵的第三个角度
其中:p11、p12、p13、p21、p22、p23、p31、p32、p33、q11、q12、q13、q21、q22、q23、q31、q32和q33已知;
⑧根据所述欧拉角
优选的,所述步骤⑦中:
所述p13=q23的允许误差为0.001~0.01;
所述p23=q13的允许误差为0.001~0.01;
所述p33=q33的允许误差为0.001~0.01。
优选的,所述步骤⑤的具体操作为:
将待旋转的欧拉角对应的矩阵与所述欧拉角
优选的,所述步骤(4)中旋转的具体操作为:
将待旋转的欧拉角对应的矩阵与所述旋转矩阵m相乘;得到的矩阵所对应的欧拉角即为旋转之后的欧拉角。
优选的,将步骤(1)~步骤(4)编入计算机程序。
与最接近现有技术相比,本发明的有益效果在于:
1、本发明提供的技术方案,利用板材的对称型,通过计算实现了对欧拉角可能存在的偏差进行自动校正,避免了手动调整带来的误差,大大提高了校正的精度,可以控制在0.1°一下;克服了大角度等难以调整的数据的调整难题;而且对欧拉空间内的所有欧拉角统一调整,大大提高了效率。
附图说明
下面结合附图对本发明作进一步说明:
图1:本发明提供的方法的流程图;
具体实施方式
下面结合本申请实施例中的附图,对本申请实施例中的技术方案进行清楚、完整地描述。
如图1所示,本发明人提供了一种修正板材织构欧拉角的方法,所述方法包括下述步骤:
(1)统计欧拉角在欧拉空间内的峰值中心的欧拉角;
(2)将所述步骤(1)中所得的欧拉角进行两两匹配后计算两个所述欧拉角间的位相差矩阵x;
(3)根据对称位相差矩阵x’求解所述旋转矩阵m;
(4)利用所述旋转矩阵m旋转所述欧拉空间内的所有欧拉角。
所述步骤(1)中的所述峰值为取向分布密度较大的欧拉空间坐标。
所述峰值的数量大于等于4。
所述步骤(2)中所述位相差矩阵x的计算方法为:
x=a1-1*a2
其中:a1、a2分别为所述步骤(1)中统计的欧拉角对应的矩阵。
所述步骤(3)中的对称位相差矩阵x’是从所述步骤(2)所得的所有位相差矩阵x中筛选得到的。
所述筛选包括:
①检查所述位相差矩阵是否为对称矩阵;
②检查所述位相差矩阵是否与矩阵λ相似;
其中
所述步骤①和所述步骤②中对矩阵的检查允许误差存在,选择所述误差最小的位相差矩阵x作为对称位相差矩阵x’。
所述步骤(3)中求解旋转矩阵m的具体方法为:
①设所述旋转矩阵
②根据公式
③根据公式φ=arccos(l)和l求出所述旋转矩阵m对应的欧拉角的第二个角度φ;
④根据公式
⑤将步骤(1)中所有欧拉角按照
⑥从所述步骤⑤中旋转后的欧拉角对应的矩阵中筛选出一组镜面对称的矩阵p和q;
所述
p13=q23、
p23=q13、
p33=q33和
(p32+q32)/(p31+q31)=-(p31-q31)/(p32-q32)=(p21+q22)/(p11+q21)=-(p11-p21)/(p12-q22)=(p22+q12)/(p21+q11)=-(p22-q12)/(p21-q11);
⑦根据公式
求出所述旋转矩阵的第三个角度
其中:p11、p12、p13、p21、p22、p23、p31、p32、p33、q11、q12、q13、q21、q22、q23、q31、q32和q33已知;
⑧根据所述欧拉角
所述步骤⑦中:
所述p13=q23的允许误差为0.001~0.01;
所述p23=q13的允许误差为0.001~0.01;
所述p33=q33的允许误差为0.001~0.01。
优选的,所述步骤⑤的具体操作为:
将待旋转的欧拉角对应的矩阵与所述欧拉角
所述步骤(4)中旋转的具体操作为:
将待旋转的欧拉角对应的矩阵与所述旋转矩阵m相乘;得到的矩阵所对应的欧拉角即为旋转之后的欧拉角。
将步骤(1)~步骤(4)编入计算机程序。
最后应该说明的是:所描述的实施例仅是本申请一部分实施例,而不是全部的实施例,基于本申请中的实施例,本领域技术人员在没有做出创造性劳动的前提下所获得的所有其它实施例,都属于本申请保护的范围。