边界几何参数对空间平面张拉膜结构固有频率影响研究-项目案例-污水池加盖-反吊膜|膜加盖-除臭加盖-膜结构公司-上海华喜膜结构工程有限公司
网站首页 解决方案 项目案例 新闻动态 膜材介绍 关于华喜 联系方式 EN
首页 > 新闻动态 > 公司动态

边界几何参数对空间平面张拉膜结构固有频率影响研究

发布时间:2019年9月30日 点击数:2939

大型空间平面张拉膜结构因质量轻、包装折叠效率高、易于展开等特点,广泛用于太空可展开结构,如太阳帆、空间望远镜防护罩、可展开反射阵列天线等。 周边索网张拉薄膜结构,见图1、图2。该结构由支撑杆、张拉悬索、平面薄膜组成,其中,张拉悬索包括外悬索、牵连索、内悬索。

针对该类型张拉薄膜结构,Sakamoto[1]理论上证明薄膜外形为圆弧状时膜内应力为最佳张紧状态,能有效避免褶皱产生; 并研究该结构在膜、索给定应变条件下内外悬索截面积、索膜总质量随圆弧数目变化的函数关系。Park等[2]研究该结构的鲁棒性、分析膜面褶皱对固有频率影响认为,在角点支撑发生面内位移情况下外索能有效抑制角点位移向膜面内传播。Leifer等[3]用ANSYS分析薄膜边界剪切变形下褶皱出现临界条件及褶皱幅值、波数。Tessler等[4]研究由薄膜结构中索、膜大变形产生的几何非线性,并用薄壳屈曲理论分析褶皱幅值。汪有伟等[5]研究力与边界弧长数奇偶性质关系,分别给出优缺点。

图1 索网张拉阵面 Fig.1 The photo of the membrane array antenna

图1 索网张拉阵面 Fig.1 The photo of the membrane array antenna   下载原图

图2 周边索网张拉薄 Fig.2 Web-cable girded membrane structure

图2 周边索网张拉薄 Fig.2 Web-cable girded membrane structure   下载原图

以上对周边索网张拉薄膜结构研究多集中于几何参数优化及薄膜在轨工作表面精度( 褶皱) ,而有关结构参数对固有频率影响研究较少,薄膜固有频率为反映结构鲁棒性及保证形面精度的重要指标。

Kang等[6,7,8]提出归一化动态响应函数法,建立任意凸边界形状薄膜固有频率求解模型及多域法模型, 可求解凹陷边界或区域内部不连续薄膜固有频率并推广到任意边界形状板固有频率求解。然而,该方法假设薄膜任意点受力相同,均为二轴等值拉伸应力状态, 理论计算不存在任何问题,而工程应用中并非任何边界形状膜可实现此应力状态。

为此,本文提出分段圆弧状边界膜,通过边界索管道式张拉,理论上可实现均匀应力状态; 并针对周边索网张拉薄膜结构,研究索膜在一定应力水平下结构参数对薄膜固有频率影响。

1任意边界形状薄膜振动模型

1. 1凸单连通区域薄膜振动模型

1. 1. 1归一化动态响应函数

研究任意固定边界形状薄膜固有频率,设薄膜: 1处于均匀二轴等值拉伸应力状态; 2不承受任何弯矩、剪切; 3材料均匀各向同性,且面内变形较小,满足Hooker定律; 4面内张拉刚度远大于横向几何刚度,忽略面内振动。

考虑任意形状均匀薄膜,单位面积质量为 ρ,单位长度张力为T,变形前完全处于xoy平面,见图3虚线。

图3任意边界形状薄膜Fig. 3 Arbitrarily shaped membrane

图3任意边界形状薄膜Fig. 3 Arbitrarily shaped membrane   下载原图

对该膜边界 Γ 进行小位移u( rΓ,t) 谐波激励,其中rΓ为边界任意点位置矢量。用w( r,t) 表示薄膜在t时刻横向位移,其中r为薄膜任意点位置矢量,则其横向振动偏微分方程及边界条件[9]分别为

 
 

对谐波激励,用分离变量法将w( r,t) 、u( rΓ,t) 写为

 
 

式中: ω 为激励源角频率。

代入式( 1) 、( 2) 得

 
 

式中: 为系统特征值( 波数) 。

式( 5) 为Helmholtz方程,坐标变换后可得Bessel方程,基本解为Bessel函数,α 阶Bessel函数表达式[10]

 

引入归一化动态响应函数为

 

式中: ω 为激励源角频率; r,rk分别为点P,Pk位置矢量; J0( x) 为零阶第一类Bessel函数。

该函数描述薄膜边界点Pk做单位振幅、频率为 ω 的简谐振动时点P的位移响应。( Λ r - rk) 量纲为1, 称F( r,ω) 为归一化动态响应函数。

式( 6) 为连续边界条件,若将薄膜边界离散化,即用诸多离散点Pi( i = 1,2,…,N) 表示边界,则可得对应式( 6) 的离散边界条件为

 

