/filename,ex6-6                    		!趨
/title,The Static Analysis of the RPC    	!趨
/prep7                            		!ǰģ
!***************************
!嵥Ԫ͡ʵͲϲ
!***************************
et,1,solid45                        	!嵥ԪSolid45
et,2,link8                          	!嵥ԪLink8
R,1,137.4e-6,2.5e-4 ,                	!ֽʵͳʼλƣԤӦ
!ϲ
mp,ex,1,4e10                     	!ˮܵϲΪ1
mp,prxy,1,0.17 
mp,dens,1,2500
mp,ex,2,2.0e11                    	!ˮܵԤӦֽϲΪ2
mp,prxy,2,0.3 
mp,dens,2,7800
!***************************
!ģ
!***************************
csys,1                           	!ѵǰϵתΪϵ
!ؼڵ
n,1,0.58                          	!ھƽڵ
n,2,0.63                          	!ֽƽڵ
n,3,0.7                          	!⾶ƽڵ
NGEN,2,3,1,3,1,0,18,0,1,           	!ƹؼƽڵ1~3
NGEN,2,6,1,6,1,0,0,0.05,1,          	!ƹؼƽڵ1~6
!ܵԪ
type,1                            	!Ԫͱ1
mat,1                                   !ϲ1
e,1,2,5,4,7,8,11,10                 	!Ԫ1
e,2,3,6,5,8,9,12,11                 	!Ԫ2 
eplot 
egen,20,12,1,2,1,,,,,,,18,,            	!ƵԪˮܵ浥Ԫ
nummrg,all                       	!Ԫڵϲ
numcmp,all                              !ѹ
egen,120,120,1,40,1,,,,,,,,0.05        	!ƵԪˮܵȫԪ
nummrg,all                      	!Ԫڵϲ
numcmp,all                              !ѹ
!ֽԪ
type,2                         	!øֽԪͺ
mat,2                          	!øֽԪϱ
!19Ԫ
*dim,nOne,array,21              	!ʼnone
*do,i,1,19                      	!iΪѭDOѭʼ
*do,j,1,21                      	!jΪѭǶDOѭʼ
nsel,s,loc,x,0.63                 	!ѡx=0.63Ľڵ
nsel,r,loc,z,0.3*i                 	!ӵǰѡѡz=0.3*iĽڵ
nsel,r,loc,y,18*(j-1)              	!ӵǰѡѡy=18*(j-1)Ľڵ
*get,nOne(j),node,,num,max       	!ڵŴnOne
*enddo                        	        !jΪѭǶDOѭ
*do,j,1,20                      	!jΪѭǶDOѭʼ
e,nOne(j),nOne(j+1)              	!ڵɵԪ
*enddo                        	        !jΪѭǶDOѭ
*enddo                         	        !iΪѭǶDOѭ
!10ݽԪ
*dim,nTwo,array,121              
*do,i,1,10 
*do,j,1,121 
nsel,s,loc,x,0.63 
nsel,r,loc,y,36*(i-1) 
nsel,r,loc,z,0.05*(j-1) 
*get,nTwo(j),node,,num,max 
*enddo 
*do,j,1,120 
e,nTwo(j),nTwo(j+1) 
*enddo 
*enddo 
Alls                         !ѡеͼԪ
nummrg,all                   !ϲغϵĲ
numcmp,all                   !ѹеĲ
!***************************
!ؼ
!***************************
/solu                            	!ģ
Antype,0                         	!趨Ϊ
!ʩλԼ
Nsel,s,loc,z,0,0.3
Nsel,r,loc,x,0.7                    	!ѡ˵Ľڵ
D,all,all                          	!ԽڵʩȫԼ
Allsel
Nsel,s,loc,z,5.7,6
Nsel,r,loc,x,0.7
D,all,all
Allsel
!ʩѹغ
Nsel,s,loc,x,0.58                    	!ѡڱĽڵ
Sf,all,PRES,1e5                         !ʩ1e5
Allsel
Nsel,s,loc,y,0,180                   		
Nsel,r,loc,x,0.7                        !ѡĽڵ
Sf,all,PRES,3e5                        	!ʩ3e5
Csys,0                                	!תǰϵΪֱϵ
Acel,,10                          	!ʩٶ
Allsel                                 	!ѡȫͼԪ
Solve                            	!ʼ
!***************************
!鿴
!***************************
/post1                                
Csys,1                            			!лǰϵΪϵ
PLNSOL, U,X, 0,1.0                 		!ʾXλͼ
PLNSOL, U,Y, 0,1.0                 		!ʾYλͼ
PLNSOL, U,Z, 0,1.0                 		!ʾZλͼ
PLNSOL, U,SUM, 0,1.0              		!ʾλͼ
PLNSOL, S,X                      		!ʾXӦͼ
PLNSOL, S,Y                      		!ʾYӦͼ
PLNSOL,S,1                    		        !ʾһӦͼ
PLNSOL,S,3                    		        !ʾӦͼ
/EXIT,ALL                            		!ȫݲ˳
