24小时热门版块排行榜    

查看: 1454  |  回复: 34

11050209

银虫 (正式写手)

[求助] 接触问题求助 已有1人参与

两块板通过螺栓连接,拉伸问题的求助,搞不懂什么原因导致算出来的应力特别小!由于第一次接触到接触分析,所以不是很懂,我觉得问题可能出在接触这一块,不知道大神能不能帮看下,还有一个可能就是我施加的约束不够?如果是约束的问题又该怎么加约束呢
附上命令:
/prep7
/TITLE, lashen  
!定义单元类型及材料属性
et,1,solid185
et,2,plane182
mp,ex,1,190e6
mp,ex,2,184e6
mp,ex,3,69e9
mp,prxy,1,0.34
mp,prxy,2,0.3
mp,prxy,3,0.33
!定义参数
l1=100/1000    !试样长度
l2=34/1000    !夹头长度
w1=24/1000    !试样宽度
w2=33.8/1000    !夹头宽度
h=1.5/1000    !试样厚度
l3=48/1000    !试样交叠部分长度
r=4.9/1000      !连接圆弧半径
r1=6.25/2/1000     !圆台下表面半径
r2=8.75/2/1000     !圆台上表面半径
r3=4/1000     !螺母半径
pload1=100000000*2*h*w2
pload2=-pload1
!建立几何模型(TA32)
block,-l3/2,-l3/4,0,w1/2,0,h
block,-l3/4,0,0,w1/2,0,h
block,-l3/4,0,-w1/2,0,0,h
block,-l3/2,-l3/4,-w1/2,0,0,h
block,0,l3/4,0,w1/2,0,h
block,l3/4,l3/2,0,w1/2,0,h
block,l3/4,l3/2,-w1/2,0,0,h
block,0,l3/4,-w1/2,0,0,h
wpoff,-12/1000,0,0  
cone,r1,r2,0,h
wpoff,24/1000,0,0
cone,r1,r2,0,h
vsel,s,,,1,4,1
allsel
vsbv,all,9
allsel
vsel,s,,,5,8,1
allsel
vsbv,all,10
numcmp,all
CSYS,0  
WPAVE,0,0,0
CSYS,0  
block,-37.1/1000,-24/1000,-w1/2,w1/2,0,h
k,73,-42/1000,w1/2+r,h
k,74,-76/1000,w1/2+r,h
k,75,-76/1000,-(w1/2+r),h
k,76,-42/1000,-(w1/2+r),h
larc,72,73,54,r
l,73,74
l,74,75
l,75,76
larc,76,69,54,r
lsel,s,,,125,129,1
lsel,a,,,120
al,all
allsel
vdrag,63,,,,,,121
block,-76/1000,-42/1000,-w2/2,w2/2,-h,0
numcmp,all
nummrg,all
a,73,82,79,76
vsba,10,8
numcmp,all
allsel
vglue,all
!建立几何模型(7050)
block,-l3/2,-l3/4,0,w1/2,-h,0
block,-l3/4,0,0,w1/2,-h,0
block,-l3/4,0,-w1/2,0,-h,0
block,-l3/2,-l3/4,-w1/2,0,-h,0
block,0,l3/4,0,w1/2,-h,0
block,l3/4,l3/2,0,w1/2,-h,0
block,l3/4,l3/2,-w1/2,0,-h,0
block,0,l3/4,-w1/2,0,-h,0
wpoff,-12/1000,0,0
cylind,0,r1,0,-h
wpoff,24/1000,0,0
cylind,0,r1,0,-h
vsel,s,,,14,16,1
vsel,a,,,9
allsel
vsbv,all,21
vsel,s,,,17,20,1
allsel
vsbv,all,22
CSYS,0  
WPAVE,0,0,0
CSYS,0  
block,24/1000,37.1/1000,-w1/2,w1/2,-h,0

k,139,42/1000,w1/2+r,0
k,140,76/1000,w1/2+r,0
k,141,76/1000,-(w1/2+r),0
k,142,42/1000,-(w1/2+r),0

