CN110866360B - 一种旋转陶瓷基复合材料叶片的模态计算方法 - Google Patents
一种旋转陶瓷基复合材料叶片的模态计算方法 Download PDFInfo
- Publication number
- CN110866360B CN110866360B CN201911130769.4A CN201911130769A CN110866360B CN 110866360 B CN110866360 B CN 110866360B CN 201911130769 A CN201911130769 A CN 201911130769A CN 110866360 B CN110866360 B CN 110866360B
- Authority
- CN
- China
- Prior art keywords
- matrix
- centrifugal force
- displacement
- stiffness
- stress
- 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.)
- Active
Links
Images
Abstract
一种旋转陶瓷基复合材料叶片的模态计算方法,包括如下步骤:1、对陶瓷基复合材料叶片划分有限元网格,并获得单元、结点信息;2、计算单元动能、应变能,获得离心力载荷下有限元形式的运动控制方程;3、基于损伤组装结构的弹性刚度矩阵、初应力刚度矩阵、离心力刚度矩阵,进行离心载荷下的静力分析获得预应力;4、将预应力加载到结构上进行模态分析,求解广义特征值问题获得离心力载荷下的模态。本发明综合考虑了由弹性模量下降引起的结构软化和离心力引起的结构刚度加强现象,运用牛顿‑拉斐逊法,并在每个迭代步中进行刚度的实时折减,能够给出离心力载荷下陶瓷基复合材料叶片的模态响应和应力、应变响应。
Description
技术领域
本发明属于复合材料力学分析技术领域,具体涉及一种旋转陶瓷基复合材料叶片的模态计算方法。
背景技术
陶瓷基复合材料因具有耐高温、密度小等优点,成为未来航空材料的候选者。转子叶片作为航空发动机典型热端部件,运用陶瓷基复合材料可以降低重量并承受更高的涡轮前进口温度。转子叶片在实际工况下,需要承受很大的离心力。一方面,离心力增强了叶片的弹性恢复力,使得叶片各阶固有频率提高。另一方面,陶瓷基复合材料的非线性本构造成旋转过程中的损伤,使得旋转叶片的模态计算区别于普通金属材料。
叶片受离心力作用而产生径向拉应力,在模态分析中将其作预应力处理。事实上,复合材料目前更广泛地被应用在风力机叶片上。Yanbin Cui(Cui Y,Shi L,Zhao F.ModalAnalysis of Wind Turbine Blade Made of Composite laminated plates[C]//2010Asia-Pacific Power and Energy Engineering Conference.IEEE,2010.)、Zhou,LiQun(Zhou L Q,Xing S H,Li Y P.Dynamic Analysis for Wind Turbine CompositeBlade[J].Applied Mechanics and Materials,2013,364:102-106.)、XING ShuaiHeng(Shuaiheng X,Liqun Z,Yuping L I,et al.WIND TURBINE COMPOSITE BLADE DYNAMICCHARACTER RESEARCH[J].Journal of Mechanical Strength,2014.)等人均对复合材料风力机叶片进行了模态计算,并分析了铺层角度的影响。ZHENG Guang-hua(Guang-Hua Z,Jia-Lin L I,Xue T,et al.Vibration Modal Analysis of Compressor Blade based onANSYS[J].International Journal of Plant Engineering and Management,2017,22(2).)基于ANSYS软件对压气机叶片进行了强度和模态计算,分析了离心力对模态参数的影响。Qiang Yang (Yang Q,Fang K Q,Shao C.Vibration Characteristics Analysis ofAeroengine Composite Blade[J].Advanced Materials Research,2012,583:57-61.)使用非接触式激光扫描振动计测量了航空发动机复合材料叶片模态响应,并与FEM结果进行了对比分析。目前的工作大多基于有限元软件进行预应力模态计算,研究离心力对叶片模态的影响,但是均未定性、定量地考虑复合材料由于弹性模量下降引起的结构软化现象。因此,确有必要对现有技术进行改进以解决现有技术之不足。
发明内容
本发明针对现有技术中的不足,提供一种旋转陶瓷基复合材料叶片的模态计算方法。
为实现上述目的,本发明采用以下技术方案:
一种旋转陶瓷基复合材料叶片的模态计算方法,其特征在于,包括如下步骤:
步骤1、对陶瓷基复合材料叶片划分有限元网格,并获得单元、结点信息;
步骤2、计算单元动能、应变能,获得离心力载荷下有限元形式的运动控制方程;
步骤3、基于损伤组装结构的弹性刚度矩阵、初应力刚度矩阵、离心力刚度矩阵,以及步骤2计算的运动控制方程,进行离心载荷下的静力分析获得预应力;
步骤4、将预应力加载到结构上进行模态分析,求解广义特征值问题获得离心力载荷下的模态;模态包括固有频率ω和固有振型向量φ。
为优化上述技术方案,采取的具体措施还包括:
进一步地,步骤1中,对陶瓷基复合材料叶片进行有限元网格划分,并提取有限元程序可识别的单元、结点信息;基于有限元原理,组装与应力状态无关的形函数矩阵[N]、质量矩阵[M]、等效离心力载荷{Qc}等。
进一步地,步骤2中,计算单元动能Te和应变能Ue分别为
代入拉格朗日方程
式中,T表示总体动能,是系统用广义坐标表示的动能,t表示时间,{Q}表示广义力,{δ}表示广义位移,即叶片上各节点在离心力作用下的位移;
化简获得有限元形式的运动控制方程
式中,{Qc}表示离心力荷载向量;反对称陀螺阵[MG]与旋转角速度Ω一次方成正比,对前几阶频率影响较小,可忽略;[Ke]是系统弹性刚度矩阵,包括小位移刚度矩阵[K0]和大变形刚度矩阵[Kn],分别为
式中,[B0]是线性几何矩阵,与节点位移无关;[Bn]是非线性几何矩阵,是节点位移的函数;e表示单元数;因为是用有限元求解,所以单元的各矩阵需要组装成结构的总体矩阵;所以,这边的方程表示在每个单元内积分获得单元的各矩阵,然后组装成总体矩阵;
初应力刚度矩阵[Kσ]是由离心载荷产生的附加刚度矩阵,导致应力刚化现象,由下式给出
式中,矩阵[G]是形函数对坐标的偏导组装的矩阵,与位移无关;矩阵[M]是9×9由六个应力分量组成的矩阵,与应力直接相关;
离心力刚度矩阵[Kc]是由旋转引起的对称阵,与旋转角速度Ω的二次方成正比,由下式给出
式中,[S]是由单位旋转角速度向量在x、y、z轴上的投影组成的矩阵;
([Ke]+[Kσ]-[Kc]){δ}={Qc} (9)
运用牛顿-拉斐逊迭代法求解上式,获得预应力。
进一步地,步骤3中,运用牛顿-拉斐逊迭代法求解上式,具体的迭代过程分成四步:
3-1)选取初始位移{δ}0=[K0]-1{Qc},并计算出离心力刚度矩阵[Kc];
32)根据初始位移{δ}0,依次计算出当前应力状态下的初应力刚度矩阵[Kσ]和大变形刚度矩阵[Kn];预应力则是由离心载荷引起的,并以初应力刚度矩阵[Kσ]的形式,作为附加刚度矩阵出现在控制方程中;
3-3)运用下式进行迭代求解
{δ}i=([Ke]i-1+[Kσ]i-1-[Kc])-1{Qc} (10)
式中,i表示迭代序号,第i步的位移是由上一步的计算结果得到的;
3-4)规定{δ}i≈{δ}i-1时计算收敛,{δ}i即为旋转态叶片平衡位置的位移;
其中,每一迭代步均考虑弹性模量降低引起的结构软化,区别于一般的线性材料;根据每一步的位移依次计算出当前状态下的应力、应变;进一步获得割线刚度,用于计算各刚度矩阵中反映材料性能的弹性矩阵[D],实现刚度的实时折减。
进一步地,步骤4中,根据步骤3的结果获得当前转速下的离心应力分布,并以预应力形式对叶片进行模态计算
求解广义特征值问题
([Ke]+[Kσ]-[Kc]-ω2[M]){φ}=0 (12)
式中,ω是模态中的固有频率,{φ}是模态中的固有振型向量,即可获得离心力载荷下陶瓷基复合材料叶片的模态。
本发明的有益效果是:本发明提供的一种旋转陶瓷基复合材料叶片的模态计算方法,综合考虑了由弹性模量下降引起的结构软化和离心力引起的结构刚度加强现象。运用牛顿-拉斐逊法,并在每个迭代步中进行刚度的实时折减,能够给出离心力载荷下陶瓷基复合材料叶片的模态响应和应力、应变响应。相比于现有的计算方法,本发明考虑了陶瓷基复合材料的非线性本构对叶片模态的影响。
附图说明
图1本发明中陶瓷基复合材料叶片模型。
图2本发明中陶瓷基复合材料叶片第一阶固有频率与转速关系。
具体实施方式
现在结合附图对本发明作进一步详细的说明。
实施例:对如图1所示的陶瓷基复合材料叶片施加离心载荷,转速Ω从0增大到2000rad/s。
步骤1:对陶瓷基复合材料叶片进行有限元网格划分,并提取有限元程序可识别的单元、结点等信息。基于有限元原理,组装与应力状态无关的形函数矩阵[N]、质量矩阵[M]、等效离心力载荷{Qc}等。
其中叶片绕着Y轴旋转,离心力的等效节点载荷由下式计算
步骤2:推导有限元形式的运动控制方程,并获得方程中各刚度矩阵的计算公式,如式(4)-(8)所示。
步骤3:进行离心载荷下的静力分析获得预应力。静力学问题的控制方程如式(9)所示。运用牛顿-拉斐逊迭代法求解上式,具体的迭代过程分成四步:
3-1)选取初始位移{δ}0=[K0]-1{Qc},并计算出离心力刚度矩阵[Kc];
3-2)根据初始位移{δ}0,依次计算出当前应力状态下的初应力刚度矩阵[Kσ]和大变形刚度矩阵[Kn];
3-3)运用式(10)进行迭代求解;
3-4)规定{δ}i≈{δ}i-1时计算收敛,{δ}i即为旋转态叶片平衡位置的位移。
上述步骤中,每一迭代步均考虑弹性模量降低引起的结构软化,区别于一般的线性材料。根据每一步中的位移可依次计算出当前状态下的应力、应变。进一步获得割线刚度,用于计算各刚度矩阵中反映材料性能的弹性矩阵[D],实现刚度的实时折减。
步骤4:根据步骤3中的结果获得当前转速下的离心应力分布,并以预应力形式对叶片进行模态分析。求解广义特征值问题获得离心力载荷下陶瓷基复合材料叶片的模态。
将考虑材料非线性的计算结果和只考虑离心力刚化作用的计算结果绘制在图2中。横坐标为转速;纵坐标为第一阶固有频率。由图2可知,随着离心载荷的增大,第一阶固有频率ω1均增大,表明固有频率受叶片旋转刚度加强效应的影响很大。在较小转速下,材料的损伤不明显,两条曲线几乎无差别。随着转速的增大材料损伤程度加剧,会发现固有频率增大的程度被减弱,是非线性本构造成的结构软化现象的表现。
需要注意的是,发明中所引用的如“上”、“下”、“左”、“右”、“前”、“后”等的用语,亦仅为便于叙述的明了,而非用以限定本发明可实施的范围,其相对关系的改变或调整,在无实质变更技术内容下,当亦视为本发明可实施的范畴。
以上仅是本发明的优选实施方式,本发明的保护范围并不仅局限于上述实施例,凡属于本发明思路下的技术方案均属于本发明的保护范围。应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明原理前提下的若干改进和润饰,应视为本发明的保护范围。
Claims (2)
1.一种旋转陶瓷基复合材料叶片的模态计算方法,其特征在于,包括如下步骤:
步骤1、对陶瓷基复合材料叶片进行有限元网格划分,并提取有限元程序可识别的单元、结点信息;基于有限元原理,组装与应力状态无关的形函数矩阵[N]、质量矩阵[M]、等效离心力载荷{Qc};
步骤2、计算单元动能、应变能,获得离心力载荷下有限元形式的运动控制方程;计算单元动能Te和应变能Ue分别为
代入拉格朗日方程
式中,T表示总体动能,是系统用广义坐标表示的动能,t表示时间,{Q}表示广义力,{δ}表示广义位移,即叶片上各节点在离心力作用下的位移;
化简获得有限元形式的运动控制方程
式中,{Qc}表示离心力荷载向量;反对称陀螺阵[MG]与旋转角速度Ω一次方成正比;[Ke]是系统弹性刚度矩阵,包括小位移刚度矩阵[K0]和大变形刚度矩阵[Kn],分别为
式中,[B0]是线性几何矩阵,与节点位移无关;[Bn]是非线性几何矩阵,是节点位移的函数;e表示单元数;
初应力刚度矩阵[Kσ]是由离心载荷产生的附加刚度矩阵,由下式给出
式中,矩阵[G]是形函数对坐标的偏导组装的矩阵,与位移无关;矩阵[M]是9×9由六个应力分量组成的矩阵,与应力直接相关;
离心力刚度矩阵[Kc]是由旋转引起的对称阵,与旋转角速度Ω的二次方成正比,由下式给出
式中,[S]是由单位旋转角速度向量在x、y、z轴上的投影组成的矩阵;
步骤3、基于损伤组装结构的弹性刚度矩阵、初应力刚度矩阵、离心力刚度矩阵,以及步骤2计算的运动控制方程,进行离心载荷下的静力分析获得预应力;其中,进行离心载荷下的静力分析,对于静力学问题,位移{δ}与时间无关,则其对时间的一阶和二阶偏导均为0,即代入式(4)得
([Ke]+[Kσ]-[Kc]){δ}={Qc} (9)
运用牛顿-拉斐逊迭代法求解上式,获得预应力;具体的迭代过程分成四步:
3-1)选取初始位移{δ}0=[K0]-1{Qc},并计算出离心力刚度矩阵[Kc];
3-2)根据初始位移{δ}0,依次计算出当前应力状态下的初应力刚度矩阵[Kσ]和大变形刚度矩阵[Kn];预应力则是由离心载荷引起的,并以初应力刚度矩阵[Kσ]的形式,作为附加刚度矩阵出现在控制方程中;
3-3)运用下式进行迭代求解
{δ}i=([Ke]i-1+[Kσ]i-1-[Kc])-1{Qc} (10)
式中,i表示迭代序号,第i步的位移是由上一步的计算结果得到的;
3-4)规定{δ}i≈{δ}i-1时计算收敛,{δ}i即为旋转态叶片平衡位置的位移;
其中,每一迭代步均考虑弹性模量降低引起的结构软化;根据每一步的位移依次计算出当前状态下的应力、应变;进一步获得割线刚度,用于计算各刚度矩阵中反映材料性能的弹性矩阵[D],实现刚度的实时折减;
步骤4、将预应力加载到结构上进行模态分析,求解广义特征值问题获得离心力载荷下的模态。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911130769.4A CN110866360B (zh) | 2019-11-19 | 2019-11-19 | 一种旋转陶瓷基复合材料叶片的模态计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911130769.4A CN110866360B (zh) | 2019-11-19 | 2019-11-19 | 一种旋转陶瓷基复合材料叶片的模态计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110866360A CN110866360A (zh) | 2020-03-06 |
CN110866360B true CN110866360B (zh) | 2021-08-24 |
Family
ID=69655189
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201911130769.4A Active CN110866360B (zh) | 2019-11-19 | 2019-11-19 | 一种旋转陶瓷基复合材料叶片的模态计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110866360B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114970033B (zh) * | 2022-06-07 | 2023-09-19 | 中国核工业二三建设有限公司 | 一种大型薄壁设备吊装过程快速有限元求解方法及系统 |
CN116050027B (zh) * | 2023-03-30 | 2023-07-07 | 陕西空天信息技术有限公司 | 叶轮叶片结构静态分析方法、存储介质和电子设备 |
Citations (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101620642A (zh) * | 2009-04-09 | 2010-01-06 | 中国人民解放军国防科学技术大学 | 一种不同形状悬臂梁压电振子振动发电分析有限元方法 |
JP2010258188A (ja) * | 2009-04-24 | 2010-11-11 | Kyocera Corp | 圧電磁器の物理量算出方法 |
CN102570368A (zh) * | 2012-01-16 | 2012-07-11 | 南京航空航天大学 | 一种基于面内或面外模态的行波型压电材料振动防除冰装置及方法 |
CN102579142A (zh) * | 2005-01-25 | 2012-07-18 | 奥美科股份有限公司 | 用于成形生坯体的方法和由该方法制成的制品 |
CN104850713A (zh) * | 2015-05-28 | 2015-08-19 | 西北工业大学 | 机械结构随机振动动态应力高精度计算方法 |
US9274036B2 (en) * | 2013-12-13 | 2016-03-01 | King Fahd University Of Petroleum And Minerals | Method and apparatus for characterizing composite materials using an artificial neural network |
CN106706241A (zh) * | 2016-12-30 | 2017-05-24 | 中国华能集团清洁能源技术研究院有限公司 | 一种风力机叶片损伤主动自检装置及方法 |
US10118696B1 (en) * | 2016-03-31 | 2018-11-06 | Steven M. Hoffberg | Steerable rotating projectile |
CN109858171A (zh) * | 2019-02-14 | 2019-06-07 | 南京航空航天大学 | 编织陶瓷基复合材料应力-应变响应和强度的预测方法 |
CN109920495A (zh) * | 2019-03-28 | 2019-06-21 | 南京航空航天大学 | 一种编织陶瓷基复合材料强度的多尺度预测方法 |
JP2019126749A (ja) * | 2018-01-22 | 2019-08-01 | Towa株式会社 | 加工装置、及び製品の製造方法 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
ITTV20130048A1 (it) * | 2013-04-10 | 2014-10-11 | Dario Toncelli | Macchina per il taglio di materiale in lastra |
US9460259B2 (en) * | 2014-08-22 | 2016-10-04 | Samsung Electronics Co., Ltd. | Methods of generating integrated circuit layout using standard cell library |
-
2019
- 2019-11-19 CN CN201911130769.4A patent/CN110866360B/zh active Active
Patent Citations (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102579142A (zh) * | 2005-01-25 | 2012-07-18 | 奥美科股份有限公司 | 用于成形生坯体的方法和由该方法制成的制品 |
CN101620642A (zh) * | 2009-04-09 | 2010-01-06 | 中国人民解放军国防科学技术大学 | 一种不同形状悬臂梁压电振子振动发电分析有限元方法 |
JP2010258188A (ja) * | 2009-04-24 | 2010-11-11 | Kyocera Corp | 圧電磁器の物理量算出方法 |
CN102570368A (zh) * | 2012-01-16 | 2012-07-11 | 南京航空航天大学 | 一种基于面内或面外模态的行波型压电材料振动防除冰装置及方法 |
US9274036B2 (en) * | 2013-12-13 | 2016-03-01 | King Fahd University Of Petroleum And Minerals | Method and apparatus for characterizing composite materials using an artificial neural network |
CN104850713A (zh) * | 2015-05-28 | 2015-08-19 | 西北工业大学 | 机械结构随机振动动态应力高精度计算方法 |
US10118696B1 (en) * | 2016-03-31 | 2018-11-06 | Steven M. Hoffberg | Steerable rotating projectile |
CN106706241A (zh) * | 2016-12-30 | 2017-05-24 | 中国华能集团清洁能源技术研究院有限公司 | 一种风力机叶片损伤主动自检装置及方法 |
JP2019126749A (ja) * | 2018-01-22 | 2019-08-01 | Towa株式会社 | 加工装置、及び製品の製造方法 |
CN109858171A (zh) * | 2019-02-14 | 2019-06-07 | 南京航空航天大学 | 编织陶瓷基复合材料应力-应变响应和强度的预测方法 |
CN109920495A (zh) * | 2019-03-28 | 2019-06-21 | 南京航空航天大学 | 一种编织陶瓷基复合材料强度的多尺度预测方法 |
Non-Patent Citations (4)
Title |
---|
《Design of Distributed Support Decision-Making System Model Based on Multi-agent Technique for Sporting Event》;Song Yingdong;《2014 Fifth International Conference on Intelligent Systems Design and Engineering Applications》;20141231;356-359 * |
《Numerical simulation of turbine cascade flow with blade-fluid heat exchange》;Petr Louda等;《Applied Mathematics and Computation》;20131231;7206-7214 * |
《Predicting the thermal conductiivty and temperature distribution in 3D braided composites Composite Structures》;Jiang L L等;《Composite Structures》;20141231;578-583 * |
《编织复合材料结构与材料一体化优化设计》;孙杰;《中国博士学位论文全文数据库 工程科技Ⅰ辑》;20120715;B020-11 * |
Also Published As
Publication number | Publication date |
---|---|
CN110866360A (zh) | 2020-03-06 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Chen et al. | Vibration and damping analysis of the bladed disk with damping hard coating on blades | |
Liu et al. | Aeroelastic analysis of a floating offshore wind turbine in platform‐induced surge motion using a fully coupled CFD‐MBD method | |
Riziotis et al. | Aeroelastic stability of wind turbines: the problem, the methods and the issues | |
Kumar et al. | Free vibration analysis of Al 2024 wind turbine blade designed for Uttarakhand region based on FEA | |
Xiong et al. | Dynamic response analysis of the rotating blade of horizontal axis wind turbine | |
Ebrahimi et al. | Transient response of the flexible blade of horizontal-axis wind turbines in wind gusts and rapid yaw changes | |
CN110866360B (zh) | 一种旋转陶瓷基复合材料叶片的模态计算方法 | |
Tsai | Rotating vibration behavior of the turbine blades with different groups of blades | |
Lee et al. | Aeroelastic analysis of wind turbine blades based on modified strip theory | |
Chetan et al. | Flutter predictions in the design of extreme-scale segmented ultralight morphing rotor blades | |
Nyssen et al. | Experimental and numerical characterization of a ceramic matrix composite shroud segment under impact loading | |
CN114357847B (zh) | 一种带冠叶片非线性模态分析方法、装置和设备 | |
Wang et al. | A study of nonlinear aeroelastic response of a long flexible blade for the horizontal axis wind turbine | |
Prasad et al. | Analysis of classical flutter in steam turbine blades using reduced order aeroelastic model | |
Kwon et al. | Structural dynamic modeling and stability of a rotating blade under gravitational force | |
Vanti et al. | An integrated numerical procedure for flutter and forced response assessment of turbomachinery blade-rows | |
Guruswamy | Computational-fluid-dynamics-and computational-structural-dynamics-based time-accurate aeroelasticity of helicopter rotor blades | |
Gantasala et al. | Aeroelastic simulations of wind turbine using 13 DOF rigid beam model | |
Das et al. | Large amplitude free vibration analysis of a rotating beam with non-linear spring and mass system | |
Manolas et al. | Cross comparison of aeroelastic state-of-the-art design tools on a 10 MW scale wind turbine | |
Sayed et al. | Aero-elastic analysis and classical flutter of a multi-megawatt slender bladed horizontal-axis wind turbine | |
Naung | Computational method for aerodynamic and aeromechanical analysis of offshore wind turbines | |
CN104199998A (zh) | 一种风力发电机组桨叶动态特性分析方法 | |
Möhren et al. | Prediction and Comparison of Dynamic Loads of Propellers | |
Lian et al. | Insight on aerodynamic damping of the civil transonic fan blade |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |