24小时热门版块排行榜    

查看: 5740  |  回复: 14

Ampare1987

新虫 (小有名气)

[求助] ANSYS约束不充分的问题 已有3人参与

我想让这个圆柱体转动起来,在圆柱体轴线节点上施加转Z向转动荷载,约束轴的节点的其他自由度,但是出现错误:
There is at least 1 small equation solver pivot term (e.g., at the UY   
degree of freedom of node 328).  Please check for an insufficiently     
constrained model.   
请问为什么约束不足呢,对于这样的实体单元模型应该施加多少个约束呢?
我的命令流如下:
/com   **** preprocessor ****
/PREP7
!******Element Type******
ET,1,SOLID185
ET,3,SHELL181
!******Material Props******
MP,DENS,1,7800  
MP,EX,1,2.06e11
MP,NUXY,1,0.3
!************Modeling************
r0=0.009                   !radius of the rigid is 9mm
h0=0.006                   !heigth of the rigid is 6mm
cyl4,,i,r0,,,,h0
HPTCREATE,AREA,1,0,COORD,0,0,0,
!************meshing************
!Area
asel,s,loc,z,0
AATT,1,,3,0
Esize,0.002
Amesh,all
!volume
VATT,       1, ,   1,       0
Esize,0.005
EXTOPT,ACLEAR,1
vsweep,all
!rotate the model
Nsel,s,loc,x,0
Nsel,r,loc,y,0
nplot
D,ALL, ,0, , , ,UX, , , , ,
D,ALL, ,0, , , ,UY, , , , ,
D,ALL, ,0, , , ,UZ, , , , ,
D,ALL,ROTX,0
D,ALL,ROTY,0
D,ALL,ROTZ,3.1415926/4
allsel,all
save
/com   **** Solution ****
/solu
solve
回复此楼

» 猜你喜欢

结构工程,泥石流拦挡结构
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

liyaobo0509

专家顾问 (著名写手)

实践检验过很多次,当然靠谱~可以试下MPC

发自小木虫Android客户端
7楼2018-06-06 19:48:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

kingspin

木虫 (著名写手)

【答案】应助回帖

感谢参与,应助指数 +1
隐式稳态计算应该是不可能让一个结构自由转动的。对于隐式稳态解,结构的三个位移和三个转动必须要确定且唯一才可能计算。
欢迎加入Digimat技术交流讨论群366061054,了解复合材料多尺度仿真技术
2楼2018-06-04 18:10:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Ampare1987

新虫 (小有名气)

结构工程,泥石流拦挡结构
3楼2018-06-04 18:44:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

shgao20

专家顾问 (职业作家)

【答案】应助回帖

感谢参与,应助指数 +1
ANSYS  静力分析也可以实现转动,只是你的做法有问题。
  看你的命令流,你是先将圆柱的一个端面划分 shell181 单元,再扫描成为 solid185 单元,同时清除了壳体单元,只留下实体单元。
  不过,实体单元的节点只有三个平移自由度,没有转动自由度。因此你对节点的转动自由度施加的约束 (0 或非 0 约束) 都是无效的。
  计算时出现的错误说明约束不足,经过分析,缺少约束与转动自由度无关,是对平移自由度的约束不足。分析如下:
    你选取用来施加约束的节点是圆柱的轴线,约束了所有平移自由度,不过,由于轴线是一条直线,因此不能约束圆柱绕轴线的转动,即少了一个约束条件。
  为了约束绕轴线的转动为某个指定值 (比如 PI/4, PI = 3.1416),需要将整个圆柱所有节点的节点坐标系转换为圆柱坐标系,然后对所有节点分别施加位移载荷:
      Uyi = Ri * PI/4   其中  i - 节点号
    节点很多是,可以用循环语句来实现。
4楼2018-06-05 08:44:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

liyaobo0509

专家顾问 (著名写手)

【答案】应助回帖

感谢参与,应助指数 +1
实体单元的节点只有三个平动自由度,可查看帮助文档对单元的说明。

ANSYS隐式求解完全可以。
如果要实现圆柱体转动,可另定义一个自由节点,建立节点与圆柱体的耦合关系,最后约束节点的平动自由度、施加转动约束,即可让圆柱转动起来。

