Shell idealization ACEL Fillet symmetry flat end hemispherical end-Ass #3

jinqinghai

hereinafter attached assignment #3 FEM for mechanical engineering E563
please check it and give your comment!

jinqinghai

The first part is APDL commands file
ATTACHMENT—APDL file
1. APDL for the tank with flat end

!! setting the background of the window
WPSTYLE,,,,,,,,0
/RGB,INDEX,100,100,100, 0   
/RGB,INDEX, 80, 80, 80,13   
/RGB,INDEX, 60, 60, 60,14   
/RGB,INDEX, 0, 0, 0,15  

!recovery the original background color
! /RGB,INDEX, 0, 0, 0, 0  
! /RGB,INDEX, 60, 60, 60,13   
! /RGB,INDEX, 80, 80, 80,14   
! /RGB,INDEX,100,100,100,15   

/PMETH,OFF,0
KEYW,PR_SET,1   
KEYW,PR_STRUC,1

/UNITS,SI

/title, Assignment #3 E563

WPSTYLE,,,,,,,,0

FINISH  

/PREP7  
Diff=0.2  !!! tolerance for exact location

!ET,1,SOLID45  !! solid95 node limit
!Et,2,solid95
!et,1,plane42
et,1,shell93

MPTEMP,,,,,,,,  
MPTEMP,1,0  
MPDATA,EX,1,,200000
MPDATA,PRXY,1,,0.3  
MPDATA,DENS,1,,7.85E-6

FINISH  

/PREP7  
/REPLOT,RESIZE  
!K,1,500,,,  
K,1,502.5,1252.5,,  
k,3,502.5,1252.5-502.5,,
K,2,0,1252.5,,

!LFILLT,2,1,25, ,

!LSEL,ALL
KSEL,ALL

KSYMM,Y,ALL, , , ,0,0  

LSTR,1,2  
LSTR,1,3
lstr,3,6

LSTR,6,4
lstr,5,4

LFILLT,1,2,25, ,
!*  
LFILLT,4,5,25, ,

lsel,all

aROTAT,all, , , , , ,2,5 ,180,1,

allsel,all
/VIEW,1,1,1,1   
/ANG,1  

wpro,,,-90
!wpro,90
WPRO,180
wpro,180
CSWPLA,11,1,
cswpla,12,0

csys,12

k,200,0,-152.5,0
k,201,0,-152.5,-680
k,202,0,0,-680
k,203,0,0,0

L,200,201
L,201,202
!l,200,203

FLST,2,2,4,ORDE,2   
FITEM,2,21  
FITEM,2,-22
FLST,8,2,3  
FITEM,8,202
FITEM,8,203
AROTAT,P51X, , , , , ,P51X, ,180,1,

allsel,all  
aptn,all

asel,s,area,,10,13,3
ADELE,all, , ,1

/VIEW,  1, -0.424968791673    , -0.295430929554    , -0.855641333718
/ANG,   1,  -93.4021764054  

lplot
/image,save,model_line,bmp,

erase
allsel,all
aplot
/image,save,model_area,bmp,

R,1,5,  
R,2,50,

ASEL, , , ,       9
AATT,       1,       2,   1,       0,   

allsel,all
asel,u,area,,9
AATT,       1,       1,   1,       0,   

CSYS,0

LSEL,S,LOC,Z,-300+DIFF,-1250-5-DIFF
lesize,all,,,36

type,1
mat,1

SMRT,4
MSHAPE,1,2D
MSHKEY,0
asel,all
amesh,all

erase
eplot
/image,save,model_element,bmp,

/SOL
!*  
ANTYPE,0

!!! ############ APPLYING CONSTRAINTS

csys,0
lsel,s,loc,z,-diff,diff
dl,all,,symm

CSYS,0
kSEL,s,LOC,y,2500/2-500-5,2500/2-500+5
kSEL,R,LOC,X,-500-DIFF,-505+DIFF
nslk,s
D,all,ux,0

