CN108388699B - 中心刚体-fgm楔形梁***末端动力学响应计算方法 - Google Patents
中心刚体-fgm楔形梁***末端动力学响应计算方法 Download PDFInfo
- Publication number
- CN108388699B CN108388699B CN201810082783.0A CN201810082783A CN108388699B CN 108388699 B CN108388699 B CN 108388699B CN 201810082783 A CN201810082783 A CN 201810082783A CN 108388699 B CN108388699 B CN 108388699B
- Authority
- CN
- China
- Prior art keywords
- fgm
- rigid body
- wedge
- central rigid
- beam system
- 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
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
-
- 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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/06—Power analysis or power optimisation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Geometry (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Evolutionary Computation (AREA)
- Computer Hardware Design (AREA)
- General Engineering & Computer Science (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Computational Mathematics (AREA)
- Rod-Shaped Construction Members (AREA)
- Buildings Adapted To Withstand Abnormal External Influences (AREA)
Abstract
本发明公开了一种中心刚体‑FGM楔形梁***末端动力学响应计算方法。运用弧长坐标来描述柔性FGM梁的几何位移关系,分别使用倾角和拉伸应变变量描述柔性梁的横向弯曲、纵向拉伸变形及剪切角;采用假设模态法离散变形场,运用第二类拉格朗日方程进行方程推导,得到中心刚体‑FGM楔形梁***的刚柔耦合动力学模型;运用C++编写FGM梁末端响应计算程序,通过读入中心刚体‑FGM楔形梁***几何参数、功能梯度参数、材料组成等参数,得出梁***末端响应随大范围转动时间的变化值。采用本发明方法,计算精度、效率较高。
Description
技术领域
本发明属于多体***动力学领域,具体地说,是一种中心刚体-FGM楔形梁***末端动力学响应计算方法。
背景技术
针对做大范围运动的中心刚体-柔性梁***的动力学响应问题,选择合适的坐标系建立精确的动力学模型,并在完整模型的基础上进行简化,在满足计算精度的情况下,得到计算效率较高的动力学模型,成为求解这类动力学问题的关键所在。
Librescu首次建立了中心刚体-FGM梁模型,并在此基础上对这一模型进行了振动分析。2005年,Librescu在以往建模方法的基础上,对圆柱薄壁梁在大范围运动下的动力学特性做出了研究。2012年,张伟将高阶剪切理论运用在建模过程中,考虑到离心力的作用,建立了旋转运动下FGM板的动力学方程。黎亮首次提出了倾角坐标,基于细长梁假设,对FGM梁***在大范围运动下的动力学问题进行了研究。目前已有的工作主要是运用传统的直角坐标系计算柔性梁的动力学问题,计算效率较为低下,选用合理的建模方法,建立满足计算精度,计算效率较高的模型,成为这类问题研究中的重点。
发明内容
本发明的目的是,针对大范围旋转运动下中心刚体-FGM楔形梁***末端动力学响应问题,提供一种数值仿真的计算方法,将FGM楔形梁几何参数、功能梯度参数、材料组成分别进行设置,得到FGM楔形梁的末端横向变形与轴向变形。
实现本发明目的地技术解决方案为:一种中心刚体-FGM楔形梁***末端动力学响应计算方法,包括以下步骤:
(1)设定中心刚体-FGM楔形梁***相关参数:中心刚体转动惯量、楔形梁几何尺寸、FGM梁组成材料组成、功能梯度指数,并给出大范围运动角速度规律;
(2)选用弧长坐标中心刚体-FGM楔形梁***进行建模,运用几何关系描述中心刚体-FGM楔形梁***的变形场,得出柔性梁末端位移表达式;
(3)取中心刚体-FGM楔形梁***的一段微元进行分析,写出柔性梁***在大范围转动下的动能和势能表达式;
(4)运用假设模态法对每段微元的横向弯曲角、纵向拉伸量及剪切角进行离散,并将动能与势能带入第二类Lagrange方程,并将方程中二次以上项舍去,得到中心刚体-柔性梁***的刚柔耦合动力学方程;
(5)针对中心刚体-FGM楔形梁***,运用梁高比Rh,梁宽比Rb描述楔形梁几何形状;运用梁悬臂端与自由端材料参数及功能梯度参数描述FGM梁材料组成;
(6)根据步骤(4)中动力学方程和步骤(5)给定的参数,得出FGM楔形梁末端横向变形及轴向变形随时间变化规律数据。
步骤(1)中,大范围运动角速度规律为:
式中,ω为转动角速度,ω0为初始转动角速度,T为大范围转动计算时长;
步骤(2)中,柔性梁末端位移表达式为:
式中,u(t)为柔性梁末端轴向变形,v(t)为柔性梁末端横向变形,α(s,t)为弧长坐标s处横截面的弯曲角度,ε(s,t)为弧长坐标s处轴向拉伸量,l为柔性梁长度。
步骤(3)中,柔性梁***的动能表达式为:
式中,Joh为中心刚体转动惯量,θ0为中心刚体角位移,ρ(s)为柔性梁沿轴向密度函数,A(s)为柔性梁沿轴向横截面积函数,x0、y0为梁轴线上一点处坐标分量,γ(s,t)为弧长坐标s处横截面的剪切角。
柔性梁势能表达式为:
式中,E(s)为柔性梁沿轴向弹性模量函数,G(s)为柔性梁沿轴向剪切模量函数,k为剪切修正系数。
步骤(4)中,采用假设模态法对柔性梁的变形进行描述,将倾角、纵向拉伸量、剪切角进行离散处理:
其中,φi(s)为一端固支一端自由杆试函数行向量,A(t)、B(t)、C(t)为与时间相关项列向量。将上式带入第二类Lagrange方程,并舍去部分高阶项,得到非惯性系下中心刚体-柔性梁***的动力学方程:
式中各项分别为:
步骤(5)中梁高比Rh,梁宽比Rb的取值范围分别为0≤Rh≤1,0≤Rb≤1,材料参数分别需设定梁固定端及自由端材料密度与弹性模量。
步骤(6)中中心刚体-FGM楔形梁***需设置的参数分别为:柔性梁长度、悬臂端截面积及惯性矩、梁高比及梁宽比、中心刚体转动惯量、悬臂端及固定端材料密度与弹性模量、功能梯度指数。
本发明与现有技术相比,其显著优点为:
(1)在建模过程中,采用倾角坐标,可通过描述梁中轴线上点的位置,对柔性梁的变形问题进行描述,在建模过程中较传统的坐标系更加方便简洁;
(2)在公式推导中,基于Timoshenko梁假设,考虑到了横向弯曲、轴向拉伸及剪切效应;计算精度较高;
(3)在计算末端响应时,在满足精度的前提下,省略了部分高阶项,计算效率较高;
(4)FGM楔形梁参数设置时,可对中心刚体转动惯量、几何尺寸、材料组成、功能梯度参数等进行设置,可运用于多种形式的FGM楔形梁。
附图说明
图1是柔性梁变形示意图。
图2是FGM楔形梁几何尺寸示意图。
图3是文件“0.txt”数据。
图4是C++程序运行过程。
图5是计算完成输出文件“v.txt”及“u.txt”。
图6是本发明方法的实现流程图。
具体实施方式
下面结合实施例及附图对本发明进行进一步介绍。
实例:一种中心刚体-FGM楔形梁***末端动力学响应计算方法,如图6所示,该方法包括以下步骤:
(1)设定旋转速度规律:
(2)设定中心刚体-FGM梁***的相关参数,悬臂端材料为铝,自由端材料为陶瓷。具体数值在表1中给出。
表1 中心刚体-FGM梁***的几何参数
表2 中心刚体-FGM梁***的材料参数
(3)在文件“0.txt”中输入数据,如图3所示;
(4)编写C++程序求解动力学方程(4),运行程序并读入文件“0.txt”中数据,计算过程如图4所示;
(5)程序运行完成后,输出FGM楔形梁***在大范围运动下横向变形与轴向变形随时间变化规律文件“v.txt”及“u.txt”(图5),文件中数据可用于生成曲线来研究FGM楔形梁末端变形随时间变化规律。
Claims (6)
1.一种中心刚体-FGM楔形梁***末端动力学响应计算方法,其特征在于包括以下步骤:
(1)设定中心刚体-FGM楔形梁***相关参数:中心刚体转动惯量、楔形梁几何尺寸、FGM梁组成材料组成、功能梯度指数,并给出大范围运动角速度规律;
(2)选用弧长坐标对中心刚体-FGM楔形梁***进行建模,运用几何关系描述中心刚体-FGM楔形梁***的变形场,得出梁末端位移表达式;
取中心刚体-FGM楔形梁***的一段微元进行分析,得出柔性梁***在大范围转动下的动能和势能表达式;柔性梁***的动能表达式为:
式中,Joh为中心刚体转动惯量,θ0为中心刚体角位移,ρ(s)为柔性梁沿轴向密度函数,A(s)为柔性梁沿轴向横截面积函数,x0、y0为梁轴线上一点处坐标分量,γ(s,t)为弧长坐标s处横截面的剪切角;
柔性梁***的势能表达式为:
(3)式中,E(s)为柔性梁沿轴向弹性模量函数,G(s)为柔性梁沿轴向剪切模量函数,k为剪切修正系数;
(4)运用假设模态法对每段微元的横向弯曲角、纵向拉伸量及剪切角进行离散,并将动能与势能带入第二类Lagrange方程,并将方程中二次以上项舍去,得到中心刚体-柔性梁***的刚柔耦合动力学方程;
(5)针对中心刚体-FGM楔形梁***,运用梁高比Rh,梁宽比Rb描述楔形梁几何形状;运用梁悬臂端与自由端材料参数及功能梯度参数描述FGM梁材料组成;
(6)根据步骤(4)中动力学方程和步骤(5)给定的参数,得出FGM楔形梁末端横向变形及轴向变形随时间变化规律数据。
5.根据权利要求1所述的中心刚体-FGM楔形梁***末端动力学响应计算方法,其特征在于:步骤(5)中梁高比Rh、梁宽比Rb的取值范围分别为0≤Rh≤1,0≤Rb≤1,材料参数分别需设定梁固定端及自由端材料密度与弹性模量。
6.根据权利要求1所述的中心刚体-FGM楔形梁***末端动力学响应计算方法,其特征在于:步骤(6)中中心刚体-FGM楔形梁***需设置的参数分别为:柔性梁长度、悬臂端截面积及惯性矩、梁高比及梁宽比、中心刚体转动惯量、悬臂端及固定端材料密度与弹性模量、功能梯度指数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810082783.0A CN108388699B (zh) | 2018-01-29 | 2018-01-29 | 中心刚体-fgm楔形梁***末端动力学响应计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810082783.0A CN108388699B (zh) | 2018-01-29 | 2018-01-29 | 中心刚体-fgm楔形梁***末端动力学响应计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108388699A CN108388699A (zh) | 2018-08-10 |
CN108388699B true CN108388699B (zh) | 2022-02-18 |
Family
ID=63074172
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810082783.0A Active CN108388699B (zh) | 2018-01-29 | 2018-01-29 | 中心刚体-fgm楔形梁***末端动力学响应计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108388699B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109940613B (zh) * | 2019-03-08 | 2022-04-01 | 南京理工大学 | 一种计算含压电材料机械臂动力学响应及控制的仿真方法 |
CN110008543B (zh) * | 2019-03-21 | 2022-09-13 | 南京理工大学 | 一种考虑梁中性轴对旋转梁动力学响应影响的仿真方法 |
CN113312775A (zh) * | 2021-06-01 | 2021-08-27 | 扬州大学 | 变温场中fgm梁的动力学仿真模型、其建立方法及仿真方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104808512A (zh) * | 2015-03-03 | 2015-07-29 | 北京空间飞行器总体设计部 | 一种航天器多级驱动刚柔耦合响应的获取方法 |
CN106777918A (zh) * | 2016-11-29 | 2017-05-31 | 河南理工大学 | 基于导波和模糊算法的功能梯度结构材料特性的反演方法 |
-
2018
- 2018-01-29 CN CN201810082783.0A patent/CN108388699B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104808512A (zh) * | 2015-03-03 | 2015-07-29 | 北京空间飞行器总体设计部 | 一种航天器多级驱动刚柔耦合响应的获取方法 |
CN106777918A (zh) * | 2016-11-29 | 2017-05-31 | 河南理工大学 | 基于导波和模糊算法的功能梯度结构材料特性的反演方法 |
Non-Patent Citations (4)
Title |
---|
Dynamic analysis of rotating axially FG tapered beams based on a new rigid-flexible coupled dynamic model using the B-spline method;Li L et al.;《Composite Structures》;20151230;第124卷;第357-367页 * |
中心刚体-功能梯度材料梁***的动力学特性;黎亮 等;《机械工程学报》;20130715;第49卷(第13期);第77-84页 * |
刚—柔耦合复合结构的动力学建模理论研究;黎亮;《中国博士学位论文全文数据库工程科技II辑》;20160415;第C028-4页 * |
旋转叶片刚柔耦合***动力学分析;韩广才 等;《哈尔滨工程大学学报》;20110615;第32卷(第6期);第736-741页 * |
Also Published As
Publication number | Publication date |
---|---|
CN108388699A (zh) | 2018-08-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108388699B (zh) | 中心刚体-fgm楔形梁***末端动力学响应计算方法 | |
CN109940613B (zh) | 一种计算含压电材料机械臂动力学响应及控制的仿真方法 | |
CN112560167B (zh) | 机翼结构力学高保真降阶仿真方法、电子设备、存储介质 | |
CN106407582B (zh) | 一种用于柔性机构静力学和动力学分析的高效计算方法 | |
CN107976908B (zh) | 一种飞行器耦合动稳定性特征分析方法 | |
CN109271655B (zh) | 一种基于非对称有限元算法的材料尺度效应分析方法 | |
Calhoun et al. | Nonlinear finite element analysis of clamped arches | |
CN108182330A (zh) | 一种基于b样条计算柔性矩形薄板刚柔耦合动力学响应的方法 | |
Ni et al. | Analysis of parasitic motion in parallelogram compliant mechanism | |
CN104850683B (zh) | 基于弱形式求积元法计算材料裂纹尖端应力场系数的方法 | |
CN107220678B (zh) | 多自由度梁式结构非线性类型确定方法 | |
CN114091153B (zh) | 考虑接缝界面动力效应的新型接缝模型及其构建方法和数值实现方法 | |
Raju et al. | Nonlinear structural analysis using integrated force method | |
CN110008543B (zh) | 一种考虑梁中性轴对旋转梁动力学响应影响的仿真方法 | |
CN109684766B (zh) | 含转角的大变形柔性梁单元建模方法 | |
CN108334688B (zh) | 基于matlab的旋转功能梯度厚板动力学行为计算方法 | |
CN208707550U (zh) | 一种精密缩小柔性微动*** | |
CN109977498A (zh) | 一种基于hoc计算fgm梯形梁动力学响应的方法 | |
Zheng et al. | An Optimized 3D Probe Using Sensitivity and Compliance Analysis | |
Kosmatka | The use of cross-section warping functions in composite rotor blade analysis | |
CN110879553B (zh) | 一种基于输出状态可用的微陀螺仪的控制方法及*** | |
CN117171915A (zh) | 一种基于显示动力学理论的梁段吊装变形快速计算方法 | |
Huang et al. | Design and Optimization of Amplifying Mechanism in Gas-Liquid Thermoelectric Power Device | |
Xu et al. | Static and dynamic analysis of flexure hinge mechanisms using the weak-form quadrature element method | |
CN114282303A (zh) | 一种基于高斯积分的旋翼桨叶结构有限元建模方法 |
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 |