24小时热门版块排行榜    

查看: 3366  |  回复: 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 !保存

谢谢各位啦!
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

TorresLIN

金虫 (小有名气)

朋友,你知道这个问题怎么解决了吗?我也在研究怎么做MPC,谢谢你
2楼2019-09-12 16:52:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 achillesue 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[论文投稿] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +5 ASdOkHsho7FD 2026-08-28 6/300 2026-08-30 17:25 by hZiFeudZyoGR
[基金申请] 面上意见出来了 +9 黄鸟于飞Chao 2026-08-29 18/900 2026-08-30 16:47 by 黄鸟于飞Chao
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 jCd0dEvKHShX 2026-08-29 6/300 2026-08-30 13:54 by GkCVyODlx4mR
[找工作] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +6 ASdOkHsho7FD 2026-08-28 9/450 2026-08-30 11:19 by l0VvVHGBGRLv
[基金申请] 29号明天会评吗 +3 笨笨唐 2026-08-28 3/150 2026-08-30 09:29 by cww8181
[论文投稿] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 G6APbkg8SA6w 2026-08-29 4/200 2026-08-30 07:39 by ZPa0EcMwuECS
[考博] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +6 ASdOkHsho7FD 2026-08-28 10/500 2026-08-30 06:11 by ZPa0EcMwuECS
[教师之家] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +5 ASdOkHsho7FD 2026-08-28 8/400 2026-08-30 05:48 by ZPa0EcMwuECS
[考研] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 G6APbkg8SA6w 2026-08-29 5/250 2026-08-30 02:40 by ZPa0EcMwuECS
[基金申请] 为什么到现在没收到通知? +4 tannykie 2026-08-29 4/200 2026-08-29 17:14 by lmz0216
[教师之家] 导师吐槽:我怎么摊上了这么个极品研究生! +9 苏东坡二世 2026-08-23 9/450 2026-08-29 14:41 by hustersqt
[基金申请] 基金不中,共勉 +11 eulota 2026-08-26 11/550 2026-08-28 14:22 by 火星超人xi
[基金申请] 基金系统什么内容也没有 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 zhaosm1982 2026-08-27 6/300 2026-08-27 16:00 by lfy8008
[文学芳草园] 梦想 +3 myrtle 2026-08-26 3/150 2026-08-27 10:01 by angelyueyi
[基金申请] 国合里面能看到了 +7 一怀馨秋 2026-08-26 7/350 2026-08-26 11:23 by zhaosm1982
[基金申请] 系统进不去 +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
信息提示
请填处理意见