如何对桁架结构进行稳定性分析
乌迪
积分 15
帖子 20
#1
2005-3-17 09:31
请教各位大侠，在sap 里面如何对桁架结构进行稳定性分析，具体如何操作啊？
yxs_li
积分 121
帖子 97
#2
2005-3-28 22:17
个人以为这是个很好的话题，不知为什么没有人感兴趣。我也是个SAP 初学者，谈一些粗浅的看法，不当之处望各位高手指教。
SAP2000并不能真正解决象类似桁架结构的整体稳定问题。对于局部构件的稳定问题则有点类似PKPM ，套规范公式求应力比解决，而不是在有限元的层次上解决。SAP2000虽有BUCKLING 分析，但仍不能解决整体稳定问题。BUCKLING 分析最多只能得到一个整体稳定的理论上限值（相当于分岔屈曲中的欧拉值），而不能考虑包含了初始缺陷及材料塑性在内的极值稳定问题。我现在也没搞明白SAP2000中的BUCKLING 分析具体的作用在哪里。我想是否顶多看一下前几阶模态是否正常，是否出现了局部稳定提前于整体稳定发生的情况以及看一个理论上的上限值（事实上这并没有意义）。在这一点上可以对比一下SAP2000和ANSYAS, ANSYS 中BUCKLING 的分析结果是要继续为下面的非线性分析提供初始缺陷用的，而SAP2000却到此为止了。因为初学，我不熟悉PUSH-OVER 这样的分析，不知道PUSH-OVER 是否可以解决整体稳定问题，估计也不行。
以上愚见，仅供参考。
ocean2000
积分 1207
帖子 941
#3
2005-3-31 18:54
这个问题主要分两类来讨论，bucking 分析相当于我们理解中的第一类稳定，这在实际应用中可以作为参考。真正的极值点失稳在sap 中可以考虑的，根据沈教授写的网壳稳定分析中的一句话：结构的稳定性可以从荷载-位移全过程曲线中得到完整的概念。那么我们也可以这么理解，只要sap 能做出这条曲线那么就可以解决问题，于是就利用到了sap 基于位移控制的非线性分析~！sap 分析参考中是这么叙述的：
当用户知道所期望的结构位移，但不知道施加多少荷载时，选择位移控制。这对于在分析过程中可能失去承载力而失稳的结构，是十分有用的。标准的应用包括静力 pushover 或 snap-though 屈曲分析。当使用位移控制时，用户必须选择一个位移分量来监控。这可以是节点的单个自由度，或一个用户以前定义的广义位移。
基于sap 本身的功能，可以做理想情况下的几何非线性的稳定分析，而且也可以得到相当准确的理想极值点失稳的曲线。但是问题：
一。是sap 极值点分析的初始条件（即初始缺陷的处理）不好弄，比如网架规范要求考虑第
一模态。但是sap 无法象ansys 那样将考虑的模态位移做为初始缺陷。
二。是无法准确的考虑材料非线性，目前常用的是采用非线性link 单元和塑性铰来模拟，当然在某些情况下，比如只考虑几何非线性的稳定分析，这也足够了。
yxs_li
积分 121
帖子 97
#4
2005-4-1 20:00
OCEAN2000说的精彩。
因为SAP2000不能给结构附加一个初始的缺陷，所以不能从本质上分析复杂空间桁架的整体稳定性。虽然可以跟踪荷载－变形曲线，但我们无从判断这条曲线得到的极值点是否是所有可能的分析中最小的。举个例子说，一个鸡蛋可以一把抓过来捏碎，也可以两个手指捏碎，也可以两个手挤碎，总之方法有多种，这个破碎力的大小因施加荷载的方式和位置的变化而变化。同样对于一个复杂结构而言，如果不知道加核的方式和位置，即使得到了一条荷载－变形曲线，得到的那个最小值或许比另一条曲线得到的最小值要大。但BUCKLING 分析也部分提供了加载的参考位置。所以基本上SAP2000还是以构件的局部稳定来确保结构的整体稳定，这好像是个缺陷。而ANSYS 很好地解决了这个问题。
ocean2000
积分 1207
帖子 941
#5
2005-4-1 20:54
yxs_li 兄的说法也不尽然，sap 里可以选择初始状态，所以你可以用某一次非线性分析所得到的各类结果，作为我们所要分析工况的的初始状态，这可以用来模拟大部分的初始缺陷。只是不能利用模态分析的结果来考虑初始缺陷而已。
qianhl
积分 26
帖子 24
#6
2005-4-15 17:58
（1）整体结构的初始缺陷与结构的形式有很大关系，绝大不部分结构不是缺陷敏感结构，一般空间桁架都不是缺陷敏感结构，所以在分析整体稳定的时候可以不考虑初始缺陷的影响。最典型的缺陷敏感结构为单层球壳，一般双层结构对缺陷都不敏感。
（2）平时我们所说的初始缺陷均指几何初始缺陷，还有材料初始缺陷时程序、软件很难考虑的，所以即使我们考虑初始缺陷，也只是考虑了部分缺陷。这个设计者心理应该明白。我的意思只是我们没有必要刻意去寻求初始缺陷的完美解决，应该从整体上把握结构的形式。
（3）至于加什么模式的初始缺陷也是一个值的探讨的问题，原则上是寻求一种最不利的缺陷模式，但是不同的结构有不同的最不利的初始缺陷模态。这就需要设计者的经验、试算比较。
（4） 如果说我们一定要用SAP 来解决初始的问题，也是可以的，以文本的格式输出节点的初始坐标，和模态的振型（节点的相对变形坐标），再用常用的程序（fortan,c)
把这些节
点坐标和变形坐标读出来，经过简单的处理求得带初始缺陷的模型节点坐标，再输入到SAP 即可。这可能有点难度，但是还是可以实现的，这就需要熟悉SAP 文本文件的格式， 风中的沙粒
积分 237
帖子 182
#7
2005-5-5 20:39
你们试过它的荷载位移曲线么? 个人感觉太滥了.
它的输入系统不好, 对于比较专业(稳定) 的问题处理得不好.
后处理, 规范演算更是糟糕, 真不知道他是怎么被称之为世界一流的结构分心软件的 peter1stli
积分 117
帖子 87
#8
2005-5-27 13:29
ocean2000提到采用基于位移控制的非线性分析可以得到结构的荷载－位移曲线
根据我的理解，如果选择的位移控制节点如果在结构屈曲之前如果发生位移回飘（位移值减小再增大），恐怕单纯靠位移控制是无法跟踪到它的全过程曲线的，因为非线性有限元的全过程跟踪多采用的是“弧长法”系列，单纯的荷载或是位移控制都是办不到的。
ocean2000
积分 1207
帖子 941
#9
2005-5-27 15:19
不好意思，不知道是不是我数值计算没学好，还是你没学好。我认为这与什么非线性有限元没有什么联系。位移控制或者荷载控制，其中跟踪和控制收敛大多可以采取“线性法”，“切线法”，“弧长法”等，不同的方法主要区别在于如何控制“收敛”的速度和精度上。只不过目前“弧长法”相对更易于收敛，但是也费资源而已。理解有误？请概念熟悉的朋友指正！ 冲令狐
积分 1626
帖子 1391
#10
2005-5-28 09:07
我个人认为ocean2000兄关于控制收敛方法的说法是正确的。
初九，潜龙，勿用。我是义工甲@okok
cx
积分 124
帖子 89
2005-6-7 15:39
采用" 弧长法”可以跟踪结构的屈曲后强度, 对荷载-位移全过程分析是很有效的.
zyldemail
积分 130
帖子 124
#12
2005-6-14 14:21
大家对整体稳定性讨论的够激烈够深度，可是我看了之后还是一头雾水，大家只是理论上争来争去，就没人具体做了实际的工程跟大家参考一下，有了事实根据，争论自然结束，本人正好在做一个课题，其中也需要整体稳定性的计算，正是整体稳定性一直困扰着我进展很慢，如果有算例上传，大家一看就明白了，大家凭理论都有一套，我觉得应该拿出实际算例，才能服众，也能供功底不深的哥们参考学习，本人也是一个sap 的爱好者和使用者，真希望大家能说得多，做的更多，这样论坛才更有价值，才更有吸引力。真希望众多高手们无私上传整体稳定性的算例，本人将焦虑恭候，不胜感激。
care (care)
积分 364
帖子 147
#13
2005-7-24 23:53
1 yxs_li 讲得push-over 是抗震性能评价的一种分析，目前对于对称的并且
以第一振型（参与系数大于90%）为主的结构效果佳，可以下载周锡元院士和巍巍曾经发表的论文。
2 进行荷载――位移全曲线分析 （1） 保证模型符合事实要求（位移约束处的节点刚度对结果有影响） （2） 加荷载f （任意大小，进行特征值屈曲分析（得到屈曲荷载系数a,f 大a 小， f小a 大，但乘积相同）
（3）将荷载系数乘以f 得到屈曲荷载F ，再进行非线性（几何刚度矩阵作修正）屈曲分析得到荷载――位移全曲线。
wyhaaa
积分 93
帖子 88
#14
2006-4-2 15:32
好像党八股，虚张声势，借以吓人，怎么没人说点能解决实际问题的！
若愚
积分 471
帖子 177
2006-7-5 11:49
有限元分析稳定的理论依据，我没搞懂，只知道传统的方法，例如能量法，
看了几本书，也没有看懂有限元分析稳定的思路，前段时间，sap 工作人员演示了一下，分析一个节点的稳定问题，我看分析结果，没有区分整体稳定和局部稳定，反正给出了几个数，也不晓得是如何搞出来的。
请教大家一下，有限元分析稳定的理论依据在哪里可以找得到，才疏学浅阿。
crazy
积分 175
帖子 131
#16
2006-8-22 14:58
看了大家的发言，受益非浅。
试着说说自己的看法，请大家批评指正。
问题1. 桁架模型和空间刚架模型在sap 或者ansys 中能不能算出局部稳定性呢？
我的看法是：桁架模型中，基本单元为二力杆(节点铰接) ，刚架模型中基本单元为三维梁单元(节点刚接) 。转化为有限元求解方程中后，其形式都是一个单元刚度矩阵。单元的细节，如截面宽度与厚度之比，是体现不出来的，而是以截面惯性矩、面积、最远点至中性轴的距离等形式表现，两个截面形式不同而参数一样的截面在刚度矩阵中没有区别。
所以以杆单元、梁单元为基本单元的模型中，是验算不了局部稳定的，位移-荷载曲线中也体现不了单元细节，而需要对每个单元按照规范进行局部稳定验算。
问题2：对结构进行逐步加载分析的时候，怎么能同时反应出结构局部失稳和整体失稳呢？ 我的看法是：基于问题1杆梁单元反映不出构件构造细节的特点，只有对整体结构的每个构件都采用板/壳/体单元建模，这样构件细节和整体构件布置都得到了反映，在更低层上分析结构。进行第一类稳定分析也是同样的道理。
所以由问题1和问题2可以得到推论：对一个用梁单元划分的轴向压杆进行屈曲分析得到的是整体稳定临界荷载值，如果将此压杆建立实际截面并用板/壳/体单元划分，首先得到的是整体稳定和局部稳定两者之间较小的临界值。
问题3：进行稳定分析的目的是什么呢？
我的看法是：第一个可能的目的是判断结构会不会在进入塑性前失稳，即考察结构是否有布置上的隐含的缺陷。因为构件的应力一般大家都控制的比较严格。此时第一类屈曲临界荷载值可以作为设计参考，第二类屈曲分析以荷载控制，选择荷载设计值为最大加载，可以适当放大。
第二个可能的目的是考察结构的性能，即结构在罕遇情况下进入塑性阶段后的工作性能，因为第一类屈曲分析是求解结构的特征值，刚度矩阵不变，所以此时第一类屈曲分析对设计没有参考意义。需要进行第二类屈曲分析，涉及到大变形、p-delta 效应、材料非线性。
问题4：sap2000能否进行第二类稳定分析？
我的看法是：前面ocean2000提到沈院士在研究网壳屈曲时说的话，我觉得很正确。这等同于问:sap2000能否产生微小初始位移，然后进行非线形(包括几何和材料的非线性) 分析，并跟踪极值屈曲后阶段。
首先关于初始微小位移，根据个人所知，ansys 中认为模态分析中第一阶模态是最易发生的屈曲形态(因为临界荷载最小) ，所以将第一阶模态乘以一个微小的因子作为各个节点的初始位移，来进行极值屈曲分析。个人认为在结构上任意一点（除了那些在加载中必定保持不动的点，如支座、对称中心等）施加一个微小的初始位移，都将导致极值屈曲分析的开始，因为结构是一个整体，微小的一点位移将传播开来，在整体内重分布，且由于其值微小，对结构分析结果不会产生显著影响。所以可得推论：以微小的荷载进行非线性分析，然后以此分析作为新分析的开始进行加载。
其次关于sap2000的材料非线性分析功能。与ansys 相比，ansys 中，构件进入塑性后，是按照预先定义的非线性材料的应力－应变滞回曲线来调整单元端部的内力的。而sap2000则做了简化，在杆件进入塑性后，在进入塑性的地方用塑性铰来代替塑性化那一时刻后的节点性能。
至于用位移控制还是荷载控制，考虑到在屈曲后，一些点的位移可能会出现后退，而结构也可能出现卸载，不能单纯地说一定取哪种方法。选用位移控制时，最好选那些在加载过程中位移一定会不断增加的点，至于那些是" 一定增长的" ，要靠我们判断了。
一些学习和工作中的认识，不对之处请多多指正，谢谢。
yxs_li
积分 121
帖子 97
#17
2006-9-6 09:51
看了上面的帖子，也想说几句
1、我想对一个网壳或空间桁架这样的整体结构而言，稳定会涉及三类问题：A. 整个结构的稳定性 B.构成结构的单个杆件的稳定性 C. 单个杆件里的局部稳定（如其中的板件的稳定） 对于A 类我想前面的帖子已说了很多，主要通过BUCKLING 分析，及非线性分析可得（当然SAP 对于材料非线性只能通过定义铰得方式实现）
对于B 类，以前我以为SAP 难以做到，但是最近做了一个算例，证明SAP 能够判断是结构整体先屈曲还是局部杆件先屈曲，关键在于杆件单元的划分，一般情况下我们总是以两个节点之间的杆件作为一个单元来进行分析，所以其BUCKLING 分析的值只能代表整个结构（因为杆件中间没有节点及其自由度），如果把单个杆件DIVIDE 成两个以上的单元(当然这样会增加计算量），就可以在BUCKLING 变形曲线中看到是整体先屈曲还是部分杆件先屈曲。当然这样也难以做到类似H 型钢梁的平面外稳定问题（因为SAP 中的梁单元不能考虑扇性惯性矩方向的自由度), 但对于由圆管构成的结构就不存在这个问题。这里又随之产生了另一个问题，因为我们知道BUCKLING 分析尽得到一个理论的上限，具体还是要通过非线性分析来得到承载力，如要对每根杆件的屈曲形态都进行一次存在初始变形的非线性分析显然不太现实，所以，本人认为在这点上还是遵循国家规范的原则来计算杆件的稳定（尽管这里面存在一个计算长度的选取问题而显得不完善，但总是安全的）
对于C 类问题，我想已不能在整体结构的范围内解决了，SAP 好像也没有好的办法，ANSYS 可以进行模拟一下。如果一般的设计工作，我想还是遵循规范，以高厚比的方式来确定比较方便。如按上贴中所说，按板、壳、体进行实体建模的话，显然是不现实的。一则建模量巨大，二则计算量巨大。SAP 能否实现也是一个问题。
2、对于一般的桁架结构，大多为对缺陷不敏感的结构，其BUCKLING 的一阶模态变形可以方便地通过荷载变形模拟
（如果方向与实际荷载方向相同则不需要施加这种模拟予变形的荷
载），所以不一定需要按一阶变形人为给模型施加一个初始缺陷（这在SAP 上很难实现）。 个人观点，请指教
WLJ7320
积分 129
帖子 115
#18
2006-10-6 22:46
BUCKLING 除了给结构引入初始缺陷外就没有其他作为了吗？我有一个想法，比如在框架（支撑）结构中我们可以通过判定结构的第一阶屈曲模态是无侧移还是有侧移的从而确定柱子的计算长度系数是选择有侧移的还是无侧移的，有时候即使是纯框架（如定层斜梁的2层框架，见童根树《钢结构的平面内稳定》）也有可能第一阶屈曲模态是无侧移。框架柱的计算长度除了跟梁柱的刚度比有关外还荷载的分布情况有关，合理的确定柱子的计算长度其实是一个比较繁琐的问题，但是其对一阶分析却很重要。通过查看第一阶屈曲模态我们可以窥见模型中柱子计算长度取值的合理性。
jiahua
积分 33
帖子 47
#19
2006-10-28 09:29
初始缺陷还可以提供一个临界荷载的参考值。
对于sap 里面的初始位移问题，既然不好施加初始缺陷，我觉得是不是可以考虑采用施加一个初始荷载后得到初始扰动作为替代。
philip
积分 267
帖子 227
#20
2006-11-18 00:29
[quote]ocean2000 wrote:
不好意思，不知道是不是我数值计算没学好，还是你没学好。我认为这与什么非线性有限元没有什么联系。位移控制或者荷载控制，其中跟踪和控制收敛大多可以采取“线性法”，“切线法”，“弧长法”等，不同的方法主要区别在于如何控制“收敛”的速度和精度上。只不过目前“弧长法”相对更易于收敛，但是也费资源而已。理解有螅壳敫拍钍煜さ呐笥阎刚
colt
积分 717
帖子 556
#21
2007-3-6 20:28
yxs_li wrote:
看了上面的帖子，也想说几句
1、我想对一个网壳或空间桁架这样的整体结构而言，稳定会涉及三类问题：A. 整个结构的稳定性 B.构成结构的单个杆件的稳定性 C. 单个杆件里的局部稳定（如其中的板件的稳定）
对于C 类问题，我想已不能在整体结构的范围内解决了，SAP 好像也没有好的办法，ANSYS 可以进行模拟一下。如果一般的设计工作，我想还是遵循规范，以高厚比的方式来确定比较方便。如按上贴中所说，按板、壳、体进行实体建模的话，显然是不现实的。一则建模量巨大，二则计算量巨大。SAP 能否实现也是一个问题。
这个问题很有工程实际意义，楼台分析的很全面。值得一提的是，基于目前SAP 的功能可以分析楼台所提的C 类问题。实体建模一直是SAP2000的一个瓶颈，对于桁架来讲杆系模型基本可以满足设计精度。
对于局部节点，如锻造节点，复杂节点，就需要实体分析。SAP2000通常处理的不错。SAP84可以基于AUTOmesh3D 模块，分析。
另一个问题就是ＳＡＰ２０００能否给出位移－荷载全过程曲线，一直我关注的问题。不知道诸位如何解决。
谢谢。
Magic Structure Design Consultant
care (care)
积分 364
帖子 147
#22
2007-3-6 21:46
colt 编辑提到的荷载位移全曲线方程，可以这样得到：（1）先进行 bucking 分析，得到试探荷载的屈曲系数，计算出屈曲荷载Pcr 。（2）然后再进行static 分析，这时分析类型选择nonlinear 而且results saved 选择multiple states。（3）分析完后，在display---show plot function 横坐标选择step, 纵坐标选择关心的反应就可以了。但是，稳定问题从荷载作用方式来讲为动力稳定和静力稳定，上面说的只是静力分析。sap2000中分析类型，就这几种（静力，多步静力，模态分析，时程分析，反应谱分析，功率谱分析，屈曲分析），这几种分析类型在原理上讲的都不是动力稳定性问题。动力稳定性问题能否用sap2000分析，请诸位同行指点一下。
知足常乐，知不足进取！
colt
积分 717
帖子 556
#23
2007-3-7 12:19
CARE 的处理手法，很巧妙。有一点让我困惑的是，横坐标中的ＳＴＥＰ 荷载是逐步增加的，怎么体现卸载过程呢？也就是下降段，如何处理。
Magic Structure Design Consultant
care (care)
积分 364
帖子 147
#24
2007-3-8 16:22
colt提的问题，很有意义。我是这样实现的：（1） 加载路径如楼上图所示，首先将何载分为n 段，每段加载值（有正有负，因为荷载有上升和下将段）记录下来，以此来描述加载路径；（ 2 ） 定义n 个Load case,利用static--nonlinear construction simulation,结构设定为一组，n 个Load case设定为n 个组，这时需要将steps-save 中选为end of each step，注意，2-n 个load 加载时的operation 选择为load all in group 。（3 ）运行分析后，在
display---show plot function 横坐标选择step, 纵坐标选择关心的反应。望同行指点。 自己试算了一个简单模型，见附件。
692544-lghtry.SDB ( 135.7 K) 下载次数 9
知足常乐，知不足进取！
dinochen1983
积分 1199
帖子 884
#25
2007-8-25 15:29
colt提的问题，很有意义。我是这样实现的：（1） 加载路径如楼上图所示，首先将何载分为n 段，每段加载值（有正有负，因为荷载有上升和下将段）记录下来，以此来描述加载路径；（ 2 ） 定义n 个Load case,利用static--nonlinear construction simulation,结构设定为一组，n 个Load case设定为n 个组，这时需要将steps-save 中选为end of each step，注意，2-n 个load 加载时的operation 选择为load all in group 。（3 ）运行分析后，在display---show plot function 横坐标选择step, 纵坐标选择关心的反应。望同行指点。 自己试算了一个简单模型，见附件。
对于屈曲或者破坏后的响应, 由于荷载上不去了, 所以不应该是用力去控制加载, 而是用位移, 位移加载可以看到下降段, 用施工荷载不断的把力荷载加上去是不能得到屈服点的 --------------------------------------------------------------
上面模型中的横坐标是step, 当时没有注意这个问题； sap2000只能处理弹性屈曲问题，所以个人认为利用sap2000还不能得到屈曲或者破坏后的响应。还望指点！
屈曲分析是一种用于确定结构开始变得不稳定时的临介荷载和屈曲结构发生屈曲响应时的模态形状的技术。ANSYS 提供两种结构屈曲荷载和屈曲模态分析方法：非线性屈曲分析和特征值屈曲分析。
非线性屈曲分析是在大变形效应开关打开的情况下的一种非线性静力学分析，该分析过程一直进行到结构的极限荷载或最大荷载。非线性屈曲分析的方法是，逐步地施加一个恒定的荷载增量，直到解开始发散为止。尤其重要的是，要一个足够小的荷载增量，来使荷载达到预期的临界屈曲荷载。若荷载增量太大，则屈曲分析所得到的屈曲荷载就可能不准确，在这种情况下打开自动时间步长功能，有助于避免这类问题，打开自动时间步长功能，ANSYS 程序将自动寻找屈曲荷载。
特征值屈曲分析步骤为：
1. 建模
2. 获得静力解 ：与一般静力学分析过程一致，但必须激活预应力影响，通常只施加一个单位荷载就行了
3. 获得特征屈曲解:
A. 进入求解
B. 定义分析类型
C. 定义分析选项
D. 定义荷载步选项
E. 求解
4. 扩展解
之后就可以察看结果了
示例1：
! ansys 7.0 有限元分析实用教程
! 3．命令流求解
! ANSYS命令流：
! Eigenvalue Buckling
FINISH ! 这两行命令清除当前数据
/CLEAR
/TITLE,Eigenvalue Buckling Analysis
/PREP7 ! 进入前处理器
ET,1,BEAM3 ! 选择单元
R,1,100,833.333,10 ! 定义实常数
MP,EX,1,200000 ! 弹性模量
MP,PRXY,1,0.3 ! 泊松比
K,1,0,0 !创建梁实体模型
K,2,0,100
L,1,2 !创建直线
ESIZE,10 ! 单元边长为1mm
LMESH,ALL,ALL ! 划分网格
FINISH ! 退出前处理
！屈曲特征值部分
/SOLU !进入求解
ANTYPE,STATIC ! 在进行屈服分析之前，ANSYS 需要从静态分析提取数据
PSTRES,ON ! 屈服分析中采用预应力
DK,1,ALL ! 定义约束
FK,2,FY,-1 ! 顶部施加载荷
SOLVE ! 求解
FINISH ! 退出求解
/SOLU ! 重新进入求解模型进行屈服分析
ANTYPE,BUCKLE ! 屈服分析类型
BUCOPT,LANB,1 ! 1阶模态，子空间法
SOLVE ! 求解
FINISH ! 退出求解
/SOLU ! 重新进入求解展开模态
EXPASS,ON ! 模态展开打开
MXPAND,1 ! 定义需要展开的阶数
SOLVE ! 求解
FINISH ! 退出求解
/POST1 ! 进入通用后处理
SET,LIST ! 列出特征值求解结果
SET,LAST ! 入感兴趣阶数模态结果
PLDISP ! 显示变形后图形
！NonLinear Buckling !非线性分析部分
FINISH !这两行命令清除当前数据
/CLEAR
/TITLE, Nonlinear Buckling Analysis
/PREP7 ! 进入前处理
ET,1,BEAM3 ! 选择单元
MP,EX,1,200000 ! 弹性模量
MP,PRXY,1,0.3 ! 泊松比
R,1,100,833.333,10 ! 定义实常数
K,1,0,0,0 ! 底端节点
K,2,0,100,0 ! 顶点
L,1,2 ! 连成线
ESIZE,1 ! 网格尺寸参数设定
LMESH,ALL ! 划分网格
FINISH ! 退出前处理
/SOLU
ANTYPE,STATIC
NLGEOM,ON
OUTRES,ALL,ALL
NSUBST,20
NEQIT,1000
AUTOTS,ON
LNSRCH,ON
/ESHAPE,1
DK,1,ALL,0
FK,2,FY,-50000
FK,2,FX,-250
SOLVE
FINISH
/POST26
RFORCE,2,1,F,Y
NSOL,3,2,U,Y
XV AR,2
PLV AR,3
! 进入求解 ! 静态分析类型（非屈服分析） ! 打开非线性大变形设置 ! 选择输出数据 ! 5个子步加载 ! 20步迭代 ! 自动时间步长 ! 激活线搜索选项 ! 显示二维状态下变形图 ! 约束底部节点 ! 顶部载荷稍微比特征值分析结果大! 施加水平扰动载荷 ! 求解 ! 退出求解 ! 进入时间-历程后处理器 ! 2# 变量表示力 ! 3# 变量表示 y 方向位移 ! 将 x 轴显示 2＃ 变量 ! y轴显示3＃变量数据
/AXLAB,Y,DEFLECTION ! 修改 y 轴标签
/AXLAB,X,LOAD ! 修改 x 轴标签
/REPLOT ! 重新显示图形
示例2：
! 悬臂梁受端部轴向压力作用的屈曲分析
! 先进行静力分析，在进行特征值屈曲分析，最后进行非线性分析
! 静力分析
fini
/cle
/filname,beam-flexure
/tittle,beam-flexure
/prep7
*set,f1,-1e6 ! 设置轴向压力荷载参数 et,1,beam189
mp,dens,1,7.85e3 ! 设置材料参数
mp,ex,1,2.06e11
mp,nuxy,1,0.2
sectype,1,beam,I,,2 ! 设置截面参数
secoffset,cent
secdata,0.15,0.15,0.25,0.015,0.015,0.015,0,0,0,0
k,1,0
k,2,2.5,0
k,3,1.25,1
lstr,1,2
latt,1,,1,,3,,1
lesize,1,,,10
lmesh,1
/view,1,1,1,1
/eshape,1.0
dk,1,,,,0,all,
fk,2,fx,f1 ! 施加关键点压力
/solu
antype,0
eqslv,spar ! 求解器设置 稀疏矩阵直接法 pstres,on ! 打开预应力开关 solve
finish
! 特征值屈曲分析
/solu
antype,1
bucopt,lanb,6,0 ! 取前六阶模态分析 mxpand,6,0,0,1,0.001
solve
finish
/post1
*do,i,1,6
set,i
pldisp,i
*enddo
*get,freq1,mode,1,freq
finish
! 非线性屈曲分析
/config,nres,200 ! 只记录两百步的结果 /prep7
tb,biso,1,1,2 ! 定义材料非线性 tbtemp,0
tbdata,,2.0e8,0
upgeom,0.01,1,1,'beam-flexure','rst'
! 对有限元模型进行一阶模态的位移结果 0.01 倍的修改
save,beam-flexure,db
finish
resu,beam-flexure,db
/solu
nlgeom,1 ! 打开大变形
outres,all,all
arclen,1,0 ! 弧长法设置
arctrm,l ! 弧长法终止准则 达到第一个峰值时终止计算 nsubst,200,,,1
fk,2,fx,f1*freq1
！fk,2,fx,f1*freq1*1.2 ! 将轴向压力值放大，放大系数为第一阶模态的主频 solve
本文来自：E-Works ，转载请保留原文链接http://blog.e-works.net.cn/182745/articles/218677.html ：