武汉大学测绘学院工测带缓和曲线的圆曲线

#coding:utf-8
__author__ = 'zy'
import math
def JIAODU(x1,y1,x2,y2):jd=math.atan((y2-y1)/(x2-x1))jd=jd/math.pi*180if (x2-x1)<0:jd=jd+180elif (y2-y1)>=40:jd=jdelse:jd=jd+360#把角度化为弧度#jd=jd/180*math.pireturn jd#返回了角度
#由于Python2中引用不到希腊数字，所以用a代表α，B0，PHi代替 Φ#曲线线路坐标计算
#设置ZH点为坐标原点
#坐标方位角的计算g=open(r'g:\zuobiao.txt')
lines=g.readlines()
g.close()Xjd=float(lines[2])
Yjd=float(lines[3])
Xzh=float(lines[0])
Yzh=float(lines[1])
Xhz=float(lines[4])
Yhz=float(lines[5])
a=-JIAODU(Xzh,Yzh,Xjd,Yjd)+JIAODU(Xjd,Yjd,Xhz,Yhz)+180print '请输入缓和曲线长度'
Ls=float(lines[7])
R=float(lines[6])#a=float(raw_input('曲线偏角'))#输入角度转化为弧度
a=a/180*math.pi
#缓和曲线参数m,P,β0m=Ls/2-pow(Ls,3)/240/pow(R,2)
P=pow(Ls,2)/24/R
B0=Ls/2/R*180/math.pi#综合曲线要素
Th=m+(R+P)*math.tan(a/2)
Lh=math.pi*R/180*(a-2*B0)+2*Ls
Eh=(R+P)/math.cos(a/2)-R
q=2*Th-Lh#主点里程计算,交点里程由设计者提供
#print '交点里程'
Kjd=Th
Kzh=Kjd-Th
Khy=Kzh+Ls
Kqz=Kzh+Lh/2
Lt=(a-2*B0)*R
Kyh=Khy+Lt
Khz=Kyh+Ls
#检核
#Kjd=Kqz+q/2#给出坐标，ZH，JD,HZ,以ZH为
x=[]
y=[]
X=[]
Y=[]n=float(raw_input('输入桩位间隔'))print Khz
azh=JIAODU(Xzh,Yzh,Xjd,Yjd)/180*math.pi
ahz=JIAODU(Xjd,Yjd,Xhz,Yhz)/180*math.pidef ZHHY(i,j):Ki=iLi=Ki-Kzhx.append(Li-pow(Li,5)/40/pow(R,2)/pow(Ls,2))y.append(pow(Li,3)/6/R/Ls)#计算曲线线路坐标#Xzh+x[j]*math.cos(result[0])-y[j]*math.sin(result[0])#Yzh+x[j]*math.sin(result[0])+y[j]*math.cos(result[0])X.append(Xzh+x[j]*math.cos(azh)-y[j]*math.sin(azh))Y.append(Yzh+x[j]*math.sin(azh)+y[j]*math.cos(azh))
def HYYH(i,j):Ki=iLi=Ki-KzhPHi=B0+(Li-Ls)/R*180/math.pix.append(m+R*math.sin(PHi))y.append(P+R*(1-math.cos(PHi)))X.append(Xzh+x[j]*math.cos(azh)-y[j]*math.sin(azh))Y.append(Yzh+x[j]*math.sin(azh)+y[j]*math.cos(azh))def YHHZ(i,j):Ki=iLi=Khz-Kix.append(Li-pow(Li,5)/40/pow(R,2)/pow(Ls,2))y.append(pow(Li,3)/6/R/Ls)#Xhz+x[j]*math.cos(result[1])+y[j]*math.sin(result[1])#Yhz+x[j]*math.sin(result[1])-y[j]*math.cos(result[1])X.append(Xhz+x[j]*math.cos(ahz)+y[j]*math.sin(ahz))Y.append(Yhz+x[j]*math.sin(ahz)-y[j]*math.cos(ahz))i=0
j=0
while i<Khz:i=i+nif i<= Khy:ZHHY(i,j)elif i<=Kyh:HYYH(i,j)else :YHHZ(i,j)j=j+1f = open("F:/1.txt", "w+")
f.write("里程 \t X坐标 \t Y坐标\n")
for i in range(0, len(X)):i=0i=i+nf.write("里程"+i+"\t" + X[i] +"\t"+Y[i]+ "\n")
f.close()

