钢筋混凝土结构非线性动力反应分析

钢筋混凝土结构非线性动力反应分析
钢筋混凝土结构非线性动力反应分析

1 问题描述 (1)

2 模型建立 (2)

2.1 轴网布置 (2)

2.2 材料定义 (2)

2.3 截面定义 (3)

2.4 楼板定义 (8)

2.5 构件布置............................................................... 错误!未定义书签。

2.6 构件剖分 (10)

2.7 节点束缚指定 (11)

2.8 先单元插入点指定 (11)

2.9 刚域指定 (12)

2.10 整体模型 (12)

2.11 荷载模式定义 (13)

2.12 质量源定义 (14)

2.13 荷载定义 (14)

2.14 荷载工况定义 (15)

2.14.1 模态分析工况 (15)

2.14.1 反应谱工况 (16)

2.14.3 COMBINE工况(计算柱轴力) (17)

2.14.3 时程反应分析工况 (18)

2.15 定义组 (21)

3 振型分解反应谱法结果 (22)

3.1 楼层剪力及水平地震作用 (22)

3.2 内力图 (23)

3.2.1 弯矩图 (23)

3.2.2 剪力图 (24)

3.2.3 轴力图 (25)

4 塑性铰的计算与定义 (26)

4.1 梁铰柱铰的计算 (26)

4.1.1 梁铰的计算 (26)

4.1.2 柱铰的计算 (29)

4.2 铰属性设置 (31)

4.2.1 梁铰属性设置 (31)

4.2.2 柱铰的设置 (32)

4.3 梁柱铰的指定 (35)

5 时程反应分析结果 (36)

5.1 层位移 (36)

5.2 层间位移 (38)

5.3 层间位移角 (38)

5.4 层剪力 (39)

5.5 出铰情况及破坏机制 (40)

1 问题描述

已知结构为一栋四层框架结构。结构布置如图1-1所示,构件尺寸及配筋情况如图1-2。混凝土强度等级为:1~5层采用C40;6、7层采用C30,恒载按实际梁、板、柱实际重量计算,不考虑装饰荷载,活荷载按2KN/m2考虑,不考虑风荷载。对N23波(峰值加速度350.13cm/s2)进行调整,满足八度(0.2g)地震的加速度幅值。计算下列问题:

1)采用振型分解法,计算该结构在小震时的地震作用,并画出地震内力图2)结构在罕遇地震作用下的时程反应性能(包括层位移、层间位移、层间位移角、层剪力、基底剪力及结构的出铰情况和破坏机制)。

图1-1

图1-2 梁柱配筋图

2 模型建立

2.1 轴网布置

根据题目要求的轴网尺寸、结构布置等结构信息,在sap2000的编辑轴网对话框中依次输入数据,建立结构的几何模型。

建模过程中,采用的统一单位制为(N,mm,C),建模过程如下所示:

图2.1.1轴网的编辑

2.2 材料定义

依题目要求,采用快速定义材料,定义混凝土材料与钢筋材料:

混凝土:1-3层采用C40,4层采用C30;

钢筋:HPB300,HRB335两种钢筋。

具体定义过程以C40和HRB335为例:

图2.2.1定义材料

图2.2.2材料属性

2.3 截面定义

这里的截面包括不同尺寸,不同配筋,不同混凝土强度,不同位置的各种截面。根据1、2、4层与3层的混凝土强度等级不同,边梁和主梁的楼板加强作用不同,主梁和次梁的截面尺寸不同,故定义10种梁截面形式,8种柱截面形式,具体如下:

“2FLG1”(二层主梁,C40混凝土,截面300mmx600mm);

“2,4FLG2”(二层、四层边次梁,C40混凝土,截面300mmx300mm);“2,4FLG3”(二层、四层中次梁,C40混凝土,截面300mmx400mm);“3FLG1”(三层主梁,C30混凝土,截面300mmx600mm);

“3FLG2”(三层边次梁,C30混凝土,截面300mmx300mm);

“3FLG3”(三层中次梁,C30混凝土,截面300mmx400mm);

“4,RFLG1”(顶层主梁,C40混凝土,截面300mmx600mm);“RFLG2”(顶层边次梁,C40混凝土,截面300mmx300mm);“RFLG3”(顶层中次梁,C40混凝土,截面300mmx400mm);

“1FLC1”(一层边柱,C40混凝土,截面500mmx500mm);

“1FLC2”(一层中柱,C40混凝土,截面500mmx500mm);

“2FLC1”(二层边柱,C40混凝土,截面500mmx500mm);

“2FLC2”(二层中柱,C40混凝土,截面500mmx500mm);

“3FLC1”(三层边柱,C40混凝土,截面500mmx500mm);

“3FLC2”(三层中柱,C40混凝土,截面500mmx500mm);

“4FLC1”(四层边柱,C40混凝土,截面500mmx500mm);

“4FLC2”(四层中柱,C40混凝土,截面500mmx500mm);

图2.3.1框架截面定义

“1FLC1”为例,说明柱截面的定义过程,如下:

2.3.2框架柱截面定义

图2.3.3框架柱截面布筋

以“2FLG1”为例,说明梁截面的定义过程,如下:

2.3.4框架梁截面定义

图2.3.5框架柱截面布筋

图2.3.6框架梁截面惯性矩放大系数

2.4 楼板定义

由于楼板的跨度相对比较大,板厚为130mm,剪切作用不是很明显,这里选用薄壳模型。由于结构的1,2,4和3层所使用的混凝土标号不同,这里需要定义两种板截面形式:

“S1,2,4”(一,二,四层楼板,C40混凝土,130mm厚)

“S3”(三层楼板,C30混凝土,130mm厚)

图2.4.1面截面定义

下面3层的楼板定义为例:

图2.4.2壳截面定义

图2.4.3混凝土壳参数定义

2.5构件剖分

将梁柱构件以及板布置到结构的对应位置上,如下图:

为计算精度的要求,这里需要对单元进行剖分。对于框架线单元,这里选择自动剖分框架单元。对于楼板,划分为1*2个单元。如下图:

图2.5.1指定自动框架剖分

图2.5.2指定自动面剖分

剖分只是利用原来的面,但是计算精度更高。若不进行剖分,则会出现楼板与主梁跨中变形不协调的情况。一般一个面至少要剖分为四个单元以上结果才会比较精确,综合考虑计算时间随剖分数量成倍增长以及本框架结构的布局和楼板

大小,采取按数目剖分为2个面。

2.6 节点束缚指定

对建筑结构使用刚性隔板,避免了由于较大平面的楼板用膜单元模拟(楼板的刚度模拟为膜的刚度)所产生的数值准确性问题。这在建筑的侧向(水平)动力分析中是十分有用的,因为其可以显著减少所求解特征值问题的计算时间。选中Z轴即表示按此束缚轴束缚节点隔板的运动,从而实现了刚性楼板的假定。依此指定不同楼层的节点束缚。

