接触问题的模态分析实例
刚出炉的：）一块蒙皮在工作条件下的模态分析结果。
四个角固定在骨架上，考虑蒙皮跟骨架的接触、摩擦系数，蒙皮的自重，工作载荷
建模和后处理都是在ANSYS中完成的，求解器是自己的。
nlp000.png是模型的网格。
nlp003.png是自重作用下的静态接触分析，变形放得有点大。
最后是用响应谱的数据做的比较。
[[i] 本帖最后由 wsjxaut 于 2006-9-22 09:47 编辑 [/i]]
图片
图片图片
图片图片
图片图片
图片 牛人 哈哈，被发现了。 真厉害啊，羡慕死了，我什么时候能达到这样啊。一定好好学。很不错!
张见识了啊 好好学习 这个是上面那个问题的线性模态，不考虑接触，四个角固定，在ANSYS中计算的。前面几阶跟非线性计算的结果比较相似。[[i] 本帖最后由 wsjxaut 于 2006-9-18 14:41 编辑 [/i]] 0Hz的图片是接触结构的冲击响应谱，峰值对应的频率与我们的非线性求解器NLMA计算得到的模态频率吻合的非常好，说明NLMA的计算结果还是比较靠谱的 :)
[[i] 本帖最后由 wsjxaut 于 2006-9-20 23:52 编辑 [/i]] 还有两个相关的非线性模态分析的东西：
机床的非线性模态分析：[[url]http://www.simwe.com/forum/viewthread.php?tid=740835&highlight=%C4%A3%CC%AC%2Bwsjxaut[/url]
悬臂梁的接触模态分析：[url]http://www.simwe.com/forum/viewthread.php?tid=740495&highlight=%C4%A3%CC%AC%2Bwsjxaut[/url]
[[i] 本帖最后由 wsjxaut 于 2006-9-24 17:50 编辑 [/i]] 下面这个是MARC做模态接触的讨论：
[url]http://www.simwe.com/forum/thread-215057-1-107.html[/url] 关于NASTRAN的非线性模态分析：
[url]http://www.chinavib.com/forum/thread-24082-2-1.html[/url]
NASTRAN的非线性模态分析也是预应力的模态分析：先算非线性的静态模型，然后把非线性因素产生的附加刚度考虑进去，在不考虑非线性因素的情况下用通常的线性模态的提取技术获得模态。对于模态分析这一步，模型是线性的，大约是因为在一定程度上考虑了非线性的因素，所以称为非线性模态，真正的非线性分析其实只是前面静态分析这一步是非线性的（当然也可以是线性的）。这还不是真正意义上的非线性模态。其它的一些软件，如ALGOR和ADINA说的非线性模态也都跟NASTRAN相同，MARC可以这么做，但是没有明确说这个叫非线性模态，ANSYS把这种分析叫预应力的模态分析，也没叫非线性模态分析。这种方法正如ADINA帮助里说的，只是在非线性程度非常小的情况下才可以使用（感觉就是废话）。
下面是一个关于NASTRAN非线性能力的介绍，里面专门强调其非线性模态是指预应力结构，比较诚实，呵呵：
发信人: bergkamp (至尊玉), 信区: Civil
标 题: MSC.Nastran for Windows简介(4)
发信站: 飘渺水云间 (Tue Apr 9 23:45:07 2002), 转信
四、非线性分析（Nonlinear Analysis）
非线性模块是MSC.NASTRAN for Windows的一个可选择模块，能使你对结构的非线
性特性进行分析。 结构中的非线性特性是由于材料的非线性，大应变以及（由接
触而引起的）边界条件的变化。所有这些特性对结构的部件或整个系统都有巨大的
影响。非线性分析包括静态和瞬态的非线性分析。
什麽是非线性分析？
绝大多数的物理现象都有非线性效应。分析人员经常用线性关系简化地建立模型。
但有的时候，这种简化是不合适的。
非线性效应由于以下情况产生：
非线性材料
大变形
接触
MSC.NASTRAN for Windows也可以解决随时间变化（瞬态）的非线性问题。
在非线性分析中，MSC/NASTRAN用非线性单元或材料特性来表达这些特性。要得到
非线性解，经常需要一些特殊技术。MSC/NASTRAN for Windows的用户界面使得生
成和设定这些特性并进而获得非线性解更加容易。
材料非线性特性
塑性（小应变）
Von Miss和Tresca屈服准则
Mohr-Coulmb和Drucker-Prager屈服准则
各向同性，运动和组合硬化
双线性或多点应力-应变曲线
大应变的超弹性材料
小应变的非线弹性
热弹性
粘弹性（蠕变）
粘弹性与塑性的组合
几何非线性特性
大位移和大转动（更新的Lagrangian）
大应变（超弹性材料的总体Lagrangian）
跳跃（snap-through）分析
追随力
位移控制
非线性边界条件
结构-结构接触（滑移线）
缝隙的开合，有无摩擦
强迫变形
非线性分析再启动
非线性曲屈
非线性模态分析（预应力结构）
自动求解方法
可选择的BFGS法和线性搜索法
自适应二分法和载荷增量的恢复
载荷与位移的控制法
组合法（弧长法，Crisfield,Riks法和修正的Riks法）
自适应弧长法
路径反转的自适应校正
卸载的刚度更新
自适应间隙单元
隐式的Newmark-Beta时间步长积分
自动二分法和自动刚度更新的自适应时间步长
容易的线性分析至非线性分析的转换
仅在线性模型上填加一些非线性特性
输入与输出的格式与线性分析类似
有效的数值求解过程
稀疏矩阵求解器可提高计算速度，有效利用磁盘空间
自适应的刚度更新
--
\__/ \__/ \__/ \__/ \__/ \__/ \__/
(oo) (o-) (@@) (xx) (--) ( ) (OO)
//||\\ //||\\ //||\\ //||\\ //||\\ //||\\ //||\\
bug bug bug/w dead bug blind bug after
winking hangover bug sleeping bug seeing a
female bug
※ 来源:·飘渺水云间 Freecity.dhs.org·[FROM: bergkamp]
[ 本帖最后由 wsjxaut 于 2006-10-4 00:59 编辑 ]
唉，文也不行，武也不行，你不做山贼，你想做状元啊？
我有想过…… ADINA, ALGOR, NASTRAN等软件的非线性模态，以及ANSYS的预应力模态分析，其实都是一回事，傲雪的一个帖子说得比较明白：
〖讨论〗结构预应力改变结构刚性问题 zz
2006-10-07 22:26:19
大中小
--------------------------------------------------------------------------------
旋转叶片会使叶片刚性增加，从而是结构模态发生改变，请问大虾该现象怎么样由力学基础理论来分析，进而怎么样从ansys算得负载时的模态？是不是先做一次预分析，加上pstres,on命令，然后再作模态分析？
普通模态分析的动力学方程为M*a+K*x=0，在这个方程中第一项是惯性力，第二项是由于叶片振动时偏离平衡位置的应力产生的回复力（叶片内力）。这两个力是平衡的。
当叶片旋转时，在叶身内产生离心力，其方向与叶片的径向相同。而叶片振动一般是弯曲型的，其运动方向是横向的与径向垂直，两者正交。对于线性小变形理论，径向的力应该不会对横向变形有影响。
但这只是小变形理论的一个假设。在叶片振动时由于偏离平衡位置一段距离x，径向的离心力与此距离的乘积会形成一个力矩，其作用与叶片的内力相同，使结构回到平衡位置。力矩的大小也与变形x成正比，可以表示成Ks*x。这时的动力学方程表示为：M*a+（K+Ks）*x=0，其中的Ks称为初应力矩阵。
由此可见由于初应力矩阵的影响使结构的刚度增加，因此其振动频率也会增加。
还可以举一个更简单的例子来说明：当我们弹动琴弦时，会由于振动发出声音。如果将弦拉紧，则弦的振动频率就会升高，音调就会高一些。由于拉紧力“加速”了弦回复到平衡位置的速度。
ANSYS在分析这类问题时，先加离心力进行静力分析，但分析时要将预应力选项打开以形成初应力矩阵，然后再进行模态分析。其他的有限元软件也是这么做的。
有许多讲振动方面的教材都讲过这个问题。比如我上研究生时的教材，西交大倪振华的《振动力学》 ，在梁的振动一章中就有这方面的内容。
[url]http://www.aoxue.org/cgi-bin/bbs/topic_show.cgi?id=73838&h=1#571575[/url] 非线性的模态分析软件NLMA简介
传统的模态分析技术是基于线性模型的线性模态分析。通常情况下，机器是有若干个零件构成的，零件与零件之间通过各种方式联接起来，接触是零件联接的基本特征。在力学分析中，零件间的接触状态跟载荷有关，表现出强烈的非线性。机器工作过程中产生的振动也表现出非线性特征。传统的模态分析技术无法有效地处理这种机器的非线性振动的模态问题。为此，人们采取了很多近似的分析方法。首先，可以不考虑零件间的接触特性，将若干个零件视为一体，这个时候，问题变为线性，可以采用传统的模态分析技术获得机器整机的振动模态。其次，可以将零件间的联接简化成线性弹簧，弹簧的刚度由联接刚度决定。这样，问题也变成线性的，可以采用传统的模态分析技术获得机器的整机模态。第三，比较理想的情况是把零件之间的联接特性用非线性的弹簧来反映，但是这个时候没办法采用线性的模态分析技术获得机器整机的模态。此外，还有一种做模态分析的方法，就是在考虑接触特性的情况下做一次静态非线性分析，获得在静态载荷作用下的结构应力，然后把得到的应力以附加刚度的形式叠加到机器刚度上，最后在不考虑接触的条件下对这种具有附加刚度的机器结构做线性模态分析，获得机器的整机模态。该方法在一定程度上考虑了接触特性对模态的影响，但是对机器整机模态的处理仍然是不全面的，仍然是线性模态，并不能反映模态的非线性特征，如模态的时变特征和瞬态载荷对模态的影响。该方法是目前在商业软件中进行非线性模态分析采用的主流方法。
NLMA是基于全新的非线性模态分析技术编写的第一个真正意义上的非线性模态分析软件。NLMA的基础是非线性模态分析技术，能够充分考虑非线性模态的特征，甚至可以用作非线性模态特征的分析工具应用于非线性模态的分析研究中。NLMA目前利用ANSYS软件完成前后处理，即模型的建立是在ANSYS的前处理器中完成并将模型数据提交给NLMA，NLMA的计算结果则由ANSYS的后处理器读入并完成阵型的显示。
NLMA可以用于包含接触的机器整机动态性能分析中，完成机器整机的非线性模态分析，可以帮助用户分析机器的振动模态与动态载荷之间的关系，工作条件下机器的振动模态与时间的变化关系。由于NLMA是基于非线性的模态分析技术，因此模态分析的范围并不仅限于接触问题的模态分析，也可以完成其它类型的非线性模态分析，如大位移、塑性等类型的非线性问题的模态分析。在这方面一个典型的应用是多关节机械手在大位移运动过程中的非线性模态分析。此外，NLMA也可以用作一种研究非线性模态问题的数值分析工具，帮助了解非线性模态的特点和性质。
关键词：非线性模态 接触模态 有限元 NLMA ANSYS
KEYWORD: Nonlinear modal, Contact modal, Finite element, NLMA, ANSYS
[[i] 本帖最后由 wsjxaut 于 2006-10-12 12:45 编辑 [/i]] NLMA计算结果在blog上的链接：
[url]http://blog.sina.com.cn/u/4705a544010005ey[/url]
[url]http://blog.sina.com.cn/u/4705a544010006tq[/url]
[url]http://blog.sina.com.cn/u/4705a544010006tp[/url]
软件的介绍：
[url]http://blog.sina.com.cn/u/4705a544010006n5[/url] 怎么才能获得这个软件呢？我想试试nlma的计算功能。 这考虑接触有意义吗？
字数不够啊.... 应该这样说：对多个零件构成的机器，整机的模态分析中不考虑接触，算出来的模态有意义么？ 嗯，可以合作
原帖由 [i]bulingshan[/i] 于 2006-10-25 11:14 发表
怎么才能获得这个软件呢？我想试试nlma的计算功能。 [/quote] 嗯，访问量超过1000！
..................................................... 精华置顶，版主英明，呵呵
——————————————————————————————————————————
[[i] 本帖最后由 wsjxaut 于 2007-1-30 23:02 编辑 [/i]] 楼主，不是一般的厉害啊。中国有限元的明天，还需楼主这样的人潜心研究。。。 做得真好啊，什么时候也能这么厉害就好了 牛人**********佩服 这个软件是LZ自己编的啊?
太强了
wsjxaut,您的意思是说也可以进行材料为非线性的模态分析?比如橡胶
能让我试试您的软件吗?
[[i] 本帖最后由 tongji_lu 于 2007-4-15 21:47 编辑 [/i]]
回复 #28 tongji_lu 的帖子
NLMA是基于有限元模型的非线性模态分析，基本上有限元法能建立起来的的非线性问题都可以做的 ANSYS在作模态分析的时候一般会把接触等非线性因素忽略掉吧？感觉作接触分析同节点结合在效果上基本相差不大。
回复 #5 wsjxaut 的帖子
分析的对吗？ 边界都穿透了！！:@ :@你第一个图片结果就不收敛,后面的计算就更没法看了,
先计算收敛再说吧? 边界穿透更说不过去
[attach]118804[/attach]
[[i] 本帖最后由 ansyslove 于 2007-4-20 18:33 编辑 [/i]] 非线性模态的性质跟线性模态有很大不同，穿透的问题，主要是存在几个模态的耦合，振动的时候，几个振型是叠加在一起的，总体的效果是不会穿透的，线性的模态，一个很重要的性质就是不耦合，非线性的模态，通常不具有这个性质，给出的振型图只是表示有这样一种振型存在，以线性模态的概念去理解非线性模态是不对的。
除此以外，非线性模态还有一些很有意思的性质，有时间会慢慢发上来 王老师来点APDL享受一下哈,
呵呵,你的建模水平很高,为什么不做个专贴,
嘻嘻,现在学ANSYS的越来越多放弃的越很多,就是因为ANSYS建模操作不便
看了哪么多贴子,多数为单一体的模形计算分析,真是有点打击ANSYS的疼处!!! [quote]原帖由 [i]ansyslove[/i] 于 2007-4-18 13:08 发表 [url=http://www.simwe.com/forum/redirect.php?goto=findpost&pid=982643&ptid=737234][img]http://www.simwe.com/forum/images/common/back.gif[/img][/url]
分析的对吗？ 边界都穿透了！！:@ :@
你第一个图片结果就不收敛,后面的计算就更没法看了,
先计算收敛再说吧? 边界穿透跟说不过去
118804 [/quote]
那个是后处理的bug，用脑子想想，有可能算那么多步么？ANSYS里这类小毛病很多，再举一例，下面这个云图，有可能是计算出来的么？
[[i] 本帖最后由 wsjxaut 于 2007-4-19 21:04 编辑 [/i]]
回复 #34 wsjxaut 的帖子
不要把问题归结为bug, sub=999999是计算不收敛的标志 好吧，按照您的意愿，承认我们算了999，999步，具体花了多长时间记不得了，反正是最后得到了一个不收敛的结果，后面的结果当然也是看不成的了。[[i] 本帖最后由 wsjxaut 于 2007-4-20 20:38 编辑 [/i]]
回复 #36 wsjxaut 的帖子
你说的不对, ANSYS出来的sub=999999并不是算了999999步的意思, 只是不收敛标志而已. 那个模型算的时候挺好的，后来做后处理的时候，读进去的结果里，子步数变成了6个9，其它数据倒是还正常，所以认为那是后处理显示的问题。刚才把从前的模型找出来，在11.0上又跑了一下，最大Z向变形跟那个图的结果只有1.5%的误差，所以结果还是靠得住的。模型的具体参数：
弹性模量：2.0E+11
泊松比：0.3
密度：7800
板厚：0.001
骨架和蒙皮都是壳单元，骨架边宽2.5E-02，骨架和蒙皮外形一致，都是正方形，边长0.6
骨架全部固定，蒙皮在四个角固定，蒙皮和骨架间加接触单元
蒙皮表面加面载荷50
另外，应该注意到，那个图里除了子步数有问题以外，DMX，SMN，SMX显示的数值也很诡异
[[i] 本帖最后由 wsjxaut 于 2007-4-21 11:01 编辑 [/i]] 嗯，不错，谢谢分享哦！ 太强悍了！太羡慕拉
页:
[1]
2