约束处理方法

fhyd212

18# bbssbb
楼主能否详细解释一下过约束 及如何释放， 新手上路 还请指教哈

pasuka

个人观点：
静力问题，楼主说的3种方法都可用。

模态分析：
1  用消去法最好。
2  如果用置0赋1法，则应该有所变通：假设约束自由度为 i ，将刚度和质量矩阵的非对角元素置0，即
       k_{ij}=k_{ji}=0
   ...
refeihc 发表于 2010-5-19 00:34

消去法、置0赋1和置大数，在动力计算中同样可以使用，至少我编程计算中都用过，没出现什么不良反应。不要误导lz
btw，lz到底是自己编程在算东西，还是用商业软件或者别人写都软件在算东西？还时不时蹦出GUI操作？真有点看不明白。。。

studentxj

22# pasuka

同意楼上说法！
     这三种方法本质上是在消除总刚和质量阵的奇异性，都属于数值计算方法，在对刚度阵的处理理论上没有太大的差别，有也是数值精度上的差别，静力计算中三种方法自己编程都实验过，没有差别。
     至于以上讨论的动力问题，是不是质量阵处理也如此没有做过，待考虑。
     楼主讨论的问题，估计是用商用软件做的，这个需要具体看看软件自身的某些设置！

fhyd212

22# pasuka
我用的是有限元软件 comsol 不是编程哈

fhyd212

23# studentxj
是的，确实是要看软件是怎么设置的，而且理论上来说三种方法的结果应该是相同或相近的。现在的问题就是这是个必选项，默认的是消去法，在处理某些问题（比如三维问题）的时候需要使用另外两种方法，主要是置大数法，所以还请多多指教

refeihc

以前和别人讨论过“大质量法”和“大刚度法”，好久不关心，快忘了，不记得是不是与约束处理有关。

静力分析3种处理方法都可以，这没有问题。

对于动力分析问题，刚才又考虑了一下，觉得置大数法可以是置大质量，或是置大刚度，也可以2个同时置。

caoer

I often use the penalty method to set boundary condition, which is easy to be coded. However the main drawback is that BC should be imposed weakly if using penalty method.

fhyd212

27# caoer
这种情况我倒是没怎么处理，我没有使用弱约束，是不是这里出问题了呢。或者说三种约束处理条件与是不是弱约束有关呢，关系是什么呢

fhyd212

26# refeihc
对，理论上应该是一致的。不过是不是处理方法应该和约束处理有关，请具体讲解一下。

bbssbb

18# bbssbb  
楼主能否详细解释一下过约束 及如何释放， 新手上路 还请指教哈
fhyd212 发表于 2010-5-19 08:35

过约束问题非常常见，尤其是在处理与约束相关的问题，比如本质边界条件，不可压缩条件，摩擦接触条件等等。这个问题本质上是约束太多了，导致拉格朗日场（约束作用反力场，接触力场，压力场等）与位移场相比过于丰富，违背了所谓的LBB条件，从而出现了所谓的失稳现象，表现为拉格朗日场出现了与网格疏密无关的高频震荡。处理方法也很简单，
1。增加位移场的丰富程度，比如用高阶单元，bubble mode，hughes变分多尺度方法等等。
2。减少拉格朗日约束的数量，比如缩减积分之类的方法。
3。增加稳定项，用所谓的稳定有限元方法，比如压力投影算法之类的。

refeihc

26# refeihc  
对，理论上应该是一致的。不过是不是处理方法应该和约束处理有关，请具体讲解一下。
fhyd212 发表于 2010-5-19 19:06

印象中是用在抗震分析中，我没有钻研过，以前在Nastran的帮助文件中见过介绍，建议你查一下。

fhyd212

31# refeihc
好的 谢谢哈

fhyd212

30# bbssbb
多谢版主支持哈