图2.6.1指定节点束缚

2.7 先单元插入点指定

对已有模型构件指定插入点,实现了

构件的偏心,构件的几何位置并不改变,

模型中的节点位置并没有改变,而力的传

递作用位置发生了变化。SAP2000中默认

的插入点的偏移值为零,而且结构的方位

基准点在楼层高程处,故为了保证梁的表

面与楼层高程平齐,框架梁的控制点要改

为8 Top Center。

2.8 刚域指定

刚域系数表示指定了端部偏移的刚

域部分,在弯曲和剪切变形时刚域程度。

在刚域系数一栏可输入0到1之间的一个

系数,表示刚域的刚性程度从完全柔性到

完全刚性。完全柔性表示刚域的弯曲和剪

切刚度由线性单元的实际刚度确定;完全

刚性即刚域内部没有弯曲和剪切变形。这

里取为0.5。

2.9 整体模型

将视图切换到x-y平面视图,z=0.0;选择所有节点,第一层柱的下端节点指定为固端。

图2.9.1节点约束

图2.9.2整体模型

2.10 荷载模式定义

这里定义两种荷载模式DEAD、LIVE以及屋面活荷载ROOF LIVE。这里定义一个屋面活荷载其目的再于后续的质量源组合中不能组合屋面活荷载。自重系数表示软件将自动计算结构中所有构件的的自重,将自重乘以这个自重系数施加在结构上。只需要定义一个为1的自重系数,即可以考虑结构的全部自重,否则会造成重复考虑结构自重的错误。

图2.10.1定义荷载模式

2.11 质量源定义

为了满足抗震规范重力荷载代表值的相关要求,这里选用“来自荷载“这一项。按抗震规范的5.1.3条规定,自重,附加恒载的系数为1.0,活荷载的系数取为0.5。

图2.11.1定义质量源

2.12 荷载定义

将非结构荷载转化为各个楼层的均布面荷载,以屋顶为例,楼层上的面荷载指定后如下:

图2.12.1定义楼面荷载(6.75kN/m)

2.13 荷载工况定义

2.1

3.1 模态分析工况

计算结构振动的模态,从而获得各个振型的周期。

图2.13.1模态分析工况定义

2.1

3.2 反应谱工况

通过定义反应谱函数定义振型分解反应谱所需的函数,然后通过荷载工况定义反应谱工况,定义时注意比例系数取9.8;所施加的反应谱为根据中国规范要求的场地类别为Ⅱ类,设计地震分组为第二组,特征周期为0.4s,水平地震影响系数最大值为0.16(八度的多遇地震)。

图2.13.2反应谱曲线

图2.13.3反应谱工况定义

2.1

3.3 COMBINE工况(计算柱轴力)

根据抗规要求,结构等效总重力荷载代表值,多自由度体系取总重力荷载代表值的85%,按照等效重力荷载下的设计值计算柱轴力是的荷载或者比例系数为:活载

=

+

.0+

S。

85

=

?

5.0

活载)

恒载

