结构动力学方程常用数值解法教学文案.doc

ID:57163441

大小:627.50 KB

页数:17页

时间:2020-08-04

结构动力学方程常用数值解法教学文案.doc_第1页
结构动力学方程常用数值解法教学文案.doc_第2页
结构动力学方程常用数值解法教学文案.doc_第3页
结构动力学方程常用数值解法教学文案.doc_第4页
结构动力学方程常用数值解法教学文案.doc_第5页
资源描述:

《结构动力学方程常用数值解法教学文案.doc》由会员上传分享,免费在线阅读,更多相关内容在工程资料-天天文库

1、结构动力学方程常用数值解法精品文档结构动力学方程常用数值解法对于一个实际结构,由有限元法离散化处理后,动力学方程可写为:从数学角度看,这是一个常系数的二阶线性常微分方程组,计算数学领域,常微分数值算法常用的有两大类:-、针对一阶微分方程数值积分法发展的欧拉法,中点法,Rugge-kutta(龙格—库塔)方法。二、直接基于二阶动力学方程发展的方法。对结构动力学问题的数值求解,常用的有两大类:一是坐标变换法,它是对结构动力方程式,在求解之前,进行模态坐标变换,实际上就是一种Rize变换,即把原物理空间的动力方程变换到模态空间中去求解。现在,普遍使用的方

2、法是模态(振型)迭加法。二是直接积分法,它是对结构动力方程式在求解之前不进行坐标变换,直接进行数值积分计算。这种方法的特点是对时域进行离散,然后将该时刻的加速度和速度用相邻时刻的各位移线性组合而成。通常又称为逐步积分法。模态迭加方法,比较常用,但如下情况通常使用直接积分方法(即求解之前不进行模态分析)一、非比例阻尼,非线性情况。二、有冲击作用,激起高频模态,力作用持续时间较短,模态迭加计算量太大。一振型迭加法与Duhamel积分数值解按照有限单元法的一般规则,经过边界条件的约束处理,结构在强迫振动时多自由度体系的运动平衡方程可以表示为:(1)收集于

3、网络,如有侵权请联系管理员删除精品文档其中,是体系的质量矩阵,C是体系的阻尼矩阵,而K则是刚度矩阵.R为外荷载向量.、和则分别是体系单元节点的位移、速度和加速度向量.上述动力平衡方程实质上是与加速度有关的惯性力和与速度有关的阻尼力及与位移有关的弹性力在时刻t与荷载的静力平衡。振型叠加法是把多自由度体系的结构的整体振动分解为与振型次数相对应的单自由度体系,求得各个单自由度体系的动力响应后,再进行叠加得出结构整体响应.振型叠加法原理是利用结构无阻尼自由振动的振型矩阵作为变换矩阵,将结构动力方程式(1)式变换成一组非耦合的微分方程.逐个地求解这些方程后,

4、将解叠加即可得到动力方程的解。将体系单元节点的位移向量表示为如下的变换形式:(2)式中的变换矩阵是由动力方程对应的无阻尼自由振动方程解出的前m阶振型矩阵.即;是与时间有关的m阶向量,的各分量称为广义位移。将式(2)代入动力方程(1)并左乘以,则可得广义位移为未知数的方程:(3)式中,,,(4)现在进一步考察式(4).考虑到特征向量的正交性,可得,(5)于是对应于振型的广义位移的平衡方程(3)可改写为收集于网络,如有侵权请联系管理员删除精品文档(6)其中,为特征值(7)将式(2)稍加运算可得广义位移用有限元位移表示的形式(8)在(6)式中,当忽略了阻

5、尼的影响,平衡方程为互不耦合的,可以对每个方程逐个地进行时间积分.出于相同的考虑,在对有阻尼的体系进行分析时仍然希望采用相同的计算过程去求解互不耦合的平衡方程式.问题是式(6)中的阻尼阵C通常不能象体系的质量阵和刚度阵那样由单元的刚度阵和质量阵装配而成.但当假定阻尼与固有频率成比例,即假定(9)式中,是振型阻尼参数;是Kronecker符号(当时,=1.当时,=0)。这时式(6)可简化为如下形式的若干个方程式(10)其中的初始条件为下式收集于网络,如有侵权请联系管理员删除精品文档,(11)式(10)表示了一个具有单位质量,刚度为的自由度体系当阻尼比

