24小时热门版块排行榜    

查看: 3363  |  回复: 1

achillesue

新虫 (初入文坛)

[求助] 有关MPC多点约束的实例的一个问题

最近在研究ANSYS中beam-solid的固接的连接方法,分别用MPC接触法和约束方程法做了一个简单的例子,但是得到的结果完全不一样,有大神能帮我看看是什么问题吗?

以下是MPC接触方法的命令流:

FINI
/CLE
/FILNAME,BEAM_AND_SOLID_ELEMENTS_CONNECTION !定义工作文件名
/TITLE,COUPLE_AND_CONSTRAINT_EQUATION !定义工作名
/PREP7
ET,1,SOLID95 !定义实体单元类型为SOLID95
ET,2,BEAM4 !定义梁单元类型为BEAM4
Et,3,170
Et,4,175
Keyopt,4,2,2
Keyopt,4,12,5
Keyopt,4,4,1
Keyopt,3,110111
MP,EX,1,3E4 !定义材料的弹性模量
MP,PRXY,1,0.3 !定义泊松比
R,1 !定义实体单元实常数
R,2,10.0,10/12.0,1000/12.0,10.0,1.0 !定义梁单元实常数
BLC4,,,20,7,10 !创建矩形块为实体模型
WPOFFS,0,3.5 !将工作平面向Y方向移动3.5
WPROTA,0,90 !将工作平面绕X轴旋转90度
VSBW,ALL !将实体沿工作平面剖开
WPOFFS,0,5 !将工作平面向Y方向移动5
WPROTA,0,90 !将工作平面绕X轴旋转90度
VSBW,ALL !将实体沿工作平面剖开
WPCSYS,-1 !将工作平面设为与总体笛卡儿坐标一致
K,100,20,3.5,5 !创建关键点
K,101,120,3.5,5 !创建关键点
L,100,101 !连接关键点生成梁的线实体
LSEL,S,LOC,X,21,130 !选择梁线
LATT,1,2,2 !指定梁的单元属性
LESIZE,ALL,,,10 !指定梁上的单元份数
LMESH,ALL !划分梁单元
VSEL,ALL !选择所有实体
VATT,1,1,1 !设置实体的单元属性
ESIZE,1 !指定实体单元尺寸
MSHAPE,0,2D !设置实体单元为2D
MSHKEY,1 !设置为映射网格划分方法
VMESH,ALL !划分实体单元
ALLS !全选
FINI !退出前处理
/SOLU !进入求解器
ASEL,S,LOC,X,0 !选择实体的端面
DA,ALL,ALL !约束实体端面
ALLS !全选
FK,101,FY,-3.0 !在两端施加Y向压力
R,3
Real,3
Type,3
TSHAP,PILO  
E,1
Esel,s,type,,1
Nsle,s
Nsel,r,loc,x,20
Type,4
Real,3
Esurf,all
Alls
Solve

以下是约束方程法:
FINI

/CLE

/FILNAME,BEAM_AND_SOLID_ELEMENTS_CONNECTION !定义工作文件名

/TITLE,COUPLE_AND_CONSTRAINT_EQUATION !定义工作名

/PREP7

ET,1,SOLID95 !定义实体单元类型为SOLID95

ET,2,BEAM4 !定义梁单元类型为BEAM4

MP,EX,1,3E4 !定义材料的弹性模量

MP,PRXY,1,0.3 !定义泊松比

R,1 !定义实体单元实常数

R,2,10.0,10/12.0,1000/12.0,10.0,1.0 !定义梁单元实常数

BLC4,,,20,7,10 !创建矩形块为实体模型

WPOFFS,0,3.5 !将工作平面向Y方向移动3.5

WPROTA,0,90 !将工作平面绕X轴旋转90度

VSBW,ALL !将实体沿工作平面剖开

WPOFFS,0,5 !将工作平面向Y方向移动5

WPROTA,0,90 !将工作平面绕X轴旋转90度

VSBW,ALL !将实体沿工作平面剖开

WPCSYS,-1 !将工作平面设为与总体笛卡儿坐标一致

K,100,20,3.5,5 !创建关键点

K,101,120,3.5,5 !创建关键点

L,100,101 !连接关键点生成梁的线实体

LSEL,S,LOC,X,21,130 !选择梁线