2.1(

恒载595

.0

02

.1

图2.13.4COMBINE工况定义

2.1

3.3 时程反应分析工况

(1)定义时程反应分析函数

第一章 非线性动力学分析方法

第一章非线性动力学分析方法(6学时) 一、教学目标 1、理解动力系统、相空间、稳定性的概念; 2、掌握线性稳定性的分析方法; 3、掌握奇点的分类及判别条件; 4、理解结构稳定性及分支现象; 5、能分析简单动力系统的奇点类型及分支现象。 二、教学重点 1、线性稳定性的分析方法; 2、奇点的判别。 三、教学难点 线性稳定性的分析方法 四、教学方法 讲授并适当运用课件辅助教学 五、教学建议 学习本章内容之前,学生要复习常微分方程的内容。 六、教学过程

本章只介绍一些非常初步的动力学分析方法,但这些方法在应用上是十分有效的。 1.1相空间和稳定性 一、动力系统 在物理学中,首先根据我们面对要解决的问题划定系统,即系统由哪些要素组成。再根据研究对象和研究目的,按一定原则从众多的要素中选出最本质要素作为状态变量。然后再根据一些原理或定律建立控制这些状态变量的微分方程,这些微分方程构成的方程组通常称为动力系统。研究这些微分方程的解及其稳定性以及其他性质的学问称为动力学。 假定一个系统由n 个状态变量1x ,2x ,…n x 来描述。有时,每个状态变量不但是时间t 的函数而且也是空间位置r 的函数。如果状态变量与时空变量都有关,那么控制它们变化的方程组称为偏微分方程组。这里假定状态变量只与时间t 有关,即X i =X i (t),则控制它们的方程组为常微分方程组。 ),,,(2111 n X X X f dt dX ???=λ ),,,(2122 n X X X f dt dX ???=λ (1.1.1) … ),,,(21n n n X X X f dt dX ???=λ 其中λ代表某一控制参数。对于较复杂的问题来说,i f (i =l ,2,…n)一般是{}i X 的非线性函数,这时方程(1.1.1)就称为非线性动力系统。由于{}i f 不明显地依赖时间t ,故称方程组(1.1.1)为自治动力系统。若{}i f 明显地依赖时间t ,则称方程组(1.1.1)为非自治动力系统。非自治动力系统可化为自治动力系统。 对于非自治动力系统,总可以化成自治动力系统。 例如:)cos(t A x x ω=+

动力学分析方法

1动力学分析方法 结构动力学的研究方法可分为分析方法(结构动力分析)和试验方法(结构动力试验)两大类。[7-10] 分析方法的主要任务是建模(modeling),建模的过程是对问题的去粗取精、去伪存真的过程。在结构动力学中,着重研究力学模型(物理模型)和数学模型。建模方法很多,一般可分为正问题建模方法和反问题建模方法。正问题建模方法所建立的模型称为分析模型(或机理模型)。因为在正问题中,对所研究的结构(系统)有足够的了解,这种系统成为白箱系统。我们可以把一个实际系统分为若干个元素或元件(element),对每个元素或元件直接应用力学原理建立方程(如平衡方程、本构方程、汉密尔顿原理等),再考虑几何约束条件综合建立系统的数学模型。如果所取的元素是一无限小的单元,则建立的是连续模型;如果是有限的单元或元件,则建立的是离散模型。这是传统的建模方法,也称为理论建模方法。反问题建模方法适用于对系统了解(称黑箱系统——black box system)或不完全了解(称灰箱系统——grey box system)的情况,它必须对系统进行动力学实验,利用系统的输入(载荷)和输出(响应——response)数据,然后根据一定的准则建立系统的数学模型,这种方法称为试验建模方法,所建立的模型称为统计模型。 在动力平衡方程中,为了方便起见一般将惯性力一项隔离出来,单独列出,因此通常表达式为: u I M&& (2) = - +P 其中M为质量矩阵,通常是一个不随时间改变的产量;I和P是与位移和速度有关的向量,而与对时间的更高阶导数无关。因此系统是一个关于时间二级导数的平衡系统,而阻尼和耗能的影响将在I和P中体现。可以定义: + = (3) I&& C u Ku 如果其中的刚度矩阵K和阻尼矩阵C为常数,系统的求解将是一个线性的问题;否则将需要求解非线性系统。可见线性动力问题的前提是假设I是与节点位移和速度是线性相关的。 将公式(2)代入(1)中,则有

非线性动力学和混沌理论

非线性动力学和混沌理论 非线性动力学 随着科学技术的发展,非线性问题出现在许多学科之中,传统的线性化方法已不能满足解决非线性问题的要求,非线性动力学也就由此产生。 非线性动力学联系到许多学科,如力学、数学、物理学、化学,甚至某些社会科学等。非线性动力学的三个主要方面:分叉、混沌和孤立子。事实上,这不是三个孤立的方面。混沌是一种分叉过程,孤立子有时也可以和同宿轨或异宿轨相联系,同宿轨和异宿轨是分叉研究中的两种主要对象。 经过多年的发展,非线性动力学已发展出了许多分支。如分叉、混沌、孤立子和符号动力学等。然而,不同的分支之间又不是完全孤立的。非线性动力学问题的解析解是很难求出的。因此,直接分析非线性动力学问题解的行为(尤其是长时期行为)成为研究非线性动力学问题的一种必然手段。 混沌理论是谁提出的? 混沌理论,是系统从有序突然变为无序状态的一种演化理论,是对确定性系统中出现的内在“随机过程”形成的途径、机制的研讨。 美国数学家约克与他的研究生李天岩在1975年的论文“周期3则乱七八糟(Chaos)”中首先引入了“混沌”这个名称。 美国气象学家洛伦茨在2O世纪 6O年代初研究天气预报中大气流动问题时,揭示出混沌现象具有不可预言性和对初始条件的极端敏感依赖性这两个基本特点,同时他还发现表面上看起来杂乱无章的混沌,仍然有某种条理性。 1971年法国科学家罗尔和托根斯从数学观点提出纳维-斯托克司方程出现湍流解的机制,揭示了准周期进入湍流的道路,首次揭示了相空间中存在奇异吸引子,这是现代科学最有力的发现之一。 1976年美国生物学家梅在对季节性繁殖的昆虫的年虫口的模拟研究中首次揭示了通过倍周期分岔达到混沌这一途径。 1978年,美国物理学家费根鲍姆重新对梅的虫口模型进行计算机数值实验时,发现了称之为费根鲍姆常数的两个常数。这就引起了数学物理界的广泛关注。 与此同时,曼德尔布罗特用分形几何来描述一大类复杂无规则的几何对象,使奇异吸引子具有分数维,推进了混沌理论的研究。20世纪70年代后期科学家们在许多确定性系统中发现混沌现象。作为一门学科的混沌学目前正处在研讨之中,未形成一个完整的成熟理论。混沌的理论 要弄明白不可预言性如何可以与确定论相调和,可以来看看一个比整个宇宙次要得多的系统——水龙头滴下的水滴。这是一个确定性系统,原则上流入水龙头中的水的流量是平稳、均匀的,水流出时发生的情况完全由流体运动定律规定。但一个简单而有效的实验证明,这一显然确定性的系统可以产生不可预言的行为。这使我们产生某种数学的“横向思维”,它向我们解释了为什么此种怪事是可能的。 假如你很小心地打开水龙头,等上几秒钟,待流速稳定下来,通常会产生一系列规则的水滴,这些水滴以规则的节律、相同的时间间隔落下。很难找到比这更可预言的东西了。但假如你缓缓打开水龙头,使水流量增大,并调节水龙头,使一连串水滴以很不规则的方式滴落,这种滴落方式似乎是随机的。只要做几次实验就会成功。实验时均匀地转动水龙头,别把龙头开大到让水成了不间断的水流,你需要的是中速滴流。如果你调节得合适,就可以在好多分钟内听不出任何明显的模式出现。 1978年,加利福尼亚大学圣克鲁斯分校的一群年青的研究生组成了一个研究动力学系统的小组。他们开始考虑水滴系统的时候,就认识到它并不像表现出来的那样毫无规则。他们用话筒记录水滴的声音,分析每一滴水与下一滴水之间的间隔序列。他们所发现的是短期的可预言性。要是我告诉你3个相继水滴的滴落时刻,你会预言下一滴水何时落下。例如,假如水滴之间最近3个间隔是0.63秒、1.17秒和0.44秒,则你可以肯定下一滴水将在0.82秒后落下这些数只是为了便于说明问题。事实上,如果你精确地知道头3滴水的滴落时刻,你就可以预言系统的全部未来。 那么,拉普拉斯为什么错了? 问题在于,我们永远不能精确地测量系统的初始状态。我们在任何物理系统中所作出的最精确的测量,对大约10位或12位小数来说是正确的。 但拉普拉斯的陈述只有在我们使测量达到无限精度即无限多位小数,当然那是办不到的时才正确。 在拉普拉斯时代,人们就已知道这一测量误差问题,但一般认为,只要作出初始测量,比如小数点后10位,所有相继的预言也将精确到小数点后10位。误差既不消失,也不放大。 不幸的是,误差确实放大,这使我们不能把一系列短期预言串在一起,得到一个长期有效的预言。例如,假设我知道精确到小数点后10位的头3滴水的滴落时刻,那么我可以精确到小数点后9位预言下一滴的滴落时刻,再下一滴精确到8位,以此类推。 误差在每一步将近放大10倍,于是我对进一步的小数位丧失信心。所以,向未来走10步,我对下一滴水的滴落时刻就一无所知

非线性动力学数据分析

时间序列分析读书报告与数据分析 刘愉 200921210001 时间序列分析是利用观测数据建模,揭示系统规律,预测系统演化的方法。根据系统是否线性,时间序列分析的方法可分为线性时间序列分析和非线性时间序列分析。 一、 时间序列分析涉及的基本概念 1、 测量 对于一个动力系统,我们可以用方程表示其对应的模型,如有限差分方程、微分方程等。如果用t X 或)(t X 表示所关心系统变量的列向量,则系统的变化规律可表示成 )(1t t X f X =+或)(X F dt dX = 其中X 可以是单变量,也可以是向量,F 是函数向量。通过这类方程,我们可以研究系统的演化,如固定点、周期、混沌等。 在实际研究中,很多时候并不确定研究对象数据何种模型,我们得到的是某类模型(用t X 或)(t X 表示)的若干观测值(用t D 或)(t D 表示),构成观测的某个时间序列,我们要做的是根据一系列观测的数据,探索系统的演化规律,预测未来时间的数据或系统状态。 2、 噪声 测量值和系统真实值之间不可避免的存在一些误差,称为测量误差。其来源主要有三个方面:系统偏差(测量过程中的偏差,如指标定义是否准确反映了关心的变量)、测量误差(测量过程中数据的随机波动)和动态噪音(外界的干扰等)。 高斯白噪声是一类非常常见且经典的噪声。所谓白噪声是指任意时刻的噪声水平完全独立于其他时刻噪声。高斯白噪声即分布服从高斯分布的白噪声。这类噪声实际体现了观测数据在理论值(或真实值)周围的随机游走,它可以被如下概率分布刻画: dx M x dx x p 2222)(exp 21 )(σπσ--= (1) 其中M 和σ均为常数,分别代表均值和标准差。 3、 均值和标准差 最简单常用的描述时间序列的方法是用均值和标准差表示序列的整体水平和波动情况。 (1)均值 如果M 是系统真实的平均水平,我们用观测的时间序列估计M 的真实水平方法是:认为N 个采样值的水平是系统水平的真实反映,那么最能代表这些观测值(离所有观测值最近)的est M 即可作为M 的估计。于是定义t D 与est M 的偏离为2 )(est t M D -,所以,使下面E 最小的M 的估计值即为所求: 21)(∑=-=N t est t M D E (2)

《从非线性动力学到复杂系统》

《从非线性动力学到复杂系统》 段法兵 系统理论博士生课程

第一讲动态系统的发展 系统是一些相互关联的客体组成的集合,动态(动力dynamical)系统是系统状态变量,比如温度、位移、价格、信号幅值等,随着时间变化的。它的描述可以用微分方程或者离散方程。 微分方程历史悠久,可追溯到牛顿、伽利略、欧拉、雅克比等人,用以描述行星的运动轨迹。研究中发现即使满足牛顿引力定律的三体运动也非常复杂,其微分方程是非线性的,非线性是指不满足叠加定律的方程,解无法利用已知函数进行描述,如果能够描述的我们称为显式解。因此,庞加莱在1880年-1910年期间,试图利用解的拓扑几何性质来解释动态系统的运动规律,发现即使确定性系统,其运动规律也会出现随机性态,非常复杂(确定性系统是指其外力是确定的不随机,只要知道初始条件和演化方程,其运动是可预先确定的)。 非线性系统运动的复杂性:李雅普诺夫研究了系统平衡点?的稳定性?问题,随后本迪尔松等发现系统的解包含(1)平衡态(静止不动);(2)周期运动(比如行星)(3)拟周期,就是几个频率不可公约周期之和。 接着1975年Li和Yorke提出了混沌的概念,即系统的解是非周期的一种类似随机运动的现象,这其中就包含了洛伦兹提出的“蝴蝶效应”,根源在于这类非线性动力系统对于初始条件的极其敏感性,初始条件的微小变化导致了系统状态的巨大改变,从此有关非线性科学的发展异常迅速,形成了现代动力学理论,其最重要的贡献是揭示了一个简单的模型可能蕴含了无比复杂的动力学性态。 例子:Van der Pol(范德波尔)方程 1920年Van der Pol利用电子震荡管研究心脏的跳动问题,比如人工心脏起

简单非线性电阻电路的分析

第五章 简单非线性电阻电路的分析 5-1 含一个非线性元件的电阻电路的分析 一、含一个非线性元件的电阻电路都可用电源等效定理来等效 N 为含源线性网络。 二、非线性电路的一般分析方法 1、图解法 2、代数法 3、分段分析法 4、假定状态分析法 1、图解法 设非线性电阻的V AR 为 在如上图所示u 和i 的参考方向如下,线形部分的V AR 为 将 代入上式得 通常,用图解法求解u 和i 如图5-2 两曲线的交点Q 是所求解答。直线称为负载线 在求出端口电压 u Q 和 i Q 后。就 可用置换定理求出线性单口网络内部的电 ) (u f i =i R u u oc 0-=)(u f i =oc oc u u u f R u f R u u =+-=)()(00

压电流。 例5-1 电路如图5-3(a)所示,二极管特性曲线如图(d)所示,输入电压随时间变化。 (1)试求所示电路输出电压u0对输入电压u i的曲线,即u0-u i转移特性; (2)若输入电压的波形如图(e)所示,试求输出电压u0的波形。 解戴维南等效电路 由电路可知 2 i oc u u= i u u30 0 + =

若 u i 变化时(交流),戴维南等效电压源也是时变的。但Ro 是定值,所以 线性网络的负载线具有不变的斜率 -1/Ro ,在 u-i 平面上作平行移动,每一时 刻负载线在电压轴的截距总是等于等效电压源在该时刻的瞬时值,负载线与二极管特性曲线的交点也在移动,即二极管的电压、电流都随时间而变。 求u 0-u i 转移特性曲线 由图(a )可得 当 时,0u 由 确定。 当 时,0i =, 可得转移特性曲线如图5-4所示 2、代数法 若i=f(u)中的f(u)可用初等函数表示,那么可利用节点法或回路法求解。 例5-2 如图5-5所示电路中,已知非线性电阻的V AR 为 试求电流i 。 030u u i =+0>i u i u u o 30+=0

第一章 非线性动力学分析方法

第一章非线性动力学分析方法(6学时) 一、教学目标 1、理解动力系统、相空间、稳定性得概念; 2、掌握线性稳定性得分析方法; ?3、掌握奇点得分类及判别条件; ?4、理解结构稳定性及分支现象; 5、能分析简单动力系统得奇点类型及分支现象. 二、教学重点 1、线性稳定性得分析方法; ?2、奇点得判别。 三、教学难点 ?线性稳定性得分析方法 四、教学方法 讲授并适当运用课件辅助教学 五、教学建议 ?学习本章内容之前,学生要复习常微分方程得内容。 六、教学过程 本章只介绍一些非常初步得动力学分析方法,但这些方法在应用上就是十分有效得。 1、1相空间与稳定性 ?一、动力系统 在物理学中,首先根据我们面对要解决得问题划定系统,即系统由哪些要素组成。再根据研究对象与研究目得,按一定原则从众多得要素中选出最本质要素作为状态变量。然后再根据一些原理或定律建立控制这些状态变量得微分方程,这些微分方程构成得方程组通常称为动力系统。研究这些微分方程得解及其稳定性以及其她性质得学问称为动力学. 假定一个系统由n个状态变量,,…来描述。有时,每个状态变量不但就是时间t得函数而且也就是空间位置得函数。如果状态变量与时空变量都有关,那么控制它们变化得方

程组称为偏微分方程组.这里假定状态变量只与时间t有关,即X =X i(t),则控制它们 i 得方程组为常微分方程组。 ?????(1。1.1) … 其中代表某一控制参数.对于较复杂得问题来说,(i=l,2,…n)一般就是得非线性函数,这时方程(1.1.1)就称为非线性动力系统。由于不明显地依赖时间t,故称方程组(1。1.1)为自治动力系统。若明显地依赖时间t,则称方程组(1、1、1)为非自治动力系统.非自治动力系统可化为自治动力系统. 对于非自治动力系统,总可以化成自治动力系统。 例如: 令,,上式化为 上式则就是一个三维自治动力系统。 又如: 令,则化为 它就就是三微自治动力系统、 对于常微分方程来说,只要给定初始条件方程就能求解。对于偏微分方程,不但要给定初始条件而且还要给定边界条件方程才能求解。 能严格求出解析解得非线性微分方程组就是极少得,大多数只能求数值解或近似解析解。 二、相空间 ,X2,…Xn)描述得系统,可以用这n个状态变量为坐标轴支由n个状态变量=(X 1 起一个n维空间,这个n维空间就称为系统得相空间。在t时刻,每个状态变量都有一个确定得值,这些值决定了相空间得一个点,这个点称为系统状态得代表点(相点),即它代表了系统t时刻得状态。随着时间得流逝,代表点在相空间划出一条曲线,这样曲线称为相轨道或轨线.它代表了系统状态得演化过程。 三、稳定性 把方程组(1。1.1)简写如下

单摆非线性动力学

单摆的非线性动力学分析 亚兵 (交通大学车辆工程专业,,730070) 摘要:研究单摆的运动,从是否有无阻尼和驱动力方面来分析它们对单摆运动的影响。对于小角度单摆的运动,从单摆的动力学方程入手,借助雅普诺夫一次近似理论,推导出单摆的运动稳定性情况。再借助绘图工具matlab,对小角度和大角度单摆的运动进行仿真,通过改变参数,如阻尼大小、驱动力大小等绘出单摆运动的不同相图,对相图进行分析比较,从验证单摆运动的稳定性情况。关键词:单摆;振动;阻尼;驱动力 Abstract:The vibration of simple pendulum is studied by analyzing whether or not damp and drive force its influence of the simple pendulum. For small angle pendulum motion, pendulum dynamic equation from the start, with an approximate Lyapunov theory of stability of motion is derived pendulum situation. Drawing tools with help from matlab, small angle and wide-angle pendulum motion simulation, by changing the parameters, such as damping size, drive size draw simple pendulum of different phase diagram, analysis and comparison of the phase diagram, from the verification the stability of the situation pendulum movement. Key words: simple pendulum; vibration; damp; drive force 1 引言 单摆是一种理想的物理模型[1],单摆作简谐振动(摆角小于5°)时其运动微分方程为线性方程,可以求出其解析解,而当单摆做大幅度摆角运动时,其运动微分方程为非线性方程,我们很难用解析的方法讨论其运动,这个时候可以用MATLAB软件对单摆的运动进行数值求解,并可以模拟不同情况下单摆的运动。 θ=时, 随着摆角的减小,摆球的运动速率将越来越大,而加速度将单调下降,至0 加速度取极小值。本文从动力学的角度详细考察了这一过程中摆球的非线性运,得出了在运动过程中.,t θθθ --的关系。

24第20章_非线性动力分析_李永双概论

第二十章非线性动力分析 本书前面已经介绍了使用SAP2000进行线性动力分析的基本内容,线性动力分析主要任务是处理结构在多遇地震及一般动力荷载作用下的效应问题,在这阶段结构并没有进入到塑性发展阶段,因此结构的响应控制在线弹性的范围。 根据我国规范提出的结构抗震设计中“小震不坏、中震可修,大震不倒”三个设防水准,以及弹性阶段承载力设计和弹塑性阶段变形验算的两阶段设计理论,进入到大震状态(罕遇地震)是允许结构构件出现塑性发展的,并且需要程序能够进行一定深度的弹塑性分析并给出相关的效应结果。此外,目前很多实际工程中已经开始使用隔振器、阻尼器等复杂保护装置,这些装置一般需要使用非线性连接单元去模拟,而线性时程分析不能够考虑非线性连接单元的非线性属性。综上所述,特定工程需要进行相关条件下结构的非线性动力分析,也要求程序能够完成这一分析。 在SAP2000中可以进行非线性时程分析,在这一分析中可以考虑结构构件的塑性发展(塑性铰),可以考虑复杂的隔振器、阻尼器等非线性连接单元,也可以完成冲击、爆炸等复杂的动力荷载作用下结构效应分析,本章将结合这些非线性时程分析的具体问题阐述其定义方式及相关需要注意的问题。另外,需要注意的是,非线性时程分析本质上仍然是一种动力时程分析,不同之处在于它可以综合考虑结构中的非线性属性,因此部分参数选择和设置方式与线性时程分析是相同的,对于这类问题由于在线性时程分析中已经进行阐述,因此本章不会重复描述,本章的重点在于使用SAP2000进行非线性时程分析时所能够考虑的非线性属性及其意义。 20.1非线性时程分析工况的定义及相关概念 本章将分别介绍非线性时程分析的相关概念、快速非线性模态积分方法和几种常见的非线性分析类型。下面从非线性时程分析工况的定义出发,阐述非线性时程分析所涉及的几个基本概念。 20.1.1时程函数的定义 与线性时程分析相同,非线性时程分析首先需要定义时程函数曲线,定义方式与线性时程分析是相同的。如果需要进行罕遇地震作用下结构的非线性分析,需要选择地震波曲线,可以使用程序联机带有的常用地震波形式以及我国规范常用的几种场地状态下地震波曲线,可以通过峰值控制来得到罕遇地震的地震时程曲线。 除了罕遇地震作用以外,作用于结构更复杂的动力荷载一般需要提供该作用的数据形式,或工程师根据荷载特征构建荷载作用的数据形式,比如一定的冲击荷载作用或爆炸荷载作用。对于这类荷载数据形式的形成和使用方式与线性时程分析中所描述的时程曲线形成的方式相同,对于几种典型动力作用的时程曲线我们在本章后面相关专题将会再次涉及到。 20.1.2时程工况的定义 与线性时程分析相同,完成时程函数曲线定义之后,需要定义非线性时程分析工况。当选择添加新工况并在分析工况类型下拉菜单中选择Time History,可以弹出时程分析工况定义对话框。非线性分析工况定义对话框与线性时程分析对话框是相同的,见图20-1。

非线性力学和混沌简介

非线性力学和混沌简介 非线性科学是一门研究非线性现象共性的基础学科。它是自本世纪六十年代以来,在各门以非线性为特征的分支学科的基础上逐步发展起来的综合性学科,被誉为本世纪自然科学的“第三次革命”。非线性科学几乎涉及了自然科学和社会科学的各个领域,并正在改变人们对现实世界的传统看法。科学界认为:非线性科学的研究不仅具有重大的科学意义,而且对国计民生的决策和人类生存环境的利用也具有实际意义。由非线性科学所引起的对确定论和随机论、有序与无序、偶然性与必然性等范畴和概念的重新认识,形成了一种新的自然观,将深刻地影响人类的思维方法,并涉及现代科学的逻辑体系的根本性问题。 一线性与非线性的意义 线性”与“非线性”是两个数学名词。所谓“线性”是指两个量之间所存在的正比关系。若在直角坐标系上画出来,则是一条直线。由线性函数关系描述的系统叫线性系统。在线性系统中,部分之和等于整体。描述线性系统的方程遵从叠加原理,即方程的不同解加起来仍然是原方程的解。这是线性系统最本质的特征之一。“非线性”是指两个量之间的关系不是“直线”关系,在直角坐标系中呈一条曲。 最简单的非线性函数是一元二次方程即抛物线方程。简单地说,一切不是一次的函数关系,如一切高于一次方的多项式函数关系,都是非

线性的。由非线性函数关系描述的系统称为非线性系统。 线性与非线性的区别 定性地说,线性关系只有一种,而非线性关系则千变万化,不胜枚举。线性是非线性的特例,它是简单的比例关系,各部分的贡献是相互独立的;而非线性是对这种简单关系的偏离,各部分之间彼此影响,发生偶合作用,这是产生非线性问题的复杂性和多样性的根本原因。正因为如此,非线性系统中各种因素的独立性就丧失了:整体不等于部分之和,叠加原理失效,非线性方程的两个解之和不再是原方程的解。因此,对于非线性问题只能具体问题具体分析。 线性与非线性现象的区别一般还有以下特征: (1)在运动形式上,线性现象一般表现为时空中的平滑运动,并可 用性能良好的函数关系表示,而非线性现象则表现为从规则运动向不规则运动的转化和跃变; (2)线性系统对外界影响的响应平缓、光滑,而非线性系统中参数的极微小变动,在一些关节点上,可以引起系统运动形式的定性改变。在自然界和人类社会中大量存在的相互作用都是非线性的,线性作用只不过是非线性作用在一定条件下的近似。 非线性问题研究的历史概况

三维柔性多体梁系统非线性动力响应分析

振 动 与 冲 击 第25卷第1期 J OURNAL OF V IBRATI ON AND SHOCK Vo.l 25No .12006  三维柔性多体梁系统非线性动力响应分析 收稿日期:2005-01-09 第一作者吴国荣男,博士,副教授,1966年10月生 吴国荣1  钟伟芳2  梁以德 2 (1.福建工程学院土木系,福州 350007; 2.华中科技大学力学系,武汉 430074) 摘 要 研究了三维柔性多体梁系统的非线性动力响应问题。将空间柔性梁的变形分解为轴向变形以及在x -y 平面的弯曲变形和在x -z 平面的弯曲变形,引用各自的精确振动模态描述变形场,利用拉格朗日乘子法建立起柔性多体梁系统约束非线性动力学方程。结合N e wma rk 直接积分法和N e w t on -R aphson 迭代法,导出了求解该非线性代数-微分 方程组的数值方法。仿真算例证明了该方法的正确性和有效性。 关键词:多体系统,三维柔性梁,非线性,动力响应 中图分类号:O 32,TB1 文献标识码:A 0 引 言 随着现代机械系统向高速、轻质和高精度发展,使 得在机械系统设计和动力学分析时,应用传统的多刚体系统动力学分析方法不能得到正确的结果,必须考虑构件在大范围刚性运动时柔性变形的影响。这类计及部件变形的多体系统称为柔性多体系统。近年来,考虑刚柔耦合的柔性多体系统的研究已成为多体系统中最主要的研究方向,各国学者对此作了大量的研究工作,提出了不少新的方法处理柔性多体系统动力学 问题,如有限元方法[1-5],有限段方法[6-8] ,子结构方法[9-10],模态综合法[11-12]。 研究了三维多体梁系统非线性动力响应问题。将空间柔性梁的变形分解为轴向变形和弯曲变形,并用各自的解析模态来近似描述。使用带乘子的拉格朗日方程导出了柔性多体系统动力学方程。基于N e wm ark 直接积分法结合N e w ton -Raphson 迭代法[13] ,给出了求 解该约束非线性方程的数值方法。数值模拟结果表明了本文方法的有效性 。 图1空间柔性梁的弹性位移示意图 1 运动学 建立如图1所示的坐标体系来描述作空间运动的柔性梁的构形,其中OXYZ 为惯性坐标系,O i X i Y i Z i 为随体坐标系。柔性体的变形在随体坐标系中描述。随体坐标系的方位由三个相互独立的欧 拉角 i , i 及ψi 确定。考虑其上任一点P i ,它的空间位置可表示为: r i P =R i +A i u i (1)上式中:R i 为随体坐标系的坐标原点的位置矢量,u i 为点P i 变形后的局部坐标,矩阵A i 为坐标变换矩阵,其定义为: A i =C i C ψi -C i S i S ψi -S ψi C i -C i S i C ψi S i S i S i C ψi + C i C i S ψi -S ψi S i +C i C i C ψi -S i C i S i S ψi S i C ψi C i (2) 其中S =si n ()及C =cos ()。 可把矢量u i 分解为: u i =u i o +u i f (3) 其中u i o 是在变形前点P i 的局部坐标,而且u i f 为随体坐标下该点的形变位移矢量。 把式(3)代入式(1),得: r i p =R i +A i (u i o +u i f )(4) 方程(4)的两边对时间求导,并注意到·u i 0=0,得:·r i p =·R i +A i (ωi ×(u i o +u i f ))+A u i i f (5)这里ωi 为随体坐标系下的角速度,它可表示为: ωi =G i θ·i (6) 其中θi =[ i i ψi ],G i 为如下定义的矩阵: G i =S i S ψi C ψi 0 -S i C ψi -S ψi 0C i 1 2 空间梁的变形描述 假设柔性构件为空间Ber noulli -Euler 梁,其弹性变形分 DOI 牶牨牥牣牨牫牬牰牭牤j 牣cn ki 牣jvs 牣牪牥牥牰牣牥牨牣牥牥牱

结构非线性动力分析方法综述_周文峰

·自然科学研究· 结构非线性动力分析方法综述 周文峰 郭 剑 (攀枝花学院土木工程学院,四川攀枝花 617000) 摘 要 时程分析法是一种计算机模拟分析方法,其优势在于能模拟出结构进入非弹性阶段的受力性能。该 方法主要包括结构分析模型、单元模型和恢复力模型三个重要方面。本文从这三个方面简单介绍了结构非线 性动力反应分析方法。 关键词 非线性;动力分析;模型 结构抗震设计方法经历了静力阶段、反应谱阶段和动力阶段。从本质上说,前二者所采用的方法均为静力法,且只能进行弹性分析。动力阶段的形成建立在计算机的普及和数值分析方法的出现基础之上,其分析方法称为时程分析法。时程分析法本质上是一种计算机模拟分析方法,能够计算出结构地震反应的全过程,该方法的突出优势在于能模拟出结构进入非弹性阶段的受力性能。 时程分析法的出现促进了结构非线性地震反应分析的发展。它主要包括结构分析模型、单元模型和恢复力模型三个重要方面,下面从这三个方面进行简单介绍。 1 结构分析模型 结构的模型化是非线性动力反应分析的第一步,结构模型的模拟应着重于其动力特性的模拟。因此体系恢复力、质量、阻尼模型的准确性是模拟精度的前提。目前的结构分析模型可分为以下几类: 1.1 层间模型 考虑到框架结构质量的分布规律,很容易形成以楼层为单元的多质点体系的思路,故将这种模型称之为层间模型。在研究框架结构动力反应时,层间模型中采用得最多的是层间剪切型模型。该模型假定框架结构层间变形以剪切变形为主,忽略其它形式变形的影响,故而比较适用于高跨比不大、层数不多的框架。为了进一步拓宽此模型的适用范围,在此模型基础上又发展了层间剪弯型模型,使之能适用于层数较多和高跨比较大的框架。 但是层间模型在实际使用中却存在比较大的困难,这主要反映在如何具体确定层间的剪切刚度及弯曲刚度的问题上,而且这二者之间又是耦合在一起的。这一问题层间模型自身是无法解决的。目前,层间模型只是对于常见的层数不多且平面布置十分简单、规则、对称并且能简化为平面结构的框架有一定的实用性,也就是说对于这类框架通常能根据经验进行适当的假设后进行简单推导得到层间单元刚度。 1.2 杆系模型 杆系模型是将整体结构离散为梁、柱单元进行分析。杆系分析模型的出现不仅解决了层间模型所面临的层间刚度无法确定的困难,而且它还解决了层间模型所固有的另外两个缺陷。其一,如果说层间模型从宏观(层单元)角度展示了结构总体动力反应规律,那么由于框架各杆进入非弹性阶段的先后次序不同所造成的整个框架动力反应规律的不同,则是层间模型所不能解释、反映的。其二,无论从抗震研究还是设计角度来看,框架结构的梁、柱构件在地震作用下的反应规律到底如何也是人们所关心的,因为结构的设计最终要落实到构件的设计。如柱端弯矩增大系数应如何取值等,这些问题采用层间模型是无法回答的,从这个角度看也必须将框架结构细化到至少是构件层次才有可能解决这些问题。 杆系分析模型分为两大类,平面杆系分析模型与空间杆系分析模型。目前,平面杆系分析模型的研究相对较为成熟,国内外已开始将注意力转向空间杆系分析模型的研究上。 2 单元模型 对于杆系分析模型,目前用于模拟单元滞回性能的模型已有很多,这些单元分析模型可采取分类的方式加以比较考察。这些模型大致可分为两大类若干小类。 2.1 集中塑性铰模型 单分量模型是集中塑性铰模型中最简单的一类,该模型将杆单元的非弹性性能用非线性弹簧反映,而不对非弹性变· 109·第23卷第4期 攀枝花学院学报 2006年8月V o l .23.N o .4 J o u r n a l o f P a n z h i h u a U n i v e r s i t y A u g .2006

电路的分析方法电子教案

第2章 电路的分析方法 本章要求: 1. 掌握支路电流法、叠加原理和戴维宁定理等电路的基本分析方法。 2. 理解实际电源的两种模型及其等效变换。 3. 了解非线性电阻元件的伏安特性及静态电阻、动态电阻的概念,以及简单非线性电阻电路的图解分析法。 重点: 1. 支路电流法; 2. 叠加原理; 3.戴维宁定理。 难点: 1. 电流源模型; 2. 结点电压公式; 3. 戴维宁定理。 2.1 电阻串并联联接的等效变换 1.电阻的串联 特点: 1)各电阻一个接一个地顺序相联; 2)各电阻中通过同一电流; 3)等效电阻等于各电阻之和; 4)串联电阻上电压的分配与电阻成正比。 两电阻串联时的分压公式: 2.电阻的并联 特点: 1)各电阻联接在两个公共的结点之间; 2)各电阻两端的电压相同; 3)等效电阻的倒数等于各电阻倒数之和; 4)并联电阻上电流的分配与电阻成反比。 U R R R U 2111+=U R R R U 2 122+=

