CN106934147B - 一种空间桁架索网天线等几何找形方法 - Google Patents
一种空间桁架索网天线等几何找形方法 Download PDFInfo
- Publication number
- CN106934147B CN106934147B CN201710138544.8A CN201710138544A CN106934147B CN 106934147 B CN106934147 B CN 106934147B CN 201710138544 A CN201710138544 A CN 201710138544A CN 106934147 B CN106934147 B CN 106934147B
- Authority
- CN
- China
- Prior art keywords
- cable
- antenna
- node
- truss
- net
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Aerials With Secondary Devices (AREA)
- Details Of Aerials (AREA)
Abstract
本发明公开了一种空间桁架索网天线等几何找形方法,能够在提高精确度的同时简单有效的对桁架索网天线进行找形;该方法为:根据设定的天线设计参数,基于天线拓扑关系,获得支撑桁架和索网当前节点坐标;根据设定的索段张力值和索网当前节点坐标,计算索段初始长度,获得索段当前力密度系数,在不考虑桁架变形的基础上,基于索段初始力密度系数,采用渐近迭代力密度法对索网进行找形分析;并对桁架和索网有限元模型进行静力平衡迭代;判断天线是否满足静力平衡误差,完成桁架索网找形。
Description
技术领域
本发明属于空间索网天线找形技术领域,具体涉及一种空间桁架索网天线等几何找形方法。
背景技术
近年来,航天科技的发展迫切需要大口径桁架索网天线,进而满足移动通信、射电天文、对地观测和军事侦察等需求。桁架索网天线主要由支撑桁架、主网、副网、张力阵、金属网组成。在天线结构设计与分析中,可将金属网对反射面的影响等效到主网中考虑,即将这类天线结构视为仅由支撑桁架和索网组成的系统。在桁架索网天线的反射面设计中,需要面临众多的问题与挑战,其中包括天线的结构刚度问题和天线的形面误差问题。桁架索网天线与一般结构的主要差异表现为柔性索网所使用的材料自身不具有刚度和形状,即在自然状态下不具有形状保持和承载能力,只有对索网中各索段施加预张力,才能使索网形成需要的形状,并具有一定刚度。
因此在天线设计阶段,需要对天线反射面进行找形,以提高反射面的形面精度,同时使得天线具有足够的刚度。目前主要有三种索网找形方法,第一种:Linkwitz和Schek提出的力密度法是迄今应用最广泛的索网找形方法之一。该方法的基本思想是通过预先给定各索段中的力和索段长度的比值(力密度),将几何非线性问题转换为线性问题,联立求解线性方程组,得到索网中各节点的坐标,由于力密度法没有考虑节点坐标变化对节点平衡的影响,得到的初始位形误差较大,天线反射面形面精度较低;第二种:基于力密度法,Morterolle等提出了迭代的力密度法。该方法的主要思路是在每个迭代步更新力密度系数,调整预张力,从而得到均匀的索段预张力分布。但该方法对反射面的网格构型有很大限制,无法获得满足工程精度需求的反射面形面精度;第三种:第三种找形方法是极小范数法,该方法将索段张力的平均值作为设计变量,索段最大拉力比作为目标函数,求解静不定平衡方程的极小范数解,从而得到所有索段预张力分布。与第一种力密度法相比,极小范数法得到的索段张力均匀性要差很多,并且极小范数法的求解区间有限,对于一些复杂索网结构会出现负拉力,得到脱离实际的结果,使得天线反射面找形任务失败。此外,如何提高桁架索网天线机构反射面的形面精度,便成为亟待解决的技术问题。
发明内容
有鉴于此,本发明提供了一种空间桁架索网天线等几何找形方法,提高了桁架索网天线机构反射面的形面精度。
实施本发明的具有步骤如下:
天线包括支撑桁架和索网,索网由索网节点和索段组成,索网节点分为前索网节点和后索网节点,该方法的具体步骤如下:
步骤1:根据设定的天线设计参数,基于天线拓扑关系,获得支撑桁架和索网当前节点坐标;
步骤301:根据索段当前力密度系数和索网节点当前坐标,采用力密度法建立索网节点线性平衡方程组,并对该方程组求解,获得所有索网节点的坐标;
步骤302:根据步骤301获得的索网节点当前坐标,计算索网索段当前长度和索段当前张力值基于步骤301获得的前索网节点坐标,采用抛物面天线公式znew=(x2+y2)/(4f),对前索网节点的z坐标进行修正;其中,(x,y,z)为修正前的前索网节点坐标,(x,y,znew)为修正后的前索网节点坐标,f为抛物面天线焦距;
步骤303:将修正前与修正后的前索网节点坐标误差进行计算,判断误差是否满足设计精度要求;若不满足设计精度要求,将修正后的(x,y,znew)作为前索网节点当前坐标,重复步骤301-303;若满足设计精度要求,根据修正后的索网节点坐标重新计算索段张力值,以索网节点当前坐标和索段当前张力值作为索网的渐近迭代力密度找形结果。
进一步包括:
步骤4:根据步骤303输出的索网节点当前坐标和索段当前张力值,更新索段当前长度和索段当前力密度系数以索段当前张力值索段当前长度和索段当前力密度系数作为迭代的初始值,依据步骤3的找形结果,建立桁架和索网非线性有限元模型,并对桁架和索网有限元模型进行静力平衡迭代,k的初始值取0;
步骤5:根据步骤4的迭代结果判断天线是否满足静力平衡误差:||Δu||<tole,Δu为静力平衡迭代获得的索段节点位移增量,tole为设定的收敛误差,||·||为向量的模,若满足静力平衡误差,直接输出步骤3的找形结果,完成桁架索网找形;若不满足静力平衡误差,根据步骤4静力平衡迭代结果,更新索段当前张力值和索段当前长度继而计算索段当前力密度系数并将更新后的索段力密度系数代入步骤3中,重复步骤3-5。
进一步地,步骤1的具体处理过程如下:
步骤101:根据设定的天线设计参数,采用最小二乘法对抛物面天线反射面进行圆弧拟合,获得拟合球的球面方程;
步骤102:基于天线拓扑关系,根据步骤101获得的拟合球的球面方程确定天线桁架下表面N个节点坐标P_down,该N个节点分别为桁架中心点和周边S个端点;
步骤103:设定桁架高度H,基于步骤102确定的天线桁架下表面N个节点坐标P_down,获得桁架上表面对应的N个节点坐标P_up,计算公式如下:
式(4)中,r为桁架下表面节点和拟合球心之间连线的矢量;
步骤104:基于步骤102确定的天线桁架下表面N个节点坐标P_down,将后索网设计成环向索网格式,前索网设计成三向索网格式,根据三向索网天线反射面形面精度估算公式,确定估算精度小于或等于设定精度的索网分段数M,根据桁架上表面节点坐标对后索网进行等分,获得后索网节点坐标P_rear;
步骤105:根据步骤104获得的后索网节点坐标P_rear,基于后索网节点和拟合球心之间的连线与抛物面天线的交点,建立如下非线性方程组:
通过对式(5)进行求解,获得所述交点坐标(xP_fornt,yP_fornt,zP_fornt),将所述交点作为前索网节点P_front,其中,(xo,yo,zo)为步骤101获得的拟合球面的球心坐标,(xP_rear,yP_rear,zP_rear)为步骤104获得的后索网节点坐标;最终获得支撑桁架当前节点坐标P_down、P_up和索网当前节点坐标P_rear、P_front。
进一步地,步骤4具体步骤如下:
步骤401:根据步骤3找形得到的索网当前节点坐标和步骤1获得的支撑桁架当前节点坐标,建立天线空间几何模型,对桁架的天线空间几何模型进行网格划分,生成绝对节点坐标单元,对索网进行控制点插入,生成等几何索单元;
步骤402:根据天线空间几何模型,建立天线有限元模型,根据连续介质力学理论计算单元刚度矩阵,并对单元刚度矩阵进行组装,获得有限元模型的非线性刚度矩阵K;
步骤403:根据步骤3找形得到的索段张力对索网施加预应力F,对天线有限元模型施加边界条件;
步骤404:建立系统静力学平衡方程KΔU=F,求解静力学平衡方程获得单元节点位移增量ΔU,其中K为步骤402获得有限元模型的非线性刚度矩阵,F为步骤403获得的预应力;
步骤405:根据步骤404获得的节点位移增量ΔU,判断节点位移增量ΔU是否收敛,若节点位移增量不收敛,根据单元节点位移增量更新索网当前节点坐标和支撑桁架当前节点坐标,重复步骤402到405,若节点位移增量收敛,则完成静力平衡迭代。
进一步地,索网天线为六棱柱模块索网天线或三棱柱模块索网天线;
进一步地,N=7,S=6。
1)本发明在不考虑桁架变形的基础上,基于索段当前的力密度系数,对索网进行找形分析,提高了桁架索网天线机构反射面的形面精度;本发明的找形方法简单易于实现,适用范围更广,不会出现负拉力。
2)本发明在找形过程中建立了天线索网等几何非线性有限元模型,为后续天线模态分析和展开动力学分析奠定了基础。
3)本发明通过建立桁架和索网非线性有限元模型,并对桁架和索网有限元模型进行静力平衡迭代,在迭代过程中,考虑了桁架弹性变形的桁架索网天线机构反射面进行找形,可以得到均匀的索段张力分布,获得较高的反射面形面精度,提高了天线反射面找形精确度,且方法简单有效。
附图说明
图1是空间桁架索网天线等几何找形方法的主流程图。
图2是设计的某模块化索网天线初始构型图。
图3是空间桁架索网天线支撑桁架机构和索网设计流程图。
图4是设计的某模块化索网天线支撑桁架结构示意图。
图5是采用渐近迭代力密度法对索网进行找形的流程图。
图6是对天线有限元模型进行静力平衡迭代的流程图。
图7是本发明找形方法应用于某模块化索网天线进行仿真结果示意图。
具体实施方式
下面结合附图并举实施例,对本发明进行详细描述。
本发明提供了一种空间桁架索网天线等几何找形方法,如图1所示,该方法包括如下步骤:
步骤101:针对空间桁架索网天线,设定天线设计参数、索网及桁架的材料参数、几何参数和拓扑关系;
步骤102:根据设定的天线设计参数,采用最小二乘法对抛物面天线反射面进行球面拟合,基于天线拓扑关系,获得支撑桁架和索网当前节点坐标;
本步骤102包含子步骤201-205,具体为:
步骤201:根据设定的天线设计参数,采用最小二乘法对抛物面天线反射面进行圆弧拟合,获得拟合球的球面方程;
步骤202:基于天线拓扑关系,根据步骤201获得的拟合球的球面方程确定天线桁架下表面N个节点坐标P_down,该N个节点分别为桁架中心点和周边S个端点;
步骤203:设定桁架高度H,基于步骤202确定的天线桁架下表面N个节点坐标P_down,获得桁架上表面对应的N个节点坐标P_up,计算公式如下:
式(4)中,r为桁架下表面节点和拟合球心之间连线的矢量;
步骤204:基于步骤202确定的天线桁架下表面N个节点坐标P_down,将后索网设计成环向索网格式,前索网设计成三向索网格式,根据三向索网天线反射面形面精度估算公式确定索网分段数M,根据桁架上表面节点坐标对后索网进行等分,获得后索网节点坐标P_rear;其中,δrms为索网估算的三向索网形面精度,l为分段后前索网最长索段的长度,f为抛物面天线焦距,C=7.872;
步骤205:根据步骤204获得的后索网节点坐标P_rear,基于后索网节点和拟合球心之间的连线与抛物面天线的交点,建立如下非线性方程组:
通过对式(5)进行求解,获得所述交点坐标(xP_fornt,yP_fornt,zP_fornt),将所述交点作为前索网节点P_front,其中,(xo,yo,zo)为步骤201获得的拟合球面的球心坐标,(xP_rear,yP_rear,zP_rear)为步骤204获得的后索网节点坐标;最终获得支撑桁架节点坐标(P_down,P_up)和索网节点坐标(P_rear,P_front)。
步骤103:设定索段当前张力值根据索网节点坐标计算索段当前长度基于索段当前张力值和索段当前长度获得索段当前力密度系数 其中,m和n分别为索段节点编号,j为索段编号,k为迭代次数,k取整数,k的初始值取0;
本步骤104包含子步骤401-404,具体为:
步骤401:根据步骤103获得的力密度系数和步骤205获得的前索网节点坐标,采用力密度法建立索网节点线性平衡方程组:
式(6)中,qj为前索网第j个索段的力密度系数,(xn,yn,zn)为第j个索段的一个端点坐标,(xm,ym,zm)为第j个索段的另一个端点坐标,s为任一节点连接的索段数;求解式(6),获得所有索网节点坐标;
步骤402:根据步骤401获得的前索网节点坐标,计算前索网索段长度lj和索段张力值Tj;
步骤403:根据步骤401获得的前索网节点的坐标,采用抛物面天线公式znew=(x2+y2)/(4f),对前索网连接点的z坐标进行修正;其中,(x,y,z)为前索网节点坐标,(x,y,znew)为修正后的前索网坐标;
步骤404:将修正前与修正后的前索网节点坐标误差进行计算,判断误差是否满足设计精度要求;若不满足设计精度要求,将修正后的(x,y,znew)替换步骤401中的前索网节点坐标,重复步骤401-404;若满足设计精度要求,执行步骤404,输出索网所有节点坐标和索段张力值,完成对索网的渐近迭代力密度找形。
步骤105:以索段当前张力值索段当前长度和索段当前力密度系数作为迭代的初始值,依据步骤104的找形结果,采用绝对节点坐标法和等几何分析法分别建立桁架和索网非线性有限元模型,并对桁架和索网有限元模型进行静力平衡迭代。
本步骤105包含子步骤501-505,具体为:
步骤501:根据步骤104找形得到的索网节点坐标和桁架节点坐标,建立天线空间几何模型,对桁架的天线空间几何模型进行网格划分,生成绝对节点坐标单元,对索网进行控制点插入,生成等几何索单元;
步骤502:根据天线空间几何模型,建立有限元模型,根据连续介质力学理论计算单元刚度矩阵,并对单元刚度矩阵进行组装,获得有限元模型的非线性刚度矩阵K;
步骤503:根据104找形得到的索段张力对索网施加预应力F,对天线有限元模型施加边界条件和载荷;
步骤504:建立系统静力学平衡方程KΔU=F,求解静力学平衡方程获得单元节点位移增量ΔU,其中K为步骤502获得有限元模型的非线性刚度矩阵,F为步骤503获得的预应力;
步骤505:根据步骤504获得的节点位移增量ΔU,判断节点位移增量ΔU是否收敛,若节点位移增量不收敛,根据单元节点位移增量更新索网节点坐标重复步骤502到505,若节点位移增量收敛,则完成静力平衡迭代。
本步骤105中,根据位移增量ΔU是否收敛确定完成静力平衡迭代,K代表有限元模型的非线性刚度矩阵,这两个量结合进行静力平衡迭代,反应了对桁架弹性变形的考虑。
步骤106:根据步骤105的迭代结果判断天线是否满足静力平衡误差:||Δu||<tole,Δu为索段节点位移增量,tole为设定的收敛误差,||·||为向量的模,若满足静力平衡误差,直接输出步骤104的找形结果,完成桁架索网找形;若不满足静力平衡误差,根据步骤105静力平衡迭代结果,更新索段当前张力值和索段当前长度继而计算索段当前力密度系数并将更新后的索段力密度系数代入步骤104中,重复步骤104-106。其中,j为索段编号,k为迭代次数,k取整数,k的初始值取0。
下面以六棱柱模块索网天线为对象,进行几何找形方法的描述。具体实施方式如下:
步骤101:针对六棱柱模块索网天线如图2所示,给定天线设计参数,如天线焦距f=13.76m,偏心距h=3.44m,索网为凯夫拉绳,杨氏模量为2e10Pa,泊松比为0,线径为1mm;桁架为碳纤维杆,杨氏模量为1.6e11Pa,泊松比为0.3,外径为22毫米,内径为20毫米;
步骤102:根据天线设计参数采用最小二乘法对抛物面天线反射面进行球面拟合,流程如图3所示,得到拟合球心坐标(0,-2.787,37.12),球面半径r=37.577m。设计得到支撑桁架如图4所示,其中桁架上表面和下表面七个节点坐标如表1所示:
表1桁架上表面和下表面七个节点坐标
步骤104:在不考虑桁架变形的基础上,采用渐近迭代力密度法对索网进行找形分析,具体流程如图5所示;
步骤105:依据步骤104的找形结果,采用绝对节点坐标法和等几何分析法分别建立桁架和索网非线性有限元模型,并对桁架和索网有限元模型进行静力平衡迭代,具体流程如图6所示;
步骤106:根据步骤105的迭代结果判断天线是否满足静力平衡误差:||Δu||<tole,一般设置tole=1e-6,若满足静力平衡误差,根据步骤104的找形结果输出索网节点位置(x,y,z)和索段张力值T,完成桁架索网找形,对图2所示的六棱柱模块索网天线进行找形得到的结果如图7所示;若不满足静力平衡误差,根据步骤105静力平衡迭代结果更新索段张力索段长度和索段力密度系数转至步骤104,进行下一次找形迭代;Δu为索段节点位移增量,tole为设定的收敛误差,||·||为向量的模;k为迭代次数。
步骤102具体实施方法如下:
步骤201:采用最小二乘法对抛物面天线进行圆弧拟合;
步骤202:根据天线拓扑关系,在步骤201得到的拟合球面上确定天线桁架下表面七个节点坐标P_down;
步骤203:对桁架高度H进行设计,H=0.72m,从而得到桁架上表面相应的七个节点坐标P_up;
步骤204:根据设计要求形面精度δrms<0.5mm,需要对索网索段进行5等分;
步骤205:求解后索网节点和拟合球心的连线与抛物面天线的交点作为前索网节点;
步骤104具体实施方法如下:
步骤401:根据步骤103计算得到的力密度系数和步骤102设计得到的索网节点坐标,采用力密度法建立索网节点线性平衡方程组,通过联立求解这些线性平衡方程组即可得到所有索网连接点的坐标。
步骤402:根据步骤301求解得到的自由索网节点坐标和与桁架相连的索网节点坐标可以求得索网索段长度和索段张力值;
步骤403:根据步骤301计算得到的前索网节点xi和yi坐标,采用抛物面天线公式z=(x2+y2)/(4f)对前索网节点z坐标进行修正;
步骤404:将修正前与修正后的前索网节点坐标误差进行计算,取设计误差为1e-12,判断误差是否满足设计精度要求;若不满足设计精度要求,转至步骤401;若满足设计精度要求,执行步骤404,输出索网所有节点坐标和索段张力值,完成对索网的渐近迭代力密度找形。
步骤105具体实施方法如下:
步骤501:根据步骤104找形得到的索网节点和桁架节点坐标,建立天线空间几何模型,对桁架的天线空间几何模型进行网格划分,生成绝对节点坐标单元,对索网进行控制点插入,生成等几何索单元;
步骤502:根据天线空间几何模型,建立天线有限元模型,根据连续介质力学理论计算单元刚度矩阵,并对单元刚度矩阵进行组装,获得有限元模型的非线性刚度矩阵K;
步骤503:根据104找形得到的索段张力对索网施加预应力F,对天线有限元模型施加边界条件;
步骤504:建立系统静力学平衡方程KU=F,求解静力学平衡方程获得单元节点位移增量U,其中K为步骤502获得的系统非线性刚度矩阵,F为步骤502获得的预应力;
步骤505:根据步骤504获得的位移增量,判断节点位移增量是否收敛,取收敛误差为1e-8,若节点位移增量不收敛,重复步骤502到505,根据单元节点位移增量更新单元节点坐标,若节点位移增量收敛,则完成静力平衡迭代。
综上所述,以上仅为本发明的较佳实施例而已,并非用于限定本发明的保护范围。凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (7)
1.一种空间桁架索网天线等几何找形方法,天线包括支撑桁架和索网,索网由索网节点和索段组成,索网节点分为前索网节点和后索网节点,其特征在于,该方法的具体步骤如下:
步骤1:根据设定的天线设计参数,基于天线拓扑关系,获得支撑桁架和索网当前节点坐标;
步骤301:根据索段当前力密度系数和索网节点当前坐标,采用力密度法建立索网节点线性平衡方程组,并对该方程组求解,获得所有索网节点的坐标;
步骤302:根据步骤301获得的索网节点当前坐标,计算索网索段当前长度和索段当前张力值基于步骤301获得的前索网节点坐标,采用抛物面天线公式znew=(x2+y2)/(4f),对前索网节点的z坐标进行修正;其中,(x,y,z)为修正前的前索网节点坐标,(x,y,znew)为修正后的前索网节点坐标,f为抛物面天线焦距;
步骤303:将修正前与修正后的前索网节点坐标误差进行计算,判断误差是否满足设计精度要求;若不满足设计精度要求,将修正后的(x,y,znew)作为前索网节点当前坐标,重复步骤301-303;若满足设计精度要求,根据修正后的索网节点坐标重新计算索段张力值,以索网节点当前坐标和索段当前张力值作为索网的渐近迭代力密度找形结果。
2.如权利要求1所述一种空间桁架索网天线等几何找形方法,其特征在于,进一步包括:
步骤4:根据步骤303输出的索网节点当前坐标和索段当前张力值,更新索段当前长度和索段当前力密度系数以索段当前张力值索段当前长度和索段当前力密度系数作为迭代的初始值,依据步骤3的找形结果,建立桁架和索网非线性有限元模型,并对桁架和索网有限元模型进行静力平衡迭代,k的初始值取0;
3.如权利要求1所述一种空间桁架索网天线等几何找形方法,其特征在于,步骤1的具体处理过程如下:
步骤101:根据设定的天线设计参数,采用最小二乘法对抛物面天线反射面进行圆弧拟合,获得拟合球的球面方程;
步骤102:基于天线拓扑关系,根据步骤101获得的拟合球的球面方程确定天线桁架下表面N个节点坐标P_down,该N个节点分别为桁架中心点和周边S个端点;
步骤103:设定桁架高度H,基于步骤102确定的天线桁架下表面N个节点坐标P_down,获得桁架上表面对应的N个节点坐标P_up,计算公式如下:
式(4)中,r为桁架下表面节点和拟合球心之间连线的矢量;
步骤104:基于步骤102确定的天线桁架下表面N个节点坐标P_down,将后索网设计成环向索网格式,前索网设计成三向索网格式,根据三向索网天线反射面形面精度估算公式,确定估算精度小于或等于设定精度的索网分段数M,根据桁架上表面节点坐标对后索网进行等分,获得后索网节点坐标P_rear;
步骤105:根据步骤104获得的后索网节点坐标P_rear,基于后索网节点和拟合球心之间的连线与抛物面天线的交点,建立如下非线性方程组:
通过对式(5)进行求解,获得所述交点坐标(xP_fornt,yP_fornt,zP_fornt),将所述交点作为前索网节点P_front,其中,(xo,yo,zo)为步骤101获得的拟合球面的球心坐标,(xP_rear,yP_rear,zP_rear)为步骤104获得的后索网节点坐标;最终获得支撑桁架当前节点坐标P_down、P_up和索网当前节点坐标P_rear、P_front。
4.如权利要求2所述一种空间桁架索网天线等几何找形方法,其特征在于,步骤4具体步骤如下:
步骤401:根据步骤3找形得到的索网当前节点坐标和步骤1获得的支撑桁架当前节点坐标,建立天线空间几何模型,对桁架的天线空间几何模型进行网格划分,生成绝对节点坐标单元,对索网进行控制点插入,生成等几何索单元;
步骤402:根据天线空间几何模型,建立天线有限元模型,根据连续介质力学理论计算单元刚度矩阵,并对单元刚度矩阵进行组装,获得有限元模型的非线性刚度矩阵K;
步骤403:根据步骤3找形得到的索段张力对索网施加预应力F,对天线有限元模型施加边界条件;
步骤404:建立系统静力学平衡方程KΔU=F,求解静力学平衡方程获得单元节点位移增量ΔU,其中K为步骤402获得有限元模型的非线性刚度矩阵,F 为步骤403获得的预应力;
步骤405:根据步骤404获得的节点位移增量ΔU,判断节点位移增量ΔU是否收敛,若节点位移增量不收敛,根据单元节点位移增量更新索网当前节点坐标和支撑桁架当前节点坐标,重复步骤402到405,若节点位移增量收敛,则完成静力平衡迭代。
5.如权利要求1所述一种空间桁架索网天线等几何找形方法,其特征在于,索网天线为六棱柱模块索网天线或三棱柱模块索网天线。
6.如权利要求3所述一种空间桁架索网天线等几何找形方法,其特征在于,N=7,S=6。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710138544.8A CN106934147B (zh) | 2017-03-09 | 2017-03-09 | 一种空间桁架索网天线等几何找形方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710138544.8A CN106934147B (zh) | 2017-03-09 | 2017-03-09 | 一种空间桁架索网天线等几何找形方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106934147A CN106934147A (zh) | 2017-07-07 |
CN106934147B true CN106934147B (zh) | 2020-03-31 |
Family
ID=59433840
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710138544.8A Active CN106934147B (zh) | 2017-03-09 | 2017-03-09 | 一种空间桁架索网天线等几何找形方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106934147B (zh) |
Families Citing this family (17)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107491594B (zh) * | 2017-07-24 | 2020-07-14 | 北京空间飞行器总体设计部 | 一种索网天线形面精度计算方法 |
CN108021758B (zh) * | 2017-12-07 | 2021-02-26 | 上海宇航系统工程研究所 | 基于有限元的索网天线等张力找形方法 |
CN108168504B (zh) * | 2017-12-26 | 2020-05-22 | 中国人民解放军战略支援部队信息工程大学 | 一种基于拟合球心法的天线旋转中心标校方法 |
CN108388698A (zh) * | 2018-01-26 | 2018-08-10 | 东南大学 | 一种改进的力密度平衡形态分析方法 |
CN110096719A (zh) * | 2018-01-31 | 2019-08-06 | 上海二十冶建设有限公司 | 一种大面积焊接球网架结构的变形预控制方法 |
CN109101747B (zh) * | 2018-08-29 | 2023-04-07 | 西安空间无线电技术研究所 | 一种环形网状反射器索网结构的确定方法及装置 |
CN109244678B (zh) * | 2018-09-27 | 2023-10-17 | 中国科学院国家天文台 | 一种索系基础望远镜反射面结构 |
CN110909435B (zh) * | 2019-11-28 | 2023-10-10 | 上海宇航系统工程研究所 | 一种网状天线索网形态分析方法 |
CN111241731B (zh) * | 2020-01-07 | 2023-08-18 | 国网浙江省电力有限公司电力科学研究院 | 一种依据水平张力的输电导线找形方法 |
CN111814231B (zh) * | 2020-06-23 | 2022-06-24 | 同济大学 | 一种基于迭代找形的网壳形状优化方法 |
CN112699462B (zh) * | 2020-12-29 | 2024-04-09 | 中国航空工业集团公司西安飞机设计研究所 | 一种大跨度开剖面舱体空间桁架结构设计方法 |
CN112800521A (zh) * | 2021-01-28 | 2021-05-14 | 北京市建筑设计研究院有限公司 | 一种适用于索网结构找形的低索夹滑移力找形方法 |
CN112800520A (zh) * | 2021-01-28 | 2021-05-14 | 北京市建筑设计研究院有限公司 | 一种索网结构的找形方法 |
CN113779828B (zh) * | 2021-08-30 | 2023-10-03 | 航宇救生装备有限公司 | 降落伞张满形状计算方法 |
CN114497999B (zh) * | 2022-01-14 | 2023-02-14 | 西安电子科技大学 | 一种面向热稳定性的多材料复合索网设计与制作方法 |
CN114707116B (zh) * | 2022-03-24 | 2024-06-07 | 西安电子科技大学 | 一种基于代理模型的索网天线制造误差灵敏度分析方法 |
CN115270540B (zh) * | 2022-05-26 | 2024-04-12 | 青岛科技大学 | 一种索膜桁架天线优化设计方法 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2006136843A1 (en) * | 2005-06-23 | 2006-12-28 | Bae Systems Plc | Improvements in or relating to antennas |
CN103399981A (zh) * | 2013-07-03 | 2013-11-20 | 上海同磊土木工程技术有限公司 | 一种索杆梁混合体系找形的力密度有限元混合法 |
CN103761368A (zh) * | 2014-01-02 | 2014-04-30 | 西安电子科技大学 | 一种索网反射面天线形面精度与张力状态的同时设计方法 |
CN104504284A (zh) * | 2015-01-04 | 2015-04-08 | 西安电子科技大学 | 一种基于悬链线单元的松弛索网找形方法 |
CN104795640A (zh) * | 2015-03-30 | 2015-07-22 | 西安电子科技大学 | 一种预防周边桁架式可展开天线索网热致松弛的方法 |
CN105740517A (zh) * | 2016-01-25 | 2016-07-06 | 西安电子科技大学 | 一种考虑在轨热环境的星载网状可展开天线结构优化方法 |
CN106354906A (zh) * | 2016-08-19 | 2017-01-25 | 西安电子科技大学 | 一种基于三角形面片误差的网状可展开天线找形方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8018079B2 (en) * | 2009-02-23 | 2011-09-13 | Tetraheed Llc | Reciprocating system with buoyant aircraft, spinnaker sail, and heavy cars for generating electric power |
-
2017
- 2017-03-09 CN CN201710138544.8A patent/CN106934147B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2006136843A1 (en) * | 2005-06-23 | 2006-12-28 | Bae Systems Plc | Improvements in or relating to antennas |
CN103399981A (zh) * | 2013-07-03 | 2013-11-20 | 上海同磊土木工程技术有限公司 | 一种索杆梁混合体系找形的力密度有限元混合法 |
CN103761368A (zh) * | 2014-01-02 | 2014-04-30 | 西安电子科技大学 | 一种索网反射面天线形面精度与张力状态的同时设计方法 |
CN104504284A (zh) * | 2015-01-04 | 2015-04-08 | 西安电子科技大学 | 一种基于悬链线单元的松弛索网找形方法 |
CN104795640A (zh) * | 2015-03-30 | 2015-07-22 | 西安电子科技大学 | 一种预防周边桁架式可展开天线索网热致松弛的方法 |
CN105740517A (zh) * | 2016-01-25 | 2016-07-06 | 西安电子科技大学 | 一种考虑在轨热环境的星载网状可展开天线结构优化方法 |
CN106354906A (zh) * | 2016-08-19 | 2017-01-25 | 西安电子科技大学 | 一种基于三角形面片误差的网状可展开天线找形方法 |
Non-Patent Citations (2)
Title |
---|
可展天线的柔性索网结构找形分析方法;李团结等;《宇航学报》;20080530;全文 * |
索网式天线结构预拉力优化的新方法;杨东武等;《西安电子科技大学学报》;20080420;全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN106934147A (zh) | 2017-07-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106934147B (zh) | 一种空间桁架索网天线等几何找形方法 | |
CN109766656B (zh) | 一种基于拓扑优化的梯度点阵结构设计方法 | |
CN109145427A (zh) | 一种基于三周期极小曲面的多孔结构设计与优化方法 | |
CN110909435B (zh) | 一种网状天线索网形态分析方法 | |
CN108629140B (zh) | 一种基于测地线距离的带孔复合材料结构设计优化方法 | |
Nie et al. | Form finding and design optimization of cable network structures with flexible frames | |
CN107590325A (zh) | 一种基于Shepard插值的纤维增强复合材料结构优化方法 | |
CN110955941B (zh) | 基于向量场的复合材料结构优化设计方法及设备 | |
CN109902350B (zh) | 对变截面梁的截面惯性矩进行模型修正中克服模态交换的方法 | |
JP2020504859A (ja) | しわを抑制する空間薄膜構造治具形状を最適化する設計方法 | |
CN113158510B (zh) | 一种索网结构的优化找形方法 | |
CN111488656A (zh) | 基于apdl和遗传算法的天线结构温度保型优化设计方法 | |
CN110096798B (zh) | 一种多状态有限元模型修正的方法 | |
CN103065015B (zh) | 一种基于内力路径几何形态的承载结构低碳节材设计方法 | |
CN110704950A (zh) | 一种消除自由飞配平载荷下飞机变形中刚性位移的方法 | |
CN108875132B (zh) | 一种多工况下应力约束类桁架连续体结构的拓扑优化方法 | |
CN104834795B (zh) | 包带连接结构接触摩擦非线性特性模拟方法及系统 | |
US11514210B2 (en) | Component-wise reduced-order model design optimization such as for lattice design optimization | |
CN110147571B (zh) | 一种组件结构的拓扑优化方法及装置 | |
CN110083946B (zh) | 一种基于无约束优化模型的多状态模型修正的方法 | |
CN108319800B (zh) | 一种基于功能目标的张力结构通用体型设计方法 | |
CN114528735A (zh) | 索膜反射面可展开天线的结构优化方法和装置 | |
CN113326582B (zh) | 一种基于应力分布的变密度点阵结构及其设计方法 | |
CN113139313A (zh) | 基于有限元理论的等效轻量化结构分析方法 | |
CN110110396B (zh) | 有限元模型修正过程中克服模态交换的方法 |
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 |