6、为时的运动平衡控制方程。这个平衡方程的求解可通过计算Duhamel积分求得。(12)式中(13)当利用式(9)来考虑阻尼的影响时意味着假设结构的总阻尼是每个振型的阻尼之和,而每个振型上的阻尼是能够量测的,况且在大多数情况下结构的阻尼比更易于量测。因而便于用来近似地反映结构体系的阻尼特性。同时在计算上也避免计算阻尼阵而只需计算刚度阵和质量阵。积分递推公式对以上方程式(10),考虑某一模态的振动,并略去下标可写为(14)在初始条件,(15)下的定解为收集于网络,如有侵权请联系管理员删除精品文档(16)式中,,将上式对时间求导,得(17)将上式中的,分别

7、代以,(其中为时间步长),并按抛物线法则计算式中的卷积,有下式:(18)收集于网络,如有侵权请联系管理员删除精品文档(19)以上两式以矩阵表示为:(20)记,,,,,,(21)则上述A矩阵元素为,,,,,,,应用上述递推公式,以前一时刻来求后一时刻的结果。计算不重复。当求出后,以后在时域中的步进求解只是一些简单的数组相乘。计算速度很快。二Newmark类方法下面讨论直接积分法(1)(2)收集于网络,如有侵权请联系管理员删除精品文档(3)将(3)代入(1),(2)得:1、可以直接略去高阶项2、用变权来调节然后假设在时刻近似满足运动方程通过变换将速度和

8、加速度用位移表示,代入运动方程,只剩n+1时刻位移一个未知数,得参数不同选取包含着三个经典算法(1)Newmark平均加速

当前文档最多预览五页,下载文档查看全文

此文档下载收益归作者所有

当前文档最多预览五页,下载文档查看全文
正文描述:

《结构动力学方程常用数值解法教学文案.doc》由会员上传分享,免费在线阅读,更多相关内容在工程资料-天天文库

1、结构动力学方程常用数值解法精品文档结构动力学方程常用数值解法对于一个实际结构,由有限元法离散化处理后,动力学方程可写为:从数学角度看,这是一个常系数的二阶线性常微分方程组,计算数学领域,常微分数值算法常用的有两大类:-、针对一阶微分方程数值积分法发展的欧拉法,中点法,Rugge-kutta(龙格—库塔)方法。二、直接基于二阶动力学方程发展的方法。对结构动力学问题的数值求解,常用的有两大类:一是坐标变换法,它是对结构动力方程式,在求解之前,进行模态坐标变换,实际上就是一种Rize变换,即把原物理空间的动力方程变换到模态空间中去求解。现在,普遍使用的方

2、法是模态(振型)迭加法。二是直接积分法,它是对结构动力方程式在求解之前不进行坐标变换,直接进行数值积分计算。这种方法的特点是对时域进行离散,然后将该时刻的加速度和速度用相邻时刻的各位移线性组合而成。通常又称为逐步积分法。模态迭加方法,比较常用,但如下情况通常使用直接积分方法(即求解之前不进行模态分析)一、非比例阻尼,非线性情况。二、有冲击作用,激起高频模态,力作用持续时间较短,模态迭加计算量太大。一振型迭加法与Duhamel积分数值解按照有限单元法的一般规则,经过边界条件的约束处理,结构在强迫振动时多自由度体系的运动平衡方程可以表示为:(1)收集于

3、网络,如有侵权请联系管理员删除精品文档其中,是体系的质量矩阵,C是体系的阻尼矩阵,而K则是刚度矩阵.R为外荷载向量.、和则分别是体系单元节点的位移、速度和加速度向量.上述动力平衡方程实质上是与加速度有关的惯性力和与速度有关的阻尼力及与位移有关的弹性力在时刻t与荷载的静力平衡。振型叠加法是把多自由度体系的结构的整体振动分解为与振型次数相对应的单自由度体系,求得各个单自由度体系的动力响应后,再进行叠加得出结构整体响应.振型叠加法原理是利用结构无阻尼自由振动的振型矩阵作为变换矩阵,将结构动力方程式(1)式变换成一组非耦合的微分方程.逐个地求解这些方程后,

