/filename,ex5-1                           	!趨
/title,The Transient Analysis of the pipe       !趨
/prep7                                   	!ǰģ
*set,dia,139.7                             	!
*set,thick,7.72
*set,hole_dia,12.7
*set,hole_length,25.4
*set,hole_num,10
*set,hole_ang,18
!
et,1,solid95                               	!ѡԪSOLID95 
mp,ex,1,2.1e5                             	!嵯ģ
mp,prxy,1,0.3                             	!岴ɱ 
mp,dens,1,0.00785                         	!ܶ
wpoff,0,0,-hole_length                     	!-zƫƹƽhole_length
cyl4,0,0,dia/2-thick,,dia/2,,hole_length        !һԲ
wpoff,0,0,hole_length/2                     	!+zƫƹƽhole_length
wpro,,,90                                 	!yתƽ90
cyl4,0,0,0,,hole_dia/2,,100                    	!Բ
!
csys,1                                    	!ϵ
*do,i,1,hole_num                           	!doѭʼ
vgen,2,2,,,,i*360/hole_num,0                  	!ѭԲ
*enddo                                   	!doѭ
!
vovlap,all                                 	!едӲ
allsel,all                            		!ѡȫͼԪ
/replot
*get,max_number,volu,0,num,max                  !ͼԪڲmax_number
vdele,1,max_number-1,1,1              		!ɾΪ1~max_number-1ͼԪ
numcmp,all                          		!ѹ
!
vgen,5,1,,,,-hole_ang,-hole_length,0       	!Բ
vglue,all                             		!ճе
/replot 
!mesh the model
lsel,s,loc,z,0                          	!ѡz=0 
lsel,a,loc,z,-5*hole_length               	!ͬʱѡz=-5*hole_length
lesize,all,,,10                          	!趨ߵķָ10
smrtsize,8                            		!趨ֵܻȼ8
allsel                                		!ѡȫͼԪ
mshkey,0                             		!趨񻮷ַʽ
mshape,1,3d                          	        !趨Ϊ
vmesh,all                             	        !񻮷
!node component
nsel,s,loc,z,0                           	!ѡz=0Ľڵ
cm,zuoduanmian,node                   	        !ɽڵ
nsel,s,loc,z,-5*hole_length                	!ѡz=-5*hole_lengthĽڵ
cm,youduanmian,node                   	        !ɽڵ
!
!apply the displacement
nsel,s,loc,z,-hole_length               !ѡz=-hole_lengthĽڵ
nsel,a,loc,z,-4*hole_length             !ͬʱѡz=-4*hole_lengthĽڵ
nsel,r,loc,x,dia/2                      !ڵǰѡѡx=dia/2Ľڵ
nsel,r,loc,y,210,330                    !ڵǰѡѡyλ210~330Ľڵ
d,all,all,0                             !ԵǰѡнڵɶʩԼ
/solu                                  	!ģ
antype,trans                            !趨
*dim,forcez1,table,13,1,1,,,            !ѭغɱforcez1
!time
*set,forcez1(1,0,1) , 0                 !趨ѭغɱʱֵ
*set,forcez1(2,0,1) , 0.0025 
*set,forcez1(3,0,1) , 0.0075
*set,forcez1(4,0,1) , 0.01  
*set,forcez1(5,0,1) , 0.0125
*set,forcez1(6,0,1) , 0.0175 
*set,forcez1(7,0,1) , 0.02
*set,forcez1(8,0,1) , 0.0225  
*set,forcez1(9,0,1) , 0.0275
*set,forcez1(10,0,1) , 0.03
*set,forcez1(11,0,1) , 0.0325   
*set,forcez1(12,0,1) , 0.0375 
*set,forcez1(13,0,1) , 0.04   
!+z 
*set,forcez1(1,1,1) , 0                       !趨ѭغɱֵ
*set,forcez1(2,1,1) , 191.82 
*set,forcez1(3,1,1) , 191.82  
*set,forcez1(4,1,1) , 0 
*set,forcez1(5,1,1) , -191.82 
*set,forcez1(6,1,1) , -191.82 
*set,forcez1(7,1,1) , 0
*set,forcez1(8,1,1) , 191.82 
*set,forcez1(9,1,1) , 191.82 
*set,forcez1(10,1,1) , 0   
*set,forcez1(11,1,1) , -191.82    
*set,forcez1(12,1,1) , -191.82    
*set,forcez1(13,1,1) , 0 
!the second force -z
*dim,forcez2,table,13,1,1,,,                !ѭغɱforcez2
!time
*set,forcez2(1,0,1) , 0                     !趨ѭغɱʱֵ
*set,forcez2(2,0,1) , 0.0025 
*set,forcez2(3,0,1) , 0.0075
*set,forcez2(4,0,1) , 0.01  
*set,forcez2(5,0,1) , 0.0125
*set,forcez2(6,0,1) , 0.0175 
*set,forcez2(7,0,1) , 0.02
*set,forcez2(8,0,1) , 0.0225  
*set,forcez2(9,0,1) , 0.0275
*set,forcez2(10,0,1) , 0.03
*set,forcez2(11,0,1) , 0.0325   
*set,forcez2(12,0,1) , 0.0375 
*set,forcez2(13,0,1) , 0.04   
!-z
*set,forcez2(1,1,1) , 0                     !趨ѭغɱֵ
*set,forcez2(2,1,1) , -191.82 
*set,forcez2(3,1,1) , -191.82  
*set,forcez2(4,1,1) , 0 
*set,forcez2(5,1,1) , 191.82 
*set,forcez2(6,1,1) , 191.82 
*set,forcez2(7,1,1) , 0
*set,forcez2(8,1,1) , -191.82 
*set,forcez2(9,1,1) , -191.82 
*set,forcez2(10,1,1) , 0   
*set,forcez2(11,1,1) , 191.82    
*set,forcez2(12,1,1) , 191.82    
*set,forcez2(13,1,1) , 0 
!
timint,on                             	!ʱ
time,0.04                             	!趨ʱ0.04
deltim,0.0005                         	!趨Ӳʱ0.0005
f,zuoduanmian,fz,%forcez1%             	!ʩ+zѭغ
f,youduanmian,fz,%forcez2%             	!ʩ-zѭغ
allsel                                 	!ѡͼԪ 
outres,all,all                          !趨ÿӲȫ 
solve                                 	!
finish                                	!ģ
/post26                               	!ʱ̺ 
nsol,2,2,u,x,                           !ڵ2xλƶΪ2
nsol,3,2,u,y,                           !ڵ2yλƶΪ3
nsol,4,2,u,z,                           !ڵ2zλƶΪ4
nsol,5,2580,u,x,                        !ڵ2580xλƶΪ5 
nsol,6,2580,u,y,                        !ڵ2580yλƶΪ6
nsol,7,2580,u,z,                        !ڵ2580zλƶΪ7 
plvar,2,3,4,5,6,7,                      !ͼʾ2~7
prvar,2                               	!ͼʾ2
finish                                	!ʱ̺ģ
/post1                                	!ͨú
set,,,1,,0.01,,,                        !ʱ0.01ʱĽ
pldisp,2                              	!ƽṹͼ
plnsol,s,eqv                           	!ͼʾƽЧӦ
finish                                	!ͨúģ
/exit,all                               !ݲ˳
