| 查看: 4560 | 回复: 6 | |||
[交流]
ANSYS梁单元与实体单元的耦合与约束方程
|
|||
|
ANSYS梁单元与实体单元的耦合与约束方程 1 概述 在ANSYS计算过程中，有时候会遇到不同单元之间进行连接，由于不同的单元自由度不同，连接时通常需要通过耦合和约束方程建立节点自由度的关系，保证结果的准确性。 耦合可以理解成是将耦合的对象某个自由度作相等处理，而约束方程则不局限于相等这个关系，其可以描述具有某种关系的自由度。如图1所示，为梁单元与平面单元的连接。如果不采用约束方程，力矩的传递无法完成，因为平面单元没有转动自由度。  图1 梁单元与平面单元连接 为使节点2具有力矩传递的能力，要求1、2、3节点之间的自由度满足以下关系： ROTZ2 = (UY3 - UY1)/10 再通过CE命令，即可将此关系通过约束方程的形式施加给1、2、3节点。 2 命令 查看ANSYS的帮助文档，查询CE命令的解释，如图2所示。  图2 ANSYS的CE命令解释 CE, NEQN, CONST, NODE1, Lab1, C1, NODE2, Lab2, C2, NODE3, Lab3, C3 其中，NEQT表示常数，用于区别约束方程，一般可以用数字1、2、3表示即可，表示第几个约束方程； CONST表示方程的常数项，一般为0； NODE1，表示第一个节点； Lab1，表示自由度标签，对于结构而言，就是三个平移和三个转动自由度； C1，表示该自由度的系数； 同理，后面的也一样。 对于方程ROTZ2 = (UY3 - UY1)/10 稍作变形，0 = UY3 - UY1 - 10*ROTZ2 由此式即可直接写出对应的ANSYS命令流： CE,1,0,3,UY,1,1,UY,-1,2,ROTZ,-10 3 模型 本次为梁单元与实体单元的连接，建立模型，如图3所示，梁单元与实体单元有一个节点位置重合，为使位移和力矩能够传递，需要耦合两个节点的三个平移自由度，同时还需要用约束方程限制梁的三个转动自由度。  图3 梁单元与实体单元 4 约束方程 节点自由度耦合比较好操作，采用CP命令，重合位置处的两个节点分别为节点1（梁）和节点21（实体），自由度耦合如下： CP,1,UX,1,21 !耦合节点1和节点21X方向自由度 CP,2,UY,1,21 !耦合节点1和节点21Y方向自由度 CP,3,UZ,1,21 !耦合节点1和节点21Z方向自由度 为约束转动自由度，由CE的参数项可知，需要先写出转动约束方程，对照图4分别写出三个转动自由度的约束方程，图4中红圈的四个节点分布在中心节点周围，将这几个节点进行约束即可限制梁单元和实体单元的转动自由度，自由度方程如下：  图4 节点分布 ROTY(1)=(UX(626)-UX(2328))/ABS(NZ(626)-NZ(2328)) !Y轴转动 ROTZ(1)=(UX(67)-UX(4283))/ABS(NY(67)-NY(4283)) !Z轴转动 ROTX(1)=(UZ(67)-UZ(4283))/ABS(NY(67)-NY(4283)) !X轴转动 稍作变换如下： 0=(UX(626)-UX(2328))- ROTY(1)*ABS(NZ(626)-NZ(2328)) 0=(UX(67)-UX(4283))- ROTZ(1)* ABS(NY(67)-NY(4283)) 0=(UZ(67)-UZ(4283))- ROTX(1)* ABS(NY(67)-NY(4283)) 对应的CE命令为： CE,1,0,626,UX,1,2328,UX,-1,1,ROTY,-ABS(NZ(626)-NZ(2328)) CE,2,0,67,UX,1,4283,UX,-1,1,ROTZ,-ABS(NY(67)-NY(4283)) CE,3,0,67,UZ,1,4283,UZ,-1,1,ROTX,-ABS(NY(67)-NY(4283)) 施加完成后如图5所示。  图5 施加耦合、约束方程之后  |
回复此楼» 猜你喜欢
关于filecode，很负责任的告诉大家
已经有7人回复
-
【2027博士申请】纳米药物递送方向
已经有5人回复
-
关于豆爷回答的JTJC与%2F数量
已经有7人回复
-
2026国自然放榜时间
已经有3人回复
-
国基金的申报应该改成非等额制，评价高的钱多评价低的钱少，但是增加资助率
已经有7人回复
-
好奇怪的filecode
已经有4人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有4人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急
已经有5人回复
» 抢金币啦！回帖就可以得到:
-
一边期待，一边忐忑：发榜前的科研人正在想什么？
+5/120
-
别研究filecode了，祝都有一个好结果：）
+2/98
-
三(三甲基硅基)甘油醚（CAS：6787-10-6 ）
+1/88
-
1,3,2-Dioxaphospholane, 2-(2,2,2-trifluoroethoxy)-
+1/87
-
三氟甲基亚磺酸锂
+1/87
-
3,4-二(三甲基硅氧基)-1-丁烯（CAS：89543-84-0）
+1/87
-
乙二醇二甲基环硅醚（CAS：57915-66-9）
+1/87
-
2-氯-1,3,2-二氧磷杂环戊烷（CAS：822-39-9）
+1/86
-
香港城市大学郑星课题组诚聘博士后 (生态学、环境科学、植物学、大气科学背景)
+1/78
-
哈工大（深圳）何自开团队诚招化学工程与技术专业博士生
+1/27
-
中国石油大学（华东）材料科学与工程学院/新颖能源材料智能设计和开发方向 特聘副教授
+2/13
-
澳大利亚 Murdoch University 全奖博士招生（3个名额）地质化工冶金领域
+1/12
-
苏州大学能源学院管大秦课题组诚聘博士后及招收博士/硕士研究生
+1/11
-
中国石油大学（华东）计算机科学与技术学院/智能化新能源系统设计和开发方向 特聘副教
+2/10
-
澳大利亚 Murdoch University 全奖博士招生（3个名额）地质化工冶金领域
+1/9
-
中山大学李国强课题组招收博士生（柔性电子方向）
+1/8
-
香港理工大学殷骏课题组诚聘1 名博士后以及2 名博士生
+1/5
-
北京理工大学-集成电路与电子学院杰青团队-招博士后
+1/2
-
氨氮空白偏高！求指教！
+1/2
-
上海交通大学化学化工学院张智涛课题组诚聘博士后
+1/1
3楼2016-04-02 08:32:52
4楼2016-04-02 09:25:32
5楼2016-08-23 17:07:44
6楼2016-11-03 22:18:48
7楼2018-11-14 14:54:06
简单回复
2016-04-02 08:02
回复
小左儿(金币+1): 谢谢参与