CSYS,0
kSEL,s,LOC,y,-(2500/2-500-5),-(2500/2-500+5)
kSEL,R,LOC,X,-500-DIFF,-505+DIFF
nslk,s

D,all,all,0

!! gravity applied

ACEL,9.8,0,0,  

!!! ######### APPLYING LOADS
!!! selecting inside surfaces to be loaded with pressure

allsel,all
asel,u,area,,7,14,7

SFA,all,1,PRES,2  !!! not 20 MPa!  

asel,s,area,,7,14,7
sfa,all,1,pres,-2

allsel,all

/PSF,PRES,NORM,2,0,1
/PBF,TEMP, ,1   
/PIC,DEFA, ,1   

aplot
/image,save,model_loadconst,bmp,

solve  !!! success!

!!!____________________________________
!!!  ##### STARTING POSTPROCESS  #####

*GET,total_Node_number,NODE,,NUM,MAX ! CURRENT MAXIMUM NODE NUMBER
*get,total_element_number,elem,,num,max !! CURRENT MAXIMUM elem NUMBER

!#### starting the process of postprocess  ####

allsel,all

/post1

/DSCALE,1,1
pldisp,1

/image,save,deflection_contour,bmp,

/DSCALE,1,1
PLNSOL, U,SUM, 0,1.0

/image,save,maxdisp,bmp, !!! displacement vector sum

erase
/post1
plnsol,s,1,
/image,save,principal1,bmp,

erase
/post1
plnsol,s,2,
/image,save,principal2,bmp,

erase
/post1
plnsol,s,3,
/image,save,principal3,bmp,

erase
/post1
plnsol,s,int,
/image,save,maxsint,bmp,

erase
/post1
plnsol,s,eqv,
/image,save,maxseqv,bmp,

!!! exploring local plot for nozzle

FLST,5,3,5,ORDE,3   
FITEM,5,9   
FITEM,5,14  
FITEM,5,-15
ASEL,S, , ,P51X
aplot   
nsla,s  
esla,s  

plnsol,s,int
/image,save,localsint,bmp,

plnsol,s,eqv,
/image,save,localseqv,bmp,

PLNSOL, U,SUM, 0,1.0
/image,save,localUsum,bmp,

pldisp,1
/image,save,localdisp,bmp,

!! ### define the max. value and the points of max. quantities
allsel,all

nsort,s,eqv
*get,Max_SEQV_value,sort,,max
*get,node_number_max_Seqv,sort,0,imax  !!! the node number of  max. Von Mises stress

nsort,s,int
*get,max_Sint_value,sort,,max
*get,node_number_max_Sint,sort,0,imax !!! the node number of max. stress intensity

nsort,u,sum
*get,max_Disp_value,sort,,max
*get,node_number_max_Disp,sort,0,imax !!! the node number of max. Displacement

*get, maxdisp_x,node,node_number_max_disp,loc,x
*get, maxdisp_y,node,node_number_max_disp,loc,y
*get, maxdisp_z,node,node_number_max_disp,loc,z
*get, maxseqv_x,node,node_number_max_seqv,loc,x
*get, maxseqv_y,node,node_number_max_seqv,loc,y
*get, maxseqv_z,node,node_number_max_seqv,loc,z
*get, maxsint_x,node,node_number_max_sint,loc,x
*get, maxsint_y,node,node_number_max_sint,loc,y
*get, maxsint_z,node,node_number_max_sint,loc,z

!!_______________________________________________________

!!! #### graph the von-mises stress distribution of the plate

allsel,all

/POST1  
CSYS,0
!!____________________________________________
!! ### point of maximum stress over the plate
!!! %%% checking -->list-->other-->parameters %%%

*create, report

*CFOPEN,reports,txt

  *VWRITE, TOTAL_ELEMENT_NUMBER
  ('THE TOTAL NUMBER OF ELEMENTS:', F9.0,TL1, ' ')

  *VWRITE, TOTAL_NODE_NUMBER
  ('THE TOTAL NUMBER OF NODES:', F9.0,TL1, ' ')

