结构动力学思考题
made by 云屹
思考题一
1、结构动力学与静力学的主要区别是什么？结构的运动方程有什么不同？
主要区别为：
(1)动力学考虑惯性力的影响，静力学不考虑惯性力的影响；
(2)动力学中位移等量与时间有关，静力学中位移等量不随时间变化；
(3)动力学的求解方法通常与荷载类型有关，静力学一般无关。
运动方程的不同：
动力学的运动方程包括位移项、速度项和加速度项；静力学的平衡方程只包括位移项。
2、什么是动力自由度？什么是静力自由度？区分动力自由度和静力自由度的意义是什么？动力自由度：确定结构体系质量位置的独立参数；
静力自由度：确定结构体系在空间中的几何位置的独立参数。
意义：通过适当的假设，当静力自由度数大于动力自由度数时，使用动力自由度可以减少未知量，简化计算，提高计算效率。
3、采用集中质量法、广义坐标法和有限元法都可以使无限自由度体系简化为有限自由度体
4、在结构振动的过程中引起阻尼的原因有哪些？
(1)材料的摩擦或材料变形引起的热耗散；
(2)构件连接处或结构构件与非结构构件之间的摩擦；
(3)结构外部介质的阻尼。
5、在建立结构运动方程时，如考虑重力的影响，动位移的运动方程有无改变？
如果满足条件：
(1)线性问题；
(2)重力的影响预先被平衡；
则动位移的运动方程不会改变，否则会改变。
思考题二
1、刚度系数k ij和质量系数m ij的直接物理意义是什么？如何直接用m ij的物理概念建立梁单元的质量矩阵[M]？
k ij：由第j自由度的单位位移所引起的第i自由度的力；
m ij：由第j自由度的单位加速度所引起的第i自由度的力。
依次令第j（j=1,2,3,4）自由度产生单位加速度，而其他的广义坐标处保持静止，使用平衡方程解出第i自由度上的力，从而得到m ij，集成得到质量矩阵[M]。
2、如何用刚度矩阵和质量矩阵，以矩阵的形式表示多自由度体系的势能和动能？
{}[]{}1=
2T
T u M u {}[]{}1=2T
V u K u
3、建立多自由度体系运动方程的直接动力平衡法和拉格朗日方程法的优缺点是什么？ (1)直接动力平衡法：
优点：概念直观，易于通过各个结构单元矩阵建立整体矩阵，便于计算机编程。 缺点：涉及矢量计算，通常计算较繁琐；涉及叠加原理，不易处理非线性问题。 (2)拉格朗日方程法：
优点：仅涉及标量计算；求解不限于线性问题，适用围广。 缺点：不便计算机编程，不适用于大规模问题。
4、什么是几何刚度，几何刚度主要与什么量有关，几何刚度对结构动力特性有什么影响？ 几何刚度：表示结构在变形状态下的刚度变化。（轴力引起的附加弯矩的影响） 几何刚度主要与轴力的大小及构件的几何形状与尺寸有关。
几何刚度会产生P-Δ效应，改变结构的动力特性。压力降低刚度，拉力增加刚度。
5、什么是结构动力问题分析中的静力凝聚法？动力自由度的概念是什么？静力凝聚法在结构动力问题分析中可起什么作用？
静力凝聚法：当静力自由度数目大于动力自由度时，消去广义质量为零或很小的广义坐标，从而缩减结构体系自由度数目的方法。
动力自由度：确定结构体系质量位置的独立参数。 作用：缩减计算规模，提高计算效率，降低计算量。
6、试证明多自由度体系的位能和动能分别为：
1111
1=21
2N N
ij i j
i j N
N
ij i j i j V k u u T m u u =====
∑∑∑∑
证明：
弹性恢复力所做的功为
()()1111111
11112
11
22
N
N
N N N N
B B
B
ij i ij j i ij j i ji i j A A A
i j i j i j N
N N N
B
ij i j ij iB jB iA jA A
i j i j W F du k u du k u du k u du k du u k u u u u =============+==-∑∑∑∑∑∑???∑∑
∑∑?
故定义弹性位能为：
11
1=2N N
ij i j i j V k u u ==∑∑
惯性力所做的功为
()()1111
1
1
1
111
1211
22N
N N N
N
N
B B
B
ij i ij j i ij j i ji i j A
A
A i j i j i j N
N
N
N
B
ij i j ij iB jB iA jA A i j i j W F du m u du m u du m u du m du u m u u u u =============+==-∑∑∑∑∑∑
?
??∑∑
∑∑?
注意到=
j j i i i j du u du du u du dt
=
故定义动能为：
11
12N N
ij i j i j T m u u ===∑∑
7、如何充分论证，当多自由度体系的动力自由度不能充分确定体系的几何位置时，初始建立的运动方程组中一定含有非动力自由度的静力自由度？ 证明：
假设初始建立的运动方程组不含非动力自由度的静力自由度，则质量矩阵[M]满秩，则动力自由度可以充分确定体系的几何位置，与前提条件矛盾。
8、在推导拉格朗日方程时，给出了以下几个基本表达式： 位移：()12,,...,;i i n u u q q q t =
(1) 动能：()1212,,...,;,,...,n n T T q q q q q q = (2) 势能：()12,,...,n V V q q q =
(3)
问题：
(1)式为什么显含时间t ？
有时位移中可以存在已知的显含的时间t 的函数，比如地基运动问题。 (2)式(2)中是否应显含时间t ？
如果位移显含t ，由于动能是位移的函数，也应该显含时间t 。
(3)难道广义坐标及速度完全确定后，体系的动能还与时间t 有关系？ 可以有关系，比如地基运动问题中，体系的动能就与时间t 有关系。 (4)势能中是否应显含时间t ？
如果位移显含时间t ，由于势能是相对位移的函数，也可能会显含时间t 。 (5)为什么在变分运算时，不对显含的时间t 进行运算？
因为显含时间的函数随时间的变化规律是已知的，它的变分为零，即δt=0。
(6)若式(2)、(3)中显含时间t ，对拉格朗日方程的推导是否有影响？ 由于时间t 的变分为零，对拉格朗日方程的推导没有影响。（实质上是方程的边界条件）
思考题三
1、在振动过程中产生阻尼的原因有哪些？什么叫临界阻尼？什么叫阻尼比？怎样测量结构
振动过程中的阻尼比？一般建筑结构的阻尼比是多少？ 产生阻尼的原因：
(1)材料的摩擦或材料变形引起的热耗散； (2)构件之间或构件与非构件之间的摩擦； (3)结构外部介质的阻力。
临界阻尼：使体系自由振动反应中不出现往复振动所需的最小阻尼值。 阻尼比：体系中实际阻尼系数与临界阻尼系数的比值。 测量结构阻尼比的方法： (1)对数衰减率法； (2)共振放大法； (3)半功率带宽法； (4)等效粘性阻尼法。 对于钢结构，ζ=0.01左右；
对于混凝土结构，脉动荷载下ζ=0.03左右，地震下ζ=0.05左右。
2、分析临界阻尼体系自由振动的可能运动形式及其满足的条件。
()()00n u u ω≥-时，位移不会变号； ()()00n u u ω<-时，位移会变号。
3、阻尼对结构的自振频率有什么影响？阻尼变大，结构的自振周期如何变化？
由2
1D ωωζ=-当ζ<1时，自振频率会变小，但当阻尼比较小时（ζ<0.2），这一影响可忽略不计。
当ζ<1时，阻尼变大，结构的自振周期变大。
4、为什么说自振周期是结构的固有特性？它与结构哪些固有量有关？
动荷载及初始条件确定后，结构的动力响应就仅由结构的自振周期（自振频率）控制。自振频率与结构的质量、刚度及阻尼比有关。
5、什么是动力放大系数？动力放大系数的大小与哪些因素有关？单自由度体系位移的动力
放大系数与力的动力放大系数是否一样？
动力放大系数：动荷载引起的响应幅值与动荷载幅值作为静荷载所引起的结构静响应的比值。
简谐荷载下的动力放大系数与频率比（自振频率、荷载频率）、阻尼比有关：
d R =
当惯性力与动荷载作用线重合时，位移动力系数与力动力系数相等，否则不相等。原因是：
当把动荷载换成作用于质量的等效荷载时，引起的质量位移相等，但力并不等效，根据动力系数的概念可知不会相等。
6、根据动力放大系数分析，什么时候动力放大系数Rd->1，如何理解下述结论：“随时间变化很慢的动荷载实际上可看作静荷载”。这里“很慢”的标准是什么？ 当ω->0时，Rd->1。
根据上述结论，当动荷载的频率很小时，动力放大系数趋于1，动荷载可以看作静荷载。 “很慢”的标准是惯性力相对于总荷载可忽略不计。
7、单自由度体系动荷载作用点不在体系的集中质量上时，动力计算如何进行？此时，体系中的动力放大系数是否仍然一样？
通过动平衡方程或虚位移原理，将原动荷载用沿自由度方向作用于质量上的等效动荷载代替。集中质量位移的动力放大系数仍然一样，但体系其他部位的位移以及力的动力系数通常不再相同，即不能采用统一的动力系数。
8、简谐荷载作用下有初始条件影响的无阻尼单自由度体系动力反应的瞬态反应项中
()
02
/sin 1/n
n n p t k ωωωωω- 一项是如何产生的，它与外荷载和初始条件的关系如何？ 是由外荷载产生的伴生自由振动，作用是使求得的解满足初始条件，它与外荷载的幅值和频率有关，与初始条件无关。
9、什么是共振？什么是共振频率？结构位移反应、速度反应和加速度反应的共振频率是否相同？ 定义一：共振是指体系在动荷载作用下振幅最大的情形，相应的动荷载的频率称为共振频率。 定义二：共振是指体系自振频率与动荷载频率相同而使振幅变得很大的一种现象。 当体系无阻尼时，结构位移反应、速度反应和加速度反应的共振频率相同；当体系有阻尼时，结构位移反应、速度反应好加速度反应的共振频率不同。
10、无阻尼体系和有阻尼体系的自振频率和共振频率是否相同？分别为多少？ 不相同，分别为：
思考题四
1、在结构动力反应分析中采用的阻尼理论有哪几种？各有什么特点？ (1)粘滞阻尼：大小与速度成正比 (2)摩擦阻尼：大小为常数
(3)滞变阻尼：大小与位移成正比
(4)流体阻尼：大小与速度的平方成正比
2、加速度计和位移计的设计原理是什么？如何设计速度计？
加速度计：在所量测的频段（低频段，0.5n ωω≤）动力放大系数接近常数。 速度计：在所量测的频段（高频段，2n ωω>）动力放大系数接近常数。
速度计：在所量测的频段（中频段）动力放大系数接近常数。
3、用拟静力试验（往复加载的静力试验）测量结构构件阻尼比的原理是什么？如何实现？ 原理是阻尼耗能与加载频率关系不大。
实现方法是通过拟静力试验测出一个周期的阻尼耗能E D ，从而计算出等效粘滞阻尼比：
20
2D
eq n E ku ζωπω=
4、测量结构阻尼比的方法有几种？每一方法的优点和缺点是什么？ (1)对数衰减率法
优点：测量一阶振型的阻尼比比较容易。
缺点：确定高阶振型的阻尼比时，要能够激发出相应振型的自由振动，这一点比较困难。 (2)共振放大法
优点：方法简单，且可处理任意类型的多自由度体系。 缺点：等效静位移较难确定。 (3)半功率带宽法 优点：可操作性强。
缺点：对多自由度体系要求共振频率稀疏。 (4)能量等效阻尼比 优点：可操作性强。
缺点：对于共振频率之外的其他频率，物理概念不一定正确，只是一个近似。
5、简谐荷载作用下，在结构的一个振动循环中，外力、阻尼力、弹性恢复力和惯性力做工
及其关系如何？
在稳定状态下，外力和阻尼力所做的功和为零，弹性恢复力和惯性力所做的功为零。
6、结构中阻尼的来源以摩擦型阻尼为主，为什么实际结构动力反应分析中采用的结构阻尼是滞变阻尼而不采用经典的库伦摩擦阻尼？仅仅是出于计算上的方便？ (1)在简谐反应分析（频域分析）中，滞变阻尼理论与试验结果相符； (2)函数连续可微，便于计算；
(3)结构的等效阻尼比可以通过试验测得，而摩擦系数较难通过试验确定。
7、滞变阻尼（复阻尼）的三种形式在复数域是完全等价的，但在一个振动循环的耗能确不相同，原因是什么？
复数域上取模后会丢失一部分信息（不再一一对应），做围道积分时会产生差异，故第一种形式与后两种有所不同。
8、第1和第3种形式的滞变阻尼在实数域的定义不尽相同，但在复数域则完全等价，若分别采用这两种形式的滞变阻尼进行结构动力反应分析，是否预示着采用实数域分析和复数域分析会获得矛盾的结果？
因为实数域与复数域并不是一一对应的关系，所以可能会获得矛盾的结果。（扯的）
思考题五
1、在杜哈梅积分中时间变量τ和t 有什么区别？怎样用杜哈梅积分求解任意动荷载作用下的动力位移问题？简谐荷载下的动位移可以用杜哈梅积分求解吗？
积分上限t 是原函数的自变量，是动力响应发生的时刻；τ是积分变量，是瞬时冲量作用的时刻。
对于无阻尼体系：
()()()0
1sin t
n
n
u t p t d m τωττω=
-?
???? 对于阻尼体系：
()()()()0
1
exp sin t
n
D
D
u t p t t d m τζωτωττω=
---?
???????? 简谐荷载可以使用杜哈梅积分。
2、采用连续傅里叶变换和离散傅里叶变换研究非周期荷载作用下体系动力反应问题时的最主要差别是什么？在采用离散傅里叶变换分析时都应注意哪些问题？
采用连续傅里叶变换时仍为非周期函数，而离散傅里叶变换会将非周期函数周期化。 注意事项：
(1)离散傅里叶变换将非周期时间函数周期化；
(2)对荷载P(t)要增加足够多的零点以增大持续时间Tp ，以保证在所计算的时间段[0,Tp]体系位移能衰减到0；
(3)频谱上限频率（奈奎斯特频率）为1/2Nyquist f t =?；
(4)频谱的分辨率为1/p f T ?=； (5)频谱的下限11/p f T =。
3、对比基于傅里叶级数和离散傅里叶变换得到的结构动力反应的解式，分析两者之间的异同。
傅里叶级数：()()()exp j
j
j u t H i i t ωω∞
=-∞
=
∑
离散傅里叶变换：()()1
1
2exp N k j j p kj u t U i T N πω-=??
=
- ???∑
相同点：
(1)计算的结果都具有周期性； 不同点：
(1)从本质上讲，傅里叶级数处理的是周期性荷载，而离散傅里叶变换处理的是任意荷载； (2)傅里叶级数得到的结果是连续的，离散傅里叶变换得到的结果是离散的； (3)傅里叶级数是无穷级数，离散傅里叶变换是有限项求和。
4、在离散傅里叶变换中，从傅里叶谱离散化给出的离散频率点看，最大频率点为1/N f t =?，但理论上给出的上限频率却仅为1/2Nyquist f t =?，为什么？
如果采样频率高于奈奎斯特频率，将会发生频率混叠现象，不能真实还原被测信号。
5、什么是奈奎斯特频率？为什么称奈奎斯特频率为折叠频率？它有什么作用？为保证离散数值分析的精度，最大有效频率应如何取值？
奈奎斯特频率是离散信号系统采样频率的一半，即1/2Nyquist f t =?
因为超过奈奎斯特频率的信号与其关于奈奎斯特频率对称的频率的信号相同，故称奈奎斯特频率为折叠频率。
作用是确定了采样频率的上限。
由于接近奈奎斯特频率的信号在采样和重建过程中可能会产生畸变，最大有效频率一般取奈奎斯特频率的三分之二。
6、在采用离散傅里叶变换方法进行结构动力反应问题分析时，将导致“周期化”，什么是非周期问题的“周期化”？如何避免“周期化”对结构动力反应的影响？ “周期化”是指离散傅里叶变换会将非周期函数变为周期函数。
对荷载P(t)要增加足够多的零点以增大持续时间Tp ，以保证在所计算的时间段[0,Tp]体系位移能衰减到0。
7、什么是数值信号处理问题中的分辨率？如何提高分辨率？分辨率对结构阻尼比的计算有什么影响？
分辨率是指采样间隔，表达式为1/p f T ?=
可以通过增大持续时间Tp 来提高分辨率。
分辨率越高，解的精度越高，计算阻尼比时的精度通常也越高。
8、什么是反应谱，它与那些物理量有关？什么是地震影响系数，它与反应谱有什么关系？什么是动力系数，它与反应谱和地震影响系数有什么关系，与动力放大系数有什么关系？ 反应谱是指在给定的地震作用下，结构的最大（相对）位移反应和最大（绝对）加速度反应随自振周期变化的曲线，它与给定的地震动加速度时程、结构的阻尼比以及自振周期有关。 地震影响系数是以重力加速度g 为单位的反应谱， 与反应谱的关系为：/a S g α= 地震动力系数是归一化的反应谱， 与反应谱的关系为：/a g S u β=
与地震影响系数的关系为：/g g u βα=? 与动力放大系数的关系为：/d st g R u u β=?
思考题六
1、什么是多自由度体系的振型，用振型对结构的位移进行展开，即采用振型叠加法进行结构动力反应分析有什么优点？
振型是体系上所有质量按相同频率做自由振动时的振动形状。
优点是将耦联的N 个自由度问题解耦为N 个独立的单自由度问题，避免求解联立方程。
2、什么是振型的正交性？振型关于刚度阵正交的物理意义是什么？振型关于质量阵正交的物理意义是什么？
振型的正交性：当m n ≠时，
{}[]{}0T
m n M φφ= {}[]{}0T
m n K φφ=
物理意义：某一振型产生的弹性恢复力（惯性力）在另一振型上做功为零，即不同振型间能量不会传递。
3、如何证明振型的完备性？如何证明结构振型之间是线性无关的？
由于振型{}i φ为矩阵[][]1
K M -?的特征向量，而矩阵[K]、[M]都是正定矩阵，故矩阵
[][]1K M -?也为正定矩阵，则其特征向量{}i φ是线性无关的。
对于N 维空间，由于向量{}()1,2,...,i i N φ=是线性无关的，故它们构成了该线性空间的一组基，故振型具有完备性。
4、什么是振型质量Mn 和阵型刚度Kn ？它们与自振频率ωn 有什么关系？对应于结构某阶振型，振兴质量Mn 和阵型刚度Kn 是否为固定常数？
{}[]{}T
n n n M M φφ= {}[]{}T
n n n K K φφ=
n ω=
它们不是固定常数，但它们的比值是固定常数2
n ω。
5、对于单自由度体系通过自由振动分析可以获得结构的无阻尼自振频率ωn 和有阻尼自振频率ωD ，对于多自由度有阻尼体系，如何获得结构的自振频率和振型？ 通过阻尼的解耦假定，仍可得到与单自由度体系类似的结论：
D ω=
振型与无阻尼体系一致。
6、振型叠加法用到了叠加原理，什么情况下能用这个方法？什么情况下不能用？ 只有线性体系才能使用叠加原理。
如果出现材料非线性、几何非线性或运动非线性，都不能使用叠加原理。
7、在多自由度体系振型阻尼比的现场动力测量时，可以采用自由振动试验法，此时需要使结构按不同振型做自由衰减振动，如何使多自由度体系只按某个特定的振型振动？ 各质点之间的初位移和初速度的比值应具有该振型的比值关系。
8、N 个自由度的体系有多少发生共振的可能性？为什么？ 有N 个，因为有N 种不同的振型。
9、多自由度体系的频率方程存在重根时，体系自振频率个数、振型个数与自由度数关系如何？各振型之间的关系如何？ 自振频率个数少于自由度数。 振型个数等与自由度数。 各振型之间仍然保持正交性。
10、如何判断频率方程是否存在重根及其为几重根？ 记()[][]2f
K M ωω=-
则若()
()()01,2,...,j i f
j k ω==且()()10k i f ω+≠
便在ωi 处存在k 重根。
11、什么是矩阵的正定条件？体系刚度矩阵和质量矩阵的正定条件是否能保证频率方程不出现重根？
对于对称矩阵A ，若对于任意非零向量x ，均满足不等式0T
x Ax >，则称矩阵A 为正定矩阵。
不能保证。
思考题七
1、为什么阻尼会对结构振型的正交性产生影响，什么时候阻尼称为经典阻尼？什么时候称为非经典阻尼？
因为推到正交性时没有考虑阻尼。（而使用复模态时就考虑了阻尼） 经典阻尼是指满足振型正交条件的阻尼。（阻尼矩阵为对角阵） 非经典阻尼是指不满足振型正交条件的阻尼。（阻尼矩阵为非对角阵）
2、当结构的阻尼为非经典阻尼时，采用振型叠加法计算结构动力反应时，避免求解联立方程组的两种基本分析方法是什么？各有什么优缺点？ (1)迭代法
优点：形式简单直观；
缺点：需要计算多步，计算量可能较大。 (2)复模态方法
优点：可以得到解耦的独立方程；
缺点：矩阵的阶数增大，且自振频率和振型都是复数，要花费更多的时间。
3、Rayleigh 阻尼的概念和特点，确定Rayleigh 阻尼公式中两参数的原则是什么？ Rayleigh 阻尼的定义：假设结构的阻尼矩阵是质量矩阵好刚度矩阵的组合，即 [C]=a 0[M]+a 1[K]
特点：满足振型正交条件。
原则：所有感兴趣的频率都在频段,i j ωω????。
4、扩展的Rayleigh 阻尼（Caughey 阻尼）的概念；用多个自振频率和振型阻尼比确定扩展Rayleigh 阻尼的常数时，在自振频率个数的选取上应注意的基本原则是什么？ Caughey 阻尼的概念：[][]
[][]()
11
L l
l
l C M a M K --==∑
基本原则：频率点应取偶数个，防止出现负阻尼。
5、为什么高阶振型对结构动力反应的影响小？ 对于地震荷载，外荷载的频率较小，故对于高阶振型，
n
ω
ω较小，则动力放大系数趋于1。
6、当结构不同部分的阻尼比存在明显差异时，如何较高精度地实现结构地震反应的振型分解反应谱分析？
(1)先将结构分为几个部阻尼接近的子结构；
(2)各子结构阻尼是经典的，可分别建立阻尼矩阵； (3)集成后得到结构总体阻尼阵。
7、构造结构阻尼矩阵的目的是什么？为什么采用Rayleigh 阻尼假设？当结构由阻尼相差较大的几部分构成时，结构体系的阻尼矩阵如何建立？ 目的：在有阻尼体系的动力分析问题中使用振型叠加法。
采用Rayleigh 阻尼假设是为了构造出既满足振型正交条件又具有一定精度的阻尼矩阵，从而简化计算。 建立方法：
(1)先将结构分为几个部阻尼接近的子结构；
(2)各子结构阻尼是经典的，可分别建立阻尼矩阵； (3)集成后得到结构总体阻尼阵。
8、Rayleigh 阻尼是一种经典阻尼，满足振型正交条件，用振型叠加法分析经典阻尼结构的动力反应问题时，是否需要采用Rayleigh 阻尼假设并构造阻尼矩阵？
对于经典阻尼结构，阻尼矩阵已经满足了振型正交条件，无需再使用Rayleigh 阻尼假设。
9、什么是振型阻尼比？实际工程中不同阶振型阻尼比的变化规律如何？数值计算时一般如何选取？
振型阻尼比：某一特定振型下的阻尼比，即
2n
n n n
C M ξω=
实际工程中振型阻尼比随振幅的增大而增大，且因振型的不同而不同。 数值计算中一般取为与振型阶数无关的定值。 10、什么是振型加速度？什么是振型加速度法？什么是静力修正法？两种分析方法有什么异同？
振型加速度：振型位移q n (t)关于时间的二阶导。
振型加速度法：在叠加公式中使用振型加速度法（包括振型速度）。 静力修正法：对于高阶振型，采用简化的静力分析方法。
除数值计算引起的误差外，静力修正法好振型加速度法给出的结果是相同的。
思考题八
1、建立杜哈梅积分公式和Newmark-β法时域逐步积分公式时最主要的区别是什么？如果在
用杜哈梅积分求解任意动荷载作用下的反应问题时，先将时间τ等步距离散化，然后采用数值积分，可否用杜哈梅积分求解非线性问题？
最主要区别是杜哈梅积分公式得到的是精确的解析解；而Newmark-β法得到的是近似的数值解。
不能，因为杜哈梅积分公式建立在叠加原理的基础上，不能求解非线性问题。
2、用中心差分逐步分析方法计算结构非线性动力反应问题时，在每一步的计算中能否像Newmark-β法那样实现对非平衡力的迭代修正计算？
不能，中心差分逐步分析方法使用的是原始的运动方程，而非增量运动平衡方程，不涉及位移增量的计算，运动状态是一个整体，无法进行迭代修正。
3、结构阻尼的存在有助于减小结构的动力反应，阻尼是否有助于控制时域逐步积分法引起的结构振荡失稳，即提高动力问题数值算法的稳定性？
不一定，比如在Clough 格式的显式中心差分法中，当阻尼比超过0.5时，稳定性随阻尼比的增大而变差。
4、Wilson-θ法将引起结构动力反应振幅的进一步衰减，这称为算法阻尼，试讨论算法阻尼存在的利弊。
利：通常能提高算法的稳定性。 弊：会降低计算结果的精度。
5、时域逐步积分法的计算精度和稳定性的概念及两者之间的关系。 计算精度：截断误差与时间步长之间的关系。
稳定性：随时间步数的增大，数值解是否会远离精确解。 一般而言，计算精度与稳定性正相关，但也有例外。
6、建立一种时域逐步积分算法，要进行包括收敛性、计算精度、稳定性和计算效率四方面基本问题的分析，如何进行算法的收敛性和计算精度的分析？ 按照定义进行，分析截断误差及计算结果随时间步长的变化规律。
7、显示算法不需要求解联立方程组，而隐式算法需要求解联立方程组。如何证明“所有的显示算法都是有条件稳定的，而隐式算法可以是有条件稳定或无条件稳定”这一论断？ 对于显示算法，所有的系数矩阵均为对角阵，故计算公式可化简为：
{}()[]{}{}1i i i u f t A u b +=?+
式中，A 为对角阵，若要求算法稳定，则
()[]()()[]()1f t A f t A ρρ?=?≤
显式地存在t?，必然对t?有所要求，故显示算法都是有条件稳定的。
f t?可能被消掉，故可以是无条件稳定的。
对于隐式算法，()
8、结构阻尼常用于描述结构线弹性动力反应时的耗能效应，对于结构弹塑性反映问题，塑形反应将引起结构振动能量的耗散，称为塑性耗能。在结构弹塑性反应问题分析中是否还应继续考虑阻尼耗能？如果不考虑，能否出现问题？如果要考虑，如何考虑？
要考虑，
(1)塑性耗能只是耗能的一种，还有摩擦耗能、环境阻力耗能等的影响；
(2)结构并不是一直处在塑性状态中，在弹性状态不考虑阻尼耗能会产生较大误差。
思考题九
1、什么是欧拉-伯努利梁？什么是铁木辛柯梁？剪切变形和转动惯量对梁的自振频率和变形有什么影响？
欧拉-伯努利梁：仅考虑弯曲变形的梁。
铁木辛柯梁：除了弯曲变形外，还考虑了剪切变形和转动惯量的影响。
剪切变形和转动惯量会使结构的自振频率降低，并会改变阵型，且阶数越高，影响越大。
2、欧拉-伯努利梁和铁木辛柯梁在梁的变形基本假设上的异同。
相同点：都使用了平截面假定；
不同点：欧拉-伯努利梁仅考虑了弯曲变形，而铁木辛柯梁还考虑了剪切变形和转动惯量的影响。
3、均匀梁有多少个自振频率？对于实际建筑结构，粗略推断作为结构构件的单梁（柱）自振频率和整体结构自振频率之比值的可能围。
有无穷多个自振频率。
可能会差很多。
4、轴力对弯曲梁和剪切梁的自振频率有什么影响？轴力对杆的自振频率是否有影响？轴力是否影响梁的振型？
拉力会使得梁的自振频率增大，压力会使得梁的自振频率减小。
对于杆，轴力是外荷载，对其自振频率没有影响。
轴力不影响梁的振型。
5、梁的每一个端点分别有位移、转角、弯矩、剪力四个条件，为什么仅能提供两个边界条件？
因为位移与剪力不独立，转角与弯矩不独立。
6、欧拉-伯努利梁为什么被称为纯弯曲梁？是梁中不存在剪力吗？
因为欧拉-伯努利梁中只有弯曲变形。
存在剪力，但不存在剪切变形。
7、铁木辛柯梁是否是完全合理、精确地描述均直梁变形关系的理论？
不是，铁木辛柯梁仍然用到了平截面假定。
8、试叙述建立铁木辛柯梁时的几何方程、物理方程和平衡方程。
见教材。
9、横向变形梁的平截面假设有几种，各有什么特点？对应着什么梁？
欧拉-伯努利梁：平面截面变形后仍为平面，且变形后的截面与变形后的轴线垂直。
铁木辛柯梁：同上。
均匀剪切梁：平面截面变形后仍为平面，且变形后的截面与变形前的轴线垂直。
轴向杆：同上。
10、如果要更深入、系统地研究梁截面的变化、变形规律及平截面假设带来的影响，都可以采用什么方法开展分析？
可以采用弹塑性力学数值计算的方法来开展分析。
11、不同梁理论（轴向，剪切，弯曲，弯剪）的关键不同点在何处？
关键不同点在于对变形形态的假设不同，轴向杆只发生轴向变形，均匀剪切梁只发生剪切变形，弯曲梁只发生弯曲变形，弯剪梁发生弯曲变形和剪切变形。
思考题十
1、当均匀梁中有一集中质量M时，梁的运动方程有什么变化，梁的振型如何推导？
运动方程中会多出一项惯性力。
在推导振型时，可以将集中质量当作边界条件来处理，即通解不变，但系数会有所变化。
2、如果在进行结构模态分析时未考虑梁上集中质量M的影响，得到的振型是否可以用于梁的动力反应分析（用振型叠加法）？如果可以，会出现什么问题？
可以，此时将集中质量的惯性力视为外荷载。
运动方程会变成耦联的。
3、当均匀简支梁的梁端支座产生竖向运动时，如何求解支座运动引起的简支梁的动力反应？可将支座运动等效为线性分布的动力荷载。
4、如何将一根梁的研究成果推广用于具有分布参数的框架结构的自振频率和振型的分析？可以使用动力直接刚度法。将梁单元的动力刚度矩阵组装成结构的动力刚度矩阵，并进行求解。
5、一根梁的自振频率ωn与n或n2成正比，为什么由梁组成的大型复杂结构体系，例如大
跨悬索桥的自振频率不按这一规律变化，而表现出频率密集的性质？
我！不！知！道！
思考题十一
1、有限元法和一般广义坐标法的试函数（有限元法中称插值函数或形函数）及广义坐标有
2、集中质量和一致质量有限元法的差异和优缺点，采用这两种有限元模型给出的结构自振频率与实际结构自振频率的关系。
3、在结构有限元模型中，动力自由度和静力自由度是否相同？
对于一致质量法，是相同的；
对于集中质量法，是不同的。
4、对于一致质量和集中质量有限元模型，如果采用Rayleigh阻尼，与质量矩阵成正比的阻尼比有什么不同？与刚度阵成正比的阻尼比是否相同？
Rayleigh阻尼不仅与质量阵、刚度阵有关，还和自振频率有关，故与质量矩阵成正比的阻尼比完全不同，与刚度阵成正比的阻尼比只是系数a1不同。
5、当采用中心差分法对无阻尼结构有限元模型的运动方程进行时域逐步积分求解时，集中质量法和一致质量法的逐步积分公式的最主要差别是什么？
集中质量法的质量阵是对角阵，公式为显式；一致质量法的质量阵是非对角阵，公式为隐式。
6、是否可以从理论上证明一致质量有限元模型的基本自振频率必定高于相应连续结构的基本自振频率？
可以，一致质量模型相当于对体系增加了约束，提高了结构刚度，故计算的自振频率不小于理论值。
思考题十二
1、结构动力反应问题中引起不确定性的来源有几种？指出工程结构动力分析中的两类不确定性问题。
有3中：
(1)结构确定，输入不确定；
(2)结构不确定，输入确定；
(3)结构和输入都不确定。
主要是前两种。
2、什么是随机变量？什么是随机过程？两者之间有什么关系？
随机变量：一个量不能预先确定，但其取值满足一定的分布。
随机过程：作为时间函数的随机变量。
3、什么是平稳随机过程？强平稳和弱平稳的定义？
平稳随机过程：统计特征不随时间的推移而发生变化的随机过程。 强平稳：随机过程的各阶矩均与时间t 无关。
弱平稳：随机过程的一阶矩和二阶矩与时间t 无关。
4、什么是宽带随机过程？什么是窄带随机过程？把平稳随机过程划分为窄带过程和宽带过程在工程上有什么实际意义？
宽带随机过程：功率谱密度函数在相当宽的频带上取有意义的数量级。 窄带随机过程：功率谱密度函数只在较窄的频带具有有意义的数量级。
实际意义：如果结构的输入谱具有窄带特征，则结构设计时应优先选择避开该频段，避免发生共振；如果结构的输入谱具有宽带特征，则应考虑通过增大结构阻尼等方式减小结构反应。
5、什么是各态历经随机过程？它与平稳随机过程的关系？
各态历经随机过程：平稳过程的一个样本包含了其他各次取样的全部特征。 平稳是各态历经的必要条件，但不是充分条件。
6、给出三种常见的随机过程描述。 无聊不无聊！
7、功率谱密度函数和自相关函数的定义，什么是Wiener-Khintchin 定理。 功率谱密度函数：()()21
lim
x T S F T
ωω→∞=
自相关函数：()()(),R t t E x t x t ττ+=+????
Wiener-Khintchin 定理：功率谱密度函数和自相关函数互为Fourier 变换对。
8、在结构动力学课程介绍的线性结构随机反应的频域分析方法中如何反映或说如何实现结构动力反应的不确定性计算。
(1)输入的功率谱密度函数()x S ω已知，确定结构的复频响应函数()H i ω； (2)计算结构反映的功率谱密度函数()()()2
y x S H i S ωωω=；
(3)通过Fourier 逆变换计算自相关函数()()()1
exp 2y y R S i d τωωτωπ
+∞
-∞
=
?
；
(4)假设结构相应的均值，计算其方差。
思考题十三
1、推导Rayleigh-Ritz 法时，通过对Rayleigh 熵取极值得到了Rayleigh-Ritz 法的基本公式，取极值的物理背景是什么？
增加约束会使体系的频率变大，Rayleigh 法得到的是固有频率的上限，故取极小值最接近真实的固有频率。
2、Rayleigh-Ritz 法首先假设一组振型，通过在一个减缩空间的模态分析获得结构体系的一组振型和自振频率。能否利用所获得的振型继续用Rayleigh-Ritz 法进行迭代分析以获得精度
更好的一组振型和自振频率？
不能，Rayleigh-Ritz 法求得的已经是线性空间中Rayleigh 熵最小的解了，继续使用Rayleigh-Ritz 法不会改变结果。
3、在用振型叠加法分析结构动力反应问题时，可以采用自振振型，也可以采用Ritz 向量，与普通自振振型相比，荷载相关的Ritz 向量的优缺点是什么？ 优点：对于相应的外荷载收敛很快，且易于编程。
缺点：与外荷载的分布形式有关，当外荷载的分布形式变化时，需要重新计算相应的Ritz 向量。
4、Lanczos 向量与Ritz 向量的关系，用Lanczos 方法得到的振型与Lanczos 向量是否是同一概念(是否相同)？
Lanczos 向量与Ritz 向量是等价的，但物理概念不同（Lanczos 方法计算的是自振振型，Lanczos 向量则是强迫振动时的振型）。
5、子空间迭代法与Rayleigh-Ritz 法和矩阵迭代法的关系，与Rayleigh-Ritz 法和矩阵迭代法相比，子空间迭代法有什么优点？
子空间迭代法是Rayleigh-Ritz 法和矩阵迭代法的综合应用，通过矩阵迭代法修正子空间，再通过Rayleigh-Ritz 法得到子空间的频率和振型。
与Rayleigh-Ritz 法相比的优点：子空间得到修正，精度更高。
与矩阵迭代法相比的优点：迭代次数减少，且一次可计算出多个频率和振型。
6、分析用Rayleigh-Ritz 法给出的连续介质模型的特征方程和一般广义坐标法给出的体系的运动方程，讨论两者之间的异同。
不同点：Rayleigh-Ritz 法是有限自由度（一种特殊的广义坐标法），一般广义坐标法是无限自由度。
相同点：表达式类似，当Rayleigh-Ritz 法的自由度数趋于无穷时，两者等价。
7、当研究的对象为非均匀连续梁时，如何运用Rayleigh 方法，或Rayleigh-Ritz 法计算非均匀梁的自振频率和振型？如何完成相应的理论和计算公式的推导？ 只需在计算Rayleigh 熵时将[K]与[M]换成相应的积分形式：
()()()()()2
02
''l
l EI x x dx
m x x dx
ψρψψ????=
?
????
?
8、在结构自振频率的Rayleigh 方法近似分析中，选取假设振型的原则是什么？ 尽量接近所关心的真实振型。
9、用Rayleigh-Ritz 法可以得到结构的一组自振频率和振型，在使用这一方法时，如何选取一组假设振型？
振型是独立的，且满足相应的边界条件。 10、为什么矩阵逆迭代法会使结果收敛于低阶自振频率和振型，而矩阵正迭代法则使结果收
敛于高阶自振频率和振型？
矩阵正迭代法等价于幂法（计算最大特征值），矩阵逆迭代法等价于反幂法（计算最小特征值）。
11、对于多自由度体系，是否根据任意选取的两阶假设振型，采用近似分析方法得到的两阶频率一定是结构前一、二阶频率的近似值？
不能，只能得到前一、二阶频率的上界。
12、能否用荷载相关的Ritz向量做为一组假设振型，由Rayleigh-Ritz法得到结构的一组自振频率和振型？如果可以，在计算Ritz向量时，计算的每一步都要求了各振型满足正交条件，为什么还可以进行迭代？
可以。在计算Ritz向量时，只要求Ritz向量关于质量阵正交，而没有要求它们关于刚度阵正交，故仍有迭代的空间。
思考题十四
1、什么是P-Δ效应？P-Δ效应的本质是什么？如何处理高层结构动力反应问题分析中的P-Δ效应？
P-Δ效应：轴力与横向荷载导致的位移相耦合而引起附加弯矩的现象。
本质：几何非线性（相对大转动）。
处理方法：引入几何刚度矩阵。
2、什么是结构健康监测？什么是结构检测？两者的研究目标和采用的研究方法有什么异同？
结构健康监测：对工程结构实施的损伤监测和识别。
结构检测：以各种试验测试手段对结构或构件的工作性能作出评价。
相同点：研究方法上均通过静力或动力试验得到数据，进行分析。
不同点：前者的目标是确定结构的健康状态，后者的目标是对工程结构的工作性能作出评价。
3、什么是结构地震反应问题中的多维地震动输入问题，这一问题的本质是什么？
结构动力反应分析中所需要考虑的结构不同反应分量之间的耦合影响，即“二阶效应”。
4、什么是结构多点地震动输入问题，如何建立这一问题的运动方程，能否用振型叠加法分析该问题？
将位移分解为相对运动和牵连运动，牵连运动也是与位置有关的量。
对于线性问题可以使用。
5、什么是物理非线性？在结构地震反应分析研究中，大体上有几种层面的结构物理非线性模型？
物理方程的非线性。
有材料层面、单轴应力应变层面、截面层面、构件层面等。
[0729]《结构力学》 1、桁架计算的结点法所选分离体包含几个结点 A. 单个 2、固定铰支座有几个约束反力分量 B. 2个 3、从一个无多余约束的几何不变体系上去除二元体后得到的新体系是 A. 无多余约束的几何不变体系 4、两刚片用三根延长线交于一点的链杆相连组成 A. 瞬变体系 5、定向滑动支座有几个约束反力分量 B. 2个 6、结构的刚度是指 C. 结构抵抗变形的能力 7、桁架计算的截面法所选分离体包含几个结点 B. 最少两个 8、对结构进行强度计算的目的，是为了保证结构 A. 既经济又安全 9、可动铰支座有几个约束反力分量 A. 1个 10、固定支座(固定端)有几个约束反力分量 C. 3个 11、改变荷载值的大小，三铰拱的合理拱轴线不变。 A.√ 12、多余约束是体系中不需要的约束。 B.× 13、复铰是连接三个或三个以上刚片的铰 A.√ 14、结构发生了变形必然会引起位移，结构有位移必然有变形发生。 B.×
15、如果梁的截面刚度是截面位置的函数，则它的位移不能用图乘法计算。 A.√ 16、一根连杆相当于一个约束。 A.√ 17、单铰是联接两个刚片的铰。 A.√ 18、连接四个刚片的复铰相当于四个约束。 B.× 19、虚功原理中的力状态和位移状态都是虚设的。 B.× 20、带拉杆三铰拱中拉杆的拉力等于无拉杆三铰拱的水平推力。 A.√ 21、瞬变体系在很小的荷载作用下会产生很大的内力，所以不能作为结构使用。 A.√ 22、一个无铰封闭框有三个多余约束。 A.√ 23、三铰拱的水平推力不仅与三铰的位置有关，还与拱轴线的形状有关。 B.× 24、三铰拱的主要受力特点是：在竖向荷载作用下产生水平反力。 A.√ 25、两根链杆的约束作用相当于一个单铰。 B.× 26、不能用图乘法求三铰拱的位移。 A.√ 27、零杆不受力，所以它是桁架中不需要的杆，可以撤除。 B.× 28、用图乘法可以求等刚度直杆体系的位移。 A.√ 29、连接四个刚片的复铰相当于四个约束。
结构动力学大作业 对于如下结构，是研究质量块的质量变化和在简支梁上位置的变化对整个系统模态的影响。 1 以上为一个简支梁结构。集中质量块放于梁上，质量块距简支梁的左端点距离为L. 将该简支梁简化为欧拉伯努利梁，并离散为N 个单元。每个单元有两个节点，四个自由度。 单元的节点位移可表示为： ]1122,,,e v v δθθ?=? 则单元内一点的挠度可计作： 带入边界条件： 1 3 32210)(x a x a x a a x v +++=0 1)0(a v x v ===3 322102)(L a L a L a a v L x v +++===1 10 d d a x v x ===θ2 321232d d L a L a a x v L x ++===θ1 0v a =
[]12 3 4N N N N N = 建立了单元位移模式后，其动能势能均可用节点位移表示。单元的动能为： 00111()222 l l T T T ke e e e e y E dx q N Ndxq q mq t ρρ?===??? 其中m 为单元质量阵，并有： l T m N Ndx ρ=? 带入公式后积分可得： 222215622541322413354 1315622420133224l l l l l l l m l l l l l l ρ-?? ??-??= ?? -?? ---? ? 单元势能可表示为 22 200 11()()22 2 T l l T T e pe e e e q y E EI dx EI N N dxq q Kq x ?''''== =??? 其中K 为单元刚度矩阵，并有 ()l T K EI N N dx ''''=? 2 23 2212 612664621261266264l l l l l l EI k l l l l l l l -????-??=??---??-?? 以上为单元类型矩阵，通过定义全局位移矩阵，可以得到系统刚度矩阵和系统质量矩 1 1θ=a )2(1)(3211222θθ+--=L v v L a )(1)(22122133θθ++-= L v v L a 1232133222231)(θ???? ??+-+???? ??+-=L x L x x v L x L x x v 2 2232332223θ??? ? ??-+???? ??-+L x L x v L x L x 2 4231211)()()()()(θθx N v x N x N v x N x v +++=
华中科技大学土木工程与力学学院 《结构动力学》考试卷（B卷、闭卷） 2013～2014学年度第一学期成绩 学号专业班级姓名 一、简答题（每题5分、共25分） 1、刚度法和柔度法所建立的体系运动方程间有何联系？各在什么情况下使用方便？ 答：从位移协调的角度建立振动方程的方法为柔度法。从力系平衡的角度建立的振动方程的方法为刚度法。这两种方法在本质上是一致的，有着相同的前提条件。在便于求出刚度系数的体系中用刚度法方便。同理，在便于求出柔度系数的体系中用柔度法方便。在超静定结构中，一般用刚度法方便，静定结构中用柔度法方便。 2、什么叫动力系数，动力系数大小与哪些因素有关？单自由度体系位移动力系数与内力动力系数是否一样？ 答：动力系数是指最大动位移[y(t)]max与最大静位移yst的比值，其与体系的自振频率和荷载频率θ有关。当单自由度体系中的荷载作用在质量处才有位移动力系数与内力动力系数一样的结果。 3、什么叫临界阻尼？怎样量测体系振动过程中的阻尼比？若要避开共振应采取何种措施？ 答：当阻尼增大到体系在自由反应中不再引起振动，这时的阻尼称为临界阻尼。根据公式即测出第k次振幅和第k+n次振幅即可测出阻尼比。 措施：○1可改变自振频率，如改变质量、刚度等。○2改变荷载的频率。○3可改变阻尼的大小，使之避开共振。 4、振型正交的物理意义是什么？振型正交有何应用？频率相等的两个主振型互相正交吗？ 答：物理意义：第k主振型的惯性力与第i主振型的位移做的功和第i主振型的惯性力与第k主振型的静位移做的功相等，即功的互等定理。 作用：○1判断主振型的形状特点。○2利用正交关系来确定位移展开公式中的系数。 5、应用能量法求频率时，所设的位移函数应满足什么条件？其计算的第一频率与精确解相比是偏高还是偏低？什么情况下用能量法可得到精确解？ 答：所设位移函数要满足位移边界条件，同时要尽可能与真实情况相符。第一频率与精确解相比偏高。如果所假设的位移形状系数与主振型的刚好一致，则可以得到精确解。
2012学年《结构动力学》作业1 发布日期：3月9日上交日期：3月16日 1．采用牛顿第二定律推导复合摆的 运动方程，该复合摆由一根长L， 单位长度的质量为m的均质棒以 及半径为R质量为M的圆盘组成 （见图1）。 图1：复合摆示意图 2．推导图2中系统的等效弹簧常数。 图2：由弹簧通过刚性连杆支持的系统 3．承受弯曲的悬臂梁是由2个均匀段 组成，如图3所示。求对应于自由 端x＝L处施加垂直力时的等效弹 簧常数。 图3：非均匀梁作为弹簧 4．如图4，比重计质量为0.0115 kg， 用于测定某液体的密度。比重计伸 出液面部分的玻璃管直径为0.8 cm，液体比重为1.02 （即是水的 密度的1.02倍）。现将比重计轻轻 地向下按一下，比重计将作上下自 由振动，求振动周期。 图4 5．如下图所示，重量为P的小车从斜面上高h处滑下，与缓冲弹簧相撞后，随同弹簧一起做自由振动。弹簧刚度为K，斜面倾角为 ，小车与斜面间摩擦不计。求小车的振动周期和振幅。（注意：振幅为相对于弹簧静平衡位置） 6．教材习题2-1 7．教材习题2-2
8． 如教材图2-7所示单自由度系统，假设m ＝1kg ，K ＝100N/m ，初始条件x(0)=0.1m ， 0)0(=x ，a) 绘制 c ＝1 N ·s/m ，5N ·s/m ，10N ·s/m 条件下，t ＝0~10s 的响应；b ）绘制 c ＝20 N ·s/m ，30N ·s/m ，40N ·s/m 条件下, t ＝0~10s 的响应。要求用Matlab 编程计算并绘图。对结果进行分析。 9． 教材习题2-4 10． 教材习题2-5 11． 一个有粘性阻尼的弹簧质量系统，作自由振动时测得振动周期为1.8s ，相邻两振幅之比 为4.2:1。求此系统的固有频率。 12． 列出下图系统的振动微分方程。已知m ＝98 N ，K ＝9800 N/m ，r ＝9800 N s/m ，a ＝L/3, b=2L/3。（1）求系统振动时的频率（注意：不是固有频率），并与无阻尼时的固有频率作比较；（2）求系统振动时振幅的对数衰减率。 13． 一质量弹簧系统的质量块重W ＝19.6 kN ，弹簧刚度系数K ＝48.02 kN/m ，今需在此系统 中配置一粘性阻尼，使系统的相对阻尼系数1.0=?，问阻尼器的粘性阻尼系数c 应为多少？系统自由振动时的频率为多少？
结构风工程课后思考题参考答案 二、大气边界层风特性 1 对地表粗糙度的两种描述方式：指数律和对数律（将公式写上）。 2 非标准地貌下的风速换算原则（P）和方法（P公式）。1514 3 脉动风的生成： 近地风在流动过程中由于受到地表因素的干扰，产生大小不同的涡旋，这些涡旋的迭加作用在宏观上表现为速度的随机脉动。在接近地面时，由于受到地表阻力的影响，导致风速减慢并逐步发展为混乱无规则的湍流。 脉动风的能量及耗散机制：而湍流运动可以看做是能量由低频脉动向高频脉动过渡，并最终被流体粘性所耗散的过程。在低频区漩涡尺度较大，向中频区（惯性子区）、高频区（耗散区）漩涡尺度逐渐减小，小尺度涡吸收由惯性子区传递过来的能量，能量最终被流体粘性所耗散。 4 Davenport谱的特点：先写出公式 通过不同水平脉动风速谱的比较： （1）D谱不随高度变化，而其他谱（如Kaimal谱、Solari谱、Karman谱）则考虑了近地湍流随高度变化的特点；（D谱不随高度变化，在高频区符合-5/3律，没有考虑近地湍流随高度变化的特点；） （2）D谱的谱值比其它谱值偏大，会高估结构的动力反应，计算结果偏于保守。（3）S(0)=0,意味着L=0，与实际不符。uu5 湍流度随高度及地面粗糙程度的变化规律：随地面粗糙度的增大而增大，随高度的增加而减小。 积分尺度随高度及地面粗糙程度的变化规律：大量观测结果表明，大气边界层中的湍流积分尺度是地面粗糙度的减函数，而且随着高度的增加而增加。 功率谱随高度及地面粗糙程度的变化规律：随着高度增大和粗糙度的减小，能量在频率上的分布趋于集中，谱形显得高瘦；随着高度减小和粗糙度的增大，能量在频率上的分布趋于分散，谱形显得扁平。 相干函数随高度及地面粗糙程度的变化规律：随地面粗糙度的增大而减小，随高度的增加而增大。 6 阵风因子与峰值因子的区别：阵风因子G=U'/U,是最大风速与平均风速的比/ σ是最大脉动风速与脉动风速均方根的比值。g=u 值；峰值因子umax联系：二者可以相互换算：G=(U'+gσ)/U'=1+gσ/U'=1+gI。Uuu 三、钝体空气动力学理论 1 钝体绕流的主要特征有： ）粘性效应：气体粘性随温度升高而增大，液体粘性随温度升高而减小。1（（2）边界层的形成：由于粘性效应，使靠近物体表面的空气流动速度减慢，形 成气流速度从表面等于零逐渐增大到与外层气流速度相等，形成近壁面流动现象。 （3）边界层分离：如果边界层内的流体微粒速度因惯性力减小到使靠近表面的气流倒流，便出现了边界层分离。 （4）再附：在一定条件下，自建筑物前缘分离的边界层会偶然再附到建筑物表面，这时附面层下会形成不通气的空腔，即分离泡。每隔一段时间分离泡破裂产生较大的风吸值，产生一个风压脉冲。 （5）钝体尾流：对于细长钝体，漩涡脱落是在其两侧交替形成的。漩涡脱落时导致建筑物出现横向振动的主要原因。
结构动力学作业 姓名： 学号：
目录 1.力插值法 (1) 1.1分段常数插值法 (1) 1.2分段线性插值法 (4) 2.加速度插值法 (7) 2.1常加速度法 (7) 2.2线加速度法 (9) 附录 (12) 分段常数插值法源程序 (12) 分段线性插值法源程序 (12) 常加速度法源程序 (13) 线加速度法源程序 (13)
1.力插值法 力插值法对结构的外荷载进行插值，分为分段常数插值法和分段线性插值法，这两种方法均适用于线性结构的动力反应计算。 1.1分段常数插值法 图1-1为一个单自由度无阻尼系统，结构的刚度为k ，质量为m ，位移为y (t )，施加的外力为P (t )。图1-2为矩形脉冲荷载的示意图，图中t d 表示作用的时间，P 0表示脉冲荷载的大小。 图1-1 单自由度无阻尼系统示意图 图1-2 矩形脉冲荷载示意图 对于一个满足静止初始条件的无阻尼单自由度体系来说，当施加一个t d 时间的矩形脉冲荷载，此时结构在t d 时间内的位移反应可以用杜哈梅积分得到： 0()sin ()2 (1cos )(1cos ) (0) t st st d P y t t d m t y t y t t T ωττω πω=-=-=-≤≤? (1-1) 如果结构本身有初始的位移和速度，那么叠加上结构自由振动的部分，结构的位移反应为： 02()cos sin (1cos ) (0 )st d y t y t y t t y t t T πωωω =+ +-≤≤ (1-2)
图1-3 分段常数插值法微段示意图 对于施加于结构任意大小的力，将其划分为Δt 的微段，每一段的荷载都为一个常数(每段相当于一个矩形的脉冲荷载)，如图1-3所示，则将每一段的位移和速度写成增量的形式为： 1cos t sin t (1cos t)i i i i y P y y k ωωωω +=?+ ?+-? (1-3) i+1/sin t cos t sin t i i i y P y y k ωωωωω =-?+ ?+ ? (1-4) 程序流程图如下
Advanced Structural Dynamics Project The dynamic response and stability analysis of the beam under vertical excitation Instructor：Dr. Li Wei Name： Student ID：
1.Problem description and thepurpose of the project 1.1 calculation model An Eular beam subjected to an axial force. Please build thedifferential equation of motion and use a proper difference method to solve this differentialequation. Study the dynamic stability of the beam related to the frequency and amplitude of the force. As shown in the Fig 1.1. Fig1.1 1.2 purpose and process arrangement a.learninghow to create mathematical model of thecontinuous system and select proper calculation method to solve it. b.learning how to build beam vibration equation and solve Mathieu equation. https://www.doczj.com/doc/5e1501457.html,ing Floquet theory to judgevibration system’s stability and analyze the relationship among the frequency and amplitude of the force and dynamic response. This project will introduce the establishment of the mathematical model of the continuous system in section 2, the movement equation and the numerical solution of using MATLAB in section 3,Applying Floquent theory to study the dynamic stability of the beam related to the frequency and amplitude of the force in section 4. In the last of the project, we get some conclusions in section 5.
西南交《结构力学E》离线作业 一、单项选择题(只有一个选项正确，共13道小题) 1. 瞬变体系在一般荷载作用下( Ｃ) (A) 产生很小的内力 (B) 不产生内力 (C) 产生很大的内力 (D) 不存在静力解答 2. 图示体系为：Ｂ (A) 几何不变无多余约束 (B) 几何不变有多余约束； (C) 常变体系； (D) 瞬变体系。 3. 图示某结构中的AB杆的隔离体受力图，则其弯矩图的形状为( Ｂ)
(A) 图a (B) 图b (C) 图c (D) 图d 4. 图示结构：Ｂ (A) ABC段有内力； (B) ABC段无内力； (C) CDE段无内力； (D) 全梁无内力。 5. 常变体系在一般荷载作用下（Ｄ） (A) 产生很小的内力 (B) 不产生内力 (C) 产生很大的内力 (D) 不存在静力解答 6. 图示体系的几何组成为Ｄ
(A) 几何不变，无多余联系； (B) 几何不变，有多余联系； (C) 瞬变； (D) 常变。 7. 在弯矩图的拐折处作用的外力是（Ｂ）。 (A) 轴向外力 (B) 横向集中力 (C) 集中力偶 (D) 无外力 8. 对于图示结构，下面哪个结论是正确的。（Ｂ） (A) 该结构为桁架结构； (B) 该结构是组合结构，其中只有57杆是受拉或受压杆（二力杆）； (C) 只有杆34的内力有弯矩； (D) 除杆123外，其余各杆均为二力杆。
9. 在径向均布荷载作用下，三铰拱的合理轴线为：( Ａ) (A) 圆弧线； (B) 抛物线； (C) 悬链线； (D) 正弦曲线。 ： 10. 如图示各结构弯矩图的形状正确的是( Ｂ) (A) 如图a (B) 如图b (C) 如图c (D) 如图d 11. 静定结构在支座移动时，会产生：( Ｃ) (A) 内力； (B) 应力； (C) 刚体位移； (D) 变形。 12. 图示桁架，各杆EA为常数，除支座链杆外，零杆数为：(Ａ )
第七章影响线 一判断题 1. 图示梁AB与A0B0，其截面C与C0弯矩影响线和剪力影响线完全相同。（X） 题1图题2图 2. 图示结构Q E影响线的AC段纵标不为零。（X） 3. 图示梁K截面的M K影响线、Q K影响线形状如图a、b所示。 4. 图示梁的M C影响线、Q C影响线形状如图a、b所示。 5. 图示梁的M C影响线、M B影响线形状如图a、b所示。 6. 图示结构M B影响线的AB段纵标为零。 7. 图示梁跨中C截面弯矩影响线的物理意义是荷载P=1作用在截面C的弯矩图形。（X） 8. 用静力法作静定结构某量值的影响线与用机动法作该结构同一量值的影响线是不等价 的。（X） 9. 求某量值影响线方程的方法，与恒载作用下计算该量值的方法在原理上是相同的。（√） 10. 影响线是用于解决活载作用下结构的计算问题，它不能用于恒载作用下的计算。（X） 11. 移动荷载是指大小，指向不变，作用位置不断变化的荷载，所以不是静力荷载。（X） 12. 用静力法作影响线，影响线方程中的变量x代表截面位置的横坐标。（X） 13. 表示单位移动荷载作用下某指定截面的内力变化规律的图形称为内力影响线。（√） 14. 简支梁跨中截面弯矩的影响线与跨中有集中力P时的M图相同。（X） 15. 简支梁跨中C截面剪力影响线在C截面处有突变。（√） 16. 绝对最大弯矩是移动荷载下梁的各截面上最大的弯矩。（√） 17. 静定结构及超静定结构的内力影响线都是由直线组成。（X） 18. 图示结构Q C影响线的CD段为斜直线。 19. 图示结构K断面的剪力影响线如图b所示。（√） 题19图 20. 用机动法作得图a所示Q B左结构影响线如图b。 题20图题21图 21. 图示结构a杆的内力影响线如图b所示 22. 荷载处于某一最不利位置时，按梁内各截面得弯矩值竖标画出得图形，称为简支梁的弯
一．是非题（将判断结果填入括弧：以O 表示正确，X 表示错误）（本大题分4小题，共 11分） 1 . (本小题 3分) 图示结构中DE 杆的轴力F NDE =F P /3。（ ）. 2 . (本小题 4分) 用力法解超静定结构时，只能采用多余约束力作为基本未知量。 （ ） 3 . (本小题 2分) 力矩分配中的传递系数等于传递弯矩与分配弯矩之比，它与外因无关。（ ） 4 . (本小题 2分) 用位移法解超静定结构时，基本结构超静定次数一定比原结构高。 （ ） 二．选择题（将选中答案的字母填入括弧内）（本大题分5小题，共21分） 1 （本小题6分） 图示结构EI=常数，截面A 右侧的弯矩为：（ ） A ．2/M ； B ．M ； C ．0； D. )2/(EI M 。 2. （本小题4分） 图示桁架下弦承载，下面画出的杆件内力影响线，此杆件是：（ ） Ａ.ch; B.ci; C.dj; D.cj. 2
3. (本小题 4分) 图a 结构的最后弯矩图为： A. 图b; B. 图c; C. 图d; D.都不对。（ ） ( a) (b) (c) (d) 4. (本小题 4分) 用图乘法求位移的必要条件之一是： A.单位荷载下的弯矩图为一直线； B.结构可分为等截面直杆段； C.所有杆件EI 为常数且相同； D.结构必须是静定的。 ( ) 5. （本小题3分） 图示梁A 点的竖向位移为（向下为正）：( ) Ａ．F P l 3 /(24EI); B. F P l 3 /(!6EI); C. 5F P l 3 /(96EI); D. 5F P l 3 /(48EI). 三（本大题 5分）对图示体系进行几何组成分析。 F P =1
结构力学大作业——五层三跨框架结构内力计算 专业班级：土木工程XXXX班 姓名 XXXXX 学号：XXXXX 指导教师：XX
目录 一、题目 (3) 二、任务 (5) 三、结构的基本数据 (5) 1.构件尺寸： (5) 2.荷载： (5) 3.材料性质： (5) 四、水平荷载作用下的计算 (5) 1.反弯点法 (6) 2.D值法 (8) 3.求解器法 (12) 五、竖直荷载作用下的计算 (15) 1.分层法 (16) 2.求解器法 (21) 六、感想 (24)
二、题目 结构（一） 1、计算简图如图1所示。 4 . 2 m 3 . 6 m 3 . 6 m 3 . 6 m 3 . 6 m 图1
’ 图2 q’ 图3
二、任务 1、计算多层多跨框架结构在荷载作用下的内力，画出内力图。 2、计算方法： （1） 水平荷载： D 值法、反弯点法、求解器，计算水平荷载作用下的框架 弯矩； （2） 竖向荷载：迭代法、分层法、求解器，计算竖向荷载作用下框架弯矩。 3、对各种方法的计算结果进行对比，分析近似法的误差。 4、把计算过程写成计算书的形式。 三、结构的基本数据 E h =3.0×107kN/m 2 柱尺寸：400×400，梁尺寸（边梁）：250×600，（中间梁）300×400 竖向荷载：q '=17kN/m 水平荷载：F P '=15kN 构件线刚度：)12 (,3 bh I l EI i == 柱子：43-3 10133.212 400400m I ?=?= 柱 第一层：m kN i ?=???= -152382.410133.2100.33 71 第二--五层：m kN i ?=???= -177786.310133.2100.33 72 梁： 边梁：43-3105.412 600250m I ?=?=边梁 m kN i ?=???=-225006105.4100.3373 中间梁：43-3106.112 400300m I ?=?=中间梁 m kN i ?=???=-228571 .2106.1100.3374 四、水平荷载作用下的计算 水平荷载： F P =16kN ，F p '=15kN
. ... . 院（系） 建筑工程系 学号 三 明 学院 姓名 . 密封 线 内 不 要 答 题 密封……………………………………………………………………………………………………… ……………………………………………结构力学试题答案汇总 一、选择题(每小题3分,共18分) 1. 图 示 体 系 的 几 何 组 成 为 ： （ A ） A. 几 何 不 变 ， 无 多 余 联 系 ； B. 几 何 不 变 ， 有 多 余 联 系 ； C. 瞬 变 ； D. 常 变 。 （第1题） （第4题） 2. 静 定 结 构 在 支 座 移 动 时 ， 会 产 生 ： （ C ） A. 力 ； B. 应 力 ； C. 刚 体 位 移 ； D. 变 形 。 3. 在 径 向 均 布 荷 载 作 用 下 ， 三 铰 拱 的 合 理 轴 线 为： （ B ） A ．圆 弧 线 ； B ．抛 物 线 ； C ．悬 链 线 ； D ．正 弦 曲 线 。 4. 图 示 桁 架 的 零 杆 数 目 为 ： （ D ） A. 6； B. 7； C. 8； D. 9。 5. 图 a 结 构 的 最 后 弯 矩 图 为 ： （ A ） A ．图 b ； B ．图 c ； C ．图 d ； D ．都不 对 。 6. 力 法 方 程 是 沿 基 本 未 知 量 方 向 的 ： （ C ） A ．力 的 平 衡 方 程 ； B ．位 移 为 零 方 程 ； C ．位 移 协 调 方 程 ； D ．力 的 平 衡 及 位 移 为 零 方 程 。
. ... . 二、填空题（每题3分,共9分） 1.从 几 何 组 成 上 讲 ， 静 定 和 超 静 定 结 构 都 是___几何不变____ 体 系 ， 前 者___无__多 余 约 束 而 后 者____有___多 余 约 束 。 2. 图 b 是 图 a 结 构 ___B__ 截 面 的 __剪力__ 影 响 线 。 3. 图 示 结 构 AB 杆 B 端 的 转 动 刚 度 为 ___i___, 分 配 系 数 为 ____1/8 ____, 传 递 系 数 为 ___-1__。 三、简答题（每题5分，共10分） 1.静定结构内力分析情况与杆件截面的几何性质、材料物理性质是否相关？ 为什么？ 答：因为静定结构内力可仅由平衡方程求得，因此与杆件截面的几何性质无关， 与材料物理性质也无关。 2.影响线横坐标和纵坐标的物理意义是什么？ 答：横坐标是单位移动荷载作用位置，纵坐标是单位移动荷载作用在此位置时物 理量的影响系数值。 四、计算分析题，写出主要解题步骤(4小题,共63分) 1.作图示体系的几何组成分析（说明理由），并求指定杆1和2的轴力。（本题16分） （本题16分）1.因为w=0 所以本体系为无多约束的几何不变体系。（4分） F N1=- F P （6分）； F N2=P F 3 10（6分）。 2.作 图 示 结 构 的 M 图 。（本题15分）
1、结构的刚度是指 1. C. 结构抵抗变形的能力 2、 图7中图A～图所示结构均可作为图7（a）所示结构的力法基本结构，使得力法计算最为简便的 C 3、图5所示梁受外力偶作用，其正确的弯矩图形状应为（）C 4、对结构进行强度计算的目的，是为了保证结构 1. A. 既经济又安全 5、改变荷载值的大小，三铰拱的合理拱轴线不变。 1. A.√ 6、多余约束是体系中不需要的约束。 1. B.×
7、结构发生了变形必然会引起位移，结构有位移必然有变形发生。 1. B.× 8、如果梁的截面刚度是截面位置的函数，则它的位移不能用图乘法计算。 1. A.√ 9、一根连杆相当于一个约束。 1. A.√ 10、单铰是联接两个刚片的铰。 1. A.√ 11、虚功原理中的力状态和位移状态都是虚设的。 1. B.× 12、带拉杆三铰拱中拉杆的拉力等于无拉杆三铰拱的水平推力。 1. A.√ 13、瞬变体系在很小的荷载作用下会产生很大的内力，所以不能作为结构使用。 1. A.√ 14、虚位移原理中的虚功方程等价于静力平衡方程，虚力原理中虚功方程等价于变形协调方程。 1. A.√ 15、体系的多余约束对体系的计算自由度、自由度及受力状态都没有影响，故称多余约束。 1. B.× 16、力矩分配中的传递系数等于传递弯矩与分配弯矩之比，它与外因无关。 1. A.√ 17、当上部体系只用不交于一点也不全平行的三根链杆与大地相连时，只需分析上部体系的几何组成，就能确1. A.√ 18、用力法计算超静定结构时，其基本未知量是未知结点位移。
B.× 19、静定结构在非荷载外因（支座移动、温度改变、制造误差）作用下，不产生内力，但产生位移。 1. A.√ 20、力法和位移法既能用于求超静定结构的内力，又能用于求静定结构的内力。（） 1. B.× 21、静定结构在非荷载外因（支座移动、温度改变、制造误差）作用下，不产生内力，但产生位移。（）1. A.√ 22、位移法和力矩分配法只能用于求超静定结构的内力，不能用于求静定结构的内力。( ） 1. B.× 23、 图2所示体系是一个静定结构。（） 1. B.× 24、力矩分配法中的分配系数、传递系数与外来因素（荷载、温度变化等）有关。 1. B.× 25、三铰拱的水平推力不仅与三铰的位置有关，还与拱轴线的形状有关。 1. B.× 26、三铰拱的主要受力特点是：在竖向荷载作用下产生水平反力。 1. A.√ 27、两根链杆的约束作用相当于一个单铰。 B.× 28、不能用图乘法求三铰拱的位移。
目录 一、结构特性矩阵 1.1框架设计 (2) 1.2截面尺寸 (2) 1.3动力自由度 (2) 1.4结构的一致质量矩阵 (3) 1.5结构的一致刚度矩阵 (13) 二、频率与振型 2.1简化的质量矩阵 (25) 2.2简化的刚度矩阵 (25) 2.3行列式法求频率与振型 (27) 2.4Stodola法求频率与振型 (27) 三、时程分析 3.1框架资料 (31) 3.2地震波波形图 (31) 3.2瑞利阻尼 (32) 3.4操作步骤 (33) 3.5各楼层位移时程反应图 (37)
一、结构特性矩阵 1.1框架设计 框架平面图如图1所示，跨度均为6.0m,层高均为3.6m，混凝土采用C30。 图1 框架平面图 1.2截面尺寸 梁均为300mm600mm ? ?,柱均为500mm500mm 1.3动力自由度 框架结构可以理想化为在节点处相互连接的梁柱单元的集合。设梁、柱的轴向变形均忽略不计，只考虑横向平面位移，则该框架有3个平动自由度和12个角自由度，共15个自由度，并对梁柱单元分别编号，如图2所示： 图2 单元编号及自由度
将结构分成在有限个节点处相互连接的○1~○21个离散单元体系，通过计算各个单元的一致质量矩阵、一致刚度矩阵，并将相关的单元叠加求得整个单元结构的一致质量矩阵、一致刚度矩阵。 1.4结构的一致质量矩阵 在节点位移作用下框架梁和柱上所引起的变形形状采用三次Hermite 多项式，因此均布质量梁的一致质量矩阵为： ??? ???? ???????4 3 2 1 I I I I f f f f =420L m ?? ? ?? ???????------222 2432213341322221315654132254156 L L L L L L L L L L L L ???? ????????? ????? (4) .. 3 2 1 v v v v 梁：m =250060.030.0??=450kg/m, L=6m;
第一章 单自由度系统 1.1 总结求单自由度系统固有频率的方法和步骤。 单自由度系统固有频率求法有：牛顿第二定律法、动量距定理法、拉格朗日方程法和能量守恒定理法。 1、 牛顿第二定律法 适用围：所有的单自由度系统的振动。 解题步骤：（1） 对系统进行受力分析,得到系统所受的合力； （2） 利用牛顿第二定律∑=F x m ,得到系统的运动微分方程； （3） 求解该方程所对应的特征方程的特征根，得到该系统的固有频率。 2、 动量距定理法 适用围：绕定轴转动的单自由度系统的振动。 解题步骤：（1） 对系统进行受力分析和动量距分析； （2） 利用动量距定理J ∑=M θ ,得到系统的运动微分方程； （3） 求解该方程所对应的特征方程的特征根，得到该系统的固有频率。 3、 拉格朗日方程法： 适用围：所有的单自由度系统的振动。 解题步骤：（1）设系统的广义坐标为θ，写出系统对于坐标θ的动能T 和势能U 的表达式；进一步写求出拉格朗日函数的表达式：L=T-U ； （2）由格朗日方程 θ θ??-???L L dt )( =0，得到系统的运动微分方程； （3） 求解该方程所对应的特征方程的特征根，得到该系统的固有频率。 4、 能量守恒定理法 适用围：所有无阻尼的单自由度保守系统的振动。
解题步骤：（1）对系统进行运动分析、选广义坐标、写出在该坐标下系统的动能T 和势能U 的表达式；进一步写出机械能守恒定理的表达式 T+U=Const （2）将能量守恒定理T+U=Const 对时间求导得零,即0) (=+dt U T d ，进一步得到系统的运动微分方程； （3） 求解该方程所对应的特征方程的特征根，得到该系统的固有频率。 1.2 叙述用衰减法求单自由度系统阻尼比的方法和步骤。 用衰减法求单自由度系统阻尼比的方法有两个：衰减曲线法和共振法。 方法一：衰减曲线法。 求解步骤：（1）利用试验测得单自由度系统的衰减振动曲线，并测得周期和相邻波峰和波谷的幅值i A 、1+i A 。 （2）由对数衰减率定义 )ln( 1 +=i i A A δ， 进一步推导有 2 12ζ πζδ-= ， 因为ζ较小， 所以有 π δζ2= 。 方法二：共振法求单自由度系统的阻尼比。 （1）通过实验，绘出系统的幅频曲线， 如下图：
结构动力学大作业 问题描述 《建筑结构抗震设计》高振世P247 该建筑为一幢六层现浇钢筋混泥土框架房屋，屋顶有局部突出的楼梯间和水箱间。混泥土强度等级：梁为C20，柱为C25。混凝土密度为2500kg/m3 本题目将对该梁柱结构的框架房屋进行模态分析，求解出该结构的前8阶固有频率及其对应的模态振型。框架的平、剖面见图1，图2。构件尺寸参见表1、表2。 其材料为混凝土，相关参数为：杨氏模量C20为2.55e10N/m2，C25为2.8e10 N/m2。 图 1 平面视图 图 2 剖面图 表 1 梁的几何尺寸 部位断面 b×h (m×m) 跨度L （m） 屋顶梁0.25×0.60 5.7 楼层梁0.25×0.65 5.7 走道凉0.25×0.40 2.1 表 2 柱的几何尺寸 层次柱高（m）断面(m×m) 1 4 0.50×0.50 2，3，4，5，6 3.6 0.45×0.45
建模 利用有限元商业软件ANSYS8.0建模和有限元分析。 1单元类型（Element Type ）: Beam4 图 3 单元BEAM4的特征 2材料模型: 本题目中所有材料都假定是各向同性线弹性体。 表 3 MARERIAL MODEL No. 梁/柱混凝土标号弹性模量EX N/m2 泊松比PXY 密度DENS kg/ m3 1 柱C25 2.8e10 0. 2 2500 2 梁C20 2.55e10 0.2 2500
3单元实常数： 表 4 REAL CONSTANT No. 部位 横截面面积 AREA（m2） Z轴惯性矩 IZZ（m4） Y轴惯性矩 IYY（m4） 高度 TKZ（m） 宽度 TKY（m） 1 底层柱0.25 5.208e-3 5.208e-3 0.5 0.5 2 其余柱0.2025 3.417e- 3 3.417e-3 0.45 0.45 3 走道梁0.1 1.333e-3 5.208e- 4 0.4 0.25 4 楼层梁0.162 5 5.721e-3 8.464e-4 0.65 0.25 5 屋顶梁0.15 0.45e-3 7.8125e-4 0. 6 0.25 4几何模型，网格划分，施加边界条件： 柱划为个6单元，走道梁划为3个单元，楼层梁划为5个单元，合计3029个单元 图 4
结 构 动 力 学 大 作 业 姓名： 学号：
习题1 用缩法减进行瞬态结构动力学分析以确定对有限上升时间得恒定力的动力学响应。实际结构是一根钢梁支撑着集中质量并承受一个动态荷载。 钢梁长L ，支撑着一个集中质量M 。这根梁承受着一个上升时间为t τ，最大值为F1的动态荷载F(t)。梁的质量可以忽略，需确定产生最大位移响应时间max t 及响应max y 。同时要确定梁中的最大弯曲应力bend σ。 已知：材料特性：25x E E MPa =，质量M ＝0.03t ，质量阻尼ALPHAD=8; 几何尺寸：L ＝450mm I=800.64 mm h=18mm; 荷载为：F1＝20N t τ＝0.075s 提示：缩减法需定义主自由度。荷载需三个荷载步（0至加质量，再至0.075s ， 最后至1s ） ANSYS 命令如下： FINISH /CLE$/CONFIG,NRES,2000 /prep7 L=450$H=18 ET,1,BEAM3 ET,2,MASS21,,,4 R,1,1,800.6,18 R,2,30 !MASS21的实常数顺序MASSX, MASSY, MASSZ, IXX, IYY, IZZ MP,EX,1,2E5$MP,NUXY ,1,0.3 N,1,0,0,0 N,2,450/2,0,0 N,3,450,0,0 E,1,2$E,2,3 !创建单元 TYPE,2$REAL,2 E,2 M,2,UY FINISH /SOLU !进入求解层 ANTYPE,TRANS
TRNOPT,REDUC OUTRES,ALL,ALL$DELTIM,0.004 !定义时间积分步长 ALPHAD,8 !质量阻尼为8 D,1,UY$D,3,UX,,,,,UY !节点1Y方向,约束节点3X、Y方向约束 F,2,FY,0 LSWRITE,1 !生成荷载步文件1 TIME,0.075 FDELE,ALL,ALL F,2,FY,20 LSWRITE,2 !生成荷载步文件2 TIME,1 LSWRITE,3 !生成荷载步文件3 LSSOLVE,1,3,1 !求解荷载文件1,2,3 FINISH /SOLU EXPASS,ON$EXPSOL,,,0.10000 !扩展处理 SOLVE FINISH /POST26 NUMV AR,0 FILE,fdy,rdsp !注意，建立的项目名称为fdy,否则超出最大变量数200，结果无效NSOL,2,2,U,Y,NSOL PLV AR,2 !时间位移曲线 PRV AR,2 !得出在0.10000该时间点上跨中位移最大 /POST1 !查看某个时刻的计算结果 SET,FIRST PLDISP,1 !系统在0.10000秒时总变形图 ETABLE,Imoment,SMISC,6 !单元I点弯矩 ETABLE,Jmoment,SMISC,12 !单元J点弯矩 ETABLE,Ishear,SMISC,2 !单元I点剪力 ETABLE,Jshear,SMISC,8 !单元J点剪力 PLLS,IMOMENT,JMOMENT,1,0 !画出弯矩图 PLLS,ISHEAR,JSHEAR,,1,0 !画出剪力图 结果如下; 随着时间位移的大小：
结构动力学大作业（威尔逊- 法） ： 学号： 班级： 专业：
威尔逊-θ法原理及应用 【摘要】在求解单自由度体系振动方程时我们用了常加速度法及线加速度法等数值分析方法。在多自由度体系中，也有类似求解方法，即中心差分法及威尔逊-θ法。实际上后两种方法也能求解单自由度体系振动方程。对于数值方法，有三个重要要求：收敛性、稳定性及精度。本文推导了威尔逊-θ法的公式，并利用MATLAB 编程来研究单自由度体系的动力特性。 【关键词】威尔逊-θ法 冲击荷载 阻尼比 【正文】威尔逊-θ法可以很方便的求解任意荷载作用下单自由度体系振动问题。实际上，当 1.37θ>时，威尔逊-θ法是无条件收敛的。 一、威尔逊-θ法的原理 威尔逊-θ法是线性加速度法的一种拓展（当1θ=时，两者相同），其基本思路和实现方法是求出在时间段[],t t t θ+?时刻的运动，其中1θ≥，然后通过插得到i t t +?时刻的运动（见图 1.1）。 图 1.1 1、公式推导 推导由t 时刻的状态求t t θ+?时刻的状态的递推公式： 对τ积分
{}{}{}{}{}{})(623 2 t t t t t t t y y t y y y y &&&&&&&-?+++=?++θτ θτττ {}{}{}{}{})2(6)(2t t t t t t t y y t y t y y &&&&&+?+?+=?+?+θθθθ {}{}{}{}{}t t t t t t t y y t y y t y &&&&&26 )()(62-?--?=?+?+θθθθ []{}{} {}[]{}{}{}[]{}{}{})223()26)(6( )(2t t t t t t t t t t y t y y t c y y t y t m P P P R &&&&&&?++?++?+?+-+=?+θθθθθ 2、MA TLAB 源程序： clc;clear; K=input('请输入结构刚度k(N/m)'); M=input('请输入质量(kg)'); C=input('请输入阻尼(N*s/m)'); t=sym('t');%产生符号对象t Pt=input('请输入荷载); Tp=input('请输入荷载加载时长(s)'); Tu=input('请输入需要计算的时间长度(s) '); dt=input('请输入积分步长(s)'); Sita=input('请输入θ'); uds=0:dt:Tu;%确定各积分步时刻 pds=0:dt:Tp; Lu=length(uds); Lp=length(pds); if isa(Pt,'sym')%荷载为函数 P=subs(Pt,t,uds); %将荷载在各时间步离散 if Lu>Lp P(Lp+1:Lu)=0; end elseif isnumeric(Pt)%荷载为散点 if Lu<=Lp