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

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

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

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

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

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

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

对该膜边界 Γ 进行小位移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) 。

在无限大薄膜上选一条与实际薄膜有相同位置、 形状虚构边界( 图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) 获得。至此,可求解任意几何边界形状薄膜固有频率。

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

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

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

或写成

亦可写为

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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