一种应力双折射计算方法技术

一种应力双折射计算方法，包括步骤：确定样品的比例极限及临界载荷；建立样品的有限元仿真并提取节点的空间应力数据；确定光传输方向，并进行坐标系转化；提取垂直于光传播方向的正应力及切应力，计算垂直于光传输平面的主应力和次主应力；计算样品的双折射数值并判定快轴角。通过样品的比例极限及临界载荷来确定是否适用光弹性方程计算样品的双折射数值，进一步给出主应力和次主应力的具体计算公式，并给出快轴角的计算公式，及快轴角的判定方法，解决了现有应力双折射计算方法定义不清晰，计算不详细的问题。算不详细的问题。算不详细的问题。

【技术实现步骤摘要】
一种应力双折射计算方法


[0001]本专利技术涉及应力双折射
，具体涉及一种应力双折射计算方法。

技术介绍

[0002]在计算应力双折射时，常采用光弹性方程，经专利技术人研究发现，光弹性方程并不是适合所有样品的应力双折射计算，因此，仍需要对应力双折射计算方法进行深入研究，界定光弹性方程的应用范围，明确主应力的具体计算方式以及明确快轴角的具体判定方法。

技术实现思路

[0003]针对应力双折射计算过程中存在的问题，本申请提供一种应力双折射计算方法，通过明确光弹性方程的应用范围，明确主应力的具体计算方式，以及明确快轴角的具体判定方法，以解决现有应力双折射计算方法定义不清晰，计算不详细的问题。
[0004]本专利技术的技术方案如下：
[0005]本专利技术提供一种应力双折射计算方法，包括步骤：
[0006]确定样品的比例极限及临界载荷；
[0007]建立样品的有限元仿真并提取节点的空间应力数据；
[0008]确定光传输方向，并进行坐标系转化；
[0009]提取垂直于光传播方向的正应力及切应力，计算垂直于光传输平面的主应力和次主应力；
[0010]计算样品的双折射数值并判定快轴角。
[0011]进一步优选的，外加载荷小于所述比例极限及临界载荷，使应力
‑
应变
‑
位移为线性化关系，以采用光弹性方程计算样品的双折射数值。
[0012]进一步优选的，所述空间应力数据包括正应力数据和切应力数据。
[0013]进一步优选的，所述坐标系转化为空间坐标系转化，由笛卡尔坐标系转化为光传输方向坐标系。
[0014]进一步优选的，所述计算垂直于光传播方向的主应力及次主应力具体为计算垂直于光传播方向平面内的平面主应力及次主应力。
[0015]进一步优选的，所述计算垂直于光传播方向的主应力及次主应力，具体计算公式为：
[0016][0017][0018]其中，σ1为主应力，σ2为次主应力，σ
x
、σ
y
为正应力数据，τ
xy
为切应力数据。
[0019]进一步优选的，所述计算样品的双折射数值，计算公式为：Δn＝C(σ1‑
σ2)，其中，Δn为双折射数值，C为常数。
[0020]进一步优选的，所述快轴角的计算公式为：
[0021]2θ＝arctan(2τ
xy
/(σ
x
‑
σ
y
))+2kπ；
[0022]其中，θ为快轴角，k为正整数，θ的判定方式为：
[0023]当σ
x
>σ
y
，θ＝
‑
θ；
[0024]当σ
x
<σ
y
，若τ
xy
<0，则θ＝(π
‑
2θ)/2，若τ
xy
>0，则θ＝(
‑
π
‑
2θ)/2。
[0025]依据上述实施例的应力双折射计算方法，通过样品的比例极限及临界载荷来确定是否适用光弹性方程计算样品的双折射数值，进一步给出主应力和次主应力的具体计算公式，并给出快轴角的计算公式，及快轴角的判定方法，解决了现有应力双折射计算方法定义不清晰，计算不详细的问题。
附图说明
[0026]图1为应力双折射计算流程图。
[0027]图2为光学平板及外加载荷图。
[0028]图3为光学平板有限元仿真图。
[0029]图4为光学平板应力双折射仿真图。
[0030]图5为光学平板快轴角分布图。
具体实施方式
[0031]下面通过具体实施方式结合附图对本专利技术作进一步详细说明。
[0032]应力双折射：在外加力的影响下，原本各项同性的光学材料转为各项异性，产生双折射效应，并且在外力撤去之后，材料又恢复各项同性，双折射效应消失，这种效应成为应力双折射。
[0033]比例极限：弹性变形阶段，金属材料所承受的应力和应变保持正比的最大应力。
[0034]临界载荷：在载荷加载过程中产生较大的位移及较大的旋转，使得几何方程由线性转为非线性。
[0035]光弹性方程：Δn＝C(σ1‑
σ2)，其中Δn为双折射率，σ1，σ2分别为第一主应力和第二主应力。
[0036]快轴：波片中传播速度快的光矢量方向。
[0037]现有的应力双折射计算通常是根据光弹性方程进行计算，计算公式如下：
[0038]Δn＝C(σ1‑
σ2)；
[0039]tan2θ＝2τ
xy
/(σ
x
‑
σ
y
)。
[0040]但是实际计算中并没有对上述公式进行明确的定义，例如，什么情况下适用光弹性方程计算应力双折射，光弹性方程中的σ1、σ2的具体计算会受哪些因素影响，总之，目前对上述公式的应用仍比较模糊，没有明确的定义，专利技术人经过深入研究，提出了一种应力双折射计算方法，通过样品的比例极限及临界载荷来确定是否适用光弹性方程计算样品的双折射数值，进一步给出主应力和次主应力的具体计算公式，并给出快轴角的计算公式，及快轴角的判定方法，解决了现有应力双折射计算方法定义不清晰，计算不详细的问题。
[0041]本申请提供的应力双折射计算方法流程图如图1所示，具体包括以下步骤。
[0042]S100：确定样品的比例极限及临界载荷。
[0043]S200：建立样品的有限元仿真并提取节点的空间应力数据。
[0044]S300：确定光传输方向，并进行坐标系转化。
[0045]S400：提取垂直于光传播方向的正应力及切应力，计算垂直于光传输平面的主应力和次主应力。
[0046]S500：计算样品的双折射数值并判定快轴角。
[0047]在步骤S100中，通过确定样品的比例极限及临界载荷来确定是否适合采用光弹性方程计算样品的双折射数值，例如，当外加载荷小于该比例极限及临界载荷时，使应力
‑
应变
‑
位移线性化，则可以采用光弹性方程计算样品的双折射数值。
[0048]例如，如图2所示，所采用的仿真样品为直径188mm，厚度25mm的熔融石英光学平板，0
°
(y轴负方向)和270
°
(y轴正方向)指向圆心各加载一个点载荷，大小为1MPa；熔融石英的比例极限为3430MPa，仿真所得的临界载荷为1.61MPa；因此外加载荷符合要求，可以采用光弹性方程计算双折射。
[0049]在步骤S200中，建立样品的有限元仿真并提取节点的空间应力数据，以图2中的仿真样品为例，则有限元仿真图如图3的所示，所提取的空间应力数据包括：σ
x
、σ
y
、σ<本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种应力双折射计算方法，其特征在于，包括步骤：确定样品的比例极限及临界载荷；建立样品的有限元仿真并提取节点的空间应力数据；确定光传输方向，并进行坐标系转化；提取垂直于光传播方向的正应力及切应力，计算垂直于光传输平面的主应力和次主应力；计算样品的双折射数值并判定快轴角。2.如权利要求1所述的应力双折射计算方法，其特征在于，外加载荷小于所述比例极限及临界载荷，使得应力
‑
应变
‑
位移为线性化关系，以采用光弹性方程计算样品的双折射数值。3.如权利要求1所述的应力双折射计算方法，其特征在于，所述空间应力数据包括正应力数据和切应力数据。4.如权利要求1所述的应力双折射计算方法，其特征在于，所述坐标系转化为空间坐标系转化，由笛卡尔坐标系转化为光传输方向坐标系。5.如权利要求1所述的应力双折射计算方法，其特征在于，所述计算垂直于光传播方向的主应力及次主应力具体为计算垂直于光传播方向平面内的平面主应力及次主应力。6.如权利要求1所述的应力双折射计算方法，其特征在于，所述计算垂直于光传播方向的主应力及次主应力，具体计算公式为：的主应力及次主应力，具体计算公式为：其中，σ1为主应力，σ2为次...

【专利技术属性】
技术研发人员：朱玲琳，曾爱军，胡敬佩，黄惠杰，
申请(专利权)人：上海镭望光学科技有限公司，
类型：发明
国别省市：