LATT,1,2,2 !指定梁的单元属性

LESIZE,ALL,,,10 !指定梁上的单元份数

LMESH,ALL !划分梁单元

VSEL,ALL !选择所有实体

VATT,1,1,1 !设置实体的单元属性

ESIZE,1 !指定实体单元尺寸

MSHAPE,0,2D !设置实体单元为2D

MSHKEY,1 !设置为映射网格划分方法

VMESH,ALL !划分实体单元

ALLS !全选

FINI !退出前处理

/SOLU !进入求解器

ASEL,S,LOC,X,0 !选择实体的端面

DA,ALL,ALL !约束实体端面

ALLS !全选

FK,101,FY,-3.0 !在两端施加Y向压力

CP,1,UX,1,21 !耦合节点1和节点21X方向自由度

CP,2,UY,1,21 !耦合节点1和节点21Y方向自由度

CP,3,UZ,1,21 !耦合节点1和节点21Z方向自由度

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)) !设置约束方程

ALLS !全选

SOLVE !保存

谢谢各位啦!
回复此楼

» 收录本帖的淘帖专辑推荐

ANSYS命令流

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)

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

TorresLIN

金虫 (小有名气)

朋友,你知道这个问题怎么解决了吗?我也在研究怎么做MPC,谢谢你
2楼2019-09-12 16:52:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 achillesue 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 投票:  有多少人是今天查系统知道结果的? +14 爱看书的可乐 2026-08-26 16/800 2026-08-26 17:44 by skyjudy
[基金申请] 范进中举一文的中心思想 +9 炎黄贵胄 2026-08-22 10/500 2026-08-26 15:40 by semaglutide
[基金申请] 我不理解! +8 Edward_pc 2026-08-26 13/650 2026-08-26 15:37 by Tide man
[基金申请] 2026国自然函评费到账 +22 羊腰板 2026-08-21 25/1250 2026-08-26 15:00 by zuocuiping
[基金申请] 出来了 +9 trojank 2026-08-26 9/450 2026-08-26 14:25 by 宝贝虫子
[基金申请] 哪位高人中了,把查询到的截图贴出来让我看看,让我长长见识 +4 yuleib84 2026-08-26 5/250 2026-08-26 13:40 by yuleib84
[基金申请] 为什么国自然不能直接公布 +4 bjdxyxy 2026-08-26 4/200 2026-08-26 13:12 by qingmu1201
[基金申请] 系统进不去 +4 yanglien 2026-08-26 5/250 2026-08-26 11:10 by wenfengw83
[基金申请] 项目信息和经费信息在系统里都可以看到了 +6 wittyboy 2026-08-26 14/700 2026-08-26 10:55 by wittyboy
[基金申请] 国合可查了 +3 paperzjh 2026-08-26 3/150 2026-08-26 10:41 by LemmonTr
[基金申请] 明天应该可查了!? +6 chengyan1220 2026-08-23 6/300 2026-08-25 19:45 by zfd97
[基金申请] 某些机构,以效率低为荣,以效率低作为存在感 +9 yuleib84 2026-08-25 10/500 2026-08-25 17:14 by alexon
[基金申请] 没有任何消息-是不是就凉了 +9 图啦图啦 2026-08-24 10/500 2026-08-25 11:59 by 南海小哥
[教师之家] 导师吐槽:我怎么摊上了这么个极品研究生! +3 苏东坡二世 2026-08-23 3/150 2026-08-25 10:35 by shisan1313
[基金申请] 我面上完蛋了 +13 且听虎啸 2026-08-20 14/700 2026-08-25 09:10 by mrkang
[基金申请] 建议基金发布提前给出明确的时间点 +13 kulium 2026-08-21 16/800 2026-08-24 16:27 by superceng
[基金申请] 让我中一个面上吧! +13 大萍1987 2026-08-20 16/800 2026-08-24 10:23 by 太傻了
[基金申请] 时间戳今天,20号变了 +5 archvillain 2026-08-20 5/250 2026-08-22 06:12 by hui_daxiao
[基金申请] 时间戳又变了 +13 wuchongjun 2026-08-20 19/950 2026-08-21 17:21 by 紫杉醇
[基金申请] 基金啊基金 +4 longfie172 2026-08-20 4/200 2026-08-21 08:58 by mark mao
信息提示
请填处理意见