» 本帖已获得的红花(最新10朵)

5楼2018-06-05 13:12:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Ampare1987

新虫 (小有名气)

送红花一朵
引用回帖:
5楼: Originally posted by liyaobo0509 at 2018-06-05 13:12:34
实体单元的节点只有三个平动自由度,可查看帮助文档对单元的说明。

ANSYS隐式求解完全可以。
如果要实现圆柱体转动,可另定义一个自由节点,建立节点与圆柱体的耦合关系,最后约束节点的平动自由度、施加转动约 ...

在尝试了很多的解决方案之后,我觉得你的说法最靠谱的,请教如何建立节点和圆柱体的耦合关系,非常感谢!
结构工程,泥石流拦挡结构
6楼2018-06-06 16:27:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Ampare1987

新虫 (小有名气)

引用回帖:
7楼: Originally posted by liyaobo0509 at 2018-06-06 19:48:13
实践检验过很多次,当然靠谱~可以试下MPC

非常感谢您再次回复我的问题,我现在就尝试使用MPC,在未得到您的关于耦合的指导前,我尝试了使用CP的方法,命令流如下
N,1000,0,0,0.01,0,0,0
cp,1,UX,1000,1
cp,2,UY,1000,1
CP,3,UZ,1000,1
CP,4,ROTX,1000,1
CP,5,ROTY,1000,1
CP,6,ROTZ,1000,1
cp,1,UX,all
cp,2,UY,all
CP,3,UZ,all
CP,4,ROTX,all
CP,5,ROTY,all
CP,6,ROTZ,all

D,1000, ,0, , , ,UX, , , , ,
D,1000, ,0, , , ,UY, , , , ,
D,1000, ,0, , , ,UZ, , , , ,
D,1000,ROTX,0
D,1000,ROTY,0
D,1000,ROTZ,3.1415926/4
结果出现了如下警告
Coupled node set 4 has unused node 1000.                                
Coupled node set 5 has unused node 1000.                                
Coupled node set 6 has unused node 1000.                                
Specified degree of freedom constraint UX at unused node 1000.         
Specified degree of freedom constraint UY at unused node 1000.         

Specified degree of freedom constraint UZ at unused node 1000.         
Specified degree of freedom constraint ROTX at unused node 1000.        
Specified degree of freedom constraint ROTY at unused node 1000.        
Specified degree of freedom constraint ROTZ at unused node 1000.        
我在网上看到了如下解释:
一般是节点并不属于单元(此节点不需要参与计算),然而还是被约束了,会出现这种警告。需要检查模型及约束。
我觉得这个解释还是有道理的,我不太明白您说的自由节点是什么意思,该如何建呢?
结构工程,泥石流拦挡结构
8楼2018-06-07 10:14:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

liyaobo0509

专家顾问 (著名写手)

引用回帖:
8楼: Originally posted by Ampare1987 at 2018-06-07 10:14:50
非常感谢您再次回复我的问题,我现在就尝试使用MPC,在未得到您的关于耦合的指导前,我尝试了使用CP的方法,命令流如下
N,1000,0,0,0.01,0,0,0
cp,1,UX,1000,1
cp,2,UY,1000,1
CP,3,UZ,1000,1
CP,4,ROTX,1000 ...

没记错的话,CE、CP耦合集中的任何节点,不可以施加位移约束或载荷,可以查看帮助资料确认一下
9楼2018-06-07 12:18:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Ampare1987

新虫 (小有名气)

引用回帖:
4楼: Originally posted by shgao20 at 2018-06-05 08:44:50
ANSYS  静力分析也可以实现转动,只是你的做法有问题。
  看你的命令流,你是先将圆柱的一个端面划分 shell181 单元,再扫描成为 solid185 单元,同时清除了壳体单元,只留下实体单元。
  不过,实体单元的节点只 ...

谢谢您的回答,对节点施加位移不可能施加在柱坐标系上,因为ANSYS只能对节点坐标系施加位移。这是我尝试了很长时间才得出的结论。您的前半部分的回答对我还是很有帮助的,再次感谢您!
结构工程,泥石流拦挡结构
10楼2018-06-13 09:42:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 Ampare1987 的主题更新
信息提示
请填处理意见