/batch
!2006/10/26 作成
!!!!!!!!!!!!!軸力発生のための不整量の与え方を変更!!!!!!!!!!!!!!!!!!
!!!!!!!!!!!!!!!!!ペナルティー剛性は0.1!!!!!!!!!!!!!!!!



!リベットの作成開始!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

!リベット寸法　
d = 10         !軸径
L1 = 0		!軸上部
L2 = -18	!軸下部


!モデリング開始
/prep7

et,1,186       !二次六面体要素
MP,EX,1,205000
MP,PRXY,1,0.3

!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!軸部を作成

*afun,deg

wpro,,-90,
cylind,0,d/2,L1,L2,0,90
cylind,0,d/2,L1,L2,90,180
cylind,0,d/2,L1,L2,180,270
cylind,0,d/2,L1,L2,270,360

!頭部内部
cylind,0,d/2,L2,L2-7,0,90
cylind,0,d/2,L2,L2-7,90,180
cylind,0,d/2,L2,L2-7,180,270
cylind,0,d/2,L2,L2-7,270,360

cylind,0,d/2,L1,L1+7.5,0,90
cylind,0,d/2,L1,L1+7.5,90,180
cylind,0,d/2,L1,L1+7.5,180,270
cylind,0,d/2,L1,L1+7.5,270,360
wpro,,90,


vsel,none
!ボルト内部の四角柱
block,-2,2,L1+7.5,L2-7,-2,2
cm,tmp,volu
allsel
csys,5
vgen,1,tmp,,,,45,,,,1
csys,0

vsel,s,loc,y,L1+7.5,L2-7	!演算関係のボリューム選択
cmsel,u,tmp
cm,tmp2,volu
vsel,all
vsbv,tmp2,tmp,sepo,,keep
cmdele,tmp,volu
numcmp,all


wpro,,-90,
!下頭部外側
cylind,9,d/2,L2,L2-7,0,90
cylind,9,d/2,L2,L2-7,90,180
cylind,9,d/2,L2,L2-7,180,270
cylind,9,d/2,L2,L2-7,270,360

!上頭部外側
cylind,8.5,d/2,L1,L1+7.5,0,90
cylind,8.5,d/2,L1,L1+7.5,90,180
cylind,8.5,d/2,L1,L1+7.5,180,270
cylind,8.5,d/2,L1,L1+7.5,270,360

*afun,rad
wpro,,90,

btol,0.01
vglue,all
btol,defa
numcmp,all


cm,rivet,volu !リベットをコンポーネント化

!モデリング終了・メッシング開始!!!!!!!!!!!!!!!!!!!!!!!

!下側頭部について!!!!!!!!!!!!!!!!!!!
!頭上部
csys,5
lsel,s,loc,x,9
csys,0
lsel,r,loc,y,L2
LESIZE,all, , ,4, , , , ,1
!頭厚さ方向
csys,5
lsel,s,loc,x,9
csys,0
lsel,r,loc,y,L2-7+1,L2-1
LESIZE,all, , ,2, , , , ,1
!大きい円径方向
csys,5
lsel,s,loc,x,6,8
csys,0
lsel,r,loc,y,L2
LESIZE,all, , ,2, , , , ,1
!小さい円径方向
lsel,s,loc,y,L2-7
csys,5
lsel,r,loc,x,3,4.5
csys,0
LESIZE,all, , ,2, , , , ,1
allsel
!頭部外周メッシュ
type,1
smrt,off
mshkey,1
mshape,0,3D
csys,5
vsel,s,loc,x,5,10
csys,0
vsel,r,loc,y,L2,L2-7
VMESH,all

!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!上側頭部について
!頭下部
csys,5
lsel,s,loc,x,8.5
csys,0
lsel,r,loc,y,L1
LESIZE,all, , ,4, , , , ,1
!頭厚さ方向
csys,5
lsel,s,loc,x,8.5
csys,0
lsel,r,loc,y,L1+1,L1+7.5-1
LESIZE,all, , ,2, , , , ,1
!大きい円径方向
csys,5
lsel,s,loc,x,6,8
csys,0
lsel,r,loc,y,L1
LESIZE,all, , ,2, , , , ,1
!小さい円径方向
lsel,s,loc,y,L1
csys,5
lsel,r,loc,x,3,4.5
csys,0
LESIZE,all, , ,2, , , , ,1
allsel
!頭部外周メッシュ
type,1
smrt,off
mshkey,1
mshape,0,3D
csys,5
vsel,s,loc,x,5,10
csys,0
vsel,r,loc,y,L1,L1+7.5
VMESH,all



!小さい円柱中はなしの選択
type,1
csys,5
vsel,s,loc,x,2,5
csys,0
vsel,r,loc,y,L2-7+1,L2-1
vsweep,all
type,1
csys,5
vsel,s,loc,x,2,5
csys,0
vsel,r,loc,y,L1+1,L1+7.5-1
vsweep,all


!軸部のサイズ
esize,,5
type,1
csys,5
vsel,s,loc,x,2,5
csys,0
vsel,r,loc,y,L2,L1
vsweep,all

