/prep7  



*set,Q,0.001
*set,R,0.001

*set,P0,100
*set,F,1
*set,weight,1000

fini


!*DIM,E0(1),,20
*DIM,E_2(i),,20
*DIM,P_2(i),,20
*set,E0,2000



*Do,i,1,50





*Do,j,1,5


/prep7

	
et,1,LINK1  		! Element type no.1 is LINK1
mp,eX,1,E0		! Material Properties: E for material no.
mp,prxy,1,0.3		! poisson's ratio for material no.1
r,1,0.4			! Real Const.ant number 1 is 0.125
			! (Cross sect.ional area)

n,1,0,0,,   
n,2,20,0,,  
n,3,40,0,,  
n,4,60,0,,  
n,5,80,0,,  
n,5,80,0,,  
n,6,0,20,,  
n,7,20,20,, 
n,8,40,20,, 
n,9,60,20,, 
n,10,80,20,,
n,11,0,40,, 
n,12,20,40,,
n,13,40,40,,
n,14,60,40,,
n,15,80,40,,
n,16,0,60,, 
n,17,20,60,,
n,18,40,60,,
n,19,60,60,,
n,20,80,60,,
n,21,0,80,, 
n,22,20,80,,
n,23,40,80,,
n,24,60,80,,
n,25,80,80,,


!------------------Element connects nodes together------------------------------

en,1,1,2
en,2,2,3
en,3,3,4
en,4,4,5
en,5,5,10
en,6,10,15
en,7,15,20
en,8,20,25
en,9,25,24
en,10,24,23
en,11,23,22
en,12,22,21
en,13,21,16
en,14,16,11
en,15,11,6
en,16,6,1
en,17,1,7
en,18,7,13
en,19,13,19
en,20,19,25
en,21,16,22
en,22,11,17
en,23,17,23
en,24,6,12
en,25,12,18
en,26,18,24
en,27,2,8
en,28,8,14
en,29,14,20
en,30,3,9
en,31,9,15
en,32,4,10
en,33,6,7
en,34,7,8
en,35,8,9
en,36,9,10
en,37,11,12
en,38,12,13
en,39,13,14
en,40,14,15
en,41,16,17
en,42,17,18
en,43,18,19
en,44,19,20
en,45,2,7
en,46,7,12
en,47,12,17
en,48,17,22
en,49,3,8
en,50,8,13
en,51,13,18
en,52,18,23
en,53,4,9
en,54,9,14
en,55,14,19
en,56,19,24

!----------------------BC-----------------------
d,1,ux,0
d,1,uy,0
d,5,ux,0
d,5,uy,0
d,6,ux,0
d,11,ux,0
d,16,ux,0
d,21,ux,0
d,10,ux,0
d,15,ux,0
d,20,ux,0
d,25,ux,0
d,2,uy,0
d,3,uy,0
d,4,uy,0


!----------------------load applying---------------------
f,23,fy,-10
FINISH


!----------------------Solution------------------------
/SOLU

ANTYPE,4
 
TRNOPT,FULL 
LUMPM,0 
  
DELTIM,1,0,0  
TIME,5 
nlgeom,on
	
SOLVE   



	*DIM,disp2000,ARRAY,5,1,0, , ,   
	*VGET,disp2000(1.1),NODE,17,U,Y, , ,2     
	*VGET,disp2000(2.1),NODE,22,U,Y, , ,2
	*VGET,disp2000(3.1),NODE,23,U,Y, , ,2
	*VGET,disp2000(4.1),NODE,24,U,Y, , ,2     
	*VGET,disp2000(5.1),NODE,19,U,Y, , ,2 




	*DIM,DISP1000,ARRAY,5,1,0, , ,  
	*SET,DISP1000(1,1) , -0.308830   
	*SET,DISP1000(2,1) , -0.296424 
	*SET,DISP1000(3,1) , -0.858322 
	*SET,DISP1000(4,1) , -0.255089  
	*SET,DISP1000(5,1) , -0.191081  


	*DIM,E(j),,5
	
	*dim,P(j),,5
	!*dim,p1(1),,1

	!*DIM,B,,1  

        p1(1)=F*P0*F+Q
	

	m=disp2000(j)*p1(1)*disp2000(j)+R
	k=p1(1)*disp2000(j)*(1/m)
	C=1-k*disp2000(j)

	E(j)=E0+k*(-DISP1000(j)+disp2000(j))

	*set,E0,E(j)

*ENDDO

	E_2(i)=E(5)
	E0=E_2(i)
	P_2(i)=C*p1(1)*C+k*R*k
	P0=weight*P_2(i)




*ENDDO

!---------------------ploting the results from nodes 17,19,22,23,24--------------------

/post26
nsol,2,17,U,Y
nsol,3,22,U,Y
nsol,4,23,U,Y
nsol,5,24,U,Y
nsol,6,19,U,Y