CN106989694B - 一种考虑膜片褶皱的索膜电极面形面精度的估计方法 - Google Patents

一种考虑膜片褶皱的索膜电极面形面精度的估计方法 Download PDF

Info

Publication number
CN106989694B
CN106989694B CN201710347619.3A CN201710347619A CN106989694B CN 106989694 B CN106989694 B CN 106989694B CN 201710347619 A CN201710347619 A CN 201710347619A CN 106989694 B CN106989694 B CN 106989694B
Authority
CN
China
Prior art keywords
normal distribution
diaphragm
point
triangular
plane
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201710347619.3A
Other languages
English (en)
Other versions
CN106989694A (zh
Inventor
张逸群
朱日升
杜敬利
杨东武
张树新
李申
李娜
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Xian University of Electronic Science and Technology
Original Assignee
Xian University of Electronic Science and Technology
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Xian University of Electronic Science and Technology filed Critical Xian University of Electronic Science and Technology
Priority to CN201710347619.3A priority Critical patent/CN106989694B/zh
Publication of CN106989694A publication Critical patent/CN106989694A/zh
Application granted granted Critical
Publication of CN106989694B publication Critical patent/CN106989694B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
    • G01B11/00—Measuring arrangements characterised by the use of optical techniques
    • G01B11/24—Measuring arrangements characterised by the use of optical techniques for measuring contours or curvatures
    • G01B11/2433—Measuring arrangements characterised by the use of optical techniques for measuring contours or curvatures for measuring outlines by shadow casting

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Aerials With Secondary Devices (AREA)
  • Measurement Of Length, Angles, Or The Like Using Electric Or Magnetic Means (AREA)

Abstract

本发明提供了一种考虑膜片褶皱的索膜电极面形面精度的估计方法，其步骤如下：首先，将索网节点、膜片测量点和三角形膜片投影至XOY平面，并在要求测量区域均匀布置离散点，删除不必要的离散点；其次，以各个膜片测量点离其三角形膜片平面的距离为褶皱波峰，生成方形二维正态分布；然后判断离散点是否落在膜片测量点的正态分布内，确定三角形膜片平面点和正态分布点的三维坐标；最终以索网节点、膜片测量点、三角形膜片平面点和正态分布点的三维坐标拟合索膜电极面的最佳拟合抛物面，计算其形面精度。本发明针对实际工程中含有褶皱膜片的索膜电极面的形面精度进行了合理地估算，其更能反映实际索膜电极面的形面精度，较准确地指导其形面调整工作。

Description

一种考虑膜片褶皱的索膜电极面形面精度的估计方法
技术领域
本发明属于计算形面精度技术领域，具体涉及一种考虑膜片褶皱的索膜电极面形面精度的估计方法。
背景技术
静电成形薄膜反射面天线是一种新型的星载可展开天线，其主要结构由周边桁架、基础电极、薄膜反射面、Astromesh索网、薄膜支撑环和边界拉索组成，基础电极铺设于Astromesh索网前索网上并施加高电势，镀有金属层的薄膜反射面为等效零势面，由于电势差产生静电力来张拉薄膜，从而使薄膜形成具有一定焦径比的高精度反射面。然而，在实际工程中，铺设在Astromesh索网前索网的基础电极三角形膜片常常由于粘贴误差导致的应力不均匀而发生褶皱。
通常，考虑膜片误差的索膜电极面形面精度是以膜片逼近抛物面的逼近误差来描述的，通过在膜片上分布的高斯积分点，对三角形膜片进行面积积分求解其逼近误差。实际工程中，若仍以此方法描述含膜片褶皱的索膜电极面的形面精度，势必不能真实反映电极面的形面精度，难以指导基础电极面的调整。目前，仍没有相关文献或者方法去描述含膜片褶皱的索膜电极面的形面精度，提出一种考虑膜片褶皱的索膜电极面形面精度的估计方法是亟待解决的问题。
发明内容
本发明的目的是解决上述含膜片褶皱的索膜电极面形面精度估算问题，提供一种考虑膜片褶皱的索膜电极面形面精度的估计方法，能够根据由摄影测量技术得到的索网节点和膜片测量点的三维坐标，假设三角形膜片的褶皱分布满足二维正态分布，进而以每个膜片测量点作为二维正态分布的波峰点，将落在三角形膜片投影面上的插值点分为正态分布点和三角形膜片平面点，最终以索网节点、膜片测量点、正态分布点和三角形膜片平面点估算电极面的形面精度。
本发明的技术方案是：一种考虑膜片褶皱的索膜电极面形面精度的估计方法，该方法具体包括以下步骤：
步骤1)首先，利用非接触式摄影测量技术，测量得到索膜电极面的索网节点、膜片测量点的三维坐标，其次将索网节点、膜片测量点和三角形膜片投影至投影面XOY内，然后利用MATLAB的meshgrid命令在索膜电极面要求测量的a×a m2方形区域内均匀生成离散点；其中，离散点在笛卡尔坐标系的x、y方向上的间隔都是bmm，b表示离散点的间距，a表示方形测量区域的边长；
步骤2)判断步骤1)均匀布置的离散点是否落在三角形膜片的投影面内，将不落在三角形膜片投影面内的离散点删去；相反，将落在三角形膜片投影面内的离散点保留，并记录该离散点落于哪一个三角形膜片投影面内；
步骤3)根据步骤1)测量得到的膜片测量点的三维坐标，可计算各个膜片测量点到其所在的三角形膜片平面的距离h，以h作为褶皱的波峰，生成边长为50mm的方形二维正态分布；
步骤4)判断投影至三角形膜片投影面内的离散点是否落在膜片测量点的正态分布范围内；若不在，将该离散点定义为三角形膜片平面点，根据三角形膜片的平面方程AX+BY+CZ+D＝0，然后将该三角形膜片平面点的Xplane、Yplane坐标代入其三角形膜片的平面方程，可求解出该三角形膜片平面点的Zplane坐标；其中A、B、C、D表示该三角形膜片平面方程的系数，可由该三角形膜片的三个顶点确定，(Xplane,Yplane,Zplane)为该离散点的三维坐标；若离散点落在膜片测量点的正态分布范围内，则将该离散点定义为正态分布点，根据正态分布函数和其膜片的平面方程确定该离散点的Zgauss坐标，Zgauss表示正态分布点的笛卡尔Z轴坐标；
步骤5)根据步骤4)得到的三角形膜片平面点和正态分布点的三维坐标，结合步骤1)测量得到的索网节点、膜片测量点的三维坐标，可以拟合确定索膜电极面的最佳拟合抛物面，最终计算索膜电极面的形面精度。
上述步骤3)中假设膜片的褶皱区域分布服从二维正态分布，然后分别以各个膜片测量点到其三角形膜片平面的距离h作为正态分布的波峰值，其中，h＝[h1,h2,…,hj]T,j＝1,2,…,n，n为膜片测量点的个数，j为膜片测量点的编号，hj表示第j个膜片测量点到其三角形膜片平面的距离，T表示矩阵的转置；最终，生成n个边长为50mm的方形二维正态分布，其具体步骤如下：
步骤3.1)将n个膜片测量点定义为其二维正态分布的波峰点，利用MATLAB软件的正态分布函数mvnpdf计算所有膜片测量点的二维正态分布波峰值pmax，得到n个边长为50mm的方形二维正态分布；
步骤3.2)根据公式
Figure GDA0002176710680000031
计算比例系数列向量α，其中，α＝[α1,α2,…αj]T，αj表示的意义为第j个膜片测量点其褶皱值与正态分布函数波峰值的比值。
上述步骤4)中处理离散点的方法，具体步骤如下：
步骤4.1)判断投影至三角形膜片投影面内的离散点是否落于膜片测量点边长为50mm的方形二维正态分布范围内，将落于二维正态分布内的离散点定义为正态分布点，记正态分布点的个数为Nguass；反之，将该离散点定义为三角形膜片平面点，记三角形膜片平面点的个数为Nplane；
步骤4.2)确定正态分布点的褶皱高度H＝[H1,H2,…,Hk]T,k＝1,2,…Nguass，k为正态分布点的编号，Hk表示第k个正态分布点的褶皱高度；首先计算正态分布点与其二维正态分布波峰点的x坐标差值列向量ΔX＝[Δx1,Δx2,…,Δxk]T和y坐标差值列向量ΔY＝[Δy1,Δy2,…,Δyk]T，其中，Δxk表示第k个正态分布点与其二维正态分布波峰点的x坐标差值，Δyk表示第k个正态分布点与其二维正态分布波峰点的x坐标差值；其次，通过MATLAB软件的正态分布函数mvnpdf计算该正态分布点的概率值p＝[p1,p2,…,pk]T，pk表示第k个正态分布点的概率值；然后根据步骤3.2)确定对应的比例系数α＝[α1,α2,…,αj]T，根据
Figure GDA0002176710680000041
确定正态分布点的褶皱高度；
步骤4.3)确定每一三角形膜片的平面方程；记三角形膜片的个数为m，假设第l个三角形膜片的平面方程可以表示为AlX+BlY+ClZ+Dl＝0，平面方程的Al、Bl、Cl、Dl系数由该三角形膜片的三个顶点确定，l为三角形膜片的编号；然后将三角形膜片平面点的Xplane、Yplane坐标代入其三角形膜片的平面方程，可求解该三角形膜片平面点的Zplane坐标；
步骤4.4)基于步骤4.3)求解得到的三角形膜片平面方程，即第l个三角形膜片的平面方程为AlX+BlY+ClZ+Dl＝0，进而，落在第l个三角形膜片的正态分布点的Zguass坐标可表示为Zguass＝(AlXguass+BlYguass+Dl)/(-Cl)+Hl，其中，(Xguass,Yguass,Zguass)为正态分布点的三维坐标，Hl表示为落在第l个三角形膜片的正态分布点的褶皱高度，l＝1,2,…m；最终，依次求解出m个三角形膜片内正态分布点的Zguass坐标。
上述步骤5)中根据索网节点、膜片测量点、三角形膜片平面点和正态分布点的三维坐标，拟合确定索膜电极面的最佳拟合抛物面，具体步骤如下：
步骤5.1)建立求解最佳拟合抛物面的优化模型：以分别绕x、y、z方向旋转的θx、θy、θz角度和沿x、y、z方向平移的平移量u、v、w为设计变量，即优化模型的设计变量s＝[θx,θy,θz,u,v,w]T，目标是使得以最佳拟合抛物面为评价基准的变形后抛物面的形面精度RMS(root-mean-square)最小，优化模型如下：
Find s＝[θx,θy,θz,u,v,w]T
Min
Figure GDA0002176710680000051
S.t. |u|,|v|,|w|≤D
|θx|,|θy|,|θz|≤θ
其中，N为坐标点数目；[newXi,newYi,newZi]T为最佳拟合抛物面第i个坐标点的三维坐标，i为坐标点的编号，F为天线的焦距，D为天线的物理口径，|u|,|v|,|w|分别为平移量u、v、w的绝对值，|θx|,|θy|,|θz|分别为旋转角θx、θy、θz的绝对值，θ为允许抛物面绕轴转动的角度上限值；
步骤5.2)确定最佳拟合抛物面的所有节点坐标[newX,newY,newZ]T的表达式；
步骤5.3)利用MATLAB软件的全局优化函数GrobalSearch求解优化模型，得到该优化模型的最优解。
本发明的技术优点在于：
1.本发明针对含有褶皱膜片的索膜电极面，提出了一种描述其真正形面精度的估算方法，该方法操作简单，可行性强，且相比传统的形面精度计算方法，更能真实反映出索膜电极面的形面精度，从而弥补了该方面的空缺。
2.本发明基于估算的形面精度，可以指导实际工程中含褶皱膜片的索膜电极面的形面调整工作，使得电极面的形面精度真正地满足设计要求。
以下将结合附图对本发明做进一步详细说明。
附图说明
图1为静电成形薄膜反射面天线结构示意图；
图2为索膜电极面的理论初始构型俯视图；
图3为索膜电极面的理论初始构型正视图；
图4为本发明的流程图；
图5为索膜电极面的索网节点、膜片测量点和膜片所处位置侧视图；
图6为索膜电极面的索网节点、膜片测量点和膜片所处位置俯视图；
图7为离散点的布置范围图；
图8为某一膜片测量点生成的正态分布图；
图9为三角形膜片平面点的分布图；
图10为正态分布点分布图；
图11为索膜电极面的变形后抛物面、理想抛物面及最佳拟合抛物面。
具体实施方式
静电成形薄膜反射面天线是一种新型的星载可展开天线，其主要结构由周边桁架、基础电极、薄膜反射面、Astromesh索网、薄膜支撑环和边界拉索组成，如图1所示。三角形的基础电极铺设在Astromesh索网的三角形网格上，形成索膜电极面，其结构的俯视图如图2所示，正视图如图3所示。
本发明提出一种考虑膜片褶皱的索膜电极面形面精度的估计方法，假设三角形膜片的褶皱区域分布满足二维的正态分布，进而，根据由摄影测量技术得到的索网节点和膜片测量点，以每个膜片测量点作为二维正态分布的波峰点，将落在三角形膜片投影面内的插值点分为正态分布点和膜片平面点，最终以索网节点、膜片测量点、正态分布点和三角形膜片平面点估算电极面的形面精度，其流程图如图4所示。其中，索膜电极面的索网节点、膜片测量点和三角形膜片分布的侧视图如图5所示，索网节点、膜片测量点分布的俯视图如图6所示，估计方法的详细步骤如下：
步骤1)利用非接触式摄影测量技术，测量得到索膜电极面的索网节点、膜片测量点的三维坐标，然后将索网节点、膜片测量点和三角形膜片投影至投影面XOY内，利用MATLAB的meshgrid命令在索膜电极面要求测量的a×a m2方形区域内均匀生成离散点，如图7所示；
步骤2)利用MATLAB编写程序，判断均匀布置的离散点是否落在三角形膜片投影面内，将不落在三角形膜片投影面上的离散点删去；相反，将落在三角形膜片投影面上的离散点保留，并记录该离散点落于哪一个三角形膜片投影面内；
步骤3)假设三角形膜片的褶皱区域分布服从二维正态分布，进而以各个膜片测量点到其三角形膜片平面的距离h＝[h1,h2,…,hj]T作为二维正态分布的波峰值，其中，j＝1,2,…,n。最终生成n个边长为50mm的方形二维正态分布，其具体步骤如下：
步骤3.1)分别将n个膜片测量点定义为其二维正态分布的波峰点，利用MATLAB的正态分布函数mvnpdf计算所有膜片测量点的二维正态分布波峰值pmax，生成n个边长为50mm的方形二维正态分布，图8为某一膜片测量点生成的二维正态分布图；
步骤3.2)计算得到比例系数列向量
Figure GDA0002176710680000071
步骤4)判断投影至三角形膜片投影面内的离散点是否落在膜片测量点生成的二维正态分布范围内；若不在，将该离散点定义为三角形膜片平面点，根据三角形膜片的平面方程，确定该离散点的Zplane坐标；若在二维正态分布范围内，将该离散点定义为正态分布点，根据正态分布函数和其三角形膜片的平面方程确定该离散点的Zgauss坐标，其具体的步骤如下：
步骤4.1)利用MATLAB编写程序，判断投影至三角形膜片投影面内的离散点是否落在步骤3)中膜片测量点生成的正态分布范围内，将落于正态分布内的离散点定义为正态分布点，记正态分布点的个数为Nguass，如图10所示的“+”点为正态分布点；反之，将该离散点定义为三角形膜片平面点，记膜片平面点的个数为Nplane，如图9所示的“*”点为三角形膜片平面点；
步骤4.2)确定步骤4.1)的正态分布点的褶皱高度H＝[H1,H2,…,Hk]T，其中，k＝1,2,…Nguass。首先，分别计算正态分布点到其二维正态分布波峰点的x、y坐标差值列向量ΔX＝[Δx1,Δx2,…,Δxk]T、ΔY＝[Δy1,Δy2,…,Δyk]T，通过MATLAB正态分布函数mvnpdf计算正态分布点的概率值p＝[p1,p2,…,pk]T；并由步骤3.2)确定的对应比例系数α＝[α1,α2,…,αj]T，根据
Figure GDA0002176710680000082
计算出正态分布点的褶皱高度；
步骤4.3)确定每一片三角形膜片的平面方程；记三角形膜片的个数为m，假设第l个三角形膜片的平面方程可以表示为AlX+BlY+ClZ+Dl＝0，平面方程的Al、Bl、Cl、Dl系数可以由该三角形膜片的三个顶点确定；然后将三角形膜片平面点的Xplane、Yplane坐标代入三角形膜片的平面方程，求解该三角形膜片平面点的Zplane坐标；
步骤4.4)基于步骤4.3)求解得到的m个三角形膜片的平面方程，落在第l个三角形膜片投影面内的正态分布点的Zguass坐标可表示为
Zguass＝(AlXguass+BlYguass+Dl)/(-Cl)+Hl
其中，Hl表示为落在第l个三角形膜片投影面内的正态分布点褶皱高度，l＝1,2,…m。最终，依次求解出m个三角形膜片投影面内正态分布点的Zguass坐标。
步骤5)根据索网节点、膜片测量点、三角形膜片平面点和正态分布点的三维坐标，拟合索膜电极面的最佳拟合抛物面，索膜电极面的最佳拟合抛物面与其理想抛物面、变形后抛物面的位置关系如图11所示，最终，计算索膜电极面的形面精度RMS，其具体步骤如下：
步骤5.1)建立求解最佳拟合抛物面的优化模型；以绕x、y、z方向旋转θx、θy、θz角度和x、y、z方向的平移量u、v、w为设计变量，目标是使得以最佳拟合抛物面为评价基准的变形后抛物面的形面精度RMS最小，优化模型如下：
Find s＝[θx,θy,θz,u,v,w]T
Min
S.t. |u|,|v|,|w|≤D
|θx|,|θy|,|θz|≤θ
步骤5.2)确定最佳拟合抛物面的所有节点坐标[newX,newY,newZ]T的表达式。设索网节点、膜片测量点、三角形膜片平面点和正态分布点的三维坐标为[X,Y,Z]T，则最佳拟合抛物面的所有节点的三维坐标[newX,newY,newZ]T可表示为
[newX,newY,newZ]T＝([X,Y,Z]T-[u,v,w]T)·Rz·Ry·Rx
其中，
Figure GDA0002176710680000091
Rx、Ry、Rz分别表示绕x、y、z旋转的旋转矩阵，[u,v,w]T为平移矩阵。
步骤5.3)利用MATLAB的全局优化函数GrobalSearch求解优化模型，可得到该优化模型的最优解s*为
Figure GDA0002176710680000092
进而，以最佳拟合抛物面为评价基准，变形后抛物面的形面精度为
Figure GDA0002176710680000093
综上，本发明具备如下技术优点：
1.本发明针对含有褶皱膜片的索膜电极面，提出了一种描述其真正形面精度的估算方法，该方法操作简单，可行性强，且相比传统的形面精度计算方法，更能真实反映出索膜电极面的形面精度，弥补了该方面的空缺。
2.本发明基于估算的形面精度，可以指导实际工程中含褶皱膜片的索膜电极面的形面调整工作，使得电极面的形面精度真正地满足设计要求。
本实施方式中没有详细叙述的部分属本行业的公知的常用手段，这里不一一叙述。以上例举仅仅是对本发明的举例说明，并不构成对本发明的保护范围的限制，凡是与本发明相同或相似的设计均属于本发明的保护范围之内。