两电阻并联时的分流公式: 2.3 电源的两种模型及其等效变换 1.电压源 电压源是由电动势 E 和内阻 R 0 串联的电源的电路模型。若 R 0 = 0,称为理想电压源。 特点: (1) 内阻R 0 = 0; (2) 输出电压是一定值,恒等于电动势(对直流电压,有 U ≡ E ),与恒压源并联的电路电压恒定; (3) 恒压源中的电流由外电路决定。 2.电流源 电流源是由电流 I S 和内阻 R 0 并联的电源的电路模型。若 R 0 = ∞,称为理想电流源。 特点: (1) 内阻R 0 = ∞ ; (2) 输出电流是一定值,恒等于电流 I S ,与恒流源串联的电路电流恒定; (3) 恒流源两端的电压 U 由外电路决定。 3.电压源与电流源的等效变换 等效变换条件: E = I S R 0 0 R E I = S 注意: ① 电压源和电流源的等效关系只对外电路而言,对电源内部则是不等效的。 ② 等效变换时,两电源的参考方向要一一对应。 ③ 理想电压源与理想电流源之间无等效关系。 ④ 任何一个电动势 E 和某个电阻 R 串联的电路,都可化为一个电流为 I S 和这个电阻并联的电路。 4.电源等效变换法 (1) 分析电路结构,搞清联接关系; (2) 根据需要进行电源等效变换; (3) 元件合并化简:电压源串联合并,电流源并联合并,电阻串并联合并; I R R R I 2121+=I R R R I 2 112+=

