结构非线性作业参考-方志老师
结构非线性分析
结构非线性分析理论1.结构设计方法结构设计方法从传统的容许应力设计法发展到了基于概率统计的极限状态设计法。
传统的容许应力设计法是基于线弹性理论，依照经验选取一定的安全系数，以构件危险截面某一点的计算应力不超过材料的容许应力为准则，目前在某些领域仍在使用。
安全系数，是一个单一的根据经验确定的数值，没有考虑不同结构之间的差异，不能保证不同结构具有同等的安全水平。
此外，容许应力设计法以弹性理论计算内力，对那些发展塑性变形能提高承载力的构件或结构(如受弯构件)，比那些发展塑性变形不能提高承载力的构件或结构(如轴心受力构件)具有较大的安全储备。
概率极限状态设计法是采用数理统计方法按照一定概率确定荷载或材料的代表值，并给出结构的功能函数，用结构失效概率或可靠指标度量结构的可靠性。
《建筑结构可靠度设计统一标准》将极限状态分为两类：(1)承载能力极限状态，是指结构或结构构件达到最大承载能力或不适于继续承载的变形；(2)正常使用极限状态，是指结构或结构构件达到正常使用或耐久性能的某项规定限值。
结构按极限状态设计应符合下列要求：()0,21≥n X X X g （1.1）式((1.1)中g(X i )为结构功能函数，X i (i =1， 2……n)为基本变量，是指影响该结构功能的各种作用、材料性能、几何参数等。
目前我国结构设计规范基本都是采用以概率理论为基础的极限状态设计方法，用分项系数设计表达式进行计算。
美国的钢结构设计采用了两种设计方法：ASD(Allowable Stress Design)和LRFD(Load and Resistance Factor Design)，即容许应力设计法和分项系数设计法，McCormac 指出LRFD 相比ASD ，并不一定节省材料，虽然在很多情况下可以取得这样的效果，而在不同荷载作用下能给结构提供等同的可靠性，对于活载和恒载，ASD 采用的安全系数是一样的，而LRFD 对恒载则采用了一个较小的荷载系数(恒载比活载能更准确的确定)，也就是说如果恒载大于活载，LRFD 比ASD 节省材料。
第2章-结构几何非线性
尤拉描述法（Eulerian Formulatlon） 独立变量是质点P当前时刻的位置xn+1与时 间tn+1。 现在用得最为广泛的是Lagriangian列式， 因此下面主要讲述T.L和U.L下几何非线性 有限元方程的建立。
第2章
6
三. 线弹性杆单元刚度矩阵建立的步骤 1.选择插值函数描述体内任一点的位移 u
i xi ui
第2章 16
Lagrange 描述－Green应变张量
x3
Q( xi dxi )
Q '(i di )
ds
P' (i )
ds0
P( xi )
O
x2
x1
第2章 17
图示两相邻点P和Q，变形前其坐标分别为xi 和xi dxi，
单元PQ的长度为ds0；变形后两点分别变位至P ' 和Q '， 其坐标分别为i 和i di，P ' Q '的长度变为ds。点P的 位移向量为ui。 很明显，若ds ds0，则单元没有变形，仅发生刚体平
第2章 20
若计轴x1 , x2 , x3分别对应x, y, z轴，u , v, w分别 为x, y, z方向的位移分量，则有工程符号表达 的 xx 和 xy为： u 1 u 2 v 2 w 2 xx 11 ( ) ( ) ( ) x 2 x x x u v u u v v w w xy 2 xy 212 ( ) y x x y x y x y
第2章
4
全拉格朗日列式法（ T.L列式法－ Total Lagrangian Formulation）。选取to=0时刻 未变形物体的构形Ao作为参照构形进行分 析。 修正的拉格朗日列式法（ U.L列式法－ Updated Lagrangian Formulation）。选取 tn时刻的物体构形An为参照构形。由于An随 计算而变化，因此其构形和坐标值也是变 化的，即与t有关。tn为非线性增量求解时 增量步的开始时刻。
03非线性分析要点
第三部分非线性分析第一章非线性有限元概述1.1非线性行为1、 非线性结构的基本特征是结构刚度随载荷的改变而变化。
如果绘制一个非线 性结构的载荷一位移曲线，则 力与位移的关系是非线性函数。
2、 引起结构非线性的原因：a 几何非线性：大应变，大位移，大旋转 （例如钓鱼竿的变形）b 材料非线性：塑性，超弹性，粘弹性，蠕变c 状态改变非线性：接触，单元死活3、 非线性行为一一分析方法特点A 不能使用叠加原理！B 结构响应与路径有关，也就是说加载的顺序可能是重要的。
C 结构响应与施加的载荷可能不成比例。
1.2非线性分析的应用1、 一些典型的非线性分析的应用包括： 非线性屈曲失稳分析金属成形研究碰撞与冲击分析制造过程分析（装配、部件接触等）材料非线性分析 （塑性材料、聚合物）2、 橡胶底密封：一个包含几何非线性（大应变与大变形），材料非线性（橡胶）， 及状态非线性（接触）的例子。
2.1非线性方程组的解法1、求解一个结构的平衡问题通常等于求解结构的总位能的驻值 问题。
结构总位能n ： 口 "3弋门心 2、 增量法：就是将荷载分成一系列的荷载增量，即 ANSYS 中的荷载步或荷载子 步。
A 要点：在每一个荷载增量求解完成后，继续进行下一个荷载增量之前， 刚度矩阵以反映结构刚度的变化。
B 增量法的优点：可以追踪结构变形历程，这对于材料或几何非线性（特别是 极限值屈曲分析）十分有用。
C 增量法的缺点：随着荷载步增量的增加而产生积累误差，导致荷载-位移曲 线飘移。
D 对飘移进行平衡修正，可以大大提高增量法的精度。
应用最广的就是在每一 级载荷增量上用Newton-Raphsor 或其变形的迭代法。
3、 迭代法：割线刚度法：收敛性差，因此很少应用切线刚度法Newto n-Ra phsor 迭代法：切向刚度法中 2.2 Newto n-Ra phsor 迭代法 1、 优点：对于一致的切向刚度矩阵有 二次收敛速度。
结构非线性
10.1 非线性问题分类及求解 10.2 非线性问题求解方法 10.3 材料非线性 10.4 几何非线性 10.5 边界非线性 10.6 非线性弹性稳定性问题 10.7非线性分析特点 10.8 ANSYS非线性结构计算示例 10.9ANSYS稳定性计算示例
第十章 结构非线性分析的有限单元法简介
图10-3 N—R迭代法的几何意义
图10-4 修正牛顿法迭代几何意义
第十章 结构非线性分析的有限单元法简介
10.2.3 载荷增量法
, K T P 0
为载荷因子，用来描述载荷变化的参数， 对应于 ， 对应于 ，则 , 0
第十章 结构非线性分析的有限单元法简介
10.4.2 几何非线性有限元分析
由虚功原理
T * e
F dxdydz dV
t
Dt mD 1 mDP
0＜m＜1
其中，m为加权因子．当m＝1时为完全弹性；m＝0 为完全塑性。m值的物理意义见图10-7。
图10-7 m值的物理意义
第十章 结构非线性分析的有限单元法简介
(2) 弹塑性有限元解法
弹塑性问题求解常用切线刚度法、初应力法或切线 刚度法等增量法。 同样，弹塑性问题的平衡方程可以表示为
式中
s
——屈服应力，
H ——塑性强化模量。
岩石在承受 的荷载超过 一定值时， 如较高的围 岩压力时表 现出理想塑 性特性。
第十章 结构非线性分析的有限单元法简介
弹塑性变形时总应变包括 两部分。
弹塑 性 元件 应力足够大 式中 e ——弹性应变， 时的金属、 岩石、土壤。 p ——塑性应变。 加载时使用增量理论。
结构力学模拟中的三类非线性问题
1. 线性分析外加载荷与系统的响应之间为线性关系。
例如线性弹簧，结构的柔度阵（将刚度阵集成并求逆）只需计算一次。
通过将新的载荷向量乘以刚度阵的逆，可得到结构对其它载荷情况的线性响应。
此外，结构对各种载荷情况的响应，可以用常数放大和/或相互叠加，以确定它对一种全新载荷情况的响应，所提供的新载荷情况是前面各种载荷的叠加（或相乘）。
这种载荷的叠加原理假定所有的载荷情况采用了相同的边界条件。
2. 非线性分析非线性结构问题是指结构的刚度随其变形而改变。
所有的物理结果均是非线性的。
线性分析只是一种近似，它对设计来说通常已经足够了。
但是，对于许多结构包括加工过程的模拟（诸如锻造或者冲压）、碰撞分析以及橡胶部件的分析（诸如轮胎或者发动机支座），线性分析是不够的。
一个简单例子就是具有非线性刚度响应的弹簧。
线性弹簧，刚度是常数非线性弹簧，刚度不是常数由于刚度依赖于位移，所以不能再用初始柔度乘以外加载荷的方法来计算任意载荷时弹簧的位移。
在非线性隐式分析中，结构的刚度阵在整个分析过程中必须进行许多次的生成和求逆，分析求解的成本比线性隐式分析昂贵得多。
在显式分析中，非线性分析增加的成本是由于稳定时间增量减小而造成的。
非线性系统的响应不是所施加载荷的线性函数，因此不能通过叠加来获得不同载荷情况的解答。
每种载荷情况都必须作为独立的分析进行定义和求解。
3. 非线性的来源在结构的力学模拟中有三种：材料非线性、边界非线性（接触）、几何非线性。
(1) 材料非线性大多数金属在低应变值时都具有良好的线性应力/应变关系；但是在高应变时材料发生屈服，此时材料的响应成为了非线性和不可恢复的。
橡胶材料等也是一种非线性、可恢复（弹性）响应的材料。
材料的非线性也可能与应变以外的其它因素有关。
应变率相关材料数据和材料失效都是材料非线性的形式。
材料性质也可以是温度和其它预先定义的场变量的函数。
(2) 边界非线性如果边界条件在分析过程中发生变化，就会产生边界非线性问题。
非线性本构理论及方程
非线性本构理论及方程非线性本构理论及方程是构成工程力学和材料科学的重要组成部分，它反映了物质的力学特性，是了解材料的自然行为的关键概念。
本文将介绍非线性本构理论及其相关方程，包括非线性本构模型、非线性本构方程、压缩圆柱模型、等因式能量函数等。
首先，介绍非线性本构模型。
非线性本构模型是描述材料性质的基本概念，它涉及材料物理本质，模型可以用来研究材料在加载过程中的全局响应，以及材料力学和结构力学性质。
常见的非线性本构模型有弹性-塑性模型、扭转模型、粘弹性模型等。
其次，介绍非线性本构方程。
非线性本构方程是描述材料性质的基本方程，它涉及材料物理本质，可以用来研究材料在加载过程中响应的性质和行为规律。
常见的非线性本构方程有Jaumann函数、等因式能量函数、Rice-Salamon函数等。
再次，介绍压缩圆柱模型。
压缩圆柱模型是用来描述材料性质的一种模型，它是一种压缩材料的流变特性模型，可以用来描述材料在压缩方向的性质，同时也可以用来分析材料的非线性行为。
压缩圆柱模型的一般形式为：σ=K_0*[1+e~(-K~2*ε)]^(-n)其中，K_0是已知的参数，e~(-K~2*ε)是可以计算的，n是未知的参数，σ是应力，ε是压缩应变。
最后，介绍等因式能量函数。
等因式能量函数是用来描述材料性质的常用方程，它是建立材料屈服条件的重要函数，可以用来表征材料在上下线性段之间的行为规律。
等因式能量函数的一般形式为：W=K_1ε^2*(1+K_2ε^n)其中，K_1、K_2和n是未知参数，W是能量，ε是应变。
综上所述，非线性本构理论及其相关方程是工程力学和材料科学的重要组成部分，它反映了物质的力学特性，是了解材料的自然行为的关键概念。
本文介绍了非线性本构模型、非线性本构方程、压缩圆柱模型、等因式能量函数等。
将本构理论和方程应用到工程设计中，将有助于更好地使用材料以解决工程问题。
5 桥梁结构的材料几何非线性分析
几何非线性理论一般可分为大位移小应变即有限位移理论和大 位移大应变即有限应变理论两种。
桥梁工程中的几何非线性问题一般都是有限位移问题。 一些简单几何非线性问题可以找到解析解，如压弯杆稳定问题
，拱圈刚度按一定规律变化的拱桥大挠度问题，悬索桥在简单荷载 作用下的大挠度问题等。
(1) 材料非线性问题的平衡方程
以钢材和混凝土为主要材料的桥梁结构，所涉及的材料非线性主 要是弹塑性问题。
以有限元分析桥梁结构时，所建立的平衡方程为
[B]T { }dV {F}
{ } [B]{ }
由于并未放弃小变形假定，对桥梁的材料非线性问题，上列两式仍
然成立，但物理方程是非线性的，可以写成
方程的非线性，则称其为材料非线性问题。如第13章所介绍的混凝 土本构关系中，大多本构模型为非线性模型，必将引起平衡方程的 非线性。
在桥梁工程问题中: 混凝土的徐变、收缩、结构弹塑性等都属于材料非线性问题 桥梁结构中常用的低碳钢在承载力的后期亦进入弹塑性阶段， 呈现出材料非线性本质。
材料非线性问题可以分为非线性弹性问题和弹塑性问题两大类， 前者在卸载后无残余应变存在，后者会存在残余变形。但两者的本 质是相同的，求解方法亦完全一样。
但多数问题还需借助有限元及其它数值法求解 （3）接触问题
若受力后的边界条件在求解前是未知的，即不满足理想约束 假定而引起的边界约束方程的非线性问题称其为接触问题。
如:悬索桥主缆与鞍座的接触状态问题 支架上预应力梁在张拉后的部分落架现象 等均属此类，此问题在桥梁工程上表现不多，但不应忽视。
（4）桥梁结构非线性 材料非线性问题在混凝土桥中表现最为突出，由于混凝土材料
非线性分析实例
ansys非线性分析实例(1) 设定分析参数，在ANSYS顶部菜单Parameters->Scalar Parameters，在弹出的Scalar Parameters窗口中输入，FORCE=100，OFFSET=0.1。
(2) 建立模型，在本算例中，我们将接触两种新的单元，三维梁单元Beam 4和三维索单元Link 10。
Beam 4单元在参数定义等方面比前面介绍的Beam 188单元简单，对于常见的矩形弹性截面，也是一个很实用的单元。
Link 10单位为ANSYS提供的空间索单元，用户可以控制该单元只能受压或者只能受拉，默认该单元只能受拉。
(3) 在ANSYS主菜单Preprocessor->Element type->Add/Edit/Delete中选择Beam 4单元和Link 10单元(4) Beam 4单元和Link 10单元都可以通过实参数来设置截面属性。
首先设置Beam 4单元的截面。
我们要建立的截面是一个0.1m×0.12m的矩形截面。
Beam 4单元如果没有指定截面主轴方向，则截面局部坐标系的Y轴方向将和整体坐标系的X-Y平面平行。
在ANSYS 主菜单Preprocessor->Real Constants->Add/Edit/Delete，选择Add，指定实参数的关联单元类型为Beam 4单元，输入截面参数为AREA: 0.1*0.12, IZZ: 0.12*0.1**3/12, IYY:0.1*0.12**3/12, TKZ: 0.12, TKY: 0.1, 如下图所示(5) 继续添加第二个实参数类型，指定第二个实参数与Link 10单元相关。
取Link 10单元的截面积为4×10-6m2，初始应变为2×10-3，这里正的表示初始应变为拉应变。
(6) 下面设定材料属性，在本例子中为了简单起见，所有材料都设定为钢材。
机械结构的非线性响应分析
机械结构的非线性响应分析随着科技的不断发展，机械结构的设计和分析变得日益重要。
在工程实践中，准确预测机械结构的响应对于保证结构的安全性和可靠性至关重要。
然而，由于机械结构在受力过程中可能会出现非线性响应，传统的线性分析方法已经不能满足实际需求。
因此，研究机械结构的非线性响应分析方法成为了一个热门的课题。
首先，我们来了解一下什么是机械结构的非线性响应。
机械结构的非线性响应是指结构在受到外力作用时，结构的变形与所受力的关系不是简单的比例关系。
这是由于材料的非线性特性导致的。
材料的非线性特性即应力和应变之间的关系不是线性的，而是随着应力的增大而变化的。
因此，在分析机械结构的响应时，需要考虑材料的非线性特性对结构的影响。
非线性响应分析可以分为几类，其中最常见的是几何非线性。
几何非线性指的是在结构变形过程中，结构本身的几何形状发生变化的情况。
例如，当结构受到大的变形作用时，结构的刚度和应力分布会随之发生变化。
这种几何非线性将导致结构的响应与线性分析结果存在差异。
另一类非线性响应分析是材料非线性。
材料非线性是由于实际材料的本构关系不满足线性弹性假设引起的。
在加载过程中，材料的应力和应变之间的关系将会发生变化。
例如在低应力范围内，材料可能表现出线性弹性行为；但当加载超过一定应力阈值时，材料可能呈现出塑性变形或破坏行为。
要进行机械结构的非线性响应分析，首先需要建立合适的数学模型。
数学模型是对于机械结构的抽象表示，可以用来描述结构的受力、变形和相互作用等情况。
建立数学模型需要考虑结构的几何形状、材料的力学性质以及结构的约束条件等影响因素。
常用的非线性分析方法包括有限元法、边界元法和解析法等。
其中最常用的是有限元法。
有限元法是一种将连续体分割成有限个离散单元进行分析的方法。
通过在每个单元上建立适当的数学模型，将问题转化为求解一组代数方程。
然后通过求解代数方程来得到结构的响应。
在进行非线性响应分析时，还需要考虑加载路径的影响。
结构非线性是什么
结构非线性是什么结构非线性，又称为非线性结构或非线性系统，是指系统内部的元素或组成部分之间的关系不遵循线性规律，即输入和输出之间的关系不是简单的比例关系。
相对于线性结构的输出随着输入的变化而以恒定速率变化，非线性结构的输出则可能呈现出复杂的非线性特征。
在现实生活中，大部分系统都具有一定的非线性特征，包括天气系统、交通系统、经济系统、生态系统等。
非线性结构的特点是具有复杂的动力学行为和多样的稳定态，在控制和预测方面具有挑战性。
非线性结构的一个重要特征是存在非线性耦合作用。
线性系统中，不同元素之间的相互作用是简单的叠加效应，而在非线性系统中，不同元素之间的相互作用可能产生非线性的引力或阻力，使得系统的行为变得复杂多样。
另外，非线性结构还表现出以下几个方面的特点：1.非线性响应：当系统受到外部激励时，其输出不是简单的线性响应，可能出现分岔、周期倍增、混沌等现象。
这意味着微小的扰动或变化可能对系统产生不可预测的巨大影响。
2.多重稳定态：非线性结构能够在不同的输入条件下呈现出不同的稳定态。
这意味着系统有可能同时存在多个平衡点，且在不同的输入条件下转变之间可能具有突变或不连续性。
3.尺度依赖性：非线性结构的行为在不同的空间尺度上可能呈现出不同的规律。
这使得我们在分析和预测非线性结构时需要考虑多个尺度的因素，以获得更准确的结果。
4.遗传性：非线性结构的特性可能通过遗传方式传递给后代系统。
这意味着非线性结构具有一定的内在记忆，历史变化对当前和未来行为有一定的影响。
非线性结构的研究对于我们理解和控制复杂系统、优化各种工程和科学问题都具有重要意义。
例如，在天气和气候预测中，非线性结构的天气系统具有混沌特性，微小的初始条件变化可能导致巨大的结果差异。
因此，我们需要开发出具有高分辨率和高精度的模型和方法来预测和控制天气系统。
另外，在经济领域，非线性结构也被广泛应用于金融市场预测、经济波动的研究和货币政策制定等方面。
为了研究非线性结构，科学家和工程师们提出了各种数学模型和方法。
非线性结构
Root
B
A
第一层
D
C
第二层 第三层 第四层
E
F
G
H
I
J
（1）叶子
没有后继的结点称为叶子（或终端端结点），图1.23 中的结点D、E、F、G、H、J为终端结点； （2）分支结点 非叶子结点称为分支结点 （或：非终端结点）。 （3）结点的度 一个结点的子树数目称为该结点的度。B的度为2, 结点C的度为3；D、J的度为0； （4）树的度 树中各结点的度的最大值称为该树的度，上图所示 的树的度为3。
中序遍历的递归定义为：
若二叉树为空，遍历结束，否则： （1）按中序遍历方式遍历根结点的左子树； （2）访问根结点； （3）按中序遍历方式遍历根结点的右子树；
二叉树的中序遍历算法
void inorder (bnode *BT) { if (BT= =NULL) return; else { if (BT->LC ！＝ NULL) inorder (BT->LC); visite (BT); /*访问BT指向的根结点*/ if (BT->RC ！＝ NULL) inorder (BT->RC); } }
§1.3 非线性结构
非线性结构的逻辑特征是一个结点元素可能 有多个直接前趋和多个直接后继。最主要的非线 性结构是树结构和图结构。 1．3．1树结构及其基本概念 树结构是一类重要的非线性结构。树结构是 结点之间有分支、层次关系的结构，在客观世界 中，树结构是大量存在的，例如家谱、行政组织 机构都可用树形象地表示。在计算机领域中，树 结构也被广泛应用，如计算机磁盘文件的管理， 是一种从根目录到各级子目录的分层结构i在数据 库系统中，常采用树来组织数据信息。
例6 已知一棵二叉树的先序遍历序列为ABCDEFGHIJ．中 序遍历序列为CBDEAFHIGJ，试构造这棵二叉树。构造过 程由图示如下： A A
非线性结构的动力特性与自振频率分析
非线性结构的动力特性与自振频率分析非线性结构是指在受力作用下，其应变与应力之间的关系不遵循线性规律的一类结构。
与线性结构相比，非线性结构具有丰富的动力特性和振动行为。
研究非线性结构的动力特性和自振频率，对于工程设计和结构安全分析至关重要。
本文将介绍非线性结构的动力特性和自振频率分析的方法和应用。
一、非线性结构的动力特性分析非线性结构的动力特性是指在受力作用下，结构发生振动时具有的特定性质和行为。
与线性结构相比，非线性结构的动力特性更为复杂，其中包括非线性振动、非线性耗能和非线性共振等现象。
1. 非线性振动非线性振动是指结构在受到激励作用下，产生的振动不符合线性规律。
这种振动可能表现为周期性振动、分岔现象、倍周期振动等。
非线性振动的出现使得结构的动力响应更为丰富，需要通过数值模拟或试验手段来分析和研究。
2. 非线性耗能非线性耗能是指结构在振动过程中由于摩擦、塑性变形等因素导致的能量损耗。
非线性耗能可以有效减小结构的振动幅值，提高结构的抗震性能。
因此，对于非线性结构的动力特性进行分析时，需要考虑非线性耗能的影响。
3. 非线性共振非线性共振是指结构在受到周期激励时，振动频率与激励频率之间存在非线性关系，导致结构响应出现共振放大现象。
非线性共振的出现可能引发结构的失稳和破坏，因此对于非线性结构的动力特性分析应重点研究非线性共振的机理和特征。
二、非线性结构的自振频率分析自振频率是指结构在无外界激励下，由自身固有刚度和质量决定的振动频率。
对于非线性结构的自振频率分析，需要考虑非线性因素对结构刚度的影响。
1. 线性刚度法线性刚度法是非线性结构自振频率分析的一种常用方法，它将非线性结构视为由各个线性小段组成的多自由度系统。
通过将非线性系统离散化为多个线性系统，可以计算出每个分段结构的自振频率，然后将其合并得到整个非线性结构的自振频率。
2. 近似解法对于复杂的非线性结构，无法直接应用线性刚度法进行自振频率分析。
此时，可以使用一些近似解法，如变分法、贝塞尔函数法、有限元法等。
结构非线性分析ABAQUS
《工程结构非线性》作业学院：土木工程学院专业：结构工程姓名：**学号：S********教师：方志（教授）目录作业 (3)1 偏压柱的跨中最大挠度的解析解 (3)2用有限元软件ABAQUS建立题中所给的弯压柱的力学模型，并计算跨中最大挠 (4)2.1 给出一个实例 (4)2.2 确定材料的本构模型 (4)2.3 建立有限元模型 (5)2.4 模拟结果分析对比 (12)3 ABAQUS有限元软件分析的理论背景（来自ABAQUS帮助文件） (14)4 对结构几何非线性和稳定的关系进行讨论 (24)结构非线性作业一（1） 求出荷载——柱中点侧移的解析解。
（2） 以具体的实例给出几何非线性效应得数值解（可用有限元程序计算），并与解析解结果对比。
（3） 给出有限元程序理论背景的详细描述。
（4） 对结构几何非线性和稳定的关系进行讨论。
1 偏压柱的跨中最大挠度的解析解图1 计算简图1.1跨中弯矩为： ()M P e y =+ （1）1.2由材料力学中梁挠曲线的近似微分方程可以得到： 22d y Mdx EI=- 将（1）式代入其中得''()P e y P Pey y EI EI EI+=-=+解微分方程得：[]csc sin ()csc sin sin ()sin csc l x x y e l l x e l x l l e l x x l e ααααααα-⎡⎤⎡⎤=--+-⎢⎥⎢⎥⎣⎦⎣⎦=-+- 其中P EIα=1.3 求跨中侧移：当2lx =时 max2sin csc (sec 1)22l ly e l e e ααα=-=-2 用有限元软件ABAQUS 建立题中所给的弯压柱的力学模型，并计算跨中最大挠度2.1 给出一个实例：假设题中所给弯压柱所受荷载P=10KN, 偏心距e=0.1m ，柱高为L=2m ，采用屈服强度为345MP 的钢材，弹性模量E=206000MP, 柱的截面尺寸如所示：图1 计算简图2.2 确定材料的本构模型采用韩林海（2007）中的二次塑性流模型来模拟钢材，其应力-应变关系曲线,分为弹性段(Oa)、弹塑性段(ab)、塑性段 (bc)、强化段(cd)和二次塑流(de)等五个阶段，如图1所示。
非线性规划例题
1、用梯度法（最速下降法）求下述函数的极小点：解：取初始点T X )0,0()0(=。
TXf )0,0()()1(=∇，故)1(X 为极小点。
其极小值0)()1(=X f 。
2、用梯度法（最速下降法）求函数22215)(x x X f +=的极小点，取允许误差7.0=ε。
解：取初始点TX)1,2()0(=。
⎪⎪⎭⎫ ⎝⎛-=⎪⎪⎭⎫ ⎝⎛-⎪⎪⎭⎫ ⎝⎛==⎪⎪⎭⎫ ⎝⎛⎪⎪⎭⎫ ⎝⎛⎪⎪⎭⎫ ⎝⎛=⎪⎪⎭⎫⎝⎛=∇=∇=∇1240.05504.11041124.0121124.010410002)10,4(104)10,4(10002)()10,4()(,)10,2()()1(02)0(21XX f X f x x X f T T λ。
其海赛矩阵ε>=∇⎪⎪⎭⎫ ⎝⎛-=∇1526.11)(,2400.11008.3)(2)1()1(Xf Xf ελελελ<=∇⎪⎪⎭⎫ ⎝⎛=∇⎪⎪⎭⎫⎝⎛=⎪⎪⎭⎫ ⎝⎛--⎪⎪⎭⎫ ⎝⎛-==⎪⎪⎭⎫ ⎝⎛-⎪⎪⎭⎫ ⎝⎛-⎪⎪⎭⎫ ⎝⎛--=>=∇⎪⎪⎭⎫ ⎝⎛-=∇⎪⎪⎭⎫⎝⎛-=⎪⎪⎭⎫ ⎝⎛-⎪⎪⎭⎫ ⎝⎛==⎪⎪⎭⎫ ⎝⎛⎪⎪⎭⎫ ⎝⎛⎪⎪⎭⎫ ⎝⎛=>=∇⎪⎪⎭⎫ ⎝⎛=∇⎪⎪⎭⎫⎝⎛=⎪⎪⎭⎫ ⎝⎛--⎪⎪⎭⎫ ⎝⎛-==⎪⎪⎭⎫ ⎝⎛-⎪⎪⎭⎫ ⎝⎛-⎪⎪⎭⎫ ⎝⎛--=6685.0)(,759.0304.0)(0759.0152.03419.08542.03221.003419.04271.03221.03419.08542.010002)3419.0,8542.0(3419.08542.0)3419.0,8542.0(8466.0)(,3419.08542.0)(03419.04271.0757.2102.11124.02757.05510.01124.0757.2102.110002)757.2,102.1(757.2102.1)757.2,102.1(815.8)(,757.2102.1)(2757.05510.02400.11008.33223.01240.05504.13223.02400.11008.310002)2400.1,1008.3(2400.11008.3)2400.1,1008.3(2)4()4()4(32)3()3()3(22)2()2()2(1X f X f X X f X f X X f X f X故以TXf )0759.0,152.0()()4(=为近似极小点，此时的函数值0519.0)()4(=X f 。
结构非线性反应PPT
m(t ) c(t )y(t ) k (t )y(t ) P(t ) ----增量方程 y
2.线加速度法
m(t ) c(t )y(t ) k (t )y(t ) P(t ) y
----增量方程
(t ) y(t ) y 2 y (t ) y (t )t (t ) (t )3 2! 3! y(t ) y (t ) (t )t y (t ) 2 2! y ( 4) (t ) (t ) y(t )t y (t ) 2 ~ k (t ) k (t ) 2!
§3.12 结构非弹性地震反应分析
一、结构的非弹性性质 滞回曲线：结构或构件在反复荷载作用下的力与非弹 性变形间的关系曲线。
受弯钢筋凝土构 件的滞回曲线
滞回模型：描述结构或 构件滞回关系的数学模 型。
双线性模型
双线性模型一般适 用于钢结构梁、柱、节 点域构件。 钢筋混凝土梁、柱、 墙等一般采用退化三线 性模型。
1 (ti ) [ P(ti ) f D (ti ) f s (ti )] y y(0) 0.0744 2. t=0.1s m ~ ~ y(0.1) y(0) y(0) 0.00248 y(0.1) y(0) y(0) 0.0744 P(ti )、k （ti ) （2）求 ~ 6 3 y(0.1) 0.05 弹性阶段 k (0.1) 60 k (t ) k (t ) m c(t )
3.计算步骤 已知ti-1时刻的状态向量及 y(ti 1 )、y(ti 1 ) 求ti时刻的状态向量及增量。
（1）求ti时刻的状态向量 y (ti ) y (ti 1 ) y (ti 1 )
y (ti ) y (ti 1 ) y (ti 1 ) (ti ) y
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
《工程结构非线性》作业
学院：土木工程学院
专业：桥梁工程
姓名：刘万事
学号：S140110021
教师：方志（教授）
结构非线性作业一
（1） 求出荷载—柱中点侧移的解析解及第一类失稳荷载；
（2） 以具体的实例给出几何非线性效应的数值解（可用有限元程序计算,但应给出有限元程
序理论背景的详细描述），并与解析解结果对比；
（3） 对结构几何非线性和稳定的关系进行讨论。
1、偏压柱的跨中最大挠度的解析解
图1 计算简图
1.1跨中弯矩为： ()M P e y =+ （1）
1.2由材料力学中梁挠曲线的近似微分方程可以得到： 22d y M
dx EI
=- 将（1）式代入其中得
''()P e y P Pe y y EI EI EI
+=-=+
解微分方程得：
[]csc sin ()csc sin sin ()sin csc l x x y e l l x e l x l l e l x x l e
ααααααα-⎡⎤⎡
⎤=--+-⎢⎥⎢⎥⎣⎦⎣⎦=-+-
其中α=
1.3 求跨中侧移：当2l
x =
时 max
2sin csc (sec 1)22
l l
y e l e e ααα=-=-
2、用有限元软件ansys 建立题中所给的弯压柱的力学模型，并计算跨中最大挠度
2.1 给出一个实例：
假设题中所给弯压柱所受荷载P=10KN, 偏心距e=0.05m ，柱高为L=5m ，采用屈服
强度为345MP 的钢材，弹性模量E=2.06x105
MP, 柱的截面尺寸如所示：
图1 计算截面
2.2 确定材料的本构模型
采用韩林海（2007）中的二次塑性流模型来模拟钢材，
其应力-应变关系曲线,分为弹性段(Oa)、弹塑性段(ab)、塑性段 (bc)、强化段(cd)和二次塑流(de)等五个阶段，如图1所示。
图1中的点划线为钢材实际的应力-应变关系曲线，实线所示为简化的应力-应变关系曲线，模型的数学表达式如式(3-1)。
其中：
e1e3e1e2e e1s y e 100,10,5.1,/8.0εεεεεεε====E f ； f p 、f y 和 f u 分别为钢材的比例极限、屈服极限和抗拉强度极限。
⎪⎪⎪⎪⎪
⎩
⎪⎪⎪⎪⎪
⎨⎧>≤<⎥⎦⎤
⎢⎣
⎡--+≤<≤<++-≤=e3
s y e3s e2e2e3e2s y e2s e1y
e1
s e s 2s e s s s s 6.16.01εεεεεεεεεεεεεεεεεεεεσf f f C B A E (3-1)
由该本构模型计算出材料的应力—应变关系
表1 计算的钢管的力学参数与应力——应变曲线
应
变 分
段
点 0ε
0.001729612
e ε 0.001383689 1
e ε
0.002075534 2e ε 0.02075534 3e ε
0.218143 系 数
A
1.48877E+14 B
6.18E+11
C -285040000
2.3 建立有限元模型
建立一个新的文件，将材料、截面形状、截面特性进行设置，支座为两端铰支，梁划分成100个单元，具体模型如下图。
施加偏心荷载，将荷载划分为轴向力和一个弯矩，加载图如下。
计算前将进行如下设置，按大位移计算，打开预应力开关，具体如图。
计算结果，位移图如下。
计算所得柱中心位移为1.54mm。
3.5结果分析
ansys模拟计算出来跨中最大挠度为：1.54mm 利用前面建立的理论公式计算跨中挠度：
max
2sin csc 12l y e l e e αα⎡⎤⎢⎥=-=-⎢⎥⎢⎥⎣⎦
10000P N =
206000000000E Pa =
3
12
bh I = b=0.1m h=0.1m t=0.01m
I =（0.1x 0.13-0.08x 0.083）/12= 4.92x 10-6 l =5m e=0.05m
所以：
=1.582mm
对比有限元软件结果1.54mm 与理论公式结果1.58mm ，差异度只有2%左右。