CN111310269B - 一种考虑边界位移的固支梁结构非线性动特性分析方法 - Google Patents
一种考虑边界位移的固支梁结构非线性动特性分析方法 Download PDFInfo
- Publication number
- CN111310269B CN111310269B CN202010114613.3A CN202010114613A CN111310269B CN 111310269 B CN111310269 B CN 111310269B CN 202010114613 A CN202010114613 A CN 202010114613A CN 111310269 B CN111310269 B CN 111310269B
- Authority
- CN
- China
- Prior art keywords
- displacement
- vibration
- formula
- section
- equation
- 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
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16C—COMPUTATIONAL CHEMISTRY; CHEMOINFORMATICS; COMPUTATIONAL MATERIALS SCIENCE
- G16C60/00—Computational materials science, i.e. ICT specially adapted for investigating the physical or chemical properties of materials or phenomena associated with their design, synthesis, processing, characterisation or utilisation
Landscapes
- Engineering & Computer Science (AREA)
- Computing Systems (AREA)
- Theoretical Computer Science (AREA)
- Life Sciences & Earth Sciences (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Bioinformatics & Computational Biology (AREA)
- Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
Abstract
本发明公开了一种考虑边界位移的固支梁结构非线性动特性分析方法,包括如下步骤:根据边界位移和结构材料参数,通过非线性分析计算固支梁的挠度曲线的解析解;建立考虑边界位移的固支梁结构非线性振动方程,通过矩阵传递法和逐步搜索法求解变系数振动微分方程,进而得到考虑边界位移的梁结构非线性振动模态频率和模态振型解析解。本发明的动特性分析方法考虑了边界位移导致的几何非线性梁的动特性的影响,并考虑了梁挠度和振动模态的耦合作用。本发明能够有效提高复杂环境下梁结构动特性分析精度,指导工程结构设计。
Description
技术领域
本发明涉及一种梁结构动特性分析方法,尤其涉及一种考虑边界位移的固支梁结构非线性动特性分析方法。
背景技术
梁结构被广泛运用于工程结构中,并在服役过程中面临复杂的载荷环境。结构特性分析能够为结构设计提供了重要指导。边界位移是梁结构在工程实际中的一种典型工况。目前针对梁结构非线性特性的研究,主要通过线性叠加的形式,借助有限元软件进行分析,未见考虑静态挠度和振动位移耦合作用的分析方法。同时,分析方法相比于有限元方法具有更高的分析精度。
因此,为了更好的指导结构设计和发挥结构性能,亟待发展一种考虑边界位移的固支梁结构非线性动特性分析方法。
发明内容
发明目的:本发明的目的是提供一种有效提高复杂环境下梁结构动特性分析精度的考虑边界位移的非线性动特性分析方法。
技术方案:为实现以上目的,本发明公开了一种考虑边界位移的固支梁结构非线性动特性分析方法,包括如下步骤:
(1)、根据初始边界位移和结构材料参数,基于假定函数,计算固支梁结构大变形的挠度曲线解析解;
(2)、考虑边界位移下挠度和振动位移间的耦合作用,建立边界位移下的固支梁结构横向运动的动力学方程;根据结构动力学方程建立变系数振动微分方程,再通过矩阵传递法和逐步搜索法得到考虑边界位移的梁结构非线性动特性分析结果。
其中,所述步骤(1)中计算固支梁结构大变形的挠度曲线解析解的具体步骤为:
首先根据固支梁边界约束条件,假定挠度曲线位移函数w0,如
式中n为根据计算精度需求选取的表达式项数,m表示第m阶表达式,am为第m阶表达式的未知系数,x为梁轴向坐标,L为梁长度;计算结构应变能Up0和外力功Wp0,计算公式为:
Wp0=Fsw0|x=0
式中E为弹性模量,A为梁截面面积,I为转动惯量,上标“'”表示对x的一阶偏导,“″”表示对x二阶偏导,Fs为约束端未知剪力;
然后计算梁结构总势能Π0=Up0-Wp0,根据变分原理建立求解方程组
与边界约束方程联立,求解挠度曲线位移函数的系数am,进而得到考虑边界位移的结构挠度曲线解析解。
优选地,所述步骤(2)中梁结构非线性动特性分析的具体步骤为:
首先计算振动下结构应变能Up、外力功Wp和动能Tk:
然后根据哈密顿原理,建立边界位移下的固支梁结构横向运动的动力学方程,
式中上标“(4)”表示对x四阶偏导;
固支梁自由振动方程为:
再由分离变量法获取变系数振动微分方程:
式中φ为梁结构模态振型,ω为模态频率;
再采用传递矩阵法求解固支梁模态振型与频率,将固支梁分为ni段,第i段梁的振动微分方程为:
式中φi为第i段梁结构振型;定义第i段梁的系数δi和β为
式中,xi和xi+1为第i段梁的两端坐标,得到第i段梁的方程解为:
φi(x)=Ai sinβi1(x-xi)+Bi cosβi1(x-xi)+Ci coshβi2(x-xi)+Di sinhβi2(x-xi)
式中,xi≤x≤xi+1,(i=1,2,3,…,ni),Ai、Bi、Ci、Di是第i段梁解的待定系数,第i段梁中βi1和βi2表达式为:
由第i段梁和第i+1段梁在xi点处位移、转角、弯矩和剪力连续,有以下关系:
φi(xi+1)=φi+1(xi+1)
φ′i(xi+1)=φ′i+1(xi+1)
EIφ″i(xi+1)=EIφ″i+1(xi+1)
式中上标“(3)”表示对x三阶偏导,令
ψi=[Ai Bi Ci Di]T
θi(x)=[sinβi1(x-xi) cosβi1(x-xi) coshβi2(x-xi) sinhβi2(x-xi)]
Φi(xj)=[θi(xj) θ′i(xj) θ″i(xj) θ″′i(xj)]T
上标“T”表示矩阵转置,则有:
上标“-1”表示矩阵求逆,令
进而,边界约束条件表示为:
式中
最后通过逐步搜索法求解下述方程,得到结构模态频率,进而得到每段梁解的待定系数和梁的模态振型:
再者,所述固支梁为铝合金梁。
进一步,所述固支梁的横截面为矩形。
有益效果:与现有技术相比,本发明具有以下显著优点:本发明的动特性分析方法考虑了边界位移导致的几何非线性梁的动特性的影响,并考虑了梁挠度和振动模态的耦合作用,本发明能够有效提高复杂环境下梁结构动特性分析精度,指导工程结构设计。
附图说明
图1为本发明的流程图;
图2为本发明中固支梁及边界位移示意图;
图3为本发明中梁结构第1阶振型图;
图4为本发明中梁结构第2阶振型图;
图5为本发明中梁结构第3阶振型图;
图6为本发明中梁结构第4阶振型图;
图7为本发明中梁结构第5阶振型图;
图8为本发明中梁结构第6阶振型图。
具体实施方式
下面结合附图对本发明的技术方案作进一步说明。
如图1所示,本发明提供一种考虑边界位移的固支梁结构非线性动特性分析方法,包括如下步骤:
(1)、根据初始边界位移和结构材料参数,基于假定函数,计算固支梁结构大变形的挠度曲线解析解;
其中计算固支梁结构大变形的挠度曲线解析解的具体步骤为:
首先根据固支梁边界约束条件,假定挠度曲线位移函数w0,如
式中n为根据计算精度需求选取的表达式项数,m表示第m阶表达式,am为第m阶表达式的未知系数,x为梁轴向坐标,L为梁长度;计算结构应变能Up0和外力功Wp0,计算公式为:
Wp0=Fsw0|x=0
式中E为弹性模量,A为梁截面面积,I为转动惯量,上标“'”表示对x的一阶偏导,“″”表示对x二阶偏导,Fs为约束端未知剪力;
然后计算梁结构总势能Π0=Up0-Wp0,根据变分原理建立求解方程组
与边界约束方程联立,求解挠度曲线位移函数的系数am,进而得到考虑边界位移的结构挠度曲线解析解;
(2)、考虑边界位移下挠度和振动位移间的耦合作用,建立边界位移下的固支梁结构横向运动的动力学方程;根据结构动力学方程建立变系数振动微分方程,再通过矩阵传递法和逐步搜索法得到考虑边界位移的梁结构非线性动特性分析结果;
其中梁结构非线性动特性分析的具体步骤为:
首先计算振动下结构应变能Up、外力功Wp和动能Tk:
然后根据哈密顿原理,建立边界位移下的固支梁结构横向运动的动力学方程,
式中上标“(4)”表示对x四阶偏导;
固支梁自由振动方程为:
再由分离变量法获取变系数振动微分方程:
式中φ为梁结构模态振型,ω为模态频率;
再采用传递矩阵法求解固支梁模态振型与频率,将固支梁分为ni段,第i段梁的振动微分方程为:
式中φi为第i段梁结构振型;定义第i段梁的系数δi和β为
式中,xi和xi+1为第i段梁的两端坐标,得到第i段梁的方程解为:
φi(x)=Ai sinβi1(x-xi)+Bi cosβi1(x-xi)+Ci coshβi2(x-xi)+Di sinhβi2(x-xi)
式中,xi≤x≤xi+1,(i=1,2,3,…,ni),Ai、Bi、Ci、Di是第i段梁解的待定系数,第i段梁中βi1和βi2表达式为:
由第i段梁和第i+1段梁在xi点处位移、转角、弯矩和剪力连续,有以下关系:
φi(xi+1)=φi+1(xi+1)
φ′i(xi+1)=φ′i+1(xi+1)
EIφ″i(xi+1)=EIφ″i+1(xi+1)
式中上标“(3)”表示对x三阶偏导,令
ψi=[Ai Bi Ci Di]T
θi(x)=[sinβi1(x-xi) cosβi1(x-xi) coshβi2(x-xi) sinhβi2(x-xi)]
Φi(xj)=[θi(xj) θ′i(xj) θ″i (xj)θi″′(xj)]T
上标“T”表示矩阵转置,则有:
上标“-1”表示矩阵求逆,令
进而,边界约束条件表示为:
式中
最后通过逐步搜索法求解下述方程,得到结构模态频率,进而得到每段梁解的待定系数和梁的模态振型:
实施例1
以两端固支梁在边界强制位移下的静力大变形分析为例,如图2所示。采用铝合金梁,弹性模量E=70000MPa,泊松比μ=0.3,密度ρ=2.7×10-9t/mm3,梁长度L=320mm,梁为矩形截面,宽b=20mm,高h=2mm,左端强制位移-10mm。
首先假定挠度曲线位移函数w0,以梁左端为坐标原点,边界约束条件可表示为w0|x=0=-10
w0|x=L=0
w′0|x=0=w′0|x=L=0
式中上标“'”表示对x的一阶偏导;根据边界约束条件,w0可假定为
式中n为根据计算精度需求选取的表达式项数,此实施例中n取为5时,挠度曲线位移函数w0有足够计算精度。m表示第m阶表达式,am为第m阶表达式的未知系数,x为梁轴向坐标,L为梁长度;计算结构应变能Up0和外力功Wp0,计算公式为:
Wp0=Fsw0|x=0
式中E为弹性模量,A为梁截面面积,I为转动惯量,上标“″”表示对x二阶偏导,Fs为约束端未知剪力;
然后计算梁结构总势能Π0=Up0-Wp0,根据变分原理建立求解方程组
联立求解方程组及约束方程,可计算得挠度曲线位移函数w0为
进一步,将挠度曲线位移函数w0代入下述变系数振动微分方程:
式中φ为梁结构模态振型,ω为模态频率;
再采用传递矩阵法求解固支梁模态振型与频率,将固支梁分为ni段,第i段梁的振动微分方程为:
式中φi为第i段梁结构振型;定义第i段梁的系数δi和β为
式中,xi和xi+1为第i段梁的两端坐标,得到第i段梁的方程解为:
φi(x)=Ai sinβi1(x-xi)+Bi cosβi1(x-xi)+Ci coshβi2(x-xi)+Di sinhβi2(x-xi)
式中,xi≤x≤xi+1,(i=1,2,3,…,ni),Ai、Bi、Ci、Di是第i段梁解的待定系数,第i段梁中βi1和βi2表达式为:
由第i段梁和第i+1段梁在xi点处位移、转角、弯矩和剪力连续,有以下关系:
φi(xi+1)=φi+1(xi+1)
φ′i(xi+1)=φ′i+1(xi+1)
EIφ″i(xi+1)=EIφ″i+1(xi+1)
式中上标“(3)”表示对x三阶偏导,令
ψi=[Ai Bi Ci Di]T
θi(x)=[sinβi1(x-xi) cosβi1(x-xi) coshβi2(x-xi) sinhβi2(x-xi)]
Φi(xj)=[θi(xj) θ′i(xj) θ″i(xj) θ″′i(xj)]T
上标“T”表示矩阵转置,则有:
上标“-1”表示矩阵求逆,令
进而,边界约束条件表示为:
式中
最后通过逐步搜索法求解下述方程,得到结构模态频率,进而得到每段梁解的待定系数和梁的模态振型:
通过方程求解,可得梁的前6阶模态频率如表1所示,前6阶结构振型如图3~图8所示。
表1前6阶模态频率
本发明的动特性分析方法考虑了边界位移导致的几何非线性梁的动特性的影响,并考虑了梁挠度和振动模态的耦合作用。本发明能够有效提高复杂环境下梁结构动特性分析精度,指导工程结构设计。
Claims (4)
1.一种考虑边界位移的固支梁结构非线性动特性分析方法,其特征在于,包括如下步骤:
(1)、根据初始边界位移和结构材料参数,基于假定函数,计算固支梁结构大变形的挠度曲线解析解;其中计算固支梁结构大变形的挠度曲线解析解的具体步骤为:
首先根据固支梁边界约束条件,假定挠度曲线位移函数w0,如
式中n为根据计算精度需求选取的表达式项数,m表示第m阶表达式,am为第m阶表达式的未知系数,x为梁轴向坐标,L为梁长度;计算结构应变能Up0和外力功Wp0,计算公式为:
Wp0=Fsw0|x=0
式中E为弹性模量,A为梁截面面积,I为转动惯量,上标“'”表示对x的一阶偏导,“"”表示对x二阶偏导,Fs为约束端未知剪力;
然后计算梁结构总势能Π0=Up0-Wp0,根据变分原理建立求解方程组
与边界约束方程联立,求解挠度曲线位移函数的系数am,进而得到考虑边界位移的结构挠度曲线解析解;
(2)、考虑边界位移下挠度和振动位移间的耦合作用,建立边界位移下的固支梁结构横向运动的动力学方程;根据结构动力学方程建立变系数振动微分方程,再通过矩阵传递法和逐步搜索法得到考虑边界位移的梁结构非线性动特性分析结果。
2.根据权利要求1所述的一种考虑边界位移的固支梁结构非线性动特性分析方法,其特征在于:所述步骤(2)中梁结构非线性动特性分析的具体步骤为:
首先计算振动下结构应变能Up、外力功Wp和动能Tk:
然后根据哈密顿原理,建立边界位移下的固支梁结构横向运动的动力学方程,
式中上标“(4)”表示对x四阶偏导;
固支梁自由振动方程为:
再由分离变量法获取变系数振动微分方程:
式中φ为梁结构模态振型,ω为模态频率;
再采用传递矩阵法求解固支梁模态振型与频率,将固支梁分为ni段,第i段梁的振动微分方程为:
式中φi为第i段梁结构振型;定义第i段梁的系数δi和β为
式中,xi和xi+1为第i段梁的两端坐标,得到第i段梁的方程解为:
φi(x)=Ai sinβi1(x-xi)+Bi cosβi1(x-xi)+Ci coshβi2(x-xi)+Di sinhβi2(x-xi)
式中,xi≤x≤xi+1,i=1,2,3,…,ni,Ai、Bi、Ci、Di是第i段梁解的待定系数,第i段梁中βi1和βi2表达式为:
由第i段梁和第i+1段梁在xi点处位移、转角、弯矩和剪力连续,有以下关系:
φi(xi+1)=φi+1(xi+1)
φ′i(xi+1)=φ′i+1(xi+1)
EIφ″i(xi+1)=EIφ″i+1(xi+1)
式中上标“(3)”表示对x三阶偏导,令
ψi=[Ai Bi Ci Di]T
θi(x)=[sinβi1(x-xi) cosβi1(x-xi) coshβi2(x-xi) sinhβi2(x-xi)]
Φi(xj)=[θi(xj) θ′i(xj) θ″i(xj) θ″′i(xj)]T
上标“T”表示矩阵转置,则有:
上标“-1”表示矩阵求逆,令
进而,边界约束条件表示为:
式中
最后通过逐步搜索法求解下述方程,得到结构模态频率,进而得到每段梁解的待定系数和梁的模态振型:
3.根据权利要求1所述的一种考虑边界位移的固支梁结构非线性动特性分析方法,其特征在于:所述固支梁为铝合金梁。
4.根据权利要求1所述的一种考虑边界位移的固支梁结构非线性动特性分析方法,其特征在于:所述固支梁的横截面为矩形。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010114613.3A CN111310269B (zh) | 2020-02-25 | 2020-02-25 | 一种考虑边界位移的固支梁结构非线性动特性分析方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010114613.3A CN111310269B (zh) | 2020-02-25 | 2020-02-25 | 一种考虑边界位移的固支梁结构非线性动特性分析方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111310269A CN111310269A (zh) | 2020-06-19 |
CN111310269B true CN111310269B (zh) | 2021-01-26 |
Family
ID=71160180
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010114613.3A Active CN111310269B (zh) | 2020-02-25 | 2020-02-25 | 一种考虑边界位移的固支梁结构非线性动特性分析方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111310269B (zh) |
Families Citing this family (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111783198B (zh) * | 2020-06-21 | 2022-08-12 | 西北工业大学 | 一种基于双梁模型的悬索桥精细化动力分析方法 |
CN112287488B (zh) * | 2020-11-06 | 2023-03-28 | 清华大学 | 一种用于微机电器件的非线性弹性约束结构及约束方法 |
CN112989486B (zh) * | 2021-03-10 | 2023-03-24 | 北京航空航天大学 | 一种构建起落架模型的方法和装置 |
CN113656984B (zh) * | 2021-08-31 | 2023-05-02 | 北京建筑大学 | 一种基于蒙特卡洛法的梁结构固有频率和振型计算方法 |
CN114065429B (zh) * | 2021-11-18 | 2023-04-18 | 哈尔滨工业大学 | 一种单对称截面风力机叶片固有特性求解方法 |
CN116844664B (zh) * | 2023-03-10 | 2023-11-17 | 大连理工大学 | 考虑不确定性的数据驱动的动力学序列线性规划方法 |
CN117271949B (zh) * | 2023-11-23 | 2024-03-12 | 华东交通大学 | 考虑弹性边界和大垂度影响的悬索自振分析方法及系统 |
CN118378383B (zh) * | 2024-04-26 | 2024-10-01 | 西安理工大学 | 框架结构参数振动的分析方法 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101114308A (zh) * | 2007-08-24 | 2008-01-30 | 西安交通大学 | 一种三类变量区间b样条小波梁单元构造方法 |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102012953B (zh) * | 2010-11-04 | 2013-05-08 | 西北工业大学 | Cfd/csd耦合求解非线性气动弹性仿真方法 |
CN104359754B (zh) * | 2014-10-10 | 2016-11-30 | 中国科学院工程热物理研究所 | 一种梁构件弯扭组合加载试验装置及方法 |
CN104778377B (zh) * | 2015-05-04 | 2017-11-10 | 中国矿业大学 | 一种组合梁弯曲振动的固有频率分析方法 |
US10365254B2 (en) * | 2015-10-14 | 2019-07-30 | The Regents Of The University Of California | Assessment of blood coagulation using an acoustic radiation force based optical coherence elastography (ARF-OCE) |
-
2020
- 2020-02-25 CN CN202010114613.3A patent/CN111310269B/zh active Active
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101114308A (zh) * | 2007-08-24 | 2008-01-30 | 西安交通大学 | 一种三类变量区间b样条小波梁单元构造方法 |
Non-Patent Citations (1)
Title |
---|
An efficient transient analysis method for time-varying structures based on statistical energy analysis;yang xuan等;《Mechanics Research Communications》;20181231;第93-99页 * |
Also Published As
Publication number | Publication date |
---|---|
CN111310269A (zh) | 2020-06-19 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111310269B (zh) | 一种考虑边界位移的固支梁结构非线性动特性分析方法 | |
Chen et al. | Generalized equations for estimating stress concentration factors of various notch flexure hinges | |
Kwon et al. | A nonlinear elastic spline finite strip analysis for thin-walled sections | |
Li et al. | A tapered Timoshenko–Euler beam element for analysis of steel portal frames | |
JP2019056220A (ja) | 床構造に用いる鉄骨梁の設計方法、床構造 | |
Banerjee | Explicit frequency equation and mode shapes of a cantilever beam coupled in bending and torsion | |
Rasmussen et al. | Beam-element-based analysis of locally and/or distortionally buckled members: Theory | |
Rahai et al. | Buckling analysis of non-prismatic columns based on modified vibration modes | |
Mahar et al. | A computational study on buckling behavior of cold-formed steel built-up columns using compound spline finite strip method | |
Bui | Semi-analytical finite strip method based on the shallow shell theory in buckling analysis of cold-formed sections | |
Wong et al. | Kriging-based timoshenko beam elements with the discrete shear gap technique | |
JP7453551B2 (ja) | 座屈応力度の推定装置、座屈応力度の推定方法、及び座屈応力度の推定プログラム | |
JP7211287B2 (ja) | 座屈応力度の推定装置、座屈応力度の推定方法、及び座屈応力度の推定プログラム | |
Zhou et al. | Distortional buckling calculation method of steel-concrete composite box beam in negative moment area | |
JP7211281B2 (ja) | 座屈応力度の推定装置、座屈応力度の推定方法、及び座屈応力度の推定プログラム | |
JP7211282B2 (ja) | 座屈応力度の推定装置、座屈応力度の推定方法、及び座屈応力度の推定プログラム | |
Akgun et al. | Nonlinear Static analysis of laminated composite hollow beams with super-elliptic cross-sections | |
JP7211285B2 (ja) | 座屈応力度の推定装置、座屈応力度の推定方法、及び座屈応力度の推定プログラム | |
Al-Bermani et al. | Stability of cold-formed members | |
He et al. | A high order finite strip transfer matrix method for buckling analysis of single-branched cross-section thin-walled members | |
Pouladkhan | Finite Element Analysis of Linear Buckling of Circular Thin Plates for Different Boundary Conditions | |
JP7211286B2 (ja) | 座屈応力度の推定装置、座屈応力度の推定方法、及び座屈応力度の推定プログラム | |
Zhang et al. | Mathematical modeling in plane strain bending | |
JP2021006788A (ja) | 変位の推定装置、変位の推定方法、及び変位の推定プログラム | |
WO2020262306A1 (ja) | 変位の推定装置、座屈応力度の推定装置、変位の推定方法、座屈応力度の推定方法、変位の推定プログラム、及び座屈応力度の推定プログラム |
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 |