当离散点个数N趋于无穷时,式( 8) 则会收敛为式( 6) 。

1. 1. 2任意固定边界薄膜方程

在无限大薄膜上选一条与实际薄膜有相同位置、 形状虚构边界( 图3) ,沿边界分配N个点。若P1, P2,. . . ,PN做简谐振动,设振幅为A1,A2,…,AN,通过线性叠加,点P总位移响应计算式为

 

利用式( 8) 可得

 

写成矩阵形式为

 

式中: SM( Λ) 为N × N对称系统矩阵,其中每项表达式为SMik= J0( Λ| ri- rk|) ; A为向量,表示每个边界点在模态坐标中的比值。

等式( 11) 称为系统矩阵方程。对于周边固定薄膜的自由振动分析,式( 11) 中令U = 0,得

 

求解det[SM( Λ) ]= 0可得系统特征值即结构固有频率,薄膜相应模态由式( 9) 获得。至此,可求解任意几何边界形状薄膜固有频率。

1. 2复杂边界形状薄膜振动模型

用单域法求解具有凹陷或含孔洞薄膜固有频率时计算结果往往不收敛,此因归一化动态响应函数为波形函数,考虑无限区域时会从一点沿全方位传播,而薄膜凹陷或孔洞区域会隔断波直线传播,因此单区域法求解得不到正确结果。

引入归一化动态响应函数多域法求解具有凹陷或含孔洞薄膜固有频率。该法将薄膜整个区域在缺陷处分成若干子区域,对每个子区域用单域法建立系统矩阵方程组,各子系统矩阵方程组通过公共边界条件相互耦合建立整体系统矩阵,从而获得系统特征值及相应模态。图4为凹陷薄膜被一条虚线公共边界分割成两子区域D、D,分别由边界 围成, 处于相同位置 。

图4凹陷薄膜被分成两子区域Fig. 4 Concavely shaped membrane subdivided into two domains

图4凹陷薄膜被分成两子区域Fig. 4 Concavely shaped membrane subdivided into two domains   下载原图

对两子区域分别用归一化动态响应函数单区域法获得系统矩阵方程分别为

 
 

或写成

 
 

亦可写为

 
 
 
 

式中: 向量U1,U2为薄膜边界 Γ12离散点位移; 向量为公共边界 点未知位移; 向量A1,A2为边界 Γ12上每个离散点在模态坐标中的比值; 向量为公共边界 上每个离散点在模态坐标中的比值。

对固定边界条件, U1= U2= 0 ,由式( 16 ) 、 ( 18 ) 知,向量A1,A2可分别由 表示,即

 
 

分别代入式( 17) 、( 19) ,且在边界 Γa上满足

 
 

式中: n为公共边界法线方向,该凹陷薄膜在公共边界法线方向上C1连续。

将式( 17) 、( 19) 代入式( 22) 、( 23) 得

 

系统矩阵SM( Λ) 、向量A表达式为

 
 

式中: SM'( )为∂SM( )/ ∂n。系统特征值可令矩阵行列式等于零获得,即det[SM( Λ) ] = 0。

同理,可导出归一化动态响应函数三、四等多区域法,据薄膜边界形状合理划分区域,求出具有凹陷或含孔洞薄膜固有频率及模态。

2计算实例

2. 1算例1: 圆形薄膜固有频率

圆形边界薄膜固有频率存在理论解,本文分别用归一化动态响应函数单区域法及有限元法求解,并与理论解对比。选半径为单位1的均匀圆形薄膜,见图5,在边界上分别等间距取8,12,16点。

图5 圆形薄膜 Fig.5 The circular membrane

图5 圆形薄膜 Fig.5 The circular membrane   下载原图

归一化动态响应函数法中,系统矩阵行列式值的对数随系统特征值 Λ 变化关系见图6,极值点横坐标为薄膜固有频率。

三种方法结果见表1。由表1看出,用归一化动态响应函数法求解结果非常精确,随边界点数增加收敛于理论解。取8或12个边界点时系统方程阶数较低, 无法获取前八阶全部固有频率。而有限元法将无限自由度结构离散成有限自由度结构,增大结构刚度,故求解的固有频率值大于理论解。随单元划分细化,有限元解收敛于理论解。该算例表明归一化动态响应函数法对固定边界薄膜固有频率求解收敛性、精确性良好。

图6圆形薄膜特征值Fig. 6 The eigenvalues of circular membrane

图6圆形薄膜特征值Fig. 6 The eigenvalues of circular membrane   下载原图

表1圆形薄膜固有频率Tab. 1 The natural frequencies of circular membrane with fixed boundary     下载原表

表1圆形薄膜固有频率Tab. 1 The natural frequencies of circular membrane with fixed boundary

2. 2算例2: L - 形薄膜固有频率

图7 L-形薄膜( 单位: m) Fig. 7 The L-shape membrane