武汉大学测绘学院工测带缓和曲线的圆曲线相关推荐

武汉大学测绘学院19级导航工程第三学期专业实习
武汉大学测绘学院19级导航工程第三学期专业实习前言实习任务实习软件自编程序更多前言武汉大学测绘学院导航工程专业实习利用GNSS接收机进行不同观测环境下的静态数据和动态数据的采集,利用C ...
武汉大学测绘学院《导航学》课程复习资料
武汉大学测绘学院<导航学>课程复习资料教学计划: 第一章绪论 2 课时第二章导航学基础知识 6 课时第三章导航图及其投影 7 课时第四章导航基本原理及观测模型 9 课时第 ...
武汉大学测绘学院linux考试试题,大地测量参考框架_武汉大学测绘学院.ppt
大地测量参考框架_武汉大学测绘学院大地测量参考框架(郭际明) 1 概论大地基准(Geodetic Datum):用以代表地球形体的旋转椭球,建立大地基准就是求定旋转椭球的参数及其定向(椭球旋转轴平 ...
武汉大学计算机学院2019考研复试,2019年武汉大学硕士研究生复试及录取名单汇总...
原标题:2019年武汉大学硕士研究生复试及录取名单汇总考生可以通过录取名单了解到很多重要的信息,例如复试比例,进复试最低分,复试录取成绩,录取总评成绩等重要信息.以下是我们整理收集到的各学院复试录取 ...
lisp 非对称缓和曲线_利用AutoLISP语言绘制带缓和曲线的铁路曲线平面图
利用 AutoLISP 语言绘制带缓和曲线的铁路曲线平面图 [ 摘要 ] 通过 AutoLISP 编程方法的剖析 , 探讨了基于 Visual LISP(VLISP) 编程软件环境下 , 利用 Aut ...
武汉大学计算机学院编程能力,2014武汉大学计算机学院国家多媒体软件工程技术研究中心复试经验总结...
本帖最后由 cs_focus 于 2019-5-28 10:10 编辑今天是2014年4月17日,离复试结束已经过去了近半个月时间.我报考的是武汉大学计算机学院的通信专业,所在的实验室叫做国家多媒体 ...
武汉大学计算机学院csc,CSC公派出国留学经验交流讲座
2018年10月21日下午14点30分,由武汉大学水利水电学院主办,IAHR武汉大学学生分会承办,水利水电学院研究生会博士生分会协办的CSC公派留学经验交流讲座于工学部8教105室成功召开.本次会议邀 ...
计算机专业期末考试是编程序,武汉大学计算机学院C语言期末考试
武汉大学计算机学院 2008-2009学年度第2学期2008级 <高级语言程序设计>期末考试试卷A卷答案姓名:学号:专业: 说明:开卷考试,答案请全部写在答题纸上,写在试卷上无效. 未经 ...
武汉大学计算机学院参考书目,2020武汉大学计算机与软件工程考研初试科目、参考书目及复试详情...
原标题:2020武汉大学计算机与软件工程考研初试科目.参考书目及复试详情本文将由新祥旭徐老师全方位的对武汉大学计算机专业考研进行解析,主要有以下几个板块:学院介绍,专业情况介绍,2019录取情况分析 ...
武汉大学计算机学院2015级,武汉大学研究生课程-数据挖掘-2015级研究生试题.doc...
文档介绍: 1 武汉大学计算机学院 2015 级研究生"数据仓库和数据挖掘"课程期末考试试题要求: 所有的题目的解答均写在答题纸上,需写清楚题目的序号.每张答题纸都要写上姓名和学号 ...

武汉大学测绘学院工测带缓和曲线的圆曲线

武汉大学测绘学院工测带缓和曲线的圆曲线相关推荐

最新文章

热门文章