Claims (2)

1.一种考虑膜片褶皱的索膜电极面形面精度的估计方法，其特征在于，该方法具体包括以下步骤：
步骤1)首先，利用非接触式摄影测量技术，测量得到索膜电极面的索网节点、膜片测量点的三维坐标，其次将索网节点、膜片测量点和三角形膜片投影至投影面XOY内，然后利用MATLAB软件的meshgrid命令在索膜电极面要求测量的a×a m2方形区域内均匀生成离散点；其中，离散点在笛卡尔坐标系的x、y方向上的间隔都是b mm，b表示离散点的间距，a表示方形测量区域的边长；
步骤2)判断步骤1)均匀布置的离散点是否落在三角形膜片的投影面内，将不落在三角形膜片投影面内的离散点删去；相反，将落在三角形膜片投影面内的离散点保留，并记录该离散点落于哪一个三角形膜片投影面内；
步骤3)根据步骤1)膜片测量点的三维坐标，计算膜片测量点到其所在的三角形膜片平面的距离hj，j为膜片测量点的编号，hj表示第j个膜片测量点到其所在三角形膜片平面的距离；将hj作为以膜片测量点为中心，边长为50mm的正方形区域的波峰值，生成第j个膜片测量点的方形二维正态分布；记
h＝[h1,h2,…,hj,…,hn-1,hn]T,(j＝1,2,…,n)为n个膜片测量点到其所在三角形膜片平面的距离列向量，该向量维度为n×1，n为膜片测量点的个数，T表示矩阵的转置；最终，生成n个边长为50mm的方形二维正态分布，其具体步骤如下：
步骤3.1)将n个膜片测量点定义为其二维正态分布的波峰点，利用MATLAB软件的正态分布函数mvnpdf计算该n个所有膜片测量点的二维正态分布波峰值，记第j个膜片测量点的二维正态分布波峰值为pmax,j，则n个膜片测量点的二维正态分布波峰值组成的列向量表示为pmax＝[p1,p2,…,pj,…,pn-1,pn]T，得到n个边长为50mm的方形二维正态分布；
步骤3.2)根据公式
Figure FDA0002176710670000021
计算比例系数列向量α，其中，α＝[α1,α2,…αj,…αn-1,αn]T，αj表示的意义为第j个膜片测量点其褶皱值与正态分布函数波峰值的比值；
步骤4)判断投影至三角形膜片投影面内的离散点是否落在膜片测量点的正态分布范围内；若不在，将该离散点定义为三角形膜片平面点，根据三角形膜片的平面方程AX+BY+CZ+D＝0，然后将该三角形膜片平面点的Xplane、Yplane坐标代入其三角形膜片的平面方程，求解出该三角形膜片平面点的Zplane坐标；其中A、B、C、D表示该三角形膜片平面方程的系数，由该三角形膜片的三个顶点确定，(Xplane,Yplane,Zplane)为该离散点的三维坐标；若离散点落在膜片测量点的正态分布范围内，则将该离散点定义为正态分布点，根据正态分布函数和其膜片的平面方程确定该离散点的Zgauss坐标，Zgauss表示正态分布点的笛卡尔Z轴坐标；
步骤5)根据步骤4)得到的三角形膜片平面点和正态分布点的三维坐标，结合步骤1)测量得到的索网节点、膜片测量点的三维坐标，拟合确定索膜电极面的最佳拟合抛物面，最终计算索膜电极面的形面精度；
所述步骤5)中根据索网节点、膜片测量点、三角形膜片平面点和正态分布点的三维坐标，拟合确定索膜电极面的最佳拟合抛物面，具体步骤如下：
步骤5.1)建立求解最佳拟合抛物面的优化模型：以分别绕x、y、z方向旋转的θx、θy、θz角度和沿x、y、z方向平移的平移量u、v、w为设计变量，即优化模型的设计变量s＝[θx,θy,θz,u,v,w]T，目标是使得以最佳拟合抛物面为评价基准的变形后抛物面的形面精度RMS最小，优化模型如下：
Find s＝[θx,θy,θz,u,v,w]T
Min
Figure FDA0002176710670000031
S.t.|u|,|v|,|w|≤D
|θx|,|θy|,|θz|≤θ
其中，N为坐标点数目；[newXi,newYi,newZi]T为最佳拟合抛物面第i个坐标点的三维坐标，i为坐标点的编号，F为天线的焦距，D为天线的物理口径，|u|,|v|,|w|分别为平移量u、v、w的绝对值，|θx|,|θy|,|θz|分别为旋转角θx、θy、θz的绝对值，θ为允许抛物面绕轴转动的角度上限值；
步骤5.2)确定最佳拟合抛物面的所有节点坐标[newX,newY,newZ]T的表达式；
步骤5.3)利用MATLAB软件的全局优化函数GrobalSearch求解优化模型，得到该优化模型的最优解。
2.根据权利要求1所述的一种考虑膜片褶皱的索膜电极面形面精度的估计方法，其特征在于，步骤4)中处理离散点的方法，具体步骤如下：
步骤4.1)判断投影至三角形膜片投影面内的离散点是否落于膜片测量点边长为50mm的方形二维正态分布范围内，将落于二维正态分布内的离散点定义为正态分布点，其中，落在同一个二维正态分布范围内的离散点，其褶皱值与正态分布函数波峰值的比值相同，记正态分布点的个数为Nguass，则正态分布点褶皱值与其正态分布函数波峰值的比值列向量表示为
Figure FDA0002176710670000041
反之，将该离散点定义为三角形膜片平面点，记三角形膜片平面点的个数为Nplane；
步骤4.2)确定正态分布点的褶皱高度
k为正态分布点的编号，Hk表示第k个正态分布点的褶皱高度；首先计算正态分布点与其二维正态分布波峰点的x坐标差值列向量
Figure FDA0002176710670000043
和y坐标差值列向量
Figure FDA0002176710670000044
其中，Δxk表示第k个正态分布点与其二维正态分布波峰点的x坐标差值，Δyk表示第k个正态分布点与其二维正态分布波峰点的y坐标差值；其次，通过MATLAB软件的正态分布函数mvnpdf计算该正态分布点的概率值pk表示第k个正态分布点的概率值；然后根据步骤4.1)确定对应的比例系数
Figure FDA0002176710670000051
根据
Figure FDA0002176710670000052
确定正态分布点的褶皱高度；
步骤4.3)确定每一三角形膜片的平面方程；记三角形膜片的个数为m，假设第l个三角形膜片的平面方程表示为AlX+BlY+ClZ+Dl＝0，平面方程的Al、Bl、Cl、Dl系数由该三角形膜片的三个顶点确定，l为三角形膜片的编号；然后将三角形膜片平面点的Xplane、Yplane坐标代入其三角形膜片的平面方程，求解该三角形膜片平面点的Zplane坐标；
步骤4.4)基于步骤4.3)求解得到的三角形膜片平面方程，即第l个三角形膜片的平面方程为AlX+BlY+ClZ+Dl＝0，进而，落在第l个三角形膜片的正态分布点的Zguass坐标表示为Zguass＝(AlXguass+BlYguass+Dl)/(-Cl)+Hl，其中，(Xguass,Yguass,Zguass)为正态分布点的三维坐标，Hl表示为落在第l个三角形膜片的正态分布点的褶皱高度，l＝1,2,…m；最终，依次求解出m个三角形膜片内正态分布点的Zguass坐标。
CN201710347619.3A 2017-05-17 2017-05-17 一种考虑膜片褶皱的索膜电极面形面精度的估计方法 Active CN106989694B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710347619.3A CN106989694B (zh) 2017-05-17 2017-05-17 一种考虑膜片褶皱的索膜电极面形面精度的估计方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710347619.3A CN106989694B (zh) 2017-05-17 2017-05-17 一种考虑膜片褶皱的索膜电极面形面精度的估计方法