*VWRITE, MAX_DISP_VALUE, NODE_NUMBER_MAX_DISP
('THE VALUE OF MAXIMUM DISPLACEMENT:', E15.6, '   AT THE NODE:', F9.0,TL1,'  ')

   *vwrite, maxdisp_x,maxdisp_y,maxdisp_z
   ('THE LOCATION OF MAXIMUM DISPLACEMENT (X=', F7.1, ' , Y=', F7.1,' , z=', F7.1, ')')

*VWRITE, MAX_SEQV_VALUE, NODE_NUMBER_MAX_SEQV  
('THE VALUE OF MAXIMUM VON MISES STRESS:', F7.1, '       AT THE NODE:', F9.0,TL1,'  ')

  *vwrite, maxSEQV_x,maxSEQV_y,maxseqv_z
   ('THE LOCATION OF MAXIMUM VON MISES STRESS (X=', F7.1, ' , Y=', F7.1,' , z=', F7.1, ')')

*VWRITE, MAX_SINT_VALUE, NODE_NUMBER_MAX_SINT
('THE VALUE OF MAXIMUM STRESS INTENSITY:', F7.1, '       AT THE NODE:', F9.0,TL1,'  ')

  *vwrite, maxSINT_x,maxSINT_y,maxsint_z
   ('THE LOCATION OF MAXIMUM STRESS INTENSITY (X=', F7.1, ' , Y=', F7.1,' , z=', F7.1, ')')


*CFCLOS

*end  !!! must match "*creat"

!!/output, results, txt

*use,report

node1p1=2
node2p1=263

nsel,s,node,,node1p1
nsel,a,node,,node2p1

path,p1,2,30,20 !!!node1p1, node2p1
!!!path,p1,2,node1p1,node2p1
PPATH,1,node1p1
ppath,2,node2p1
PATH,STAT   
PATH,P1
AVPRIN,0, ,
PDEF, ,S,EQV,AVG
PLSECT,S,X, ,0  
/image,save,p1_path,bmp,

2. APDL for the tank with hemispherical end

!! setting the background of the window
WPSTYLE,,,,,,,,0
/RGB,INDEX,100,100,100, 0   
/RGB,INDEX, 80, 80, 80,13   
/RGB,INDEX, 60, 60, 60,14   
/RGB,INDEX, 0, 0, 0,15  

!recovery the original background color
! /RGB,INDEX, 0, 0, 0, 0  
! /RGB,INDEX, 60, 60, 60,13   
! /RGB,INDEX, 80, 80, 80,14   
! /RGB,INDEX,100,100,100,15   

/PMETH,OFF,0
KEYW,PR_SET,1   
KEYW,PR_STRUC,1

/REPLO  

/UNITS,SI

/title, Assignment #3 E563

WPSTYLE,,,,,,,,0

FINISH  

/PREP7  
Diff=0.2  !!! tolerance for exact location

!ET,1,SOLID45  !! solid95 node limit
!Et,2,solid95
!et,1,plane42
et,1,shell93

MPTEMP,,,,,,,,  
MPTEMP,1,0  
MPDATA,EX,1,,200000
MPDATA,PRXY,1,,0.3  
MPDATA,DENS,1,,7.85E-6
fINISH  

/PREP7  
K,1,502.5,1252.5,,  
k,2,502.5,1252.5-502.5,,
k,3,502.5,-1252.5,
K,4,0,1252.5,,  
k,5,0,-1252.5,,
k,6,502.5,-1252.5+502.5,,

CIRCLE,4,502.5 , , ,90,1,
allsel,all  
kplot   

FLST,3,1,4,ORDE,1   
FITEM,3,1   
LSYMM,Y,P51X, , , ,0,0  
lplot   
LSTR,       3,       6  
LSTR,       6,       2  
LSTR,       2,       1  

allsel,all

AROTAT,1, , , , , ,4,5 ,180,1,

AROTAT,1, , , , , ,4,5 ,180,1,

