CN106055764B - 基于三维壳有限元-梁模型的风力机叶片位移计算方法 - Google Patents
基于三维壳有限元-梁模型的风力机叶片位移计算方法 Download PDFInfo
- Publication number
- CN106055764B CN106055764B CN201610356352.XA CN201610356352A CN106055764B CN 106055764 B CN106055764 B CN 106055764B CN 201610356352 A CN201610356352 A CN 201610356352A CN 106055764 B CN106055764 B CN 106055764B
- Authority
- CN
- China
- Prior art keywords
- blade
- model
- finite element
- section
- displacement
- 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.)
- Expired - Fee Related
Links
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]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/17—Mechanical parametric or variational design
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Geometry (AREA)
- General Physics & Mathematics (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Computer Hardware Design (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Aerodynamic Tests, Hydrodynamic Tests, Wind Tunnels, And Water Tanks (AREA)
Abstract
本发明公开了一种基于三维壳有限元‑梁模型的风力机叶片位移计算方法,在该计算方法中将叶片三维壳有限元模型等效为梁模型,梁模型的节点与截面参考中心重合,对叶片叶根进行完全约束,在壳有限元模型的自由端分别施加6个独立的单位载荷,基于梁模型柔度矩阵的节点位移计算公式,只需要求解一次各梁模型柔度矩阵,在载荷已知情况下,通过做简单的矩阵乘法,可以直接计算出叶片各气动中心的位移,极大提高在叶片气弹耦合、弯扭耦合特性分析计算中叶片弹性变形的计算效率。
Description
技术领域
本发明涉及一种基于三维壳有限元-梁模型的风力机叶片位移计算方法,用于计算叶片等效梁模型的节点位移,同时也用于气弹分析中的位移计算及弯扭耦合特性计算。
背景技术
叶片变形是风力机叶片进行结构设计和气弹耦合分析的重要参数。叶片刚度设计需要准确地计算变形;气弹耦合分析中,需要反复迭代计算叶片变形及气动性能。因此,准确、快速地计算出叶片变形是叶片设计中的重要问题。
风力机叶片具有长厚比、长宽比大的特点,常被简化为梁模型,叶片变形计算采用梁模型时自由度可以很好适应基于叶素动量理论得到的叶片气动载荷,同时梁模型由于具有单元数量少、计算时间短的特点,是目前应用最广泛的描述叶片的模型之一。在风力机的设计过程中,为了分析叶片的动力学行为,需要进行风力机多体系统的气弹分析,采用梁模型使得具有大自由度的多体系统分析成为可能。
现代大型风力机叶片具有复杂的截面形状、多种类的铺层材料、丰富的主梁结构等特点。采用三维壳单元建立叶片模型比采用传统的一维梁单元建立模型能更准确更真实地涵盖大型叶片的这些特点。此外,三维分析结果能够提供叶片更丰富的变形、应变等信息,但是叶片三维模型建立的效率低。
综上所述,在叶片三维壳有限元模型的建模过程中,由于未引入任何假设,因此具有很高的建模精度。但是,叶片壳模型建模过程较为繁琐,特别是需要花费大量的时间设置壳单元的铺层信息,这导致壳有限元模型的建模效率很低。此外,叶片壳有限元模型的自由度总数很大,无法高效地用于风力机多体系统的气弹分析中。因此,如何从三维壳有限元模型数据直接得到多体系统分析所需的梁模型数据非常值得研究。
发明内容
本发明要解决的技术问题是提供一种基于三维壳有限元-梁模型的风力机叶片位移计算方法,该方法基于梁模型柔度矩阵的节点位移计算公式,在已知叶片各叶素上的载荷后,可以直接计算出叶片各参考中心的位移,极大提高在叶片气弹耦合、弯扭耦合特性计算中叶片弹性变形计算效率。
为解决上述技术问题,本发明所采取的技术方案是:1、一种基于三维壳有限元-梁模型的风力机叶片位移计算方法,其特征在于包括如下计算步骤:
1)依据叶片的铺层、翼型和材料数据,建立叶片三维壳有限元模型,壳有限元模型包含n个横向的截面,将叶片三维壳有限元模型等效为梁模型,梁模型的节点与截面参考中心重合,对叶片叶根进行完全约束,在壳有限元模型的自由端分别施加6个独立的单位载荷,其中包括3个相互垂直的力载荷和3个绕力载荷逆时针方向的力矩载荷,形成6种工况;
2)计算每种载荷工况下,叶片各个截面参考中心的位移,
式中,(u,v,w)为叶片壳有限元模型上任意节点(x,y,z)的位移;(U,V,W)为截面参考中心的线位移,(ψx,ψy,ψz)为截面参考中心的角位移;
将各个截面参考中心的位移代入下式,得到每个截面的位移矩阵Di;
3)根据壳有限元模型的自由端施加的单位载荷,计算每个横向截面的载荷矩阵Fi
式中,L为叶片长度,Li为第i个截面距离叶片约束端的距离,i=1~n;
4)根据位移矩阵Di和载荷矩阵Fi,计算出叶片各悬臂梁模型的6×6柔度矩阵Ci
Ci=DiFi -1
5)根据每种工况下作用于叶片叶素的载荷fi、叶片梁模型的柔度矩阵Ci和载荷矩阵Fi,计算得到作用在第j(j=1,2,…,i-1)个截面的气动载荷在i截面引起的刚体位移和作用在第j(j=i,i+1,…,m)个截面的气动载荷作用在i截面引起的变形
式中,Mij为转换矩阵,其计算公式如下
式中,Lj为第j个截面距离叶片约束端的距离,j=1~n;
6)根据刚体位移及变形计算梁模型各节点的位移di
7)重复步骤5)和6),可计算不同载荷下的叶片梁模型节点位移。
对上述方法作进一步补充,所述的叶片三维壳有限元模型包括叶片坐标信息,其中坐标原点位于叶根截面的参考中心,沿叶片长度方向上、各截面参考中心连线为z轴,且沿叶片展向为正,xoz平面位于旋转平面内,y为挥舞方向,x为摆振方向。
对上述方法作进一步补充,所述的步骤4)中,柔度矩阵Ci可表示为
式中:各元素为梁模型在叶片坐标系下的柔度,C11和C22为剪切柔度,C33为拉伸柔度,C44为x向弯曲柔度,C55为y向弯曲柔度,C66为扭转柔度;非对角线元素表示耦合柔度。
采用上述技术方案所产生的有益效果在于:本发明用于风力机叶片位移计算方法,从叶片三维壳有限元模型出发,采用梁模型柔度定义,已知叶片各叶素上的载荷后,可以直接计算出叶片各参考中心的位移,极大提高在叶片气弹耦合、动态分析计算中叶片弹性变形计算效率;在计算过程中,本发明提出方法中,只需要求解一次各悬臂梁模型的柔度矩阵,即可计算不同载荷作用下的梁模型节点位移,计算得到的叶片气动中心变形与三维壳有限元计算得到的叶片气动中心的变形具有相同精度,为叶片的气弹耦合分析中位移的计算或叶片弯扭耦合特性的计算提供便利。
附图说明
下面结合附图和具体实施方式对本发明作进一步详细的说明。
图1是本发明中风力机叶片三维壳有限元模型;
图2是本发明中风力机叶片梁模型节点载荷示意图;
图3是本发明中风力机叶片梁模型节点位移示意图;
图4是风力机叶片叶素载荷示意图;
图5是各向同性箱型叶片主梁的三维有限元模型;
图6是图5中叶片主梁的自由端局部放大图;
图7是主梁梁模型刚度结果对比示意图,其中图7a为扭转刚度结果,图7b为拉伸刚度结果,图7c为弯曲刚度结果;
图8是叶片主梁的载荷简图;
图9是叶片主梁变形对比,其中图9a为挠度对比,图9b为角位移对比;
图10是5MW风力机叶片壳有限元模型;
图11是5MW风力机叶片壳有限元模型在额定风速工况下的变形;
图12是5MW风力机叶片梁模型在额定风速工况下变形,其中图12a为角位移,图12b为线位移。
具体实施方式
本发明具体为一种基于三维壳有限元-梁模型的风力机叶片位移计算方法,用于计算叶片等效梁模型的节点位移,包括如下计算步骤:
1)依据叶片的铺层、翼型和材料数据,建立叶片三维壳有限元模型,壳有限元模型分为n个横向的截面,如附图1所示。xoz平面位于旋转平面内,y为挥舞方向,x为摆振方向,叶片各截面参考中心的连线为z轴,且沿叶片展向为正。叶片壳有限元模型的叶根固定,壳模型上任意节点(x,y,z)的位移为(u,v,w)。将叶片三维壳有限元模型等效为梁模型,如图1中黑色粗实线所示,即梁模型的节点与截面参考中心重合,梁模型位于z轴上。
为了得到叶片等效梁模型第i个节点(i截面参考中心)的位移,i节点到叶根的距离为Li,将从叶根至截面i的叶片壳有限元模型段截取出来,如图2和3所示。那么,该叶片段可以简化为一个悬臂梁模型,见图中黑色粗实线,该悬臂梁模型由一个梁单元组成,且梁单元的左节点固定(位移为0)。
对壳有限元模型的自由端分别施加6个独立的单位载荷,其中包括3个相互垂直的力载荷和3个绕力载荷逆时针方向的力矩载荷,形成6种工况;
假设悬臂梁模型i的右节点(i截面参考中心)的载荷f和位移向量d如图2和3所示,为
fi=[fx fy fz mx my mz]T
di=[U V W ψx ψy ψz]T
式中,(U,V,W)为截面参考中心的线位移,(ψx,ψy,ψz)为参考中心的角位移;(fx,fy,fz)为截面力,(mx,my,mz)为截面力矩。
平衡方程为
di=Cifi
式中,Ci为6×6的梁模型柔度矩阵,
式中:各元素为梁模型在叶片坐标系下的柔度,如C11和C22为剪切柔度,C33为拉伸柔度,C44为x向弯曲柔度,C55为y向弯曲柔度,C66为扭转柔度;非对角线元素表示耦合柔度,如C46为弯曲扭转耦合柔度。
2)假设叶片在变形前后,截面均为平面,且形状保持不变,即忽略面内、面外翘曲。则得到由截面参考中心线位移(U,V,W)和角位移(ψx,ψy,ψz)表示的叶片三维壳有限元模型各截面节点的位移。计算每种载荷工况下,叶片各个截面等效节点的位移,其表达式如下
式中,(u,v,w)为叶片壳有限元模型上任意节点(x,y,z)的位移;(U,V,W)为截面参考中心的线位移,(ψx,ψy,ψz)为截面参考中心的角位移;
在上式中,给出了叶片壳有限元模型某截面节点位移和该截面参考中心位移之间的关系。该方程有6个未知变量而方程数目为3n个,为超静定方程。采用最小二乘法,可求得6个未知变量。
分别计算在这六个独立工况下叶片壳有限元模型的变形,并将各工况下的壳有限元模型节点位移代入上式,计算得到这六个工况下悬臂梁模型的i截面参考中心的位移列向量,组装位移列向量,得到每个截面的位移矩阵Di;
3)根据叶片的三维壳有限元模型,在叶尖施加沿叶片坐标系X、Y和Z轴方向的单位力或单位力矩,则在这六个独立工况下,由静力平衡方程计算得到第i个梁模型自由端的载荷矩阵Fi
式中,L为叶片长度;Li为第i个悬臂梁模型的长度,即第i个截面到叶片根部的距离。
4)根据位移矩阵Di和载荷矩阵Fi,计算出叶片各悬臂梁模型的6×6柔度矩阵Ci
Ci=DiFi -1
5)由叶素动量理论,可以得到叶片在相应载荷工况下的各叶素的气动载荷,如附图4所示,
根据每种工况下作用于叶片叶素的载荷、叶片梁模型的柔度矩阵Ci和载荷矩阵Fi,计算得到作用在第j(j=1,2,…,i-1)个截面的气动载荷在i截面引起的刚体位移和作用在第j(j=i,i+1,…,m)个截面的气动载荷作用在i截面引起的变形
式中,Mij为转换矩阵,其计算公式如下
式中,Lj为第j个截面距离叶片约束端的距离,j=1~n;
6)由叠加原理,根据刚体位移及变形计算梁模型各节点的位移di
7)重复步骤5)和6),可计算不同载荷下的梁模型节点位移。
由以上求解梁模型节点位移的计算步骤可以看出,基于三维壳有限元-梁模型节点位移计算方法,只需要求解一次各梁模型柔度矩阵,在载荷已知情况下,通过做简单的矩阵乘法,即可求得叶片变形。该方法可以快速计算不同载荷下的变形,可用于气弹耦合等分析中。
在气弹分析过程中,需要反复迭代计算叶片在不同载荷作用下的变形,采用给出的方法,只需要求解一次各悬臂梁模型的柔度矩阵(SSM),然后重复步骤5)和6),即可计算不同载荷作用下的梁模型节点位移。因此,本发明的方法具有较高的计算效率。此外,从本发明方法的推导过程可以看出,该方法建立了由叶片三维壳模型等效为叶片梁模型的方法,没有引入假设直接建立梁模型,因此,计算精度较高。此方法定义的叶片变形计算模型与叶片动量叶素理论中的自由度相适应,没有涉及到载荷从一维到三维的转换,非常方便。
下面给出两个计算实例,来说明本发明中计算方法的准确性和计算精度。
实例一为各向同性材料箱型叶片主梁。
各向同性材料箱型叶片主梁的宽度和高度分别为1.5m和3m,厚度为0.05m,长度为30m。箱型叶片主梁材料为钢材,弹性模量为200GPa,泊松比为0.3,密度为7800kg/m3。
采用有限元ANSYS软件建立的根部固定端约束的箱型叶片主梁三维壳有限元模型,如图5所示。其中,为了减小边界效应的影响,在自由端处耦合了一段梁单元,用于施加载荷,如图6所示。
按照本发明的方法,计算得到各梁模型的柔度矩阵C,代入式
中得到梁模型的刚度矩阵S。
为了验证基于三维壳有限元-梁模型风力机叶片位移计算方法的正确性,将刚度矩阵S的各系数与经典梁的刚度矩阵公式进行比较。
经典梁的刚度矩阵公式为
式中,EA为截面拉伸刚度;EIx和EIy为弯曲刚度,GJ为扭转刚度,各刚度的计算公式如下,
EA=E∫dA,EIx=E∫y2dA,EIy=E∫x2dA,
式中,Ω为箱型梁中面所围面积,δ为壁厚。
将本发明的方法计算得到的刚度矩阵与经典梁的刚度矩阵计算得到的结果相比较,得到
将本发明的方法计算得到的刚度矩阵与经典梁的刚度矩阵计算得到S11、S22、S33和S66进行比较,如图7所示。
从图7可以看出,采用本发明计算方法与经典梁的刚度矩阵计算得到的拉伸、扭转刚度几乎相同,误差均小于5%。其中由于边界效应的影响,弯曲刚度在约束部位附近误差大于其他截面。
箱型叶片主梁受到载荷如图8所示,采用本发明的计算步骤,计算得到叶片各截面变形。将计算结果与梁模型、壳有限元模型、材料力学解析解进行比较如图9所示。
从图9中可以看出,各种计算方法得到的挠度和弯曲角的变化趋势是相同的。采用该方法计算得到的叶片气动中心变形与三维壳有限元模型计算得到的叶片气动中心的变形具有相同精度,验证了计算方法的正确性。
实例二为5MW风力机叶片变形计算。
建立的5MW风力机叶片模型,5MW叶片由蒙皮、梁帽和腹板组成,叶片长度为61.5m,由8种不同的翼型组成,在20%附近的位置达到最大弦长,为4.7m。
表1给出了5MW叶片所用复合材料的力学特性。采用ANSYS软件的shell99单元建立的三维叶片壳有限元模型如图10所示,计算了5MW叶片在额定风速为11.4m/s工况下的叶片壳有限元模型的位移如图11所示。
表1铺层材料的力学特性
按照位移计算步骤,计算了5MW叶片在额定风速为11.4m/s工况下的叶片梁模型的位移如图12所示。从图12可以看出,在11.4m/s风速下,叶片拉伸变形及扭转变形均很小,通常忽略不计;叶尖挥舞方向变形为4.6m,达到了合理安全距离,避免了与塔架发生碰撞。
由上述两个实例可知,本发明的计算方法定义的叶片变形计算模型与叶片动量叶素理论中的自由度相适应,没有涉及到载荷从一维到三维转换,非常方便。而且该方法基于叶片三维模型,未引入各种假设直接建立梁模型,因此,计算结果比传统方法精度更高。
本发明的计算方法兼具了三维壳有限元模型的高精度和叶片梁模型的高效率。在已知叶片各叶素上的载荷的情况下,利用本计算方法可以直接计算出叶片梁模型节点的位移。为叶片的气弹耦合分析中气动载荷的计算或叶片弯扭耦合特性的计算提供便利。
Claims (3)
1.一种基于三维壳有限元-梁模型的风力机叶片位移计算方法,其特征在于包括如下计算步骤:
1)依据叶片的铺层、翼型和材料数据,建立叶片三维壳有限元模型,壳有限元模型包含n个横向的截面,将叶片三维壳有限元模型等效为梁模型,梁模型的节点与截面参考中心重合,对叶片叶根进行完全约束,在壳有限元模型的自由端分别施加6个独立的单位载荷,其中包括3个相互垂直的力载荷和3个绕力载荷逆时针方向的力矩载荷,形成6种工况;
2)计算每种载荷工况下,叶片各个截面参考中心的位移,
式中,(u,v,w)为叶片壳有限元模型上任意节点(x,y,z)的位移;(U,V,W)为截面参考中心的线位移,(ψx,ψy,ψz)为截面参考中心的角位移;
将各个截面参考中心的位移代入下式,得到每个截面的位移矩阵Di;
3)根据壳有限元模型的自由端施加的单位载荷,计算每个横向截面的载荷矩阵Fi
式中,L为叶片长度,Li为第i个截面距离叶片约束端的距离,i=1~n;
4)根据位移矩阵Di和载荷矩阵Fi,计算出叶片各悬臂梁模型的6×6柔度矩阵Ci
Ci=DiFi -1
5)根据每种工况下作用于叶片叶素的载荷fi、叶片梁模型的柔度矩阵Ci和载荷矩阵Fi,计算得到作用在第j(j=1,2,…,i-1)个截面的气动载荷在i截面引起的刚体位移和作用在第j(j=i,i+1,…,m)个截面的气动载荷作用在i截面引起的变形
式中,Mij为转换矩阵,其计算公式如下
式中,Lj为第j个截面距离叶片约束端的距离,j=1~n;
6)根据刚体位移及变形计算梁模型各节点的位移di
7)重复步骤5)和6),可计算不同载荷下的叶片梁模型节点位移。
2.根据权利要求1所述的基于三维壳有限元-梁模型的风力机叶片位移计算方法,其特征在于所述的叶片三维壳有限元模型包括叶片坐标信息,其中坐标原点位于叶根截面的参考中心,沿叶片长度方向上、各截面参考中心连线为z轴,且沿叶片展向为正,xoz平面位于旋转平面内,y为挥舞方向,x为摆振方向。
3.根据权利要求2所述的基于三维壳有限元-梁模型的风力机叶片位移计算方法,其特征在于所述的步骤4)中,柔度矩阵Ci可表示为
式中:各元素为梁模型在叶片坐标系下的柔度,C11和C22为剪切柔度,C33为拉伸柔度,C44为x向弯曲柔度,C55为y向弯曲柔度,C66为扭转柔度;非对角线元素表示耦合柔度。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610356352.XA CN106055764B (zh) | 2016-05-26 | 2016-05-26 | 基于三维壳有限元-梁模型的风力机叶片位移计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610356352.XA CN106055764B (zh) | 2016-05-26 | 2016-05-26 | 基于三维壳有限元-梁模型的风力机叶片位移计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106055764A CN106055764A (zh) | 2016-10-26 |
CN106055764B true CN106055764B (zh) | 2019-04-16 |
Family
ID=57174501
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610356352.XA Expired - Fee Related CN106055764B (zh) | 2016-05-26 | 2016-05-26 | 基于三维壳有限元-梁模型的风力机叶片位移计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106055764B (zh) |
Families Citing this family (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107490345A (zh) * | 2017-07-26 | 2017-12-19 | 中国电建集团西北勘测设计研究院有限公司 | 一种基于三维激光扫描的高耸塔筒弯曲度检测方法 |
CN109299492B (zh) * | 2018-05-23 | 2023-04-18 | 中国航空工业集团公司沈阳飞机设计研究所 | 一种验证柔性双圆管翼梁扭转变形的方法 |
US11714938B2 (en) | 2019-08-16 | 2023-08-01 | General Electric Company | Method for distortion prediction and compensation for sintering of binder jet printed parts |
CN111814265B (zh) * | 2020-05-27 | 2023-12-22 | 明阳智慧能源集团股份公司 | 一种风电叶片质量计算方法 |
CN111814225B (zh) * | 2020-06-04 | 2022-12-23 | 中国铁路设计集团有限公司 | 一种基于广义柔度矩阵的外荷载墩顶位移计算方法 |
CN111832211B (zh) * | 2020-07-27 | 2023-07-07 | 内蒙古工业大学 | 一种复合纤维风力机叶片的刚度优化方法 |
CN112052524B (zh) * | 2020-09-25 | 2022-09-06 | 中国直升机设计研究所 | 一种直升机吊挂柔性机身建模方法 |
CN113219926A (zh) * | 2021-05-13 | 2021-08-06 | 中国计量大学 | 基于数字孪生系统的人机共融制造单元安全风险评估方法 |
CN114199500A (zh) * | 2021-07-14 | 2022-03-18 | 南京工业大学 | 大型风力机叶片实尺寸模拟风力加载测试系统 |
CN113536648B (zh) * | 2021-09-06 | 2021-12-17 | 广东工业大学 | 一种海上风电平台的模拟仿真可对中匹配计算方法及系统 |
CN115758830A (zh) * | 2022-11-23 | 2023-03-07 | 中国核电工程有限公司 | 一种管道力学计算方法及装置 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105181318A (zh) * | 2015-09-23 | 2015-12-23 | 华北电力大学(保定) | 一种风力发电机叶片弯扭耦合向量测量装置 |
-
2016
- 2016-05-26 CN CN201610356352.XA patent/CN106055764B/zh not_active Expired - Fee Related
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105181318A (zh) * | 2015-09-23 | 2015-12-23 | 华北电力大学(保定) | 一种风力发电机叶片弯扭耦合向量测量装置 |
Non-Patent Citations (4)
Title |
---|
5 MW 风力机叶片弯扭耦合特性分析;赵鹤翔等;《华北电力大学学报》;20130731;全文 |
Warping Effect to Bend-twist Coupling Characteristics for Box Spar of Wind Turbine Blade;Xing-Yin ZHOU等;《International Conference on Material Science and Application (ICMSA 2015)》;20151231;全文 |
航空发动机压气机叶片的静力学及模态分析;曲文浩;《中国优秀硕士学位论文全文数据库 工程科技II辑》;20140715;全文 |
销钉连接叶片固有特性分析有限元模型研究;陈军;《计算机仿真》;20091130;全文 |
Also Published As
Publication number | Publication date |
---|---|
CN106055764A (zh) | 2016-10-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106055764B (zh) | 基于三维壳有限元-梁模型的风力机叶片位移计算方法 | |
Tsushima et al. | Geometrically nonlinear static aeroelastic analysis of composite morphing wing with corrugated structures | |
Guo | Aeroelastic optimization of an aerobatic aircraft wing structure | |
Yu et al. | Generalized Timoshenko theory of the variational asymptotic beam sectional analysis | |
Lee et al. | Multi‐flexible‐body dynamic analysis of horizontal axis wind turbines | |
CN112560167B (zh) | 机翼结构力学高保真降阶仿真方法、电子设备、存储介质 | |
Buckney et al. | Structural efficiency of a wind turbine blade | |
Rezaei et al. | Development of a reduced order model for nonlinear analysis of the wind turbine blade dynamics | |
Xie et al. | Static aeroelastic analysis including geometric nonlinearities based on reduced order model | |
CN112580241B (zh) | 一种基于结构降阶模型的非线性气动弹性动响应分析方法 | |
La Magna | Bending-active plates: strategies for the induction of curvature through the means of elastic bending of plate-based structures | |
CN113886967A (zh) | 多巡航工况的大型飞机机翼气动弹性优化方法 | |
Zhou et al. | Geometrically nonlinear beam analysis of composite wind turbine blades based on quadrature element method | |
Hodges | Unified approach for accurate and efficient modeling of composite rotor blade dynamics The Alexander A. Nikolsky honorary lecture | |
Invernizzi et al. | A fully consistent linearized model for vibration analysis of rotating beams in the framework of geometrically exact theory | |
Meddaikar et al. | Optimization, manufacturing and testing of a composite wing with maximized tip deflection | |
CN116305589B (zh) | 一种直升机桨叶结构降阶分析方法、系统、设备及介质 | |
CN116484671A (zh) | 风力机叶片全耦合动力响应精细化分析方法 | |
Fulton et al. | Aeroelastic stability of composite hingeless rotor blades in hover—Part I: theory | |
Navadeh et al. | Approximate Mode-Based Simulation of Composite Wind Turbine Blade Vibrations Using a Simplified Beam Model | |
Lin et al. | Combined analytical and finite element beam model for wind turbine blades | |
Barnes et al. | Evaluation of simplified loading models for finite element analysis of composite wind turbine blades | |
Tsushima et al. | Structural and aerodynamic models for aeroelastic analysis of corrugated morphing wings | |
Otero et al. | Structural analysis of complex wind turbine blades: flexo-torsional vibrational modes | |
Datta et al. | Integrated Aeromechanics with Three-Dimensional Solid-Multibody Structures |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20190416 |