!中心の直方体メッシュ
type,1
csys,5
vsel,s,loc,x,0,2
csys,0
vsel,r,loc,y,L2,L1
vsweep,all	
allsel
!メッシング終了

!!!!!!!!!!!!!コンポーネント化!!!!!!

!上頭部下面
asel,s,loc,y,L1
csys,5
asel,r,loc,x,d/2,8.5
csys,0
cm,up_bearing,area
allsel

!下頭部上面
asel,s,loc,y,L2
csys,5
asel,r,loc,x,d/2,9
csys,0
cm,down_bearing,area
allsel
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!


parsav
allsel
cdopt,anf
cdwrite,comb,rivet
finish
/clear,start
!リベットの作成終了!!!!!!!!!!!!!!!!!!!!!!!!!!!!!



!被締結物の作成開始!!!!!!!!!!!!!!!!!!!!!!!!!!!
!plate1（上）の作成
!モデリング開始
parres,new
/prep7
et,1,186
thick=9      !板厚

k,1000,0,0,0
k,1001,-40/2,0,40/2
k,1002,    0,0,40/2
a,1000,1001,1002

k,1003,0,0,0
k,1004,    0,0,40/2
k,1005, 40/2,0,40/2
a,1003,1004,1005

k,1006,0,0,0
k,1007,-40/2,0,-40/2
k,1008,0,0,-40/2
a,1006,1007,1008

k,1009,0,0,0
k,1010,0,0,-40/2
k,1011, 40/2,0,-40/2
a,1009,1010,1011

k,1012,0,0,0
k,1013,-40/2,0,-40/2
k,1014,-40/2,0,0
a,1012,1013,1014

k,1015,0,0,0
k,1016,-40/2,0,0
k,1017,-40/2,0,40/2
a,1015,1016,1017

k,1018,0,0,0
k,1019, 40/2,0,-40/2
k,1020, 40/2,0,0
a,1018,1019,1020

k,1021,0,0,0
k,1022, 40/2,0,0
k,1023, 40/2,0,40/2
a,1021,1022,1023

csys,5
asel,s,loc,x,12,20
csys,0
asel,r,loc,y,0 
vext,all,,,0,-thick,0
cm,tmp,volu
vsel,none

wprota,,-90,
cylind,0,12,-thick,0,0,360
wprota,,90,
cm,tmp2,volu

vsel,all
vsbv,tmp,tmp2,sepo,,dele
cmdele,tmp
cmdele,tmp2
numcmp,all

wpro,,-90,
cylind,5.5,12,-thick,0,-45,45
cylind,5.5,12,-thick,0,45,135
cylind,5.5,12,-thick,0,135,225
cylind,5.5,12,-thick,0,225,315
wpro,,90

!拘束用端部の作成
k,1024,-20,0,0
k,1025,-20,0,20
k,1026,-25,0,20
k,1027,-25,0,0
a,1024,1025,1026,1027
k,1028,-20,0,-20
k,1029,-20,0,0
k,1030,-25,0,0
k,1031,-25,0,-20
a,1028,1029,1030,1031
asel,s,loc,x,-21,-24
vext,all,,,0,-thick,0
allsel

BLOCK,-25,-125,0,-9,-20,0,
BLOCK,-25,-125,0,-9,20,0,

vglue,all
numcmp,all

!モデリング終了・メッシング開始!!!!!!!!!!!!
!!!!!!!!!!!!!!!!!!メッシュサイズ定義
!外周
lsel,s,loc,x,20
lsel,a,loc,x,-20
lsel,a,loc,x,-25
lsel,a,loc,z,20
lsel,a,loc,z,-20
lsel,u,loc,y,-1,-thick+1
lsel,u,loc,x,-125+1,-25-1
LESIZE,all, , ,3, , , , ,1

lsel,s,loc,x,-125+1,-25-1
LESIZE,all, , ,20, , , , ,1


!筒の対角線
lsel,s,loc,y,0
lsel,a,loc,y,-thick
csys,5
lsel,r,loc,x,7,10
csys,0
LESIZE,all, , ,5, , , , ,1

!四角の対角線など
lsel,s,loc,y,0
lsel,a,loc,y,-thick
csys,5
lsel,r,loc,x,13,22
csys,0
LESIZE,all, , ,2, , , , ,1

!板厚方向
lsel,s,loc,x,20
lsel,a,loc,x,-20
lsel,a,loc,x,-25
lsel,a,loc,x,-125
lsel,a,loc,z,20
lsel,a,loc,z,-20
lsel,r,loc,y,-1,-thick+1
LESIZE,all, , ,3, , , , ,1

!筒の周方向
lsel,s,loc,y,0
lsel,a,loc,y,-thick
csys,5
lsel,r,loc,x,5.5
csys,0
LESIZE,all, , ,6, , , , ,1
allsel


type,1        !二次要素
vsweep,all
numcmp,all
!メッシング終了

