| 查看: 3524 | 回复: 3 | ||||
求索者807新虫 (正式写手)
|
[求助]
接触分析，在求解后接触面脱离？？？？？？ 已有2人参与
|
|
模拟球压在面上，然后求接触面的接触应力，用二维分析，对称约束，求解后接触面脱离，一直没明白，求大神指导？ 球 目标面 TARGE169 盘 接触面 CONTA172 附命令流： /CLEAR /FILNAME, EXAMPLE11.12 /PREP7 ET, 1, PLANE182 Keyopt,1,3,1 !第三个数从0开始 ET, 2, TARGE169 ET, 3, CONTA172,,,,,1 R, 1,,,0.1 MP, EX, 1, 2E11 !球的弹性模量 MP, PRXY, 1, 0.3 MP, EX, 2, 1.74E11 !块的弹性模量 MP, PRXY, 2, 0.34 CYL4,0,0,3,-90,0,0 !建立一个圆 Rectng,0,6,-3,-8 Asel,,,,1 AATT, 1, 1, 1 !指定面的单元属性,材料1、实常数号是1，单元类型是1 Asel,,,,2 AATT, 2, 1, 1 ALLS,ALL ESIZE, 0.5 MSHAPE, 0 !指定单元形状 0 四边形（2D），六面体（3D） MSHKEY, 0 !0 自由网格划分 AMESH, ALL NREFINE,80, , , 3 ! 将节点附近的单元网格细化 NREFINE,3, , , 3 ! 将节点附近的单元网格细化 DL, 4, 1,UX,0 DL, 4, 1,UY,0 DL, 3, 1,UX !3,1分别为受约束线段及线段所属面积的号码 DL, 7, 2,UX !节点号12 F, 12, FY,-10 !在指定线上加荷载 ALLSEL,ALL ! *** 创建目标目单元 !命令流CM Lsel,s,,,1 !球上的线 NSLL, S, 1 !线:球上L1，块L6 Nsel,S,NODE,,1,10 Cm,N-contact,node Lsel,s,,,6 !块上的线 NSLL, S, 1 !选择线上的节点 Nsel,S,NODE,,200,208 Nsel,A,NODE,,69 Cm,N-target,node Alls !选择所有 Eplot !显示单元 !手动进行接触对的设置 /SOLU Nlgeom,1 Deltim,0.01,0.001,0.2 ALLS SOLVE FINISH /post1 Set,last Plnsol,s,eqv,0,1.0  结果图片.PNG |
回复此楼» 猜你喜欢
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有7人回复
-
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急
已经有10人回复
-
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急
已经有6人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有8人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有5人回复
-
售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急
已经有10人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有8人回复
-
国社科又开始会评了，不知道这次命运如何
已经有14人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有8人回复
-
麻烦专家们看看评委们的意见（F口面上）
已经有13人回复
» 本主题相关商家推荐: (我也要在这里推广)
2楼2015-11-14 01:10:29
shgao20
专家顾问 (职业作家)
-
专家经验: +809 - 仿真EPI: 1
- 应助: 1311 (讲师)
- 金币: 13381.4
- 红花: 267
- 帖子: 3531
- 在线: 473.4小时
- 虫号: 3644138
- 注册: 2015-01-14
- 专业: 动力学与控制
- 管辖: 仿真模拟
【答案】应助回帖
感谢参与，应助指数 +1
|
检查了你的模型，发现一些问题。最主要的问题是：虽然你的命令流中说：”手动进行接触对的设置”，但是你的命令流中没有该部分的内容，即实际上你没有创建接触，因而计算时上面的球受力后直接穿过下面的盘，去向远方了。 除此以外，你的模型还有其它问题： 1 你准备设置接触对用的两个 component，分开很远，即使设置接触对，也不可能接触上； 2 你的网格划分并不合理，小球的网格加密区与盘的网格加密区分开了一段距离，并没有对准； 3 载荷为节点 12 处的集中力，对收敛不利，一般不建议这样做，较好的办法是对小球上边界线施加压力载荷； 4 接触属性基本没有设置。对你的问题可能需要设置：初始调整 - 减少穿透和闭合间隙；将接触行为设置为不分离，以有利于收敛；此外，可以设置一定的摩擦系数，比较符合实际情况。 5 求解设置过于简单； 修改了你的命令流，可以收敛。命令流如下，从 !*** 开始为修改的部分： fini /CLEAR /FILNAME, EXAMPLE11.12 /PREP7 ET, 1, PLANE182 Keyopt,1,3,1 !第三个数从0开始 ET, 2, TARGE169 ET, 3, CONTA172,,,,,1 R, 1,,,0.1 MP, EX, 1, 2E11 !球的弹性模量 MP, PRXY, 1, 0.3 MP, EX, 2, 1.74E11 !块的弹性模量 MP, PRXY, 2, 0.34 CYL4,0,0,3,-90,0,0 !建立一个圆 Rectng,0,6,-3,-8 Asel,,,,1 AATT, 1, 1, 1 !指定面的单元属性,材料1、实常数号是1，单元类型是1 Asel,,,,2 AATT, 2, 1, 1 ALLS,ALL ESIZE, 0.5 MSHAPE, 0 !指定单元形状 0 四边形（2D），六面体（3D） MSHKEY, 0 !0 自由网格划分 AMESH, ALL NREFINE,80, , , 3 ! 将节点附近的单元网格细化 NREFINE,3, , , 3 ! 将节点附近的单元网格细化 DL, 4, 1,UX,0 DL, 4, 1,UY,0 DL, 3, 1,UX !3,1分别为受约束线段及线段所属面积的号码 DL, 7, 2,UX !节点号12 F, 12, FY,-10 !在指定线上加荷载 ALLSEL,ALL ! *** 创建目标目单元 !命令流CM Lsel,s,,,1 !球上的线 NSLL, S, 1 !线:球上L1，块L6 Nsel,S,NODE,,1,10 Cm,N-contact,node Lsel,s,,,6 !块上的线 NSLL, S, 1 !选择线上的节点 Nsel,S,NODE,,200,208 Nsel,A,NODE,,69 Cm,N-target,node Alls !选择所有 Eplot !显示单元 !手动进行接触对的设置 !**** 以下重新选择节点，创建 component，准备设置接触对 LSEL,S, , , 1 NSLL,S,1 !* FLST,5,21,1,ORDE,7 FITEM,5,1 FITEM,5,3 FITEM,5,-11 FITEM,5,433 FITEM,5,-437 FITEM,5,441 FITEM,5,-446 NSEL,R, , ,P51X CM,N_contact,NODE !* LSEL,S, , , 6 NSLL,S,1 !* FLST,5,26,1,ORDE,6 FITEM,5,72 FITEM,5,79 FITEM,5,-83 FITEM,5,146 FITEM,5,192 FITEM,5,-210 NSEL,R, , ,P51X CM,N_target,NODE ALLSEL,ALL !* !**** 以下创建接触对 ! /COM, CONTACT PAIR CREATION - START CM,_NODECM,NODE CM,_ELEMCM,ELEM CM,_KPCM,KP CM,_LINECM,LINE CM,_AREACM,AREA CM,_VOLUCM,VOLU ! /GSAV,cwz,gsav,,temp MP,MU,1,0.2 MAT,1 MP,EMIS,1,7.88860905221e-031 R,3 REAL,3 ET,4,169 ET,5,172 R,3,,,1.0,0.1,0, RMORE,,,1.0E20,0.0,1.0, RMORE,0.0,0,1.0,,1.0,0.5 RMORE,0,1.0,1.0,0.0,,1.0 KEYOPT,5,3,0 KEYOPT,5,4,2 KEYOPT,5,5,3 KEYOPT,5,7,0 KEYOPT,5,8,0 KEYOPT,5,9,0 KEYOPT,5,10,2 KEYOPT,5,11,0 KEYOPT,5,12,0 KEYOPT,5,2,0 ! Generate the target surface NSEL,S,,,N_TARGET CM,_TARGET,NODE TYPE,4 ESLN,S,0 ESURF CMSEL,S,_ELEMCM ! Generate the contact surface NSEL,S,,,N_CONTACT CM,_CONTACT,NODE TYPE,5 ESLN,S,0 ESURF ALLSEL ESEL,ALL ESEL,S,TYPE,,4 ESEL,A,TYPE,,5 ESEL,R,REAL,,3 ! /PSYMB,ESYS,1 ! /PNUM,TYPE,1 ! /NUM,1 ! EPLOT ESEL,ALL ESEL,S,TYPE,,4 ESEL,A,TYPE,,5 ESEL,R,REAL,,3 CMSEL,A,_NODECM CMDEL,_NODECM CMSEL,A,_ELEMCM CMDEL,_ELEMCM CMSEL,S,_KPCM CMDEL,_KPCM CMSEL,S,_LINECM CMDEL,_LINECM CMSEL,S,_AREACM CMDEL,_AREACM CMSEL,S,_VOLUCM CMDEL,_VOLUCM ! /GRES,cwz,gsav CMDEL,_TARGET CMDEL,_CONTACT ! /COM, CONTACT PAIR CREATION - END ! /MREP,EPLOT ALLSEL,ALL !* /AUTO,1 EPLOT !* FINISH !* /SOL !* ANTYPE,0 !**** 设置求解参数 Nlgeom,1 DELTIM,0.01,0.00001,0.2 LNSRCH,1 NCNV,2,0,0,0 cnvtol, f, , 0.001, L2, -1 !* OUTPR,BASIC,ALL, /GST,1 OUTRES,ALL,ALL, kbc,0 !* ALLS SAVE !* SOLVE SAVE !* FINISH /POST1 !* /EFACET,1 PLNSOL, U,SUM, 0,1.0 /wait,3 !* /EFACET,1 PLNSOL, S,EQV, 0,1.0 /wait,3 !* FINISH |
3楼2015-11-15 11:24:14
求索者807
新虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 3610.1
- 散金: 880
- 红花: 3
- 帖子: 898
- 在线: 148.1小时
- 虫号: 3044207
- 注册: 2014-03-12
- 性别: GG
- 专业: 机械摩擦学与表面技术
太感谢了，好心人
4楼2017-04-05 08:21:51
30