bolted conections

[Mambuscay]

hi, i have a problem with a simple problem of contact, is a simulation of a test of bolt parallel to grain embedment

! ************** Macro: P.Mac *******************************
! Descripcion :
! Todas las unidades en el SI
!*********************************************************************
/title,"Calculo de parametros"
/prep7
!Element Attributes
ET,1,SOLID45,
KEYOPT,1,1,0 !0 - Include extra displacement shapes
KEYOPT,1,2,0 !0 - Full integration with or without extra displacement shapes, depending on the setting of KEYOPT(1)
KEYOPT,1,4,0 !0 - Element coordinate system is parallel to the global coordinate system
KEYOPT,1,5,2 !2 -- Nodal Stress Solution
KEYOPT,1,6,3 !0 -- Include nonlinear solution at each integration point
keyopt,1,9,0,!0 -- No user subroutine to provide inicial stress (default)
ET,2,170
!***Configuaraciones de elementos de contacto
KEYOPT,2,1,0
KEYOPT,2,2,0
KEYOPT,2,3,0
KEYOPT,2,4,0
KEYOPT,2,5,0
KEYOPT,2,6,000
ET,3,174
KEYOPT,3,1,0 !Selects degrees of freedom 0-- UX,UY,UZ
KEYOPT,3,2,2
KEYOPT,3,3,0
KEYOPT,3,4,0
KEYOPT,3,5,2
KEYOPT,3,6,0
KEYOPT,3,7,0
KEYOPT,3,8,0
KEYOPT,3,9,4
KEYOPT,3,10,2
KEYOPT,3,11,0 !Beam/Shell thickness effect: 0 -- Exclude
KEYOPT,3,12,4

R1 = 3.05/1000 ! target sphere/cylinder/cone radius
R2 = 3.05/1000 ! target cone 2nd radius
FKN = 1 ! * normal penalty stiffness
FTOLN = 0.001 ! * penetration tolerance
ICONT = 0.0 ! * initial contact closure
PINB = 0.001 ! * pinball region
PMAX = 0.0 ! * initial allowable penetration < PMAX
PMIN = 0.0 ! * initial allowable penetration > PMIN
TAUMAX = 1.0E+20 ! maximum friction stress
CNOF = 0.0 ! contact surface offset (+ closes, - opens)
FKOP = 1.0 ! * contact opening stiffness
FKT = 0.0 ! * tangent penalty stiffness
COHE = 0.0 ! contact cohesion

R, 1, R1, R2, FKN, FTOLN, ICONT, PINB
RMORE, PMAX, PMIN, TAUMAX, CNOF, FKOP, FKT

EL=32271000000
ET=361350000
ER=ET

vLT=0.25
vRL=0.05*vLT
vRT=vLT
vLR=vLT
vTL=vRL
vTR=vRT

GLT=((EL*ET)**0.5)/(2*(1+(vTL*vLT)**0.5))
GTL=GLT
GRL=GLT
GRT=(ET)/(2*(1+vRT))
GLR=GRL
GTR=GRT

Tan_EL=27931000000
Tan_ET=252300000
Tan_ER=Tan_ET
Tan_GRL=0.1*GRL
Tan_GLT=0.1*GLT
Tan_GRT=0.1*GRT


EsfFy_R=4400000
EsfFy_L=29830000
EsfFy_T=EsfFy_R

Dey_TL=((EsfFy_T)/(2*(ET-Tan_ET)))*(ET/GTL)**0.5
Dey_RT=((EsfFy_T)/(2*(ET-Tan_ET)))*(ET/GRT)**0.5
Dey_LT=Dey_TL
Dey_RL=Dey_LT
Dey_TR=Dey_RT
Dey_LR=Dey_RL

EsfFy_RT=GRT*Dey_RT
EsfFy_LT=GLT*Dey_LT
EsfFy_RL=GRL*Dey_RL
EsfFy_TR=GTR*Dey_TR
EsfFy_TL=GTL*Dey_TL
EsfFy_LR=GLR*Dey_LR


!Parametros para definir el modelo generalizado de Hill


!*********Definicion de dimensiones
D_dim=3.05/1000 !m CAMBIAR
L_dim=50/1000 !m CAMBIAR
W_dim=40/1000 !m CAMBIAR
T_dim=30/1000 !m CAMBIAR
F_zapl=2 !CAMBIAR
!*******Defincion de tiempos y carga aplicada
Carga=-0.429*9.81*(10**3)
NStep=1

Despl_y=-0.5*D_dim
! ************** Macroefinicion de materiales.Mac *******************************
! Descripcion :
! Todas las unidades en el SI
!*********************************************************************

