CN111324933B - 抗振型光学参考腔振动敏感度分析及设计参考腔的方法 - Google Patents
抗振型光学参考腔振动敏感度分析及设计参考腔的方法 Download PDFInfo
- Publication number
- CN111324933B CN111324933B CN202010103216.6A CN202010103216A CN111324933B CN 111324933 B CN111324933 B CN 111324933B CN 202010103216 A CN202010103216 A CN 202010103216A CN 111324933 B CN111324933 B CN 111324933B
- Authority
- CN
- China
- Prior art keywords
- reference cavity
- support
- model
- vibration
- cavity
- 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.)
- Active
Links
Images
Classifications
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
Abstract
本发明公开了抗振型光学参考腔振动敏感度分析及设计参考腔的方法,以第一主加速度方向为一维坐标轴,坐标数值大的支撑物进行三维固定约束,坐标数值小的支撑物施加挤压力,计算第一参考腔腔长形变量;建立与第一参考腔模型完全相同的第二参考腔模型,只改变主加速度载荷值,其他设置均不变,重复以上步骤进行仿真计算,计算第二参考腔腔长形变量;最后,计算拟分析的参考腔的振动敏感度,根据上述计算步骤,计算各个参数化变量下的振动敏感度,选出参考腔最佳的设计参数,本发明能准确的反应参考腔的实际受力情况,充分考虑了各个参数化变量的情况,具备最佳的设计参数。
Description
技术领域
本发明属于超稳光学参考腔技术领域,具体属于抗振型光学参考腔振动敏感度分析及设计参考腔的方法。
背景技术
具有极高频率稳定度的窄线宽激光器作为高精密测量的一种手段,在基本物理常数测量、引力波探测、测地学、原子光钟等科学与技术领域有着广泛的应用。实现窄线宽激光器的方法主要有两种,一种是利用原子或分子能级跃迁谱线的稳定性,如饱和吸收谱等,将激光频率锁定在跃迁谱线中心频率上,达到激光稳频的目的,该方法得到的激光频率稳定度约为10-13量级;另一种是基于光学参考腔的Pound-Drever-Hall(PDH)稳频技术,该方法将相位调制光谱技术与光外差探测技术结合,把激光频率精密锁定在光学参考腔的共振频率上,具有鉴频信号强、中心频率处斜率大、控制范围宽等特点,广泛用于超稳激光的产生。使用PDH稳频技术进行激光稳频时,是将激光频率锁定在参考腔的共振频率上,该共振频率的稳定性取决于参考腔的长度稳定性,因此参考腔的长度稳定性决定了稳频激光的频率稳定性。而振动导致的参考腔长度稳定性变化已经成为限制超稳激光性能进一步提升的关键因素之一。因此计算参考腔振动敏感度成为了研制超稳激光的关键一环。
目前大多数超稳激光都是仅工作在实验室环境中,参考腔的一般固定方式是将参考腔放置三个或者四个支撑面上,在重力约束下,使其保持平衡状态。2006年Chen等人首次提出利用有限元分析软件通过数值模拟的方法分析计算参考腔的振动敏感度(2006,PHYSICAL REVIEW A,74,053801,)。随后世界各研究小组均采用数值模拟的方法分析参考腔的振动敏感度,从而优化设计参考腔结构及支撑位置。具体操作方法是:第一,建立参考腔的几何模型,主要包括参考腔腔体、反射镜(由于反射镜表面镀膜很薄,一般不需要单独对镀膜层建模);第二,选取参考腔的支撑面,并对所有支撑面进行三维固定设置,同时增加外部载荷——主加速度;第三,对模型进行网格划分;第四,运行有限元分析软件,计算参考腔形变;第五,以两反射镜内表面中心位置的距离表示参考腔长度,仿真前后该距离沿光轴方向的变化量ΔL与参考腔原长度L的比值ΔL/L,即为该参考腔的振动敏感度,单位为g-1,g为主加速度。
随着航天技术的发展,美国、欧盟、日本等世界发达国家和组织争相开展在太空领域的光钟研究计划,以争取获得时间频率的制高点,达到导航、卫星、国防更高精度的应用,抗振型超稳激光越来越受到各研究小组的青睐,其中抗振型超稳光学参考腔成为该研究领域的热点。2011年Webster等人提出了一种基于立方腔体四点对称挤压支撑的光学参考腔(2011,Opt.Lett.36,3572),如图1所示,该立方体参考腔八个顶角均沿体对角线进行切割,切割成八个正三角形面,对其中构成正四面体的四个顶角面使用特氟龙半球进行支撑,由于该参考腔进行了三维约束,具有较好的抗振效果,作者使用有限元分析软件分析了四个顶面挤压力对参考腔长度变化的影响,并从实验上测得该参考腔最大的振动敏感度为2.5×10-11/g,该振动敏感度为已知参考腔中振动敏感度最低的参考腔。2018年华东师范大学研究了类似100mm边长的立方体参考腔,并且在8个顶角面的正三角形平面中心挖圆锥孔,类似的于其中四个圆锥孔进行固定支撑,在进行振动敏感度的有限元分析时,作者将四个支撑锥面进行三维固定设置,并施加主加速度,仿真分析腔长的变化。
以立方体形状的抗振型参考腔为例,当前抗振型参考腔在进行振动敏感度有限元仿真分析时,主要存在以下两个方面的不足之处:第一,将参考腔所有支撑面都加载三维固定约束不能较为准确的体现参考腔受力情况,如立方体参考腔,所有支撑面与腔体之间应为挤压力,而若将四个锥面固定,在重力作用下局部支撑面对参考腔的作用力则为拉力;第二,将预紧力导致的参考腔腔长形变叠加于主加速度导致的腔长形变,无法体现仅由主加速度导致的参考腔腔长形变量。
发明内容
为了解决现有技术中存在的问题,本发明提供一种抗振型光学参考腔振动敏感度分析方法,解决现有技术中的超稳光学参考腔振动敏感度的有限元分析方法不能准确反映光学参考腔的受力问题;
为实现上述目的,本发明提供如下技术方案:
一种抗振型光学参考腔振动敏感度分析方法,首先,对第一参考腔模型施加第一主加速度,以第一主加速度方向为一维坐标轴,坐标数值大的支撑物进行三维固定约束,坐标数值小的支撑物施加挤压力约束,仿真计算得到第一参考腔腔长形变量ΔL1;仿真计算得到第一参考腔腔长形变量ΔL1;
然后,建立第二参考腔模型,所述第二参考腔模型与第一参考腔模型完全相同,对第二参考腔模型施加第二主加度,所述第二主加速度的载荷值与第一主加速度的载荷值不同但方向相同,第二参考腔模型上支撑物的约束方式与所述第一参考腔模型上支撑物约束方式相同,仿真计算得到第二参考腔腔长形变量ΔL2;
最后,计算振动敏感度,单位为g-1:
(ΔL2-ΔL1)/L/(a2-a1),
式中:L为第一参考腔模型或第二参考腔模型未发生形变的参考腔腔长,a1为第一主加速度,a2为第二主加速度。
进一步的,具体包括以下步骤:
1)根据拟分析的参考腔建立分析建模,作为第一参考腔模型;
2)对第一参考腔模型施加第一主加速度,以第一主加速度方向为一维坐标轴,坐标数值大的支撑物进行三维固定约束,坐标数值小的支撑物施加挤压力约束,或以第一主加速度方向为一维坐标轴,坐标数值大的支撑物仅在垂直于支撑面方向进行一维固定约束;
3)根据第一参考腔模型的尺寸对第一参考腔模型进行单元格划分;
4)通过步骤3)单元格划分后的第一参考腔模型,计算第一参考腔腔长形变量ΔL1;
5)建立与步骤1)完全相同的参考腔模型为第二参考腔模型;
6)对第二参考腔模型施加第二主加速度,所述第二主加速度的载荷值与第一主加速度的载荷值不同但方向相同,第二参考腔模型上支撑物的约束方式与第一参考腔模型上支撑物约束方式相同;
7)根据第二参考腔模型的尺寸对第二参考腔模型进行单元格划分;
8)通过步骤7)单元格划分后的第二参考腔模型,计算得出第二参考腔腔长形变量ΔL2;
计算拟分析的参考腔振动敏感度,单位为g-1:
(ΔL2-ΔL1)/L/(a2-a1),
式中:L为第一参考腔模型或第二参考腔模型未发生形变的参考腔腔长,a1为第一主加速度,a2为第二主加速度。
进一步的,所述步骤3)中对模型进行单元格划分时,对支撑物和腔镜进行网格细化。
进一步的,所述步骤3)和步骤7)中所述单元格的尺寸进行均匀过渡设置。
进一步的,所述挤压力包括对所述支撑物的外端面的法向施加挤压力。
本发明还提供了一种抗振型光学参考腔振动敏感度的设计参考腔的方法,计算若干种参数化变量下的参考腔振动敏感度;
拟合出参考腔振动敏感度随各个参数化变量变化的曲线,根据所述曲线中振动敏感度最接近零或者为零的数值选出对应的各个参数化变量的数值为参考腔最佳设计参数。
进一步的,所述参数化变量包括参考腔体的切割尺寸,支撑物和支撑面的接触形状和接触面的面积,支撑物厚度以及支撑物材料属性中至少一项。
进一步的,所述参考腔体的切割尺寸包括沿立方腔体对角线切割顶角的深度或参考腔围脖的厚度及内径和外径尺寸、腔镜的直径、腔镜的厚度、参考腔补偿环厚度及内径和外径。
与现有技术相比,本发明至少具有以下有益效果:
参考腔在支撑物的挤压下,才能实现整个参考腔的三维约束,本发明提供了一种在有限元软件上对抗振型光学参考腔振动敏感度分析方法,改变了现有技术中将所有支撑面都固定的约束方式,此种固定约束方式,在主加速度的作用下,被固定约束的支撑物对参考腔的作用表现为拉力,而不是实际中的挤压力,本发明充分考虑了参考腔实际的受力情况,对一部分支撑面施加挤压力约束,同时考虑到参考腔在主加速度作用下,会有沿主加速度运动的趋势,为了防止参考腔运动,以第一主加速度方向为一维坐标轴,坐标数值大的支撑物进行了三维固定约束,即对坐标数值大的支撑物施加拉力,本发明对参考腔的约束方式与实际的参考腔受力情况更加接近,能够更加准确的仿真出参考腔的形变,从而准确的评估加速度对参考腔形变的影响,即振动敏感度;
进一步的,本发明还提出前后两次施加不同主加速度仿真参考腔腔长变化,并对两次结果做差,消除预紧力导致的参考腔腔长变化,进而得出参考腔腔长对主加速度的敏感度,即得到所谓的参考腔振动敏感度。本发明能够更为准确的仿真计算抗振型参考腔振动敏感度,抗振型超稳光学参考腔的设计,节省抗振型超稳参考腔研制的经济成本与时间成本,同时本发明也可以应用于光学干涉仪、超窄线宽激光器、光学原子钟、引力波探测等高科技领域;
进一步的,本发明通过腔体结构的切割尺寸,支撑物和腔体支撑面的接触形状和接触面面积,支撑物厚度和支撑物材料属性进行参数化设置,生成不同参数情况,经过计算不同参数情况下的参考腔振动敏感度,选出最佳的振动敏感度为参考腔的最佳设计参数,通过该设计方法,本发明设计的参考腔充分考虑了各个参数化变量的情况,具备最佳的设计参数,整体结构更加稳固,性能更好。
附图说明
图1为本发明某一实施例的立方体光学参考腔结构示意图;
图2为本发明某一实施例的梭形光学参考腔结构示意图。
附图中:1-第一支撑物,2-第二支撑物,3-第三支撑物,4-第四支撑物,5-补偿环,6-腔镜,7-围脖。
具体实施方式
下面结合附图和具体实施方式对本发明作进一步的说明。
本发明使用有限元软件分析抗振型光学参考腔振动敏感度的方法并基于该方法选出参考腔最佳的设计参数,具体步骤如下:步骤一,根据拟分析的参考腔建立分析建模,一般支撑物的尺寸较小,建立模型时可以考虑忽略支撑物,而在腔体支撑面上根据支撑物与腔体的接触面积截取一定的面积作为支撑约束面,也可以将支撑物作为分析模型的一部分,进行多体仿真分析。
步骤二,建立适当的约束条件和外部载荷;对于抗振型超稳光学参考腔,为实现三维约束,参考腔的支撑面通常大于等于三个,对参考腔施加外部载荷-主加速度,为了准确体现参考腔的受力情况,应根据主加速度方向,以第一主加速度方向为一维坐标轴,坐标数值大的支撑物进行三维固定约束,坐标数值小的支撑物施加挤压力约束,或以第一主加速度方向为一维坐标轴,坐标数值大的支撑物仅在垂直于支撑面方向进行一维固定约束。
步骤三,参考腔结构参数化设置;参数化变量包括以下至少一项,腔体结构的切割尺寸(如沿立方腔体对角线切割顶角的深度、参考腔围脖的厚度及内径和外径尺寸、腔镜的直径、腔镜的厚度、参考腔补偿环厚度及内径和外径)、支撑物和腔体支撑面的接触形状和接触面面积,支撑物厚度和支撑物材料属性,仿真分析每个变量对参考腔长度变化的影响。
步骤四,模型网格划分;根据建立参考腔模型的尺寸对参考腔模型进行单元格划分,便于得出参考腔的形变量;在载荷施加面的支撑物、腔镜关键部分进行更加精细的网格细化;不同网格的尺寸需要进行均匀过渡设置,即网格的边长逐渐增加。
步骤五,算例运行,提取计算结果;仿真计算各种参数情况下参考腔腔长形变,即参考腔长度变化;两个腔镜内表面中心位置的距离表示参考腔长度,提取两个腔镜内表面中心沿光轴方向位移变化,以其相对变化代表此条件下参考腔长度变化。
步骤六,重新建立与步骤一完全相同的参考腔模型,在重新建立的参考腔模型的同一位置施加与步骤一载荷数值不同但方向相同的主加速度,重新建立的参考腔模型上支撑物的约束方式与步骤二的参考腔模型上支撑物的约束方式相同,继续进行仿真计算,计算本步骤下的参考腔腔长形变,通过步骤五的参考腔腔长形变和步骤六的参考腔腔长形变进行运算,得出拟分析的参考腔长度对主加速度的敏感程度,即振动敏感度,拟合出参考腔振动敏感度随各个参数化变量变化的曲线,根据所述曲线中振动敏感度最接近零,最好为零的数值选出对应的各个参数化变量的数值为参考腔最佳设计参数。
实施例1
以立方体光学参考腔为例,包括一个超低热膨胀玻璃(ULE)参考腔腔体、一对(或两对)熔融石英(FS)超高反射率反射镜、一对(或两对)ULE补偿环5、两个(或一个)通气孔和四个圆柱形支撑物;
如图1所示,本实施例中的立方体光学参考腔包括一个超低热膨胀玻璃(ULE)参考腔腔体、一对熔融石英(FS)超高反射率反射镜、一对ULE补偿环5、两个通气孔和四个圆柱形支撑物;
当光学参考腔水平放置,即主加速度方向与参考腔光轴垂直,光轴沿Y轴方向,主加速度沿负Z轴方向(或沿正Z轴方向,亦可沿负X轴方向或沿正X轴方向):
步骤一,根据拟分析的立体参考腔建立分析建模;设置参考腔每个部分的材料属性及每个部分间的接触方式。
步骤二,沿负Z轴方向施加1g主加速度,并将第三支撑物3的外端面和第四支撑物4的外端面进行三维固定约束,同时对第一支撑物1的外端面和第二支撑物2外端面施加法向挤压力F,
步骤三,对参考腔体的八个顶角的切割深度、四个支撑物和腔体支撑面的接触形状和接触面面积,支撑物厚度和支撑物材料属性,进行参数化设置。
步骤四,在计算资源允许的条件下,将模型进行精细化网格划分,特别是四个支撑物和两个腔镜6部分,需要进行更加精细的网格划分,同时不同网格尺寸间需要进行均匀过渡设置。
步骤五,算例运行,仿真计算各种参数情况下参考腔形变,提取两个腔镜6内表面中心沿光轴方向位移变化,以其相对变化代表此条件下参考腔长度变化;提取两个腔镜6内表面中心处单元节点沿Y轴方向的位置变化,并将两位移相减,得到各个参数化变量条件下对应的各个参考腔腔长的变化量ΔL1。
步骤六,重新建立与步骤一完全相同的参考腔模型,只将步骤二中主加速度载荷数值由1g改为2g,其他设置不变,重复以上步骤进行仿真计算,得到2g主加速度情况时的各个参数化变量条件下对应的各个参考腔腔长的变化量ΔL2,为消除预紧力导致的参考腔长度变化,计算各个参数化变量下参考腔长度对主加速度的敏感程度(ΔL2-ΔL1)/L/(a2-a1)即振动敏感度,单位为g-1,式中:L为第一参考腔模型或第二参考腔模型未发生形变的参考腔腔长,a1为第一主加速度,a2为第二主加速度,然后整理每个参数化变量下的参考腔振动敏感度,拟合出参考腔振动敏感度随参数化变量的曲线,根据曲线中振动敏感度最接近零或者为零的数值选出对应的参数化变量的数值为参考腔最佳设计参数。
实施例2
如图1所示,以立方体光学参考腔为例;
当光学参考腔竖直放置,即主加速度方向与参考腔光轴平行,光轴沿Y轴方向,加速度沿负Y轴方向(或沿正Y轴方向):
步骤一,根据图1所示参考腔的建立分析建模;设置参考腔各部分的材料属性及各部分间的接触方式。
步骤二,沿负Y轴方向施加1g主加速度,并将第一支撑物1的外端面和第三支撑物3的外端面进行三维固定约束,同时对第二支撑物2的外端面和第四支撑物4外端面施加法向挤压力F。
步骤三,参考腔体的八个顶角的切割深度、四个支撑物和腔体支撑面的接触形状和接触面面积,支撑物厚度和支撑物材料属性,进行参数化设置。
步骤四,在计算资源允许的条件下,将模型进行精细化网格划分,特别是四个支撑物和两个腔镜6部分,需要进行更加精细的网格划分,同时不同网格尺寸间需要进行均匀过渡设置。
步骤五,算例运行,仿真计算各种参数情况下参考腔形变,提取两个腔镜6内表面中心沿光轴方向位移变化,以其相对变化代表此条件下参考腔长度变化;提取两个腔镜6内表面中心处单元节点沿Y轴方向的位置变化,并将两位移相减,得到各个参数化变量条件下对应的各个参考腔腔长的变化量ΔL1。
步骤六,重新建立于步骤一完全相同的参考腔模型,只将步骤二中主加速度载荷数值由1g改为2g,其他设置不变,重复以上步骤进行仿真计算,得到2g主加速度情况时各个参数化变量条件下对应的各个参考腔腔长的变化量ΔL2,为消除预紧力导致的参考腔长度变化,计算各个参数化变量下参考腔长度对主加速度的敏感程度(ΔL2-ΔL1)/L/(a2-a1)即振动敏感度,单位为g-1,式中:L为第一参考腔模型或第二参考腔模型未发生形变的参考腔腔长,a1为第一主加速度,a2为第二主加速度,然后整理每个参数化变量下的参考腔振动敏感度,拟合出参考腔振动敏感度随参数化变量的曲线,根据曲线中振动敏感度最接近零或者为零的数值选出对应的参数化变量的数值为参考腔最佳设计参数。
实施例3
如图2所示,以竖直放置带围脖7的梭形光学参考腔为例(在本实施例中,梭形光学参考腔的围脖上有六个用于安装固定螺栓的通孔,共计12个缓冲支撑物);
当光学参考腔竖直放置,即主加速度方向与参考腔光轴平行,光轴沿Z轴方向,主加速度沿负Z轴方向:
步骤一,根据图2所示参考腔的建立分析建模;设置参考腔各部分的材料属性及各部分间的接触方式。
步骤二,沿负Z轴方向施加1g主加速度,并将围脖7下方六个支撑物的外端面进行沿Z轴的一维固定约束,同时将围脖7上方六个支撑物的外端面施加法向挤压力F。
步骤三,将围脖7厚度、内径、外径、腔镜6的直径、腔镜6的厚度、参考腔补偿环厚度及内外径,梭形腔体的上、下端面直径和十二个支撑物和腔体支撑面的接触形状和接触面面积,支撑物厚度和支撑物材料属性进行参数化设置。
步骤四,在计算资源允许的条件下,将模型进行精细化网格划分,特别是支撑物和两个腔镜6部分,需要进行更加精细的网格划分,同时不同网格尺寸间需要进行均匀过渡设置,即单元格的尺寸逐渐递增。
步骤五,算例运行,仿真计算各种参数情况下参考腔形变,提取两个腔镜内表面中心沿光轴方向位移变化,以其相对变化代表此条件下参考腔长度变化;提取两个腔镜内表面中心处单元节点沿Y轴方向的位置变化,并将两位移变化相减,得到各个参数化变量条件下对应的各个参考腔腔长的变化量ΔL1。
步骤六,将步骤二中主加速度载荷数值由1g改为2g,其他设置不变,重复以上步骤进行仿真计算,得到2g主加速度情况时的各个参数化变量条件下对应的各个参考腔腔长的变化量ΔL2,为消除预紧力导致的参考腔长度变化,计算各个参数化变量下参考腔长度对主加速度的敏感程度(ΔL2-ΔL1)/L/(a2-a1)即振动敏感度,单位为g-1,式中:L为第一参考腔模型或第二参考腔模型未发生形变的参考腔腔长,a1为第一主加速度,a2为第二主加速度,然后整理每个参数化变量下的参考腔振动敏感度,拟合出参考腔振动敏感度随参数化变量的曲线,根据曲线中振动敏感度最接近零或者为零的数值选出对应的参数化变量的数值为参考腔最佳设计参数。
在本发明的某一实施例中,梭形光学参考腔的围脖上有三个用于安装固定螺栓的通孔,共计6个缓冲支撑物;并且与实施例3中六个用于安装固定螺栓的通孔,共计12个缓冲支撑物的分析及设计方法一致;
在本发明中,带围脖的圆柱体的分析及设计方法与梭形光学参考腔一致。
参考腔在支撑物的挤压下,才能实现整个参考腔的三维约束,本发明提供了一种在有限元软件上对抗振型光学参考腔振动敏感度分析方法,改变了现有技术中将所有支撑面都固定的约束方式,此种固定约束方式,在主加速度的作用下,被固定约束的支撑物对参考腔的作用表现为拉力,而不是实际中的挤压力,本发明充分考虑了参考腔实际的受力情况,对一部分支撑面施加挤压力约束,同时考虑到参考腔在主加速度作用下,会有沿主加速度运动的趋势,为了防止参考腔运动,以第一主加速度方向为一维坐标轴,坐标数值大的支撑物进行了三维固定约束,即对坐标数值大的支撑物施加拉力,本发明对参考腔的约束方式与实际的参考腔受力情况更加接近,能够更加准确的仿真出参考腔的形变,从而准确的评估加速度对参考腔形变的影响,即振动敏感度;
最后应说明的是:以上所述实施例,仅为本发明的具体实施方式,用以说明本发明的技术方案,而非对其限制,本发明的保护范围并不局限于此,尽管参照前述实施例对本发明进行了详细的说明,本领域的普通技术人员应当理解:任何熟悉本技术领域的技术人员在本发明揭露的技术范围内,其依然可以对前述实施例所记载的技术方案进行修改或可轻易想到变化,或者对其中部分技术特征进行等同替换;而这些修改、变化或者替换,并不使相应技术方案的本质脱离本发明实施例技术方案的精神和范围,都应涵盖在本发明的保护范围之内。因此,本发明的保护范围应所述以权利要求的保护范围为准。
Claims (7)
1.抗振型光学参考腔振动敏感度分析方法,其特征在于,首先,对第一参考腔模型施加第一主加速度,以第一主加速度方向为一维坐标轴,坐标数值大的支撑物进行三维固定约束,坐标数值小的支撑物施加挤压力约束,仿真计算得到第一参考腔腔长形变量ΔL1;
然后,建立第二参考腔模型,所述第二参考腔模型与第一参考腔模型完全相同,对第二参考腔模型施加第二主加速度,所述第二主加速度的载荷值与第一主加速度的载荷值不同但方向相同,第二参考腔模型上支撑物的约束方式与所述第一参考腔模型上支撑物约束方式相同,仿真计算得到第二参考腔腔长形变量ΔL2;
最后,计算振动敏感度,单位为g-1:
(ΔL2-ΔL1)/L/(a2-a1),
式中:L为第一参考腔模型或第二参考腔模型未发生形变的参考腔腔长,a1为第一主加速度,a2为第二主加速度;
1)根据拟分析的参考腔建立分析建模,作为第一参考腔模型;
2)对第一参考腔模型施加第一主加速度,以第一主加速度方向为一维坐标轴,坐标数值大的支撑物进行三维固定约束,坐标数值小的支撑物施加挤压力约束,或以第一主加速度方向为一维坐标轴,坐标数值大的支撑物仅在垂直于支撑面方向进行一维固定约束;
3)根据第一参考腔模型的尺寸对第一参考腔模型进行单元格划分;
4)通过步骤3)单元格划分后的第一参考腔模型,计算第一参考腔腔长形变量ΔL1;
5)建立与步骤1)完全相同的参考腔模型为第二参考腔模型;
6)对第二参考腔模型施加第二主加速度,所述第二主加速度的载荷值与第一主加速度的载荷值不同但方向相同,第二参考腔模型上支撑物的约束方式与第一参考腔模型上支撑物约束方式相同;
7)根据第二参考腔模型的尺寸对第二参考腔模型进行单元格划分;
8)通过步骤7)单元格划分后的第二参考腔模型,计算得出第二参考腔腔长形变量ΔL2;
计算拟分析的参考腔振动敏感度,单位为g-1:
(ΔL2-ΔL1)/L/(a2-a1),
式中:L为第一参考腔模型或第二参考腔模型未发生形变的参考腔腔长,a1为第一主加速度,a2为第二主加速度。
2.根据权利要求1所述的抗振型光学参考腔振动敏感度分析方法,其特征在于,所述步骤3)中对模型进行单元格划分时,对支撑物和腔镜进行网格细化。
3.根据权利要求1所述的抗振型光学参考腔振动敏感度分析方法,其特征在于,所述步骤3)和步骤7)中所述单元格的尺寸进行均匀过渡设置。
4.根据权利要求1所述的抗振型光学参考腔振动敏感度分析方法,其特征在于,所述挤压力包括对所述支撑物的外端面的法向施加挤压力。
5.根据权利要求1-4任一项所述的抗振型光学参考腔振动敏感度分析方法的设计参考腔的方法,其特征在于,计算若干种参数化变量下的参考腔振动敏感度;
拟合出参考腔振动敏感度随各个参数化变量变化的曲线,根据所述曲线中振动敏感度最接近零或者为零的数值选出对应的各个参数化变量的数值为参考腔最佳设计参数。
6.根据权利要求5所述的抗振型光学参考腔振动敏感度的设计参考腔的方法,其特征在于,所述参数化变量包括参考腔体的切割尺寸,支撑物和支撑面的接触形状和接触面的面积,支撑物厚度以及支撑物材料属性中至少一项。
7.根据权利要求6所述的抗振型光学参考腔振动敏感度的设计参考腔的方法,其特征在于,所述参考腔体的切割尺寸包括沿立方腔体对角线切割顶角的深度或参考腔围脖的厚度及内径和外径尺寸、腔镜的直径、腔镜的厚度、参考腔补偿环厚度及内径和外径。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010103216.6A CN111324933B (zh) | 2020-02-19 | 2020-02-19 | 抗振型光学参考腔振动敏感度分析及设计参考腔的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010103216.6A CN111324933B (zh) | 2020-02-19 | 2020-02-19 | 抗振型光学参考腔振动敏感度分析及设计参考腔的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111324933A CN111324933A (zh) | 2020-06-23 |
CN111324933B true CN111324933B (zh) | 2023-03-21 |
Family
ID=71163757
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010103216.6A Active CN111324933B (zh) | 2020-02-19 | 2020-02-19 | 抗振型光学参考腔振动敏感度分析及设计参考腔的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111324933B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111912608B (zh) * | 2020-06-24 | 2022-01-18 | 中国科学院国家授时中心 | 一种可搬运光学参考腔振动敏感度的测试方法及装置 |
CN114122888A (zh) * | 2021-10-22 | 2022-03-01 | 中国科学院国家授时中心 | 一种用于光钟的频率可调谐、可搬运稳频激光系统 |
CN114740581B (zh) * | 2022-04-12 | 2023-03-17 | 中国科学院国家授时中心 | 一种光学参考腔零膨胀温度点调节装置 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104361631A (zh) * | 2014-10-24 | 2015-02-18 | 哈尔滨工业大学 | 参数化建立季节冻土区铁路路基振动反应三维数值模型的方法 |
CN104898403A (zh) * | 2015-05-19 | 2015-09-09 | 中国科学院国家授时中心 | 一种空间应用暨可搬运超稳光学参考腔 |
CN107862119A (zh) * | 2017-10-31 | 2018-03-30 | 无锡市五十五度科技有限公司 | 一种振动可靠性测试方法 |
CN109828342A (zh) * | 2019-02-28 | 2019-05-31 | 中国科学院国家授时中心 | 一种多功能超稳光学参考腔及其安装方法 |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101510229B (zh) * | 2009-03-20 | 2011-09-21 | 西安电子科技大学 | 基于机电热三场耦合的电子设备机箱结构优化设计方法 |
CN104008250B (zh) * | 2014-06-10 | 2016-01-20 | 广东工业大学 | 基于主频能量时域最优分布的非对称变加速度规划方法 |
US10931077B2 (en) * | 2018-02-13 | 2021-02-23 | Government Of The United States Of America, As Represented By The Secretary Of Commerce | Optical reference cavity |
-
2020
- 2020-02-19 CN CN202010103216.6A patent/CN111324933B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104361631A (zh) * | 2014-10-24 | 2015-02-18 | 哈尔滨工业大学 | 参数化建立季节冻土区铁路路基振动反应三维数值模型的方法 |
CN104898403A (zh) * | 2015-05-19 | 2015-09-09 | 中国科学院国家授时中心 | 一种空间应用暨可搬运超稳光学参考腔 |
CN107862119A (zh) * | 2017-10-31 | 2018-03-30 | 无锡市五十五度科技有限公司 | 一种振动可靠性测试方法 |
CN109828342A (zh) * | 2019-02-28 | 2019-05-31 | 中国科学院国家授时中心 | 一种多功能超稳光学参考腔及其安装方法 |
Non-Patent Citations (2)
Title |
---|
叶森 ; 龙云 ; 王如泉 ; 曹强 ; 肖三强 ; 赵宏鸣 ; 杜渐 ; .一种优化超稳光学腔的数值方法.2016,(06),全文. * |
郭海凤 ; 吕莎莎 ; 刘芳 ; 王春 ; 沈辉 ; 李刘锋 ; 陈李生 ; .法布里珀罗腔振动敏感度的分析.2012,(12),全文. * |
Also Published As
Publication number | Publication date |
---|---|
CN111324933A (zh) | 2020-06-23 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111324933B (zh) | 抗振型光学参考腔振动敏感度分析及设计参考腔的方法 | |
CN114364508A (zh) | 高刚度重量比反射性光学器件的3d打印 | |
Hecht et al. | Three-dimensional simulations and analysis of the nonlinear stage of the Rayleigh-Taylor instability | |
Zhikharev | A Sierpiński triangle geometric algorithm for generating stronger structures | |
Atkins et al. | Lightweighting design optimisation for additively manufactured mirrors | |
Wang et al. | Frequency splitting of hemispherical resonators trimmed with focused ion beams | |
Yong et al. | Numerical simulations of instabilities in the implosion process of inertial confined fusion in 2D cylindrical coordinates | |
Olhoff | On optimum design of structures and materials | |
Xie et al. | The influence of finite element meshing accuracy on a welding machine for offshore platform’s modal analysis | |
Wang et al. | An investigation on the precision mounting process of large-aperture potassium dihydrogen phosphate converters based on the accurate prediction model | |
Yuan et al. | A New Approach to Nonlinear Dynamic Modeling and Vibration Analysis of Tensegrity Structures | |
Pollard et al. | The structural analysis of a light-weight aluminum foam core mirror | |
Aghaie et al. | Birefringence analysis of photonic-bandgap fibers using the hexagonal Yee's cell | |
Smekalin et al. | Finite element analysis of the structural mechanics of the telescope objective for the lira-B space experiment | |
CN112926220A (zh) | 一种基于模态频率约束的结构件制备方法 | |
Guyomarc’h et al. | Some aspects of simulation and realization of an optical reference cavity | |
Roche et al. | Opto-mechanical analyses for performance optimization of lightweight grazing-incidence mirrors | |
Agresti et al. | Design and construction of a prototype of a flat top beam interferometer and initial tests. | |
Catalucci et al. | Point cloud processing techniques and image analysis comparisons for boat shapes measurements | |
Qian et al. | Modeling and Finite Element Analysis of Frame Type High Precision Adjustment Structure | |
Jedamzik et al. | Design and fabrication of a 3m class light weighted mirror blank for the E-ELT M5 | |
CN118049979A (zh) | 一种半球谐振陀螺的谐振子去重调平和平面度检测方法 | |
Wang et al. | Design and Analysis of Antenna Tower Tooling with High Stiffness for Complex Vibration Environments | |
CN108549772B (zh) | 一种结构优化的偏振片轮及其设计方法 | |
Canestrari et al. | Lightweight optical segment prototype for adaptive optics manufactured by hot slumping |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |