|
fini
/cle
C************************************************************
C*** GRAPHICS SETTINGS
C************************************************************
/vie,1,-1,-2,-1
/vup,1,z
/pnu,mat,1
/num,1
/esha,1
/dsc,1,1
C************************************************************
C*** PARAMETERS
C************************************************************
a=10 ! EDGE LENGTH, mm
t=0.1 ! PLATE THICKNESS, mm
r1=4 ! RADIUS OF SPHERICAL PUNCH
r2=5 ! RADIUS OF SPHERICAL DIE
E=2e5 ! ELASTIC MODULUS, MPa
rho=7800e-9 ! DENSITY (kg/mm^3)
sy=200 ! YIELD STRESS, MPa
e2=2e3 ! TANGENT MODULUD AFTER YIELD (MPa)
n_int=5 ! # OF INTEGRATION PTS THRU SHELL THICKNESS (5 MAX)
esz=a/10 ! ELEMENT SIZE
dz=-3.5 ! DISTANCE TRAVELED BY PUNCH
t_umax=0.0005 ! TIME WHEN PUNCH REACHES MAX EXCURSION
gap_0=t/2 ! INITIAL GAP (mm)
/fil,explct
/prep7
C************************************************************
C*** MODEL WORKPEICE (PLATE)
C************************************************************
rect,0,a,0,a ! PLATE
cm,plate_a,area
aatt,1,1,1
C************************************************************
C*** MODEL DIE INTO WHICH PLATE WILL BE PUSHED
C************************************************************
wpof,0,0,-t ! FORM
cswp,11,1
cswp,12,2
sphere,r2,,0,90
vsbw,all
csys
vsel,s,loc,z,-t,r2
vdel,all,,,1
alls
vdel,all
asel,s,loc,x
asel,a,loc,y
asel,a,loc,z,-t
adel,all,,,1
csys,11
k,,2*a
k,,2*a,90
a,kp(r2,0,0),kp(2*a,0,0),kp(2*a,90,0),kp(r2,90,0)
alls
cmse,u,plate_a
cm,die_a,area
aatt,2,2,2
C************************************************************
C*** MODEL PUNCH
C************************************************************
wpcs,-1,0
wpof,0,0,r1+t+gap_0
cswp,21,1
cswp,22,2
sphere,r1,,0,90
vsbw,all
csys,21
vsel,s,loc,z,0,r1
vdel,all,,,1
alls
vdel,all
csys
asel,s,loc,x
asel,a,loc,y
csys,21
asel,a,loc,z,0
adel,all,,,1
csys,21
k,,2*a
k,,2*a,90
a,kp(r1,0,0),kp(2*a,0,0),kp(2*a,90,0),kp(r1,90,0)
alls
cmse,u,plate_a
cmse,u,die_a
cm,punch_a,area
aatt,3,3,3
C************************************************************
C*** ATTRIBUTES
C************************************************************
et,1,163,10
et,2,163,1
et,3,163,1
mp,ex,1,E
mp,dens,1,rho
mp,nuxy,1,0.3
tb,biso,1
tbdata,1,sy
tbdata,2,e2
mp,ex,2,E
mp,dens,2,rho
mp,nuxy,2,0.3
mp,ex,3,E
mp,dens,3,rho
mp,nuxy,3,0.3
edmp,rigid,2,7,7
edmp,rigid,3,4,7
r,1,,n_int,t
r,2,,,t
r,3,,,t
C************************************************************
C*** MESH
C************************************************************
esiz,esz
!mshk,1
alls
ames,all
csys,1
nsel,s,loc,x,0,1.5*r2
esln
erefine,all,,,3
C************************************************************
C*** BOUNDARY CONDITIONS
C************************************************************
csys
lsel,s,loc,x
dl,all,,ux
dl,all,,roty
dl,all,,rotz
lsel,s,loc,y
dl,all,,uy
dl,all,,rotx
dl,all,,rotz
C************************************************************
C*** CONTACT
C************************************************************
alls
edcg,ag
C************************************************************
C*** MOVE PUNCH DOWN AND UP
C************************************************************
alls
edpart,create
cbc_rmp5,0,t_umax,dz,100,1,''tm'',''u_punch'',2.5*t_umax,0
edload,add,rbuz,,3,''tm'',''u_punch''
fini
C************************************************************
C*** SOLVE
C************************************************************
/solu
time,2.5*t_umax
solve
fini
C************************************************************
C*** PLOT DISPLACEMENT LAST RESULTS SET
C************************************************************
/post1
set,last
pldi
C************************************************************
C*** GET PLATE DISPLACEMENT FOR CONSTRAINTS TO BE IMPOSED
C*** FOR IMPLICIT RUN
C************************************************************
csys
asel,s,mat,,1
alls,belo,area
nd=node(a,a,0)
uz_nd=uz(nd)
fini
C************************************************************
C*** SAVE EXPLICIT DATABASE
C************************************************************
save
C************************************************************
C*** SWITCH TO IMPLICIT FOR SPRINGBACK: CHANGE JOBNAME
C************************************************************
/fil,implct
C************************************************************
C*** MODIFY ATTRIBUTES
C************************************************************
/prep7
et,1,0 $et,2,0 $et,3,0
et,1,181
r,1,t
tbde,biso,1
C************************************************************
C*** DELETE DISPLACEMENTS ON "NULLED" ELEMENTS COMPRISING
C*** PUNCH AND DIE
C************************************************************
asel,s,mat,,2,3,
alls,belo,area
dlde,all,all
C************************************************************
C*** UPDATE GEOMETRY
C************************************************************
alls
shpp,off
upgeom,1,,,explct,rst
C************************************************************
C*** PREVENT PLATE RIGID BODY MOTION
C************************************************************
alls
d,nd,uz,uz_nd
fini
C************************************************************
C*** IMPORT STRESSES FROM FINAL LS DYNA RESULT
C************************************************************
/solu
rimport,dyna,stress,elem,,,explct,rst
C************************************************************
C*** OBTAIN EXPLICIT SPRINGBACK SOLUTION
C************************************************************
nlge,on
ematwrite,yes ! FLAG TO WRITE EMAT FILE
pstress,on ! FLAG TO WRITE PRESTRESS MATRIX
save
solve
fini
C************************************************************
C*** PLOT NORMAL COMPONENT OF DISPLACEMENT
C************************************************************
/post1
set,last
plns,u,z
fini
C************************************************************
C*** PERFORM MODAL ANALYSIS OF FORMED STRUCTURE
C************************************************************
/fil,modal
/sys,copy implct.emat modal.emat
/sys,copy implct.esav modal.esav
/solu
lscle,lsopt
anty,modal ! PERFORM MODAL ANALYSIS
upcoord,1,on ! ADD PREVIOUS DISPLACEMENTS TO EXISTING COORDS
modo,lanb,5 ! USE BLOCK LANCHOS TO EXTRACT 5 MODES
mxpand,5,,,yes ! EXPAND 5 MODES, INCLUDE ELEM RESULTS
ddel,nd,uz ! DELETE CONSTRAINT ON CORNER NODE
psolve,eiglanb ! CALCULATE EIGEN VALUES AND EIGEN VECTORS
fini
/dsc
save
/solu
expass,on ! EXPAND EIGEN VECTOR SOLUTION
psolve,eigexp
fini
/eof
fini
/cle
C************************************************************
C*** GRAPHICS SETTINGS
C************************************************************
/vie,1,-1,-2,-1
/vup,1,z
/pnu,mat,1
/num,1
/esha,1
/dsc,1,1
C************************************************************
C*** PARAMETERS
C************************************************************
a=10 ! EDGE LENGTH, mm
t=0.1 ! PLATE THICKNESS, mm
r1=4 ! RADIUS OF SPHERICAL PUNCH
r2=5 ! RADIUS OF SPHERICAL DIE
E=2e5 ! ELASTIC MODULUS, MPa
rho=7800e-9 ! DENSITY (kg/mm^3)
sy=200 ! YIELD STRESS, MPa
e2=2e3 ! TANGENT MODULUD AFTER YIELD (MPa)
n_int=5 ! # OF INTEGRATION PTS THRU SHELL THICKNESS (5 MAX)
esz=a/10 ! ELEMENT SIZE
dz=-3.5 ! DISTANCE TRAVELED BY PUNCH
t_umax=0.0005 ! TIME WHEN PUNCH REACHES MAX EXCURSION
gap_0=t/2 ! INITIAL GAP (mm)
/fil,explct
/prep7
C************************************************************
C*** MODEL WORKPEICE (PLATE)
C************************************************************
rect,0,a,0,a ! PLATE
cm,plate_a,area
aatt,1,1,1
C************************************************************
C*** MODEL DIE INTO WHICH PLATE WILL BE PUSHED
C************************************************************
wpof,0,0,-t ! FORM
cswp,11,1
cswp,12,2
sphere,r2,,0,90
vsbw,all
csys
vsel,s,loc,z,-t,r2
vdel,all,,,1
alls
vdel,all
asel,s,loc,x
asel,a,loc,y
asel,a,loc,z,-t
adel,all,,,1
csys,11
k,,2*a
k,,2*a,90
a,kp(r2,0,0),kp(2*a,0,0),kp(2*a,90,0),kp(r2,90,0)
alls
cmse,u,plate_a
cm,die_a,area
aatt,2,2,2
C************************************************************
C*** MODEL PUNCH
C************************************************************
wpcs,-1,0
wpof,0,0,r1+t+gap_0
cswp,21,1
cswp,22,2
sphere,r1,,0,90
vsbw,all
csys,21
vsel,s,loc,z,0,r1
vdel,all,,,1
alls
vdel,all
csys
asel,s,loc,x
asel,a,loc,y
csys,21
asel,a,loc,z,0
adel,all,,,1
csys,21
k,,2*a
k,,2*a,90
|