CN102663229A - 一种索杆张力结构的找形方法 - Google Patents
一种索杆张力结构的找形方法 Download PDFInfo
- Publication number
- CN102663229A CN102663229A CN2012100600976A CN201210060097A CN102663229A CN 102663229 A CN102663229 A CN 102663229A CN 2012100600976 A CN2012100600976 A CN 2012100600976A CN 201210060097 A CN201210060097 A CN 201210060097A CN 102663229 A CN102663229 A CN 102663229A
- Authority
- CN
- China
- Prior art keywords
- length
- node
- cable
- partiald
- weight coefficient
- 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.)
- Granted
Links
Images
Abstract
本发明公开了一种索杆张力结构的找形方法,属于空间结构的建筑设计和结构设计领域。本发明方法步骤是:设定索杆张力结构的拓扑关系,并给出其约束条件和压杆长度;确定索杆张力结构找形分析中的目标函数;设定每根拉索和压杆的长度权重系数以及每个边界条件的权重系数,并给出梯度法优化中的初始步长;计算每个节点坐标变量在xk处的梯度向量。本发明从能量的角度出发,并将索杆张力结构的边界条件等约束写入梯度法优化方法的功能函数中,从而较容易地解决了具有约束条件的张力结构的找形问题。传统找形方法中加入了过多的约束条件以及几何限制,本发明初始限制条件较少,可以找到大量非规则体系。
Description
技术领域
本发明公开了一种基于梯度法优化的索杆张力结构的找形方法,属于空间结构的建筑设计和结构设计领域。
背景技术
索杆张力结构的刚度是由预应力提供的,它们需要一个过程达到自平衡状态,这个状态一般称为预应力态。而结构的预应力态取决于结构的形状,所以达到预应力态的过程也称之为找形。而索杆张力结构的找形同时也是一个找力的过程,在具体的找形方法中可以以形状参数为变量,也可以以内力参数为变量。因此,可以将找形方法划为“找力”和“找形”两大类,前者主要指搜索可行预应力或预应力优化,后者的代表性方法有力密度法、非线性有限元法和动力松弛法。力密度法最早由Schek提出,用于索网等只有受拉构件结构的找形分析中。随后,众多学者对力密度法进行了改进,使之能够运用于含有拉压杆的张拉整体结构的找形。但在运用力密度法对索杆张力结构找形时,确定力密度后,结构的几何构形并不唯一,而这一问题一直没有很好的解决。
在一个三维索杆结构中,假定自由节点i与节点j、节点k相连,如图1所示。索杆体系的节点平衡方程可以表述为:
对于任意一个连接节点i、j的单元(i,j),其内力为fij,长度为lij;作用于节点i的外荷载为pi,在三个方向的分量分别为:pix,piy,piz。力密度为单元内力和长度的比值,即为qij=fij/lij。式1可整理为:
假定索杆张力结构有b个单元,ns个节点,其中n个自由节点,nf个固定节点,引入b×n的拓扑矩阵:
在Cs中将自由节点排列在固定节点之前,Cs又可以拆分为自由节点拓扑矩阵和固定节点拓扑矩阵:Cs=[C Cf],则自由节点的力平衡关系可表示为(以x方向为例):
CTQCx+CTQCfxf=Px 式4
式中:x为自由节点坐标向量;xf为固定节点坐标向量;Px为节点外荷载x方向向量;C为自由节点拓扑矩阵;Cf为固定节点拓扑矩阵;Q为力密度对角矩阵。
在初始预应力确定的过程中,不考虑外荷载的作用,并且如果仅考虑结构是自平衡的(没有固定点),则结构的平衡可以简化为:
其中D为力密度矩阵,可以表示为:
D=CTQC 式6
图2所示为平面张拉整体结构,由4根拉索和2根压杆组成,压杆通过拉索相连,其本身互不相连。对于这种自平衡结构,需满足平衡条件式5。如果式5中力密度矩阵D是满秩矩阵,则节点坐标需满足:
x=y=z=0 式7
这样图2所示的结构就缩为1个点了。如果力密度矩阵D是非满秩矩阵,则可以通过奇异值分解等数学手段得到其零向量空间。假定4根拉索和2根压杆的力密度比值为:1∶1∶1∶1∶-1∶-1,则可以获得多组节点坐标解。如图3和图4所示的两种几何构形都满足式5。该结构的力密度矩阵D和节点坐标可以表示为:
式中a,b,c为系数,可取任意值。
如果4根拉索和2根压杆的力密度比值不是:1∶1∶1∶1∶-1∶-1,则节点坐标无解。结构的力密度矩阵D和节点坐标可以表示为:
结构的所有节点都集中于一点上。
发明内容
为了克服力密度法找形过程中几何构形不唯一的缺点,本发明从能量的角度出发,提供了一种基于优化方法的索杆结构找形方法。另外,现有的找形方法大多针对自平衡的索杆体系,而土木工程中应用的结构大多有约束条件(支座等边界条件)。本发明将这些约束条件写入梯度法优化方法的功能函数中,从而较容易地解决了该问题。
一种索杆张力结构的找形方法,其特征在于,包括以下步骤:
第一步、设定索杆张力结构的拓扑关系,并给出其约束条件和压杆长度;
第二步、确定索杆张力结构找形分析中的目标函数,
其中,x为所有节点坐标向量,表示为:{x1,x2,x3,…,xp,…,xd*n}T;n为结构的节点数;d为结构的维数,当体系为二维结构时,d=2;当体系为三位结构时,d=3;lj(x)为第j根拉索的计算长度,是坐标向量x的函数;Lm为第m根压杆的计算长度,是坐标向量x的函数;lm为第m根压杆的给定长度;Xi是步骤1给出第i个约束条件的节点坐标值;b为索杆张力结构的杆件数;ns是索杆张力结构的压杆数量;nf是约束条件的数量;wj为第j根拉索的长度权重系数;w_strm为第m根压杆的长度权重系数;w_cori为第i个边界条件的权重系数;
第三步、设定每根拉索和压杆的长度权重系数以及每个边界条件的权重系数,并给出梯度法优化中的初始步长step;
第四步、输入所有节点坐标向量的初始值x0,并带入目标函数,计算得到目标函数为f0;
第五步、在第k步循环过程中,其对应的节点坐标向量为xk,其目标函数值为fk,计算每个节点坐标变量在xk处的梯度向量为:
判断梯度向量中第p个元素是否大于零,如果大于零,xp k+1=xp k-step;否则,xp k+1=xp k+step。从而可以得到下一个循环步的节点坐标向量xk+1;
第六步、计算xk+1对应的目标函数值fk+1,如果fk+1>fk,则step=step/4;否则,step不变,重复第五步,直至完成所有循环步,输出最终节点坐标向量x。
本发明的有益效果是,本发明提出的索杆体系“找形”的方法适合于绝大多数索杆体系,无论是简单的平面结构,还是复杂的空间结构。在给定的初始条件下,该发明可以找到同时满足能量最小和一定形态要求的索杆张力体系。该发明的另外一个优点是可以在目标函数中直接引入边界约束条件,从而可以很方便地求解考虑约束条件的索杆张拉体系的几何形状。传统找形方法中加入了过多的约束条件以及几何限制,因此适合规则索杆体系的找形,本发明初始限制条件较少,可以找到大量非规则体系。
附图说明
下面结合附图及实施例对本发明进一步详细说明:
图1索杆张力结构节点平衡示意图。
图2两杆四索平面张力结构示意图。
图3两杆四索平面张力结构构形一。
图4两杆四索平面张力结构构形二。
图5是本发明计算流程图。
图6是平面索杆张力结构示意图。
图7是平面索杆张力结构计算构形。
图8是四压杆索穹顶结构示意图。
图9是四压杆索穹顶结构计算构形。
图10是五压杆索穹顶结构示意图。
图11是五压杆索穹顶结构计算构形。
具体实施方式
下面结合具体实施案例对本发明进行更为详细的表述:
组成索杆张力结构的拉索和压杆,都只受到轴力的作用,其应变能V可以表述为:
其中F为杆件轴力,l为杆件长度,EA为杆件的轴向刚度。式(1)还可以表述为:
式中q为杆件的力密度,w为权重系数。
和力密度法相等价的能量函数П(x)可以表示为:
其中,qj为第j根杆件的力密度,lj(x)为第j根杆件的计算长度,向量x为节点的x,y,z向的坐标,可以表述为:{x1,x2,x3,…,xp,…,xd*n}T。
考虑式(3)的驻值条件为:
式(4)和节点的平衡方程等价,所以式(3)为力密度法所对应的能量函数。
假定压杆的长度为定值,则可以利用Lagrange乘子法将能量函数П(x,λ)修改为:
其中,λm为第m根压杆所对应的Lagrange乘子;Lm(x)为第m根压杆的计算长度,是坐标向量x的函数;lm为第m根压杆的给定长度。对于图2所示的索杆结构,如果我们将所有拉索的力密度比值定为:1∶1∶1∶1,两根压杆的长度也相等,则图3和图4所示的两种几何构形都符合式(5)所对应的驻值条件。如果将式(5)的功能函数П(x,λ)假定为其它形式:
其中wj为第j根拉索的长度权重系数;w_strm为第m根压杆的长度权重系数。式(6)给出的函数没有实际的物理意义。这个时候将图2所示结构的拉索的力密度比值定为:1∶1∶1∶1,两根压杆的长度也相等,根据式(3)的驻值条件,可以得到图3所示的几何构形是满足条件的唯一解。这主要是因为,式(3)给出的功能函数和长度的平方相关,而根据勾股定理可知,在直角三角形中,斜边长度的平方永远等于两直角边长度的平方和。所以满足式(3)驻值条件的几何构形有无穷多组。而式(6)中没有考虑索杆张力结构的边界条件,确定考虑边界条件约束的索杆张力结构初始平衡态包含如下步骤:
(1)设定索杆张力结构的拓扑关系,并给出其约束条件和压杆长度;
(2)确定索杆张力结构找形分析中的目标函数,
其中,x为所有节点坐标向量,表示为:{x1,x2,x3,…,xp,…,xd*n}T,n为结构的节点数,d为结构的维数(当体系为二维结构时,d=2;当体系为三位结构时,d=3);lj(x)为第j根拉索的计算长度,是坐标向量x的函数;Lm(x)为第m根压杆的计算长度,是坐标向量x的函数;lm为第m根压杆的给定长度;Xi是步骤1给出第i个约束条件的节点坐标值;b为索杆张力结构的杆件数;ns是索杆张力结构的压杆数量;nf是约束条件的数量;wj为第j根拉索的长度权重系数;w_strm为第m根压杆的长度权重系数;w_cori为第i个边界条件的权重系数。
(3)设定每根拉索和压杆的长度权重系数以及每个边界条件的权重系数,并给出梯度法优化中的初始步长step。
(4)输入所有节点坐标向量的初始值x0,并带入目标函数,计算得到目标函数为f0。
(5)在第k步循环过程中,其对应的节点坐标向量为xk,其目标函数值为fk。计算每个节点坐标变量在xk处的梯度向量为:
判断梯度向量中第p个元素是否大于零,如果大于零,xp k+1=xp k-step;否则,xp k+1=xp k+step。从而可以得到下一个循环步的节点坐标向量xk+1。
(6)计算xk+1对应的目标函数值fk+1。如果fk+1>fk,则step=step/4;否则,step不变,重复第5步,直至完成所有循环步,输出最终节点坐标向量x。
利用图5所示的流程编制相应程序,结构具体算例对本发明的具体实施方式作出更为详细的说明:
算例1图6所示的二维索杆结构是由2根压杆和6根拉索组成。其已知条件为:压杆(杆件编号分别为7,8)的长度均为4;已知节点5的坐标为(0,0),节点6的坐标为(0,7);压杆的长度权重系数为1000,边界条件的权重系数为5000。将6根拉索的权重系数矩阵w设定为[1 1 1 1 1 1]T,则6个节点的坐标为:x=[2.1692 4.8308 4.8308 2.1692 0.00007.0000]T;y=[-1.9995 -1.9995 1.9993 1.9993 -0.0000 -0.0000]T。其几何形状如图7所示。
算例2四压杆索穹顶结构是由4根压杆和16根拉索组成,如图8所示。其已知条件为:压杆的长度矩阵均为1;固定支座节点1~4坐标分别为(0,0,0),(10,0,0),(10,10,0)(0,10,0);压杆的长度权重系数为1000,边界条件的权重系数为5000。
若将16根拉索的权重系数矩阵w设定为[1 1 1 1 1 1 1 1 2 2 2 2 2 2 2 2]T,则8个自由节点的坐标见表1。其几何形状如图9所示。
表1四压杆索穹顶节点坐标
算例3五压杆索穹顶结构是由5根压杆和20根拉索组成,如图10所示。其已知条件为:压杆的长度矩阵均为2;固定支座节点1~5坐标分别为(-10.0000,0,0),(-3.0902,-9.5106,0),(8.0902,-5.8779,0),(8.0902,5.8779,0),(-3.0902,9.5106,0);压杆的长度权重系数为1000,边界条件的权重系数为5000。
若将16根拉索的权重系数矩阵w设定为[1 1 1 1 1 1 1 1 1 2 2 2 2 2 2 2 2 2]T,则10个自由节点的坐标如表2所示。其几何形状如图11所示。
表2五压杆索穹顶节点坐标
Claims (1)
1.一种索杆张力结构的找形方法,其特征在于,包括以下步骤:
第一步、设定索杆张力结构的拓扑关系,并给出其约束条件和压杆长度;
第二步、确定索杆张力结构找形分析中的目标函数,
其中,x为所有节点坐标向量,表示为:{x1,x2,x3,…,xp,…,xd*n}T;n为结构的节点数;d为结构的维数,当体系为二维结构时,d=2;当体系为三位结构时,d=3;lj(x)为第j根拉索的计算长度,是坐标向量x的函数;Lm为第m根压杆的计算长度,是坐标向量x的函数;lm为第m根压杆的给定长度;Xi是步骤1给出第i个约束条件的节点坐标值;b为索杆张力结构的杆件数;ns是索杆张力结构的压杆数量;nf是约束条件的数量;wj为第j根拉索的长度权重系数;w_strm为第m根压杆的长度权重系数;w_cori为第i个边界条件的权重系数;
第三步、设定每根拉索和压杆的长度权重系数以及每个边界条件的权重系数,并给出梯度法优化中的初始步长step;
第四步、输入所有节点坐标向量的初始值x0,并带入目标函数,计算得到目标函数为f0;
第五步、在第k步循环过程中,其对应的节点坐标向量为xk,其目标函数值为fk,计算每个节点坐标变量在xk处的梯度向量为:
判断梯度向量中第p个元素是否大于零,如果大于零,xp k+1=xp k-step;否则,xp k+1=xp k+step。从而可以得到下一个循环步的节点坐标向量xk+1;
第六步、计算xk+1对应的目标函数值fk+1,如果fk+1>fk,则step=step/4;否则,step不变,重复第五步,直至完成所有循环步,输出最终节点坐标向量x。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210060097.6A CN102663229B (zh) | 2012-03-08 | 2012-03-08 | 一种索杆张力结构的找形方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210060097.6A CN102663229B (zh) | 2012-03-08 | 2012-03-08 | 一种索杆张力结构的找形方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102663229A true CN102663229A (zh) | 2012-09-12 |
CN102663229B CN102663229B (zh) | 2014-11-05 |
Family
ID=46772720
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201210060097.6A Expired - Fee Related CN102663229B (zh) | 2012-03-08 | 2012-03-08 | 一种索杆张力结构的找形方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102663229B (zh) |
Cited By (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103953116A (zh) * | 2014-05-22 | 2014-07-30 | 中国建筑西南设计研究院有限公司 | 大开口索承网格结构的找形方法 |
CN104765925A (zh) * | 2015-04-14 | 2015-07-08 | 温州大学 | 基于误差敏感性分析的索杆预张力结构张拉成形优选方法 |
CN105205266A (zh) * | 2015-09-23 | 2015-12-30 | 东南大学 | 一种基于优化算法的预应力索杆结构初始构形的设计方法 |
CN105426597A (zh) * | 2015-11-09 | 2016-03-23 | 中国电力科学研究院 | 一种架空输电线路大截面导线的找形确定方法 |
CN105956256A (zh) * | 2016-04-28 | 2016-09-21 | 河海大学 | 一种既有混凝土桥梁复杂区域的承载力评估与加固计算方法 |
CN106934155A (zh) * | 2017-03-13 | 2017-07-07 | 北京市建筑设计研究院有限公司 | 一种索桁结构的找形方法 |
CN107016179A (zh) * | 2017-03-29 | 2017-08-04 | 北京市建筑设计研究院有限公司 | 一种张弦结构的找形方法 |
CN108229054A (zh) * | 2018-01-26 | 2018-06-29 | 东南大学 | 一种基于群论的对称张拉整体结构找形方法 |
CN108319800A (zh) * | 2018-03-20 | 2018-07-24 | 中国建筑西南设计研究院有限公司 | 一种基于功能目标的张力结构通用体型设计方法 |
CN108345729A (zh) * | 2018-01-26 | 2018-07-31 | 东南大学 | 一种基于群论的对称索杆结构找形方法 |
CN108388698A (zh) * | 2018-01-26 | 2018-08-10 | 东南大学 | 一种改进的力密度平衡形态分析方法 |
CN110502810A (zh) * | 2019-08-05 | 2019-11-26 | 中国建筑西南设计研究院有限公司 | 一种车辐式索网的直接找形方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20040134141A1 (en) * | 2001-05-10 | 2004-07-15 | Sven Hebbelinck | Flexible tensioned structure and method of calculating such a structure |
CN101582095A (zh) * | 2009-06-19 | 2009-11-18 | 东南大学 | 确定索杆系静力平衡状态的非线性动力有限元法 |
-
2012
- 2012-03-08 CN CN201210060097.6A patent/CN102663229B/zh not_active Expired - Fee Related
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20040134141A1 (en) * | 2001-05-10 | 2004-07-15 | Sven Hebbelinck | Flexible tensioned structure and method of calculating such a structure |
CN101582095A (zh) * | 2009-06-19 | 2009-11-18 | 东南大学 | 确定索杆系静力平衡状态的非线性动力有限元法 |
Non-Patent Citations (2)
Title |
---|
HOANG CHI TRAN: "Geometric and material nonlinear analysis of tensegrity structures", 《ACTA MECHANICA SINICA》 * |
蔡建国: "半刚性结构找形分析的实用方法", 《土木建筑与环境工程》 * |
Cited By (20)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103953116B (zh) * | 2014-05-22 | 2015-12-30 | 中国建筑西南设计研究院有限公司 | 大开口索承网格结构的找形方法 |
CN103953116A (zh) * | 2014-05-22 | 2014-07-30 | 中国建筑西南设计研究院有限公司 | 大开口索承网格结构的找形方法 |
CN104765925A (zh) * | 2015-04-14 | 2015-07-08 | 温州大学 | 基于误差敏感性分析的索杆预张力结构张拉成形优选方法 |
CN105205266B (zh) * | 2015-09-23 | 2018-03-20 | 东南大学 | 一种基于优化算法的预应力索杆结构初始构形的设计方法 |
CN105205266A (zh) * | 2015-09-23 | 2015-12-30 | 东南大学 | 一种基于优化算法的预应力索杆结构初始构形的设计方法 |
CN105426597A (zh) * | 2015-11-09 | 2016-03-23 | 中国电力科学研究院 | 一种架空输电线路大截面导线的找形确定方法 |
CN105426597B (zh) * | 2015-11-09 | 2019-08-16 | 中国电力科学研究院 | 一种架空输电线路大截面导线的找形确定方法 |
CN105956256A (zh) * | 2016-04-28 | 2016-09-21 | 河海大学 | 一种既有混凝土桥梁复杂区域的承载力评估与加固计算方法 |
CN106934155B (zh) * | 2017-03-13 | 2020-04-17 | 北京市建筑设计研究院有限公司 | 一种索桁结构的找形方法 |
CN106934155A (zh) * | 2017-03-13 | 2017-07-07 | 北京市建筑设计研究院有限公司 | 一种索桁结构的找形方法 |
CN107016179A (zh) * | 2017-03-29 | 2017-08-04 | 北京市建筑设计研究院有限公司 | 一种张弦结构的找形方法 |
CN107016179B (zh) * | 2017-03-29 | 2020-06-16 | 北京市建筑设计研究院有限公司 | 一种张弦结构的找形方法 |
CN108229054A (zh) * | 2018-01-26 | 2018-06-29 | 东南大学 | 一种基于群论的对称张拉整体结构找形方法 |
CN108345729A (zh) * | 2018-01-26 | 2018-07-31 | 东南大学 | 一种基于群论的对称索杆结构找形方法 |
CN108388698A (zh) * | 2018-01-26 | 2018-08-10 | 东南大学 | 一种改进的力密度平衡形态分析方法 |
CN108345729B (zh) * | 2018-01-26 | 2021-07-27 | 东南大学 | 一种基于群论的对称索杆结构找形方法 |
CN108319800A (zh) * | 2018-03-20 | 2018-07-24 | 中国建筑西南设计研究院有限公司 | 一种基于功能目标的张力结构通用体型设计方法 |
CN108319800B (zh) * | 2018-03-20 | 2021-10-26 | 中国建筑西南设计研究院有限公司 | 一种基于功能目标的张力结构通用体型设计方法 |
CN110502810A (zh) * | 2019-08-05 | 2019-11-26 | 中国建筑西南设计研究院有限公司 | 一种车辐式索网的直接找形方法 |
CN110502810B (zh) * | 2019-08-05 | 2023-05-30 | 中国建筑西南设计研究院有限公司 | 一种车辐式索网的直接找形方法 |
Also Published As
Publication number | Publication date |
---|---|
CN102663229B (zh) | 2014-11-05 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102663229B (zh) | 一种索杆张力结构的找形方法 | |
Munk et al. | Topology and shape optimization methods using evolutionary algorithms: a review | |
CN104504284B (zh) | 一种基于悬链线单元的松弛索网找形方法 | |
CN106934155B (zh) | 一种索桁结构的找形方法 | |
CN105426594A (zh) | 一种基于时空场和条件插值的平稳均质风场快速模拟方法 | |
CN104933528B (zh) | 一种基于稀疏矩阵技术快速形成电力系统潮流计算中雅可比矩阵的方法 | |
CN105426599A (zh) | 一种简化的拉线塔主柱有限元模型的计算方法 | |
CN105113628B (zh) | 一种预应力索杆结构的对称找力方法 | |
CN102663228B (zh) | 一种索杆张力结构初始预应力分布确定方法 | |
CN106599509A (zh) | 一种模拟弦支穹顶结构施加预应力的方法 | |
CN101261651A (zh) | 主梁断面气动自激力的全过程数值模拟方法 | |
CN113792440B (zh) | 计算浮式结构物在非稳态载荷作用下结构力响应的方法 | |
CN107016179B (zh) | 一种张弦结构的找形方法 | |
CN105888068A (zh) | 一种柔性建筑的建造方法 | |
Lu et al. | Form-finding analysis of irregular tensegrity structures by matrix iteration | |
González et al. | Construction of a unit cell Tensegrity structure | |
CN103399981A (zh) | 一种索杆梁混合体系找形的力密度有限元混合法 | |
CN106021186A (zh) | 一种高效求解大规模非线性随机结构系统状态的多尺度迭代方法 | |
Wang et al. | Regular hexagonal grid single-layer cylindrical latticed shell parametric design | |
CN110704894A (zh) | 斜拉桥桥塔地震响应的计算方法 | |
Qi et al. | Analysis on 500kV Compact and Narrow Basal Transmission Tower in ANSYS Modeling and Dynamic Properties | |
Xue et al. | Research status and prospects of cable domes | |
Yang et al. | Analysis of weak parts of narrow transmission towers and the reinforcement measure | |
Ibragimov | Convertible mesh dome | |
Xie et al. | Stability analysis and study of space grid structure |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C53 | Correction of patent of invention or patent application | ||
CB02 | Change of applicant information |
Address after: 211161 18, Nanjing Road, Jiangning street, Jiangning District, Jiangsu, China Applicant after: Southeast University Address before: Four pailou Nanjing Xuanwu District of Jiangsu Province, No. 2 210096 Applicant before: Southeast University |
|
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20141105 Termination date: 20210308 |
|
CF01 | Termination of patent right due to non-payment of annual fee |