CN105115998A - 一种粉末衍射择优取向参数确定方法 - Google Patents
一种粉末衍射择优取向参数确定方法 Download PDFInfo
- Publication number
- CN105115998A CN105115998A CN201510516202.6A CN201510516202A CN105115998A CN 105115998 A CN105115998 A CN 105115998A CN 201510516202 A CN201510516202 A CN 201510516202A CN 105115998 A CN105115998 A CN 105115998A
- Authority
- CN
- China
- Prior art keywords
- preferred orientation
- crystal
- angle
- parameter
- function
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Landscapes
- Analysing Materials By The Use Of Radiation (AREA)
- Developing Agents For Electrophotography (AREA)
Abstract
本发明涉及一种粉末衍射择优取向参数确定方法,本发明以粉末样品定向压实时产生择优取向为模型,确定择优取向函数,利用所确定的函数计算每个晶面的择优取向函数值;然后根据每个晶面的择优取向函数值得到晶体结构的计算谱图;比较计算谱图与实际谱图之间的误差,调整所设定的择优取向分布参数,直至误差小于设定值或不再减小,以调整后的择优取向分布参数作为该晶体的择优取向分布参数;最后利用该分布参数计算得到的择优取向值即为本发明所要确定的择优取向。本发明充分考虑了衍射强度与晶面取向分布概率成正比和受择优取向影响的总体衍射强度恒定个因素,所确定的择优取向函数更加符合粉末衍射的实际情况,所得到的择优取向也更加精准。
Description
技术领域
本发明涉及一种粉末衍射择优取向参数确定方法,属于技术领域。
背景技术
在一般多晶体中,每个晶粒有不同于邻晶的结晶学取向,从整体看,所有晶粒的取向是任意分布的;某些情况下,晶体的晶粒在不同程度上围绕某些特殊的取向排列,就称为择优取向或简称织构。择优取向会影响粉末衍射谱图的形态,在有关粉末衍射谱图的晶体学反演计算中,正确地描述择优取向的影响是非常的必要,择优取向对晶体学反演的影响如下:
Ip(hkl)=P(hkl)*I(hkl)
其中I(hkl)为原始计算衍射强度,Ip(hkl)为考虑择优取向影响后的衍射强度,P(hkl)为择优取向函数,择优取向函数描述了择优取向对晶面hkl衍射强度的影响。
P(hkl)=f(a,G)
其中为晶面hkl与衍射面的夹角,为锐角;G为择优取向分布状态参数,其意义主要包括择优取向的集中程度或强度。目前择优函数具体形式有以下6种,分别为D1~D6.
D1=exp(-Ga^2)
D2=exp(G(pi/2-a)^2)
D3=exp(-Gsin^2(a))
D4=exp(-G(1-cos^2(a))
D5=c+(1-c)exp(-Ga^2)
D6=(G^2cos^2(a)+sin^2(a)/G)^-3/2
Bragg-Brentano几何衍射面就是平行于样品台的平面,根据研究可知,只有平行衍射面的晶面才可以产生有效衍射,因此,在研究粉末样品择优取向时可只考虑法线包含于光路所在平面的晶面。在无择优取向时,晶面均匀分布,
P(hkl)=f(a,G)=1,-pi/2≤a≤pi/2
对上式进行积分可得
只考虑择优取向的影响,对于一个无择优取向的样品,假设一种特定情况,原来取向为a1的部分晶面,在有择优取向时,取向改变为a2,那么,在a1处P(hkl)应该有所降低,而在a2处P(hkl)应该有等量的增加,即择优取向函数的积分值应该恒定,如图1所示。而目前给出的择优取向函数均不满足上述要求,如图2所示,该图中给出了W.A.DOLLASE总结的5个单参数函数a>0右半边图形。目前的给出的择优取向函数均没有考虑受择优取向影响的总体衍射强度恒定的因素,导致所确定的择优取向不够准确。
发明内容
本发明的目的是提供一种粉末衍射择优取向参数确定方法,以解决目前由于没有考虑受择优取向影响的总体衍射强度恒定的因素导致所确定的择优取向不准确的问题。
本发明为解决上述技术问题提供了一种粉末衍射择优取向参数确定方法,该确定方法包括以下步骤:
1)测量某种晶体的单项衍射谱图,根据该衍射谱图确定晶体的结构;
2)根据晶体的结构分别计算晶体每个晶面与衍射面的夹角;
3)以粉末样品定向压实时产生择优取向为模型,确定择优取向函数,利用该函数和每个晶面与衍射面的夹角以及设定的择优取向分布参数计算每个晶面的择优取向函数值P,
P=(Gcos^2(a)+sin^2(a)/G)^-1
其中G为择优取向分布状态,0<G≤1,a为压实后晶面与衍射面的夹角;
4)根据每个晶面的择优取向函数值P得到晶体结构的计算谱图;
5)计算谱图与实际谱图之间的误差,判断误差是否小于设定值,若小于设定值,以此时所选定的择优取向分布参数作为该晶体的择优取向分布参数;
6)调整所设定的择优取向分布参数,重复步骤3)-5),直至误差小于设定值,以调整后的择优取向分布参数作为该晶体的择优取向分布参数;
7)以晶体的择优取向分布参数计算晶体衍射择优取向,该结果即为所要确定的择优取向值。
所述步骤3)中择优取向函数P的确定过程如下:
a)测量未压实时的晶面与衍射面的夹角a和定向压实时晶面与衍射面的夹角b;
b)以压实的比例参数H作为择优取向分布参数G,0<H≤1,H为0时表示最强压实,H为1时表示无压实;
c)根据夹角a和夹角b与比例参数H之间的位置关系计算夹角b概率密度函数,所计算出的概率密度函数即为所求的择优取向函数。
所述步骤4)中误差采用均方误差法计算得到。
本发明的有益效果是:本发明以粉末样品定向压实时产生择优取向为模型,重新确定择优取向函数,利用所确定的函数和每个晶面与衍射面的夹角以及设定的择优取向分布参数计算每个晶面的择优取向函数值P;然后根据每个晶面的择优取向函数值P得到晶体结构的计算谱图;比较计算谱图与实际谱图之间的误差,调整所设定的择优取向分布参数,直至误差小于设定值,以调整后的择优取向分布参数作为该晶体的择优取向分布参数;最后根据所确定的晶体择优取向分布参数计算该晶体的择优取向值即为本发明所要确定的择优取向。本发明充分考虑了受择优取向影响的总体衍射强度恒定的因素,所确定的择优取向函数更加符合粉末衍射的实际情况,所得到的择优取向更加精准。
附图说明
图1是局部择优取向改变示意图;
图2是现有技术中5个择优函数的分布示意图;
图3是定向施压产生择优取向示意图;
图4是本发明所确定的择优函数的分布示意图。
具体实施方式
下面结合附图对本发明的具体实施方式做进一步的说明。
本发明以粉末样品定向压实产生择优取向为模型,得到一个新的择优取向函数,与以往的择优取向函数相比,具备函数与晶面取向分布概率成正比和函数积分值恒定这两个特征,使根据该函数计算得到择优取向值更加精准。下面以某种晶体为例来说明本发明的具体实施过程。
1.测量某种晶体的单项衍射谱图,根据该衍射谱图确定晶体的结构。
2.根据晶体的结构分别计算晶体每个晶面与衍射面的夹角。
3.以粉末样品定向压实时产生择优取向为模型,重新确定择优取向函数,利用该函数和每个晶面与衍射面的夹角以及设定的择优取向分布参数计算每个晶面的择优取向函数值。
众所周知,在有片状、针状晶粒,如云母存在时,粉末样品定向压实会产生择优取向,如图3所示,只考虑其法线包含于光路所在平面的晶面,径向的直线表示片状晶面,垂直方向压缩了以后,会产生水平方向上的择优取向,是一种常见的择优取向分布状态。本发明以上述择优取向的分布状态为模型,以片状择优取向为例来说明择优取向函数的确定过程。
首先测量未压实时的晶面hkl与衍射面的夹角a和定向压实时晶面hkl与衍射面的夹角b,a和b均为锐角;
然后以压实的比例参数H作为择优取向分布参数G,G=H,0<G≤1,G为0时表示最强压实,G为1时表示无压实。在本实施例中G表示颓垣短轴长度与长轴长度的比值。
根据夹角a和夹角b与比例参数G之间的位置关系计算夹角b概率密度函数,所计算出的概率密度函数即为所求的择优取向函数。
b=atan(G*tan(a))
已知a在[0,pi/2]区间均匀分布,求解b的概率密度函数,以求解得到的b的概率密度函数作为择优取向函数P
P=(Gcos^2(a)+sin^2(a)/G)^-1
其中G为择优取向分布状态,0<G≤1,a为晶面与衍射面的夹角。
其图形如图4所示,由一个平台和叠加的一个对称峰形组成,G值增大,则平台升高,峰形降低,以蒙特-卡罗方法模拟100000个数据点统计验证函数,结论符合。
4.根据每个晶面的择优取向函数值P得到晶体结构的计算谱图,该谱图计算过程为行业内众所周知的,属于常规方法,这里不再详述。
5.比较计算谱图与实际谱图之间的误差,这里的误差计算采用均方误差,判断误差是否小于设定值,若小于设定值,以此时所选定的择优取向分布参数作为该晶体的择优取向分布参数;
6.调整所设定的择优取向分布参数,重复步骤3-5,直至误差小于设定值或不再减小,以调整后的择优取向分布参数作为该晶体的择优取向分布参数。
以上实施过程针对的是片状择优取向,对于线状择优取向也可以采用上述方法,只需将线状择优取向理解为是关于线轴法平面的片状择优取向,只需将片状择优取向函数中的夹角a替换为其余角pi/2-a即可。
Claims (3)
1.一种粉末衍射择优取向参数确定方法,其特征在于,该确定方法包括以下步骤:
1)测量待检测晶体的单项衍射谱图,根据该衍射谱图确定晶体的结构;
2)根据晶体的结构分别计算晶体每个晶面与衍射面的夹角;
3)以粉末样品定向压实时产生择优取向为模型,确定择优取向函数,利用该函数和每个晶面与衍射面的夹角以及设定的择优取向分布参数计算每个晶面的择优取向函数值P,
P=(Gcos^2(a)+sin^2(a)/G)^-1
其中G为择优取向分布状态,0<G≤1,a为压实后晶面与衍射面的夹角;
4)根据每个晶面的择优取向函数值P得到晶体结构的计算谱图;
5)计算谱图与实际谱图之间的误差,判断误差是否小于设定值,若小于设定值,以此时所选定的择优取向分布参数作为该晶体的择优取向分布参数;
6)调整所设定的择优取向分布参数,重复步骤3)-5),直至误差小于设定值,以调整后的择优取向分布参数作为该晶体的择优取向分布参数;
7)以晶体的择优取向分布参数计算晶体衍射择优取向,该结果即为所要确定的择优取向值。
2.根据权利要求1所述的粉末衍射择优取向参数确定方法,其特征在于,所述步骤3)中择优取向函数P的确定过程如下:
a)测量未压实时的晶面与衍射面的夹角a和定向压实时晶面与衍射面的夹角b;
b)以压实的比例参数H作为择优取向分布参数G,0<H≤1,H为0时表示最强压实,H为1时表示无压实;
c)根据夹角a和夹角b与比例参数H之间的位置关系计算夹角b概率密度函数,所计算出的概率密度函数即为所求的择优取向函数。
3.根据权利要求2所述的粉末衍射择优取向参数确定方法,其特征在于,所述步骤4)中误差采用均方误差法计算得到。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510516202.6A CN105115998B (zh) | 2015-08-20 | 2015-08-20 | 一种粉末衍射择优取向参数确定方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510516202.6A CN105115998B (zh) | 2015-08-20 | 2015-08-20 | 一种粉末衍射择优取向参数确定方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105115998A true CN105115998A (zh) | 2015-12-02 |
CN105115998B CN105115998B (zh) | 2017-11-28 |
Family
ID=54664038
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510516202.6A Active CN105115998B (zh) | 2015-08-20 | 2015-08-20 | 一种粉末衍射择优取向参数确定方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105115998B (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107590296A (zh) * | 2016-07-08 | 2018-01-16 | 深圳大学 | 一种小角x射线散射的全谱拟合方法及系统 |
CN108896591A (zh) * | 2018-05-11 | 2018-11-27 | 深圳晶泰科技有限公司 | 一种晶体粉末衍射的指标化方法 |
CN109916928A (zh) * | 2019-03-15 | 2019-06-21 | 浙江大学 | 一种判断各向异性纳米晶体择优取向性组装结果的方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CA822213A (en) * | 1969-09-02 | Esso Research And Engineering Company | Elimination of the preferred orientation distortions of an x-ray diffraction pattern | |
US3596092A (en) * | 1968-05-10 | 1971-07-27 | Rolls Royce | X-ray diffraction method for determining the value of a mechanical property or properties of a fiber |
EP1365230A2 (en) * | 2002-05-23 | 2003-11-26 | Rigaku Corporation | Method for estimating preferred orientation of polycrystalline material |
WO2011053692A1 (en) * | 2009-10-29 | 2011-05-05 | Glaxo Group Limited | Sample plate for an x-ray powder diffraction apparatus |
CN102721712A (zh) * | 2012-04-27 | 2012-10-10 | 中国航空工业集团公司北京航空材料研究院 | 一种用于测量金属管材织构的方法和装置 |
-
2015
- 2015-08-20 CN CN201510516202.6A patent/CN105115998B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CA822213A (en) * | 1969-09-02 | Esso Research And Engineering Company | Elimination of the preferred orientation distortions of an x-ray diffraction pattern | |
US3596092A (en) * | 1968-05-10 | 1971-07-27 | Rolls Royce | X-ray diffraction method for determining the value of a mechanical property or properties of a fiber |
EP1365230A2 (en) * | 2002-05-23 | 2003-11-26 | Rigaku Corporation | Method for estimating preferred orientation of polycrystalline material |
WO2011053692A1 (en) * | 2009-10-29 | 2011-05-05 | Glaxo Group Limited | Sample plate for an x-ray powder diffraction apparatus |
CN102721712A (zh) * | 2012-04-27 | 2012-10-10 | 中国航空工业集团公司北京航空材料研究院 | 一种用于测量金属管材织构的方法和装置 |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107590296A (zh) * | 2016-07-08 | 2018-01-16 | 深圳大学 | 一种小角x射线散射的全谱拟合方法及系统 |
CN108896591A (zh) * | 2018-05-11 | 2018-11-27 | 深圳晶泰科技有限公司 | 一种晶体粉末衍射的指标化方法 |
CN108896591B (zh) * | 2018-05-11 | 2020-10-16 | 深圳晶泰科技有限公司 | 一种晶体粉末衍射的指标化方法 |
CN109916928A (zh) * | 2019-03-15 | 2019-06-21 | 浙江大学 | 一种判断各向异性纳米晶体择优取向性组装结果的方法 |
Also Published As
Publication number | Publication date |
---|---|
CN105115998B (zh) | 2017-11-28 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US7937247B2 (en) | Virtual metering | |
CN104040378B (zh) | 气象预测装置以及气象预测方法 | |
CN103336891A (zh) | 一种用于配电网状态估计的伪量测生成方法 | |
CN109933852A (zh) | 预测车辆尺寸偏差的方法、装置、存储介质及电子设备 | |
CN108548833B (zh) | 一种混凝土脱空检测方法、装置及终端设备 | |
CN105115998A (zh) | 一种粉末衍射择优取向参数确定方法 | |
JP2016103972A (ja) | 電力網内の負荷をモデリングするためのシステムおよび方法 | |
CN112487732A (zh) | 一种冲积性河道冲淤模拟新方法及系统 | |
CN115854999A (zh) | 基于场景自适应的h-adcp断面平均流速自校正方法 | |
CN106556877B (zh) | 一种地磁通化方法及装置 | |
CN101988957B (zh) | 变电站合并单元电量采样准确性的动态检测方法及系统 | |
CN111998918A (zh) | 一种误差校正方法、误差校正装置及流量传感系统 | |
CN104777356A (zh) | 一种基于神经网络的实时高精度谐波检测方法 | |
CN110186533A (zh) | 一种高精度的河口短期潮位预报方法 | |
CN110715632B (zh) | 强潮河口海湾桥墩局部冲刷深度预测及预警的专用终端 | |
CN109783934A (zh) | 一种基于h-adcp的断面平均流速拟合率定方法 | |
CN104655496B (zh) | 自重对离轴反射镜面形影响的测试方法 | |
Gerritsen et al. | A modelling study of tidally induced equilibrium sand balances in the North Sea during the Holocene | |
CN104614144B (zh) | 泄洪消能诱发场地振动的预测方法 | |
CN103605323B (zh) | 化工生产的离散控制方法及装置 | |
CN105588569A (zh) | 一种利用影子定位的方法 | |
Wang et al. | An improved algorithm for state estimator based on maximum normal measurement rate | |
KR20150015200A (ko) | 제로-크로싱 방식의 전원 주파수 검출장치 및 제로-크로싱 방식의 전원 주파수 검출방법 | |
CN114353891A (zh) | 超声水流流量的计量方法、装置、电子设备及存储介质 | |
CN109387806B (zh) | 基于simulink仿真的数字化电能表溯源系统误差影响分析方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |