输入动力荷载后网格发生了畸变！（附命令代码）

zy7788250

边坡开挖后，施加动力荷载，以zvel的形式施加，模型底部是quiet边界，周围采用ff边界，进行动力计算时，一开始就出错，网格发生了畸变。在边坡顶部或者是底部输入都会出现这样的现象。
畸变的网格在边界上，会不会和动力边界条件有关系？请达人指教。

[ 本帖最后由 zy7788250 于 2007-10-23 20:30 编辑 ]

zy7788250

经过了修改以后，不再出现图中的那种畸变。
但是 计算时一开始还是出错，
math error:
invalid
并且 plo sur
出现 exaggerated grid distortion.
请达人指教。
谢谢。

chunhuo420

可以用小变形模式试试

zy7788250

非常感谢您的关注。
您说的小变形模式，是指不要写入set large命令吗？

zy7788250

没有用set large，还是出现了上面的情况。
请达人明示。

zy7788250

一直在调试
但是还是存在那样的问题
1.改变荷载量级，减小了一个数量级
2.改变荷载作用点，把荷载改成作用在边坡的侧面，边界条件也随之改成quiet边界
3.改变frequ，原先是选取的荷载的卓越频率，现在用坡体的自振频率
做了以上调试还是出现exaggerater grid distortion
因为是个开挖问题，发现开挖后台阶上的位移有0.35m，这个数值是不是偏大？
因此打算把坡体的静力边界由先前的
fix x y z rang z -72.1 -71.9
fix x range x -200.1 -199.9
fix x range x -0.1 0.1
fix y range y -90.1 -89.9
fix y range y 89.9 90.1
改为
fix x y z rang z -72.1 -71.9
fix x range x -200.1 -199.9
fix x range x -0.1 0.1
fix y
再试算一下
希望达人们能为小女出谋划策，不甚感激

zy7788250

结果还是出现上图中的问题

sch

不知道你的代码！看你介绍有可能节点不耦合！

zy7788250

网格重新调整过了
但是 动力计算还是进行不了
一到solve age 就出现 divided by zero或者是 invalid number
赋上命令流，请达人指教。
小女不胜感激。

taoyhrich

看你的命令流好象没有set dyna on 啊

zy7788250

加上去，还是不对
divided by zero

zy7788250

又重新进行了尝试
就是把模型直接建成开挖后的状态，（就是说，不存在开挖过程）
初始应力计算后，施加动力荷载，运行正常！

说明 开挖导致动力计算无法运行。
可不可以这样理解
m n range group ×× 后，动力计算认为group ××没有赋予材料属性，从而出现divided by zero或者是invalid number等错误提示？

请达人明示！

zy7788250

又将命令流进行了一些修改，前两步开挖后施加动荷载，动力计算可以进行，到了第四步开挖后，进行动力计算，出现了error，plo sur时可以看到exaggerated grid distoration，可不可以理解为，此时网格变形太大，是因为实际上是坡体发生了严重的破坏呢？

dubing222

apply nquiet squiet dquiet range z -72.1 -71.9
apply zvel 1.0 hist wave range x -97 -94 y -50 50 z 48

好像设置quiet边界就不能用速度输入地震波！

dubing222

不知道对不对

chensha

很多网格都没有连续起来，为什么要这样做呢？

sch

确实是网格节点不耦合。

chensha

这个网格估计静力计算都有问题，建议采用同一尺寸连续的网格，一样可以构造出想要得边坡形状。

zy7788250

我是在受动力荷载影响比较大的区域划分比较细的网格，粗网格和细网格用attach face 命令连接，为什么会不耦合呢？
另外，我看到书上说，可能是由于在边坡的坡面上出现了太多三角形网格，造成施加动力荷载发生了畸变。（见刘波那本书，P.227）

采用同一尺寸划分网格，就是要细化很多，这样计算必然会耗费大量的时间。

sch

你太固执了！无语。