钢筋混凝土结构非线性动力反应分析

1 问题描述 (1) 2 模型建立 (2) 2.1 轴网布置 (2) 2.2 材料定义 (2) 2.3 截面定义 (3) 2.4 楼板定义 (8) 2.5 构件布置............................................................... 错误!未定义书签。 2.6 构件剖分 (10) 2.7 节点束缚指定 (11) 2.8 先单元插入点指定 (11) 2.9 刚域指定 (12) 2.10 整体模型 (12) 2.11 荷载模式定义 (13) 2.12 质量源定义 (14) 2.13 荷载定义 (14) 2.14 荷载工况定义 (15) 2.14.1 模态分析工况 (15) 2.14.1 反应谱工况 (16) 2.14.3 COMBINE工况(计算柱轴力) (17) 2.14.3 时程反应分析工况 (18) 2.15 定义组 (21) 3 振型分解反应谱法结果 (22) 3.1 楼层剪力及水平地震作用 (22) 3.2 内力图 (23) 3.2.1 弯矩图 (23) 3.2.2 剪力图 (24) 3.2.3 轴力图 (25) 4 塑性铰的计算与定义 (26) 4.1 梁铰柱铰的计算 (26) 4.1.1 梁铰的计算 (26) 4.1.2 柱铰的计算 (29) 4.2 铰属性设置 (31) 4.2.1 梁铰属性设置 (31) 4.2.2 柱铰的设置 (32)