larc,85,139,19,r
l,139,140
l,140,141
l,141,142
larc,82,142,24,r
lsel,s,,,186,187,1
lsel,a,,,192,193,1
lsel,a,,,167,197,30
al,all
vdrag,99,,,,,,178
block,42/1000,76/1000,-w2/2,w2/2,0,h
a,117,121,114,100
vsba,18,143
vsel,s,,,14,17,1
vsel,a,,,23,26,1
vsel,a,,,19,21,1
vsel,a,,,9
vglue,all
numcmp,all
!建立几何模型(1s627)
wpoff,-12/1000,0,0  
cone,r1,r2,0,h
cylind,0,r1,0,-(h+h)
cylind,r1,r3,-h,-(h+h)
allsel
vsel,s,,,25,27,1
vglue,all
wpoff,24/1000,0,0
cone,r1,r2,0,h
cylind,0,r1,0,-(h+h)
cylind,r1,r3,-h,-(h+h)
allsel
vsel,s,,,27,30,3
vsel,a,,,26
vglue,all
CSYS,0  
WPAVE,0,0,0
CSYS,0  
allsel
!划分网格(TA32)
MSHAPE,0,2D
MSHKEY,1
mat,1
lsel,s,,,7,8
lesize,all,,,6
lsel,s,,,74,76,2
lesize,all,,,4
ALLSEL,ALL  
AMAP,39,4,33,32,5
lsel,s,,,17,18,1
lesize,all,,,6
AMAP,43,5,32,31,9  
allsel
AMAP,47,9,31,34,13
AMAP,50,13,34,33,4  
lsel,s,,,4,6,2
lesize,all,,,4
ALLSEL,ALL  
AMAP,7,9,41,40,18  
AMAP,14,18,40,39,21
AMAP,22,21,39,42,25
AMAP,34,25,42,41,9  
lsel,s,,,35
lesize,all,,,3
vsel,s,,,1,8,1
vsweep,all
asel,s,,,39,43,4
asel,a,,,47,50,3
asel,a,,,7,14,7
asel,a,,,22,34,12
aclear,all
allsel
vsweep,13
lsel,s,,,93
lesize,all,,,2,2
allsel
lsel,s,,,97
lesize,all,,,2,0.5
vsweep,11
lsel,s,,,94,96,2
lesize,all,,,17
vsweep,12
vsweep,10
allsel
!划分网格(7050)
mat,3
AMAP,96,62,101,100,63     
AMAP,121,63,100,99,76   
AMAP,132,76,99,102,68     
AMAP,129,68,102,101,62  
AMAP,73,76,109,108,77
AMAP,114,77,108,107,83
AMAP,117,83,107,110,88
AMAP,111,88,110,109,76
lsel,s,,,125
lesize,all,,,3  
vsel,s,,,17,21,1
vsel,a,,,23,24,1  
vsel,a,,,9
vsweep,all  
asel,s,,,73,96,23
asel,a,,,111,117,3
asel,a,,,121,129,8
asel,a,,,132
aclear,all
allsel
lsel,s,,,218,220,2
lesize,all,,,7
vsweep,22
lsel,s,,,199,200,1
lesize,all,,,2,2
vsweep,14
lsel,s,,,185,187,2
lesize,all,,,17
vsweep,16
vsweep,15
!划分网格(1s627)
mat,2
lsel,s,,,231,234,1
lesize,all,,,12
vsweep,25
lsel,s,,,259,262,1
lesize,all,,,12
vsweep,28
lsel,s,,,255,258,1
lesize,all,,,12
vsweep,29
lsel,s,,,241,244,1
LESIZE,all, , ,12
vsweep,26
lsel,s,,,291,294,1
lesize,all,,,12
vsweep,31
lsel,s,,,287,290,1
lesize,all,,,12
vsweep,32
allsel
!建立接触对
mp,mu,4,0.3
mat,4
r,4
real,4
et,3,targe170
et,4,conta173
!generate the target surface
asel,s,,,135,136,1
cm,left_up_bolt,area
cmsel,s,left_up_bolt
type,3
nsla,s,1
esln,s,0
esurf,all
allsel
!generate the contact surface
asel,s,,,37,45,4
asel,a,,,48
cm,left_up_plate,area
cmsel,s,left_up_plate
type,4
nsla,s,1
esln,s,0
esurf,all
allsel
!generate the target surface
mat,4
r,5
real,5
et,5,targe170
et,6,conta173
asel,s,,,147,148,1
cm,left_down_bolt,area
cmsel,s,left_down_bolt
type,5
nsla,s,1
esln,s,0
esurf,all
allsel
!generate the contact surface
asel,s,,,94,97,3
asel,a,,,99,100,1
cm,left_down_plate,area
cmsel,s,left_down_plate,area
type,6
nsla,s,1
esln,s,0
esurf,all
allsel
!generate the target surface
mat,4
r,6
real,6
et,7,targe170
et,8,conta173
asel,s,,,140,141,1
cm,right_up_bolt,area
cmsel,s,right_up_bolt
type,7
nsla,s,1
esln,s,0
esurf,all
allsel
!generate the contact surface
asel,s,,,3,9,6
asel,a,,,19,32,13
cm,right_up_plate,area
cmsel,s,right_up_plate
type,8
nsla,s,1
esln,s,0
esurf,all
allsel
!generate the target surface
mat,4
r,7
real,7
et,9,targe170
et,10,conta173
asel,s,,,160,161,1
cm,right_down_bolt,area
cmsel,s,right_down_bolt
type,9
nsla,s,1
esln,s,0
esurf,all
allsel
!generate the contact surface
asel,s,,,56,74,18
asel,a,,,77,87,10
cm,right_down_plate,area
cmsel,s,right_down_plate
type,10
nsla,s,1
esln,s,0
esurf,all
allsel
!generate the target surface
mat,4
r,8
real,8
et,11,targe170
et,12,conta173
asel,s,,,6,13,7
asel,a,,,20,33,13
asel,a,,,38,46,4
asel,a,,,49
cm,left_plate,area
cmsel,s,left_plate
type,11
nsla,s,1
esln,s,0
esurf,all
allsel
!generate the contact surface
asel,s,,,129,132,3
asel,a,,,96,121,25
asel,a,,,111,117,6
asel,a,,,73,114,41
cm,right_plate,area
cmsel,s,right_plate
type,12
nsla,s,1
esln,s,0
esurf,all
allsel
/sol
antype,static
nlgeom,on
autots,on
lnsrcch,on
pred,on
outres,all,all
time,1
nsubst,80,1000,20
nsel,s,loc,x,76/1000,76/1000
cp,1,ux,all
d,7377,uy,0
d,7377,uz,0
nsel,s,loc,x,76/1000,76/1000
nsel,r,loc,y
nsel,r,loc,z
f,all,fx,pload1
allsel
nsel,s,loc,x,-76/1000,-76/1000
cp,1,ux,all
nsel,s,loc,x,-76/1000,-76/1000
nsel,r,loc,y
nsel,r,loc,z
f,all,fx,pload2
allsel
solve
finish
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