Publications (2)

Publication Number Publication Date
CN106989694A CN106989694A (zh) 2017-07-28
CN106989694B true CN106989694B (zh) 2020-01-14

Family

ID=59420461

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710347619.3A Active CN106989694B (zh) 2017-05-17 2017-05-17 一种考虑膜片褶皱的索膜电极面形面精度的估计方法

Country Status (1)

Country Link
CN (1) CN106989694B (zh)

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP4119817B2 (ja) * 2003-10-31 2008-07-16 三菱重工業株式会社 パラメータ実長展開装置、方法及びそのプログラム
CN103678810A (zh) * 2013-12-17 2014-03-26 西安电子科技大学 一种静电成形薄膜天线的电极布局方法
CN104123421A (zh) * 2014-07-31 2014-10-29 西安电子科技大学 一种基于机电场耦合的静电成形薄膜反射面形态设计方法
CN104748714A (zh) * 2014-08-15 2015-07-01 西安电子科技大学 一种星载索网可展开天线形面精度测量方法
CN104866666A (zh) * 2015-05-20 2015-08-26 西安电子科技大学 一种基于有限元模型修正的索网反射面形面精度的鲁棒调整方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP4119817B2 (ja) * 2003-10-31 2008-07-16 三菱重工業株式会社 パラメータ実長展開装置、方法及びそのプログラム
CN103678810A (zh) * 2013-12-17 2014-03-26 西安电子科技大学 一种静电成形薄膜天线的电极布局方法
CN104123421A (zh) * 2014-07-31 2014-10-29 西安电子科技大学 一种基于机电场耦合的静电成形薄膜反射面形态设计方法
CN104748714A (zh) * 2014-08-15 2015-07-01 西安电子科技大学 一种星载索网可展开天线形面精度测量方法
CN104866666A (zh) * 2015-05-20 2015-08-26 西安电子科技大学 一种基于有限元模型修正的索网反射面形面精度的鲁棒调整方法

