{"id":314791,"date":"2023-11-03T15:33:55","date_gmt":"2023-11-03T15:33:55","guid":{"rendered":"\/forum\/forums\/topic\/how-to-select-elements-in-elliptical-shape-in-a-cylinder-using-ansys-apdl\/"},"modified":"2023-11-06T12:54:04","modified_gmt":"2023-11-06T12:54:04","slug":"how-to-select-elements-in-elliptical-shape-in-a-cylinder-using-ansys-apdl","status":"closed","type":"topic","link":"https:\/\/innovationspace.ansys.com\/forum\/forums\/topic\/how-to-select-elements-in-elliptical-shape-in-a-cylinder-using-ansys-apdl\/","title":{"rendered":"how to select elements in elliptical shape in a cylinder using ansys apdl?"},"content":{"rendered":"

Dear all,<\/p>\n

I am working on damage detection on a composite cylinder. the idea is to insert damage with an elliptical shape in a cylinder in a different location and use the buckling response as a damage indicator that later will be inserted into a CNN-based algorithm for damage detection. I have modeled the cylinder and now want to insert the elliptical shape damage (which is selecting elements in a defined elliptical shape and assigning them a low rigidity). unfortunately, I didn’t succeed till now.<\/p>\n

Any hints or ideas are welcome
Ps: the code I did is below<\/p>\n

FINISH<\/p>\n

FINISH<\/p>\n

\/CLEAR, NOSTART ! NOSTART avoids reading again START.ANS
\/CWD,’D:\\Damage_Detection_CNN\\NewLastCode’
\/RGB,INDEX,100,100,100,0
\/RGB,INDEX,0,0,0,15
\/PREP7
\/TRIAD,LBOT
\/VIEW,1,1,1,1
\/ANGLE,1
\/VUP,,Z
\/UNITS,SI
\/ESHAPE,0
*AFUN,DEG ! Use degrees for input and output of parameter ANGular functions
DmageON = 1 ! 1: Damage is activated
! 0: Damage is Desactivated
C*** Input data
C*** RUN<\/p>\n

*SET, PI, 3.14159265359<\/p>\n

! Ler variaveis externas
PARRES, change, variaveis.txt<\/p>\n

Dc= 0.05 ! [m] Diameter of cylinder
Lc= 0.3 ! [m] Length of the cylinder
Th= 0.00012 ! [m] Wall thickness
ElSiz = 2e-3 ! [mm] Element Size<\/p>\n

C*** Material proprties ************************************<\/p>\n

*SET,N_camadas,10
*SET,ANG1,45.0000
*SET,ANG2,0.0000
*SET,ANG3,-45.0000
*SET,ANG4,0.0000
*SET,ANG5,45.0000
*SET,ANG6,45.0000
*SET,ANG7,0.0000
*SET,ANG8,-45.0000
*SET,ANG9,0.0000
*SET,ANG10,45.0000<\/p>\n

*SET,E1,101860000000.0000
*SET,E2,3410000000.0000
*SET,E3,3410000000.0000
*SET,v12,0.2900
*SET,v23,0.2900
*SET,v13,0.2900
*SET,G12,7560000000.0000
*SET,G23,1321705426.3566
*SET,G13,7560000000.0000
*SET,rho,1510.0000
*SET,E_Layer,Th<\/p>\n

C*** Force ************************************
*SET,Fcomp,-1<\/p>\n

*SET,RedictioRatio,2000<\/p>\n

C*** Failure Criteria ************************************
FC,1,S,XTEN,1363.49e6
FC,1,S,YTEN,5.86e6
FC,1,S,ZTEN,5.86e6
FC,1,S,XCMP,-572.27e6
FC,1,S,YCMP,-102e6
FC,1,S,ZCMP,-102e6
FC,1,S,XY,200.61e6
FC,1,S,YZ,200.61e6
FC,1,S,XZ,200.61e6
FC,1,S,XYCP,-1
FC,1,S,YZCP,-1
FC,1,S,XZCP,-1
FC,1,EPEL,XTEN,
FC,1,EPEL,YTEN,
FC,1,EPEL,ZTEN,
FC,1,EPEL,XCMP,
FC,1,EPEL,YCMP,
FC,1,EPEL,ZCMP,
FC,1,EPEL,XY,
FC,1,EPEL,YZ,
FC,1,EPEL,XZ,<\/p>\n

FC,2,S,XTEN,1363.49e6
FC,2,S,YTEN,5.86e6
FC,2,S,ZTEN,5.86e6
FC,2,S,XCMP,-572.27e6
FC,2,S,YCMP,-102e6
FC,2,S,ZCMP,-102e6
FC,2,S,XY,200.61e6
FC,2,S,YZ,200.61e6
FC,2,S,XZ,200.61e6
FC,2,S,XYCP,-1
FC,2,S,YZCP,-1
FC,2,S,XZCP,-1
FC,2,EPEL,XTEN,
FC,2,EPEL,YTEN,
FC,2,EPEL,ZTEN,
FC,2,EPEL,XCMP,
FC,2,EPEL,YCMP,
FC,2,EPEL,ZCMP,
FC,2,EPEL,XY,
FC,2,EPEL,YZ,
FC,2,EPEL,XZ,<\/p>\n

C*** Assigne Material Proprties ************************************<\/p>\n

MP,EX,1,E1
MP,EY,1,E2
MP,EZ,1,E3
MP,PRXY,1,V12
MP,PRYZ,1,V23
MP,PRXZ,1,V13
MP,GXY,1,G12
MP,GYZ,1,G23
MP,GXZ,1,G13
MP,DENS,1,rho<\/p>\n

E11=RedictionRatio*E1
E22=RedictionRatio*E2
E33=RedictionRatio*E3
G122=RedictionRatio*G12
G232=RedictionRatio*G23
G132=RedictionRatio*G13<\/p>\n

MP,EX,2,E11
MP,EY,2,E22
MP,EZ,2,E33
MP,PRXY,2,V12
MP,PRYZ,2,V12
MP,PRXZ,2,V12
MP,GXY,2,G122
MP,GYZ,2,G232
MP,GXZ,2,G132
MP,DENS,2,RedictionRatio*2100<\/p>\n

C*** Material Composite ************************************<\/p>\n

ET,1,SHELL281
SECT,1,SHELL,,
SECDATA, E_Layer,1,ANG1,3
SECDATA, E_Layer,1,ANG2,3
SECDATA, E_Layer,1,ANG3,3
SECDATA, E_Layer,1,ANG4,3
SECDATA, E_Layer,1,ANG5,3
SECDATA, E_Layer,1,ANG6,3
SECDATA, E_Layer,1,ANG7,3
SECDATA, E_Layer,1,ANG8,3
SECDATA, E_Layer,1,ANG9,3
SECDATA, E_Layer,1,ANG10,3
secoffset,MID
seccontrol,,,, , , ,
SECPLOT, 1,,,<\/p>\n

ET,2,SHELL281
SECT,2,SHELL,,
SECDATA, E_Layer,2,ANG1,3
SECDATA, E_Layer,2,ANG2,3
SECDATA, E_Layer,2,ANG3,3
SECDATA, E_Layer,2,ANG4,3
SECDATA, E_Layer,2,ANG5,3
SECDATA, E_Layer,2,ANG6,3
SECDATA, E_Layer,2,ANG7,3
SECDATA, E_Layer,2,ANG8,3
SECDATA, E_Layer,2,ANG9,3
SECDATA, E_Layer,2,ANG10,3
secoffset,MID
seccontrol,,,, , , ,<\/p>\n

C*** Cylinder Geomtry ************************************<\/p>\n

K,1,0.,0.,0. ! CENT Point
K,2,0.,0.,Lc ! AXIS Point
K,3,0.,Dc\/2,0. ! ZERO Point
CIRCLE,1,Dc\/2,2,3,360 ! CIRCLE, PCENT, RAD, PAXIS, PZERO, ARC, NSEG
L,1,2 ! Line for Path in following ADRAG –> Line N.5 (4 Lines used for the CIRCLE)
ADRAG,ALL,,,,,,5<\/p>\n

C*** Meshing (Shells) ************************************<\/p>\n

ASEL,ALL
ESIZE,ElSiz
AMESH,ALL<\/p>\n

C*** Damage ************************************<\/p>\n

!pos_z = Lc\/3
!a = 0.02
!b = 0.01
!theta = 10
!r = 1<\/p>\n

*IF,DmageON,EQ,1,then
LOCAL,11,1,0,Dc\/2,pos_z,,,theta,r
CSYS, 11
ESEL, S, CENT, 0.01, 0.00001, b, , , , 0
!ESEL, S, CENT, Y, 0, (b\/(10*r)),,,,0
EMODIF,ALL,MAT,2,,,,,,
EMODIF,ALL,Type,2,,,,,
EMODIF,ALL,SECT,2,,,,,
ESEL,ALL
!CSYS,0
*EndIF
Allsel,all
Finish<\/p>\n

C*** Analysis ************************************<\/p>\n

\/SOL
Lsel,s,line,,1,4
DL,ALL, ,ALL,
Allsel,all
!! Analysis
ANTYPE,0
Lsel,s,line,,6
Lsel,a,line,,9
Lsel,a,line,,11
Lsel,a,line,,13
NSLL, S, 1
F,all,FZ,Fcomp<\/p>\n

Allsel,All<\/p>\n

PSTRES,1
SOLVE
FINISH
\/SOLUTION
ANTYPE,1
BUCOPT,LANB,20,0,0,CENTER
SOLVE
FINISH<\/p>\n

\/POST1
\/output,
\/output,’resultado’,txt<\/p>\n

SET,LIST
\/out<\/p>\n

\"\"<\/a><\/p>\n","protected":false},"template":"","class_list":["post-314791","topic","type-topic","status-closed","hentry","topic-tag-ansys-mechanical-apdl","topic-tag-cylinder","topic-tag-damage-1","topic-tag-material-damage"],"aioseo_notices":[],"acf":[],"custom_fields":[{"0":{"_bbp_subscription":["39239","20905"],"_bbp_author_ip":["23.206.193.146"]," _bbp_last_reply_id":["0"]," _bbp_likes_count":["0"],"_btv_view_count":["1270"],"_edit_lock":["1699274805:221391"],"_bbp_topic_status":["unanswered"],"_bbp_status":["publish"],"_bbp_topic_id":["314791"],"_bbp_forum_id":["27791"],"_bbp_engagement":["20905","39239"],"_bbp_voice_count":["2"],"_bbp_reply_count":["6"],"_bbp_last_reply_id":["315296"],"_bbp_last_active_id":["315296"],"_bbp_last_active_time":["2023-11-07 16:48:14"]},"test":"ikram-djabrouhouuniv-dbkm-dz"}],"_links":{"self":[{"href":"https:\/\/innovationspace.ansys.com\/forum\/wp-json\/wp\/v2\/topics\/314791","targetHints":{"allow":["GET"]}}],"collection":[{"href":"https:\/\/innovationspace.ansys.com\/forum\/wp-json\/wp\/v2\/topics"}],"about":[{"href":"https:\/\/innovationspace.ansys.com\/forum\/wp-json\/wp\/v2\/types\/topic"}],"version-history":[{"count":0,"href":"https:\/\/innovationspace.ansys.com\/forum\/wp-json\/wp\/v2\/topics\/314791\/revisions"}],"wp:attachment":[{"href":"https:\/\/innovationspace.ansys.com\/forum\/wp-json\/wp\/v2\/media?parent=314791"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}