**
***  file ellipse.inp
*HEADING
Elliptic hole in an elastic plate
full plate section 20x20 mm2, one quarter: 10x10 mm2
*SYSTEM
***** units: mm, N, kg, s, MPa
***** quadratic elements with 2 planes of symmetry
*****
*node,nset=n1-21
*node,nset=n21
*node,nset=n41
****** nodes 21 to 41 on line y = 0.
*nfill,nset=n21-41,bias=0.75
n21,n41,20,1
***** nodes 2001 to 2041 on line y = 10.
*node
*ngen,nset=nn1
*******************
*nset,nset=nn0
n1-21,n21-41
**** all nodes in xy plane
*nfill,nset=nall,bias=0.75
nn0,nn1,20,100
****
******************************
*element,type=cps8
**** quadratic element with 8 nodes
*elgen,elset=eall
****************************
*elset,elset=etop,generate
*Surface, type=ELEMENT, name=top
etop,S3
***********  for BC's
*nset,nset=nx,generate
?????????
*nset,nset=ny,generate
?????????
************  for results:
*nset,nset=nmes
*elset,elset=emes
*******
*SOLID SECTION, ELSET=eall, MATERIAL=alu
****
*MATERIAL,NAME=alu
*density
***** en kg/mm3
2.7e-6
*elastic
75.e3,0.3
***************************************************
*STEP,nlgeom=no
Stresses in plate with an elliptical hole: loading 1 MPa
******
*static
****dt0, step time, min, max
0.01,1.
**** BC's
*boundary
????
????
**** pressure
*dsload
top,p,-1
*** surface,trcvec,amplitude,x-coord,y-coord
*** top,trvec,1,0.,1.
****************************
*el print,position=nodes,elset=emes,freq=100
mises,s11,s22
*node print,nset=nn0,freq=100
coord,u1,u2
************************************************
*OUTPUT,FIELD,freq=1,var=preselect
*output,history,freq=1,var=preselect
************************************
*monitor,node=1,dof=2
*END STEP
*************