探讨一下对于刚体对时间步的影响的

jasonxiao523 · 发表于 2015-11-7 20:42:54

本帖最后由 jasonxiao523 于 2015-11-8 09:58 编辑

查看d3hsp文件，part id10是刚体，刚体划分的网格不是很规整，可能有的单元是毫米级，土体单元是厘米级，两者密度在一个数量级。发现整个计算时间步是由刚体的1.9e-8决定的，但是最小时间步的单元又是土体单元，增大刚体的密度，发现整个计算时间能够节省，对于刚体节省计算时间的用法有点疑惑。同时对这几个time step也有点疑惑，改模型试了下图2和图3的time step和图1中的current minimum time step0.9倍数的关系，tssfac我设置的0.9。

roczhp · 发表于 2015-11-7 20:42:55

本帖最后由 roczhp 于 2015-11-8 21:49 编辑

jasonxiao523 发表于 2015-11-8 21:36
我变现体网格是一个立方体的土体每个网格都一样的，只要我有刚体的模型smallest里面的单元都不会是，但实 ...

和你的接触设置有关，你是type 14，*CONTACT_ERODING_SURFACE_TO_SURFACE，这个接触可以设置为接触时间步长控制整个计算时间步长，当然你可以设置为实际计算时间步长取决于单元尺寸的情况（此时刚体网格不影响实际计算时间步长）。在你现有的设置下，刚体网格确实会影响计算时间！！接触时间步长和接触对材料参数、网格大小等有关。

roczhp · 发表于 2015-11-8 16:04:42

本帖最后由 roczhp 于 2015-11-8 16:12 编辑

图1中surface timestep:接触稳定需要的时间步长；一般要保证实际时间计算步长（global timestep）不大于这个时间步长（contact surface timestep），而实际上global timestep稍大或稍小于 surface timestep都不会出现接触不稳定，但是每次计算程序都会有这个提醒。这个步长和接触对材料参数、接触对网格大小以及接触类型和设置有关。
图1中smallest timesteps：dyna按材料参数和网格大小计算的100个最小时间步长；在材料密度、模量一定时，和网格大小有关；这100个最小时间步长实际对应着100个最小尺寸单元。
图2中time step是计算一开始的初始时间步长。
图3中dt后的时间步长：实际计算时的时间步长（global timestep），只有这个时间步长影响你计算所需要的时间。

增大刚体密度，实际相当于selective mass scaling，mass scaling用来将动态计算转化为准静态计算，需要自己校验计算结果的可靠性。

一般情况下，刚体网格大小不会影响实际计算时间步长。最终的求解所需时间和很多因素有关。

jasonxiao523 · 发表于 2015-11-8 21:09:02

本帖最后由 jasonxiao523 于 2015-11-8 21:17 编辑

roczhp 发表于 2015-11-8 16:04
图1中surface timestep:接触稳定需要的时间步长；一般要保证实际时间计算步长（global timestep）不大 ...

谢谢你的解答

这几天我也一直在试，也不知道我探索的对不对，发现接触时间步长由于刚体一般网格比较细，弹性模量比较大大，所以刚体那个面步长一般比变形体小，在没用质量缩放的时候图三那个步长基本就是图一步长乘以0.9或者0.7 这个根据自己设置的来。使用质量缩放的时候，输出matsum发现刚体是不会出现质量增加的，所以只要增加不是太多，added mass那一项会是0，但是计算时间会成倍的减少很多，当然准确性我现在没验证，因为不使用质量缩放所需要的计算时间根本不能忍受！搜索下论坛发现也有两个帖子是一样的问题http://forum.simwe.com/forum.php ... 6%E9%97%B4%E6%AD%A5
http://forum.simwe.com/forum.php ... 3%E7%BD%91%E6%A0%BC
再弱弱问一句实际的时间步怎么远远小于smallest那100个单元的时间步。

roczhp · 发表于 2015-11-8 21:24:33

本帖最后由 roczhp 于 2015-11-8 21:29 编辑

jasonxiao523 发表于 2015-11-8 21:09
谢谢你的解答
这几天我也一直在试，也不知道我探索的对不对，发现接触时间步长由于刚体一般网格比较细， ...

"实际的时间步怎么远远小于smallest那100个单元的时间步"，一般不会，看仔细了，看100 smallest里面的最小值，这个最小值决定实际计算时间，除非你的单元划分是均匀的。链接那个帖子中只是提醒你注意时间步长不能太大，否则出现接触不稳定，就是无法实现有效的接触分析，如果接触没有问题，你可以不理会这个提醒。

jasonxiao523 · 发表于 2015-11-8 21:36:36

roczhp 发表于 2015-11-8 21:24
"实际的时间步怎么远远小于smallest那100个单元的时间步"，一般不会，看仔细了，看100 smallest里面的最小 ...

我变现体网格是一个立方体的土体每个网格都一样的，只要我有刚体的模型smallest里面的单元都不会是，但实际时间步是系统自己根据接触面的确定的时间步。只要刚体网格细就会很小，上午看了一个帖子提到刚体和变形体网格尽量匹配。

jasonxiao523 · 发表于 2015-11-8 21:56:40

roczhp 发表于 2015-11-8 21:47
和你的接触设置有关，你是type 14，*CONTACT_ERODING_SURFACE_TO_SURFACE，这个接触可以设置为接触时间步 ...

学习了，我再看看帮助手册，确实设置接触没太在意把card1几个号填了后面就没在意了

账号		自动登录	找回密码
密码			注册

账号		自动登录	找回密码
密码			立即注册

[ANSYS/LS-DYNA] 探讨一下对于刚体对时间步的影响的

最佳答案

评分