AROTAT,2, , , , , ,4,5 ,180,1,
AROTAT,3, , , , , ,4,5 ,180,1,
AROTAT,4, , , , , ,4,5 ,180,1,
AROTAT,5, , , , , ,4,5 ,180,1,

allsel,all
/VIEW,1,1,1,1   
/ANG,1  
/REP,FAST   
wpro,,,-90
!wpro,90
WPRO,180
wpro,180
CSWPLA,11,1,
cswpla,12,0

csys,12

k,200,0,-152.5,0
k,201,0,-152.5,-680
k,202,0,0,-680
k,203,0,0,0

L,200,201
L,201,202
!l,200,203

FLST,2,2,4,ORDE,2   
FITEM,2,21  
FITEM,2,-22
FLST,8,2,3  
FITEM,8,202
FITEM,8,203
AROTAT,P51X, , , , , ,P51X, ,180,1,

allsel,all  
aptn,all

asel,s,area,,9,12,3
ADELE,all, , ,1

/VIEW,  1, -0.424968791673    , -0.295430929554    , -0.855641333718
/ANG,   1,  -93.4021764054  

lplot
/image,save,model_line,bmp,

erase
allsel,all
aplot
/image,save,model_area,bmp,

R,1,5,  
R,2,50,

ASEL, , , ,       8
AATT,       1,       2,   1,       0,   

allsel,all
asel,u,area,,8
AATT,       1,       1,   1,       0,   

CSYS,0

LSEL,S,LOC,Z,-300+DIFF,-1250-5-DIFF
lesize,all,,,36

type,1
mat,1

SMRT,4
MSHAPE,1,2D
MSHKEY,0
asel,all
amesh,all

erase
eplot
/image,save,model_element,bmp,

/SOL
ANTYPE,0

!!! ############ APPLYING CONSTRAINTS

csys,0
lsel,s,loc,z,-diff,diff
dl,all,,symm

CSYS,0
kSEL,s,LOC,y,2500/2-500-5,2500/2-500+5
kSEL,R,LOC,X,-500-DIFF,-505+DIFF
nslk,s
D,all,ux,0

CSYS,0
kSEL,s,LOC,y,-(2500/2-500-5),-(2500/2-500+5)
kSEL,R,LOC,X,-500-DIFF,-505+DIFF
nslk,s

D,all,all,0

!! gravity applied

ACEL,9.8,0,0,  

!!! ######### APPLYING LOADS
!!! selecting inside surfaces to be loaded with pressure

allsel,all
asel,s,area,,3,8,5

SFA,all,1,PRES,2  !!! not 20 MPa!  

allsel,all
asel,u,area,,3,8,5
sfa,all,1,pres,-2

allsel,all

/PSF,PRES,NORM,2,0,1
/PBF,TEMP, ,1   
/PIC,DEFA, ,1   

aplot
/image,save,model_loadconst,bmp,

solve  !!! success!

!!!  ##### STARTING POSTPROCESS  #####

*GET,total_Node_number,NODE,,NUM,MAX ! CURRENT MAXIMUM NODE NUMBER
*get,total_element_number,elem,,num,max !! CURRENT MAXIMUM elem NUMBER

!#### starting the process of postprocess  ####

allsel,all

/post1

/DSCALE,1,1
pldisp,1
/image,save,deflection_contour,bmp,

/DSCALE,1,1
PLNSOL, U,SUM, 0,1.0

/image,save,maxdisp,bmp, !!! displacement vector sum

erase
/post1
plnsol,s,1,
/image,save,principal1,bmp,

erase
/post1
plnsol,s,2,
/image,save,principal2,bmp,

erase
/post1
plnsol,s,3,
/image,save,principal3,bmp,

erase
/post1
plnsol,s,int,
/image,save,maxsint,bmp,

erase
/post1
plnsol,s,eqv,
/image,save,maxseqv,bmp,

!!! exploring local plot for nozzle

FLST,5,3,5,ORDE,3   
FITEM,5,8   
FITEM,5,13  
FITEM,5,-14
ASEL,S, , ,P51X
aplot   
nsla,s  
esla,s  

