基于有限元分析的林木风致随机振动

基于有限元分析的林木风致随机振动

随机噪动是自然界普遍存在的现象。风一到，树就随风飘动。风一停不下来，房间就会缩小。只有树枝和树叶才能随时随地旋转。枝干上几乎没有线条。它像铺石。风停不下来。这是树木随机噪动的现象。而强风来临时,树木会产生很大的摆幅,甚至倾倒、折断。一直以来强风的造成的森林灾害对社会产生了重大的经济损失,严重影响了当地的国民经济增长,同时对森林生态系统的稳定也造成了很大的影响。近年来风害发生的强度和频率越来越高,引起越来越多的国内外学者注意,国内学者李秀芬、侯倩分析了风害产生的类型及因素,并为森林管理人员提供了防灾减灾措施。林木风害是由于其风致随机振动平衡的破坏引起,为防风害有必要了解林木的风致随机振动,进行风振时程分析是最有效的方法之一。为此在本研究中讨论了产生随机振动的脉动风荷载数值模拟,并建立了林木的有限元模型,通过导入模拟的风荷载得到了其在脉动风荷下的随机动力响应。1自然风作用下的弹性变形风是空气分子的运动,由分子从气压大的地方流向气压小的地方而产生。科学研究者依据得到大量风的实测资料,通过风的顺风向时程曲线,得出自然风是由平均风和脉动风这两部分所组成,平均风周期较长,通常在10分钟以上,远远大于一般的结构自振周期,相当于静力作用在结构上。脉动风周期短,一般只有几秒左右,相当于动力作用,激起结构振动。因此在自然风作用下,林木结构的风荷载可认为是由两部分组成:一是平均风作用下的静风荷载,二是由于自然风的紊流成分诱发导致的随机动荷载。林木在随机振动平衡未被破坏时,主要产生线弹性变形,依据叠加原理,林木的在风载作用下的变形可分为平均风和脉动风分别作用下的叠加。也即:式(1)中:为平均风作用下产生的变形,可用静力学求解。s(t)为脉动风作用下的变形,需要应用动力学求解。林木的随机振动是由脉动风作用产生,本研究主要讨论的就是模拟脉动风荷导致的林木随机振动。1.1平均风速沿高度变化的规律由于风在流动中与地表摩擦,使得风速随距离地表高度的增加而增大。Davenport等根据分析实测得到的结果,提出平均风速沿高度变化的规律如下:式(2)中:表示标准高度为10m处的平均风速(m/s),α为地面粗糙度系数。1.2平均风速谱的确定Davenport通过分析总结世界上不同地点,不同高度处测得的90多次强风记录,于1961年提出风速谱的经验公式:式(3)中:K为地面粗糙度系数,表示标准高度为10m处的平均风速(m/s),n表示脉动风频率,,其风速谱峰值不随高度而变化。除此之外还有其他常用风速谱如Kaimal谱和Simiu谱,它们的风速谱都随高度变化,在本研究中脉动风模拟采用Da-venport谱。1.3功率谱的分解脉动风是一种随机干扰,在模拟之前必须确定其概率分布及功率谱密度,功率谱密度是脉动风最重要的统计特征,能反映出在某一频率上脉动风的能量大小,又因功率与荷载成平方关系,因此在模拟中必然会存在如何分解功率谱,从而得到荷载表达式的问题。国内外大量研究表明脉动风可作为均值为零的高斯过程及平稳随机过程来考虑,具有很明显的各态历经性,对于n个具有零均值的平稳高斯过程,其谱密度函数矩阵可表示为:式(4)中s11(ω),s22(ω)等可由式(2)求得,且对于Davenport谱有:s11(ω)=s22(ω)=…=snn(ω);对S(ω)做Cholesky分解有:S(ω)=H(ω)·H*(ω)T,式中H*(ω)T为H(ω)的共轭转置。1.4林木结构风压的空间相关性脉动风作用下,在林木结构表面风向及风速在不同位置处并不同步,有的甚至是相互独立无关的,因此,模拟林木结构的脉动风压须考虑其空间相关性。根据Shiotani的建议,在本次模拟研究中仅考虑竖向相关性,相干函数表达式为1.5节点处的脉动风荷载根据空气动力学原理,风压与风速的基本关系为:根据每个节点的受风荷面积Si,乘以节点处的脉动风压时程wi(t)就可得到脉动风荷载的时程曲线也即:Fi=wi(t)·Si。2基于模态分析方法求解非物质劳动过程中的动力特性对于建立林木动力学模型主要有两种,一种是基于理论分析的SIA方法,这种方法用非常简单的图形表示法描述林木的形状,具有应用广泛、不需要电脑辅助计算等优点,缺点是不能适用于所有树种,不能使研究者了解其内部动力机理。另外一种基于数学描述语言精确描述几何形状的方法,称之为L-system方法,它能描述出复杂几何形状的树种,可用有限元方法分析,缺点是求解困难甚是得不到解。鉴于以上两种方法的优缺点,在本研究中采用的是模态分析方法求解,它是介于以上二者之间的方法。模态分析主要用于确定结构体系的振动特性,同时也是其他动力学分析的基础,如谐响应分析、瞬态动力学分析、谱分析等。进行林木结构模态分析的有限元软件采用国内外广泛应用的Ansys。3示例模拟3.1新疆杨树模型的初生长模型本研究模拟实例为北京郊区某公路旁的进行过现场试验的新疆杨,新疆杨树干和树枝结构细长,主干直立向上,呈圆柱形。树根入土深,与土壤有很强的锚固力,抗风性好。常作为风景树、行道树及绿化树。根据新疆杨以上特性,在Ansys模型中其树干及第一主分支采用Beam188模拟,Beam188基于Timoshenko梁结构理论,考虑了剪切变形的影响,适用于分析从细长到中等粗短的梁结构。同时为研究方便,模型中暂不考虑树叶和其它次分枝的影响,根部和土壤之间考虑为固结。树干、树枝质量密度为均匀分布,全树弹性模量为定值。树枝作为悬臂梁固结于树干,而树干也同样作为悬臂梁固结于地面。新疆杨有限元模型及主要参数见表1。3.2基于matlab的模型分析采用线性滤波法原理模拟脉动风,线性滤波法又称白噪声滤波法,它将随机过程抽象为满足一定条件的白噪声,然后经某一假定系统进行适当变换而拟合出具有随机性、时间相关性、空间相关性的风速时程模型。线性滤波法中的自回归(Auto—Regressive,AR)模型因其计算量小、速度快的特点,被广泛应用于随机振动和时间系列分析中。依据文献,[13-15]并结合林木结构的特点,编制了基于Matlab的模拟程序。本研究中脉动风AR模拟使用的主要参数如表2所示。因篇幅限制,只选取了新疆杨模型中2、3结点风速时程曲线,从图2、图3中可以看出不同高度处的脉动风速大小不同,在互不相同的时刻达到各自峰值,体现出空间相关性的影响,其值的变化范围在0附近波动,表现出平稳随机过程的特性。图4中,将节点3的Davenport脉动风速谱与通过模拟获得的脉动风速谱用双对数坐标形式比较,可以看出二者相当吻合,表明了对于给定风速谱密度的随机过程,采用AR模型能够很好的实现人工模拟,并且具有较高的准确性和精度。3.3速度、加速度时程曲线在Matlab中输出各节点风荷载时程数据,施加到Ansys中新疆杨有限元模型相应节点上,并进行瞬态动力学分析,就可以获得新疆杨的各节点处的位移、速度、加速度时程曲线,为使图清晰可见,速度及加速度只提取了前15s的时程曲线。从图6至图8可以看出,在脉动风荷载作用下,新疆杨风致振动的位移、速度、加速度也呈现出明显的随机性特点。4风致随机振动模拟本研究利用AR模拟方法实现了林木脉动风荷载的模拟,而修改AR模拟参数还可以对处于不同地区、不同地貌的林木进行风荷模拟,同时对新疆杨的有限元模型修改物理参数及建模的几何形状后又可以推广到其他的树种,将以上二者结合起来就可以实现不同地区、不同树种的风致随机振动模拟。在有限元模型中考虑树叶、其它次要树枝后,能使模拟更加趋于真实。对于实际为林木采用的防灾减灾措施,如除叶剪枝、施加减震器等,都能用是否建立枝,叶模型、改变林木阻尼(质量阻尼、刚度阻尼)等来实现