!Ӧʵ3-6
!===غ===
FINI
/CLE

/FILNAME,EX3-6
/TITLE,EX3-6

/prep7
et,1,45				!solid45ʵ嵥Ԫ
et,2,42				!plane42ƽ浥Ԫ
mp,ex,1,26e9			!C20
mp,prxy,1,0.2			!ɱ
mp,dens,1,2449			!ܶ

!ģ
k,1
k,2,,-9
k,3,-6.35,-9
k,4,-3.35
a,1,2,3,4

type,2
aesize,all,0.5
amap,1,1,2,3,4			!ӳ仮浥Ԫ

type,1
extopt,ESIZE,15			!ָ췽򻮷Ϊ15
extopt,ACLEAR,1			!Ϻ浥Ԫ
vext,all,,,,,8.2		!z8.2m

!ʩӲѹ
et,500,shell63			!63ŵԪʩغ
!ѡȡ
type,500
asel,s,loc,x,0			!ѡʩغɵ
amesh,all 
esla,s				!ѡеԪ
nsla,s,1			!ѡнڵ	

*get,enmax,elem,,num,max	!õѡĵԪ
*get,enmin,elem,,num,min	!õѡСĵԪ

dofsel,s,fx,fy,fz		!ѡ񽫽вɶ
fcum,add 			!ʩӷʽΪ"ۼ",ȱʡ""

*do,i,enmin,enmax		!ʼԵenminŵenmaxŵԪѭ
  *if,esel(i),eq,1,then		!ѡ䣬жϵǰԪǷѡ
    *get,ae,elem,i,area 	!õԪʵ
    xe=centrx(i) 		!ԪĵX(ѹֵ)
    ye=centry(i) 		!ԪĵY(ѹֵ)
    ze=centrz(i) 		!ԪĵZ(ѹֵ)
    !ѹ仯Ĺʽ
    p_e=1.71e4*(ye-1)		!ѹʽ
    f_tot=p_e*ae		!Ԫϵĺ
    esel,s,elem,,i		!ѡiԪ
    nsle,s,corner		!ѡԪǵϵĽڵ
    *get,nn,node,,count		!õǰѡнڵĸ
    f_n=f_tot/nn		!õÿڵʩӵ
    *do,j,1,nn			!ѭԵiԪϵÿڵʩӼ
      f,nelem(i,j),fx,f_n*cos(16.5*3.14/180) !ѹ÷ΪX
      f,nelem(i,j),fy,f_n*sin(16.5*3.14/180) !ѹ÷ΪY
    *enddo
  *endif
  esla,s			!ѡеԪ
*enddo

aclear,all			!ϵshell63Ԫ
fcum,repl 			!ʩӷʽԭΪȱʡ""
dofsel,all			!ԭΪɶ
allsel 

!Լ
asel,s,loc,y,-9
da,all,all

/solu
alls
solve