Also Published As

Publication number Publication date
CN106989694A (zh) 2017-07-28

Similar Documents

Publication Publication Date Title
US10412365B2 (en) Calibration arrangement for structured light system using a tele-centric lens
CN110335301B (zh) 一种基于激光雷达和结构光的三维扫描方法及装置
CN106199605B (zh) 风场误差修正方法
CN106229605A (zh) 一种基于数学建模的大型相控阵天线精准安装方法
CN107576275A (zh) 一种利用摄影测量技术对充气结构进行应变场测量的方法
CN107271957A (zh) 基于tdoa和toa的室内三维定位方法
Yuan et al. Direct root-mean-square error for surface accuracy evaluation of large deployable mesh reflectors
Ulmer et al. Slope measurements of parabolic dish concentrators using color-coded targets
Guo et al. An adaptive sampling methodology for measuring blades with CMM based on dominant feature points
CN114444216A (zh) 基于数值模拟的高空条件下飞行器姿态控制方法及系统
Kumar et al. A novel real-time DIC-FPGA-based measurement method for dynamic testing of light and flexible structures
CN106989694B (zh) 一种考虑膜片褶皱的索膜电极面形面精度的估计方法
Amoignon et al. Study of parameterizations in the project CEDESA
CN112525145B (zh) 一种飞机降落相对姿态动态视觉测量方法及系统
Marinić-Kragić et al. Adaptive re-parameterization based on arbitrary scalar fields for shape optimization and surface fitting
WO2019194114A1 (ja) 処理装置、ｃａｄモデルの特徴部分検出方法及びプログラム
CN111211424B (zh) 一种适用于ska天线的主面面型测量调整方法及装置
CN111366149A (zh) 一种基于月相的太阳矢量方向的提取方法
Khare et al. Flexible DIC System for Wind Tunnel Testing of Lightweight Structures
Zheng et al. A new large-scale posture measurement system based on a six-laser tracer multilateral method
CN113124821A (zh) 一种基于曲面镜和平面镜的结构测量方法
CN110426562B (zh) 基于分层搜索和距离空间投影的高精度闪电三维定位方法
Chen et al. Coordinate stitching measurement of highly steep freeform surfaces
Willis et al. A quadratic basis function, quadratic geometry, high order panel method
Arcondoulis et al. Acoustic beamforming array design methods over irregular shaped areas

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant