Wind disturbance simulation and optical performance prediction for large ground-based optical telescopes based on two-dimensional stochastic fields

基于二维随机场的地基大型光学望远镜风扰动时程模拟与性能预测 Wind disturbance simulation and optical performance prediction for large ground-based optical telescopes based on two-dimensional stochastic fields 1 2 corresp first-author E-mail： yuyan_cao@126.com E-mail： yuyan_cao@126.com 曹玉岩（1986-），男，吉林大安人，博士研究生，副研究员，2012年于西安电子科技大学获得硕士学位，主要从事地基大型光学望远镜系统集成仿真模拟方面的研究。E-mail：yuyan_cao@126.com http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=24671382&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=24671391&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=24671386&type= 22.01333237 32.00400543 曹玉岩 （1986-），男，吉林大安人，博士研究生，副研究员，2012年于西安电子科技大学获得硕士学位，主要从事地基大型光学望远镜系统集成仿真模拟方面的研究。E-mail： yuyan_cao@126.com yuyan_cao@126.com 1 2 王建立（1971-），男，山东曲阜人，研究员，博士生导师，主要从事空间探测技术和地基高分辨率成像光电望远镜总体技术的研究。E-mail：wangjianli@ciomp.ac.cn http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=24671395&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=24671401&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=24671398&type= 22.01333618 32.00399017 王建立 （1971-），男，山东曲阜人，研究员，博士生导师，主要从事空间探测技术和地基高分辨率成像光电望远镜总体技术的研究。E-mail： wangjianli@ciomp.ac.cn wangjianli@ciomp.ac.cn 1 2 1 1 2 1 1 中国科学院 长春光学精密机械与物理研究所，吉林 长春 130033 中国科学院 长春光学精密机械与物理研究所，吉林 长春 130033 Changchun Institute of Optics， Fine Mechanics and Physics， Chinese Academy of Sciences， Changchun 130033， China Changchun Institute of Optics， Fine Mechanics and Physics， Chinese Academy of Sciences， Changchun 130033， China 中国科学院大学，北京 100049 中国科学院大学，北京 100049 University of Chinese Academy of Sciences， Beijing 100049， China University of Chinese Academy of Sciences， Beijing 100049， China 随着地基光学望远镜口径的不断增大，风扰动已成为影响望远镜成像质量最为关键的因素之一。为了深入研究风扰动对望远镜系统性能的影响规律及作用机理，从时域角度对望远镜在风扰动作用下的响应进行时程模拟以及系统性能预测。首先，简要介绍了望远镜结构组成并采用有限元方法建立了结构动力学模型，通过模态变换方法将动力学模型转换到模态坐标系下，从而降低了模型维数、提高了计算效率。然后，提出了一种基于二维随机场的风速时程模拟方法，将望远镜圆顶内的风速场表达为随时间和空间位置变化的二维随机场，通过引入波数谱，克服了谱表达方法中各个离散采样点处互功率谱矩阵Cholesky分解出现数值不稳定的问题，从而可以表达空间和时间频率范围内几乎连续的随机场，并且在数值计算中引入了快速傅里叶变换FFT算法，进一步提高了模拟效率。最后，以2 m口径望远镜为例，对风速环境进行了时程模拟，并对外界平均风速为10 m/s及15 m/s的风扰动作用下，望远镜的性能进行了预测。仿真分析结果表明：作用在望远镜主镜上的风扰动分别造成主镜最大接近45 nm和70 nm的面形误差；作用在次镜及桁架上，主要造成低频的主次镜相对位置和角度偏差。 As the aperture of large ground-based optical telescopes increases， wind disturbance has become one of the most significant dynamic factors degrading the performance of telescopes. To investigate the influence and interaction principles of wind disturbance on the performance of telescopes， a time-history simulation of the wind disturbance response and performance prediction for telescopes were performed in the time domain. First， the structural parts of the telescope were simply introduced， and a dynamic model of the telescope was established using the finite element method. Using modal transformation， the dynamic model， expressed in nodal coordinates， was transformed into modal coordinates， thus greatly decreasing the model dimension and greatly improving the computational efficiency. Second， a time-history simulation method for wind speed was presented based on two-dimensional stochastic fields. The wind speed field in the doom was expressed as a two-dimensional stochastic field that varied according to the temporal and spatial frequencies. The numerical instability occurring in the Cholesky decomposition of the cross-power spectrum matrix in the spectral method was avoided by introducing the wave-number spectrum. The simulation efficiency was further improved by applying the fast Fourier transform （FFT） technique. Finally， using a ground-based telescope with a 2 m aperture as a case study， a time-history simulation of the wind speed， wind disturbance response， and system performance prediction for the telescope was performed. The simulation results showed that wind disturbances with a mean speed of 10 and 15 m/s would result in a maximum surface root mean square error of the primary mirror of approximately 45 and 70 nm， respectively. In addition， the wind disturbance acting on the secondary mirror and truss would mainly cause optical axis angle and position errors. 地基望远镜 随机场 风扰动 时程模拟 动力学模型 ground-based telescope stochastic fields wind disturbance time-history simulation dynamic structural model 1 引　言 1 地基望远镜是连接人类与未知宇宙的桥梁，帮助人们观测、探索宇宙奥秘。随着对宇宙暗弱目标研究的深入，对如深空恒星的形成与演变、星系分布与结构及演变物理过程对人类的影响、暗物质与暗能量的本质等等的探索，要求未来的望远镜不但要看的清、还要看得远，这推动了地基望远镜朝大口径、高分辨力方向发展 ［ 1 1 ］ 。 随着地基望远镜口径不断增大，风扰动成为影响望远镜成像质量最为关键的因素之一。首先，望远镜台址通常选在高海拔偏僻地区来避免杂光污染并提高大气视宁度，这些地区通常风速都很大（如夏威夷莫纳克亚山天文台风速高达20 m/s）。其次，望远镜结构刚度随口径的增大而下降，系统谐振频率随之降低。此外，风扰动功率谱中包含有较大的低频能量（0.1~1 Hz），在望远镜基频处风压能量集中，容易引起望远镜系统的共振，而且风载具有随机和时变特性 ［ 2 2 ］ 。一些望远镜往往因为风大而难以获取有效科学数据，甚至不得不放弃观测。风扰动的分析与抑制已成为地基望远镜领域的热点和难点，受到了越来越多的科学家们的关注 ［ 3 3 - 10 10 3-10 ］ 。 尽管光学望远镜外部一般有圆顶防护，能够显著降低其内部风速，但为了满足圆顶内视宁度的要求，内部要适当通风，而且光学望远镜波长较短，圆顶内的残余风仍然会对望远镜系统性能产生显著影响 ［ 3 3 ］ 。针对这一问题，MacMynowski等基于Gemini望远镜镜面风速和风压测量 ［ 5 5 ］ 、风道测试 ［ 6 6 ］ 和计算流体动力学仿真CFD ［ 7 7 ］ 得到的数据对望远镜圆顶内部的风扰动特性进行了深入研究，并提出了参数化模型 ［ 8 8 ， 9 9 ］ ，该参数化模型可根据望远镜、风速及圆顶等相关参数将风扰动转化为作用在望远镜主、次镜上载荷，给结构设计优化带来极大的方便。在国内，杨德华等 ［ 11 11 ］ 针对LAMOST反射施密特改正镜，研究了其拼接子镜及整体在风扰动作用下的静态和动态响应，采用风速功率谱密度法分析了其对镜面面形和跟踪指向精度的影响；周超等 ［ 12 12 ］ 以长春光机所研制的1.23 m自适应光学望远镜为例，研究了风扰动对主镜面形的影响及望远镜在动态风扰作用下的随机响应；潘年等 ［ 13 13 ］ 采用流体动力学仿真方法CFD研究了国内某2 m望远镜风扰特性，分析了外界风速10 m/s情况下望远镜不同俯仰角时流场中截面空气速度、压力、湍流动能以及主镜面静压力的瞬态分布情况，并与Gemini望远镜实测结果进行了比较。徐江海等 ［ 14 14 ］ 采用CFD分析圆顶形式对内部风湍流的影响。以上研究从不同的角度分析了外界风扰动对望远镜性能产生的影响。 从上述文献中可以看出，望远镜口径不断增大，系统刚度和谐振频率下降，风扰动的影响应在望远镜建造阶段予以重视。风扰动包含稳态分量和非稳态分量两部分，稳态分量（静态风载荷）是恒定的，通过系统主动控制很容易得到抑制，但抑制非稳态分量（动态风扰动）引起的振动则存在一定的困难。针对此问题，已有多位学者采用计算机仿真手段深入分析了非稳态风扰动对望远镜性能的影响程度 ［ 8 8 - 9 9 8-9 ］ 。非稳态风扰动分量具有较强的随机性，目前已有的大部分研究均是从频率角度考虑风扰动的影响，时域角度的研究则很少。然而，对于大型地基望远镜这样复杂结构来说，其结构中往往存在一定的非线性特性，风扰动的时程分析模拟是必不可少的，而且从时域角度分析与工作过程更为接近，这也更利于进行实验比较。 基于上述原因，本文以长春光机所研制的2 m口径地基光学望远镜为例，从时域角度研究风扰动对望远镜的影响，并进行性能预测。首先，详细介绍了望远镜结构组成，基于有限元方法建立结构动力学模型，并通过模态变换方法将动力学模型转换到模态坐标系下，从而降低了模型维数，提高了计算效率。然后，考虑到大型地基望远镜结构尺度相对较大，将望远镜圆顶内的风速场表达为随时间和空间位置变化的二维随机场，通过引入互相关函数或互功率谱矩阵来刻画随机场的统计特征。基于互功率谱矩阵，采用谱表达方法实现时程模拟。为了克服在离散采样点数提高时，谱表达方法中各个离散采样点处互功率谱矩阵Cholesky分解出现数值不稳定的问题，引入了文献［ 15 15 - 17 17 15-17 ］所述波数谱方法，从而可以表达空间和时间频率范围内几乎连续的随机场。最后，对2 m望远镜进行风速场模拟并在此风速场下进行性能分析。 2 望远镜结构动力学建模 1 2 m级地基光学望远镜结构简图如 图1 图1 所示，总体高度为6.4 m，宽度为3.7 m，总质量接近30吨。从望远镜结构角度考虑，主要包括地平式机架、主镜和次镜三部分，其余附属结构对本文所研究内容影响较小，在此不予考虑。 自无穷远目标发出的光束经过主镜和次镜反射后，被像面位置处的高灵敏度探测器接收，实现光信号到电信号的转换，从而获得清晰的目标图像。地平式机架是主、次镜及其他附件的载体，通过俯仰轴和方位轴的高精度回转运动实现望远镜对任意天区的精确指向和对目标的稳定跟踪测量。 望远镜主镜材料为SiC，具有比刚度大、质量轻、热稳定性好等优点，采用A-Frame型柔性侧支撑方式，并结合18点机械式whiffletree轴向支撑方式，支撑结构原理如 图1 图1 所示。同样，次镜采用具有热变形匹配的bipod型结构来消除热膨胀系数差异性对次镜面形精度的影响。此外，次镜结构整体通过Hexapod并联机构实现位置调整与补偿。 针对上述2 m望远镜简化结构，采用有限元方法建立望远镜结构动力学模型如 图2 图2 所示。在全局坐标系下，望远镜结构动力学模型可以表达为二阶微分方程，即： <math display="block" id="M1"><mi mathvariant="bold-italic">M</mi><mover accent="true"><mrow><mi mathvariant="bold-italic">q</mi></mrow><mi mathvariant="bold-italic">¨</mi></mover><mo>+</mo><mi mathvariant="bold-italic">D</mi><mover accent="true"><mi mathvariant="bold-italic">q</mi><mo>˙</mo></mover><mo>+</mo><mi mathvariant="bold-italic">K</mi><mi mathvariant="bold-italic">q</mi><mo>=</mo><mi mathvariant="bold-italic">F</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=24671429&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=24671426&type= 32.08866501 3.72533321 ， <math display="block" id="M2"><mi mathvariant="bold-italic">y</mi><mo>=</mo><msub><mrow><mi mathvariant="bold-italic">C</mi></mrow><mrow><mi mathvariant="normal">o</mi></mrow></msub><mi mathvariant="bold-italic">q</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=24671436&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=24671433&type= 12.69999981 4.57200003 ， 其中： <math id="M3"><mi mathvariant="bold-italic">M</mi><mo>,</mo><mi mathvariant="bold-italic">D</mi><mo>,</mo><mi mathvariant="bold-italic">K</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=24671442&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=24671439&type= 12.78466606 3.21733332 分别为质量矩阵，阻尼矩阵和刚度矩阵； F 为载荷向量； y 为输出位移向量； <math id="M4"><mover accent="true"><mrow><mi mathvariant="bold-italic">q</mi></mrow><mi mathvariant="bold-italic">¨</mi></mover><mo>,</mo><mover accent="true"><mi mathvariant="bold-italic">q</mi><mo>˙</mo></mover><mo>,</mo><mi mathvariant="bold-italic">q</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=24671448&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=24671445&type= 9.56733322 3.64066648 分别为节点加速度，速度和位移向量； C o 为输出矩阵，在风扰动分析中，主要关注的是主、次镜节点的位移，因此输出矩阵为与主、次镜节点位置相关的矩阵。 由于望远镜结构非常大，有限元模型节点数量非常多，上述以节点变量来表达的结构动力学模型将非常大。为了提高计算效率，降低后续风扰动模拟的计算难度，对结构动力学模型进行模态变换，即 <math display="block" id="M5"><mi mathvariant="bold-italic">q</mi><mo>=</mo><mi mathvariant="bold-italic">Φ</mi><msub><mrow><mi mathvariant="bold-italic">q</mi></mrow><mrow><mi mathvariant="normal">m</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=24671456&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=24671451&type= 13.63133430 4.57200003 ， 其中， <math id="M6"><mi mathvariant="bold-italic">Φ</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=24671461&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=24671459&type= 3.13266683 2.79399991 为振型矩阵，其表达式为： <math display="block" id="M7"><mtable><mtr><mtd><mi mathvariant="bold-italic">Φ</mi><mo>=</mo><mo stretchy="false">[</mo><mtable><mtr><mtd><msub><mrow><mi>ϕ</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow></msub></mtd><mtd><msub><mrow><mi>ϕ</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msub></mtd><mtd><mo>⋯</mo></mtd><mtd><msub><mrow><mi>ϕ</mi></mrow><mrow><mi>n</mi></mrow></msub></mtd></mtr></mtable><mo stretchy="false">]</mo><mo>=</mo></mtd></mtr><mtr><mtd><mfenced open="[" close="]" separators="|"><mrow><mtable><mtr><mtd><msub><mrow><mi>ϕ</mi></mrow><mrow><mn mathvariant="normal">11</mn></mrow></msub></mtd><mtd><msub><mrow><mi>ϕ</mi></mrow><mrow><mn mathvariant="normal">12</mn></mrow></msub></mtd><mtd><mo>⋯</mo></mtd><mtd><msub><mrow><mi>ϕ</mi></mrow><mrow><mn mathvariant="normal">1</mn><mi>n</mi></mrow></msub></mtd></mtr><mtr><mtd><msub><mrow><mi>ϕ</mi></mrow><mrow><mn mathvariant="normal">21</mn></mrow></msub></mtd><mtd><msub><mrow><mi>ϕ</mi></mrow><mrow><mn mathvariant="normal">22</mn></mrow></msub></mtd><mtd><mo>⋯</mo></mtd><mtd><msub><mrow><mi>ϕ</mi></mrow><mrow><mn mathvariant="normal">2</mn><mi>n</mi></mrow></msub></mtd></mtr><mtr><mtd><mo>⋮</mo></mtd><mtd><mo>⋮</mo></mtd><mtd><mo>⋱</mo></mtd><mtd><mo>⋮</mo></mtd></mtr><mtr><mtd><msub><mrow><mi>ϕ</mi></mrow><mrow><msub><mrow><mi>n</mi></mrow><mrow><mi>d</mi></mrow></msub><mn mathvariant="normal">1</mn></mrow></msub></mtd><mtd><msub><mrow><mi>ϕ</mi></mrow><mrow><msub><mrow><mi>n</mi></mrow><mrow><mi>d</mi></mrow></msub><mn mathvariant="normal">2</mn></mrow></msub></mtd><mtd><mo>⋯</mo></mtd><mtd><msub><mrow><mi>ϕ</mi></mrow><mrow><msub><mrow><mi>n</mi></mrow><mrow><mi>d</mi></mrow></msub><mi>n</mi></mrow></msub></mtd></mtr></mtable></mrow></mfenced></mtd></mtr></mtable></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=24671468&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=24671463&type= 41.06333542 27.77066803 ， 其中： <math id="M8"><msub><mrow><mi>ϕ</mi></mrow><mrow><mi>i</mi><mi>j</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=24671474&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=24671471&type= 3.55599999 4.57200003 表示第 i 个节点自由度上的第 j 阶振型； n d 为模型阶数，即节点的全部自由度总和； n 为提取的模态阶数。 将 式（3） 式（3） 带入 式（1） 式（1） 和 式（2） 式（2） ，并在 式（1） 式（1） 左乘 <math id="M9"><msup><mrow><mi mathvariant="bold-italic">Φ</mi></mrow><mrow><mi mathvariant="normal">T</mi></mrow></msup></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=24671480&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=24671477&type= 4.74133301 3.04800010 得： <math display="block" id="M10"><msub><mrow><mi mathvariant="bold-italic">M</mi></mrow><mrow><mi mathvariant="normal">m</mi></mrow></msub><msub><mrow><mover accent="true"><mrow><mi mathvariant="bold-italic">q</mi></mrow><mi mathvariant="bold-italic">¨</mi></mover></mrow><mrow><mi mathvariant="normal">m</mi></mrow></msub><mo>+</mo><msub><mrow><mi mathvariant="bold-italic">D</mi></mrow><mrow><mi mathvariant="normal">m</mi></mrow></msub><msub><mrow><mover accent="true"><mi mathvariant="bold-italic">q</mi><mo>˙</mo></mover></mrow><mrow><mi mathvariant="normal">m</mi></mrow></msub><mo>+</mo><msub><mrow><mi mathvariant="bold-italic">K</mi></mrow><mrow><mi mathvariant="normal">m</mi></mrow></msub><msub><mrow><mi mathvariant="bold-italic">q</mi></mrow><mrow><mi mathvariant="normal">m</mi></mrow></msub><mo>=</mo><msup><mrow><mi mathvariant="bold-italic">Φ</mi></mrow><mrow><mi mathvariant="normal">T</mi></mrow></msup><mi mathvariant="bold-italic">F</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=24671486&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=24671483&type= 48.51399994 4.82600021 ， <math display="block" id="M11"><mi mathvariant="bold-italic">y</mi><mo>=</mo><msub><mrow><mi mathvariant="bold-italic">C</mi></mrow><mrow><mi mathvariant="normal">o</mi></mrow></msub><mi mathvariant="bold-italic">Φ</mi><msub><mrow><mi mathvariant="bold-italic">q</mi></mrow><mrow><mi mathvariant="normal">m</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=24671492&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=24671489&type= 17.61066628 4.57200003 ， 其中： <math id="M12"><msub><mrow><mi mathvariant="bold-italic">M</mi></mrow><mrow><mi mathvariant="normal">m</mi></mrow></msub><mo>=</mo><msup><mrow><mi mathvariant="bold-italic">Φ</mi></mrow><mrow><mi mathvariant="normal">T</mi></mrow></msup><mi mathvariant="bold-italic">M</mi><mi mathvariant="bold-italic">Φ</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=24671499&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=24671495&type= 21.84399986 4.06400013 ， <math id="M13"><msub><mrow><mi mathvariant="bold-italic">D</mi></mrow><mrow><mi mathvariant="normal">m</mi></mrow></msub><mo>=</mo><msup><mrow><mi mathvariant="bold-italic">Φ</mi></mrow><mrow><mi mathvariant="normal">T</mi></mrow></msup><mi mathvariant="bold-italic">D</mi><mi mathvariant="bold-italic">Φ</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=24671505&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=24671502&type= 20.57399940 3.97933316 ， <math id="M14"><msub><mrow><mi mathvariant="bold-italic">K</mi></mrow><mrow><mi mathvariant="normal">m</mi></mrow></msub><mo>=</mo><msup><mrow><mi mathvariant="bold-italic">Φ</mi></mrow><mrow><mi mathvariant="normal">T</mi></mrow></msup><mi mathvariant="bold-italic">K</mi><mi mathvariant="bold-italic">Φ</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=24671511&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=24671508&type= 20.40466690 3.97933316 。 对以模态变量表达的结构动力学方程（5）和（6）进一步简化可得： <math display="block" id="M15"><msub><mrow><mover accent="true"><mrow><mi mathvariant="bold-italic">q</mi></mrow><mi mathvariant="bold-italic">¨</mi></mover></mrow><mrow><mi mathvariant="normal">m</mi></mrow></msub><mo>+</mo><mn mathvariant="normal">2</mn><mi mathvariant="bold-italic">Z</mi><mi mathvariant="bold-italic">Ω</mi><msub><mrow><mover accent="true"><mi mathvariant="bold-italic">q</mi><mo>˙</mo></mover></mrow><mrow><mi mathvariant="normal">m</mi></mrow></msub><mo>+</mo><msup><mrow><mi mathvariant="bold-italic">Ω</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup><msub><mrow><mi mathvariant="bold-italic">q</mi></mrow><mrow><mi mathvariant="normal">m</mi></mrow></msub><mo>=</mo><msubsup><mrow><mi mathvariant="bold-italic">M</mi></mrow><mrow><mi mathvariant="normal">m</mi></mrow><mrow><mo>-</mo><mn mathvariant="normal">1</mn></mrow></msubsup><msup><mrow><mi mathvariant="bold-italic">Φ</mi></mrow><mrow><mi mathvariant="normal">T</mi></mrow></msup><mi mathvariant="bold-italic">F</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=24671517&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=24671514&type= 52.15466690 4.82600021 ， <math display="block" id="M16"><mi mathvariant="bold-italic">y</mi><mo>=</mo><msub><mrow><mi mathvariant="bold-italic">C</mi></mrow><mrow><mi mathvariant="normal">m</mi></mrow></msub><msub><mrow><mi mathvariant="bold-italic">q</mi></mrow><mrow><mi mathvariant="normal">m</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=24671525&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=24671522&type= 14.98600006 4.57200003 ， 其中： <math id="M17"><mi mathvariant="bold-italic">Z</mi><mo>=</mo><mi>d</mi><mi>i</mi><mi>a</mi><mi>g</mi><mo stretchy="false">(</mo><msub><mrow><mi>ζ</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow></msub><mo>,</mo><msub><mrow><mi>ζ</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msub><mo>,</mo><mo>⋯</mo><mo>,</mo><msub><mrow><mi>ζ</mi></mrow><mrow><mi>n</mi></mrow></msub><mo stretchy="false">)</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=24671530&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=24671528&type= 35.56000137 4.82600021 为阻尼系数矩阵， <math id="M18"><mi mathvariant="bold-italic">Ω</mi><mo>=</mo><mi>d</mi><mi>i</mi><mi>a</mi><mi>g</mi><mo stretchy="false">(</mo><msub><mrow><mi>ω</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow></msub><mo>,</mo><msub><mrow><mi>ω</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msub><mo>,</mo><mo>⋯</mo><mo>,</mo><msub><mrow><mi>ω</mi></mrow><mrow><mi>n</mi></mrow></msub><mo stretchy="false">)</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=24671534&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=24671531&type= 39.11600113 4.82600021 为模态频率矩阵， <math id="M19"><msub><mrow><mi mathvariant="bold-italic">C</mi></mrow><mrow><mi mathvariant="normal">m</mi></mrow></msub><mo>=</mo><msub><mrow><mi mathvariant="bold-italic">C</mi></mrow><mrow><mi mathvariant="normal">o</mi></mrow></msub><mi mathvariant="bold-italic">Φ</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=24671538&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=24671535&type= 16.25600052 3.72533321 。 为了应用线性系统相关理论，将 式（7） 式（7） 和（8）表达的结构动力学模型转化为状态空间表达形式。定义状态变量为： <math display="block" id="M20"><mi mathvariant="bold-italic">x</mi><mo>=</mo><mfenced open="[" close="]" separators="|"><mrow><mtable columnalign="left"><mtr><mtd><msub><mrow><mi>x</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow></msub></mtd></mtr><mtr><mtd><msub><mrow><mi>x</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msub></mtd></mtr></mtable></mrow></mfenced><mo>=</mo><mfenced open="[" close="]" separators="|"><mrow><mtable columnalign="left"><mtr><mtd><mi mathvariant="normal">Ω</mi><msub><mrow><mi>q</mi></mrow><mrow><mi>m</mi></mrow></msub></mtd></mtr><mtr><mtd><msub><mrow><mover accent="true"><mi>q</mi><mo>˙</mo></mover></mrow><mrow><mi>m</mi></mrow></msub></mtd></mtr></mtable></mrow></mfenced></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=24671542&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=24671540&type= 27.17799950 9.90600014 . 将 式（9） 式（9） 带入 式（7） 式（7） 和（8）可得： <math display="block" id="M21"><mover accent="true"><mi mathvariant="bold-italic">x</mi><mo>˙</mo></mover><mo>=</mo><mi mathvariant="bold-italic">A</mi><mi mathvariant="bold-italic">x</mi><mo>+</mo><mi mathvariant="bold-italic">B</mi><mi mathvariant="bold-italic">u</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=24671546&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=24671544&type= 21.25133324 2.96333337 ， <math display="block" id="M22"><mi mathvariant="bold-italic">y</mi><mo>=</mo><mi mathvariant="bold-italic">C</mi><mi mathvariant="bold-italic">x</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=24671550&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=24671548&type= 11.34533310 3.64066648 ， 式中： <math id="M23"><mi mathvariant="bold-italic">A</mi><mo>=</mo><mfenced open="[" close="]" separators="|"><mrow><mtable><mtr><mtd><mtext>    </mtext><mn mathvariant="normal">0</mn></mtd><mtd><mi mathvariant="normal">Ω</mi></mtd></mtr><mtr><mtd><mo>-</mo><mi mathvariant="normal">Ω</mi></mtd><mtd><mo>-</mo><mn mathvariant="normal">2</mn><mi>Z</mi><mi mathvariant="normal">Ω</mi></mtd></mtr></mtable></mrow></mfenced></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=24671554&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=24671552&type= 29.97200012 8.63599968 ， <math id="M24"><mi mathvariant="bold-italic">B</mi><mo>=</mo><mfenced open="[" close="]" separators="|"><mrow><mtable><mtr><mtd><mn mathvariant="normal">0</mn></mtd></mtr><mtr><mtd><msub><mrow><mi>B</mi></mrow><mrow><mi mathvariant="normal">m</mi></mrow></msub></mtd></mtr></mtable></mrow></mfenced></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=24671558&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=24671556&type= 14.64733315 8.55133343 ， <math id="M25"><mi mathvariant="bold-italic">C</mi><mo>=</mo><mo stretchy="false">[</mo><mtable><mtr><mtd><msub><mrow><mi>C</mi></mrow><mrow><mi mathvariant="normal">m</mi></mrow></msub><msup><mrow><mi mathvariant="normal">Ω</mi></mrow><mrow><mo>-</mo><mn mathvariant="normal">1</mn></mrow></msup></mtd><mtd><mn mathvariant="normal">0</mn></mtd></mtr></mtable><mo stretchy="false">]</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=24671562&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=24671560&type= 27.26266861 3.55599999 。 与 式（1） 式（1） 和（2）相比，状态空间方程（10）和（11）表达的结构动力学模型的维数大大降低，从而可以显著提高计算效率。 3 随机风扰动建模 1 对于地基大型光学望远镜而言，风扰动的直接影响是造成望远镜视轴抖动及光学元件（主要是主镜和次镜）的振动和变形，进而导致焦面处图像抖动和模糊。风扰动具有随机性，作用在望远镜结构上产生的变形非常复杂，即使采用计算流体动力学分析方法CFD也很难准确模拟 ［ 7 7 ］ 。为了简化计算模拟过程，本文将风扰动对望远镜的作用简化为如 图3 图3 所示的三种方式：（1）载荷作用在次镜及其支撑桁架上，引起次镜抖动；（2）载荷作用在主镜及其支撑结构上引起主镜变形和位置的变化；（3）载荷作用在次镜及支撑桁架上，造成俯仰轴的扰动，进而引起视轴抖动。风扰动的综合作用效果是引起主、次镜光轴相对位置出现偏差，镜面出现面形误差以及视轴抖动，进而造成光学系统性能下降。 3.1　基于二维随机场的风速时程模拟 2 从 图1 图1 已知望远镜的横向和纵向最大跨度相对较大，分别为6.4 m和3.7 m，对风速进行时程模拟需要同时考虑时间频率和空间频率，且空间频率需要考虑多个维度。这里提出了以平面内的二维随机场来模拟作用在望远镜结构上的风速，风速在空间上任何一点的时间历程上既具有随机性又符合统计规律，其功率谱与预先已知的功率谱密度一致。首先推导一维随机场，然后再将其拓展至二维情况。 定义 X （ t ）为离散且相互独立的随机过程向量，即： <math display="block" id="M26"><mi mathvariant="bold-italic">X</mi><mo stretchy="false">(</mo><mi>t</mi><mo stretchy="false">)</mo><mo>=</mo><mo stretchy="false">[</mo><msub><mrow><mi>X</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow></msub><mo stretchy="false">(</mo><mi>t</mi><mo stretchy="false">)</mo><mo>,</mo><mo>⋯</mo><mo>,</mo><msub><mrow><mi>X</mi></mrow><mrow><mi>j</mi></mrow></msub><mo stretchy="false">(</mo><mi>t</mi><mo stretchy="false">)</mo><mo>,</mo><mo>⋯</mo><mo>,</mo><msub><mrow><mi>X</mi></mrow><mrow><mi>n</mi></mrow></msub><msup><mrow><mo stretchy="false">(</mo><mi>t</mi><mo stretchy="false">)</mo><mo stretchy="false">]</mo></mrow><mrow><mi mathvariant="normal">T</mi></mrow></msup></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=24671574&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=24671571&type= 59.85933304 4.82600021 ， 其中， <math id="M27"><msub><mrow><mi>X</mi></mrow><mrow><mi>j</mi></mrow></msub><mo stretchy="false">(</mo><mi>t</mi><mo stretchy="false">)</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=24671580&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=24671578&type= 8.38199997 4.48733330 为第 j 个单变量随机过程。 随机过程向量中的各元素均为相互独立的单变量随机过程，但作为一个整体，不同元素之间存在一定的相关关系，可以用互相关函数矩阵或互功率谱密度矩阵来表示，即： <math display="block" id="M28"><msub><mrow><mi>S</mi></mrow><mrow><mi>j</mi><mi>k</mi></mrow></msub><mo stretchy="false">(</mo><mi>ω</mi><mo stretchy="false">)</mo><mo>=</mo><mfrac><mrow><mn mathvariant="normal">1</mn></mrow><mrow><mn mathvariant="normal">2</mn><mi mathvariant="normal">π</mi></mrow></mfrac><msubsup><mo>∫</mo><mrow><mo>-</mo><mi mathvariant="normal">∞</mi></mrow><mrow><mi mathvariant="normal">∞</mi></mrow></msubsup><mrow><msub><mrow><mi>R</mi></mrow><mrow><mi>j</mi><mi>k</mi></mrow></msub><mo stretchy="false">(</mo><mi>τ</mi><mo stretchy="false">)</mo><msup><mrow><mi mathvariant="normal">e</mi></mrow><mrow><mo>-</mo><mi>i</mi><mi>ω</mi><mi>τ</mi></mrow></msup><mi mathvariant="normal">d</mi><mi>τ</mi></mrow></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=24671584&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=24671582&type= 46.31266785 7.87400007 ， <math display="block" id="M29"><msub><mrow><mi>R</mi></mrow><mrow><mi>j</mi><mi>k</mi></mrow></msub><mo stretchy="false">(</mo><mi>τ</mi><mo stretchy="false">)</mo><mo>=</mo><msubsup><mo>∫</mo><mrow><mo>-</mo><mi mathvariant="normal">∞</mi></mrow><mrow><mi mathvariant="normal">∞</mi></mrow></msubsup><mrow><msub><mrow><mi>S</mi></mrow><mrow><mi>j</mi><mi>k</mi></mrow></msub><mo stretchy="false">(</mo><mi>ω</mi><mo stretchy="false">)</mo><msup><mrow><mi mathvariant="normal">e</mi></mrow><mrow><mi>i</mi><mi>ω</mi><mi>τ</mi></mrow></msup><mi mathvariant="normal">d</mi><mi>ω</mi></mrow></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=24671588&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=24671586&type= 40.21666718 7.87400007 . 为了从时程角度模拟随机过程，首先需要对其互功率谱密度矩阵进行分解，即： <math display="block" id="M30"><mi mathvariant="bold-italic">S</mi><mo stretchy="false">(</mo><mi>ω</mi><mo stretchy="false">)</mo><mo>=</mo><mi mathvariant="bold-italic">H</mi><mo stretchy="false">(</mo><mi>ω</mi><mo stretchy="false">)</mo><msup><mrow><mi mathvariant="bold-italic">H</mi></mrow><mrow><mi mathvariant="normal">T</mi><mi mathvariant="normal">*</mi></mrow></msup><mo stretchy="false">(</mo><mi>ω</mi><mo stretchy="false">)</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=24671593&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=24671591&type= 34.88266373 3.80999994 ， 其中： H （ ω ）为下三角矩阵，T*表示复共轭转置，上述分解可以通过Cholesky方法实现。 根据 式（15） 式（15） ，随机过程向量 X （ t ）的时间历程可以表达为： <math display="block" id="M31"><mtable columnalign="left"><mtr><mtd><msub><mrow><mi>X</mi></mrow><mrow><mi>j</mi></mrow></msub><mo stretchy="false">(</mo><mi>t</mi><mo stretchy="false">)</mo><mo>=</mo><mn mathvariant="normal">2</mn><mstyle displaystyle="true"><munderover><mo largeop="true">∑</mo><mrow><mi>m</mi><mo>=</mo><mn mathvariant="normal">1</mn></mrow><mrow><mi>j</mi></mrow></munderover></mstyle><mrow><mstyle displaystyle="true"><munderover><mo largeop="true">∑</mo><mrow><mi>l</mi><mo>=</mo><mn mathvariant="normal">0</mn></mrow><mrow><mi>N</mi><mo>-</mo><mn mathvariant="normal">1</mn></mrow></munderover></mstyle><mrow><mfenced open="|" close="|" separators="|"><mrow><msub><mrow><mi>H</mi></mrow><mrow><mi>j</mi><mi>m</mi></mrow></msub><mo stretchy="false">(</mo><msub><mrow><mi>ω</mi></mrow><mrow><mi>l</mi></mrow></msub><mo stretchy="false">)</mo></mrow></mfenced><mroot><mrow><mi mathvariant="normal">Δ</mi><mi>ω</mi></mrow><mrow/></mroot></mrow></mrow></mtd></mtr><mtr><mtd><mi mathvariant="normal">c</mi><mi mathvariant="normal">o</mi><mi mathvariant="normal">s</mi><mo stretchy="false">[</mo><msub><mrow><mi>ω</mi></mrow><mrow><mi>l</mi></mrow></msub><mi>t</mi><mo>+</mo><msub><mrow><mi>θ</mi></mrow><mrow><mi>j</mi><mi>m</mi></mrow></msub><mo stretchy="false">(</mo><msub><mrow><mi>ω</mi></mrow><mrow><mi>l</mi></mrow></msub><mo stretchy="false">)</mo><mo>+</mo><msub><mrow><mi>ϕ</mi></mrow><mrow><mi>m</mi><mi>l</mi></mrow></msub><mo stretchy="false">]</mo><mtext>   </mtext><mi>j</mi><mo>=</mo><mn mathvariant="normal">1,2</mn><mo>,</mo><mo>…</mo><mo>,</mo><mi>n</mi></mtd></mtr></mtable></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=24671596&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=24671595&type= 62.31467056 15.66333294 ， 其中： <math id="M32"><msub><mrow><mi>ω</mi></mrow><mrow><mi>l</mi></mrow></msub><mo>=</mo><mi>l</mi><mi mathvariant="normal">Δ</mi><mi>ω</mi><mo>,</mo><mi>l</mi><mo>=</mo><mn mathvariant="normal">0,1</mn><mo>,</mo><mo>…</mo><mo>,</mo><mi>N</mi><mo>-</mo><mn mathvariant="normal">1</mn></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=24671600&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=24671598&type= 44.78866577 4.14866638 ； <math id="M33"><mi mathvariant="normal">Δ</mi><mi>ω</mi><mo>=</mo><mfrac><mrow><msub><mrow><mi>ω</mi></mrow><mrow><mi mathvariant="normal">u</mi></mrow></msub></mrow><mrow><mi>N</mi></mrow></mfrac></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=24671604&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=24671602&type= 14.22399998 8.04333401 ； <math id="M34"><msub><mrow><mi>θ</mi></mrow><mrow><mi>j</mi><mi>m</mi></mrow></msub><mo stretchy="false">(</mo><msub><mrow><mi>ω</mi></mrow><mrow><mi>l</mi></mrow></msub><mo stretchy="false">)</mo><mo>=</mo><mi mathvariant="normal">t</mi><mi mathvariant="normal">a</mi><msup><mrow><mi mathvariant="normal">n</mi></mrow><mrow><mo>-</mo><mn mathvariant="normal">1</mn></mrow></msup><mfenced open="{" close="}" separators="|"><mrow><mfrac><mrow><mi mathvariant="normal">I</mi><mi mathvariant="normal">m</mi><mo stretchy="false">[</mo><msub><mrow><mi>H</mi></mrow><mrow><mi>j</mi><mi>m</mi></mrow></msub><mo stretchy="false">(</mo><msub><mrow><mi>ω</mi></mrow><mrow><mi>l</mi></mrow></msub><mo stretchy="false">)</mo><mo stretchy="false">]</mo></mrow><mrow><mi mathvariant="normal">R</mi><mi mathvariant="normal">e</mi><mo stretchy="false">[</mo><msub><mrow><mi>H</mi></mrow><mrow><mi>j</mi><mi>m</mi></mrow></msub><mo stretchy="false">(</mo><msub><mrow><mi>ω</mi></mrow><mrow><mi>l</mi></mrow></msub><mo stretchy="false">)</mo><mo stretchy="false">]</mo></mrow></mfrac></mrow></mfenced></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=24671608&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=24671606&type= 49.95333481 10.58333302 ； <math id="M35"><msub><mrow><mi>ϕ</mi></mrow><mrow><mi>m</mi><mi>l</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=24671612&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=24671610&type= 4.48733330 4.57200003 为相互独立的随机相位角，在 <math id="M36"><mo stretchy="false">[</mo><mn mathvariant="normal">0,2</mn><mi mathvariant="normal">π</mi><mo stretchy="false">]</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=24671699&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=24671698&type= 11.09133339 3.55599999 区间内服从均匀分布。 尽管 式（16） 式（16） 给出了随机过程向量的时间历程显式表达式，但当采样点数 <math id="M37"><mi>N</mi><mo>→</mo><mi mathvariant="normal">∞</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=24671621&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=24671619&type= 11.34533310 2.87866688 时，功率谱密度矩阵接近奇异，无法应用Cholesky方法对其进行分解。为了解决该问题，本文引入波数谱，建立在空间和时间上几乎连续的风速模型，这可以在规定空间范围内得到无穷采样点的风速，从而避免了互功率谱矩阵的Cholesky分解。波数谱的表达式可以由自功率谱和相干函数积分得到，即： <math display="block" id="M38"><mi>S</mi><mo stretchy="false">(</mo><mi>ω</mi><mo>,</mo><mi>κ</mi><mo stretchy="false">)</mo><mo>=</mo><mfrac><mrow><mn mathvariant="normal">1</mn></mrow><mrow><mn mathvariant="normal">2</mn><mi>π</mi></mrow></mfrac><msubsup><mo>∫</mo><mrow><mo>-</mo><mi mathvariant="normal">∞</mi></mrow><mrow><mi mathvariant="normal">∞</mi></mrow></msubsup><mrow><mi>S</mi><mo stretchy="false">(</mo><mi>ω</mi><mo stretchy="false">)</mo><mo>⋅</mo><mi>γ</mi><mo stretchy="false">(</mo><mi>ξ</mi><mo>,</mo><mi>ω</mi><mo stretchy="false">)</mo><msup><mrow><mi mathvariant="normal">e</mi></mrow><mrow><mi>i</mi><mi>κ</mi><mi>ξ</mi></mrow></msup><mi mathvariant="normal">d</mi><mi>ξ</mi></mrow></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=24671625&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=24671623&type= 58.25066376 7.87400007 ， 其中， <math id="M39"><mi>S</mi><mo stretchy="false">(</mo><mi>ω</mi><mo stretchy="false">)</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=24671629&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=24671627&type= 8.12800026 3.55599999 为自谱， <math id="M40"><mi>ξ</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=24671634&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=24671631&type= 1.60866666 3.64066648 为空间距离， <math id="M41"><mi>κ</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=24671638&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=24671636&type= 1.86266661 2.79399991 为波数， <math id="M42"><mi>γ</mi><mo stretchy="false">(</mo><mi>ξ</mi><mo>,</mo><mi>ω</mi><mo stretchy="false">)</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=24671640&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=24671639&type= 11.51466751 3.97933316 为相干函数，通常选择Davenport相干函数模型。 基于波数谱（17），可以得到随机过程向量 X （ t ）的时间历程表达式为： <math display="block" id="M43"><mtable><mtr><mtd><mi>X</mi><mo stretchy="false">(</mo><mi>x</mi><mo>,</mo><mi>t</mi><mo stretchy="false">)</mo><mo>=</mo><mstyle displaystyle="true"><munderover><mo largeop="true">∑</mo><mrow><mi>l</mi><mo>=</mo><mn mathvariant="normal">1</mn></mrow><mrow><msub><mrow><mi>N</mi></mrow><mrow><mi>κ</mi></mrow></msub></mrow></munderover></mstyle><mrow><mstyle displaystyle="true"><munderover><mo largeop="true">∑</mo><mrow><mi>m</mi><mo>=</mo><mn mathvariant="normal">1</mn></mrow><mrow><msub><mrow><mi>N</mi></mrow><mrow><mi>ω</mi></mrow></msub></mrow></munderover></mstyle><mrow><mroot><mrow><mn mathvariant="normal">4</mn><mi>S</mi><mo stretchy="false">(</mo><msub><mrow><mi>κ</mi></mrow><mrow><mi>l</mi></mrow></msub><mo>,</mo><msub><mrow><mi>ω</mi></mrow><mrow><mi>m</mi></mrow></msub><mo stretchy="false">)</mo><mi mathvariant="normal">Δ</mi><mi>κ</mi><mi mathvariant="normal">Δ</mi><mi>ω</mi></mrow><mrow/></mroot></mrow></mrow><mo>⋅</mo></mtd></mtr><mtr><mtd><mo stretchy="false">[</mo><mi mathvariant="normal">c</mi><mi mathvariant="normal">o</mi><mi mathvariant="normal">s</mi><mo stretchy="false">(</mo><msub><mrow><mi>κ</mi></mrow><mrow><mi>l</mi></mrow></msub><mi>x</mi><mo>+</mo><msub><mrow><mi>ω</mi></mrow><mrow><mi>m</mi></mrow></msub><mi>t</mi><mo>+</mo><msubsup><mrow><mi>ϕ</mi></mrow><mrow><mi>m</mi><mi>l</mi></mrow><mrow><mfenced separators="|"><mrow><mn mathvariant="normal">1</mn></mrow></mfenced></mrow></msubsup><mo stretchy="false">)</mo><mo>+</mo></mtd></mtr><mtr><mtd><mi mathvariant="normal">c</mi><mi mathvariant="normal">o</mi><mi mathvariant="normal">s</mi><mo stretchy="false">(</mo><msub><mrow><mi>κ</mi></mrow><mrow><mi>l</mi></mrow></msub><mi>x</mi><mo>-</mo><msub><mrow><mi>ω</mi></mrow><mrow><mi>m</mi></mrow></msub><mi>t</mi><mo>+</mo><msubsup><mrow><mi>ϕ</mi></mrow><mrow><mi>m</mi><mi>l</mi></mrow><mrow><mfenced separators="|"><mrow><mn mathvariant="normal">2</mn></mrow></mfenced></mrow></msubsup><mo stretchy="false">)</mo><mo stretchy="false">]</mo></mtd></mtr></mtable></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=24671642&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=24671641&type= 56.13399887 22.01333237 ， 其中： <math id="M44"><msub><mrow><mi>ω</mi></mrow><mrow><mi>m</mi></mrow></msub><mo>=</mo><mi>m</mi><mi mathvariant="normal">Δ</mi><mi>ω</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=24671691&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=24671690&type= 16.76399994 3.72533321 ， <math id="M45"><mi mathvariant="normal">Δ</mi><mi>ω</mi><mo>=</mo><mfrac><mrow><mn mathvariant="normal">2</mn><mi mathvariant="normal">π</mi></mrow><mrow><mi>T</mi></mrow></mfrac></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=24671646&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=24671645&type= 14.39333248 7.87400007 ， T 为模拟时长； <math id="M46"><mi mathvariant="normal">Δ</mi><mi>κ</mi><mo>=</mo><mfrac><mrow><mn mathvariant="normal">2</mn><mi mathvariant="normal">π</mi></mrow><mrow><mi>L</mi></mrow></mfrac></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=24671648&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=24671647&type= 13.80066681 7.87400007 ， L 为空间尺度。 式（18） 式（18） 仅给出了考虑一维空间频率和时间频率的表达式。然而，对于望远镜来说，尤其是对于大口径主镜来说，风速模拟时考虑平面内二维空间频率更合理。为此，需要将 式（18） 式（18） 拓展至二维空间频率范围。 对于二维空间上的两个不同位置 <math display="block" id="M47"><mo stretchy="false">(</mo><msub><mrow><mi>x</mi></mrow><mrow><mi>i</mi></mrow></msub><mo>,</mo><msub><mrow><mi>y</mi></mrow><mrow><mi>i</mi></mrow></msub><mo stretchy="false">)</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=24671650&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=24671649&type= 10.66800022 4.91066647 和 <math display="block" id="M48"><mo stretchy="false">(</mo><msub><mrow><mi>x</mi></mrow><mrow><mi>j</mi></mrow></msub><mo>,</mo><msub><mrow><mi>y</mi></mrow><mrow><mi>j</mi></mrow></msub><mo stretchy="false">)</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=24671652&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=24671651&type= 10.58333302 4.91066647 ，其Davenport相干函数模型可以表达为： <math display="block" id="M49"><mi>γ</mi><mo stretchy="false">(</mo><msub><mrow><mi>ξ</mi></mrow><mrow><mi>x</mi></mrow></msub><mo>,</mo><msub><mrow><mi>ξ</mi></mrow><mrow><mi>y</mi></mrow></msub><mo>,</mo><mi>ω</mi><mo stretchy="false">)</mo><mo>=</mo><mi mathvariant="normal">e</mi><mi mathvariant="normal">x</mi><mi mathvariant="normal">p</mi><mfenced separators="|"><mrow><mo>-</mo><mfrac><mrow><mi>ω</mi></mrow><mrow><mn mathvariant="normal">2</mn><mi mathvariant="normal">π</mi><msub><mrow><mi>v</mi></mrow><mrow><mn mathvariant="normal">10</mn></mrow></msub></mrow></mfrac><mroot><mrow><msubsup><mrow><mi>λ</mi></mrow><mrow><mi>x</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msubsup><msubsup><mrow><mi>ξ</mi></mrow><mrow><mi>x</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msubsup><mo>+</mo><msubsup><mrow><mi>λ</mi></mrow><mrow><mi>y</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msubsup><msubsup><mrow><mi>ξ</mi></mrow><mrow><mi>y</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msubsup></mrow><mrow/></mroot></mrow></mfenced></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=24671654&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=24671653&type= 65.87065887 9.90600014 ， 其中： <math id="M50"><msub><mrow><mi>ξ</mi></mrow><mrow><mi>x</mi></mrow></msub><mo>=</mo><mfenced open="|" close="|" separators="|"><mrow><msub><mrow><mi>x</mi></mrow><mrow><mi>i</mi></mrow></msub><mo>-</mo><msub><mrow><mi>x</mi></mrow><mrow><mi>j</mi></mrow></msub></mrow></mfenced></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=24671657&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=24671655&type= 21.33600044 4.91066647 和 <math id="M51"><msub><mrow><mi>ξ</mi></mrow><mrow><mi>y</mi></mrow></msub><mo>=</mo><mfenced open="|" close="|" separators="|"><mrow><msub><mrow><mi>y</mi></mrow><mrow><mi>i</mi></mrow></msub><mo>-</mo><msub><mrow><mi>y</mi></mrow><mrow><mi>j</mi></mrow></msub></mrow></mfenced></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=24671659&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=24671658&type= 20.91266632 4.91066647 为空间两个维度上的空间距离； <math id="M52"><msub><mrow><mi>λ</mi></mrow><mrow><mi>x</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=24671664&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=24671661&type= 2.96333337 3.72533321 和 <math id="M53"><msub><mrow><mi>λ</mi></mrow><mrow><mi>y</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=24671666&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=24671665&type= 2.87866688 3.72533321 为指数衰减系数； <math id="M54"><msub><mrow><mi>v</mi></mrow><mrow><mn mathvariant="normal">10</mn></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=24671668&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=24671667&type= 3.97933316 3.72533321 为距离地面高度10 m处的平均风速。 二维空间的波数谱可以由 式（19） 式（19） 积分得到，即： <math display="block" id="M55"><mtable><mtr><mtd><mi>S</mi><mo stretchy="false">(</mo><msub><mrow><mi>κ</mi></mrow><mrow><mi>x</mi></mrow></msub><mo>,</mo><msub><mrow><mi>κ</mi></mrow><mrow><mi>y</mi></mrow></msub><mo>,</mo><mi>ω</mi><mo stretchy="false">)</mo><mo>=</mo><mfrac><mrow><mn mathvariant="normal">1</mn></mrow><mrow><msup><mrow><mo stretchy="false">(</mo><mn mathvariant="normal">2</mn><mi>π</mi><mo stretchy="false">)</mo></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup></mrow></mfrac><msubsup><mo>∫</mo><mrow><mo>-</mo><mi mathvariant="normal">∞</mi></mrow><mrow><mi mathvariant="normal">∞</mi></mrow></msubsup><mrow><msubsup><mo>∫</mo><mrow><mo>-</mo><mi mathvariant="normal">∞</mi></mrow><mrow><mi mathvariant="normal">∞</mi></mrow></msubsup><mrow><msub><mrow><mi>S</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub><mo stretchy="false">(</mo><mi>ω</mi><mo stretchy="false">)</mo><mo>⋅</mo><mi>γ</mi><mo stretchy="false">(</mo><msub><mrow><mi>ξ</mi></mrow><mrow><mi>x</mi></mrow></msub><mo>,</mo><msub><mrow><mi>ξ</mi></mrow><mrow><mi>y</mi></mrow></msub><mo>,</mo><mi>ω</mi><mo stretchy="false">)</mo><msup><mrow><mi mathvariant="normal">e</mi></mrow><mrow><mi>i</mi><mfenced separators="|"><mrow><msub><mrow><mi>κ</mi></mrow><mrow><mi>x</mi></mrow></msub><msub><mrow><mi>ξ</mi></mrow><mrow><mi>x</mi></mrow></msub><mo>+</mo><msub><mrow><mi>κ</mi></mrow><mrow><mi>y</mi></mrow></msub><msub><mrow><mi>ξ</mi></mrow><mrow><mi>y</mi></mrow></msub></mrow></mfenced></mrow></msup><mi mathvariant="normal">d</mi><msub><mrow><mi>ξ</mi></mrow><mrow><mi>x</mi></mrow></msub><mi mathvariant="normal">d</mi><msub><mrow><mi>ξ</mi></mrow><mrow><mi>y</mi></mrow></msub></mrow></mrow><mo>=</mo><mfrac><mrow><mn mathvariant="normal">1</mn></mrow><mrow><mn mathvariant="normal">2</mn><mi mathvariant="normal">π</mi><msub><mrow><mi>λ</mi></mrow><mrow><mi>x</mi></mrow></msub><msub><mrow><mi>λ</mi></mrow><mrow><mi>y</mi></mrow></msub></mrow></mfrac><mfrac><mrow><mn mathvariant="normal">1</mn></mrow><mrow><msup><mrow><mfenced separators="|"><mrow><mfrac><mrow><mn mathvariant="normal">1</mn></mrow><mrow><mn mathvariant="normal">2</mn><mi mathvariant="normal">π</mi><msub><mrow><mi>v</mi></mrow><mrow><mn mathvariant="normal">10</mn></mrow></msub></mrow></mfrac><mi>ω</mi></mrow></mfenced></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup></mrow></mfrac><mo>×</mo></mtd></mtr><mtr><mtd><mfrac><mrow><mn mathvariant="normal">1</mn></mrow><mrow><msup><mrow><mfenced separators="|"><mrow><mn mathvariant="normal">1</mn><mo>+</mo><mrow><mrow><mfenced open="[" close="]" separators="|"><mrow><msup><mrow><mfenced separators="|"><mrow><mfrac><mrow><msub><mrow><mi>κ</mi></mrow><mrow><mi>x</mi></mrow></msub></mrow><mrow><msub><mrow><mi>λ</mi></mrow><mrow><mi>x</mi></mrow></msub></mrow></mfrac></mrow></mfenced></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup><mo>+</mo><msup><mrow><mfenced separators="|"><mrow><mfrac><mrow><msub><mrow><mi>κ</mi></mrow><mrow><mi>y</mi></mrow></msub></mrow><mrow><msub><mrow><mi>λ</mi></mrow><mrow><mi>y</mi></mrow></msub></mrow></mfrac></mrow></mfenced></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup></mrow></mfenced></mrow><mo>/</mo><mrow><msup><mrow><mfenced separators="|"><mrow><mfrac><mrow><mn mathvariant="normal">1</mn></mrow><mrow><mn mathvariant="normal">2</mn><mi mathvariant="normal">π</mi><msub><mrow><mi>v</mi></mrow><mrow><mn mathvariant="normal">10</mn></mrow></msub></mrow></mfrac><mi>ω</mi></mrow></mfenced></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup></mrow></mrow></mrow></mfenced></mrow><mrow><mfrac><mrow><mn mathvariant="normal">3</mn></mrow><mrow><mn mathvariant="normal">2</mn></mrow></mfrac></mrow></msup></mrow></mfrac><mo>⋅</mo><msub><mrow><mi>S</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub><mo stretchy="false">(</mo><mi>ω</mi><mo stretchy="false">)</mo></mtd></mtr></mtable></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=24671670&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=24671669&type= 135.55133057 38.43866730 . 通常描述风扰动的功率谱选择Davenport功率谱，其表达式为： <math display="block" id="M56"><mtable><mtr><mtd><msub><mrow><mi>S</mi></mrow><mrow><mi>D</mi></mrow></msub><mo stretchy="false">(</mo><mi>ω</mi><mo stretchy="false">)</mo><mo>=</mo></mtd></mtr><mtr><mtd><mn mathvariant="normal">4</mn><msubsup><mrow><mi>v</mi></mrow><mrow><mi mathvariant="normal">m</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msubsup><msup><mrow><mfenced separators="|"><mrow><mfrac><mrow><mi>k</mi></mrow><mrow><mi mathvariant="normal">l</mi><mi mathvariant="normal">n</mi><mo stretchy="false">(</mo><mrow><mrow><mi>z</mi></mrow><mo>/</mo><mrow><msub><mrow><mi>z</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub></mrow></mrow><mo stretchy="false">)</mo></mrow></mfrac></mrow></mfenced></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup><mfrac><mrow><msup><mrow><mfenced separators="|"><mrow><mfrac><mrow><mn mathvariant="normal">1200</mn><mi>ω</mi></mrow><mrow><mn mathvariant="normal">2</mn><mi mathvariant="normal">π</mi><msub><mrow><mi>v</mi></mrow><mrow><mn mathvariant="normal">10</mn></mrow></msub></mrow></mfrac></mrow></mfenced></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup></mrow><mrow><mfrac><mrow><mi>ω</mi></mrow><mrow><mn mathvariant="normal">2</mn><mi mathvariant="normal">π</mi></mrow></mfrac><msup><mrow><mfenced open="[" close="]" separators="|"><mrow><mn mathvariant="normal">1</mn><mo>+</mo><msup><mrow><mfenced separators="|"><mrow><mfrac><mrow><mn mathvariant="normal">1200</mn><mi>ω</mi></mrow><mrow><mn mathvariant="normal">2</mn><mi mathvariant="normal">π</mi><msub><mrow><mi>v</mi></mrow><mrow><mn mathvariant="normal">10</mn></mrow></msub></mrow></mfrac></mrow></mfenced></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup></mrow></mfenced></mrow><mrow><mrow><mrow><mn mathvariant="normal">4</mn></mrow><mo>/</mo><mrow><mn mathvariant="normal">3</mn></mrow></mrow></mrow></msup></mrow></mfrac></mtd></mtr></mtable></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=24671672&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=24671671&type= 59.60533142 30.22599983 ， 其中： <math id="M57"><msub><mrow><mi>v</mi></mrow><mrow><mi mathvariant="normal">m</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=24671674&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=24671673&type= 3.55599999 3.72533321 为平均风速， k 为风阻系数。 根据二维空间的波数谱表达 式（20） 式（20） ，可得二维空间随机场的时程表达式为： <math display="block" id="M58"><mtable><mtr><mtd><mi>X</mi><mo stretchy="false">(</mo><mi>x</mi><mo>,</mo><mi>t</mi><mo stretchy="false">)</mo><mo>=</mo></mtd></mtr><mtr><mtd><mstyle displaystyle="true"><munderover><mo largeop="true">∑</mo><mrow><mi>i</mi><mo>=</mo><mn mathvariant="normal">1</mn></mrow><mrow><msub><mrow><mi>N</mi></mrow><mrow><mi>x</mi></mrow></msub></mrow></munderover></mstyle><mrow><mstyle displaystyle="true"><munderover><mo largeop="true">∑</mo><mrow><mi>j</mi><mo>=</mo><mn mathvariant="normal">1</mn></mrow><mrow><msub><mrow><mi>N</mi></mrow><mrow><mi>y</mi></mrow></msub></mrow></munderover></mstyle><mrow><mstyle displaystyle="true"><munderover><mo largeop="true">∑</mo><mrow><mi>m</mi><mo>=</mo><mn mathvariant="normal">1</mn></mrow><mrow><msub><mrow><mi>N</mi></mrow><mrow><mi>ω</mi></mrow></msub></mrow></munderover></mstyle><mrow><mroot><mrow><mn mathvariant="normal">4</mn><mi>S</mi><mo stretchy="false">(</mo><msubsup><mrow><mi>κ</mi></mrow><mrow><mi>x</mi></mrow><mrow><mfenced separators="|"><mrow><mi>i</mi></mrow></mfenced></mrow></msubsup><mo>,</mo><msubsup><mrow><mi>κ</mi></mrow><mrow><mi>y</mi></mrow><mrow><mfenced separators="|"><mrow><mi>j</mi></mrow></mfenced></mrow></msubsup><mo>,</mo><msub><mrow><mi>ω</mi></mrow><mrow><mi>m</mi></mrow></msub><mo stretchy="false">)</mo><mi mathvariant="normal">Δ</mi><msub><mrow><mi>κ</mi></mrow><mrow><mi>x</mi></mrow></msub><mi mathvariant="normal">Δ</mi><msub><mrow><mi>κ</mi></mrow><mrow><mi>y</mi></mrow></msub><mi mathvariant="normal">Δ</mi><mi>ω</mi></mrow><mrow/></mroot></mrow></mrow></mrow><mo>⋅</mo></mtd></mtr><mtr><mtd><mo stretchy="false">[</mo><mi mathvariant="normal">c</mi><mi mathvariant="normal">o</mi><mi mathvariant="normal">s</mi><mo stretchy="false">(</mo><msubsup><mrow><mi>κ</mi></mrow><mrow><mi>x</mi></mrow><mrow><mfenced separators="|"><mrow><mi>i</mi></mrow></mfenced></mrow></msubsup><mi>x</mi><mo>+</mo><msubsup><mrow><mi>κ</mi></mrow><mrow><mi>y</mi></mrow><mrow><mfenced separators="|"><mrow><mi>j</mi></mrow></mfenced></mrow></msubsup><mi>y</mi><mo>+</mo><msub><mrow><mi>ω</mi></mrow><mrow><mi>m</mi></mrow></msub><mi>t</mi><mo>+</mo><msubsup><mrow><mi>ϕ</mi></mrow><mrow><mi>i</mi><mi>j</mi><mi>m</mi></mrow><mrow><mfenced separators="|"><mrow><mn mathvariant="normal">1</mn></mrow></mfenced></mrow></msubsup><mo stretchy="false">)</mo><mo>+</mo></mtd></mtr><mtr><mtd><mi mathvariant="normal">c</mi><mi mathvariant="normal">o</mi><mi mathvariant="normal">s</mi><mo stretchy="false">(</mo><msubsup><mrow><mi>κ</mi></mrow><mrow><mi>x</mi></mrow><mrow><mfenced separators="|"><mrow><mi>i</mi></mrow></mfenced></mrow></msubsup><mi>x</mi><mo>+</mo><msubsup><mrow><mi>κ</mi></mrow><mrow><mi>y</mi></mrow><mrow><mfenced separators="|"><mrow><mi>j</mi></mrow></mfenced></mrow></msubsup><mi>y</mi><mo>-</mo><msub><mrow><mi>ω</mi></mrow><mrow><mi>m</mi></mrow></msub><mi>t</mi><mo>+</mo><msubsup><mrow><mi>ϕ</mi></mrow><mrow><mi>i</mi><mi>j</mi><mi>m</mi></mrow><mrow><mfenced separators="|"><mrow><mn mathvariant="normal">2</mn></mrow></mfenced></mrow></msubsup><mo stretchy="false">)</mo><mo>+</mo></mtd></mtr><mtr><mtd><mi mathvariant="normal">c</mi><mi mathvariant="normal">o</mi><mi mathvariant="normal">s</mi><mo stretchy="false">(</mo><msubsup><mrow><mi>κ</mi></mrow><mrow><mi>x</mi></mrow><mrow><mfenced separators="|"><mrow><mi>i</mi></mrow></mfenced></mrow></msubsup><mi>x</mi><mo>-</mo><msubsup><mrow><mi>κ</mi></mrow><mrow><mi>y</mi></mrow><mrow><mfenced separators="|"><mrow><mi>j</mi></mrow></mfenced></mrow></msubsup><mi>y</mi><mo>+</mo><msub><mrow><mi>ω</mi></mrow><mrow><mi>m</mi></mrow></msub><mi>t</mi><mo>+</mo><msubsup><mrow><mi>ϕ</mi></mrow><mrow><mi>i</mi><mi>j</mi><mi>m</mi></mrow><mrow><mfenced separators="|"><mrow><mn mathvariant="normal">3</mn></mrow></mfenced></mrow></msubsup><mo stretchy="false">)</mo><mo>+</mo></mtd></mtr><mtr><mtd><mi mathvariant="normal">c</mi><mi mathvariant="normal">o</mi><mi mathvariant="normal">s</mi><mo stretchy="false">(</mo><msubsup><mrow><mi>κ</mi></mrow><mrow><mi>x</mi></mrow><mrow><mfenced separators="|"><mrow><mi>i</mi></mrow></mfenced></mrow></msubsup><mi>x</mi><mo>-</mo><msubsup><mrow><mi>κ</mi></mrow><mrow><mi>y</mi></mrow><mrow><mfenced separators="|"><mrow><mi>j</mi></mrow></mfenced></mrow></msubsup><mi>y</mi><mo>-</mo><msub><mrow><mi>ω</mi></mrow><mrow><mi>m</mi></mrow></msub><mi>t</mi><mo>+</mo><msubsup><mrow><mi>ϕ</mi></mrow><mrow><mi>i</mi><mi>j</mi><mi>m</mi></mrow><mrow><mfenced separators="|"><mrow><mn mathvariant="normal">4</mn></mrow></mfenced></mrow></msubsup><mo stretchy="false">)</mo><mo stretchy="false">]</mo></mtd></mtr></mtable></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=24671676&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=24671675&type= 61.63733673 41.40200043 ， 其中： <math id="M59"><msubsup><mrow><mi>κ</mi></mrow><mrow><mi>x</mi></mrow><mrow><mfenced separators="|"><mrow><mi>i</mi></mrow></mfenced></mrow></msubsup><mo>=</mo><mi>i</mi><mi mathvariant="normal">Δ</mi><msub><mrow><mi>κ</mi></mrow><mrow><mi>x</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=24671679&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=24671678&type= 16.08666801 4.57200003 ， <math id="M60"><mi>i</mi><mo>=</mo><mn mathvariant="normal">1,2</mn><mo>,</mo><mo>…</mo><mo>,</mo><msub><mrow><mi>N</mi></mrow><mrow><mi>x</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=24671681&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=24671680&type= 23.19866753 4.14866638 ， <math id="M61"><mi mathvariant="normal">Δ</mi><msub><mrow><mi>κ</mi></mrow><mrow><mi>x</mi></mrow></msub><mo>=</mo><mfrac><mrow><mn mathvariant="normal">2</mn><mi mathvariant="normal">π</mi></mrow><mrow><msub><mrow><mi>L</mi></mrow><mrow><mi>x</mi></mrow></msub></mrow></mfrac></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=24671683&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=24671682&type= 14.98600006 8.80533314 ； <math id="M62"><msubsup><mrow><mi>κ</mi></mrow><mrow><mi>y</mi></mrow><mrow><mfenced separators="|"><mrow><mi>j</mi></mrow></mfenced></mrow></msubsup><mo>=</mo><mi>j</mi><mi mathvariant="normal">Δ</mi><msub><mrow><mi>κ</mi></mrow><mrow><mi>y</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=24671685&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=24671684&type= 16.17133331 5.41866684 ， <math id="M63"><mi>j</mi><mo>=</mo><mn mathvariant="normal">1,2</mn><mo>,</mo><mo>…</mo><mo>,</mo><msub><mrow><mi>N</mi></mrow><mrow><mi>y</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=24671687&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=24671686&type= 23.02933502 4.57200003 ， <math id="M64"><mi mathvariant="normal">Δ</mi><msub><mrow><mi>κ</mi></mrow><mrow><mi>y</mi></mrow></msub><mo>=</mo><mfrac><mrow><mn mathvariant="normal">2</mn><mi mathvariant="normal">π</mi></mrow><mrow><msub><mrow><mi>L</mi></mrow><mrow><mi>y</mi></mrow></msub></mrow></mfrac></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=24671689&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=24671688&type= 14.90133286 8.80533314 ； <math id="M65"><msub><mrow><mi>ω</mi></mrow><mrow><mi>m</mi></mrow></msub><mo>=</mo><mi>m</mi><mi mathvariant="normal">Δ</mi><mi>ω</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=24671691&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=24671690&type= 16.76399994 3.72533321 ， <math id="M66"><mi>m</mi><mo>=</mo><mn mathvariant="normal">1,2</mn><mo>,</mo><mo>…</mo><mo>,</mo><msub><mrow><mi>N</mi></mrow><mrow><mi>ω</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=24671693&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=24671692&type= 25.14599991 4.14866638 ， <math id="M67"><mi mathvariant="normal">Δ</mi><mi>ω</mi><mo>=</mo><mfrac><mrow><mn mathvariant="normal">2</mn><mi mathvariant="normal">π</mi></mrow><mrow><mi>T</mi></mrow></mfrac></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=24671695&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=24671694&type= 14.39333248 7.87400007 ； <math id="M68"><msubsup><mrow><mi>ϕ</mi></mrow><mrow><mi>i</mi><mi>j</mi><mi>m</mi></mrow><mrow><mfenced separators="|"><mrow><mn mathvariant="normal">1,2</mn><mo>,</mo><mn mathvariant="normal">3,4</mn></mrow></mfenced></mrow></msubsup></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=24671697&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=24671696&type= 10.66800022 4.06400013 为4组相互独立的相位角，同样在 <math id="M69"><mo stretchy="false">[</mo><mn mathvariant="normal">0,2</mn><mi mathvariant="normal">π</mi><mo stretchy="false">]</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=24671699&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=24671698&type= 11.09133339 3.55599999 区间内服从均匀分布。 3.2　FFT算法 2 式（18） 式（18） 和 式（22） 式（22） 分别给出了一维空间和二维空间情况下的风速时程表达式，具有相同的表达形式，分别为二重和三重求和表达式。随着采样点数量的增大，计算效率将会显著下降。为了改善计算效率，这里引入快速傅里叶算法FFT。 为了利用FFT，将 式（18） 式（18） 改写为： <math display="block" id="M70"><mtable><mtr><mtd><mi>X</mi><mo stretchy="false">(</mo><mi>x</mi><mo>,</mo><mi>t</mi><mo stretchy="false">)</mo><mo>=</mo></mtd></mtr><mtr><mtd><mi mathvariant="normal">R</mi><mi mathvariant="normal">e</mi><mfenced open="{" close="}" separators="|"><mrow><mstyle displaystyle="true"><munderover><mo largeop="true">∑</mo><mrow><mi>l</mi><mo>=</mo><mn mathvariant="normal">0</mn></mrow><mrow><msub><mrow><mi>N</mi></mrow><mrow><mi>κ</mi></mrow></msub><mo>-</mo><mn mathvariant="normal">1</mn></mrow></munderover></mstyle><mrow><mstyle displaystyle="true"><munderover><mo largeop="true">∑</mo><mrow><mi>m</mi><mo>=</mo><mn mathvariant="normal">0</mn></mrow><mrow><msub><mrow><mi>N</mi></mrow><mrow><mi>ω</mi></mrow></msub><mo>-</mo><mn mathvariant="normal">1</mn></mrow></munderover></mstyle><mrow><mo stretchy="false">[</mo><msubsup><mrow><mi>B</mi></mrow><mrow><mi>l</mi><mi>m</mi></mrow><mrow><mfenced separators="|"><mrow><mn mathvariant="normal">1</mn></mrow></mfenced></mrow></msubsup><msup><mrow><mi mathvariant="normal">e</mi></mrow><mrow><mi>i</mi><mfenced separators="|"><mrow><msub><mrow><mi>κ</mi></mrow><mrow><mi>l</mi></mrow></msub><mi>x</mi><mo>+</mo><msub><mrow><mi>ω</mi></mrow><mrow><mi>m</mi></mrow></msub><mi>t</mi></mrow></mfenced></mrow></msup><mo>+</mo><msubsup><mrow><mi>B</mi></mrow><mrow><mi>l</mi><mi>m</mi></mrow><mrow><mfenced separators="|"><mrow><mn mathvariant="normal">2</mn></mrow></mfenced></mrow></msubsup><msup><mrow><mi mathvariant="normal">e</mi></mrow><mrow><mi>i</mi><mfenced separators="|"><mrow><msub><mrow><mi>κ</mi></mrow><mrow><mi>l</mi></mrow></msub><mi>x</mi><mo>-</mo><msub><mrow><mi>ω</mi></mrow><mrow><mi>m</mi></mrow></msub><mi>t</mi></mrow></mfenced></mrow></msup><mo stretchy="false">]</mo></mrow></mrow></mrow></mfenced></mtd></mtr></mtable></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=24671701&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=24671700&type= 62.99200058 16.00200081 ， 其中， <math id="M71"><msubsup><mrow><mi>B</mi></mrow><mrow><mi>l</mi><mi>m</mi></mrow><mrow><mfenced separators="|"><mrow><mi>n</mi></mrow></mfenced></mrow></msubsup><mo>=</mo><mroot><mrow><mn mathvariant="normal">4</mn><mi>S</mi><mo stretchy="false">(</mo><msub><mrow><mi>κ</mi></mrow><mrow><mi>l</mi></mrow></msub><mo>,</mo><msub><mrow><mi>ω</mi></mrow><mrow><mi>m</mi></mrow></msub><mo stretchy="false">)</mo><mi mathvariant="normal">Δ</mi><mi>κ</mi><mi mathvariant="normal">Δ</mi><mi>ω</mi></mrow><mrow/></mroot><mi mathvariant="normal">e</mi><mi mathvariant="normal">x</mi><mi mathvariant="normal">p</mi><mo stretchy="false">(</mo><mi>i</mi><msubsup><mrow><mi>ϕ</mi></mrow><mrow><mi>l</mi><mi>m</mi></mrow><mrow><mfenced separators="|"><mrow><mi>n</mi></mrow></mfenced></mrow></msubsup><mo stretchy="false">)</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=24671703&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=24671702&type= 55.45666885 5.33400011 。 根据FFT理论， 式（23） 式（23） 可以进一步改写为： <math display="block" id="M72"><mtable><mtr><mtd><mi>X</mi><mo stretchy="false">(</mo><mi>x</mi><mo>,</mo><mi>t</mi><mo stretchy="false">)</mo><mo>=</mo><mi mathvariant="normal">R</mi><mi mathvariant="normal">e</mi><mo stretchy="false">{</mo><mi>I</mi><mi>F</mi><mi>F</mi><msub><mrow><mi>T</mi></mrow><mrow><mi>κ</mi></mrow></msub><mo stretchy="false">[</mo><mi>I</mi><mi>F</mi><mi>F</mi><msub><mrow><mi>T</mi></mrow><mrow><mi>ω</mi></mrow></msub><mo stretchy="false">(</mo><msubsup><mrow><mi>B</mi></mrow><mrow><mi>l</mi><mi>m</mi></mrow><mrow><mfenced separators="|"><mrow><mn mathvariant="normal">1</mn></mrow></mfenced></mrow></msubsup><mo stretchy="false">)</mo><mo stretchy="false">]</mo><mo>+</mo></mtd></mtr><mtr><mtd><mi>I</mi><mi>F</mi><mi>F</mi><msub><mrow><mi>T</mi></mrow><mrow><mi>κ</mi></mrow></msub><mo stretchy="false">[</mo><mi>F</mi><mi>F</mi><msub><mrow><mi>T</mi></mrow><mrow><mi>ω</mi></mrow></msub><mo stretchy="false">(</mo><msubsup><mrow><mi>B</mi></mrow><mrow><mi>l</mi><mi>m</mi></mrow><mrow><mfenced separators="|"><mrow><mn mathvariant="normal">2</mn></mrow></mfenced></mrow></msubsup><mo stretchy="false">)</mo><mo stretchy="false">]</mo><mo stretchy="false">}</mo></mtd></mtr></mtable></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=24671705&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=24671704&type= 61.12933731 10.07533360 ， 其中， <math id="M73"><mi>F</mi><mi>F</mi><mi>T</mi><mo stretchy="false">(</mo><mo>⋅</mo><mo stretchy="false">)</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=24671707&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=24671706&type= 11.00666618 3.55599999 和 <math id="M74"><mi>I</mi><mi>F</mi><mi>F</mi><mi>T</mi><mo stretchy="false">(</mo><mo>⋅</mo><mo stretchy="false">)</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=24671709&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=24671708&type= 12.19200039 3.55599999 分别表示快速傅里叶变换和逆变换。 同样， 式（22） 式（22） 可以变换为： <math display="block" id="M75"><mtable columnalign="left"><mtr><mtd><mi>X</mi><mo stretchy="false">(</mo><msub><mrow><mi>p</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow></msub><mi mathvariant="normal">Δ</mi><mi>x</mi><mo>,</mo><msub><mrow><mi>p</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msub><mi mathvariant="normal">Δ</mi><mi>y</mi><mo>,</mo><msub><mrow><mi>p</mi></mrow><mrow><mn mathvariant="normal">3</mn></mrow></msub><mi mathvariant="normal">Δ</mi><mi>t</mi><mo stretchy="false">)</mo><mo>=</mo><mi mathvariant="normal">R</mi><mi mathvariant="normal">e</mi><mo stretchy="false">{</mo><mstyle displaystyle="true"><munderover><mo largeop="true">∑</mo><mrow><msub><mrow><mi>n</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow></msub><mo>=</mo><mn mathvariant="normal">0</mn></mrow><mrow><msub><mrow><mi>M</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow></msub><mo>-</mo><mn mathvariant="normal">1</mn></mrow></munderover></mstyle><mrow><mstyle displaystyle="true"><munderover><mo largeop="true">∑</mo><mrow><msub><mrow><mi>n</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msub><mo>=</mo><mn mathvariant="normal">0</mn></mrow><mrow><msub><mrow><mi>M</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msub><mo>-</mo><mn mathvariant="normal">1</mn></mrow></munderover></mstyle><mrow><mstyle displaystyle="true"><munderover><mo largeop="true">∑</mo><mrow><msub><mrow><mi>n</mi></mrow><mrow><mn mathvariant="normal">3</mn></mrow></msub><mo>=</mo><mn mathvariant="normal">0</mn></mrow><mrow><msub><mrow><mi>M</mi></mrow><mrow><mn mathvariant="normal">3</mn></mrow></msub><mo>-</mo><mn mathvariant="normal">1</mn></mrow></munderover></mstyle><mrow><mroot><mrow><mn mathvariant="normal">4</mn><mi>S</mi><mfenced separators="|"><mrow><msubsup><mrow><mi>κ</mi></mrow><mrow><mi>x</mi></mrow><mrow><mfenced separators="|"><mrow><msub><mrow><mi>n</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow></msub></mrow></mfenced></mrow></msubsup><mo>,</mo><msubsup><mrow><mi>κ</mi></mrow><mrow><mi>y</mi></mrow><mrow><mfenced separators="|"><mrow><msub><mrow><mi>n</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msub></mrow></mfenced></mrow></msubsup><mo>,</mo><msub><mrow><mi>ω</mi></mrow><mrow><msub><mrow><mi>n</mi></mrow><mrow><mn mathvariant="normal">3</mn></mrow></msub></mrow></msub></mrow></mfenced><mi mathvariant="normal">Δ</mi><msub><mrow><mi>κ</mi></mrow><mrow><mi>x</mi></mrow></msub><mi mathvariant="normal">Δ</mi><msub><mrow><mi>κ</mi></mrow><mrow><mi>y</mi></mrow></msub><mi mathvariant="normal">Δ</mi><mi>ω</mi></mrow><mrow/></mroot></mrow></mrow></mrow><mo>⋅</mo></mtd></mtr><mtr><mtd><mo stretchy="false">[</mo><msubsup><mrow><mi>B</mi></mrow><mrow><msub><mrow><mi>n</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow></msub><msub><mrow><mi>n</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msub><msub><mrow><mi>n</mi></mrow><mrow><mn mathvariant="normal">3</mn></mrow></msub></mrow><mrow><mfenced separators="|"><mrow><mn mathvariant="normal">1</mn></mrow></mfenced></mrow></msubsup><mi mathvariant="normal">e</mi><mi mathvariant="normal">x</mi><mi mathvariant="normal">p</mi><mfenced separators="|"><mrow><mi>i</mi><mfrac><mrow><mn mathvariant="normal">2</mn><mi mathvariant="normal">π</mi><msub><mrow><mi>n</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow></msub><msub><mrow><mi>p</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow></msub></mrow><mrow><msub><mrow><mi>M</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow></msub></mrow></mfrac><mo>+</mo><mi>i</mi><mfrac><mrow><mn mathvariant="normal">2</mn><mi mathvariant="normal">π</mi><msub><mrow><mi>n</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msub><msub><mrow><mi>p</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msub></mrow><mrow><msub><mrow><mi>M</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msub></mrow></mfrac><mo>+</mo><mi>i</mi><mfrac><mrow><mn mathvariant="normal">2</mn><mi mathvariant="normal">π</mi><msub><mrow><mi>n</mi></mrow><mrow><mn mathvariant="normal">3</mn></mrow></msub><msub><mrow><mi>p</mi></mrow><mrow><mn mathvariant="normal">3</mn></mrow></msub></mrow><mrow><msub><mrow><mi>M</mi></mrow><mrow><mn mathvariant="normal">3</mn></mrow></msub></mrow></mfrac></mrow></mfenced><mo>+</mo><msubsup><mrow><mi>B</mi></mrow><mrow><msub><mrow><mi>n</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow></msub><msub><mrow><mi>n</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msub><msub><mrow><mi>n</mi></mrow><mrow><mn mathvariant="normal">3</mn></mrow></msub></mrow><mrow><mfenced separators="|"><mrow><mn mathvariant="normal">2</mn></mrow></mfenced></mrow></msubsup><mi mathvariant="normal">e</mi><mi mathvariant="normal">x</mi><mi mathvariant="normal">p</mi><mfenced separators="|"><mrow><mi>i</mi><mfrac><mrow><mn mathvariant="normal">2</mn><mi mathvariant="normal">π</mi><msub><mrow><mi>n</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow></msub><msub><mrow><mi>p</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow></msub></mrow><mrow><msub><mrow><mi>M</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow></msub></mrow></mfrac><mo>+</mo><mi>i</mi><mfrac><mrow><mn mathvariant="normal">2</mn><mi mathvariant="normal">π</mi><msub><mrow><mi>n</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msub><msub><mrow><mi>p</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msub></mrow><mrow><msub><mrow><mi>M</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msub></mrow></mfrac><mo>-</mo><mi>i</mi><mfrac><mrow><mn mathvariant="normal">2</mn><mi mathvariant="normal">π</mi><msub><mrow><mi>n</mi></mrow><mrow><mn mathvariant="normal">3</mn></mrow></msub><msub><mrow><mi>p</mi></mrow><mrow><mn mathvariant="normal">3</mn></mrow></msub></mrow><mrow><msub><mrow><mi>M</mi></mrow><mrow><mn mathvariant="normal">3</mn></mrow></msub></mrow></mfrac></mrow></mfenced><mo>+</mo></mtd></mtr><mtr><mtd><msubsup><mrow><mi>B</mi></mrow><mrow><msub><mrow><mi>n</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow></msub><msub><mrow><mi>n</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msub><msub><mrow><mi>n</mi></mrow><mrow><mn mathvariant="normal">3</mn></mrow></msub></mrow><mrow><mfenced separators="|"><mrow><mn mathvariant="normal">3</mn></mrow></mfenced></mrow></msubsup><mi mathvariant="normal">e</mi><mi mathvariant="normal">x</mi><mi mathvariant="normal">p</mi><mfenced separators="|"><mrow><mi>i</mi><mfrac><mrow><mn mathvariant="normal">2</mn><mi mathvariant="normal">π</mi><msub><mrow><mi>n</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow></msub><msub><mrow><mi>p</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow></msub></mrow><mrow><msub><mrow><mi>M</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow></msub></mrow></mfrac><mo>-</mo><mi>i</mi><mfrac><mrow><mn mathvariant="normal">2</mn><mi mathvariant="normal">π</mi><msub><mrow><mi>n</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msub><msub><mrow><mi>p</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msub></mrow><mrow><msub><mrow><mi>M</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msub></mrow></mfrac><mo>+</mo><mi>i</mi><mfrac><mrow><mn mathvariant="normal">2</mn><mi mathvariant="normal">π</mi><msub><mrow><mi>n</mi></mrow><mrow><mn mathvariant="normal">3</mn></mrow></msub><msub><mrow><mi>p</mi></mrow><mrow><mn mathvariant="normal">3</mn></mrow></msub></mrow><mrow><msub><mrow><mi>M</mi></mrow><mrow><mn mathvariant="normal">3</mn></mrow></msub></mrow></mfrac></mrow></mfenced><mo>+</mo><msubsup><mrow><mi>B</mi></mrow><mrow><msub><mrow><mi>n</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow></msub><msub><mrow><mi>n</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msub><msub><mrow><mi>n</mi></mrow><mrow><mn mathvariant="normal">3</mn></mrow></msub></mrow><mrow><mfenced separators="|"><mrow><mn mathvariant="normal">4</mn></mrow></mfenced></mrow></msubsup><mi mathvariant="normal">e</mi><mi mathvariant="normal">x</mi><mi mathvariant="normal">p</mi><mfenced separators="|"><mrow><mi>i</mi><mfrac><mrow><mn mathvariant="normal">2</mn><mi mathvariant="normal">π</mi><msub><mrow><mi>n</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow></msub><msub><mrow><mi>p</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow></msub></mrow><mrow><msub><mrow><mi>M</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow></msub></mrow></mfrac><mo>-</mo><mi>i</mi><mfrac><mrow><mn mathvariant="normal">2</mn><mi mathvariant="normal">π</mi><msub><mrow><mi>n</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msub><msub><mrow><mi>p</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msub></mrow><mrow><msub><mrow><mi>M</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msub></mrow></mfrac><mo>-</mo><mi>i</mi><mfrac><mrow><mn mathvariant="normal">2</mn><mi mathvariant="normal">π</mi><msub><mrow><mi>n</mi></mrow><mrow><mn mathvariant="normal">3</mn></mrow></msub><msub><mrow><mi>p</mi></mrow><mrow><mn mathvariant="normal">3</mn></mrow></msub></mrow><mrow><msub><mrow><mi>M</mi></mrow><mrow><mn mathvariant="normal">3</mn></mrow></msub></mrow></mfrac></mrow></mfenced><mo stretchy="false">]</mo><mo stretchy="false">}</mo></mtd></mtr></mtable></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=24671711&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=24671710&type= 146.98133850 33.69733429 ， 其中， <math id="M76"><msubsup><mrow><mi>B</mi></mrow><mrow><msub><mrow><mi>n</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow></msub><msub><mrow><mi>n</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msub><msub><mrow><mi>n</mi></mrow><mrow><mn mathvariant="normal">3</mn></mrow></msub></mrow><mrow><mfenced separators="|"><mrow><mi>k</mi></mrow></mfenced></mrow></msubsup><mo>=</mo><mi mathvariant="normal">e</mi><mi mathvariant="normal">x</mi><mi mathvariant="normal">p</mi><mo stretchy="false">(</mo><mi>i</mi><msubsup><mrow><mi>ϕ</mi></mrow><mrow><msub><mrow><mi>n</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow></msub><msub><mrow><mi>n</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msub><msub><mrow><mi>n</mi></mrow><mrow><mn mathvariant="normal">3</mn></mrow></msub></mrow><mrow><mfenced separators="|"><mrow><mi>k</mi></mrow></mfenced></mrow></msubsup><mo stretchy="false">)</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=24671713&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=24671712&type= 31.58066559 4.48733330 ， <math id="M77"><mi>k</mi><mo>=</mo><mn mathvariant="normal">1,2</mn><mo>,</mo><mn mathvariant="normal">3,4</mn></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=24671715&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=24671714&type= 19.47333336 3.21733332 。 同样由FFT理论， 式（25） 式（25） 可以进一步改写为： <math display="block" id="M78"><mtable><mtr><mtd><mi>X</mi><mo stretchy="false">(</mo><msub><mrow><mi>p</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow></msub><mi mathvariant="normal">Δ</mi><mi>x</mi><mo>,</mo><msub><mrow><mi>p</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msub><mi mathvariant="normal">Δ</mi><mi>y</mi><mo>,</mo><msub><mrow><mi>p</mi></mrow><mrow><mn mathvariant="normal">3</mn></mrow></msub><mi mathvariant="normal">Δ</mi><mi>t</mi><mo stretchy="false">)</mo><mo>=</mo><mi mathvariant="normal">R</mi><mi mathvariant="normal">e</mi><mfenced open="{" close="" separators="|"><mrow><mi>I</mi><mi>F</mi><mi>F</mi><msub><mrow><mi>T</mi></mrow><mrow><mi>x</mi></mrow></msub><mfenced open="{" close="}" separators="|"><mrow><mi>I</mi><mi>F</mi><mi>F</mi><msub><mrow><mi>T</mi></mrow><mrow><mi>y</mi></mrow></msub><mfenced open="[" close="]" separators="|"><mrow><mi>I</mi><mi>F</mi><mi>F</mi><msub><mrow><mi>T</mi></mrow><mrow><mi>ω</mi></mrow></msub><mfenced separators="|"><mrow><msubsup><mrow><mi>B</mi></mrow><mrow><msub><mrow><mi>n</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow></msub><msub><mrow><mi>n</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msub><msub><mrow><mi>n</mi></mrow><mrow><mn mathvariant="normal">3</mn></mrow></msub></mrow><mrow><mfenced separators="|"><mrow><mn mathvariant="normal">1</mn></mrow></mfenced></mrow></msubsup></mrow></mfenced></mrow></mfenced></mrow></mfenced><mo>+</mo><mi>I</mi><mi>F</mi><mi>F</mi><msub><mrow><mi>T</mi></mrow><mrow><mi>x</mi></mrow></msub><mfenced open="{" close="}" separators="|"><mrow><mi>I</mi><mi>F</mi><mi>F</mi><msub><mrow><mi>T</mi></mrow><mrow><mi>y</mi></mrow></msub><mfenced open="[" close="]" separators="|"><mrow><mi>F</mi><mi>F</mi><msub><mrow><mi>T</mi></mrow><mrow><mi>ω</mi></mrow></msub><mfenced separators="|"><mrow><msubsup><mrow><mi>B</mi></mrow><mrow><msub><mrow><mi>n</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow></msub><msub><mrow><mi>n</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msub><msub><mrow><mi>n</mi></mrow><mrow><mn mathvariant="normal">3</mn></mrow></msub></mrow><mrow><mfenced separators="|"><mrow><mn mathvariant="normal">2</mn></mrow></mfenced></mrow></msubsup></mrow></mfenced></mrow></mfenced></mrow></mfenced><mo>+</mo></mrow></mfenced></mtd></mtr><mtr><mtd><mfenced open="" close="}" separators="|"><mrow><mi>I</mi><mi>F</mi><mi>F</mi><msub><mrow><mi>T</mi></mrow><mrow><mi>x</mi></mrow></msub><mfenced open="{" close="}" separators="|"><mrow><mi>F</mi><mi>F</mi><msub><mrow><mi>T</mi></mrow><mrow><mi>y</mi></mrow></msub><mfenced open="[" close="]" separators="|"><mrow><mi>I</mi><mi>F</mi><mi>F</mi><msub><mrow><mi>T</mi></mrow><mrow><mi>ω</mi></mrow></msub><mfenced separators="|"><mrow><msubsup><mrow><mi>B</mi></mrow><mrow><msub><mrow><mi>n</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow></msub><msub><mrow><mi>n</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msub><msub><mrow><mi>n</mi></mrow><mrow><mn mathvariant="normal">3</mn></mrow></msub></mrow><mrow><mfenced separators="|"><mrow><mn mathvariant="normal">3</mn></mrow></mfenced></mrow></msubsup></mrow></mfenced></mrow></mfenced></mrow></mfenced><mo>+</mo><mi>I</mi><mi>F</mi><mi>F</mi><msub><mrow><mi>T</mi></mrow><mrow><mi>x</mi></mrow></msub><mfenced open="{" close="}" separators="|"><mrow><mi>F</mi><mi>F</mi><msub><mrow><mi>T</mi></mrow><mrow><mi>y</mi></mrow></msub><mfenced open="[" close="]" separators="|"><mrow><mi>F</mi><mi>F</mi><msub><mrow><mi>T</mi></mrow><mrow><mi>ω</mi></mrow></msub><mfenced separators="|"><mrow><msubsup><mrow><mi>B</mi></mrow><mrow><msub><mrow><mi>n</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow></msub><msub><mrow><mi>n</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msub><msub><mrow><mi>n</mi></mrow><mrow><mn mathvariant="normal">3</mn></mrow></msub></mrow><mrow><mfenced separators="|"><mrow><mn mathvariant="normal">4</mn></mrow></mfenced></mrow></msubsup></mrow></mfenced></mrow></mfenced></mrow></mfenced></mrow></mfenced></mtd></mtr></mtable></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=24671717&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=24671716&type= 149.09800720 16.67933464 . 对比 式（26） 式（26） 和（22）可以发现，通过引入FFT算法，空间频率和时间频率的多重求和被FFT算法代替，从而显著改善了风速模拟的计算效率。 3.3　风载荷计算方法 2 从前面推导已知风速具有随机变化的特征，风速的幅值 v 可以表达为常量或平均风速 v m 和大气湍流引起的时变量 <math id="M79"><mi mathvariant="normal">Δ</mi><mi>v</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=24671735&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=24671734&type= 4.31799984 2.79399991 之和，即： <math display="block" id="M80"><mi>v</mi><mo>=</mo><msub><mrow><mi>v</mi></mrow><mrow><mi mathvariant="normal">m</mi></mrow></msub><mo>+</mo><mi mathvariant="normal">Δ</mi><mi>v</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=24671721&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=24671720&type= 19.38866615 3.80999994 . 风载荷的表达式为： <math display="block" id="M81"><mi>F</mi><mo>=</mo><mfrac><mrow><mn mathvariant="normal">1</mn></mrow><mrow><mn mathvariant="normal">2</mn></mrow></mfrac><mi>ρ</mi><msub><mrow><mi>C</mi></mrow><mrow><mi mathvariant="normal">D</mi></mrow></msub><mi>A</mi><msup><mrow><mi>v</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=24671723&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=24671722&type= 23.28333473 7.87400007 ， 其中： <math id="M82"><mi>ρ</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=24671725&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=24671724&type= 1.86266661 3.64066648 为空气密度， <math id="M83"><msub><mrow><mi>C</mi></mrow><mrow><mi mathvariant="normal">D</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=24671727&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=24671726&type= 4.23333359 3.72533321 为风阻系数， <math id="M84"><mi>A</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=24671729&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=24671728&type= 2.87866688 2.79399991 为迎风面积。 同样，风载荷也可以分解为稳态量和非稳态量的和，即： <math display="block" id="M85"><mi>F</mi><mo>=</mo><msub><mrow><mi>F</mi></mrow><mrow><mi mathvariant="normal">m</mi></mrow></msub><mo>+</mo><msub><mrow><mi>F</mi></mrow><mrow><mi mathvariant="normal">w</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=24671731&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=24671730&type= 20.57399940 3.80999994 . 非稳态风载荷 <math id="M86"><msub><mrow><mi>F</mi></mrow><mrow><mi mathvariant="normal">w</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=24671733&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=24671732&type= 4.14866638 3.72533321 与 <math id="M87"><mi mathvariant="normal">Δ</mi><mi>v</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=24671735&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=24671734&type= 4.31799984 2.79399991 有关，通过对（28）Taylor展开可得： <math display="block" id="M88"><msub><mrow><mi>F</mi></mrow><mrow><mi mathvariant="normal">w</mi></mrow></msub><mo>=</mo><msub><mrow><mfenced open="" close="|" separators="|"><mrow><mfrac><mrow><mo>∂</mo><mi>F</mi></mrow><mrow><mo>∂</mo><mi>v</mi></mrow></mfrac></mrow></mfenced></mrow><mrow><mi>v</mi><mo>=</mo><msub><mrow><mi>v</mi></mrow><mrow><mi mathvariant="normal">m</mi></mrow></msub></mrow></msub><mi mathvariant="normal">Δ</mi><mi>v</mi><mo>=</mo><mi>ρ</mi><msub><mrow><mi>C</mi></mrow><mrow><mi mathvariant="normal">D</mi></mrow></msub><mi>A</mi><msub><mrow><mi>v</mi></mrow><mrow><mi mathvariant="normal">m</mi></mrow></msub><mi mathvariant="normal">Δ</mi><mi>v</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=24671737&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=24671736&type= 48.17533112 10.58333302 . 望远镜结构非常复杂，作用在望远镜上的风载荷也同样非常复杂。从望远镜结构的角度，主要考虑 图3 图3 所示的三方面作用，即主镜、桁架以及环梁和次镜上的风扰动作用。为了简化分析计算难度，将望远镜划分为三层，如 图4 图4 所示，以这三层的位置来离散采样点，然后根据 式（30） 式（30） 来计算风载荷，并根据结构状态方程（10）和（11）来计算风扰动响应。 4 数值算例 1 以 图1 图1 所示的2 m口径地基望远镜为例，采用本文提出的基于二维随机场的风扰动建模方法模拟望远镜的风扰动环境，并与已有文献进行对比，其中与望远镜结构相关的参数如 表1 表1 所示。此外，对随机风扰动对望远镜的影响进行了详细的仿真分析并对光学性能进行预测。 Parameters Meaning Value v m External wind speed 10 m/s and 15 m/s ρ Air density 1.293 kg/m 3 A sm Area for SM and baffle 227 500 mm 2 A truss Area for truss 1 025 088 mm 2 A spider Area for spider 809 600 mm 2 A ring Area for top ring 392 370 mm 2 C D，sm Drag coefficient for SM 1.2 C D，truss Drag coefficient for truss 1.2 C D，spider Drag coefficient for spiders 1.2 C D，pm Drag coefficient for PM 1.05 L x Spatial length in x direction 2 500 mm L y Spatial length in y direction 2 500 mm H Height 6 400 mm 4.1　随机风扰动的模拟与验证 2 根据3.2节所述二维随机场的理论，按 图4 图4 所示将作用在望远镜的随机风速场简化为三层离散网格并在这些网格点生成随机风速数据，随机风速场离散参数如 表2 表2 所示。以主镜为例，在主镜镜面位置按 图5 图5 所示随机采样6个位置，外界平均风速10 m/s情况下，风速的时程变化过程以及相应的功率谱如 图6 图6 所示。 Parameters Meaning Value z Height above ground 10 m z 0 Ground roughness coefficient 0.01 λ x Decay coefficient in x direction 7 λ y Decay coefficient in y dir 7 N x Wave-number in x dir 128 N y Wave-number in y dir 128 N ω Frequencies 1 024 M 1 Sample numbers in x dir 256 M 2 Sample numbers in y dir 256 M 3 Sample numbers in time domain 2048 T Simulation time length 100s ∆κ x Wave-number space in x dir 0.002 5 mm -1 ∆κ y Wave-number space in y dir 0.002 5 mm -1 ∆ ω Frequency space 0.062 8 s -1 ∆ x Spatial step in x dir 9.766 mm ∆ y Spatial step in y dir 9.766 mm ∆ t Time step 0.048 8 s 1;2;3;4;5;6; 从 图6 图6 所示的主镜上6个采样点位置的风速时程变化曲线以及功率谱曲线可以得出以下结论： （1）从时间历程上看，各个采样点的风速均随时间随机变化； （2）在模拟时长范围内，各个采样点风速时程的功率谱与初始给定的Davenport功率变化趋势基本一致，且与文献［ 18 18 ］给出的结论相符。 从以上结论可以看出，基于二维随机场理论而生成的风速场具有空间和时间随机性的特征，这与实际望远镜使用环境中风速的湍流特性是基本一致的，以此来模拟实际中的风速环境是合理可行的。 4.2　风扰动分析及光学性能预测 2 利用4.1节给出的随机风速场，按3.3节所述的计算方法得到平均风速10 m/s和15 m/s情况下风载荷时间历程数据，并加载到望远镜结构动力学模型中，即可模拟出望远镜上各个位置的风扰动响应。由于望远镜工作中主要关注光学系统性能，而光学系统主要受主、次镜的镜面面形以及相对位置影响，因此在性能模拟与预测时主要考虑主、次镜镜面节点位置的动态变化情况。根据离散节点位置变化，采用非线性最佳拟合的方法计算主、次镜镜面面形精度以及相对空间位置，具体计算方法可参考文献［ 19 19 ］，并以此来评估风扰动对望远镜的影响。 在风扰动模拟时，将主镜与次镜及桁架分开考虑，这样有利于详细分析主镜与次镜及桁架之间的相互耦合关系。在主镜单独作用风载荷情况下，主、次镜镜面面形精度变化如 图7 图7 所示，光轴角度变化如 图8 图8 所示，主、次镜相对位置偏差如 图9 图9 所示。 1;2; 1;2; 在主镜单独作用10 m/s和15 m/s风扰动载荷情况下，主镜的面形精度影响很大，最差面形精度分别达到了45 nm和70 nm，而次镜面形精度基本不受影响；主、次镜光轴角度偏差最差时分别达到0.15″和0.25″，对相对位置偏差几乎没有影响。 在次镜及桁架单独作用风载荷情况下，主、次镜镜面面形精度变化如 图10 图10 所示，光轴角度变化如 图11 图11 所示，主、次镜相对位置偏差如 图12 图12 所示。 1;2; 1;2; 在次镜及桁架单独作用风载荷情况下，主镜面形精度受影响相对较小，在整个仿真时段均在5 nm左右；而次镜面形精度变化同样很小，原因在于次镜口径较小；主次镜光轴角度较前一种工况增加幅度很大，初始阶段最大偏差接近8″~10″，剩余阶段在4.5″以内，但波动频率明显提高；同样主、次镜 y 向相对位置偏差增加幅度很大，分别达到了0.1 mm和0.2 mm。对 图11 图11 和 图12 图12 （b）所示数据进行FFT变换得到频谱如 图13 图13 所示。从频谱曲线中可以看出，风扰动载荷对望远镜的影响集中在低频部分，约0~5 Hz区域，高频分量的影响较小。 1;2; 从以上两种情况下的分析结果对比中可以得出以下结论： （1）由于次镜口径较小，相对刚度较大，风扰动对次镜面形精度影响较小； （2）主镜面形精度变化主要来源于作用在主镜镜面的风载荷，受次镜及桁架的影响较小，即主、次镜相互耦合较小；从望远镜结构特征角度看，主镜连同主镜室以及桁架分别连接在四通下表面和上表面，而四通本身刚度很大，因此相互耦合不明显的结果是合理的； （3）次镜及桁架单元作用风载荷时，主、次镜光轴角度偏差以及相对位置偏差明显增大，将会导致望远镜光学系统成像抖动和模糊。 5 结　论 1 本文针对大型地基光学望远镜系统性能受风扰动影响显著的问题，研究了风扰动的时程模拟方法并对望远镜在风扰动下的性能进行预测。利用有限元方法及模态变换，建立了望远镜结构动力学模型。基于二维空间随机场，提出望远镜圆顶内风速的时程模拟方法，此方法可以表达空间和时间频率范围内几乎连续的随机场。对2 m口径望远镜的仿真分析结果表明：在外界平均风速10 m/s和15 m/s情况下，随机风扰动作用在望远镜主镜上会造成主镜面形变差，分别造成最大接近45 nm和70 nm的面形误差；作用在次镜及桁架上的风扰动主要造成低频0~5 Hz光轴角度误差及主、次镜相对位置偏差；主、次镜之间相互耦合效应较小。 参考文献 BELY P Y . The Design and Construction of Large Optical Telescopes ［M］. New York ： Springer ， 2003 . doi: 10.1007/b97612 http://dx.doi.org/10.1007/b97612 The Design and Construction of Large Optical Telescopes ANDERSEN T ， ENMARK A . Integrated Modeling of Telescopes ［M］. New York， NY ： Springer New York ， 2011 . doi: 10.1007/978-1-4614-0149-0 http://dx.doi.org/10.1007/978-1-4614-0149-0 Integrated Modeling of Telescopes MACMYNOWSKI D G ， BLAUROCK C ， ANGELI G Z ， et al . Modeling wind-buffeting of the thirty meter telescope ［C］. SPIE Proceedings" ， " Modeling ， Systems Engineering ， and Project Management for Astronomy II. Orlando ， Florida ， USA. SPIE ， 2006 ： 6271 （ 01 ）： M2710 - M2710 . doi: 10.1117/12.672706 http://dx.doi.org/10.1117/12.672706 Modeling wind-buffeting of the thirty meter telescope KAN F W ， EGGERS D W . Wind vibration analyses of Giant Magellan Telescope ［C］. SPIE Proceedings" ， " Modeling ， Systems Engineering ， and Project Management for Astronomy II. Orlando ， Florida ， USA. SPIE ， 2006 . doi: 10.1117/12.672685 http://dx.doi.org/10.1117/12.672685 Wind vibration analyses of Giant Magellan Telescope ERM T M ， ANGELI G Z . TMT wind model validation with measurements on Keck and Gemini ［C］. SPIE Proceedings" ， " Modeling ， Systems Engineering ， and Project Management for Astronomy II. Orlando ， Florida ， USA. SPIE ， 2006 . doi: 10.1117/12.670002 http://dx.doi.org/10.1117/12.670002 TMT wind model validation with measurements on Keck and Gemini POTTEBAUM T ， MACMYNOWSKI D G . Buffeting of large telescopes： wind-tunnel measurements of the flow inside a generic enclosure ［J］. Journal of Fluids and Structures ， 2006 ， 22 （ 1 ）： 3 - 19 . doi: 10.1016/j.jfluidstructs.2005.08.005 http://dx.doi.org/10.1016/j.jfluidstructs.2005.08.005 Buffeting of large telescopes： wind-tunnel measurements of the flow inside a generic enclosure VOGIATZIS K ， SEGURSON A ， ANGELI G Z . Estimating the effect of wind loading on extremely large telescope performance using computational fluid dynamics ［C］. SPIE Astronomical Telescopes + Instrumentation. Proc SPIE 5497， Modeling and Systems Engineering for Astronomy， Glasgow， United Kingdom ， 2004 ， 5497 ： 311 - 320 . doi: 10.1117/12.551477 http://dx.doi.org/10.1117/12.551477 Estimating the effect of wind loading on extremely large telescope performance using computational fluid dynamics MACMYNOWSKI D G ， ANDERSEN T . Wind buffeting of large telescopes ［J］. Applied Optics ， 2010 ， 49 （ 4 ）： 625 . doi: 10.1364/ao.49.000625 http://dx.doi.org/10.1364/ao.49.000625 Wind buffeting of large telescopes MACMYNOWSKI D G ， VOGIATZIS K ， ANGELI G Z ， et al . Wind loads on ground-based telescopes ［J］. Applied Optics ， 2006 ， 45 （ 30 ）： 7912 - 7923 . doi: 10.1364/ao.45.007912 http://dx.doi.org/10.1364/ao.45.007912 Wind loads on ground-based telescopes GLAESE R M ， SHEEHAN M . Vibration mitigation for wind-induced jitter for the giant Magellan telescope ［C］. SPIE Astronomical Telescopes + Instrumentation. Proc SPIE 8444， Ground-Based and Airborne Telescopes IV ， Amsterdam， Netherlands ， 2012 ， 8444 ： 300 - 314 . doi: 10.1117/12.926337 http://dx.doi.org/10.1117/12.926337 Vibration mitigation for wind-induced jitter for the giant Magellan telescope 杨德华 ， 徐灵哲 ， 徐欣圻 . 大型光学天文望远镜风载作用分析 ［J］. 光学技术 ， 2009 ， 35 （ 3 ）： 342 - 346 . doi: 10.3321/j.issn:1002-1582.2009.03.012 http://dx.doi.org/10.3321/j.issn:1002-1582.2009.03.012 大型光学天文望远镜风载作用分析 YANG D H ， XU L Z ， XU X Q . Wind disturbance study on large astronomical telescopes ［J］. Optical Technique ， 2009 ， 35 （ 3 ）： 342 - 346 . （in Chinese） . doi: 10.3321/j.issn:1002-1582.2009.03.012 http://dx.doi.org/10.3321/j.issn:1002-1582.2009.03.012 Wind disturbance study on large astronomical telescopes 周超 ， 杨洪波 ， 吴小霞 ， 等 . 1.2 m望远镜风载作用分析 ［J］. 红外与激光工程 ， 2011 ， 40 （ 5 ）： 889 - 893 . doi: 10.3969/j.issn.1007-2276.2011.05.023 http://dx.doi.org/10.3969/j.issn.1007-2276.2011.05.023 1.2 m望远镜风载作用分析 ZHOU C ， YANG H B ， WU X X ， et al . Influence of wind loading on the 1.2 m telescope ［J］. Infrared and Laser Engineering ， 2011 ， 40 （ 5 ）： 889 - 893 . （in Chinese） . doi: 10.3969/j.issn.1007-2276.2011.05.023 http://dx.doi.org/10.3969/j.issn.1007-2276.2011.05.023 Influence of wind loading on the 1.2 m telescope 潘年 ， 马文礼 ， 黄金龙 . 地基望远镜风载数值分析 ［J］. 红外与激光工程 ， 2015 ， 44 （ 1 ）： 134 - 140 . doi: 10.3969/j.issn.1007-2276.2015.01.024 http://dx.doi.org/10.3969/j.issn.1007-2276.2015.01.024 地基望远镜风载数值分析 PAN N ， MA W L ， HUANG J L . Numerical analysis of wind load on ground-based telescope ［J］. Infrared and Laser Engineering ， 2015 ， 44 （ 1 ）： 134 - 140 . （in Chinese） . doi: 10.3969/j.issn.1007-2276.2015.01.024 http://dx.doi.org/10.3969/j.issn.1007-2276.2015.01.024 Numerical analysis of wind load on ground-based telescope 徐江海 ， 宫雪非 . 利用数值风洞实验进行四种天文圆顶的风载研究 ［J］. 光学学报 ， 2015 ， 35 （ 5 ）： 0501005 . doi: 10.3788/aos201535.0501005 http://dx.doi.org/10.3788/aos201535.0501005 利用数值风洞实验进行四种天文圆顶的风载研究 XU J H ， GONG X F . Study of wind load on four astronomical enclosures based on numerical wind tunnel experiments ［J］. Acta Optica Sinica ， 2015 ， 35 （ 5 ）： 0501005 . （in Chinese） . doi: 10.3788/aos201535.0501005 http://dx.doi.org/10.3788/aos201535.0501005 Study of wind load on four astronomical enclosures based on numerical wind tunnel experiments 宋玉鹏 ， 陈建兵 ， 彭勇波 . 二维空间脉动风场波数-频率联合功率谱表达的FFT模拟 ［J］. 振动工程学报 ， 2020 ， 33 （ 4 ）： 660 - 666 . doi: 10.16385/j.cnki.issn.1004-4523.2020.04.003 http://dx.doi.org/10.16385/j.cnki.issn.1004-4523.2020.04.003 二维空间脉动风场波数-频率联合功率谱表达的FFT模拟 SONG Y P ， CHEN J B ， PENG Y B . Simulation of fluctuating wind field in two-spatial dimensions by FFT implementation of wavenumber-frequency joint power spectrum ［J］. Journal of Vibration Engineering ， 2020 ， 33 （ 4 ）： 660 - 666 . （in Chinese） . doi: 10.16385/j.cnki.issn.1004-4523.2020.04.003 http://dx.doi.org/10.16385/j.cnki.issn.1004-4523.2020.04.003 Simulation of fluctuating wind field in two-spatial dimensions by FFT implementation of wavenumber-frequency joint power spectrum 刘章军 ， 叶永友 ， 刘增辉 . 脉动风速连续随机场的降维模拟 ［J］. 工程力学 ， 2018 ， 35 （ 11 ）： 8 - 16 . doi: 10.6052/j.issn.1000-4750.2017.07.0570 http://dx.doi.org/10.6052/j.issn.1000-4750.2017.07.0570 脉动风速连续随机场的降维模拟 LIU Z J ， YE Y Y ， LIU Z H . Simulation of fluctuating wind velocity continuous stochastic field by dimension reduction approach ［J］. Engineering Mechanics ， 2018 ， 35 （ 11 ）： 8 - 16 . （in Chinese） . doi: 10.6052/j.issn.1000-4750.2017.07.0570 http://dx.doi.org/10.6052/j.issn.1000-4750.2017.07.0570 Simulation of fluctuating wind velocity continuous stochastic field by dimension reduction approach BENOWITZ B A ， DEODATIS G . Simulation of wind velocities on long span structures： a novel stochastic wave based model ［J］. Journal of Wind Engineering and Industrial Aerodynamics ， 2015 ， 147 ： 154 - 163 . doi: 10.1016/j.jweia.2015.10.004 http://dx.doi.org/10.1016/j.jweia.2015.10.004 Simulation of wind velocities on long span structures： a novel stochastic wave based model MACMYNOWSKI D G ， ANGELI G Z ， VOGIATZIS K ， et al . Parametric modeling and control of telescope wind-induced vibration ［C］. SPIE Astronomical Telescopes+Instrumentation. Proc SPIE 5497， Modeling and Systems Engineering for Astronomy， Glasgow， United Kingdom ， 2004 ， 5497 ： 266 - 277 . doi: 10.1117/12.552282 http://dx.doi.org/10.1117/12.552282 Parametric modeling and control of telescope wind-induced vibration 曹玉岩 ， 王建立 ， 陈涛 ， 等 . 基于Hexapod平台的地基大型光学望远镜失调误差主动补偿 ［J］. 光学 精密工程 ， 2020 ， 28 （ 11 ）： 2452 - 2465 . doi: 10.37188/OPE.20202811.2452 http://dx.doi.org/10.37188/OPE.20202811.2452 基于Hexapod平台的地基大型光学望远镜失调误差主动补偿 CAO Y Y ， WANG J L ， CHEN T ， et al . Active compensation of aberration for large ground-based telescope based on Hexapod platform ［J］. Opt. Precision Eng. ， 2020 ， 28 （ 11 ）： 2452 - 2465 . （in Chinese） . doi: 10.37188/OPE.20202811.2452 http://dx.doi.org/10.37188/OPE.20202811.2452 Active compensation of aberration for large ground-based telescope based on Hexapod platform 10.37188/OPE.20223004.0442 TP273 A 1004-924X（2022）04-0442-13 曹玉岩,王建立,王志臣等.基于二维随机场的地基大型光学望远镜风扰动时程模拟与性能预测[J].光学精密工程,2022,30(04):442-454. 曹玉岩,王建立,王志臣等.基于二维随机场的地基大型光学望远镜风扰动时程模拟与性能预测[J].光学精密工程,2022,30(04):442-454. CAO Yuyan,WANG Jianli,WANG Zhichen,et al.Wind disturbance simulation and optical performance prediction for large ground-based optical telescopes based on two-dimensional stochastic fields[J].Optics and Precision Engineering,2022,30(04):442-454. CAO Yuyan,WANG Jianli,WANG Zhichen,et al.Wind disturbance simulation and optical performance prediction for large ground-based optical telescopes based on two-dimensional stochastic fields[J].Optics and Precision Engineering,2022,30(04):442-454. 国家自然科学基金资助项目（No.11803035，No.U2031126） 2021-06-22 2021-08-17 2022-02-25 2022-02-26 光学精密工程 4 30