OkapiEng
Mechanical
- Jul 12, 2007
- 2
Hi everybody
My name is Okapi and i'm new on this forum.
My english is not perfect, because i'm not native speaker...
I have a probleme to analyse a not penetration surface
Please, look the programme and modify the programme to have not penetration surface (i looked with TARGE170 et CONTA but i'dont do that)
thank you
!*************************************
! Penetration problem
!*************************************
fini
/clear
/title,toto
/gra,full
/vie,,-1,-1,1
/vup,,z
/type,,4
/prep7
/num,1
tnum,1
btol,1e-5
epsi=1e-4
!Définition du fichier de plots de sortie
!----------------------------------------
filplot='toto'
taille=0.005 !taille des elements
!Définitition des éléments
!-------------------------
ET,1,63 !plaques mécaniques structure
r,1,0.8e-3
!========================
!Définition des matériaux
!========================
!Mat1
!---------
dens,1,3000
ex,1,20000e6
nuxy,1,0.33
!Mat2
!---------
dens,2,1000
ex,2,70000e6
nuxy,2,0.33
!==============================================
! Géométrie
!==============================================
allsel
csys,0 $ wpcsys,-1
rect,0e-3,20e-3,0,200e-3
csys,0 $ cloc,11,0,0,0,2e-3,,, $ wpcsys,-1
rect,0e-3,20e-3,0,200e-3
csys,0 $ cloc,11,0,0,0,0,,90, $ wpcsys,-1
rect,0e-3,20e-3,0,2e-3
csys,0 $ cloc,11,0,0,200e-3,0,,90, $ wpcsys,-1
rect,0e-3,20e-3,0,2e-3
aplot
allsel
nummrg,kp,1e-6
!-----------
!Mesh
!-----------
csys,0 $ wpcsys,-1
asel,s,loc,z,0
esize,2e-3
type,1
mat,1
real,1
amesh,all
csys,0 $ wpcsys,-1
asel,s,loc,z,2e-3
esize,2e-3
type,1
mat,2
real,1
amesh,all
csys,0 $ wpcsys,-1
asel,s,loc,y,0
esize,2e-3
type,1
mat,2
real,1
amesh,all
csys,0 $ wpcsys,-1
asel,s,loc,y,200e-3
esize,2e-3
type,1
mat,2
real,1
amesh,all
nummrg,kp,1e-6
*get,NUM,node,,num,max
!==============
!analyse modale
!==============
csys,0 $ wpcsys,-1
esel,s,mat,,1
nsle
nsel,r,loc,y,0
d,all,all,0
csys,0 $ wpcsys,-1
esel,s,mat,,1
nsle
nsel,r,loc,y,200e-3
d,all,all,0
!==============
!Calcul de la solution
!==============
/solu
allsel
antyp,modal
modopt,lanb,10,0,2000
mxpand,10,0,2000,no
outpr,,
/out,outtoto,lst
solve
/out
/post1
set,1,3
/view,,-1,-1,1
/title,MODAL ANALYSIS - Travail
plns,u,sum
My name is Okapi and i'm new on this forum.
My english is not perfect, because i'm not native speaker...
I have a probleme to analyse a not penetration surface
Please, look the programme and modify the programme to have not penetration surface (i looked with TARGE170 et CONTA but i'dont do that)
thank you
!*************************************
! Penetration problem
!*************************************
fini
/clear
/title,toto
/gra,full
/vie,,-1,-1,1
/vup,,z
/type,,4
/prep7
/num,1
tnum,1
btol,1e-5
epsi=1e-4
!Définition du fichier de plots de sortie
!----------------------------------------
filplot='toto'
taille=0.005 !taille des elements
!Définitition des éléments
!-------------------------
ET,1,63 !plaques mécaniques structure
r,1,0.8e-3
!========================
!Définition des matériaux
!========================
!Mat1
!---------
dens,1,3000
ex,1,20000e6
nuxy,1,0.33
!Mat2
!---------
dens,2,1000
ex,2,70000e6
nuxy,2,0.33
!==============================================
! Géométrie
!==============================================
allsel
csys,0 $ wpcsys,-1
rect,0e-3,20e-3,0,200e-3
csys,0 $ cloc,11,0,0,0,2e-3,,, $ wpcsys,-1
rect,0e-3,20e-3,0,200e-3
csys,0 $ cloc,11,0,0,0,0,,90, $ wpcsys,-1
rect,0e-3,20e-3,0,2e-3
csys,0 $ cloc,11,0,0,200e-3,0,,90, $ wpcsys,-1
rect,0e-3,20e-3,0,2e-3
aplot
allsel
nummrg,kp,1e-6
!-----------
!Mesh
!-----------
csys,0 $ wpcsys,-1
asel,s,loc,z,0
esize,2e-3
type,1
mat,1
real,1
amesh,all
csys,0 $ wpcsys,-1
asel,s,loc,z,2e-3
esize,2e-3
type,1
mat,2
real,1
amesh,all
csys,0 $ wpcsys,-1
asel,s,loc,y,0
esize,2e-3
type,1
mat,2
real,1
amesh,all
csys,0 $ wpcsys,-1
asel,s,loc,y,200e-3
esize,2e-3
type,1
mat,2
real,1
amesh,all
nummrg,kp,1e-6
*get,NUM,node,,num,max
!==============
!analyse modale
!==============
csys,0 $ wpcsys,-1
esel,s,mat,,1
nsle
nsel,r,loc,y,0
d,all,all,0
csys,0 $ wpcsys,-1
esel,s,mat,,1
nsle
nsel,r,loc,y,200e-3
d,all,all,0
!==============
!Calcul de la solution
!==============
/solu
allsel
antyp,modal
modopt,lanb,10,0,2000
mxpand,10,0,2000,no
outpr,,
/out,outtoto,lst
solve
/out
/post1
set,1,3
/view,,-1,-1,1
/title,MODAL ANALYSIS - Travail
plns,u,sum