|
|
马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。
您需要 登录 才可以下载或查看，没有账号？我要加入
x
本帖最后由 Rainyboy 于 2015-6-5 10:42 编辑
============
系 列 前 言
============
mayaview同学关于波动和振动的关系的帖子振动和波动有什么区别和联系得到了大家的热烈讨论，其中有许多非常有意思的论点。我之前也在辨析“波动”和“振动”这两个概念上花了一些力气，通过一些计算对它们进行了一些理解。一直有将它们分享出来的愿望，大家在上面那个帖子里的讨论更给了我动力。
关于辨析和理解“波动”和“振动”的途径，我想解决困惑最直接的办法是进行一些计算，看看所得的结果是否与预期相同。
我想通过这个系列的帖子整理自己之前为了理解和辨析这两个概念所进行的计算，如果同志们自己也有一些有意思的算例，不妨贴上来，大家一起玩赏一下。
在开始之前，我们姑且约定：
将以模态叠加法和动力刚度法为基础的结构响应求解方法称为“振动”方法，它们用一系列频率上离散的“模态”来刻画系统的固有特性；
将以傅里叶变换法为基础的结构响应求解方法称为“波动”方法，它们用“波数空间”（wave number space）来刻画系统的固有特性。
这一系列算例其实是为那些一直以来沉浸于“振动”方法，却打算（或不得不）稍严肃地考察“波动”方法及其范畴的同志准备的。
我相信，这些算例中很有可能是存在错误的:(，所以请准备好你挑剔的眼光，我们一起在批判中进步！
Rainyboy
2014年2月14日
========================
算例（一）：在有限结构上算出“波动”
========================
作为开始的第一个算例，先让我们在最熟悉的“振动”方法领域中，用最熟悉的工具（有限元，ANSYS）算出波动的效果吧！
这个算例可以帮助我们利用手边的“振动”分析工具观察波动效应。
对应的APDL是：
- FINISH
- /CLEAR
- /PREP7
- /eshape,1
- BL = 1
- BW = 2e-2
- BH = 2e-2
- NMESH = 100
- ET,1,BEAM188
- MP,EX,1,2.1E11
- MP,PRXY,1,0.3
- MP,DENS,1,7800
- SECTYPE,1,BEAM,RECT
- SECDATA,BW,BH,2,2
- K,1,-BL,0,0
- K,3,0,0,0
- K,2,BL,0,0
- L,2,3,NMESH
- L,1,3,NMESH
- LMESH,ALL
- NSEL,S,LOC,X,-BL
- D,ALL,ALL
- ALLSEL,ALL
- NSEL,S,LOC,X,BL
- D,ALL,ALL
- ALLSEL,ALL
- !**********************
- !模态分析
- !**********************
- !/SOLU
- !ANTYPE,2!指定为模态分析
- !MODOPT,LANB,10 !指定提取特征值数及所用方法
- !MXPAND,10,,,YES!指定扩展的模态阶数
- !SOLVE
- /PREP7
- *SET,N_CIRCLE,20!计算的周期数
- *SET,NSTP_C,30!每周期的载荷步数
- *SET,FRE,10000!载荷频率
- *SET,CURR_T,0
- *SET,DT,1.0/FRE/NSTP_C
- *SET,PI,3.141592653
- /SOLU
- ANTYPE,4!指定分析类型：动力学分析
- TRNOPT,FULL!完全法
- TIMINT,ON!时间积分效应
- OUTRES,ALL,LAST !所有子步，所有记录写入DATABASE
- *DO,I,0,N_CIRCLE-1,1!外部循环：共计算N_CIRCLE个周期
- *DO,J,1,NSTP_C,1!内部循环：每个周期NSTP_C个载荷步
- CURR_T=(I*NSTP_C+J)*DT!计算当前载荷步的时间点
- *IF,I,EQ,0,THEN
- ALLSEL,ALL
- NSEL,S,LOC,X,0 !选择加载点
- F,ALL,FY,100*SIN(2*PI*FRE*CURR_T)!施加当前时间点的载荷
- ALLSEL,ALL
- *ENDIF
- ALLSEL,ALL
- TIME,CURR_T !设置当前载荷步的时间点
- AUTOTS,1 !自动时间步长
- KBC,0 !载荷步之间线性插值：递增方式
- DELTIM,DT/10,DT/12,DT/8 !本载荷步之中指定时间步长
- LSWRITE,I*NSTP_C+J !写载荷步
- *ENDDO
- *ENDDO
- LSSOLVE,1,N_CIRCLE*NSTP_C!求解所有载荷步
- FINISH
若将时序计算的位移结果绘成动态图（由于一些原因，并未从0秒开始作图），就是这样：
详细的计算参数和结果简要总结如下：
=======
Remarks
=======
这是一个在有限域上进行的瞬态计算，激振力只按照正弦形式持续了1个周期，后续的过程既可以理解为自由振动，也可以理解为波的传导。
除了结构的震动，在结果中我们发现了较为典型的“波动”现象，尤其在波动尚未传递到边界时，见动态图。
而通过最后绘制的两个离激振点不同距离的点，绘制了它们的响应时程，更是清晰地发现了“滞后”现象。
这次计算是两年前作的了，在计算之前与当时同一个实验室的同志打赌，让大家猜计算的结果大概长什么样子，猜错的请猜对的人吃饭。
结果有些同志说一开始结构的所有点就会参与响应，有些甚至断言所有的点会以相同的相位响应。
我虽然觉得应该是有一些部分先变形，另一些部分后变形，但是不确信这样的特性能否由ANSYS（有限元）这样的“振动”分析工具计算出来。
其实现在想想，有限元哪里是什么“振动”分析工具，有限元就是一个普遍的动力响应分析工具，很多基于“波动”理论的方法也是以有限元为起点的。
好了，第一个算例就说这么多吧，大家可以把APDL拿去，多算一些时间，换一个更低/高的频率，再计算看看。如果计算结果与你的预期完全一样吗？不一样的话，是哪里出了问题呢？
===========================
系列链接
===========================
“理解波动”系列算例（一）：在有限结构上算出“波动” （同时包含对这系列算例的一些说明）
“理解波动”系列算例(二)：将无限域波动解映射到有限域
理解波动系列算例（三）：Wave-mode duality，以悬臂梁为例.
|
评分
-
1
查看全部评分
-
|
发表于 2014-2-14 03:15
提升卡
变色卡
显身卡
发表于 2014-2-14 08:11
发表于 2014-2-16 22:27
楼主
谢谢楼主无私分享！！！
发表于 2016-4-25 16:56
