XANSYS Forum Index
    Register    

FAQ    Search    Memberlist    Usergroups    SmartFeedSmartFeed    Profile    Log in
Inistate not working for strain state
 
This forum is locked: you cannot post, reply to, or edit topics.   This topic is locked: you cannot edit posts or make replies.    XANSYS Forum Index -> XANSYS
Author Message
abdullah.al.amin
User


Joined: 23 Oct 2014
Posts: 12

PostPosted: Tue Jan 26, 2016 9:51 pm  Reply with quote

Hi Everone,

I have to solve a prestressed cylinder which has seven layers. Hence, I subsequently used inistate command to apply prestressed condition in azimuthal direction of a 2D axisymmetric setup. But the strain values are not right. They are non-zero even if the stress is zero.

Alternately, if I use initial strain state instead of stress; the results make sense as strain is zero where stresses are zero. But, this seems to be not working for later layers. It appears that the later cylinder layers (layer 2 and onwards) has no effect on the analysis whatsoever.

Can somebody please help me figure out what is going wrong here? What it is I am missing for the strain state?

Regards,
Abdullah Al Amin
Graduate Research Assistant
Case Western Reserve University

Code:

FINISH
/clear, nostart

!CHANGE DIRECTORY
/MKDIR, Test1
/CWD, D:\Abdullah\APDL
/DIRECTORY,,HOOPSTRESS,,
/filename, Hoop Stress, 1
/com, STRESS MODELING
/prep7
/title, STRESS MODELING


ET,1,PLANE182,,,1

MP, EX, 1, 100E9
MP, MURX, 1, 1
MP, PRXY, 1, 0.3
MP, ALPX, 1, 9.95E-6


MP, EX, 2, 200E9
MP, MURX, 2, 1
MP, NUXY, 2, 0.305
MP, ALPX, 2, 1.2E-5


LAYER = 7
HEIGHT = 1e-3
R_i = 0.90
R_o = 0.93
X_COORD = R_i

Z_1 = 0.61
Z_2 = 0.80

WIDTH = (R_o-R_i)/Layer
Height = Z_2-Z_1

MANDREL_Thick = 10e-3
MANDREL_RADIUS = R_i-MANDREL_Thick

Y_COORD = 0
BLC4, MANDREL_RADIUS,Y_COORD, MANDREL_Thick, Height

*DO, j, 1, LAYER

Y_COORD=0
BLC4, X_COORD, Y_COORD, WIDTH, HEIGHT
numcmp, all
X_COORD=X_COORD+WIDTH
*ENDDO

!formation of mandrel


AGLUE, ALL
numcmp, all


!ASSIGN MATERIAL AND ELEMENT

ASEL, S, AREA,,1
AATT,2,,1
APLOT

ASEL, s, AREA, ,2,LAYER+1
AATT,1,,1

ASEL, ALL
APLOT

!ASSIGN AREA MESH SIZE
AESIZE, ALL,1E-3

AMESH, ALL

ANTYPE, 0
NLGEOM, ON
NROPT, FULL
PSTRES, OFF

FINISH



/SOLU

SOLCONTROL, ON
CNVTOL,U,1e-7,0.001,2,0.001
NSUBST, 1

LSEL, S, LOC,Y,0,0
DL, ALL,,SYMM


ALLSEL, ALL



!Assigning Hoop Stress on Elements

ASEL, S, AREA, ,2, LAYER+1
ESLA
EPLOT
!inistate,defi,,,,,,,3.11482215e7
inistate,set,dtyp,epel
inistate,defi,,,,,,,1.68e-4


current_t=1

ASEL, S, AREA,,2,LAYER+1
APLOT
ESLA
EKILL, ALL


*DO, i, 1,LAYER
TIME, i
ASEL, S, AREA,,i+1
ESLA
EALIVE, ALL

ALLSEL, ALL
SAVE
SOLVE

UPCOORD, 1, ON
*ENDDO

/post1
PATH, CYL, 2,30,30
PPATH,1,,Mandrel_radius,(Z_2-Z_1)/2,0
PPATH,2,,R_o-width/3,(Z_2-Z_1)/2,0
PDEF,S_X,S,X
PDEF,S_z,S,Z

PDEF,e_r_tot,EPTO,x
PDEF,e_z_tot,EPTO,z

PLPATH, S_z
Back to top
View user's profile Send private message
Display posts from previous:   
This forum is locked: you cannot post, reply to, or edit topics.   This topic is locked: you cannot edit posts or make replies.    XANSYS Forum Index -> XANSYS
All times are GMT - 7 Hours
Page 1 of 1

 
Jump to:  
You cannot post new topics in this forum
You cannot reply to topics in this forum
You cannot edit your posts in this forum
You cannot delete your posts in this forum
You cannot vote in polls in this forum
You cannot attach files in this forum
You cannot download files in this forum


Powered by phpBB © 2001, 2005 phpBB Group

sleek template created by Andrew Charron