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) 以道路中线的前进方向（即里程增大的方向）区分左右；当线元往左偏时，
Q=-1；当线元往右偏时，Q=1；当线元为直线时，Q=0。
(2) 当所求点位于中线时，Z=0；当位于中线左铡时，Z取负值；当位于中线中线右
侧时，Z取正值。
(3) 当线元为直线时，其起点、止点的曲率半径为无穷大，以10的45次代替。
(4) 当线元为圆曲线时，无论其起点、止点与什么线元相接，其曲率半径均等于圆
弧的半径。
(5) 当线元为完整缓和曲线时，起点与直线相接时，曲率半径为无穷大，以10的45
次代替；与圆曲线相接时，曲率半径等于圆曲线的半径。止点与直线相接时，曲率半
径为无穷大，以10的45次代替；与圆曲线相接时，曲率半径等于圆曲线的半径。
(6) 当线元为非完整缓和曲线时，起点与直线相接时，曲率半径等于设计规定的
值；与圆曲线相接时，曲率半径等于圆曲线的半径。止点与直线相接时，曲率半径等
于设计规定的值；与圆曲线相接时，曲率半径等于圆曲线的半径。
2、输入与显示说明
输入部分：
1. SZ => XY
2. XY = > SZ
N ? 选择计算方式，输入1表示进行由里程、边距计算坐标；输入2表示由坐标反算里程和边距。
S ？正算时所求点的里程，反算时为所求点的近似里程，反算时输入的近似里程不能夸线元，若夸线元的话，需进行第二次反算，并以第一次的结果做为近似里程，否则反算的里程及边距不对，Z ？正算时所求点距中线的边距(左侧取负，值右侧取正值，在中线上取零)
X ？反算时所求点的X坐标
Y ？反算时所求点的Y坐标
显示部分：
XS=×××正算时，计算得出的所求点的X坐标
YS=×××正算时，计算得出的所求点的Y坐标
FS=×××正算时，所求点对应的中线点的切线方位角
S=×××反算时，计算得出的所求点的里程
Z=×××反算时，计算得出的所求点的边距
分两步计算：1、计算曲线要素 P=S2÷(24R);Q=S÷2-S3÷(240R2);T=(P+R)tan|0.5B|+Q;L=πRB÷
180+S ZH(C)=A-T;HY(D)=C+S;YH(H)=C+L-S;HZ(G)=C+L
式中：F----方位角；A----交点桩号；B-----路线偏角（左偏为负，右偏为正），R----半径；S---缓和曲线长
T----切线长；L-----曲线长
2、计算方位角：ZH=F;HY=F+S÷2R;HZ=F+B;YH=HZ-S÷2R;(当F＜0时S÷2R为负，大于0时为正)
线元法线路坐标正反算程序
经苦心钻研，奋战多日，终于编写出了代码短，速度快，精度高，功能全的线路坐标正反算程序，欢迎试用并提出宝贵意见。 功能简介及特点： 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Δ
线元法简介
线元法万能曲线正反算简介 我的线元法是把线形分为直线和曲线，直线就不用说了，起止点桩号，坐标和方位角就可以算了；曲线最基本的组合：是由一段缓和曲线+一段圆曲线组成，任意复杂的曲线都可以分解成缓和曲线+圆曲线或者其中之一就可以。 分析最复杂的曲线可以看到: 一般复杂线形由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/8378556.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处的平曲线， 的输入是否正确，有的图纸给的方位角数据较少，需要每隔几个线元才能检验方位角。
5800坐标程序
颜色文字Casio5800交点法与线元法（积木法）匝道坐标正反算放样程序（XUFENG 2011.2.14） 本人一直以来想找一个交点法与线元法相结合的坐标正反算程序，在网上找了很久很久，没能找到一个较为满意的，有幸在测量空间看到大歪哥的《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↙
公路坐标计算公式
一、缓和曲线上的点坐标计算 已知：①缓和曲线上任一点离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 计算过程： 五、超高缓和过渡段的横坡计算
5800简单全线坐标计算程序
5800全线任意坐标计算程序 1. 正算主程序(ZHCX) （不运行） 8→DimZ 1÷P→Z[4 ]：(P-R)÷(2HPR)→D: 180÷π→E “Z=”?Z:”YJJ=”?A:Abs(S-O)→W 0.26→Z[1 ]: 0.74→B: 0.02→K: 0.82→Z[3 ]: 1-Z[3 ]→F:1-K→Z[2 ] U+W(Z[1 ]cos(G+QEKW(Z[4 ]+KWD))+Bcos(G+Z[3 ]QEW(Z[4 ]+ Z[3 ]WD))+Bcos(G+QEFW (Z[4 ]+FWD))+ Z[1 ]cos(G+ Z[2 ]QEW(Z[4 ]+ Z[2 ]WD)))→X: V+W(Z[1 ] sin (G+QEKW(Z[4 ]+KWD))+B sin(G+ Z[3 ]QEW(Z[4 ]+ Z[3 ]WD))+B sin(G+QEFW (Z[4 ]+FWD))+ Z[1 ] sin(G+ Z[2 ]QEW(Z[4 ]+ Z[2 ]WD)))→Y： G+QEW(Z[4 ]+WD)→F：X+Zcos(F+A)→X：Y+Zsin(F+A)→Y:If F≧360:Then F-360→F:IfEnd ”X=”:X→X◢ ”Y=”:Y→Y◢ If F﹤0:Then F+360→F:IfEnd ”QX FWJ=”:F▼DMS◢ “C=1=>XX: C=2=>XZ”: ”C=”?C: ”QHJU=”?L: If C=1：Then Goto 1:Else Goto 2: IfEnd 可以计算斜交斜做或斜交正做的桥涵坐标 Lbi 1 X+L cos(F)→X：Y+Lsin(F)→Y: Goto 3 Lbi 2 X+L cos(F+A-90)→X：Y+Lsin(F+A-90)→Y: Goto 3 Lbi 3 “QH-X=”: X →X◢ “QH-Y=”: Y →Y◢ Prog “FY” 2 . 参数子程序（直接运行） M(主线) 一条线路一个名称 “S=”?S If S≦线元终点：Then 线元起点X值→U: 线元起点Y值→V:线元起点切线方位角→G:线元起点桩号→O:线元长度→H：线元起点半径→P：线元终点半径→R：（左偏-1，或右偏 1）→Q:Goto 1:IfEnd … … If S≦线元终点：Then 线元起点X值→U: 线元起点Y值→V:线元起点切线方位角→G:线元起点桩号→O:线元长度→H：线元起点半径→P：线元终点半径→R：（左偏-1，或右偏 1）→Q:Goto 1:IfEnd Lbi 1 Prog “ZBJS” 3. 放样程序（FY）（不运行） “X0=”?M:“Y0=”？N Pol((X-M, Y-N)
匝道线元法
任意里程中边桩坐标正反算(CASIO fx-4800P计算器)程序 时间:2007-12-18 13:13:15 来源:工测员网作者:未知 1. 加编数据库作为主程序 , 计算中不必逐项输入 " 线元要素 ", 提高运算速度，避免现场忙中出错 2. 将原来的主程序并入数据库 3. 计算直观 , 人性化 4. 正算直接输入里程和边距 , 反算输入近似里程便可 5. 增加了“ 计算点与测站点” 的距离和方位角计算语句，方便直接放样 6. 愿收获与大家共享 7. 核心计算程序摘自“yshf” 1. 正算子程序 (SUB1) A=0.1739274226 ： B=0.3260725774 ： K=0.0694318442 ： L=0.3300094782 ：F=1-L ： M=1-K ： X=U+W(Acos(G+57.2958QKW(1/P+KWD))+Bcos(G+57.2958QLW(1/P+LWD))+Bcos(G+ 57.2958QFW (1/P+FWD))+Acos(G+57.2958QMW(1/P+MWD))) ： Y=V+W(Asin(G+57.2958QKW(1/P+KWD))+Bsin(G+ 57.2958QLW(1/P+LWD))+Bsin(G+57.2958QFW(1/P+FWD))+Asin(G+57.2958QMW(1/ P+MWD))) ： F=G+57.2958QW(1/P+ WD)+90 ： X=X+ZcosF ： Y=Y+ZsinF 2. 反算子程序 (SUB2) T=G-90 ： W=Abs((Y-V)cosT-(X-U)sinT) ： Z=0 ： Lbl 0 ： Prog "SUB1" ：L=T+57.2958QW(1/P+ WD) ： Z=(J-Y)cosL-(I-X)sinL ： AbsZ<1E-6=>Goto1 ：≠>W=W+Z ： Goto 0Δ←┘ Lbl 1 ： Z=0 ： Prog "SUB1" ： Z=(J-Y)÷sinF 3、 . 增设数据库程序（ＳＪＫ主程序） Lb1 4 ： "1.SZ => XY" ： "2.XY => SZ" ：｛ NS ｝：S ∠ 下一线元起点里程＝＞O = 本线元起点里程： U= 本线元起点 X ： V= 本线元起点 Y ： G= 本线元起算方位角： H= 本线元长度： P= 起点曲率半径： R= 终点曲率半径： Q=0 或 1 、 -1 ： Goto0Δ←┘( 第一线元数据要素) S ∠ 下一线元起点里程＝＞O = 本线元起点里程： U= 本线元起点 X ： V= 本线元起点 Y ： G= 本线元起算方位角： H= 本线元长度： P= 起点曲率半径： R= 终点曲率半径： Q=0 或 1 、 -1 ： Goto0Δ←┘( 第二线元数据要素) S ∠ 下一线元起点里程＝＞O = 本线元起点里程： U= 本线元起点 X ： V= 本线元起点 Y ： G= 本线元起算方位角： H= 本线元长度： P= 起点曲率半径： R= 终点曲率半径： Q=0 或 1 、 -1 ： Goto0Δ←┘( 第三线元数据要素)
fx-5800p全线坐标正反算带高程计算程序(线元法)
曲线任意里程中边桩坐标正反算(CASIO fx-5800P计算器)程序 一、程序功能及原理 1.功能说明：本程序由一个主程序(TYQXJS)和五个子程——正算子程序(SUB-ZS)、反算子程序(SUB-FS)等构成，可以根据曲线段——直线、圆曲线、缓和曲线（完整或非完整型）的线元要素（起点坐标、起点里程、起点切线方位角、线元长度、起点曲率半径、止点曲率半径）及里程边距或坐标，对该曲线段范围内任意里程中边桩坐标进行正反算。本修改版程序既可实现正算全线贯通，亦可实现反算全线贯通。本程序在CASIO fx-5800P计算器运行。 2．计算原理：利用Gauss-Legendre 5点通用公式正算线路中边桩坐标、线外测点至曲线元起点和终点的垂距的符号是否相异（即Dca×Dcb<0=>该测点在其线元内）进行判断并利用该线元要素反算中桩里程、支距，最后计算出放样数据。 二、源程序 1.主程序(TYQXJS)(A) Deg:fix 3 119→DimZ “INPUT(0) Or DATA(Else)”?I Lbl 0:“1.SZ=>XY,2.XY=>SZ，3.TF=>CK,4.SD=>FY，5.TW=>FY”?N If N=1 Or N=5:Then Goto 1 Else If N=2 Or N=3 Or N=4:Then Goto 2 Else Goto 3 IfEnd:IfEnd Lbl 1:“K(m)=”?S If S<0:Then Goto 0:IfEnd “JL(m)=”?Z If Z≠0:Then “ANGLE→R(Deg)=”?M:IfEnd If I=0:Then Prog “DAT1”:Else Prog “DAT2”:IfEnd S-O→W:If W<0:Then Goto 0:Else If W>H:Then Goto 0:IfEnd:IfEnd Prog “SUB-ZS”:Prog “SUB-GC” If Z<0:Then“XL(m)=”:X◢“YL(m)=”:Y◢ If N=5:Then Prog “SUB-TW”:IfEnd Else If Z>0:Then “XR(m)=”:X◢“YR(m)=”:Y◢ If N=5:Then Prog “SUB-TW”:IfEnd Else “X(m)=”:X◢“Y(m)=”:Y◢“Hs(m)=”:L◢“FWJ=”: F?DMS◢ IfEnd:IfEnd
边界元与有限元
边界元与有限元 边界元法boundary element method 定义：将力学中的微分方程的定解问题化为边界积分方程的定解问题，再通过边界的离散化与待定函数的分片插值求解的数值方法。 所属学科：水利科技(一级学科) ；工程力学、工程结构、建筑材料(二级学科) ；工程力学（水利）(三级学科) 边界元法（boundary element method）是一种继有限元法之后发展起来的一种新数值方法，与有限元法在连续体域内划分单元的基本思想不同，边界元法是只在定义域的边界上划分单元,用满足控制方程的函数去逼近边界条件。所以边界元法与有限元相比，具有单元个数少,数据准备简单等优点.但用边界元法解非线性问题时，遇到同非线性项相对应的区域积分，这种积分在奇异点附近有强烈的奇异性，使求解遇到困难。 简介 边界元法是在有限元法之后发展起来的一种较精确有效的工程数值分析方法。又称边界积分方程-边界元法。它以定义在边界上的边界积分方程为控制方程，通过对边界分元插值离散，化为代数方程组求解。它与基于偏微分方程的区域解法相比，由于降低了问题的维数，而显著降低了自由度数，边界的离散也比区域的离散方便得多，可用较简单的单元准确地模拟边界形状，最终得到阶数较低的线性代数方程组。又由于它利用微分算子的解析的基本解作为边界积分方程的核函数，而具有解析与数值相结合的特点，通常具有较高的精度。特别是对于边界变量变化梯度较大的问题，如应力集中问题，或边界变量出现奇异性的裂纹问题，边界元法被公认为比有限元法更加精确高效。由于边界元法所利用的微分算子基本解能自动满足无限远处的条件，因而边界元法特别便于处理无限域以及半无限域问题。边界元法的主要缺点是它的应用范围以存在相应微分算子的基本解为前提，对于非均匀介质等问题难以应用，故其适用范围远不如有限元法广泛，
CASIO5800线元法程序
卡西欧5800P 线元法正反算程序 1-JIN YU(主程序) Lbl 4: “1→ZS,2→FS?N(选择计算模式，1为正算，2为反算，3为超欠挖) N=1=>Goto 1:N=2=>Goto 2 Lbl 1: “K=”?S：“P=”?Z：Prog “SJ-PM”： Abs(S-O) → W：Prog “SUB1-ZS”：“X=”：Locate4,4,X:“Y=”：Locate4,4,Y:F-90→F: “W=”: Locate4,4,F°◢ Goto 4 (正算-输入待求点里程K=、输入待求点偏距P=、显示待求点里程X=、显示待求点里程Y=、显示待求点里程设计高程H=) Lbl 2:“X=”?X：“Y=”？Y:Prog“SJ-PM”： X→ I： Y→J：Prog "SUB2-FS"：O+W→S: “K="：Locate4,4, S：“P="：Locate4,4, Z：S→ K:Prog“SJ-GC”：“H=”：Locate4,4, H:F-90→F: “W=”: Locate4,4,F°◢ Goto 4 (反算-输入实测点X=、输入实测点Y=、显示实测点里程K=、显示实测点偏距P=、显示实测点的设计高程H=) SJ-PM（子程序名-平面线形数据库） If S ≥59227.681（线元起点里程）:Then 2599818.013→U（线元起点X坐标）：496887.918→V（线元起点Y坐标）：59227.681→O（线元起点里程）：208028’39”→G（线元起点方位角）：100 →H（线元长度）：1×1045→P（线元起点曲率半径）：1×1045→R（线元终点曲率半径）：0→Q（线元左右偏标志：左-1右1直0）：IfEnd If S ≥59327.681（线元起点里程）:Then 2599730.112→U（线元起点X坐标）：496840.237→V（线元起点Y坐标）：59327.681→O（线元起点里程）：208028’39”→G（线元起点方位角）：90 →H（线元长度）：1×1045→P（线元起点曲率半径）：750→R（线元终点曲率半径）：-1 →Q（线元左右偏标志：左-1右1直0）：IfEnd (注：有多个平面线形，依照上面的依次变更, 每多一个，就增加一个。每次只需要修改以上的数据或增加一个判断。) SUB1-ZS(正算子程序)
交点法与线元法
本人一直以来想找一个交点法与线元法相结合的坐标正反算程序，在网上找了很久很久，没能找到一个较为满意的，有幸在测量空间看到大歪哥的《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:
通用显式非线性有限元程序：LS-DYNA
通用显式非线性有限元程序：LS-DYNA LS-DYNA 是世界上最著名的通用显式非线性有限元分析程序，能够模拟真实世界的各种复杂问题，特别适合求解各种二维、三维非线性结构的碰撞、金属成型等非线性动力冲击问题，同时可以求解传热、流体及流固耦合问题。在工程应用领域被广泛认可为最佳的分析软件包。与实验的无数次对比证实了其计算的可靠性。 LS-DYNA 是功能齐全的几何非线性（大位移、大转动和大应变）、材料非线性（140多种材料动态模型）和接触非线性（50多种）软件。它以Lagrange 算法为主，兼有ALE 和Euler 算法；以显式求解为主，兼有隐式求解功能；以结构分析为主，兼有热分析、流体-结构耦合功能；以非线性动力分析为主，兼有静力分析功能（如动力分析前的预应力计算和薄板冲压成型后的回弹计算）;是通用的结构分析非线性有限元程序。 特色功能 ? 显式求解为主，兼有隐式算法，适合于求解高度非线性问题； ? 具有多种求解算法，以Lagrange 算法为主，兼有ALE、Euler 算法、SPH （Smoothed Particle Hydrodynamics）光顺质点流体动力算法和边界元法BEM（Boundary Element Method）； ? 具有160多种材料模型，是材料模型非常丰富的有限元软件； ? 具有50多种接触类型，是接触类型非常齐全的有限元软件； ? 极好的并行计算能力，包括分布式并行算法（MPP）和共享内存式并行（SMP）； ? 良好的自适应网格剖分技术，包括自适应网格细分和粗化； ? 行业化的专用功能：如针对汽车行业的安全带单元、滑环、预紧器、牵引器、传感器、加速计、气囊等。 客户价值 ? 拥有显式和隐式算法，各向异性材料模型，使得板成型、回弹、预应力计算等，可以连续求解； ? 多种控制选项和用户子程序使得用户在定义和分析问题时有很大的灵活性； ? MPP 版本大幅度减少计算时间，计算效率随计算机数目增多而显著提高； ? 与大多数的CAD/CAE 软件集成并有接口。 广州有道科技培训中心 h t t p ://w w w .020f e a .c o m
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) 以道路中线的前进方向（即里程增大的方向）区分左右；当线元往左偏时，
线元法
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 ↙
边界元法发展综述
边界元法发展综述 刘娅君 学号：11080922005 从工程实际中提出的力学问题，一般可归结为数学的定解问题。但其中只有极少数简单情况可以求得解析解，而大多情况都必需借助于有效的数值方法来求解。有限元法是目前工程中应用最广泛的数值方法，已有很多通用程序和专用程序在各个工程领域投人了实际应用。然而，有限元法本身还存在一些缺点。例如，在应力分析中对于应力集中区域必须划分很多的单元，从而增加了求解方程的阶数，计算费用也就随之增加；用位移型有限元法求解出的应力的精度低于位移的精度，对于一个比较复杂的问题必须划分很多单元，相应的数据输人量就很大，同时，在输出的大量信息中，又有许多并不是人们所需要的。 边界积分方程—边界元法在有限元法之后发展起来成为工程中广泛应用的一种有效的数值分析方法。它的最大特点就是降低了问题的维数，只以边界未知量作为基本未知量，域内未知量可以只在需要时根据边界未知量求出。在弹性问题中，由于边界元法的解精确满足域内的偏微分方程，因此它相对有限元法的解具有较高的精度。同时在一些领域里，例如线弹性体的应力集中问题，应力有奇异性的弹性裂纹问题，考虑脆性材料中裂纹扩展的结构软化分析，局部进人塑性的弹塑性局部应力问题以及弹性接触问题…等，边界元法已被公认为比有限元法更为有效。正是因为这些特点，使边界元法受到了力学界、应用数学界及许多工程领域的研究人员的广泛重视。 边界元与有限元相比有很多优点：首先，它能使问题的维数降低一维，如原为三维空间的可降为二维空间，原为二维空间的问题可降为一维。其次，它只需将边界离散而不象有限元需将区域离散化，所划分的单元数目远小于有限元，这样它减少了方程组的方程个数和求解问题所需的数据，不但减少了准备工作，而且节约了计算时间。第三，由于它是直接建立在问题控制微分方程和边界条件上的，不需要事先寻找任何泛函，不像以变分问题为基础的有限元法，如果泛函不存在就难于使用。所以边界元法可以求解经典区域法无法求解的无限域类问题。最后，由于边界元法引入基本解，具有解析与离散相结合的特点，因而具有较高的精度。
交点法和线元法曲线要素输入简介
测量坐标计算程序V5 输入简介 本程序运用Office Excel 软件VBE标准模块编写，其功能基本全面集成了以往所更新的Excel程序，程序适用于公路、铁路等线路坐标计算，程序主要包括（交点法、线元法、直线坐标正反算，竖曲线计算，平面控制网“导线、高程”平差，隧道超欠挖，超高加宽，测量工具箱等，还可以全自动生成卡西欧5800、9750程序数据库，其中包括：隧道超欠挖、交点法、线元法、竖曲线一系列数据库），已知数据输入明确，操作简单易懂，是工程测量人员的好帮手！ 交点法曲线要素输入简介 一、适用平曲线类型 交点法计算坐标适用的平曲线为对称或不对称缓和曲线、圆曲线。 注意：对于非普通的三单元曲线，本程序交点法不适用。非普通的三单元曲线 体现在本程序中的《直线、曲线及转角表》内，点击“生成要素”之后，计算 值与设计图纸《直线、曲线及转角表》上的切线长和曲线主点位置等不一致， 此时只能采用线元法进行坐标计算。 例如：下表的JD18及JD19处的平曲线，经本程序交点法计算之后发现，为非普通的三单元曲线，交点法不适用该类曲线的坐标计算，故只能采用线元法进行坐标计算。
二、交点法曲线要素输入说明 本程序交点法输入的要素有7个（程序不限制输入行数）： 1、QD起点坐标： 起点坐标必须在直线段上，或填写前一交点的坐标。 2、JD交点曲线要素： （1）交点桩号K，注意：当起始平曲线上的ZH点（缓和曲线）或ZY点（圆曲线）的桩号为负数时，交点桩号K统一加上100000（即增加100Km），以避免坐标正算时出现桩号计算范围错误（但是，线元法计算坐标时可以输入负坐标，坐标正算与反算都不会出现错误）。 （2）交点桩号（X,Y） （3）曲线半径R （4）第一缓和曲线长度LS1,若为0,输入0,不能为空。 （5）第二缓和曲线长度LS2,若为0,输入0,不能为空。 3、ZD终点坐标： 终点坐标也必须在直线段上，或填写后一交点的坐标。 三、操作流程： 1、根据设计图纸《直线、曲线及转角表》输入第一个交点坐标，作为QD起点坐标。 2、依次输入各交点的曲线要素。 3、输入最后一个交点坐标，作为ZD终点坐标。 4、点击“点击进入直曲表”，然后点击“生成要素”，根据计算的转角值、曲线要素、曲线主点位置、直线长度及方向与设计图纸《直线、曲线及转角表》上的设计值进行核对，看各要素输入是否正确。 注意：交点法与线元法在计算坐标时，线元法可能存在mm毫米 级以内的计算误差，坐标计算值与设计值也可能存在mm毫米级以内