APDL PROBLEM

[SKJoe]

Hi, friends

I want to put initial deformation to my finite model by the formula written below ( parameter WW). By this formula is obtained deflection, which is consequently issign to desired nod...

Is there anyone so kind to explain my why this do-loops doesnt work ? - only the first cycle is executed...

P.S. sorry for my bad english

LX=3000 ! LENGTH X
LY=1000 ! LENGHT Y
XX=30 ! NUMBER OF ELEMENTS IN X DIRECTION
YY=10 ! NUMBER OF ELEMENTS IN Y DIRECTION

KEYW,PR_SET,1
KEYW,PR_STRUC,1

/PREP7

ET,1,SHELL181

KEYOPT,1,1,0
KEYOPT,1,3,2
KEYOPT,1,8,2
KEYOPT,1,9,0
KEYOPT,1,10,0

MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,EX,1,,210000
MPDATA,PRXY,1,,.3

sect,1,shell,,HRUBKA
secdata, 10,1,0.0,3
secoffset,MID
seccontrol,,,, , , ,

RECTNG,,LX,,LY,

FLST,5,2,4,ORDE,2
FITEM,5,2
FITEM,5,4
CM,_Y,LINE
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,S,_Y

AATT, 1, , 1, 0, 1

FLST,5,2,4,ORDE,2
FITEM,5,2
FITEM,5,4
CM,_Y,LINE
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y
!*
LESIZE,_Y1, , ,YY, , , , ,1

FLST,5,2,4,ORDE,2
FITEM,5,1
FITEM,5,3
CM,_Y,LINE
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y
!*
LESIZE,_Y1, , ,XX, , , , ,1

CM,_Y,AREA
ASEL, , , , 1
CM,_Y1,AREA
CHKMSH,'AREA'
CMSEL,S,_Y
MSHKEY,1
AMESH,_Y1
MSHKEY,0
CMDELE,_Y
CMDELE,_Y1
CMDELE,_Y2

NWRITE,'NODE_01',' ',' ',0

FINISH

/CLEAR,NOSTART

/PREP7

NRRANG,1,99999999,1,
NREAD,'NODE_01',' ',' '

!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

LX=3000 ! LENGTH X
LY=1000 ! LENGHT Y
XX=30 ! NUMBER OF ELEMENTS IN X DIRECTION
YY=10 ! NUMBER OF ELEMENTS IN Y DIRECTION

W1=950.9
W2=700
W3=1050
PI=3.141592654

*DIM,RY,TABLE,(YY-1),1,1, , ,
*SET,RY(0,1,1) , 1
*SET,RY((1YY-1)),1,1) ,((LY/YY)LY-2*(LY/YY))LY/YY))
*SET,RY((1YY-1)),0,1) ,(1YY-1))

*DIM,RX,TABLE,(XX-1),1,1, , ,
*SET,RX(0,1,1) , 1
*SET,RX((1XX-1)),1,1) ,((LX/XX)LX-2*(LX/XX))LX/XX))
*SET,RX((1XX)-1),0,1) ,(1XX-1))

*DO,DY,1,(YY-1),1

*DO,DX,1,(XX-1),1

WW=(W1*(SIN(PI*(RX(DX))/LX))+W2*(SIN(2*PI*(RX(DX))/LX))+W3*(SIN(3*PI*(RX(DX))/LX)))*(SIN(PI*(RY(DY))/LY))

NSEL,S,LOC,Y,((RY(DY))-0.3*LY/YY),((RY(DY))+0.3*LY/YY)
NSEL,R,LOC,X,((RX(DX))-0.3*LX/XX),((RX(DX))+0.3*LX/XX)

*GET,NOD,NODE,0,NUM,MAX

FLST,2,1,1
FITEM,2,NOD
FLST,3,1,8
FITEM,3,(RX(DX)),(RY(DY)),WW
NMODIF,P51X,P51X, , , , , ,

ALLSEL,ALL

*ENDDO

*ENDDO

 

[SKJoe]

...I have realized that mistake is not in do-loop, but in filling table for that do-loop. So why this prescription doesnt work ?

*DIM,RX,TABLE,100,1,
*SET,RX(0,1,1),1
*SET,RX((1:100),1,1),(10:1000:10)
*SET,RX((1:100),0,1),(1:100)