4.3 梁柱铰的指定 (35) 5 时程反应分析结果 (36) 5.1 层位移 (36) 5.2 层间位移 (38) 5.3 层间位移角 (38) 5.4 层剪力 (39) 5.5 出铰情况及破坏机制 (40)

非线性动力反应分析

目录 第一章采用底部剪力法分别计算一榀平面和三维整体框架的地震内力 (2) 1.1 底部剪力法计算B轴线一榀平面框架的地震内力 (2) 1.11 重力荷载代表值的计算 (2) 1.12 框架的侧移刚度和基本周期 (4) 1.13 各楼层的水平地震作用力和剪力 (7) 1.2 底部剪力法计算三维整体结构时B轴框架的地震内力 (13) 1.21 重力荷载代表值的计算 (13) 1.22框架的侧移刚度和基本周期 (15) 1.23各楼层的水平地震作用力和剪力 (17) 1.3三维整体结构与一榀平面框架在地震作用时内力对比图 (21) 第二章采用SAP2000程序计算结构朝X方向的前五阶振型和周期 (24) 2.1沿X方向的前五阶振型及其周期 (28) 第三章采用振型分解法计算B轴框架中产生的地震内力 (29) 3.1计算水平地震作用力 (29) 3.2地震作用下各振型相应的柱子的剪力 (31) 3.3地震作用下各振型相应在B轴框架的内力 (32) 3.4采用平方和开平方求B轴框架的最大的反应 (34) 3.5振型分解法和底部剪力法求解三维整体结构时B轴框架的内力对比图 41

