线元法参数计算程序(自编)
在网上寻找计算线元法需要的参数的5800程序找了很久没找到，后来索性自己更具《详解用线元法计算公路中线坐标》中的公式自己编了个5800程序。
为了区分字母O和数字0（字母O染红了的）
Fix 3：14→DimZ
Lbl 0
"1(YP),-1(ZP)="?→Z[14]:"QD(KM)="?A（右偏输入1，左偏输入-1）"LS="?B（线元长度）
"QD(QL)="?C（起点曲率，半径的倒数，直线输入1/1045趋近于0）"ZD(QL)="?D（终点曲率）
"QD(X)="?E（起点X坐标）
"ZD(Y)="?F（起点Y坐标）
"QD(FWJ)="?G（起点方位角）
0.1184634425→R:0.2393143352→Z[11]:0.2844444444→Z[12]:0.0469100770→S:0.2307653449→T: 0.953089923→W:0.769234655→V:0.5→U:0.2393143352→K:0.1184634425→O
Z[14]=1=?Goto 1：Z[14]=-1=?Goto 3
Lbl 1
D-C→Z
G+(D+C)90B÷∏→H
Goto 2
Lbl 3
-D-C→Z
G+(-C-D)90B÷∏→H
Goto 2
Lbl 2
Rcos((G∏+180CSB+90ZBS2)÷∏)→Z[1]:Z[11]cos((G∏+180CTB+90ZBT2)÷∏)→Z[2]:Z[12]cos((G∏+180CBU+90ZBU2)÷∏)→Z[3]:Kcos((G∏+180CBV+90ZBV2)÷∏)→Z[4]:O cos((G∏+180CBW+90ZBW2)÷∏)→Z[5]:Rsin((G∏+180CBS+90ZBS2)÷∏)→Z[6]:Z[11]sin((G∏+180CBT+90ZBT2)÷∏)→Z[7]:Z[12]sin((G∏+180CBU+90ZBU2)÷∏)→Z[8]:Ksin((G∏+180CBV+90ZBV2)÷∏)→Z[9]:O sin((G∏+180CBW+90ZBW2)÷∏)→Z[10]
E+B(Z[1]+ Z[2]+ Z[3]+ Z[4]+ Z[5]) →Z[13]
F+B(Z[6]+ Z[7]+ Z[8]+ Z[9]+ Z[10]) →N
Cls
“ZDX=”:Locate 5,1,Z[13]: Locate 1,2,”ZDY=”: Locate5,2,N: Locate1,3,”F=”: Locate4,3,H o: Locate1,4,”ZDK=”: Locate5,4,(A+B)◢
Cls
Goto 0
ZDX=待求点的X坐标
ZDY=待求点的Y坐标
F=待求点的方位角 ZDK待求点的桩号
线元法线路坐标正反算程序
经苦心钻研，奋战多日，终于编写出了代码短，速度快，精度高，功能全的线路坐标正反算程序，欢迎试用并提出宝贵意见。 功能简介及特点： 1、选用高斯-勒让德公式作计算内核，保证精度，模块化设计，便于扩充功能。 2、线元数据可自动从数据库调用，也可手工输入。 3、可管理多条线路，如里程不在线路或线元范围，将警告里程偏大、偏小。 4、边桩计算设计为导线式递推方式，可用于由一个中桩推出结构物所有角点坐标。 5、反算实现了智能化操作，只需输入线路号（或手工输线元资料）、坐标，不需近似里程，即可自动从起点向后开始试算出里程、位置，如对算出里程、位置表示怀疑，还可以让计算器从终点起再向前试算下一个可能的位置（匝道、回头曲线同一坐标可能会有一个以上结果）。第三次及以后试算才要求输入近似里程。 6、程序代码规范简洁，便于阅读、理解。 完整程序清单: ZFS %正反算主程序 B=.1739274226:C=.5-B: Lbl 1:U"0 ZS 1 FS"=0=>Prog "ZS": ≠>U=1=>Prog"FS":≠>Goto 1
ZS %正算子程序 {K}:Prog"ZZ":I=0:{I}:I"L"≠0=>"Prog"WY":≠>Prog"ZB" FS %反算子程序 {KVW}:V"XC"W"YC":Lbl 2:Prog "ZZ":I=V-S:J=W-T:Pol(I,J: J=J-F:K=K+Rec(I,J:AbsI<1m=>Prog"WZ":≠>Goto 2Δ M=0:{M}:M"0 NEXT"=0=>U=U+1:Goto 2:≠>U=1 ZZ %高斯法中桩子程序(4节点) Prog"XL":M=K-L:O=(P-R)÷2PQR: D=.0694318442:E=.3300094782:F=1:G=1-E:H=1-D: I=5:Lbl 1:C[I]=A+MrC[I](1÷P+OMC[I]:Dsz I:Goto 1: S=X+M(BcosD+CcosE+CcosG+BcosH: T=Y+M(BsinD+CsinE+CsinG+BsinH WY %外移点计算子程序 Lbl 1:J=90:{J}:J=F+J"<":F=J:S=S+Rec(I,J:T=T+J: Prog"ZB":I=0:{I}:I"L"≠0=>Goto 1 WZ %位置显示子程序 "KJ":K:Pause 1:J◢ ZB %坐标显示子程序 "XY":S:Pause 1:T◢ YC %异常处理子程序 U=1=>K=L:U=2Δ U=3=>K=M:U=4Δ
中海达七参数计算
HI-RTK道路版简易操作流程 一、架设基准站： 选择视野开阔且地势较高的地方架设基站，基站附近不应有高楼或成片密林（卫星接收不好）、大面积水塘（多路径效应严重）、高压输电线或变压器（有干扰）。基站一般架设在未知点上，后面的说明均征对这种情况。（此种情况下基站无需对中整平） 二、新建项目： 打开HI-RTK道路软件，进入“项目”，选定Unnamed，“套用”，输入项目名称后确认，（选择‘套用’而不是‘新建’的目的是为了使建立的项目里面不含任何人为参数） 然后：项目信息---坐标系统---（将坐标系统名称改为“中国-‘项目名’ ”）并确认每个选项的原始参数是否正确，需要改动的地方请改正---保存---退出---（弹出“是否更新点库”）是。 三、设置基准站： 1. GPS---接收机信息---连接GPS---连接---搜索（接收机）---（搜索到仪器后）停止---（选择仪器号）连接。 2.接收机信息---基准站设置---平滑---（采集10秒后）确认---（查看并确认另外两个选项内容是否正确）---确定---断开蓝牙连接。 四、移动站设置： 1. GPS---接收机信息---连接GPS---连接---搜索（接收机）---（搜索到仪器后）停止---（选择仪器号）连接。 2.接收机信息---移动站设置---（确认每个选项内容）---确定。 五、采集已知点并求取参数： 1.采集已知点：已知点采集的时候建议采用“平滑采集”，按钮为工具栏倒数第二个按钮。（最少采集两个已知点，计算七参数时至少需要三个已知点）
2.输入已知点理论坐标到点库：碎步测量---控制点库---添加（工具栏第一个按钮）---（输入点名，X,Y,H后确认）。 3.参数计算： （主界面）参数---坐标系统---参数计算---（选择计算类型，采集两个已知点时用‘四参数+高程拟合’）---添加---（‘源点’为外业采集的点，‘目标’为输入的已知点，按钮为调用点库信息。）---保存---（继续添加）---解算---运用---（坐标系统）保存---（是否覆盖）确定---确定---（更新点库）是---退出。（请确认点对配对正确） 4.进行碎步采集或者放样。 5.数据导出：从项目或者测量界面进入“记录点库”，点击工具栏最后一个按钮，输入导出文件名、选择导出文件类型后确定，然后手簿连接电脑拷贝出对应数据即可。 这个是最全面，最权威的说明书了。
平面三角形单元有限元程序设计
. 一、题目 如图1所示，一个厚度均匀的三角形薄板，在顶点作用沿板厚方向均匀分布的竖向载荷。已知：P=150N/m ，E=200GPa ，=0.25，t=0.1m ，忽略自重。试计算薄板的位移及应力分布。 要求： 1. 编写有限元计算机程序，计算节点位移及单元应力。（划分三角形 单元，单元数不得少于30个）； 2. 采用有限元软件分析该问题（有限元软件网格与程序设计网格必 须一致），详细给出有限元软件每一步的操作过程，并将结果与程序计算结果进行对比（任选取三个点，对比位移值）； 3. 提交程序编写过程的详细报告及计算机程序； 4. 所有同学参加答辩，并演示有限元计算程序。 有限元法中三节点三角形分析结构的步骤如下： 1）整理原始数据，如材料性质、荷载条件、约束条件等，离散结构并进行单元编码、结点编码、结点位移编码、选取坐标系。 2）单元分析，建立单元刚度矩阵。 3）整体分析，建立总刚矩阵。 4）建立整体结构的等效节点荷载和总荷载矩阵 5）边界条件处理。 6）解方程，求出节点位移。 7）求出各单元的单元应力。 8）计算结果整理。 一、程序设计 网格划分 如图，将薄板如图划分为6行，并建立坐标系，则
刚度矩阵的集成 建立与总刚度矩阵等维数的空矩阵，已变单元刚度矩阵的集成。 由单元分析已知节点、单元的排布规律，继而通过循环计算求得每个单元对应的节点序号。 通过循环逐个计算：（1）每个单元对应2种单元刚度矩阵中的哪一种； （2）该单元对应总刚度矩阵的那几行哪几列 （3）将该单元的单元刚度矩阵加入总刚度矩阵的对应行列 循环又分为3层循环：（1）最外层：逐行计算 （2）中间层：该行逐个计算 （3）最里层：区分为第 奇/偶 数个计算 单元刚度的集成：[ ][][][][][]' '''''215656665656266256561661e Z e e e Z e Z e e e e k k k K k k k k k k +?++=? =?==?==?=?????? 边界约束的处理：划0置1法 X Y P X Y P
线元法简介
线元法万能曲线正反算简介 我的线元法是把线形分为直线和曲线，直线就不用说了，起止点桩号，坐标和方位角就可以算了；曲线最基本的组合：是由一段缓和曲线+一段圆曲线组成，任意复杂的曲线都可以分解成缓和曲线+圆曲线或者其中之一就可以。 分析最复杂的曲线可以看到: 一般复杂线形由Ls1 ,R1,Ls2, R2组成，相邻的Ls1+R1,一般满足A*A=Ls1*R1,这就是一个线元法单元，即使不满足也可以作为一个线元： 当Ls1= Ls2，且R1= R2时，为单曲线 当Ls1≠ Ls2，或者R1≠R2时，为复合曲线 当Ls1= Ls2=0时，线性为圆曲线， 当圆曲线长度为0时，线性为缓和曲线+缓和曲线， 当A*A≠Ls1*R1时，为卵形曲线，需要计算虚拟起点坐标 综合以上线形，本程序正反算计算全部可以处理。结合目前流行的线元法，本程序也可以，分为缓和曲线和圆曲线录入，方法是一样的，所不同的是起点要注意，复杂曲线，是两边向中间定义数据库，缓和曲线永远是ZH点或HZ点为起点。 曲线要素说明（有9个）： 1、起点桩号：（一般为ZH点或HZ点，或ZY点或YZ点，或者卵形公切点GQ） 2~3、起点坐标：（X，Y） 4、起点方位角：FWJ 114°15′24.33″写成：114.152433 5、线性特征：直线，左偏，右偏；三个选一个 6、终点桩号：如果起点为ZH点，终点一边为YH点，QZ点，HY点，都可以，一般为YH点，缓和曲线+圆曲线。如果缓和曲线Ls=0，就是YZ点；大小不一定按路线顺序，如果起点为HZ点，终点根据缓和曲线+圆曲线的特点，和上个线元对接上就可以了。 7、缓和曲线长度Ls： 8、圆曲线半径R: 9、回旋参数A: 一般满足A*A=Ls1*R1，不满足条件的是卵形曲线。 可以处理任意数量断链。 操作流程：1、先编辑线元数据，保存后推出。 2、如果有线元断链的输以下线元断链数据 3、打开线元万能曲线计算单点计算就可以了。 目前，已有一个例子文件在里面，在安装文件目录下“ \dmfx4.0\demo\左线”，有个CAD文件，里面有校核数据，可以看到本软件处理的逐桩表和要素表，可以验证软件的数据，任意数据坐标反算可以得到桩号和距中，任意输入桩号和距中可以正算得坐标。 授权版用户，可以通过运行交点文件编辑，保存后，退出；打开线元法数据编辑，浏览正在使用的主项目文件，就可以看到一个线元数据，点击这个文件确定，保存退出。就完成交点法数据转换线元法数据过程。
线元法万能坐标计算程序
线元法万能坐标计算程序(适用于CASIO fx-9750GⅡ计算器) 论文https://www.360docs.net/doc/a03098981.html,/：本论文仅供学习交流使用，本站仅作合理转载，原作者可来邮要求删除论 文。 摘要：我国公路建设事业正处于一个高速发展的时期，在公路工程施工过程中，施工技术人员经常要使用全站仪、水准仪进行施工放样、高程测量，在测量过程中，手工计算速度慢，失误率高，工作效率极低。利用CASIO fx-9750GⅡ编程函数计算器强大的内存（可诸存63000个字符）和编程功能，编写各种计算程序，能够在2秒钟内计算出施工放样、桩点坐标等施工过程中的各项数据资料，同时也使我们有更多的时间去挑战更富有创造性的工作。 关键词：坐标放线线元测量程序 1、前言 本程序采用Gauss-Legendre（高斯-勒让德）五节点公式作内核,计算速度（太约2秒）适中，计算精度很高。在此之前，本人曾用过以下公式作内核：①积分公式simpson法②双重循环复化高斯2节点③高斯-勒让德3节点④求和公式复化simpson法⑤双重循环复化simpson法⑥高斯-勒让德4节点，⑦高斯-勒让德5节点，经过测试③计算最快，⑦代码稍长但计算速度只比③⑥稍慢，精度最高，可满足线元长小于1/2πD 的所有线形的精度要求。⑦作内核分别计算圆曲线长1/4πD、1/2πD、3/4πD、πD处的精度,1/4πD时偏差为0.001mm,1/2πD时偏差为0.55m m,3/4πD时偏差为31.63mm，πD时偏差为968mm，偏差按半径倍数增大，如线元长大于1/2πD（1/2圆周长）时，可将其拆分二个或多个线元单位，以确计算保精度。 2、程序特点 事先将所有的平曲线交点的线元要素诸存到计算器内，测量时只输桩号、边距等程序会自动寻找各类要素，一气呵成地完成施工测量任务，中途不需人工转换各类要素数据，本程序可诸存几百条线路的要素数据，计算时可按需选择线路编号进行测量。测量时不需查阅及携带图纸，仅一台CASIO fx-9750GⅡ编程函数计算器即可。 本程序含一个主程序:3XYF，五个子程序:GL（公式内核）、QD（线路选择）、XL（线路要素判断）、GF（坐标反算）、File 1 （要素存放的串列工作簿）。可以根据曲线段——直线、圆曲线、缓和曲线（完整或非完整型）的线元要素（起点坐标、起点里程、起点切线方位角、终点里程、起点曲率半径、止点曲率半径）及里程边距或坐标，对该线元段范围内任意里程中边桩坐标进行正反算。 3、计算公式及原理 如图：BC 间为一曲线元，曲线元上任一点的曲率随至B 点的弧长作线性变化。设起点B 的曲率为KA ，终点C 的曲率为KB ，R 为曲线半径。±表示曲线元的偏向，当曲线元左偏时取负号，当曲线元右偏时取正号,直线段以1的45次方代替（即半径无穷大）。 式中：αΑ＝起始方位角l ＝p 点到B的距离lS=曲线总长αp＝p 点切线方位角 R1=R5=0.118463442528095 ,R2 = R4 = 0.239314335249683 , R3 = 0.28444444444444 V1=1-V5= 0.046910070 ,V 2= 1-V4 = 1 0.2307653449 V3= 0.5 利用上面公式及CASIO fx-9750GⅡ编程函数计算器可编写下列计算程序。 4、程序清单 （1）、3XYF（主程序） "1→XY2→FS"?→V:V=1=>Goto 1:V=2=>Goto 2↙（选择计算功能） Lbl 1:File 1:”XLn”?→S:Prog “QD”↙（选择线路）
交点法线元法坐标计算
3、交点法、线元法坐标计算 坐标计算是根据图纸中“直线及曲线转角一览表”提供的数据计算道路中桩坐标，然后和图纸提供的“逐桩坐标表”比对，如果一样则说明输入平曲线参数输入正确，可以计算边桩坐标和其他结构物坐标了；如果中桩坐标不一样，一般是平曲线参数输入有误，需要重新检查输入，另一种结果是图纸有错，这种情况少见，但不代表没有。“直线及曲线转角一览表”和“逐桩坐标表”见附件1、附件2。 线元法是以路线的起点坐标、方位角、起终点桩号等节点元素来计算出要求的坐标；交点法是以路线的交点要素和路线的主要要素来求得坐标。 ①交点法 交点：路线的转折点，路线改变方向是相邻两直线的延长线相交的点。用JD表示， 有些图 纸上用 IP表示。 看下图： 交 点是针对曲线的（包含圆曲线和缓和曲线），一段曲线就有一个交点。交点参数有：坐标（X,Y）、交点桩号、转角值、圆曲线半径R、缓和曲线长度。 教学提供软件（轻松测量、双心软件、测量工具）交点法曲线要素输入说明： 1、QD起点坐标： 起点坐标必须在直线段上，或填写前一交点的坐标。
2、JD交点曲线要素： （1）交点桩号 （2）交点坐标（X,Y） （3）曲线半径R 始点的话，起始里程有时候需要校正，当然，并不是每个图纸给出的起点里程都需要校正，大多数图纸的起点里程已经被设计院校正过，我们输入平曲线的时候需要验证一下。如果我们按照图纸给出的起点里程输入，发现后面的交点里程都和图纸相差一个相同的值，这就表明我们输入的起点里程需要校正。 起始点里程正常输入，第二、三个交点输入完成后，检查第二个交点的切线长和交点
里程是否和图纸一样，如果切线长正确，交点里程不正确，说明起点里程需要校正，将第二个交点的里程与正确里程的差值，应用到起点里程中，从而使第二个交点里程和后面交点的里程与图纸吻合。 注意：交点法计算坐标适用的平曲线为对称或不对称缓和曲线、圆曲线。对于非普通的三单元曲线，交点法不适用。非普通的三单曲线例如下页的JD18及JD19处的平曲线， 的输入是否正确，有的图纸给的方位角数据较少，需要每隔几个线元才能检验方位角。
转坐标系详细步骤
转坐标系详细步骤
————————————————————————————————作者：————————————————————————————————日期: ?
“北京54坐标系”转“西安80坐标系”一、数据说明 北京54坐标系和西安8０坐标系之间的转换其实是两种不同的椭球参数之间的转换,一般而言比较严密的是用七参数布尔莎模型，即X平移,Y平移,Ｚ平移,X旋转(ＷX)，Y旋转（WY)，Ｚ旋转（WY)，尺度变化（DＭ)。若得七参数就需要在一个地区提供３个以上的公共点坐标对（即北京54坐标下x、y、z和西安80坐标系下ｘ、y、ｚ），可以向地方测绘局获取。 下面讲述利用已知的3个以上（本例采用4个点计算)的公共点计算七参数方法转换： 二、利用４个已知公共点计算转换七参数 1、数据准备 （１）将已知５4、80坐标系直角坐标拷贝到文本文档,其排列格式如下（图１、图2)：不加带号。 图1 5４直角坐标 图2 ８０直角坐标 (２）将已知５4、80坐标系直角坐标利用MAPＧIS“投影变换”转换为经纬度坐标,且坐标单位为“秒”,这样计算出的参数用来转换为80坐标系时更精确。具体操作步骤如下： 1)启动MＡＰGIS下“投影变换模块”，点击“投影变换”下“用户文件投影转换”弹出“用户数据点文件投影转换”对话框,如图3； 2)点击“打开文件”,选择已准备的“54直角坐标.ｔxｔ”文本文档,打开后选择“按指定分隔符”后弹出的对话框点击确定激活“设置分隔符”选项,点击“设置分隔符”,其设置方式为：①“Ｔａｌ键”、“空格”两个选
图3 图4
项前画勾,②修改“属性名称所在行”，点击其下拉箭头选择“无”字下面一组数据，③将“属性名称”修改为x、y，④“数据类型”修改为“5双精度”,⑤“小数位”修改为“5”或其他均可,但最好至少为“2”，其设置与最终转换出坐标的小数位数相关。设置完成后点击“确定”。如图４。 3)设置“用户投影参数”及“结果投影参数”其设置方式如图５、图6。注意：投影中心点经度一定要输入,如经度为10５°,其格式为10500０0,“用户投影参数”为“投影平面直角坐标”;“结果投影参数”为“地理坐标系”，且“比例尺分母”为“1”,“坐标单位”为“妙”,“投影中心点经度”要输入。二者“椭球参数”均为“54坐标系”。 图５用户投影参数 图6 结果投影参数 ４)以上参数设置完成后点击“投影变换”——“写到文件”,弹出对话框如图7 ，先新建“54经纬度坐标.txt”，选中后点击保存,选择替换。 5)按照上述1）—4)步骤将已知的8０直角坐标转换为以“秒”为单位的经纬度坐标。注意:在“用户投影参数”及“结果投影参数”设置时,二者“椭球参数”均为“80坐标系”,其他参数同上。 转换后的5４和80坐标系以“秒”为单位的经纬度坐标如下：图7、图8。坐标中小数点前为“６位数”的是“经度”，小数点前为“5位数”的是“纬度”。 图7 54经纬度坐标图8 80经纬度坐标
公路坐标计算公式
一、缓和曲线上的点坐标计算 已知：①缓和曲线上任一点离ZH点的长度：l ②圆曲线的半径：R ③缓和曲线的长度：l0 ④转向角系数：K(1或－1) ⑤过ZH点的切线方位角：α ⑥点ZH的坐标：x Z，y Z 计算过程： 说明：当曲线为左转向时，K=1，为右转向时，K=-1，公式中n的取值如下： 当计算第二缓和曲线上的点坐标时，则： l为到点HZ的长度
α为过点HZ的切线方位角再加上180°K值与计算第一缓和曲线时相反 x Z，y Z为点HZ的坐标 切线角计算公式： 二、圆曲线上的点坐标计算 已知：①圆曲线上任一点离ZH点的长度：l ②圆曲线的半径：R ③缓和曲线的长度：l0 ④转向角系数：K(1或－1) ⑤过ZH点的切线方位角：α ⑥点ZH的坐标：x Z，y Z 计算过程：
说明：当曲线为左转向时，K=1，为右转向时，K=-1，公式中n的取值如下： 当只知道HZ点的坐标时，则： l为到点HZ的长度 α为过点HZ的切线方位角再加上180° K值与知道ZH点坐标时相反 x Z，y Z为点HZ的坐标 三、曲线要素计算公式
公式中各符号说明： l——任意点到起点的曲线长度（或缓曲上任意点到缓曲起点的长度）l1——第一缓和曲线长度 l2——第二缓和曲线长度 l0——对应的缓和曲线长度 R——圆曲线半径 R1——曲线起点处的半径 R2——曲线终点处的半径
P1——曲线起点处的曲率 P2——曲线终点处的曲率 α——曲线转角值 四、竖曲线上高程计算 已知：①第一坡度：i1(上坡为“＋”，下坡为“－”) ②第二坡度：i2(上坡为“＋”，下坡为“－”) ③变坡点桩号：S Z ④变坡点高程：H Z ⑤竖曲线的切线长度：T ⑥待求点桩号：S 计算过程： 五、超高缓和过渡段的横坡计算
有限元编程算例(fortran)
有限元编程算例(Fortran) 本程序通过Fortran语言编写，程序在Intel Parallel Studio XE 2013 with VS2013中成功运行，程序为《计算力学》（龙述尧等编）一书中的源程序，仅作研究学习使用，省去了敲写的麻烦。 源程序为： !Page149 COMMON/X1/NJ,NE,NZ,NDD,NPJ,IND,NJ2,EO,UN,GAMA,TE,AE COMMON/X2/JM(100,3),NZC(50),CJZ(100,2),PJ(100,2),B(3,6),D(3,3),S(3,6),TKZ(200,20),EKE(6,6),P(200)
OPEN(5,FILE='DATAIN') !OPEN(6,FILE='DATAOUT',STATUS='NEW') CALL DATA IF(IND.EQ.0)GOTO 10 EO=EO/(1.0-UN*UN) UN=UN/(1.0-UN) 10 CALL TOTSTI CALL LOAD CALL SUPPOR CALL SOLVEQ CALL STRESS PAUSE !STOP END SUBROUTINE DATA COMMON/X1/NJ,NE,NZ,NDD,NPJ,IND,NJ2,EO,UN,GAMA,TE,AE COMMON/X2/JM(100,3),NZC(50),CJZ(100,2),PJ(100,2),B(3,6),D(3,3),S(3,6),TKZ(200,20),EKE(6,6),P(200) READ(5,*)NJ,NE,NZ,NDD,NPJ,IND NJ2=NJ*2 NPJ1=NPJ+1 READ(5,*)EO,UN,GAMA,TE READ(5,*)((JM(I,J),J=1,3),I=1,NE) READ(5,*)((CJZ(I,J),J = 1,2),I=1,NJ) !Page150 READ(5,*)(NZC(I),I=1,NZ) READ(5,*)((PJ(I,J),J=1,2),I=1,NPJ1) WRITE(6,10)(I,(CJZ(I,J),J=1,2),I=1,NJ) 10 FORMA T(4X,2HNO,6X,1HX,6X,1HY/(I6,2X,F7.2,F7.2)) RETURN END SUBROUTINE ELEST(MEO,IASK) COMMON/X1/NJ,NE,NZ,NDD,NPJ,IND,NJ2,EO,UN,GAMA,TE,AE COMMON/X2/JM(100,3),NZC(50),CJZ(100,2),PJ(100,2),B(3,6),D(3,3),S(3,6),TKZ(200,20),EKE(6,6),P(200)
已知七参数输入方法
已知七参数输入方法 我们在测量过程中，常常会遇到要求我们利用已知的七参数进行测量的情况，下面我们来看一下如何在仪器中输入七参数。 1、在主菜单屏幕上选择管理： 七参数：使用严格3D 经典方法产生转换的参数. 该方法使用GPS 测量点(WGS84 椭球 )的直角坐标,并将这些坐标与地 方坐标的直角坐标相比较.通过这种方法,计算出用来将坐标从一个系统转换到另一个系统中平移量,旋转量和尺度因子.经典 3D 转换方法可确定最多7个转换参数(3个平移参数,3个旋转参数,和1个尺度因子). 2、选择坐标系： 3、新建一个坐标系：
4、在名称行里输入一个坐标系统的名字： 5、将光标移至转换一行，点击回车键： 6、点击F2新建：
7、在概要界面输入一个七参数名称，然后点击参数： 8、输入已知的七参数，（也有输入四参数的，即不输旋转参数）： 9、在更多界面下选择莫洛金斯基或布沙-沃尔夫，一般选择后者，然后保存： Molodensky-Badekas ——莫洛金斯基 一种转换模型,其旋转原点是系统A 中公共点的重心. Bursa-Wolf ——布沙-沃尔夫 对系统A 来说,旋转原点为笛卡儿坐标系统原点的转换模型.
10、选择做好参数的转换文件，继续： 11、将光标移至椭球行，回车： 在大地测量中,除非特别定义,椭球是 指椭圆绕短半轴旋转形成的数学图形 (有时也称回转椭球体),两个量定义一 个椭球,它们是长半轴的长度; 扁率 f. The Flattening is one of the quantities to specify an ellipsoid. f = (a-b)/a = 1 - sqrt(1-e2) where: a ... semi-major axis b ... semi-minor axis e ... eccentricity 12、选择要用的椭球（西安-80或北京-54） 如果没有需要的椭球，请点击 SHIFT键，在点击F5键即可调 阅所有椭球 13、将光标移至投影行，回车，然后新建，选择横轴莫卡托，然后输入投影参数，保存： 假定东坐标：为避免坐标出现负值，我 国将坐标原点东坐标规定为500,000 米。 中央子午线：定义地图投影经度的中央 线。是使用在地图投影中的带常数。 带宽：投影带的宽度。 注意：投影参数一定要在开始工作前落 实清楚，否则将影响投影后坐标。
《有限单元法》编程作业
湖南大学 《有限单元法》编程大作业 专业：土木工程 姓名： 学号： 2013年12月
目录 程序作业题目： (3) 1、程序编制总说明 (3) 2、Matlab程序编制流程图 (3) 3、程序主要标示符及变量说明 (4) 4、理论基础和求解过程 (5) 4.1、构造插值函数 (5) 4.2位移插值函数及应变应力求解 (5) 5.程序的验证 (6) 附录：程序代码 (15)
程序作业题目： 完成一个包含以下所列部分的完整的有限元程序( Project) 须提供如下内容的文字材料（1500字以上）： ①程序编制说明； ②方法的基本理论和基本公式； ③程序功能说明； ④程序所用主要标识符说明及主要流程框图； ⑤ 1～3 个考题：考题来源、输出结果、与他人成果的对比结果（误差百分比）； ⑥对程序的评价和结论（包括正确性、适用范围、优缺点及其他心得等）。 须提供源程序、可执行程序和算例的电子文档或文字材料。选题可根据各自的论文选题等决定。 1、程序编制总说明 a.该程序采用平面三角形等参单元，能解决弹性力学的平面应力、平面应变问题。 b.能计算单元受集中力的作用。 c.能计算结点的位移和单元应力。 d.考题计算结果与理论计算结果比较，并给出误差分析。 e.程序采用MATLAB R2008a编制而成。 2、Matlab程序编制流程图
图1 整个程序流程图 3、程序主要标示符及变量说明 1、变量说明： Node ------- 节点定义 gElement ---- 单元定义 gMaterial --- 材料定义，包括弹性模量，泊松比和厚度 gBC1 -------- 约束条件 gNF --------- 集中力 gk------------总刚 gDelta-------结点位移 输入结构控制参数 输入其它数据 形成整体刚度阵 引入支承条件 解方程，输出位移 求应力，输出应力 形成节点荷载向量 开始 结束 1 单元面积 求弹性矩阵 单元刚度矩阵 位移-应变矩阵 6 7 8 9 10 2 3 4 5
MAPGIS中坐标转换中七参数法
MAPGIS 中坐标转换中七参数法 京54坐标系和西安80坐标系之间的转换其实是两种不同的椭球参数之间的转换，一般而言比较严密的是用七参数布尔莎模型，即X 平移，丫平移，Z平移，X旋转（WX，丫旋转（WY，Z旋转（WY，尺度变化（DM。若得七参数就需要在一个地区提供3个以上的公共点坐标对（即北京54坐标下x、y、z和西安80坐标系下x、y、z），可以向地方测绘局获取。 下面具体的步骤： 启动“投影变换模块”，单击“文件”菜单下“打开文件”命 令，将演示数据“演示数据_北京54.WT、“演示数据_北京 54.WL、“演示数据—北京54.WP打开。1、单击“投影转换” 菜单下“S坐标系转换”命令，系统弹出“转换坐标值” “话框⑴、在“输入”一栏中，坐标系设置为“北京54坐标系”，单位设置为“线类单位—米”；⑵、在“输出”一栏中，坐标系设置为“西 安80坐标系”，单位设置为“线类单位—米”；⑶、在“转换方法”一栏中，单击“公共点操作求系数”项；⑷、在“输入”一栏中， 输入北京54坐标系下一个公共点的（x、y、z）,如图2所示；⑸、在“输出”一栏中，输入西安80坐标系下对应的公共点的（x、y、z）, 如图2所示；⑹、在窗口右下角，单击“输入公共点”按钮，右边的数字变为1,表示输入了一个公共点对，如图2所示；⑺、依照相同的方法，再输入另外的2个公共点对；⑻、在“转换方法”一
栏中，单击“七参数布尔莎模型”项，将右边的转换系数项激活；⑼、 单击“求转换系数”菜单下“求转换系数”命令，系统根据输入的3个公共点对坐标自动计算出7个参数，如图3所示，将其记录下来；2、单击“投影转换”菜单下“编辑坐标转换参数”命令，系统弹出“不同地理坐标系转换参数设置”对话框，如图4所示；在“坐标系选项”一栏中，设置各项参数如下：源坐标系：北京54坐标系；目的坐标系：西安80坐标系；转换方法：七参数布尔莎模型；长度单位：米；角度单位：弧度；然后单击“添加项”按钮，则在窗口左边的“不同椭球间转换”列表中将该转换关系列出；在窗口下方的“参数设置”一栏中，将上一步得到的七个参数依次输入到相应的文本框中，如图4所示；单击“修改项”按钮，输入转换关系，并单击“确定”按钮；接下来就是文件投影的操作过程了。 3、单击“投影转换”菜单下“ MAPGI毀影转换/选转换线文件”命令，系统弹出“选择文件”对话框 选中待转换的文件“演示数据_北京54.WL',单击“确定”按 钮； 4、设置文件的Tic点，在“投影变换”模块下提供了两种方法：手工设置和文件间拷贝，这里不作详细的说明； 5、单击“投影转换”菜单下“编辑当前投影参数”命令，系统弹出 “输入投影参数”对话框，如图6所示，根据数据的实际情况来设置 其地图参数坐标系类型：大地坐标系 椭球参数：北京54投影类型：高斯-克吕格投影比例尺分母：1坐标单
arcgis七参数精确转换
在ArcGIS Desktop中进行三参数或七参数精确投影转换 Desktop, 投影, ArcGIS, 参数 ArcGIS中定义的投影转换方法，在对数据的空间信息要求较高的工程中往往不能适用，有比较明显的偏差。在项目的前期数据准备工作中，需要进行更加精确的三参数或七参数投影转换。下面介绍两种办法来在ArcGIS Desktop中进行这种转换。 方法1：在ArcMap中进行动态转换(On the fly) 假设原投影坐标系统为Xian80坐标系统，本例选择为系统预设的Projected Coordinate Systems\Gauss Kruger\Xian 1980\Xian 1980 GK Zone 20投影，中央经线为117度，要转换成Beijing 1954\Beijing 1954 GK Zone 20N。 在ArcMap中加载了图层之后，打开View-Data Frame Properties对话框，显示当前的投影坐标系统为Xian 1980 GK Zone 20，在下面的选择坐标系统框中选择Beijing 1954 GK Zone 20N，在右边有一个按钮为Transformations...
点击打开一个投影转换对话框，可以在对话框中看到Convert from和Into表明了我们想 从什么坐标系统转换到什么坐标系统。
在下方的using下拉框右边，点击New...，新建一个投影转换公式，在Method下拉框中可以选择一系列转换方法，其中有一些是三参数的，有一些是七参数的，然后在参数表中输 入各个转换参数。 输入完毕以后，点击OK，回到之前的投影转换对话框，再点击OK，就完成了对当前地图的动态投影转换。这时还没有对图层文件本身的投影进行转换，要转换图层文件本身的投影，
交点法与线元法
本人一直以来想找一个交点法与线元法相结合的坐标正反算程序，在网上找了很久很久，没能找到一个较为满意的，有幸在测量空间看到大歪哥的《Casio5800交点法程序》与《线元法（积木法）匝道坐标正反算放样程序》，根据歪哥意见“需要的自行修改结合XY框架自己修改为数据库反算程序等”，本人不才，采用最笨的办法将两个程序综合了一下，使之能既能进行交点法正反算，又能进行线元法正反算。在此特别感谢大歪哥！将程序发上来，愿与大家一同交流学习欢迎大家吐口水，只要能进步就行！ 程序由一个主程序ZBZFS和8个子程序（JS、XY-A、XY-B、JDYS、1、2、3、4）构成，运行时只需运行主程序即可！ 本程序适用于单交点对称型、不对称型、无缓和曲线单圆曲线型一个交点范围内（含交点前后有直线段时）的曲线要素核对和坐标正反算，手工输入要素，对设计图纸的“直线、曲线转角表”中交点数据进行复核验证，并能对单一线元进行坐标正反算。 1主程序名：ZBZFS（功能：进入计算主程序） 65→Dimz↙ Deg:Fix 3↙ ＂1.JD ZFS 2. ZHADAO ZFS＂? I: I→Z[61]: ＂1.ZHONG SHU JS 2. JS＂? I↙ If I=1: Then Goto1: Else Goto2:IfEnd↙ LbI 1 :If Z[61]=1: Then Prog＂JDYS＂:Else Cls:＂K0＂?A:＂KN＂?L :＂X0＂?U :＂Y0＂?V :＂F0＂?W :＂R0＂?P :＂RN＂?Q:＂ZX:-1,+1,0＂?G:IfEnd↙ LbI 2 :Prog＂JS＂ 2子程序名：JS（功能：选择正算或反算模式） Cls:＂XC＂?H:＂YC＂?Z↙ Cls:＂1.ZS 2.FS＂? I: I=2=>Goto 3↙ LbI 1 : Cls: If Z[61]=1: Then＂JD ZS KX+XXX＂?K :Prog＂4＂: Else ＂ZHADAO ZS KX+XXX＂?K :IfEnd↙ LbI 2: Cls:90→B: Cls:＂RJ Or 0 To K＂?B:B=0 =>Goto 1:＂Z＂?T↙ Prog ＂XY-A＂↙ X+Tcos(M+B)→X↙ Y+Tsin(M+B)→Y↙ 360Frac((M+360)÷360→M↙ Pol(X-H,Y-Z : 360Frac((J+360)÷360→J↙ 2→O: Prog ＂XY-B＂:Goto 2↙ LbI 3 : Cls: If Z[61]=1: Then＂JD FS KN+＂?K:＂X＂?C:＂Y＂?D:Prog＂4＂:Else Cls: ＂ZHADAO FS＂:＂X＂?C:＂Y＂?D:IfEnd↙ LbI 4 :Prog ＂XY-A＂↙ (D-Y)sin(M)+(C-X)cos(M)→H↙ If Abs(H)＞X10-3 :Then K+H→K:Goto 4:IfEnd↙ (D-Y)÷cos(M)→T↙ 3→O: Prog ＂XY-B＂:Goto 3↙ 3子程序名：XY-A（功能：坐标计算程序） 5→N: G(Q-1-P-1)÷Abs（L-A）→F: Abs（K-A）÷N→R: 90R÷π→S:
手持GPS三参数计算方法
手持GPS三参数计算方法 南方测绘石家庄工程项目部靳超 新机拿到手之后，设计方都给提供一个投影参数，这对于要求不高的一般用户来说基本可以满足工作需要，而对于一些专业用户来说，就要自己来测算参数。 一般型号的导航型手持GPS自定义坐标系统（User）投影参数设置界面都提供了五个变量（△X、△Y、△Z、△A、△F）需要设置，而实际工作中，后两个参数（△A、△F）针对某一坐标系统来说为固定参数（北京54坐标系△A=-108、△F=0.0000005西安80△A=-3、△F=0），无需改动，需要自己测算的参数主要为前三个（△X、△Y、△Z），一般称为三参数。 测算三参数的基本方法是，首先在已知控制点上测量一个稳定的WGS-84大地坐标（BLH）值，然后，运用专用测量程序既可算出一个三参数来。三参数计算出来后，将其输入GPS中再到已知控制点上观测比对，最好再到另一已知控制点上观测检校，如比对检校差值在规定允许误差范围之内，既可运用于实际工程测量工作。一般来说，只要到一新工区或工程点间距较远（数十至上百公里以外）都要到已知控制点上重新进行观测比对检校，没有问题才能进行实际工作。 三参数的求取步骤如下： 一、获取已知点的经纬度 利用手持GPS到一个已知控制点上测量一个稳定（即精度比较高）的WGS84大地坐标（即B,L,H）,也就是在手持GPS中将坐标系设置为:WGS84坐标系，显示格式为：经纬度格式。每种手持机设置的位置有所不同，请参阅说明书进行操作。 二、计算转换参数 一般手持机参数为：△X、△Y、△Z、△A、△F。△A、△F在北京54和西安80为固定值，我们主要计算：△X、△Y、△Z,即三参数。 我们使用COORD4.1(在此，感谢软件的作者：Jerry , 注意网上有其它版本的软件，某些功能可能有错误,如4.2版本)软件来自己求三参数。打开软件如图
5800线元法程序
1. 主程序 Lbl 4："1.SZ => XY"："2.XY => SZ"：？N：？S：Prog“SUB0”↙ 1÷P→C: (P-R)÷(2HPR) →D:180÷∏→E:N=1 => Goto1：Goto2↙ Lbl 1:? Z：Abs(S-O) →W：Prog "SUB1"："XS="：X◢ "YS="：Y◢F-90→F：“FS=”：F▲DMS◢：Goto4↙ Lbl 2：？X：？Y：X→I：Y→J：Prog“SUB2”：O+W→S：“S=”：S◢“Z=”：Z◢Goto4↙ 2. 正算子程序(SUB1) 0.1739274226→A：0.3260725774→B：0.0694318442→K：0.3300094782→L：1-L→F：1-K→M：U+W(Acos(G+QEKW(C+KWD))+Bcos(G+QELW(C+LWD))+Bcos(G+QEFW (C+FWD))+Acos(G+QEMW(C+MWD))) →X： V+W(Asin(G+QEKW(C+KWD))+Bsin(G+QELW(C+LWD))+Bsin(G+QEFW(C+FWD))+Asin(G+QEMW(C +MWD))) →Y： G+QEW(C+WD)+90→F：X+Zcos（F）→X：Y+Zsin（F）→Y 2. 反算子程序(SUB2) G-90→T：Abs((Y-V)cosT-(X-U)sinT) →W：0→Z：Lbl 0：Prog "SUB1"： T+QEW(C+WD) →L：(J-Y)cos（L）-(I-X)sin（L）→Z：Abs（Z）<1E-6=>Goto1：W=W+Z：Goto 0←┘Lbl 1：Z=0：Prog "SUB1"：(J-Y)÷sin（F）→Z SUB0 数据库子程序 Goto1↙同时保存多个曲线时的指针 Lbl 1：IF S<***（线元终点里程）:Then***→G（线元起点方位角）:***→O（线元起点里程）:***→U（线元起点X）:***→V（线元起点Y）:***→P（线元起点曲率半径）:***→R（线元终点曲率半径）: ***→H（线元起点至终点长度）：0或1、-1→Q：Return:IfEnd↙ Lbl 1：IF S<***（线元终点里程）:Then***→G（线元起点方位角）:***→O（线元起点里程）:***→U（线元起点X）:***→V（线元起点Y）:***→P（线元起点曲率半径）:***→R（线元终点曲率半径）: ***→H（线元起点至终点长度）：0或1、-1→Q：Return:IfEnd ……………………….. 为了便于解读，每增加一个线元增加一行语句，每增加一条曲线增加一个Lbl，每增加一个工程增加一个文件。 三、使用说明 1、规定 (1) 以道路中线的前进方向（即里程增大的方向）区分左右；当线元往左偏时，
GPS七参数的计算
通过三个或三个以上已知点求解七参数模型中的参数： 不同空间直角坐标系之间的变换，其参数有（ΔX0，ΔY0，ΔZ0，ωX，ωY，ωZ，m）七个，其中（ΔX0，ΔY0，ΔZ0）为坐标平移量，（ωX，ωY，ωZ）为坐标轴间的三个旋转角度（又称为欧拉角），m为尺度因子。七参数模型如图。 以WGS84坐标系转换为北京54坐标系为例： 为计算模型中的七个参数，至少需要三个已知点的北京54空间坐标 （X,Y,Z）BJ54和WGS-84空间坐标（X,Y,Z）WGS84，利用最小二乘法求出七参数。 然而，我们已知的三个公共控制点的坐标成果，一种是GPS观测中可直接获得的WGS84椭球下的大地坐标经纬度（B,L,H），另一种是工程测量中使用的是高斯投影后的平面直角坐标（x,y,h）。即已知的三个公共控制点的坐标成果就是这两种形式的坐标表来表示的。首先，我们要把这两种形式的坐标都转换为七参数模型中的空间直角坐标。步骤如下： 1.将WGS84椭球下的大地坐标经纬度（B,L,H），采用WGS84椭球参数，转换为WGS84的空间直角坐标（X,Y,Z） 2.将北京54投影平面直角坐标（x,y,h），采用克拉索夫斯基椭球参数，转换为大地坐标（（B,L,H）后，再转换为北京54的空间直角坐标（X,Y,Z）。 3.将转换得到的三个公共点的北京54空间坐标（X,Y,Z）BJ54和WGS-84空间坐标（X,Y,Z）WGS84代入七参数模型中，求解七个参数。 以上转换过程十分复杂，即涉及到大地坐标经纬度与空间直角坐标的换算，还涉及到空间直角坐标与平面直角坐标的投影。 通常，我也使用已有的计算程序来求解七参数的，在很多这些求解七参数的程序中，直接采用的是WGS84的大地坐标和北京54大地坐标来计算，就是你只需输入三个已知点的一套WGS84的大地坐标和一套北京54大地坐标，即可为你求解出七参数。
线元法
5800 计算程序 主程序 QXJS Fix 3:Deg:Lbl 4:“1.SZ=>XY”:“2.XY=>SZ”:? Q Lbl 4: “LICHENG= ” ?S:Prog“SUB0” ↙ Lbl 0:If Q=1:Then Goto1:IfEnd IfQ=2:ThenGoto2:IfEnd ↙ Lbl 1:”-B,0,B=”? Z: “J J右交角=”?G:Prog“SUB1”: Fix 4:Cls “X=”:N →N ◢“X=”: Locate3,1,N◢ “Y=”:E →E ◢“Y=”: Locate3,1,E◢ Prog“JI”:Goto4 “QXFWJ=”:F →F:F ▲ DMS ◢ Goto4 ↙ Lbl 2: “X=”? B: “Y=”? C:B→N: C→E:Prog“SUB2”: “LICHENG=”:S◢ “OUT JL=”:Z◢ Goto4 ↙ 说明： Q: 代表正反算，其中 1 为正算， 2 为反算； S: 代表里程； Z ：代表偏移距离； G ：代表偏移角度（以线路前进方向为 X 方向，顺时针转为正； N ： X 坐标； E ： Y 坐标； F ：切线方位角； JI Clstat Pol(N-G,-E-H):Cls If S<0:Then J+360→Y:Ease J→Y:Ifend “F W J=”:Y▲ DMS ◢黄色为计算机程序
SUB0 ( 数据库 ) Goto1 ↙ Lbl 1 IF S<157687.528:THEN 2884169.2517→U:471475.6573→V:157547.528→O:98 ° 32 ′ 43.08 ″ →A:140→L:10^45→P:10000→R: Return:IfEnd ↙ IF S<163781.879:THEN 2883008.7030→U:477458.2815→V:163641.879→O:101 ° 6 ′ 4.08 ″ →A:140→L:10^45→P:10000→R: Return:IfEnd ↙ IF S<164195.661:THEN 2882981.4268→U:477595.5984→V:163781.879→O:101 ° 30 ′ 7.93 ″ →A:413.7833→L:10000→P:10000→R: Return:IfEnd ↙ IF S<164335.661:THEN 2882890.5519→U:477999.2492→V:164195.6623→O:103 ° 52 ′ 22.82 ″ →A:140→L:10000→P:10^45→R: Return:IfEnd ↙ IF S<171831.142:THEN 2882856.3502→U:478135.0069→V:164335.6623→O:104 ° 16 ′ 26.67 ″ 说明： S ：里程；157547.528→O 为线元终点里程； 2884169.2517→U 为线元起点 X 坐标；471475.6573→V 为线元起点 Y 坐标；98 ° 32 ′ 43.08 ″ →A 线元起点切线方位角；0^45→P 线元起点半径（左转为负右转为正）；10000→R 线元终点半径（左转为负右转为正） SUB1 正算子程序 0.5 （1÷R-1÷P）÷L→D:S-O→X ↙ U+∫(cos(A+(X÷P+DX2)×180÷π,0,X)→N ↙ V+∫( sin(A+(X÷P+DX2)×180÷π,0,X)→E ↙ A+(X÷P+DX2)×180÷π→F ↙ N+Zcos(F+G) →N:E+Zsin(F+G) →E Return SUB2 反算子程序 Lbl 1:0→Z ：1→Q ：Prog“SUB0”: 0.5 （1÷R-1÷P ）÷L→D:S-O→X ↙ U+∫(cos(A+(X÷P+DX2)×180÷π,0,X)→N ↙ V+∫( sin(A+(X÷P+DX2)×180÷π,0,X)→E ↙ A+(X÷P+DX2)×180÷π→F ↙ N+Zcos(F+90) →N:E+Zsin(F+90) →E : Pol(N-B+10^(-46), E-C+10^(-46)):Isin(F-90-J) →W:S+W→S ↙