Finish
/clear ,start
/filname ,Arch-dam 123
/title ,ellipse-line model
/plopts,date,0
/TRIAD,LBOT
*AFUN,DEG
NUM=8
*DIM,ELEVATION,ARRAY, NUM
*DIM ,Tc ,ARRAY, NUM
*DIM,ARCH_A,ARRAY, NUM,2
*DIM,ARCH_B,ARRYA,NUM,2
*DIM,ARCH_ANGLE,ARRAY,NUM,2
*DIM ,Ta ,ARRAY, NUM,2
*DIM ,Barch ,ARRAY, NUM
*DIM ,Yu ,ARRAY,NUM
ELEVATION(1)=2460.0,2430.0,2400.0,2360.0,2320.0,2280.0,2240.0,2210.0
H_DAM=ELEVATION(1)-ELEVATION(NUM)
Tc(1)=8.0,16.935,23.167,28.598,32.407,36.413,42.436,49.372
ARCH_A(1,1)=381.568,304.505,251.884,208.573,181.902,156.548,117.188,89.258
ARCH_A(1,2)=346.422,286.802,244.090,205.807,179.223,153.827,119.108,86.627
ARCH_B(1,1)=461.288,362.616,297.409,247.998,221.844,197.526,153.625,94.961
ARCH_B(1,2)=394.568,326.588,279.615,239.267,211.866,183.755,141.274,91.768
ARCH_ANGLE(1,1)=44.656,47.003,46.986,48.506,48.881,43.655,37.600,34.588
ARCH_ANGLE(1,2)=40.934,45.676,49.866,50.505,51.119,48.367,40.995,36.705
Ta(1,1)=8.557,22.952,35.882,50.323,60.905,66.909,67.616,64.237
Ta(1,2)=8.557,25.037,38.268,51.347,59.828,64.386,65.695,64.949
Yu(1)=0.000,-10.887,-18.824,-25.638,-29.175,-30.558,-30.910,-31.174
*VFACT,,,0.5
*VOPER ,Barch ,Yu ,add ,Tc
/PREP7
*DO,ii,1,num
A4=arch_angle(num+1-ii,1)
A5=arch_angle(num+1-ii,2)
A6=ELEVATION(1)-ELEVATION(NUM+1-II)
*do,j,10,1,-1
A1=SQRT(ARCH_B(num+1-ii,1)**2+ ARCH_A(num+1-ii,1)**2*TAN(j*A4/10)**2)
A2=ARCH_A(num+1-ii,1)**2*tan(j*A4/10)
BB=ARCH_B(num+1-ii,1)-ARCH_B(num+1-ii,1)**2/A1+Barch(num+1-ii)
TL=Tc(num+1-ii)+(Ta(num+1-ii,1)-Tc(num+1-ii))*(j/10)**4
PX1=A2/A1+TL*SIN (j*A4/10)/2
PY1=BB-TL*COS(j*A4/10)/2
PX2=A2/A1-TL*SIN (j*A4/10)/2
PY2=BB+TL*COS (j*A4/10)/2
K,(II-1)*100+11-J,PX1,PY1,A6
K,(II-1)*100+32-J,PX2,PY2,A6
*enddo
*do,J,0,-10,-1
A1=SQRT(ARCH_B(num+1-ii,2)**2+ARCH_A(num+1-ii,2)**2*TAN(j*A5/10)**2)
A2=ARCH_A(num+1-ii,2)**2*tan(j*A5/10)
BB=ARCH_B(num+1-ii,2)-ARCH_B(num+1-ii,2)**2/A1+Barch(num+1-ii)
TR=Tc(num+1-ii)+(Ta(num+1-ii,2)-Tc(num+1-ii))*(J/10)**4
PX3=A2/A1+TR*SIN(J*A5/10)/2
PY3=BB-TR*COS(J*A5/10)/2
PX4=A2/A1-TR*SIN(J*A5/10)/2
PY4=BB+TR*COS(J*A5/10)/2
K,(II-1)*100+11-J,PX3,PY3,A6
K,(II-1)*100+32-J,PX4,PY4,A6
*ENDDO
*ENDDO
Ksel ,all
l,1,22
*repeat,8,100,100
l,21,42
*repeat,8,100,100
*do,i,1,8,1
ksel,s,,,(100*(i-1)+1),(100*(i-1)+21)
bsplin,all
allsel,all
ksel,s,,,(100*(i-1)+22),(100*(i-1)+42)
bsplin,all
allsel,all
*enddo
askin,18,20,22,24,26,28,30,32
askin,17,19,21,23,25,27,29,31
askin,1,2,3,4,5,6,7,8
askin,9,10,11,12,13,14,15,16
al,1,9,17,18
al,8,16,31,32
va,1,2,3,4,5,6
ASEL,S,,,4
VEXT,ALL,,,-4*H_DAM
ASEL,S,,,3
VEXT,ALL,,,4*H_DAM
ALLSEL
asel,s,loc,z,250
vext,all,,,,,nint(h_dam/10)*10
asel,s,,,8,13,5
asel,a,,,23,27,4
asel,a,,,21
VEXT,ALL,,,,-4*H_DAM
asel,s,,,10,25,5
asel,a,,,19,29,10
VEXT,ALL,,,,4*H_DAM |