|
补充一下关于体积粘性阻尼的资料，来自于帮助文档，9.5.1 体粘性 体粘性引入了与体积应变相关的阻尼。它的目的是改进对高速动力学事件的模拟。ABAQUS/Explicit包括体粘性的线性和二次的形式。你可以在定义分析步时修改默认的体粘性参数，尽管只有很少情况需要这么做。因为仅仅将它作为一个数值影响，在材料点的应力中不包括体粘性压力。这样，并不将它考虑成为部分的材料本构响应。 线性体粘性 在默认情况下，总是包括线性体粘性来阻尼在单元最高阶频率中的振荡。根据下面的方程，它生成一个与体积应变率成线性关系的体粘性压力：  其中 是一个阻尼系数，它的默认值为0.06， 是当前的材料密度， 是当前的膨胀波速， 是单元的特征长度，和 是体积应变速率。 二次体粘性 仅在实体单元中（除了平面应力单元CPS4R外）包括二次体粘性，并且只有当体积应变速率可压缩时才采用。根据下面的方程，体粘性压力是应变速率的二次方：  其中， 是阻尼系数，它的默认值为1.2。 二次体粘性抹平了一个仅横跨几个单元的振荡波前，引入它是为了防止单元在极端高速度梯度下发生破坏。设想一个单元的简单问题，固定单元一个侧面的节点，并且另一个侧面的节点有一个指向固定节点方向的初始速度，如图9-12所示。稳定时间增量尺度是精确地等于一个膨胀波穿过单元的瞬时时间。因此，如果节点的初始速度是等于材料的膨胀波速，在一个时间增量里，这个单元发生崩溃至体积为零。二次体粘性压力引入一个阻抗压力以防止单元压溃。 
基于体粘性的临界阻尼比 体粘性压力只是基于每个单元的膨胀模式。在最高阶单元模式中的临界阻尼比由如下方程给出：  其中， 是临界阻尼比。线性项单独代表了6％的临界阻尼，而二次项一般是个更小的量。 原载于http://www.020fea.com/a/5/130/11307.html 自己也进行了想关的设置，对于线性体粘性，我把默认的0.06改为0.6，结果结果很快进入稳态，响应时间和幅度减小。 | 
IP卡
狗仔卡
发表于 2011-4-4 21:34:38
显身卡
发表于 2013-5-5 12:51:15