12静力非线性
非线性静力分析需求谱建立方法的两点改进
+& + ( s n ) , gx =一 i () t () 6
法, 采用弹塑性反应谱代替等效高阻尼弹性反应谱 . 目前常用的有等延性强度需求谱和等强度延性需求 谱两种不同的弹塑性反应谱 . 以往估计 结构 的 目 标 位移 , 一般采用等延性强度需求谱 的统计结果反推 出结构的延性系数 , 文献 [] 4采用等强度延性谱来估
中. 非线性静 力分析 的研 究 目前有若干分析 方法, 这些方 法在 结构能 力曲线建 立 的方法上是 相 同的 , 即首先建 立力 一位移 曲 线, 然后转化为承载力谱 , 但是在 需求 曲线 的建 立上各有不 同. 通过介绍几种典 型的非 线性静力 分析 需求谱 建立方法 , 出其 指 和原有的需求谱建立方法的 不同和 改进之 处 .
用等延性强度需求谱 的统计结果反推出结构 的延性 系数 . 文献[ ] 4 认为直接采用等强度延性谱来估计地 震位移需求 , 更加方便 、 简单 , 且计算精度也更高 .
1 考虑高振 型特性的影响
1 1 多 自由度弹性 体 系高振 型分 析 . 多 自由度体 系 在地震 力作 用 下 的振 动反 应可 由 公 式 () 1表示
系数 . 得
() 2 采用等强度延性谱来估计地震位移需求 . 等延性强度需求谱和等强度延性需求谱是两种不同
的弹塑性反应谱 . 以往估计结构的 目 标位移 , 一般采
收 稿 日期 :O 8 3 l 2O —0 一l
作者简介 : 朱星平 ( 6 一 , , 1 5 ) 江西省修水县人 , 9 男 江西修水建筑设计院一级注册结构师 , 工程师 , 研究方 向为建筑结构设计
() t :∑
n I
() t
静力弹塑性分析(Pushover分析)两种方法剖析
静力弹塑性分析(Pushover 分析)■ 简介Pushover 分析是考虑构件的材料非线性特点,分析构件进入弹塑性状态直至到达极限状态时结构响应的方法。
Pushover 分析是最近在地震研究及耐震设计中经常采用的基于性能的耐震设计(Performance-Based Seismic Design, PBSD)方法中最具代表性的分析方法。
所谓基于性能的耐震设计就是由用户及设计人员设定结构的目标性能(target performance),并使结构设计能满足该目标性能的方法。
Pushover 分析前要经过一般设计方法先进行耐震设计使结构满足小震不坏、中震可修的规X 要求,然后再通过pushover 分析评价结构在大震作用下是否能满足预先设定的目标性能。
计算等效地震静力荷载一般采用如图2.24所示的方法。
该方法是通过反应修正系数(R)将设计荷载降低并使结构能承受该荷载的方法。
在这里使用反应修正系数的原因是为了考虑结构进入弹塑性阶段时吸收地震能量的能力,即考虑结构具有的延性使结构超过弹性极限后还可以承受较大的塑性变形,所以设计时的地震作用就可以比对应的弹性结构折减很多,设计将会更经济。
目前我国的抗震规X 中的反应谱分析方法中的小震影响系数曲线就是反应了这种设计思想。
这样的设计方法可以说是基于荷载的设计(force-based design)方法。
一般来说结构刚度越大采用的修正系数R 越大,一般在1~10之间。
但是这种基于荷载与抗力的比较进行的设计无法预测结构实际的地震响应,也无法从各构件的抗力推测出整体结构的耐震能力,设计人员在设计完成后对结构的耐震性能的把握也是模糊的。
基于性能的耐震设计中可由开发商或设计人员预先设定目标性能,即在预想的地震作用下事先设定结构的破坏程度或者耗能能力,并使结构设计满足该性能目标。
结构的耗能能力与结构的变形能力相关,所以要预测到结构的变形发展情况。
所以基于性能的耐震设计经常通过评价结构的变形来实现,所以也可称为基于位移的设计(displacement-based design)。
非线性有限元法综述
非线性有限元法综述摘要:本文针对非线性有限元法进行综述,分别从UL列式及TL列式、CR列式、几何精确梁、壳理论三个方面介绍其分析思路和发展动态,旨在为相关学者提供一些思路参考。
关键词:几何非线性;UL列式;TL列式;CR列式;几何精确梁、壳理论1引言几何非线性是由于位置改变引起了结构非线性响应。
进行结构几何非线性分析,实质上就是要得到结构真实的变形与受力情况。
有限元方法是进行结构几何非线性分析的最成熟的方法,也是应用最广泛的分析方法.2非线性有限元法研究思路非线性有限元法主要指UL列式法、TL列式法、CR列式法和几何精确梁、壳理论等,它们有着基本相同的思路,即利用虚功原理建立平衡方程。
方程中充分考虑了非线性因素对结构应变和应力的影响,也就是将线性应变和非线性应变都代入到表达式中,然后确定单元的本构关系并选取合适的形函数,导出单元对应的弹性刚度矩阵和几何刚度矩阵,再选取合适的增量-迭代算法进行求解,由此就完成了结构的整个几何非线性分析求解过程。
非线性有限元法将结构的变形过程划分为三个主要阶段:C0状态、C1状态和C2状态,如图1所示。
图1 单元的变形C0状态是单元的初始状态,C1状态是单元受力变形后上一次处于平衡的状态;C2状态是单元的当前状态,也就是所求的状态。
2.1UL法和TL法研究思路UL法和TL法为几何非线性问题提供了新的分析思路。
这两种方法本质上没有很大区别,但是方程建立的参考状态有所不同。
完全拉格朗日法(TL法)是以结构变形前C0状态为参考建立平衡方程的,考虑结构从C0状态到C2状态之间的变形;而更新的拉格朗日法(UL法)以结构变形后C1状态为参考建立平衡方程的[2],考虑结构从C1状态到C2状态之间的变形。
两种拉格朗日法的主要形式如下:(1)TL列式(2)UL列式从上面两式可以看出:TL法和UL法的另一个不同是TL法的增量平衡方程中考虑了初位移矩阵的影响,而UL法则忽略了其影响,只考虑了弹性刚度矩阵和初应力矩阵的影响。
有限元非线性分析
下表简要列出了线性和非线性有限元分析之间的主要不同。关于荷载-位移关系、应力-应变关系、应力-应变度量 等主要不同将在本章详细介绍。
序号 1.
特征 荷载-位移关系
2.
应力-应变关系
3.
比例缩放
4.
线性叠加
5.
可逆性
6.
求解序列
7.
计算时间
8.
用户与软件的交互
13.3 非线性的类型
2)对数应变和真实应力 对数应变/自然应变/真实应变是度量大应变的方法,计算公式如下:
它是非线性应变的度量,因此是关于最终长度的非线性函数。与线性应变相比,对数应变(或真实应变)是可加
的。考虑一个初始长度为1m的杆经过下面3步的变形: 第1步: 从1m 变形至1.2m 第2步:从1.2m 变形至1.5m 第3步:从1.5m变形至2m 在下表中我们比较了工程应变和真实应变。可以清楚地看到,只有真实应变是可加的,因此在非线性分析中应该
13.6 非线性静力分析的一般流程
一个典型的非线性静力分析项目需要以下步骤:
网格划分:有限元模型的创建是有限元分析一个非常重要的步骤,不论进行什么样的分析。在第4-7章已经讨论过对 于某些应用的如何选择适当的单元类型。FEA小组会得到零件的几何数据,需要对这些几何进行网格划分以得到零件 网格。当装配中所有的零件划分网格后,使用适当的连接单元把它们都连接在一起如CWELD或CBUSH。一般来说, 四边形单元和六面体单元优于三角形单元、楔形单元和四面体单元。应该注意模型中的关键特征,比如圆角、孔和倒 角。如果在两个平行表面之间有紧固件或焊接,应该尽量在两个面上创建相似的网格。这将有助于焊接单元或刚性单 元垂直于表面而不破坏壳单元。然而,许多有限元分析(FEA)代码支持不依赖于节点焊接,而是基于绑定接触。这 允许用户在两个焊接零件之间创建不依赖于节点的连接单元。建议首先对复杂零件进行网格划分,然后对简单或平面 几何进行网格划分以保证良好的单元质量。需要用适当的方式来模拟夹紧、铰接和焊接以在结构中正确地传递荷载。 为单元定义适当的刚度和预荷载以得到更高的精度。如果荷载从结构上的某个面传递到另一个面上,应该在两个面间 定义接触。每个FEA代码都有自己的接触参数输入格式。一个典型的接触定义需要主从节点或单元,摩擦系数,接触 面间的间隙和接触算法。
迈达斯之——静力弹塑性分析基本原理及方法
m i d a s C i v i lm i d a s C i v i lm i d a s C i v i l图2.8.38 基于位移设计法的结构抗震性能评价m i d a s C i v i l示。
m i d a s C i v i lm i d a s C i v i lm i d a s C i v i lm i d a s C i v i lm i d a s C i v i lm i d a s C i v i lm i d a s C i v i lm i d a s C i v i lm i d a s C i v i lm i d a s C i v i lm i d a s C i v i l1n λ- : 前一步骤(n-1)的荷载因子1λ : 第1荷载步的荷载因子nstep : 总步骤数i : 等差增量步骤号当前步骤的外力向量如下。
0n n λ=⋅P P(10)(3) 第3阶段: 最终步骤的荷载增量(n nstep =) 最终荷载步骤(nstep )的外力向量如下、0nstep nstep λ=⋅P P ; 1.0nstep λ= (11)图2.8.43 自动调整荷载步长的例题(荷载因子结果)m i d a s C i v i l2. 点击步长控制选项 > 增量控制函数定义步长控制函数m i d a s C i v i lm i d a s C i v i lm i d a s C i v i lm i d a s C i v i lm i d a s C i v i lm i d a s C i v i lm i d a s C i v i lm i d a s C i v i lm i d a s C i v i lm i d a s C i v i lm i d a s C i v i lm i d a s C i v i lm i d a s C i v i lm i d a s C i v i lm i d a s C i v i lm i d a s C i v i lm i d a s C i v i lm i d a s C i v i lm i d a s C i v i lm i d a s C i v i lm i d a s C i v i lm i d a s C i v i lm i d a s C i v i lm i d a s C i v i lm i d a s C i v i lm i d a s C i v i lm i d a s C i v i lm i d a s C i v i lm i d a s C i v i lm i d a s C i v i lm i d a s C i v i lATC-40中对不同结构响应类型规定了谱折减系数的下限值(参见表2.8.7)。
第9章 非线性静力平衡方程组的解法
移近似值
{δ i+1}= {δi }+ {∆δ }i+1
(9-19)
式(9-17)即表示切线刚度迭代的策略,直至满足收敛为止。
2.迭代流程
切线刚度迭代法和割线刚度法都为变刚度的迭代法。切线刚度迭代法要求
混凝土的非线性应力应变关系具备下列形式
{dσ } = [Dt (ε )]{dε }
式中[Dt(ε)]为切线弹性矩阵,为应变的函数。 实施切线刚度迭代的具体步骤如下:
3.中点刚度增量法(又称为中点龙格-库塔法) 这种方法首先以上一增量步分析结束时的切线刚度和一半的当前荷载增
量 0.5{∆Ri}计算中点节点位移增量{∆δi-1/2}以及中点位移{δi-1/2}
[ ]{ K i−1 ∆δ i−1/
2}
=
1 2
{∆Ri
}
(9-11)
{δ i−1/ 2} = {δ i} + {∆δi−1/2}
(9-1)
[K(δ )] = ∑ ∫ [B]T [D(δ )][B]dV
(9-2)
从式中可以看出结构整体刚度矩阵[K(δ)]与节点位移{δ}相关,也就是说在
分析过程中,随着单元节点位移的变化,其弹性矩阵[D(δ)]也将发生变化,进
而结构整体刚度矩阵中的元素也在不断发生着变化,以体现材料非线性所带
来结构性能的变化。针对钢筋混凝土结构而言,其具体体现为混凝土的开裂
二、切线刚度迭代法
1.方法描述
假设式的真实解为{δ*},则非平衡力为零(由虚位移原理推导出来),式
(9-13)即变为
{Ψ(δ * )}= 0
(9-15)
若有一组近似解{δi},将上式在{δi}处进行泰勒展开并舍去高阶项,则得
静力弹塑性分析原理
第4章 静力弹塑性分析原理4.1 概要4.1.1 非线性分析的目的非线性抗震分析方法可分为非线性静力分析方法和非线性动力分析方法。
非线性动力分析方法可以认为是比较准确的方法,但是因为分析时间较长并对技术人员理论水准有较高的要求,所以在实际工程上的普及应用受到了限制。
相反静力分析方法虽然在反映结构动力特性方面有所不足,但是因为计算效率较高和操作简单、理论概念清晰等原因被广大设计人员所普遍使用。
静力弹塑性分析又被称为Pushover 分析,是基于性能的抗震设计(Performance-Based Seismic Design, PBSD)中最具代表性的分析方法。
所谓基于性能的抗震设计是以某种目标性能(target performance)为设计控制目标,而不是单纯的满足规范要求的极限承载能力的设计方法。
其步骤是先按照规范要求进行抗震分析和构件设计,然后通过Pushover 分析获得结构的极限承载能力,最后通过非线性位移结果评价结构是否满足目标性能要求。
目前规范中推荐的基底剪力法和反应谱分析方法均为弹性分析方法,其评价标准是地震作用下的抗力不小于地震作用下产生的内力,这些方法也被称为基于荷载的设计方法。
而基于性能的设计方法则是使用与结构损伤直接相关的位移来评价结构的变形能力(耗能能力),所以又被称为基于位移的设计(displa cement-based design)方法。
通过Pushover 分析可得如图4.1.1所示的荷载-位移关系曲线(能力谱),根据结构耗能情况可得到非线性需求谱。
能力谱与需求谱的交点就是结构对于地震作用的性能点(performance point)。
性能点意味着结构对于地震作用所拥有的最大的非线性承载力和最大位移,该点在控制目标性能范围内则表示该结构满足了性能要求。
通过非线性分析可以了解结构具有的的极限承载能力和安全度。
Spectral DisplacementS p e c t r a l A c c e l e r a t i o nS ddesign图4.1.1 基于位移设计法的结构抗震性能评价4.1.2 静力弹塑性分析的抗震设计原理基于性能的抗震分析方法有下列四种。
PUSHOVER分析
静力非线性(Pushover)分析静力非线性(包括 pushover)分析是一个强有力的功能,仅提供在ETABS 非线性版本中。
除了为基于抗震设计性能执行 Pushover 分析外,此功能还可用于执行常规静力非线性分析和分段式(增加)构造的分析。
执行任何非线性将花费许多时间与耐性。
在执行静力非线性分析前,请仔细阅读下列全部信息。
要特别注意其中的重要事项。
非线性静力非线性分析中可以考虑几类非线性特征。
在框架/线单元中不连续的用户定义铰的材料非线性。
铰沿着任何框架单元长度指定到任何位置数上(参见线对象的框架非线性铰指定)。
非耦合弯矩、扭矩、轴力和剪力铰是有效的。
也有根据铰位置上的交互作用轴力和弯矩所屈服的耦合 P-M2-M3 铰。
在相同的位置可存在多于一种的铰类型。
例如,可以指定一个 M3(弯矩)和一个 V2(剪力)铰到框架单元的相同端部。
所提供的默认铰属性是基于 ATC-40 和 FEMA-273 标准的。
在连接单元中材料的非线性。
有效非线性特征包括沿任何自由角度的缝隙(仅压力)、hook(仅张力)、单轴塑性,以及两种基本隔震器类型(双轴塑性和双轴磨擦/摆动)(参见线对象的连接属性指定)。
连接阻尼属性在静力非线性分析中没有效应。
所有单元中的几何非线性。
可以选择仅考虑 P-△ 效应或考虑 P-△ 效应加上大位移(请参见几何非线性效应)。
大位移效应考虑变形配置的平衡,并允许用于大平移和旋转。
但是,每个单元中的应变被假设保留为小值。
分段(顺序)施工。
在每个分析工况中,可按阶段施工顺序添加或删除构件(请参见静力非线性分段施工)。
分析工况静力非线性分析可由任何数量的工况组成。
每个静力非线性工况在结构中可有不同的荷载分布。
例如:典型静力非线性分析可由三种工况组成。
第一种为结构应用重力荷载,其次为在结构的高度上应用一个横向荷载分布,第三种将在结构高度上应用另一个横向荷载分布。
静力非线性工况可从零初始状态开始,或从前一工况末的结果开始。
sap2000线性和非线性
sap2000线性和非线性sap线性与非线性的比较从某种角度来讲,索的非线性问题计算结果不收敛,就等于劳而无获。
所以一定要将结果在保证精度下,调为收敛。
从荷载作用方式角度,非线性问题分为静力非线性和动力非线性。
静力非线性收敛问题初始条件(Initial Condition):初始条件选为零,这就意味着索初始应力为零,索很柔。
这种情况下,在此基础之上的工况分析,就很难收敛;如果给索一点初应力(施加应变或降低温度),使索具有一定的初始刚度,再以此应力为初始条件进行工况分析,计算速度就大有提高。
但是会与零初始条件计算结果有一定差别,这需要变动初应力和工况荷载进行反复试算。
荷载的施加(Load Application) :这里只讨论载控制(程序默认),在这种情况下,注意节点监测位移(monitored displacement)(迭代过程中,对此节点位移容差进行判断)可以自己选择;如果在水平荷载作用下,选择结构顶点和离地面较低点,一般前者收敛速度会慢(带有索的的结构会明显感觉到这一点),如果将监测位移的节点选择在索上,那就更慢(此时,因为有索,一般把P-delta+大变形打开;在竖向荷载下,只要不选择为索的节点,差别不大。
非线性参数(Nonlinear Parameters):(1)几何非线性参数。
这个就是按照“无”、“P-delta”和“P-delta+大变形”顺序收敛速度减慢,有索的结构体系也应该选择“P-delta+大变形”。
(2)求解控制参数。
如果默认参数情况下,不能够收敛,先不要改变收敛容差,先将最大总步、最大空步和每步最大迭代次数加大,一般情况下,都可以得到收敛的结果;否则,改变收敛容差,一般不小于5.0E-3。
动力非线性直接积分法收敛问题SAP2000动力分析方法有四种:(1)线性振型叠加法,是大多数教材里讲到的振型叠加法,它只适用于线性结构。
(2)非线性振型叠加法,针对于结构中加有非线性单元(粘弹、粘滞、摩擦、阻尼单元、隔震支座)的一种快速非线性解法(FNA),它适用于结构本身必须是线性的,非线性单元里的刚度可以为非线性(例如:考虑屈服刚度的隔震支座)的结构体系。
一种非线性隔振器的静力分析与实验研究
合 正刚 度 结构 , 为 一个 低 动 刚度 结 构 。还 有 些研 成 究者将 带 有 负刚度 的结构 单元 或者用 磁 体产 生负 刚 度 与 带 有 正 刚度 的 线 性 弹 簧 并 联 , 产 生 低 动 刚 来 度 [4 2] 论 文利 用 受 压梁 结 构 的屈 曲变形 来 指 导 -。本 设 计 出一种 具有 高静 低动 非线 性 隔振 器 。
有 高 静 刚度 一小 静 变 形 , 动 刚度 一 力 的变 化 率 很 小 。 低
关键词:振动 与波 ;高静低动 ; 刚度 ; 负 薄壁梁 ; 隔振器 中图分类号: B5 5. T 3 1 文献标 识码: A DO 编码 :03 6 8i n10 .3 52 1. .1 I 1.99 .s.0 615 .0 1 60 5 s 0
Ke r s: i r t n a d wa e; LDS; e ai esi n s t i e m ; i r t n io a o ywo d vb a i v HS o n n g t t e s; h n b a v f v b ai s lt r o
隔振 对 于大 部分 的工 程 结构 是 必须 要考 虑 的 问 题 , 动 隔振是 控制 振动 最 普遍 的方 法之 一 。 被
( L s f es h rc r t s r ai tdb m tac mpr o fh ert a a d x ei narsl . HS D)t f s aat i i e l ae y:u lo ai no eh oei l n p r in c e sc a v d u s t t c e metleut s
理论力学中的非线性静力分析方法是什么?
理论力学中的非线性静力分析方法是什么?在理论力学的广袤领域中,非线性静力分析方法是一个至关重要却又颇具挑战性的研究方向。
当我们面对复杂的力学系统和结构时,线性分析往往无法准确捕捉其真实的力学行为,这时非线性静力分析方法就成为了我们揭示真相的有力工具。
那么,到底什么是非线性静力分析方法呢?简单来说,它是一种用于研究物体或结构在静力作用下,其力学响应呈现非线性特征时的分析手段。
要理解非线性静力分析方法,首先得明白什么是“非线性”。
在力学中,如果一个系统的响应与施加的荷载不成正比关系,那么我们就称其为非线性系统。
比如说,一根弹簧,在其弹性限度内,拉伸的长度与施加的力成正比,这就是线性关系。
但当拉伸超过弹性限度,弹簧的变形就不再与力成正比,这就是非线性了。
非线性静力分析方法与线性静力分析方法最大的不同在于,它需要考虑更多复杂的因素。
在线性分析中,我们可以使用简单的数学公式和方程来描述系统的行为,但在非线性分析中,这些简单的公式就不再适用了。
非线性静力分析方法的应用场景非常广泛。
在建筑结构设计中,比如大型桥梁、高层建筑等,材料的非线性特性(如混凝土的开裂、钢筋的屈服等)会对结构的承载能力和稳定性产生重要影响。
通过非线性静力分析,工程师可以更准确地评估结构在极端荷载下的性能,确保其安全性。
在机械工程领域,例如汽车零部件的设计,发动机部件在高温、高压等恶劣条件下的变形和应力分布往往是非线性的。
运用非线性静力分析方法,可以优化零部件的设计,提高其可靠性和使用寿命。
在航空航天领域,飞机机翼在气流作用下的变形、火箭发动机喷管在高温燃气作用下的热膨胀等,都需要通过非线性静力分析来精确预测和设计。
实现非线性静力分析的方法有多种。
其中一种常见的是材料非线性分析。
这意味着要考虑材料本身的非线性特性,比如材料的塑性变形、蠕变、损伤等。
通过建立合适的材料本构模型,来描述材料在不同应力和应变状态下的力学行为。
几何非线性分析也是重要的一类。
静力非线性分析概论
2.6 非线性静力分析的主要目的通过静力非线性分析主要可以得到:(1) 得到结构在水平荷载作用时内力和变形的全过程,得到结构的最大承载能力和极限变形能力,即基底最大剪力-顶点最大位移,包括楼层位移、层间位移角和顶点位移等重要指标,可以估计相对与设计荷载而言的结构承载力的安全储备大小。
(2) 得到结构首次出现塑性铰位置和各个阶段的塑性铰出现先后顺序和对应的分布状态,可以判断结构是否符合强柱弱梁、强剪弱弯等设计要求,根据塑性铰可以判断结构的薄弱层。
(3) 得到结构在不同的受力阶段下楼层侧移和层间位移角沿高度的分布,校核最大层间位移角是否满足相关规范中规定的结构在多遇地震和罕遇地震作用下的最大层间位移角限值,并结合塑性铰的分布情况可以检查是否存在薄弱层。
(4) 得到结构在不同受力阶段时各部分塑性内力重分布的情况,结合塑性铰先后出现的顺序及分布,进而检查设计的多道设防意图是否能够实现。
(5) 得到结构在每一层的层剪力-层间位移角曲线,可以作为弹塑性层模型时程分析需要的各层等效层刚度。
(6) 得到结构总承载力-顶点位移曲线,它综合表示结构在各个受力阶段的能力和性能。
经过转换,把弹塑性分析得到的结构性能曲线称为“能力曲线”。
可与“需求曲线”进行比较,地震反应谱曲线可以作为“需求曲线“。
如果“需求曲线”和“能力曲线”有交点,则表示结构可以抵抗该地震,交点称为“结构性能表现点”,交点对应的位移就是结构在地震作用下的顶点位移,对该位移处结构的各项性能进行分析,可以得到结构在地震作用下的表现。
如果“需求曲线”和“能力曲线”没有交点,则表示该结构抵抗地震的能力不足,结构不合格。
参考文献:西安科技大学(硕士学位论文)==人字形支撑钢框架结构弹塑性受力性能分析。
千米级斜拉桥汽车荷载作用下非线性静力响应研究
元; 因此研究千米级斜拉桥的非线性影响对以后 类似桥 梁的设计建设有非常重要的意义。分析 汽车荷 载作 用 下千 米级 斜拉 桥 的静 力 非 线性 响 应 , 出了斜拉 桥 主 梁、 得 主塔 和 拉 索这 三 个 不 同
构件 汽 车荷 载作 用 下非线 性静 力响 应的 大 小; 并进 一步 分析 了跨度 增 大对 非线性 响应 的影 响 。
利 的加载位 置后 , 中荷 载乘 以影 响线 值 以及均 布 荷 载乘 以影 响 线 面积 , 集 即可 得 到 活 载效 应 。计 算 影 响 线时 , 如果结 构 的 刚度 只考虑 构件 的弹性 刚度 , 称之 为活 载 效应 的线性 理 论 ; 则 如果 结构 的刚度 还计 入 了
恒 载状 态 时的几 何 刚度 , 则称 之为 活 载效 应 的线 性 二 阶 理论 。线 性 二 阶 理论 忽 略 了活 载 引起 的构形 、 内 力 的改变 , 以叠加 原理 仍然适 用 所 。
关 键词 : 斜拉 桥 ; 千米 级 ; 非线性 ; 车荷 载 汽
中图分 类号 : 4 8 2 文献标 识码 : U 4 .7 A 文 章编 号 :17 0 0 (0 9 0 0 1 0 6 4— 3 0 2 0 ) 4— 0 2— 6
1 引言
斜拉 桥 是塔 、 、 梁 拉索 三者 构成 的复 杂 的超 静定 结 构 体 系 , 着 斜 拉桥 跨 径 的增 大 , 构 的 刚度 急 速 随 结
进行 活载几 何非 线性 分析 时 , 遇 到 以下两方 面 的 问题 : 会 一是 基 于线 性 叠加 原理 的影 响线 方法 失 效 ;
二是 影响线 的形 状 和范 围与结 构 的构 形 、 内力 相关 , 计 算 活载 效 应 又 必 须事 先 已知 影 响 线 。考虑 到 结 而
千米级斜拉桥主梁挠度非线性随机静力分析
吕
1 .
嘣
÷ H
构 失效 概率 , 常用来 检验 其 它方法 的精 确性 响应面 法通 过 研究 结 构 响应 与结 构 各 种 随机 因素 之 间 的关 系来 构造 结构 响应 面方 程 , 用 响应 面方程 即可 进一 步研 究千 米级斜 拉 桥 的结 构 随机特 性 。 运
1 响 应 面 法
大 型桥梁 结构 由于 具有很 强 的非线 性 特 征 , 容易 获 得 其极 限状 态 方 程 的 具体 表 达 式 , 不 因而 无 法直 接采 用一 次二 阶矩 法 ( O M) F R 计算 其可 靠指 标 。采用 不含 交叉项 的二次 多项 式来表 示 响应 面函数 。 。
中 图分 类 号 :4 82 文献 标识码 :A 文章 编号 : 0 5—33 2 1 )2—0 1 0 U 4 .7 2 9 07 (0 10 00 —5
0 概 述
随着苏 通长 江大 桥和 香港 昂船 洲大桥 的建 成通 车 , 拉桥 已进入 千 米级 领 域 。随着 斜拉 桥 跨 径 的增 斜
足 此 目标 可靠 指标 要求 。
表 1 苏 通 长 江 大 桥 随 机 输 入 变 量 的 统 计 特 征
表 2 可靠指标计算结果对 比
2 4 可 靠指标 敏感 因素 分析 .
可靠 指标对 随机 变 量平均 值 和标准差 的敏感度 , 体计 算公 式为 具
0/ o / 0" 一 3 卢
非线性基本概念
非线性概述
... 用线性求解器求解非线性
• 如果初始构形在收敛半径外部, 如果初始构形在收敛半径外部, 帮助获得收敛解: 帮助获得收敛解:
有两种技术可
F F1 ustart u
F
ustart
u
递增加载使目标更接近初始点
用收敛增强工具扩大收敛半径
• 通常结合两种 策略获得收敛. 策略获得收敛.
ANSYS TRAINING
非线性分析
• 在非线性分析中, 响应不能直接用一组线性方程预报,如 在非线性分析中, 响应不能直接用一组线性方程预报, 何 用线性方程近似非线性响应呢? 用线性方程近似非线性响应呢? : 算法: • ANSYS运用 运用Newton-Raphson 算法 运用
– 以增量形式逐渐施加载荷。 以增量形式逐渐施加载荷。 – 在每一载荷增量中完成平衡迭代来使得增量求解达到平衡。 在每一载荷增量中完成平衡迭代来使得增量求解达到平衡。 – 求解平衡方程[KT]{∆u} = {F} - {Fnr} 求解平衡方程[K ]{∆ [KT] =切线刚度矩阵 [KT] {∆u} =位侈增量 F
ANSYS TRAINING
什么是非线性? 什么是非线性?
• 17世纪, 罗伯特•虎克发现力(F) 和位移(u) 之间 17世纪, 罗伯特•虎克发现力(F) 和位移(u) 世纪 存在一个简单的线性关系, 存在一个简单的线性关系, 称为虎克定律: F = Ku
–常数 K为结构的刚度. 常数 为结构的刚度.
非线性概述
ANSYS TRAINING
非线性概述
什么是 “非线性” 行为? 非线性” 行为?
• 如果载荷引起刚度的显著改变, 此结构就是非线性 的. 如果载荷引起刚度的显著改变, • 刚度改变的典型原因是: 刚度改变的典型原因是:
静力计算手册
静力计算手册一、引言静力计算是结构工程中的一项基本技术,主要用于分析结构在恒定载荷作用下的响应。
本手册将详细介绍静力计算的相关理论、方法和应用,旨在为工程师和研究者提供实用的参考。
二、静力学基础1.静力学的概念:静力学主要研究物体在恒定外力作用下,如何维持平衡状态。
2.力的合成与分解:介绍了力的矢量表示方法,以及力的合成与分解的基本原理。
3.力的矩:介绍了力矩的概念,以及如何计算力对点之矩。
4.平衡条件:阐述了静力学中的平衡条件,包括力的平衡和力矩平衡。
三、结构分析方法1.结构分析概述:介绍了结构分析的基本概念、目的和方法。
2.线性静力分析:阐述了线性静力分析的基本原理和步骤,包括建立模型、施加载荷和求解等。
3.非线性静力分析:介绍了非线性静力分析的特点、常见问题和解决方法。
4.有限元法:详细介绍了有限元法的基本原理、实施步骤和应用。
四、杆件内力计算1.杆件内力的计算:介绍了如何通过平衡条件计算杆件的内力。
2.弯矩和剪力:阐述了弯矩和剪力的概念,以及如何计算它们。
3.轴力和扭矩:介绍了轴力和扭矩的概念,以及如何计算它们。
五、截面承载能力计算1.截面几何特性:介绍了截面的面积、惯性矩等基本几何特性。
2.截面承载能力分析:分析了各种形式截面的承载能力,如抗弯、抗剪、抗扭等。
3.材料的强度指标:介绍了常用材料的抗拉、抗压、抗剪强度等指标。
六、结构位移计算1.结构位移的计算方法:介绍了结构位移的几种主要计算方法,如虚功原理、单位载荷法等。
2.常见结构的位移特性:分析了常见结构的位移特性,如梁、柱、拱等。
3.位移与内力的关系:阐述了位移与内力之间的关系,以及如何通过位移计算内力。
七、结构优化设计1.结构优化设计概述:介绍了结构优化设计的基本概念、方法和应用领域。
2.结构优化设计的数学模型:阐述了结构优化设计的数学模型,包括设计变量、约束条件和目标函数等。
3.结构优化算法:介绍了常见的结构优化算法,如梯度法、遗传算法等。
静力非线性方法用于桥梁高墩抗震性能分析
Vo_ 4NO 5 l . 3
非 线性 静力分 析 ( u h v r 法是 一 种简 化 的非 P s o e)
应 尤其 是非 线性 地 震 反 应 的简 单 而 有 效 的方 法 , 因 此 P s o e 分 析法 为各 国规 范所 接 受 , uh vr 中国建筑 结 构 抗震 设计 规范也 引入 了 P s o e 分 析方法 j uh v r 。
MP 法 。 A
关 键词 : 高墩 ; 震 性 能 ; 力 非 线 性 分 析 ; 效 振 型 高 度 抗 静 等 中 图 分 类 号 : U3 2 1 U4 8 2 T 5 . ; 4. 2 文献 标志码 : A 文 章 编 号 :6 44 6 ( 0 2 o —0 20 1 7 — 7 4 2 1 )50 4 —8
LIZh n y n , DANG n p n e g ig。。, Pe g e g。。 LI igt o。 n a J
( . a u t fCi i En i e rn b Ke a o a o y o w c n l g o o tu t n o te n M o n a n Ar a a F c l o v l gn e ig; . y L b r t r fNe Te h oo y f r C s r c i fCiis i u t i e , y o
Oc . 2 2 t 01
d i 1 . 9 9 j is . 6 44 6 . 0 2 0 . 0 o :0 3 6 /.s n 1 7 — 7 4 2 1 . 5 0 7
地震作用计算的方法及各自的使用范围
地震作用计算的方法及各自的使用范围1.引言地震是地球上常见的自然灾害之一,对人类社会和基础设施造成了严重的破坏。
为了准确预测和评估地震对结构物的影响,地震作用计算方法至关重要。
本文将介绍几种常见的地震作用计算方法,并详细阐述它们各自的使用范围。
2.位移法位移法是一种简化的地震作用计算方法,通过假设结构物在地震作用下发生弹性变形,计算结构体的位移响应。
该方法适用于小型结构和较小地震作用的情况,如住宅、小型商业建筑等。
然而,在大震和长周期地震作用下,位移法的精度会降低,因为它无法考虑非线性效应和耗散力的影响。
3.非线性静力法非线性静力法是一种考虑结构物非弹性变形的地震作用计算方法。
该方法通过采用非线性弹簧模型或塑性铰模型,对结构体的产生的非线性效应进行建模,从而计算结构体的应力和变形响应。
非线性静力法适用于中小型结构,可以更准确地预测和评估结构体在地震作用下的性能。
4.动力时程分析法动力时程分析法是一种基于结构体惯性力和地震激励之间相互作用的地震作用计算方法。
通过将结构体建模为质点体系,并考虑结构体和地震作用之间的相互作用力,该方法可以模拟结构体在地震波荷载下的真实动态响应。
动力时程分析法适用于大型或特殊结构,如桥梁、高层建筑等。
5.响应谱分析法响应谱分析法是一种将地震波和结构体的频率特性结合起来,评估结构体在地震作用下的响应的方法。
该方法通过使用结构体的频响函数和地震波的谱函数,计算结构体的响应谱曲线,从而评估结构体的抗震性能。
响应谱分析法广泛应用于工程设计和结构性能评估。
6.使用范围比较不同的地震作用计算方法适用于不同的结构类型和地震作用水平。
以下是各种方法的使用范围比较:-位移法:适用于小型结构和较小地震作用,计算简便,精度相对较低。
-非线性静力法:适用于中小型结构,可以考虑非线性效应,具有较高的精度。
-动力时程分析法:适用于大型或特殊结构,可以模拟真实的动态响应,精度高。
-响应谱分析法:广泛适用于各种结构类型,通过结构体的频率特性评估抗震性能。
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
目的是残差迭代至足够小, 即直到解收敛. 若收敛, 则解在允许的容差范围内平衡.
{F}
{Fnr}
Space Structure Research Center, HIT, CHINA
7 / 27
... 用线性求解器求解非线性
但 Newton-Raphson 法不能保证在所有情况下都收敛!
常数 K 代表结构刚度.
线性结构服从这种线性关系. 常见的例子是简单
弹簧:
F F
K
u
K
u
线性结构非常适合基于线性矩阵代数的有限元 分析方法.
Space Structure Research Center, HIT, CHINA
3 / 27
... 什么是非线性?
然而,当前无法显示此图像。
17 / 27
… 非线性 FEA难点
一个典型的非线性分析包括以下内容:
外载荷 载荷步 (LS) 2
LS 1
一个或更多的载荷步来施加外部 载荷以及边界条件。(这对于线 性分析也同样适用。)
多个子步来逐渐施加载荷。每一 子步代表一个载荷增量。(线性 分析每一载荷步仅需要一个子步)
子步 “时间"
u 位移
Space Structure Research Center, HIT, CHINA
5 / 27
… 用线性求解器求解非线性
Newton-Raphson 法用下列方程迭代求解: [KT]{u} = {F} - {Fnr}
式中:
[KT] = 切向刚度矩阵 {u} = 位移增量 {F} = 外载荷矢量
时间步长决定了在一个子步中的载荷增量。时间步长越大, 载荷增量越大,因此时间步长对求解的精度有直接的影响。
ANSYS具有一个自动时间步长的功能,它会在一个载荷 步的所有子步中预测并控制时间步长
F 外载荷
t
1.0
2.0
“时间“
Space Structure Research Center, HIT, CHINA
并不意味着实际的时间。
外载荷
缺省情况下,在第一载荷步 末赋与时间为1.0,在第二载 荷步末赋与时间为2.0,以此 类推。
1.0
2.0
“时间"
Space Structure Research Center, HIT, CHINA
19 / 27
… 非线性 FEA难点
两个子步间的时间增量是时间步长t
4 / 27
B.用线性求解器求解非线性
非线性分析中, 响应不能用一组线性方程式直接预报. 然而, 非线性结构可以用具有修正的线性逼近迭代系列进
行分析.
ANSYS 采用 Newton-Raphson方法 迭代. 每一次迭代
称为平衡迭代.
载荷
F 34
2
1
一个载荷增量的全 NewtonRaphson 迭代分析. (图所示 为四次迭代)
小或大变形
时间和时间步或子步 数
输出控制
施加载荷。 保存数据库。 求解。
Space Structure Research Center, HIT, CHINA
25 / 27
… 非线性分析基本过程
非线性建模有什么不同之处?
问题的属性影响非线性模型和线性模型的区别:
有的情况下没有区别! 有的情况必须包括特殊的非线性特征:
材料数据 (如塑性应力-应变数据). 特殊单元 (如接触单元).
Space Structure Research Center, HIT, CHINA
几何(大应变)
状态 (接触 / 不接触)
材料 (塑性变形)
铝挤压
Space Structure Research Center, HIT, CHINA
12 / 27
… 三类非线性
几何非线性
如果结构经历大变形, 则几何构形的改变会导致非线性行为.
该例中, 在轻微横向载荷作用下, 杆的端部是柔性的. 载荷增加时, 杆的几 何形状改变 (变弯), 力臂减小 (载荷移动), 引起杆的刚化响应.
F
载荷 收敛
F
U初始
u 位移
初始点在收敛半径外部
U初始u 位移
初始点在收敛半径内部
Space Structure Research Center, HIT, CHINA
9 / 27
... 用线性求解器求解非线性
如果初始构形在收敛半径外部, 有两种技术可帮助获得收 敛解:
F
F
F1
ustart
u
递增加载使目标更接近初始点
采用渐变加载. 采用小的时间步.
Space Structure Research Center, HIT, CHINA
11 / 27
C. 三类非线性
结构非线性由很多原因引起, 可归为三大类:
几何非线性 (大应变, 大挠度, 应力刚化) 材料非线性 (塑性, 超弹性, 蠕变) 状态改变非线性 (接触, 单元死活)
20 / 27
E. 非线性分析的基本过程
和线性分析一样, 非线性分析的基本过程包括三 个主要步骤:
建模 求解 检查结果
本节主要讨论非线性分析的特殊细节— 非线性
分析的不同 之处.
Space Structure Research Center, HIT, CHINA
21 / 27
高级接触和螺栓预紧培训手册
高级结构非线性培训手册
Space Structure Research Center, HIT, CHINA
2 / 27
A. 什么是非线性?
17世纪, Robert Hooke 发现了力(F)和位移(u)
之间的简单线性关系, 即 Hooke 定律: F = Ku
但无法确定收敛半径! 若解收敛, 初始点在收敛半径内. 若不收敛, 初始点在收敛半径外. 需要试错法获得收敛. 经验和训练可减少试错工作.
复杂的问题获得收敛需要多个载荷增量, 每个载荷增 量需多次迭代.
多次迭代会增加总的求解时间.
Space Structure Research Center, HIT, CHINA
Space Structure Research Center, HIT, CHINA
23 / 27
… 非线性分析基本过程
非线性结果分析有什么不同之处?
非线性结果能得到线性分析中没有的结果项. 例如:
塑性应力和应变, 接触状态等.
由于多次增量求解使非线性分析产生响应历史.
响应历史曲线图
15 / 27
D. 非线性 FEA难点
非线性有限元分析带来三个主要难点:
获得收敛 权衡代价与精度 验证
成功解决三个难点需要 小心和技巧!
Space Structure Research Center, HIT, CHINA
16 / 27
… 非线性 FEA难点
克服收敛困难
通常是最大的挑战. 求解必须在收敛半径内开始.
12. ANSYS软件非线性分析
非线性概述
当作用在结构上的载荷引起结构刚度的重大改变时,要进 行非线性分析。引起刚度改变的主要原因是:
应变超过弹性极限(塑性) 大变形,如受力的钓鱼杆 两个物体之间的接触
Stress
Strain
Space Structure Research Center, HIT, CHINA
22 / 27
… 非线性分析基本过程
非线性求解有什么不同之处?
多次矩阵解:
线性静态仅需通过矩阵方程求解器一次. 非线性每次迭代都进行一次新的求解.
求解控制选项:
线性静态一般不需用户控制求解选项. 非线性用求解控制选项来:
激活几何非线性. 改善收敛性. 管理结果文件大小. 优化重启动能力.
Space Structure Research Center, HIT, CHINA
24 / 27
… 非线性分析基本过程
假定几何模型与网格划分已完成,非线性分析的典型步骤 如下: 1. 确定分析类型(通常为静态)。 2. 确定求解控制Solution>Sol’n Control. 许多控制可以采 用,但是常用的是:
1 / 27
…非线性概述
本章通过以下主题介绍非线性有限元分析(FEA) 的基础:
A. 什么是非线性? B. 用线性求解器求解非线性 C. 三类非线性 D. 非线性 FEA 难点 E. 一般的非线性分析过程
目的是了解非线性 FEA 的基本性质. 如何模拟各 种非线性现象的更多细节参考下面的两本书中:
Space Structure Research Center, HIT, CHINA
13 / 27
… 三类非线性
材料非线性
非线性的应力-应变关系是非线性结构行为的常见原因.
应力
应力
应变
钢
橡胶
Space Structure Research Center, HIT, CHINA
应变
14 / 27
ustart
u
用收敛增强工具扩大收敛半径
通常结合两种策略获得收敛.
Space Structure Research Center, HIT, CHINA
10 / 27
... 用线性求解器求解非线性
一般的规律是系统任何方面的突变会导致收敛困难.
刚度突变. 载荷突变.
最佳收敛行为是把突变分成一系列很多小的递增的变化.
[KT1]
F
34
Fnr1
2
1