1.1采用底部剪力法计算B 轴线一榀平面框架的地震作用和地震内力 1.11重力荷载代表值的计算 i +E K E ik G G Q ?=∑, K G 其中为恒载标准值 ,i E ?为组合值系数,ik Q 为活载标准值 恒载标准值: 本设计无填充墙,,混凝土的容重为3 25/KN m ,为了计算方便,构件的 重叠区的质量将被重复计入, 板:0.12×25=32/K N m 次梁:0.3×0.45×25=3.375/K N m 主梁:0.3×0.5×25=3.75/KN m 柱子:0.5×0.5×25=6.25/KN m 剪力墙 0.2×25=52 /KN m 各楼层的恒载: FL1: ()()3617 3.37563 3.7517 3.75646.25 3.753245 3.753262807.375K N ??+??+?+??+ ?+÷?+?+÷??= FL2: ()()3617 3.37563 3.7517 3.75646.253324533262775.5K N ??+??+?+??+ ?+÷?+?+÷??= FL2=FL3=FL4=FL5=FL6 FL7: 3617 3.37563 3.7517 3.75646.2532453262648K N ??+??+?+??+?÷?+?÷??= 各楼面活载:

非线性系统的一些动力学与控制问题

釜七届全国非线性动力学学术会议和第九届全国非线性振动学术会议论文集南京,200410.28-29复杂非线性系统的一些动力学与控制问题。 陆启韶王士敏 f北京航空航天大学理学院北京10083) E-mail:qishaolu(岔hotraail.eom 摘要本文根据非线性动力学的研究现状和发展趋势,对复杂非线性系统动力学与控制的理论和应用研究中的一些重要问题进行探讨和展望 关键词非线性,复杂系统,动力学,控制 前言 非线性动力学研究非线性系统丰富的运动模式和演化过程,是非线性科学技术的重要理论基础。非线性动力学研究的最终目的在于深刻揭示非线性世界的复杂性和多样性。非线性系统运动的复杂性来源于多个方面,例如几何关系、本构关系、约束条件、拓扑结构、激励因素、耦合方式、时空尺度、演化机理等,它们都会带来复杂的运动模式。30多年来,尽管非线性动力学对单自由度简单振动系统和低维映射系统的研究已经取得一系列重要成果,发现了大量新的非线性现象。提出并发展了基本的理论方法,但是面对在理论和应用研究中遇到的高维复杂系统问题往往束手无策,仍然缺乏有效的分析策略和手段。因此,复杂非线性系统研究已成为当务之急。 本文根据当前非线性动力学的研究现状和发展趋势,针对复杂非线性系统动力学与控制的理论和应用研究中的一些重要问题进行探讨和展望,希望引起同行关注,共同开创该方面研究的新局面。1.多自由度非线性系统组合振动、全局分析和同步实际非线性振动系统通常是多自由度的,且存在多种外界激励,因此组合振动和模态相互作用是普遍的重要现象。对单自由度系统来说,组合共振只能在多种激励并存的情形下出现。但是对多自由度系统,由于可以存在内共振和自参数共振机理,因 ’国家自然科学基金(10172011)资助项目此在单个激励作用下也可能发生组合共振。内共振(或自参数共振)发生在其线性化系统的各模态的固有频率可以通约或接近通约的情况,其类型依赖于非线性项形式和相应的分岔类型。在没有内共振时,系统的共振响应只包含由外部激励直接激发的主共振或亚,超谐共振模态。但是内共振会引起与非线性项有关的间接激发模态,并导致多模态相互作用,产生诸如饱和、跳跃、锁相、周期调制、混沌调制等复杂现象,造成弹性结构中由高频激励引起的低频大幅共振事故。现在对多自由度系统的组合振动和模态相互作用动力学研究已经取得一些重要成果,并且扩展到梁、板、壳、弦线、悬索、传送带、流一固耦合结构等系统,涉及不同的本构关系(包括粘弹性材料、复合材料、智能材料等)、约束条件和控制方式,成为十分活跃的研究方向。但是,目前这方面的研究主要局限于具体问题,对于组合振动的一般规律和分析方法仍有待于深入探讨。 高维非线性振动系统的全局动力学分析是十分重要且难度很大的问题。目前仍然主要依靠数值模拟手段.成功地用于全局分析的理论方法不多,主要是高维Melnikov方法和Shilnikov方法。近年来,人们发现了大重新的非线性动力学现象,除了混沌激变、瞬态混沌、奇怪混沌不变集之外,还有超混沌、Wada吸引域、筛形吸引域、混沌鞍等,需要从机理上予以明确阐述。因此,当务之急是将动力系统理论、强非线性系统