!!!!!!!!!!!!コンポーネント化!!!!!!!!!!!
cm,plate1,volu

!plate1上面
ASEL,S,LOC,Y,0
CSYS,5  
ASEL,R,LOC,X,6,12
CSYS,0
cm,plate1_bearing,area
allsel

!plate1下面
CMSEL,S,PLATE1
ALLSEL,BELOW,VOLU 
ASEL,R,LOC,Y,-thick
ASEL,u,LOC,x,-125+1,-25-1
cm,plate1_slip,area
allsel

!plate1作成終了・plate2（下）作成開始!!!!!!!!!!!!
VGEN,2,all, , ,0,-9, , ,0  
VSEL,s,LOC,Y,-9,-18 
VSYMM,X,all, , , ,0,1  

cm,plate2,volu

!plate2上面
CMSEL,S,PLATE2
ALLSEL,BELOW,VOLU
ASEL,R,LOC,Y,-thick
ASEL,u,LOC,x,25+1,125-1
cm,plate2_slip,area
allsel

!plate2下面
ASEL,S,LOC,Y,-2*thick
CSYS,5  
ASEL,R,LOC,X,6,12
CSYS,0
cm,plate2_bearing,area
allsel


parsav
allsel
cdopt,anf
cdwrite,comb,plate
finish
/clear,start

!被締結物の作成終了!!!!!!!!!!!!!!!!!!!!!!!!!!!!





parres,new
/prep7
cdread,comb,rivet
cdread,comb,plate
allsel



!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!被締結物の摩擦係数0.20
!リベットの摩擦係数0.20
!接触要素の設定
et,10,targe170
et,11,conta174
keyopt,11,2,1	!ペナルティ関数
!keyopt,11,9,2	!形状的初期食い込み/ギャップとオフセットの両方をランプ荷重として含める
mp,mu,2,0.20      !リベット
mp,mu,3,0.20      !被締結物

!リベット（上）-被締結物間
R,10,,,0.1,,,,
!接触定義リアルコンスタント番号10
cmsel,s,plate1_bearing
nsla,s,1
real,10
type,10          !targe170
mat,3
esurf
allsel
cmsel,s,up_bearing
nsla,s,1
type,11          !conta174
mat,2
esurf
allsel,all

!被締結物同士間
R,11,,,0.1,,,,
!接触定義リアルコンスタント番号11
cmsel,s,plate1_slip
nsla,s,1
real,11
type,10          !targe170
mat,3
esurf
allsel
cmsel,s,plate2_slip
nsla,s,1
type,11          !conta174
mat,3
esurf
allsel,all

!リベット（下）-被締結物間
R,12,,,0.1,,,,
!接触定義リアルコンスタント番号12
cmsel,s,plate2_bearing
nsla,s,1
real,12
type,10          !targe170
mat,3
esurf
allsel
cmsel,s,down_bearing
nsla,s,1
type,11          !conta174
mat,2
esurf
allsel,all


!境界条件!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!被締結物端部をx方向にカップリング拘束
CMSEL,S,PLATE1
ALLSEL,BELOW,VOLU 
asel,r,loc,x,-125
nsla,s,1
cp,1,ux,all
allsel

CMSEL,S,PLATE2
ALLSEL,BELOW,VOLU 
asel,r,loc,x,125
nsla,s,1
cp,2,ux,all
allsel


!被締結物端部をy,z方向に拘束
CMSEL,S,PLATE1
ALLSEL,BELOW,VOLU
asel,r,loc,x,-125
nsla,s,1
d,all,uy
d,all,uz
allsel
CMSEL,S,PLATE2
ALLSEL,BELOW,VOLU
asel,r,loc,x,125
nsla,s,1
d,all,uy
d,all,uz
allsel



!
!プリテンション要素の定義
et,13,prets179
R,13, ,8000, ,  
type,13
real,13
CMSEL,S,RIVET  
ALLSEL,BELOW,VOLU 
PSMESH,1,ziku,25000,ALL, ,   0,Y,-10, , ,13,prete,pretn 

!プリテンション荷重(軸力8000N)の設定
SLOAD,1,9,LOCK,FORC,8000, 2,3 
!SLOAD,ALL,9,,DISP,0.025446, 1
allsel,all



finish



!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!解析の実行!!!!!!!!!!!!!!!!!!!!!!

/solu
antype,static


!!!!!!!!!!!!!!!
!plate1の代表接点
d,910,ux,0
!!!!!!!!!!!!!!!
!plate2の代表接点
d,6260,ux,0
!!!!!!!!!!!!!!!


!nldiag,nrre,1
!nldiag,maxf,10

solcontrol,on
nsub,20,1000,1
nlgeom,on
nropt,full,,off
outres,all,all
solve
save




antype,,restart,
solve

save

!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

antype,,restart,

!plate1の代表接点
d,910,ux,0.2
!plate2の代表接点
d,6260,ux,0.0

solcontrol,on
autots,on
nsub,35,1000,30
nlgeom,on
nropt,full,,off
outres,all,all
solve

save

FINISH
  
/exit,all