24小时热门版块排行榜    

查看: 4562  |  回复: 6

小左儿

银虫 (正式写手)


[交流] ANSYS梁单元与实体单元的耦合与约束方程

ANSYS梁单元与实体单元的耦合与约束方程
1 概述
在ANSYS计算过程中,有时候会遇到不同单元之间进行连接,由于不同的单元自由度不同,连接时通常需要通过耦合和约束方程建立节点自由度的关系,保证结果的准确性。
耦合可以理解成是将耦合的对象某个自由度作相等处理,而约束方程则不局限于相等这个关系,其可以描述具有某种关系的自由度。如图1所示,为梁单元与平面单元的连接。如果不采用约束方程,力矩的传递无法完成,因为平面单元没有转动自由度。
ANSYS梁单元与实体单元的耦合与约束方程
图1 梁单元与平面单元连接
为使节点2具有力矩传递的能力,要求1、2、3节点之间的自由度满足以下关系:
ROTZ2 = (UY3 - UY1)/10
再通过CE命令,即可将此关系通过约束方程的形式施加给1、2、3节点。
2 命令
查看ANSYS的帮助文档,查询CE命令的解释,如图2所示。
ANSYS梁单元与实体单元的耦合与约束方程-1
图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所示,梁单元与实体单元有一个节点位置重合,为使位移和力矩能够传递,需要耦合两个节点的三个平移自由度,同时还需要用约束方程限制梁的三个转动自由度。
ANSYS梁单元与实体单元的耦合与约束方程-2
图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中红圈的四个节点分布在中心节点周围,将这几个节点进行约束即可限制梁单元和实体单元的转动自由度,自由度方程如下:
ANSYS梁单元与实体单元的耦合与约束方程-3
图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所示。
ANSYS梁单元与实体单元的耦合与约束方程-4
图5 施加耦合、约束方程之后

CODE:
ANSYS梁单元与实体单元的耦合与约束方程-5

回复此楼

» 猜你喜欢

» 抢金币啦!回帖就可以得到:

查看全部散金贴

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

李晨奇

木虫 (正式写手)



小左儿(金币+1): 谢谢参与
好帖子,学习一下。

发自小木虫Android客户端
3楼2016-04-02 08:32:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

swir

金虫 (小有名气)



小左儿(金币+1): 谢谢参与
4楼2016-04-02 09:25:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

机械508

木虫 (小有名气)



小左儿(金币+1): 谢谢参与
这算是 两种方法的组合使用?
5楼2016-08-23 17:07:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

匿名

用户注销 (小有名气)



小左儿(金币+1): 谢谢参与
本帖仅楼主可见
6楼2016-11-03 22:18:48
已阅   申请仿真EPI   回复此楼   编辑   查看我的主页

runt22

新虫 (初入文坛)



小左儿(金币+1): 谢谢参与
楼主,使用CE命令的时候,要是约束方程里面有四个自由度,该怎么办啊。谢谢
7楼2018-11-14 14:54:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
简单回复
2016-04-02 08:02   回复  
小左儿(金币+1): 谢谢参与
相关版块跳转 我要订阅楼主 小左儿 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 关于豆爷回答的JTJC与%2F数量 +5 yang182083 2026-08-06 7/350 2026-08-08 18:29 by zhanghaozhu
[基金申请] 国基金的申报应该改成非等额制,评价高的钱多评价低的钱少,但是增加资助率 +7 a089 2026-08-07 7/350 2026-08-08 18:05 by gltch
[公派出国] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 qTvzQBRAHjCi 2026-08-07 4/200 2026-08-08 17:27 by oEVWOejN9taj
[基金申请] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 qTvzQBRAHjCi 2026-08-07 4/200 2026-08-08 17:22 by oEVWOejN9taj
[硕博家园] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 qTvzQBRAHjCi 2026-08-07 4/200 2026-08-08 17:19 by oEVWOejN9taj
[论文投稿] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 qTvzQBRAHjCi 2026-08-07 4/200 2026-08-08 16:47 by oEVWOejN9taj
[考博] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 Vi50GxzrFcSG 2026-08-07 4/200 2026-08-08 16:27 by oEVWOejN9taj
[基金申请] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 Vi50GxzrFcSG 2026-08-07 5/250 2026-08-08 16:27 by oEVWOejN9taj
[硕博家园] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 HEQlVqMTIA7d 2026-08-07 6/300 2026-08-08 14:47 by oEVWOejN9taj
[公派出国] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 HEQlVqMTIA7d 2026-08-07 6/300 2026-08-08 14:39 by oEVWOejN9taj
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +4 HEQlVqMTIA7d 2026-08-07 5/250 2026-08-08 14:27 by oEVWOejN9taj
[论文投稿] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 HEQlVqMTIA7d 2026-08-07 5/250 2026-08-08 14:22 by oEVWOejN9taj
[教师之家] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 KXLV3nuBVBY7 2026-08-07 5/250 2026-08-08 14:07 by oEVWOejN9taj
[论文投稿] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 KXLV3nuBVBY7 2026-08-07 6/300 2026-08-08 14:01 by oEVWOejN9taj
[基金申请] 基金中了 +14 laoda193707 2026-08-06 14/700 2026-08-08 00:23 by 实验小白ha
[基金申请] 化学口download_prp&fileCode的固定段好像这几天一直没变,有变的大神么? +3 Tide man 2026-08-07 4/200 2026-08-07 22:39 by Tide man
[基金申请] 大家散了吧,后缀研究没有意义,别浪费时间了,过好目前的每一天,不要焦虑 +5 Tide man 2026-08-06 7/350 2026-08-07 13:11 by 医学老男孩
[基金申请] 求各位大神看下 100+6 hpkpkpkp 2026-08-05 33/1650 2026-08-06 14:49 by zhiyanjiang
[基金申请] 影响面上的因素 +8 布布和一二 2026-08-05 11/550 2026-08-06 10:41 by 宝贝虫子
[基金申请] 8月时间戳变的,举个手。玩一下,释放压力 +9 archvillain 2026-08-04 11/550 2026-08-05 20:06 by wlwhappy
信息提示
请填处理意见