!**Material 2=el material NO AFECTADO
MP,EY,1,EL
MP,EX,1,ER
MP,EZ,1,ET
MP,GXY,1,GRL
MP,GYZ,1,GLT
MP,GXZ,1,GRT
MP,PRXY,1,vRL
MP,PRYZ,1,vLT
MP,PRXZ,1,vRT
MP,MU,3,0.7
! Define anisotropic material properties
tb,aniso,1
tbdata,1,EsfFy_R,EsfFy_L,EsfFy_T ! Tensile yield stress (x,y & z)
tbdata,4,Tan_ER,Tan_EL,Tan_ET ! Tangent moduli (tensile)
tbdata,7,EsfFy_R,EsfFy_L,EsfFy_T ! Compressive yield stress (x,y & z)
tbdata,10,Tan_ER,Tan_EL,Tan_ET ! Tangent moduli (compressive)
tbdata,13,EsfFy_RL,EsfFy_LT,EsfFy_RT ! Shear yield stress (xy,yz,xz)
tbdata,16,Tan_RL,Tan_LT,Tan_RT ! Tangent moduli (shear)

!**Material 2= material del sujetador
MP,EX,2,2.1*10**11 ! Young's modulus for material
MP,PRXY,2,0.28 !



! ************** MODULO:Geometria y enmallado.Mac *******************************
! Descripcion :
! Todas las unidades en el SI
!*********************************************************************
c=0.5*W_dim-0.5*D_dim*1.1
k,1,0,0,0
k,2,c,0,0
k,3,W_dim-c,0,0
k,4,W_dim,0,0
k,5,0,L_dim,0
k,6,c,L_dim,0
k,7,W_dim-c,L_dim,0
k,8,W_dim,L_dim,0
k,9,c,L_dim-1.1*0.5*D_dim,0
k,10,W_dim-c,L_dim-0.5*D_dim*1.1,0
k,11,W_dim*0.5,L_dim-0.5*D_dim,0
k,12,W_dim*0.5,L_dim-0.5*D_dim*1.1
k,13,0,L_dim-0.5*D_dim*1.1,0
k,14,W_dim,L_dim-0.5*D_dim*1.1,0

L,11,12
A,1,2,9,13
A,2,3,10,9
A,3,4,14,10
A,13,9,6,5
A,9,10,7,6
A,10,14,8,7

CYL4,0.5*W_dim,L_dim,0.5*D_dim,-180
ASBA,5,7
ASBL,8,1
VEXT,1,7,,,,T_dim
vglue,1,2,3,4,5,6,7
CYL4,0.5*W_dim,L_dim,0.5*D_dim,,,,T_dim


k,41,W_dim*0.5,L_dim+D_dim*0.5,0
k,42,W_dim*0.5,L_dim+D_dim*0.5,T_dim
l,41,42
ASBL,39,71

div=20 !PAR*2********
divh=10

TYPE,1
MAT,1
ESIZE,T_dim/divh,
MSHAPE,0,3D
MSHMID,2
EXTOPT,ESIZE,div,1
ACCAT,33,16
ACCAT,17,23
Lesize,34,,,div/4
Lesize,43,,,div/4
Lesize,55,,,div/4
Lesize,33,,,div/4
Lesize,12,,,div/4
Lesize,16,,,div/4
Lesize,15,,,div/4
Lesize,21,,,div/4
MSHKEY,1
VSEL,S,VOLU,,7,7
VMESH,7,7,1
VSEL,S,VOLU,,5,5
VMESH,5,5,1

VSEL,S,VOLU,,6,6
VMESH,6,6,1

VSEL,S,VOLU,,4,4
VMESH,4,4,1

VSEL,S,VOLU,,2,2
EXTOPT,ESIZE,10,1
VSWEEP,2,


VSEL,S,VOLU,,1,1
VMESH,1,1,1

VSEL,S,VOLU,,3,3
VMESH,3,3,1


MAT,2
ESIZE,3.1415926535897932384*D_dim/(div*2),
VSEL,S,VOLU,,8,8
EXTOPT,ESIZE,10,1
VSWEEP,8,
ALLSEL,ALL

TYPE,3
REAL,1
MAT,1
ASEL,S,AREA,,27
ASEL,A,AREA,,36
NSLA,S,1
ESURF,ALL

TYPE,2
REAL,1
MAT,1
ASEL,S,AREA,,40
NSLA,S,1
ESURF,ALL
FINISH

/SOLU
ASEL,S,AREA,,9,9,1
NSLA,S,1
D,ALL,ALL,0
ASEL,S,AREA,,14,14,1
NSLA,S,1
D,ALL,ALL,0
ASEL,S,AREA,,19,19,1
NSLA,S,1
D,ALL,all,0
ALLSEL,ALL
NROP,UNSYM
SOLCONTROL,ON
KBC,0
TIME,1
NROPT,FULL,,OFF
NSUBST,NStep,NStep,NStep
NEQIT,1000
OUTRES,BASIC,ALL,
PRED,OFF
lsel,s,line,,72,72,1
NSLL,s,71
cm,mid,node !CM, Cname, Entity Groups geometry items into a component.

*GET,nodes_num,NODE,0,count
F,all,FY,Carga/nodes_num
ALLSEL,ALL
FINISH

 

[msquared48]

And the problem is????????????

Mike McCann
MMC Engineering

 

[hokie66]

New members get some leeway, but that is ridiculous. Suggest you red-flag your post and try again, with a simple statement of your problem.