本帖最后由 VibInfo 于 2016-4-15 15:03 编辑
*elseif,tpos1,le,(L1+L2),then ! 煤水车进桥
c1 = xlth1 - tpos1
*if,tele1,eq,tele2,then !火车头在同一个单元上
c2 = xlth2 - tpos2
c3 = elth2 - c2
qq = (q1*L1+q2*c3)/elth1
sfbeam,eldk(tele1),1,pres,qq,qq,,,,,
*do,j,1,tele1-1
sfbeam,eldk(j),1,pres,q2,q2,,,,,
*enddo
*else
dtele = tele1 - tele2
*if,dtele,eq,1, then
c2 = xlth2 - tpos2
c3 = elth2 - c2
qq0 = q1*(elth1-c1)/elth1
qq = (q1*c2+q2*c3)/elth2
sfbeam,eldk(tele1),1,pres,qq0,qq0,,,,,
sfbeam,eldk(tele2),1,pres,qq,qq,,,,,
*do,j,1,tele2-1
sfbeam,eldk(j),1,pres,q2,q2,,,,,
*enddo
*else
c2 = xlth2 - tpos2
c3 = elth2 - c2
qq0 = q1*(elth1-c1)/elth1
sfbeam,eldk(tele1),1,pres,qq0,qq0,,,,,
*do,j,tele2+1,tele1-1
sfbeam,eldk(j),1,pres,q1,q1,,,,,
*enddo
qq = (q1*c2+q2*c3)/elth2
sfbeam,eldk(tele2),1,pres,qq,qq,,,,,
*do,j,1,tele2-1
sfbeam,eldk(j),1,pres,q2,q2,,,,,
*enddo
*endif
*endif
*elseif,tpos1,le,qth,then ! 旅客列车进桥
c1 = xlth1 - tpos1
*if,tele1,eq,tele2, then
c2 = xlth2 - tpos2
c3 = elth2 - c2
qq = (q1*L1+q2*c3)/elth1
sfbeam,eldk(tele1),1,pres,qq,qq,,,,,
c4 = xlth3 - tpos3
c5 = elth3 - c4
dtele1 = tele2 - tele3
*if,dtele1,eq,1,then
qq1 = (q2*c4+q3*c5)/elth3
sfbeam,eldk(tele3),1,pres,qq1,qq1,,,,,
*do,j,1,tele3-1
sfbeam,eldk(j),1,pres,q3,q3,,,,,
*enddo
*else
*do,j,tele3+1,tele2-1
sfbeam,eldk(j),1,pres,q2,q2,,,,,
*enddo
qq1 = (q2*c4+q3*c5)/(elth3)
sfbeam,eldk(tele3),1,pres,qq1,qq1,,,,,
*do,j,1,tele3-1
sfbeam,eldk(j),1,pres,q3,q3,,,,,
*enddo
*endif
*else
dtele1 = tele1 -tele2
*if,dtele1,eq,1,then
c2 = xlth2 - tpos2
c3 = elth2 -c2
qq0 = q1*(elth1-c1)/elth1
sfbeam,eldk(tele1),1,pres,qq0,qq0,,,,,
c4 = xlth3- tpos3
c5 = elth3 - c4
dtele2 = tele2 - tele3
*if,dtele2,eq,0,then
qq0 = (q1*c2+q2*L2+q3*c3)/elth2
sfbeam,eldk(tele2),1,pres,qq0,qq0,,,,,
*do,j,1,tele2-1
sfbeam,eldk(j),1,pres,q3,q3,,,,,
*enddo
*elseif,dtele2,eq,1,then
qq0 = (q1*c2+q2*c3)/elth2
sfbeam,eldk(tele2),1,pres,qq0,qq0,,,,,
qq1 = (q2*c4+c3*c5)/elth3
sfbeam,eldk(tele3),1,pres,qq1,qq1,,,,,
*do,j,1,tele3-1
sfbeam,eldk(j),1,pres,q3,q3,,,,,
*else
qq0 = (q1*c2+q2*c3)/elth2
sfbeam,eldk(tele2),1,pres,qq0,qq0,,,,,
*do,j,tele3+1,tele2-1
sfbeam,eldk(j),1,pres,q2,q2,,,,,
*enddo
qq1 = (q2*c4+q3*c5)/(elth3)
sfbeam,eldk(tele3),1,pres,qq1,qq1,,,,,
*do,j,1,tele3-1
sfbeam,eldk(j),1,pres,q3,q3,,,,,
*enddo
*endif
*else
c2 = xlth2 - tele2
c3 = elth2 - c2
qq0 = q1*(elth1-c1)/elth1
sfbeam,eldk(tele1),1,pres,qq0,qq0,,,,,
*do,j,tele2+1,tele1-1,
sfbeam,eldk(j),1,pres,q1,q1,,,,,
*enddo
qq = (q1*c2+q2*c3)/(elth2)
sfbeam,eldk(tele2),1,pres,qq,qq,,,,,
dtele2 = tele2 - tele3
*if,dtele2,eq,1,then
qq1 = (q2*c4+q3*c5)/(elth3)
sfbeam,eldk(tele3),1,pres,qq1,qq1,,,,,
*do,j,1,tele3-1
sfbeam,eldk(j),1,pres,q3,q3,,,,,
*enddo
*else
*do,j,tele3+1,tele2-1
sfbeam,eldk(j),1,pres,q2,q2,,,,,
*enddo
qq1 = (q2*c4+q3*c5)/(elth3)
sfbeam,eldk(tele3),1,pres,qq1,qq1,,,,,
*do,j,1,tele3-1
sfbeam,eldk(j),1,pres,q3,q3,,,,,
*enddo
*endif
*endif
*endif
|