4、将解叠加即可得到动力方程的解。将体系单元节点的位移向量表示为如下的变换形式:(2)式中的变换矩阵是由动力方程对应的无阻尼自由振动方程解出的前m阶振型矩阵.即;是与时间有关的m阶向量,的各分量称为广义位移。将式(2)代入动力方程(1)并左乘以,则可得广义位移为未知数的方程:(3)式中,,,(4)现在进一步考察式(4).考虑到特征向量的正交性,可得,(5)于是对应于振型的广义位移的平衡方程(3)可改写为收集于网络,如有侵权请联系管理员删除精品文档(6)其中,为特征值(7)将式(2)稍加运算可得广义位移用有限元位移表示的形式(8)在(6)式中,当忽略了阻

5、尼的影响,平衡方程为互不耦合的,可以对每个方程逐个地进行时间积分.出于相同的考虑,在对有阻尼的体系进行分析时仍然希望采用相同的计算过程去求解互不耦合的平衡方程式.问题是式(6)中的阻尼阵C通常不能象体系的质量阵和刚度阵那样由单元的刚度阵和质量阵装配而成.但当假定阻尼与固有频率成比例,即假定(9)式中,是振型阻尼参数;是Kronecker符号(当时,=1.当时,=0)。这时式(6)可简化为如下形式的若干个方程式(10)其中的初始条件为下式收集于网络,如有侵权请联系管理员删除精品文档,(11)式(10)表示了一个具有单位质量,刚度为的自由度体系当阻尼比

6、为时的运动平衡控制方程。这个平衡方程的求解可通过计算Duhamel积分求得。(12)式中(13)当利用式(9)来考虑阻尼的影响时意味着假设结构的总阻尼是每个振型的阻尼之和,而每个振型上的阻尼是能够量测的,况且在大多数情况下结构的阻尼比更易于量测。因而便于用来近似地反映结构体系的阻尼特性。同时在计算上也避免计算阻尼阵而只需计算刚度阵和质量阵。积分递推公式对以上方程式(10),考虑某一模态的振动,并略去下标可写为(14)在初始条件,(15)下的定解为收集于网络,如有侵权请联系管理员删除精品文档(16)式中,,将上式对时间求导,得(17)将上式中的,分别

7、代以,(其中为时间步长),并按抛物线法则计算式中的卷积,有下式:(18)收集于网络,如有侵权请联系管理员删除精品文档(19)以上两式以矩阵表示为:(20)记,,,,,,(21)则上述A矩阵元素为,,,,,,,应用上述递推公式,以前一时刻来求后一时刻的结果。计算不重复。当求出后,以后在时域中的步进求解只是一些简单的数组相乘。计算速度很快。二Newmark类方法下面讨论直接积分法(1)(2)收集于网络,如有侵权请联系管理员删除精品文档(3)将(3)代入(1),(2)得:1、可以直接略去高阶项2、用变权来调节然后假设在时刻近似满足运动方程通过变换将速度和

8、加速度用位移表示,代入运动方程,只剩n+1时刻位移一个未知数,得参数不同选取包含着三个经典算法(1)Newmark平均加速

显示全部收起
温馨提示:
1. 部分包含数学公式或PPT动画的文件,查看预览时可能会显示错乱或异常,文件下载后无此问题,请放心下载。
2. 本文档由用户上传,版权归属用户,天天文库负责整理代发布。如果您对本文档版权有争议请及时联系客服。
3. 下载前请仔细阅读文档内容,确认文档内容符合您的需求后进行下载,若出现内容与标题不符可向本站投诉处理。
4. 下载文档时可能由于网络波动等原因无法下载或下载错误,付费完成后未能成功下载的用户请联系客服处理。
关闭