plnsol,s,int
/image,save,localsint,bmp,

plnsol,s,eqv,
/image,save,localseqv,bmp,

PLNSOL, U,SUM, 0,1.0
/image,save,localUsum,bmp,

pldisp,1
/image,save,localdisp,bmp,

!! ### define the max. value and the points of max. quantities
allsel,all

nsort,s,eqv
*get,Max_SEQV_value,sort,,max
*get,node_number_max_Seqv,sort,0,imax  !!! the node number of  max. Von Mises stress

nsort,s,int
*get,max_Sint_value,sort,,max
*get,node_number_max_Sint,sort,0,imax !!! the node number of max. stress intensity

nsort,u,sum
*get,max_Disp_value,sort,,max
*get,node_number_max_Disp,sort,0,imax !!! the node number of max. Displacement

*get, maxdisp_x,node,node_number_max_disp,loc,x
*get, maxdisp_y,node,node_number_max_disp,loc,y
*get, maxdisp_z,node,node_number_max_disp,loc,z
*get, maxseqv_x,node,node_number_max_seqv,loc,x
*get, maxseqv_y,node,node_number_max_seqv,loc,y
*get, maxseqv_z,node,node_number_max_seqv,loc,z
*get, maxsint_x,node,node_number_max_sint,loc,x
*get, maxsint_y,node,node_number_max_sint,loc,y
*get, maxsint_z,node,node_number_max_sint,loc,z

!!! #### graph the von-mises stress distribution of the plate

allsel,all

/POST1  
CSYS,0
!!____________________________________________
!! ### point of maximum stress over the plate
!!! %%% checking -->list-->other-->parameters %%%

*create, report

*CFOPEN,reports,txt

  *VWRITE, TOTAL_ELEMENT_NUMBER
  ('THE TOTAL NUMBER OF ELEMENTS:', F9.0,TL1, ' ')

  *VWRITE, TOTAL_NODE_NUMBER
  ('THE TOTAL NUMBER OF NODES:', F9.0,TL1, ' ')

*VWRITE, MAX_DISP_VALUE, NODE_NUMBER_MAX_DISP
('THE VALUE OF MAXIMUM DISPLACEMENT:', E15.6, '   AT THE NODE:', F9.0,TL1,'  ')

   *vwrite, maxdisp_x,maxdisp_y,maxdisp_z
   ('THE LOCATION OF MAXIMUM DISPLACEMENT (X=', F7.1, ' , Y=', F7.1,' , z=', F7.1, ')')

*VWRITE, MAX_SEQV_VALUE, NODE_NUMBER_MAX_SEQV  
('THE VALUE OF MAXIMUM VON MISES STRESS:', F7.1, '       AT THE NODE:', F9.0,TL1,'  ')

  *vwrite, maxSEQV_x,maxSEQV_y,maxseqv_z
   ('THE LOCATION OF MAXIMUM VON MISES STRESS (X=', F7.1, ' , Y=', F7.1,' , z=', F7.1, ')')

*VWRITE, MAX_SINT_VALUE, NODE_NUMBER_MAX_SINT
('THE VALUE OF MAXIMUM STRESS INTENSITY:', F7.1, '       AT THE NODE:', F9.0,TL1,'  ')

  *vwrite, maxSINT_x,maxSINT_y,maxsint_z
   ('THE LOCATION OF MAXIMUM STRESS INTENSITY (X=', F7.1, ' , Y=', F7.1,' , z=', F7.1, ')')

*CFCLOS

*end  !!! must match "*creat"

!!/output, results, txt

*use,report

ASEL,S, , ,       1
nsla,s  
esla,s  

/post1  
plnsol,s,EQV
/image,save,hemisphericalendSEAV,bmp,
plnsol,s,int
/image,save,hemisphericalendSINT,bmp,
plnsol,U,SUM
/image,save,hemisphericalendUsum,bmp,

jinqinghai

report attachment

jinqinghai

image1