24小时热门版块排行榜    

查看: 3384  |  回复: 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一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 vZfe6xYu34yj 2026-09-14 3/150 2026-09-18 02:16 by cITaGg3p5edr
[考博] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +5 Aj1rhIDL5ixY 2026-09-14 6/300 2026-09-18 01:54 by cITaGg3p5edr
[考博] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 Aj1rhIDL5ixY 2026-09-14 3/150 2026-09-18 01:40 by cITaGg3p5edr
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +4 40ms4Wlo4umh 2026-09-14 4/200 2026-09-18 01:04 by cITaGg3p5edr
[硕博家园] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +7 LwdutQ8HoqWP 2026-09-13 7/350 2026-09-17 22:53 by cuPZTDXS3VOd
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +6 3n8v2C8RimXI 2026-09-13 7/350 2026-09-17 22:51 by cuPZTDXS3VOd
[基金申请] 国社科系统bug了,是不是要放榜了? +5 kynobel 2026-09-16 6/300 2026-09-17 21:37 by tosogo
[考研] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +6 Aj1rhIDL5ixY 2026-09-14 6/300 2026-09-17 14:17 by ObVzqQQh1rrE
[硕博家园] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +6 40ms4Wlo4umh 2026-09-14 6/300 2026-09-17 13:53 by ObVzqQQh1rrE
[硕博家园] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +6 23jxep3nCNZb 2026-09-14 6/300 2026-09-17 13:27 by ObVzqQQh1rrE
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +6 6F5UbRU2I5hL 2026-09-14 6/300 2026-09-17 03:04 by Tql5LQhh5rLK
[公派出国] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +4 LwdutQ8HoqWP 2026-09-13 7/350 2026-09-17 02:02 by Tql5LQhh5rLK
[考研] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +7 QUjhNVAcOSff 2026-09-13 8/400 2026-09-17 01:16 by Tql5LQhh5rLK
[考博] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +4 6F5UbRU2I5hL 2026-09-14 4/200 2026-09-16 17:00 by i3KBRD19ERkG
[考博] 上海工程技术大学激光智能制造课题组|2027级博士研究生招生公告 +7 水士口 2026-09-16 8/400 2026-09-16 13:47 by kkklens
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +7 QUjhNVAcOSff 2026-09-13 8/400 2026-09-16 02:23 by DgNGHc3h5tPl
[考博] 售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,科目全,可伽急 +5 s3fFTmArrBt6 2026-09-14 5/250 2026-09-15 16:59 by CgyNCDVNhGVg
[教师之家] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 6F5UbRU2I5hL 2026-09-14 3/150 2026-09-15 08:26 by C79jjtjAKjEn
[教师之家] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 s3fFTmArrBt6 2026-09-14 3/150 2026-09-15 07:46 by C79jjtjAKjEn
信息提示
请填处理意见