24小时热门版块排行榜    

查看: 6222  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 投票:  有多少人是今天查系统知道结果的? +17 爱看书的可乐 2026-08-26 19/950 2026-08-30 22:24 by xiangy672
[基金申请] 面上意见出来了 +9 黄鸟于飞Chao 2026-08-29 18/900 2026-08-30 16:47 by 黄鸟于飞Chao
[考博] 找导师 +6 yuanjiabao 2026-08-29 7/350 2026-08-30 14:40 by 生科新手
[基金申请] 29号明天会评吗 +3 笨笨唐 2026-08-28 3/150 2026-08-30 09:29 by cww8181
[基金申请] 麻烦专家们看看评委们的意见(F口面上) +6 gdd2018 2026-08-28 11/550 2026-08-30 08:58 by 超级无敌华子
[基金申请] 我就是申请一个面上项目而已,这评审意见是按照杰青的条件评的吧? +6 gouxfjh 2026-08-28 11/550 2026-08-30 07:57 by gouxfjh
[基金申请] 国自然评审意见 +13 wangmingqi 2026-08-28 19/950 2026-08-29 10:22 by Poppy1104
[基金申请] 2026年叶企孙基金 +4 bud_bud 2026-08-27 7/350 2026-08-29 07:23 by foolishmani
[基金申请] 为什么资助数各大高校都创新高,自己申请怎么就这么难 +12 Kittylucky 2026-08-27 13/650 2026-08-29 00:04 by superceng
[基金申请] 基金系统什么内容也没有 30+4 winsaint 2026-08-27 9/450 2026-08-28 11:06 by maolC
[基金申请] 怎么查啊 +6 huang1991js 2026-08-26 6/300 2026-08-28 08:42 by winsaint
[基金申请] 面上合作单位盖章 +5 ssyjh 2026-08-27 5/250 2026-08-27 20:50 by gdfollow
[基金申请] 怎么看青基中了没有啊 +5 叶九微 2026-08-26 5/250 2026-08-27 10:35 by l_zh2008
[基金申请] 出来了 +9 trojank 2026-08-26 9/450 2026-08-26 14:25 by 宝贝虫子
[基金申请] 系统进不去 +4 yanglien 2026-08-26 5/250 2026-08-26 11:10 by wenfengw83
[基金申请] 国合可查了 +3 paperzjh 2026-08-26 3/150 2026-08-26 10:41 by LemmonTr
[基金申请] 今天务委会开完了,明天出结果吗 +19 angus9576 2026-08-25 23/1150 2026-08-26 10:03 by zp519
[基金申请] 国合现在查不到了吗? +10 chengyan1220 2026-08-24 20/1000 2026-08-26 08:57 by peasantsprig
[基金申请] 某些机构,以效率低为荣,以效率低作为存在感 +9 yuleib84 2026-08-25 10/500 2026-08-25 17:14 by alexon
[基金申请] 如果此刻你正在为国基感到焦虑,不妨来听听这首《基金之外》 +8 scalable 2026-08-24 8/400 2026-08-25 12:52 by jnhyjjm
信息提示
请填处理意见