应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 仿真模拟 » ANSYS & ABAQUS » ANSYS梁单元与实体单元的耦合与约束方程
51/1返回列表
查看: 4742  |  回复: 6
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

小左儿

银虫 (正式写手)

[交流] 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
回复此楼

» 猜你喜欢

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

查看全部散金贴
1楼 2016-03-26 10:26:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

李晨奇

木虫 (正式写手)
★
小左儿(金币+1): 谢谢参与
好帖子，学习一下。

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

swir

金虫 (小有名气)
★
小左儿(金币+1): 谢谢参与
马克

发自小木虫Android客户端
回复此楼
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   回复此楼   编辑   查看我的主页
查看全部 7 个回答
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[博后之家] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +3 Aj1rhIDL5ixY 2026-09-14 3/150 2026-09-15 20:26 by zNcEhTcH7ihq
[考博] 上海工程技术大学激光智能制造课题组招收博士研究生 +6 两三岁ss 2026-09-14 6/300 2026-09-15 20:17 by 只吃肉葡萄
[公派出国] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +3 Aj1rhIDL5ixY 2026-09-14 3/150 2026-09-15 19:59 by CgyNCDVNhGVg
[硕博家园] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +3 40ms4Wlo4umh 2026-09-14 3/150 2026-09-15 18:59 by CgyNCDVNhGVg
[找工作] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +3 23jxep3nCNZb 2026-09-14 3/150 2026-09-15 18:47 by CgyNCDVNhGVg
[考博] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +3 6F5UbRU2I5hL 2026-09-14 3/150 2026-09-15 17:39 by CgyNCDVNhGVg
[找工作] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +3 s3fFTmArrBt6 2026-09-13 3/150 2026-09-15 15:50 by CgyNCDVNhGVg
[考博] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +6 QUjhNVAcOSff 2026-09-13 6/300 2026-09-15 14:25 by VMDqgnlKimtZ
[考研] 售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急 +6 QUjhNVAcOSff 2026-09-13 7/350 2026-09-15 14:22 by VMDqgnlKimtZ
[博后之家] 售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急 +6 QUjhNVAcOSff 2026-09-13 7/350 2026-09-15 14:13 by VMDqgnlKimtZ
[公派出国] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +7 QUjhNVAcOSff 2026-09-13 9/450 2026-09-15 14:10 by VMDqgnlKimtZ
[硕博家园] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +3 23jxep3nCNZb 2026-09-14 3/150 2026-09-15 09:38 by VMDqgnlKimtZ
[教师之家] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +3 6F5UbRU2I5hL 2026-09-14 3/150 2026-09-15 08:14 by C79jjtjAKjEn
[教师之家] 售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急 +4 s3fFTmArrBt6 2026-09-13 4/200 2026-09-15 07:02 by 5BDX0d0WFp7t
[考研] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +3 s3fFTmArrBt6 2026-09-13 5/250 2026-09-15 06:49 by 5BDX0d0WFp7t
[考研] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +4 LwdutQ8HoqWP 2026-09-13 4/200 2026-09-15 06:14 by 5BDX0d0WFp7t
[博后之家] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +6 3n8v2C8RimXI 2026-09-13 6/300 2026-09-15 05:50 by 5BDX0d0WFp7t
[教师之家] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +6 3n8v2C8RimXI 2026-09-13 6/300 2026-09-15 05:38 by 5BDX0d0WFp7t
[基金申请] 求教各位大神：2026教育部人文社科青年基金项目何时公示呀？ +5 云宴山人 2026-09-10 7/350 2026-09-14 12:11 by ecnu2013
[找工作] 浙江师范大学是怎么坑我的 +6 recruit123 2026-09-09 6/300 2026-09-14 12:00 by 萨仁其其格
信息提示
关闭
请填处理意见
关闭