CN107292054B - 一种复杂船舶推进轴系扭转振动特性分析方法 - Google Patents
一种复杂船舶推进轴系扭转振动特性分析方法 Download PDFInfo
- Publication number
- CN107292054B CN107292054B CN201710563535.3A CN201710563535A CN107292054B CN 107292054 B CN107292054 B CN 107292054B CN 201710563535 A CN201710563535 A CN 201710563535A CN 107292054 B CN107292054 B CN 107292054B
- Authority
- CN
- China
- Prior art keywords
- propulsion shafting
- complex
- vibration
- shafting
- torsional
- 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
- 238000000034 method Methods 0.000 title claims abstract description 32
- 230000007246 mechanism Effects 0.000 claims abstract description 19
- 230000004044 response Effects 0.000 claims abstract description 14
- 238000006073 displacement reaction Methods 0.000 claims description 21
- 238000013461 design Methods 0.000 claims description 5
- 238000004364 calculation method Methods 0.000 abstract description 19
- 238000010586 diagram Methods 0.000 description 8
- 230000008569 process Effects 0.000 description 6
- 230000008859 change Effects 0.000 description 4
- 230000009471 action Effects 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 3
- 230000007547 defect Effects 0.000 description 2
- 238000005381 potential energy Methods 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 238000002485 combustion reaction Methods 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 238000000556 factor analysis Methods 0.000 description 1
- 230000009347 mechanical transmission Effects 0.000 description 1
- 238000012805 post-processing Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000000926 separation method Methods 0.000 description 1
- 238000001228 spectrum Methods 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
- 238000011426 transformation method Methods 0.000 description 1
Images
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
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Pharmaceuticals Containing Other Organic And Inorganic Compounds (AREA)
- Vibration Prevention Devices (AREA)
Abstract
本发明提供一种复杂船舶推进轴系扭转振动特性分析方法,包括以下步骤:测定复杂船舶推进轴系基本参数;结合边界条件求解复杂船舶推进轴系扭转自由振动的固有频率以及固有频率所对应的复杂船舶推进轴系扭转自由振动的振型;复杂船舶推进轴系运转,在轴上施加扭矩;再计算F‑1(F(p,s));结合边界条件,利用复杂船舶推进轴系在附属机构位置处的强迫振动振型和轴承处的强迫振动振型表示轴上位置x=0处和x=L处强迫振动角速度;求解复杂船舶推进轴系扭转在附属机构位置处的强迫振动振型和轴承处的强迫振动振型;求解复杂船舶推进轴系扭转强迫振动的振型;结合留数定理求解出复杂船舶推进轴系扭转强迫振动响应。本发明分析船舶推进轴系扭振特性结果可靠度高、计算速度快。
Description
技术领域
本发明属于机械传动系统振动技术领域,涉及一种复杂船舶推进轴系扭转振动特性分析方法。
背景技术
船舶推进轴系是船舶动力装置的重要组成部分,包括从主机输出端推力轴承直到螺旋桨之间的传动轴以及轴上的附件。在船舶航行过程中,推进轴系受多种复杂载荷的作用,轴系扭转振动是影响舰船动力装置可靠性的重要因素之一。如果轴系长期工作于危险的临界转速区域内,则可能导致轴系动态性能的下降,产生变形及其疲劳失效,甚至发生扭断现象。这些都直接影响船舶航行性能和安全性,所以对轴系扭转振动计算的研究具有十分重要的意义。
在对复杂船舶推进轴系扭转特性分析的过程中,应从扭转自由振动和强迫振动两个方面入手,全面了解轴系扭转振动特性。传统的轴系扭转振动计算大多将轴系进行离散化,采用集中质量法来进行轴系建模,此种建模方法受模型离散化程度的影响,在进行轴系扭转振动固有特性分析过程中会产生误差,并且离散程度越高会导致计算速度显著降低。同样,在进行扭转强迫振动计算的过程中,普遍采用模态叠加法、纽马克法等数值解法,依然不能同时满足计算精度高、计算速度快的要求。因此,目前尚缺少一种高精度,快速求解复杂船舶推进轴系扭转振动特性的方法。
发明内容
本发明的目的在于提供一种精度更高、效率更优的复杂船舶推进轴系扭转振动特性分析方法。
本发明的目的是通过以下技术方案来实现的,包括以下步骤:
步骤一:测定复杂船舶推进轴系基本参数;
所述复杂船舶推进轴系基本参数包括:轴的半径和长度;附属机构相对轴的位置和质量;轴承相对轴的位置;
所述附属机构包括:齿轮箱、联轴器、曲柄连杆机构;
步骤二:结合边界条件求解复杂船舶推进轴系扭转自由振动的固有频率;
所涉及固有频率的计算表达式为
式中,s′为轴上附属机构个数,r为轴上轴承个数;Jj表示附属机构惯量的大小且j=1,2,3,...,s;Ki代表轴承的刚度且i=1,2,3,...,r;表示第j个附属机构的位置,表示第i个轴承的位置;Y(0)是轴上x=0处的角位移,Y′(0)是轴上x=0处的角速度;Y(L)是x=L处的角位移,Y′(L)是轴上x=L处的角速度;L为轴的长度;ρ为轴的密度,ω为固有频率,G为剪切模量,Ip为截面的极惯性矩;x为轴上位移;
所述复杂船舶推进轴系扭转自由振动的振型表达式为
步骤五:复杂船舶推进轴系运转,在轴上位置x=x′处施加扭矩M(x′,t),式中t表示施加力矩的时间;再计算F-1(F(p,s));
所涉及F-1(F(p,s))为F(p,s)的傅里叶逆变换,其中F(p,s)由M(x′,t)先进行拉氏变换再进行傅里叶变换得到;
步骤八:求解复杂船舶推进轴系扭转强迫振动的振型U(x,s);
所涉及复杂船舶推进轴系扭转强迫振动的振型U(x,s)的表达式为
式中,e为指数常数;s为复变量;符号"*"代表卷积;
步骤九:结合留数定理求解出复杂船舶推进轴系扭转强迫振动响应θ(x,t)。
本发明与现有技术相比的有益效果是:
1)本发明提供了一种高精度、快速计算复杂船舶推进轴系扭转振动特性分析方法,克服了常规算法的计算精度和快速运算不能同时兼容的缺点,同时也避免了使用有限元软件进行大规模建模、网格划分等繁琐的过程,大大地缩短了计算的时间;
2)本发明广泛适用于包含各类附属结构和支撑复杂工程实际轴系扭转振动计算;
3)本发明无需对大型轴系进行实际测试就能得到轴系的振动响应特征频谱,及多种影响因素对轴系振动响应的影响规律,大大节省了反复试验投入的人力和物力,并且本发明方法易于学习和操作。
附图说明
图1是带有多个集中惯量和弹性支承的轴系扭转强迫振动模型。
图2是本发明方法流程图。
图3是分布载荷作用下轴上不同位置点的轴系扭转强迫振动响应幅值图。
图4是分布载荷作用下确定位置的轴系扭转强迫振动响应的时间与幅值图。
图5为支撑刚度为107N/m条件下轴系扭转强迫振动响应图。
图6为支撑刚度为106N/m条件下轴系扭转强迫振动响应图。
图7为支撑刚度为105N/m条件下轴系扭转强迫振动响应图。
图8为支撑刚度为104N/m条件下轴系扭转强迫振动响应图。
具体实施方式
下面结合图1模型,对图2流程图的具体实施过程进行详细的说明。应该强调的是模型1仅为例子,而不是为了限制本发明的应用范围。
本发明的目的在于提出一种高精度、快速求解复杂船舶推进轴系扭转振动特性分析方法,弥补现有技术的不足。其具体实施过程包括以下步骤:
步骤一:复杂船舶推进轴系扭转振动特性分析前处理:
a)选择汽轮发电机组、燃气轮机推进轴系、内燃机轴系等大型复杂船舶推进轴系为计算对象,选定固有频率、振型、强迫振动计算或是影响因素分析计算。
b)计算复杂推进轴系附属结构的附加惯量。复杂推进轴系附属结构包括:齿轮箱、联轴器、曲柄连杆机构。
c)确定轴承刚度参数。
d)轴系所受到的外部扭矩如:内燃机汽缸压力产生的扭矩可以以简谐力的形式来表示。
e)确定轴系及轴承的参数:包括轴段的基本几何尺寸参数及材料属性参数、附加惯量大小及位置参数、不平衡质量大小及偏心距参数、轴承的刚度和位置参数以及不同形式扭矩的大小和作用位置参数。
步骤二:建立复杂船舶推进轴系扭转强迫振动微分方程:
基于哈密顿原理,建立了带有多个集中惯量和弹性支承的复杂推进轴系扭转振动微分方程。以图1模型为例,其中,ρ为轴的密度,R为轴的半径,G为剪切模量。θ=θ(x,t)为角位移。轴由s个集中惯量和r个弹性支撑组成,其中Jj(j=1,2,3···s)表示集中惯量的大小,Ki(i=1,2,3···,r)代表弹性支承的刚度。M(x,t)代表施加在轴上的单位长度的扭矩。应用哈密顿原理建模有:其中,T为动能,U为势能,W为非保守力作的功。具体表达形式如下:
动能的变分方程为:
势能的变分方程为:
非保守力做的功的变分方程为:
将上述变分方程带入哈密顿方程中,得到图1模型的动力学方程为:
步骤三:求解复杂船舶推进轴系扭转自由振动微分方程:
a)轴系扭转自由振动微分方程的求解过程:首先基于变量分离的方法将
时间项和位移项进行分离,假设θ(x,t)=Y(x)T(t),其中Y(x)仅是x的函数,T(t)仅是时间t的函数。将其带入到扭转自由振动微分方程中可得:
将时间项和位移项进行分离可得:
T″(t)+T(t)ω2=0
然后根据拉氏变换方法对单独的位移函数方程在拉氏域上进行求解可得:
最后根据拉氏逆变换求解出振型函数
其中,H(·)为Heaviside函数。
b)结合边界条件求解复杂船舶推进轴系扭转自由振动的固有频率;
所述的边界条件决定了复杂船舶推进轴系扭转自由振动中轴两端位置处的角位移和角速度值;在复杂船舶推进轴系扭转自由振动条件下,当轴系只有x=0受约束时,Y(0)=0;Y′(0)=1;Y(0)是轴上x=0处的角位移,Y′(0)是轴上x=0处的角速度;Y(L)是x=L处的角位移,Y′(L)是轴上x=L处的角速度。
所涉及固有频率的计算表达式为
式中,s′为轴上附属机构个数,r为轴上轴承个数;Jj表示附属机构惯量的大小且有j=1,2,3,...,s;Ki代表轴承的刚度且有i=1,2,3,...,r;表示第j个附属机构的位置,表示第i个轴承的位置;Y(0)是轴上x=0处的角位移,Y′(0)是轴上x=0处的角速度;Y(L)是x=L处的角位移,Y′(L)是轴上x=L处的角速度;ω为固有频率,G为剪切模量,Ip为截面的极惯性矩;x为轴上位移;
所述复杂船舶推进轴系扭转自由振动的振型表达式为
从而可以精确地求出轴系扭转振动的各阶固有频率对应的振型。
步骤四:求解复杂船舶推进轴系扭转强迫振动微分方程:
a)轴系扭转强迫振动微分方程的求解过程:首先将微分方程中的时间t在拉氏域上进行求解得到:
然后对上述方程进行傅里叶逆变换可以得到:
其中,F-1(F(p,s))代表F(p,s)的逆变换,符号"*"代表卷积。
b)复杂船舶推进轴系运转,在轴上位置x=x′处施加扭矩M(x′,t),式中t表示施加力矩的时间;计算F-1(F(p,s));
所涉及F-1(F(p,s))为F(p,s)的傅里叶逆变换,其中F(p,s)由M(x′,t)先进行拉氏变换再进行傅里叶变换得到;
e)求解复杂船舶推进轴系扭转强迫振动的振型U(x,s);
所涉及复杂船舶推进轴系扭转强迫振动的振型U(x,s)的表达式为
式中,e为指数常数;s为复变量;符号"*"代表卷积;
f)结合留数定理求解出复杂船舶推进轴系扭转强迫振动响应θ(x,t)。
步骤五:计算结果的后处理及输出:
a)固有频率和振型的计算结果输出:可输出任意阶数的固有频率及其相应的振型图;若输出多于两阶的固有频率,则固有频率将以列的形式来展示,振型图按照每行两个的排布原则,从上至下依次排布。此外,还可以输出集中惯量或弹性支承的大小或位置发生变化时,轴系固有频率和某阶振型的变化对比图,以集中惯量或弹性支承的大小或位置为横坐标,轴系固有频率或某阶振型为纵坐标。
b)强迫振动响应的计算结果输出:通过在不同的位置施加不同形式的激励力可以得到,轴系在受到分布载荷或是集中载荷作用下扭角随时间和位置变化的三维谱图,并且可以根据需要切换为扭角分别随时间和位置变化的平面图,便于直观的观察扭角随时间和位置的变化趋势,如图3和图4所示;并且可以输出随集中惯量、弹性支承以及激励力的大小或位置变换时的轴系扭角变化谱图;进一步的可以根据输出地扭角及轴段刚度得到轴系的应力分布图,图5~图8为不同支撑刚度条件下的轴系扭转强迫振动响应图。
Claims (4)
1.一种复杂船舶推进轴系扭转振动特性分析方法,其特征在于,包括以下步骤:
步骤一:测定复杂船舶推进轴系基本参数;
所述复杂船舶推进轴系基本参数包括:轴的半径和长度;附属机构相对轴的位置和质量;轴承相对轴的位置;
所述附属机构包括:齿轮箱、联轴器、曲柄连杆机构;
步骤二:结合边界条件求解复杂船舶推进轴系扭转自由振动的固有频率;
步骤五:复杂船舶推进轴系运转,在轴上位置x=x′处施加扭矩M(x′,t),式中t表示施加力矩的时间;再计算F-1(F(p,s));s为复变量;
所涉及F-1(F(p,s))为F(p,s)的傅里叶逆变换,其中F(p,s)由M(x′,t)先进行拉氏变换再进行傅里叶变换得到;
步骤八:求解复杂船舶推进轴系扭转强迫振动的振型U(x,s);
步骤九:结合留数定理求解出复杂船舶推进轴系扭转强迫振动响应θ(x,t)。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710563535.3A CN107292054B (zh) | 2017-07-12 | 2017-07-12 | 一种复杂船舶推进轴系扭转振动特性分析方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710563535.3A CN107292054B (zh) | 2017-07-12 | 2017-07-12 | 一种复杂船舶推进轴系扭转振动特性分析方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107292054A CN107292054A (zh) | 2017-10-24 |
CN107292054B true CN107292054B (zh) | 2020-11-20 |
Family
ID=60100582
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710563535.3A Active CN107292054B (zh) | 2017-07-12 | 2017-07-12 | 一种复杂船舶推进轴系扭转振动特性分析方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107292054B (zh) |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108827630B (zh) * | 2018-06-20 | 2019-11-05 | 武汉理工大学 | 船舶电力推进轴系扭转振动特性分析方法 |
CN112395709B (zh) * | 2020-11-30 | 2024-03-29 | 中船动力研究院有限公司 | 间接轴系扭振动态特性的修改方法、装置、设备和介质 |
CN112668094B (zh) * | 2020-12-22 | 2022-07-12 | 北京理工大学 | 一种双电机耦合驱动电动汽车动力传动系统扭振分析方法 |
CN112623143B (zh) * | 2020-12-25 | 2021-12-21 | 江苏集萃复合材料装备研究所有限公司 | 一种基于动力学的船舶动力轴系扭振疲劳参数拟合方法 |
CN113221244B (zh) * | 2021-05-17 | 2023-09-22 | 江苏集萃复合材料装备研究所有限公司 | 一种船舶螺旋桨转轴疲劳寿命的分析方法 |
CN113358308B (zh) * | 2021-06-03 | 2022-10-25 | 哈尔滨工业大学 | 基于有限测点和全局模态的组合结构横向位移确定方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101975633A (zh) * | 2010-09-16 | 2011-02-16 | 华南理工大学 | 应用连续细化分析傅立叶变换法的发动机激励力测量方法 |
CN102110186A (zh) * | 2009-12-25 | 2011-06-29 | 北京四方继保自动化股份有限公司 | 汽轮发电机组轴系扭振暂态响应实时计算方法 |
CN103177171A (zh) * | 2011-12-22 | 2013-06-26 | 中国舰船研究设计中心 | 基于解析法的船舶推进轴系纵向强迫振动计算方法 |
DE102011089895A1 (de) * | 2011-12-23 | 2013-06-27 | Fr. Lürssen Werft Gmbh & Co. Kg | Schiffsantriebswellenlager |
US9157741B2 (en) * | 2013-03-14 | 2015-10-13 | Raytheon Company | Simulator for simulating the operation of a fiber optic gyroscope |
CN105468046A (zh) * | 2016-01-08 | 2016-04-06 | 中国船舶重工集团公司第七0四研究所 | 船舶主推进系统多功率分支局部振动控制方法 |
-
2017
- 2017-07-12 CN CN201710563535.3A patent/CN107292054B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102110186A (zh) * | 2009-12-25 | 2011-06-29 | 北京四方继保自动化股份有限公司 | 汽轮发电机组轴系扭振暂态响应实时计算方法 |
CN101975633A (zh) * | 2010-09-16 | 2011-02-16 | 华南理工大学 | 应用连续细化分析傅立叶变换法的发动机激励力测量方法 |
CN103177171A (zh) * | 2011-12-22 | 2013-06-26 | 中国舰船研究设计中心 | 基于解析法的船舶推进轴系纵向强迫振动计算方法 |
DE102011089895A1 (de) * | 2011-12-23 | 2013-06-27 | Fr. Lürssen Werft Gmbh & Co. Kg | Schiffsantriebswellenlager |
US9157741B2 (en) * | 2013-03-14 | 2015-10-13 | Raytheon Company | Simulator for simulating the operation of a fiber optic gyroscope |
CN105468046A (zh) * | 2016-01-08 | 2016-04-06 | 中国船舶重工集团公司第七0四研究所 | 船舶主推进系统多功率分支局部振动控制方法 |
Non-Patent Citations (1)
Title |
---|
考虑齿轮啮合激励的齿轮传动轴系扭振特性分析;于姝雯等;《动力学与控制学报》;20161031;第14卷(第5期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN107292054A (zh) | 2017-10-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107292054B (zh) | 一种复杂船舶推进轴系扭转振动特性分析方法 | |
CN112765748B (zh) | 机理-数据异构信息融合的旋转机械数字孪生建模方法 | |
Chen | Vibration modelling and verifications for whole aero-engine | |
Wan et al. | An improved time-varying mesh stiffness algorithm and dynamic modeling of gear-rotor system with tooth root crack | |
Zou et al. | Coupled longitudinal-transverse dynamics of a marine propulsion shafting under superharmonic resonances | |
CN108776734B (zh) | 一种螺栓连接鼓筒转子结构的响应特性分析方法 | |
Huang et al. | Vibration analysis of marine propulsion shafting by the coupled finite element method | |
Xiang et al. | Experimental modal test of the spiral bevel gear wheel using the PolyMAX method | |
Su et al. | Analysis of dynamic characteristic for misalignment-spline gear shaft based on whole transfer matrix method | |
CN108534996B (zh) | 一种基于微振动瀑布图的飞轮刚度参数提取方法 | |
CN107563000B (zh) | 一种内燃机曲轴系统纵扭耦合刚度的计算方法 | |
CN109933916B (zh) | 基于轴系纵向振动响应测量反演螺旋桨纵向激励的方法及系统 | |
CN109388907B (zh) | 一种具有预设纵向振动动柔度的轴系的设计方法 | |
CN103177171A (zh) | 基于解析法的船舶推进轴系纵向强迫振动计算方法 | |
Bouslema et al. | Effects of modal truncation and condensation methods on the Frequency Response Function of a stage reducer connected by rigid coupling to a planetary gear system | |
CN113465961A (zh) | 风电机组传动链全尺寸地面加载测试系统及其控制方法 | |
HUANG et al. | Composite propeller's strain modal and structural vibration performance | |
Firouzi et al. | The effect of the shaft diameter and torsional stiffness on the whirling speed of the ship propeller shafting system | |
Liu et al. | Coupled vibration analysis of a bevel geared rotor-bearing system | |
Xu et al. | Research on vibration characteristics of propulsion shafting under unbalanced-misaligned coupling fault | |
Huang et al. | Analysis on propulsion shafting coupled torsional-longitudinal vibration under different applied loads | |
Mondragon et al. | Design, build, and test of a thrust test stand | |
Parikyan | Multi-cycle torsional vibration simulation with AVL excite designer | |
Tsyss et al. | Flexible coupling dynamic analysis as a ship shafting part | |
Kulesza et al. | Controlled deflection approach for rotor crack detection |
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 |