|
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
|