C*** PRE-PROCESSAMENTO ===========================================
/PREP7
/TITLE, EXEMPLO - PLACA


C*** DEFINICAO DAS CONSTANTES ====================================

PI=3.141593
g=9.81				! aceler da gravidade
E=200e9			! modulo de elasticidade do aco
densAco=7.8e3			! densidade do aco

t12 = 2.78e-3
t14 = 1.98e-3


C*** DEFINICAO DAS CONSTANTES GEOMETRICAS ========================

Lx = 1
Lx2 = .1
Ly = 1.


C*** DEFINICAO DOS ELEMENTOS =========================================

ET,1,SHELL63,,,

ET,2,BEAM4,,,



C*** DEFINICAO DAS PROPRIEDADES ======================================

MP,EX,1,E				! Material da Estrutura
MP,NUXY,1,0.3
MP,DENS,1,densAco

R,1,t12				! espessura da placa



C*** DEFINICAO DA GEOMETRIA ===========================================

x1=0
y1=0
x2=0
y2=Ly
x3=Lx+Lx2
y3=Ly
x4=Lx+Lx2
y4=Ly-Lx2
x5=Lx
y5=y4
x6=Lx
y6=y1


C*** DEFINICAO DOS KEYPOINTS

k,1,x1,y1
k,2,x2,y2
k,3,x3,y3
k,4,x4,y4
k,5,x5,y5
k,6,x6,y6
!k,7,x7,y7
!k,8,x8,y8


C*** GERACAO DAS AREAS =====================

A,1,2,3,4,5,6

C*** GERACAO DA MALHA =====================

TYPE,1
MAT,1
REAL,1

!ESIZE,ALL,,5
!ESIZE,,5
SMRTSIZE,6

AMESH,all


C***  CONDICOES DE CONTORNO E FORCAMENTO ============================

DK,1,ALL,0
DK,2,ALL,0
DK,3,ALL,0

NSEL,S,LOC,X,(Lx/2-0.1),(Lx/2+0.1)		! Condicao de Contorno de Engaste
NSEL,R,LOC,y,(Ly/2-0.1),(Ly/2+0.1)		
F,ALL,FZ,-400000
NSEL,ALL


C*** SOLUCAO ===========================================
!/SOLUTION

!SOLVE
!FINISH
!/POST1
!SET,1,1

!ESEL,S,TYPE,,1
!ETABLE,SeqvtopI,NMISC,5
!ETABLE,SeqvtopJ,NMISC,10
!ETABLE,SeqvtopK,NMISC,15
!ETABLE,SeqvtopL,NMISC,20

!ETABLE,SeqvbotI,NMISC,25
!ETABLE,SeqvbotJ,NMISC,30
!ETABLE,SeqvbotK,NMISC,35
!ETABLE,SeqvbotL,NMISC,40
!ESEL,ALL