thingisshit

新虫 (小有名气)

引用回帖:
21楼: Originally posted by 11050209 at 2017-03-08 16:08:18
拉伸过程中端面的节点位移本来就是一致的,另外你的意思是顶帽和板之间的接触设置有问题?
...

我的意思是你两个端面的节点都是用的cp,1,ux,all,第二组你应该用cp,2,ux,all,不然的话相当于两个端面的节点位移都耦合到了一起

[ 发自手机版 http://muchong.com/3g ]
22楼2017-03-08 16:28:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

11050209

银虫 (正式写手)

2楼2017-03-06 22:24:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

11050209

银虫 (正式写手)

引用回帖:
2楼: Originally posted by 11050209 at 2017-03-06 22:24:14

应力云图

发自小木虫Android客户端
3楼2017-03-06 22:24:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

11050209

银虫 (正式写手)

4楼2017-03-06 22:24:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

11050209

银虫 (正式写手)

5楼2017-03-07 00:14:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

thingisshit

新虫 (小有名气)

你是按照例子算出来,然后比例子小?
6楼2017-03-07 09:10:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

11050209

银虫 (正式写手)

引用回帖:
6楼: Originally posted by thingisshit at 2017-03-07 09:10:11
你是按照例子算出来,然后比例子小?

自己算的
7楼2017-03-07 11:20:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

thingisshit

新虫 (小有名气)

引用回帖:
7楼: Originally posted by 11050209 at 2017-03-07 11:20:15
自己算的...

应力看起来也不算小啊,contact pressure也有200多呢
8楼2017-03-07 13:22:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

11050209

银虫 (正式写手)

引用回帖:
8楼: Originally posted by thingisshit at 2017-03-07 13:22:42
应力看起来也不算小啊,contact pressure也有200多呢...

要是200MPa还差不多吧,200Pa太小了
9楼2017-03-07 14:32:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

buaa_hw

银虫 (正式写手)

10楼2017-03-07 16:30:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 11050209 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 关于Filecode分析方法 +5 majunge000 2026-08-10 5/250 2026-08-12 18:31 by wangze12014
[基金申请] 不应该看fileCode +5 且听虎啸 2026-08-12 6/300 2026-08-12 16:26 by Tide man
[基金申请] 好奇怪的filecode +5 布布和一二 2026-08-08 6/300 2026-08-12 16:11 by 云上清扬
[基金申请] 应该是93bebmhtak前后十一个字符比较关键 +22 Lanmanbaby 2026-08-09 35/1750 2026-08-12 14:38 by sdfapple719
[论文投稿] 职称评审,求友友推荐见刊最快的期刊 +4 工厂打螺丝 2026-08-08 4/200 2026-08-12 09:18 by hansi2025
[基金申请] 综述论文作为代表作会不会影响评审专家的印象分? +11 yufeiwaner 2026-08-09 13/650 2026-08-12 08:17 by yufeiwaner
[基金申请] 是这周出结果还是下周出结果? +3 yuleib84 2026-08-11 3/150 2026-08-11 21:53 by jnhyjjm
[硕博家园] 读博的好处 +3 lnee 2026-08-11 3/150 2026-08-11 18:10 by 希望我好好的
[基金申请] 分享一下我之前已中青C的计划书的filecode +3 布布和一二 2026-08-11 4/200 2026-08-11 15:20 by aasahr
[基金申请] 听说今天filecode变了 +26 布布和一二 2026-08-06 49/2450 2026-08-11 13:18 by WH3796
[基金申请] 我的国基提前知道中了,可是同事的操作让我实在接受不了,怎么会有这样的人 +8 家与远方 2026-08-10 13/650 2026-08-11 09:49 by 家与远方
[基金申请] 基金中了 +15 laoda193707 2026-08-06 15/750 2026-08-11 00:11 by jiafei2190
[基金申请] 关于代码变化问题,想知道的进来 +17 且听虎啸 2026-08-07 24/1200 2026-08-10 18:35 by zhangduo2008
[基金申请] 这样的filecode谁见过 +11 布布和一二 2026-08-08 22/1100 2026-08-10 11:10 by wmfsnow
[基金申请] 关于filecode,很负责任的告诉大家 +6 爱看书的可乐 2026-08-08 7/350 2026-08-08 22:13 by a_niu
[基金申请] 国基金的申报应该改成非等额制,评价高的钱多评价低的钱少,但是增加资助率 +7 a089 2026-08-07 7/350 2026-08-08 18:05 by gltch
[基金申请] 关于filecode +4 布布和一二 2026-08-07 7/350 2026-08-07 22:55 by zhanghaozhu
[基金申请] 化学口download_prp&fileCode的固定段好像这几天一直没变,有变的大神么? +3 Tide man 2026-08-07 4/200 2026-08-07 22:39 by Tide man
[基金申请] 固定端突然变了,今天 +6 archvillain 2026-08-06 10/500 2026-08-07 16:03 by 医学老男孩
[基金申请] filecode +14 等待解的谜 2026-08-06 19/950 2026-08-07 12:20 by wlwhappy
信息提示
请填处理意见