图7 L-形薄膜( 单位: m) Fig. 7 The L-shape membrane   下载原图

L - 形薄膜固有频率不存在理论解,本文分别用归一化动态响应函数多区域法及有限元法求解, 将两结果对比,表明多区域法的有效性。将L - 形薄膜分成两子区域,在公共边界取3个离散点,固定边界取24个离散点,见图7。

归一化动态响应函数两区域法中,系统矩阵行列式值的对数随系统特征值 Λ 变化关系,如图8,极小值点的横坐标表示薄膜的固有频率。

图8 L-形薄膜特征值Fig. 8 The eigenvalues of L-shape membrane

图8 L-形薄膜特征值Fig. 8 The eigenvalues of L-shape membrane   下载原图

两区域法与有限元法结果见表2。由表2看出归一化动态响应函数多区域法的有效性与准确性。与有限元法相比,多区域法计算所需节点数目较少,且计算中无需数值积分,计算量大大减少。

表2 L - 形薄膜的固有频率Tab. 2 The natural frequencies of L-shape membrane with fixed boundary     下载原表

表2 L - 形薄膜的固有频率Tab. 2 The natural frequencies of L-shape membrane with fixed boundary

2. 3算例3: 圆弧状薄膜固有频率

周边索网张拉薄膜结构中薄膜部分见图9,各边界等份若干段圆弧。其中l1= 20 m,l2= 4 m为薄膜长、 高; N1,N2为薄膜长、高边界圆弧段数; θ12为长、高上圆弧圆心角。

图9圆弧状固定边界薄膜Fig. 9 Shapes of circular segments of a membranewith fixed boundary

图9圆弧状固定边界薄膜Fig. 9 Shapes of circular segments of a membranewith fixed boundary   下载原图

薄膜受均匀二轴等值拉伸应力T = 1 440 Pa,材料参数见表3。

表3材料参数Tab. 3 Material parameters in the analysis     下载原表

表3材料参数Tab. 3 Material parameters in the analysis

保持条件 θ1= θ2= 10°,N2= 4不变前提下改变长边分段数,用多域法所得薄膜固有频率见表4。由表4看出,随长边圆弧段数增加薄膜面积变大,薄膜固有频率降低。同理,保持条件N1= 20,N2= 4,θ2= 10° 不变前提下改变长边圆弧圆心角,用多域法所得薄膜固有频率见表5。由表5看出,随圆心角增大薄膜面积减小,薄膜固有频率增大。

表4薄膜固有频率随长边弧长数目变化关系Tab. 4 The relationship between the natural frequencies of membrane of fixed boundary and number of circular segments     下载原表

表4薄膜固有频率随长边弧长数目变化关系Tab. 4 The relationship between the natural frequencies of membrane of fixed boundary and number of circular segments

表5薄膜固有频率随长边弧长圆心角变化关系Tab. 5 The relationship between the natural frequencies of membrane of fixed boundary and angel of circular segments     下载原表

表5薄膜固有频率随长边弧长圆心角变化关系Tab. 5 The relationship between the natural frequencies of membrane of fixed boundary and angel of circular segments

3结论

( 1) 利用归一化动态响应函数法求解任意形状固定边界薄膜固有频率,给出高效、精确求解薄膜固有频率方法,计算所需节点数远少于有限元方法。

( 2) 对复杂薄膜边界,所提多域法可有效弥补单域法不足。对周边索网张拉结构薄膜给出固有频率及边界结构参数( 曲率半径、弦长) 变化关系。

专题报道             more...
  • 轨道交通中膜结构的应
    ...

    查看更多

  • 膜结构建筑保温内衬技
    刚查县为青海省海北藏族自治州辖县,青海省措温波高原海滨藏城演艺中心,作为刚查县的标志性建筑,演艺中心为直径50米的圆形建...

    查看更多

  • 膜结构幕墙的应用
    膜结构幕墙是膜结构在建筑外围护结构的应用,具有膜结构的共同特性和优点:膜结构是一种非传统的全新结构方式。...

    查看更多

  • 膜结构屋面的应用
    屋盖是房屋最上部的围护结构,应满足相应的使用功能的要求,为建筑提供适宜的内部空间环境。屋盖也是房屋顶部的承重结构,受到材...

    查看更多

  • 膜结构应用于环保工程
    随着我国国民经济飞速发展和市政基础设施建设全面展开,特别是污水处理厂等环保项目日益增多,其中有相当数量的污水处理厂的厌氧...

    查看更多

  • 膜结构在污水处理厂中
    相当数量的污水处理厂的厌氧池、污泥浓缩池、生物絮凝池等建于居民区、厂区的周边,污水池的环境、风貌及污水臭味等直接影响人们...

    查看更多

关于华喜

硬件实力 质量控制 发展历程 公司简介

软件实力 经营理念  解决方案 联系方式

中国华喜建筑网站

+021-59198545 400-176-6885 dshx@hxmjg99.com www.hxmjg.com 沪ICP备08009856号 使用条款