精华内容
下载资源
问答
  • 2021-10-29 13:46:04

     

     

     

     

     

     

     

    更多相关内容
  • 非球面镜具有优良的光学性能,在现代光学系统设计中广泛采用。...顶点曲率半径是非球面镜加工、检测和光学系统装校中的一个重要参数。介绍了一种基于光线追迹原理的非球面镜顶点曲率半径的测量方法和一些实验结果。
  • 由于影响曲率半径测量精度的因素众多,因此高精度曲率半径测量一直是光学检测中的难题之一。为了实现球面光学元件曲率半径的高精度测量,提出通过环境补偿和提高对准精度的方法来减小测量误差。首先从理论的角度分析...
  • 在误差理论中,有多种处理数据的方法。文章详细介绍了用逐差法、最小二乘法、加权平均法处理牛顿环测透镜曲率半径的数据的方法和过程。
  • #资源达人分享计划#
  • 给出了一般光束等效曲率半径在湍流大气和自由空间中的传输方程,推导出等效曲率半径达到最小值时光束的传输距离公式,解析地讨论了等效曲率半径的特性,并用数值计算例加以说明。研究表明,在自由空间中,一般光束...
  • 根据高精密光学系统的需求,结合目前的研究成果,基于干涉测量法研制了长曲率半径测量系统。分析了影响长曲率半径测量系统测量精度的多种误差因素,根据误差理论和系统组成建立了长曲率半径测量系统误差分配的数学...
  • 基于立式Fizeau 型干涉仪曲率半径测量系统,提出了一种简捷的离焦补偿算法,该算法对猫眼和共焦零位置的离焦量进行补偿,达到了提高曲率半径测精度的目的。从理论的角度分析了待测面波前像差中的离焦量与待测面位置...
  • 曲面主曲率半径在许多工程问题中是一个非常重要的几何量,微分几何对此有专门的计算方法。以研究平面曲线曲率半径的方法对旋转曲面主曲率半径进行图形处理,避开建立曲面方程,求解曲面第一类、第二类基本量的常规...
  • 经典的GW(Greenwood和Williamson)模型假设微凸体均为弹性,并没有考虑到微凸体磨损后曲率半径变化的情况。在经典赫兹接触理论与GW模型基础上推导出微凸体曲率为定值时的摩擦面剪切刚度公式。根据Mohr-Coulomb准则探讨...
  • 用牛顿环测曲率半径实验报告.pdf
  • 本文讨论了激光球面干涉仪以及数字干涉仪在球面曲率半径测量中的应用,同传统的球面样板检测进行了比较,分析并计算了两者的测量精度,给出了提高球径仪测量精度的原则性方案.
  • 分析了当带式输送机在垂直凹向弯曲运行状态时,为了保证输送带不失稳即边缘不产生负拉力时所需要保证的临界弯曲曲率半径。根据理论分析,建立了计算该临界曲率半径的数学模型。
  • 在简要总结各种检测光学球面曲率半径方法优缺点的基础上,提出了利用激光跟踪仪和激光干涉仪测量光学球面曲率半径的新方法。首先,通过激光跟踪仪精确定位测量干涉仪出射球面波前的焦点和待测球面镜的曲率中心点坐标...
  • 圆螺旋线求曲率半径

    2022-04-06 20:24:28
    看论文的时候看到给出了螺旋线的曲率半径求解公式: 便想着这个公式是怎么求出来的。 上述螺旋线我们知道h(转一圈,又称为导程),螺旋的半径是r。 我们可以用参数方程来描述螺旋线: {x=rcosθy=rsinθz=h2πθ\...

    引言

    看论文的时候看到给出了螺旋线的曲率半径求解公式:
    在这里插入图片描述
    便想着这个公式是怎么求出来的。

    在这里插入图片描述
    上述螺旋线我们知道h(转一圈,又称为导程),螺旋的半径是r。

    我们可以用参数方程来描述螺旋线:

    { x = r c o s θ y = r s i n θ z = h 2 π θ \left\{\begin{aligned} x&=rcos\theta \\ y&=rsin\theta \\ z&=\frac{h}{2\pi}\theta \end{aligned}\right. xyz=rcosθ=rsinθ=2πhθ

    为什么z的表达式中系数是 h 2 π \frac{h}{2\pi} 2πh呢?想一下, θ \theta θ角从0转到 2 π 2\pi 2π,这样整个旋转了一圈,把螺旋线拉开截面方向就从0到 2 π r 2\pi r 2πr,纵向方向从0到h,两个都是线性增加,很明显
    { z ( 0 ) = 0 z ( 2 π ) = h z ( θ ) = k θ \left\{\begin{aligned} z(0)&=0 \\ z(2\pi)&=h \\ z(\theta)&=k\theta \end{aligned}\right. z(0)z(2π)z(θ)=0=h=kθ

    则有

    k = z ( 2 π ) − z ( 0 ) 2 π − 0 = h 2 π k=\frac{z(2\pi)-z(0)}{2\pi-0}=\frac{h}{2\pi} k=2π0z(2π)z(0)=2πh

    这就解释了z的表达式中系数k是 h 2 π \frac{h}{2\pi} 2πh
    写出了参数方程,参考维基百科[1]参数方程的曲率公式(曲率半径是曲率的倒数):

    对于一个以参数化形式给出的空间曲线 r ( t ) = ( x ( t ) , y ( t ) , z ( t ) )   r ( t ) = ( x ( t ) , y ( t ) , z ( t ) )   {\displaystyle {\boldsymbol {r}}(t)=(x(t),y(t),z(t))\,}{\displaystyle {\boldsymbol {r}}(t)=(x(t),y(t),z(t))\,} r(t)=(x(t),y(t),z(t))r(t)=(x(t),y(t),z(t))其曲率为在这里插入图片描述

    列出需要的表达式:
    { x ′ ( θ ) = − r s i n θ x ′ ′ ( θ ) = − r c o s θ y ′ ( θ ) = r c o s θ y ′ ′ ( θ ) = − r s i n θ z ′ ( θ ) = h 2 π z ′ ′ ( θ ) = 0 \left\{\begin{aligned} x'(\theta)=&-rsin\theta \\ x''(\theta)=&-rcos\theta\\ y'(\theta)=&rcos\theta \\ y''(\theta)=&-rsin\theta\\ z'(\theta)=&\frac{h}{2\pi}\\ z''(\theta)=&0 \end{aligned}\right. x(θ)=x(θ)=y(θ)=y(θ)=z(θ)=z(θ)=rsinθrcosθrcosθrsinθ2πh0

    代入计算
    κ = ( z ′ ′ y ′ − y ′ ′ z ′ ) 2 + ( x ′ ′ z ′ − z ′ ′ x ′ ) 2 + ( y ′ ′ x ′ − x ′ ′ y ′ ) 2 ( x ′ 2 + y ′ 2 + z ′ 2 ) 3 2 = [ 0 ⋅ r c o s θ − ( − r s i n θ ) ⋅ h 2 π ] 2 + [ − r c o s θ ⋅ h 2 π − 0 ⋅ ( − r s i n θ ) ] 2 + [ ( − r s i n θ ) ⋅ ( − r s i n θ ) − ( − r c o s θ ) ⋅ r c o s θ ] 2 [ ( − r s i n θ ) 2 + ( r c o s θ ) 2 + ( h 2 π ) 2 ] 3 2 = ( h 2 π ) 2 r 2 + r 4 ( ( h 2 π ) 2 + r 2 ) 3 2 = r ( h 2 π ) 2 + r 2 ( ( h 2 π ) 2 + r 2 ) 3 2 = r ( h 2 π ) 2 + r 2 \begin{aligned} \kappa&=\frac{\sqrt{(z''y'-y''z')^2+(x''z'-z''x')^2+(y''x'-x''y')^2}}{(x'^2+y'^2+z'^2)^{\frac{3}{2}}}\\ &=\frac{\sqrt{[0\cdot rcos\theta-(-rsin\theta)\cdot \frac{h}{2\pi}]^2+[-rcos\theta\cdot \frac{h}{2\pi}-0\cdot(-rsin\theta)]^2+[(-rsin\theta)\cdot(-rsin\theta)-(-rcos\theta)\cdot rcos\theta]^2}}{[(-rsin\theta)^2+(rcos\theta)^2+(\frac{h}{2\pi})^2]^{\frac{3}{2}}} \\&=\frac{\sqrt{(\frac{h}{2\pi})^2r^2+r^4}}{((\frac{h}{2\pi})^2+r^2)^{\frac{3}{2}}}\\&=\frac{r\sqrt{(\frac{h}{2\pi})^2+r^2}}{((\frac{h}{2\pi})^2+r^2)^{\frac{3}{2}}}\\&=\frac{r}{(\frac{h}{2\pi})^2+r^2}\end{aligned} κ=(x2+y2+z2)23(zyyz)2+(xzzx)2+(yxxy)2 =[(rsinθ)2+(rcosθ)2+(2πh)2]23[0rcosθ(rsinθ)2πh]2+[rcosθ2πh0(rsinθ)]2+[(rsinθ)(rsinθ)(rcosθ)rcosθ]2 =((2πh)2+r2)23(2πh)2r2+r4 =((2πh)2+r2)23r(2πh)2+r2 =(2πh)2+r2r

    曲率求倒数,所以圆螺旋线曲率半径为

    ρ = ( h 2 π ) 2 + r 2 r \rho=\frac{(\frac{h}{2\pi})^2+r^2}{r} ρ=r(2πh)2+r2

    和论文所述一致

    参考

    维基百科-曲率

    展开全文
  • 摆动凸轮廓线曲率半径计算的CAD动画方法.pdf
  • 地球形状与曲率半径

    2021-09-05 09:54:41
    \qquad实际的地球表面是一个凹凸不平、形状十分复杂的...近似成圆球体, 中心在地球质心,半径6371km,适用于精度要求不高的场景。 \qquad 2.近似成旋转椭球体,地球自转轴与与椭圆短半轴重合,椭圆度1300\frac{1}{300

    \qquad 实际的地球表面是一个凹凸不平、形状十分复杂的物理面,难以准确量化描述。为了方便研究,通常以平静的海平面为基准,并把它向大陆延伸形成一个封闭曲面,称之为大地水准面。它所包含的的几何形状称为大地水准体。地球的实际形状如下图所示:在这里插入图片描述
    \qquad 实际应用中总是希望可以使用相对简单的数学形式来近似拟合地球几何形状,常用的近似有三种:
    \qquad 1.近似成圆球体, 中心在地球质心,半径6371km,适用于精度要求不高的场景。
    \qquad 2.近似成旋转椭球体,地球自转轴与与椭圆短半轴重合,椭圆度 1 300 \frac{1}{300} 3001(长短半轴相差21km),椭圆绕其短半轴旋转一周构成的旋转椭球体的表面,该描述中地区吃到近似为原型。
    \qquad 3.近似成三轴椭球体,其中赤道是椭圆的,赤道椭圆度约为 1 100000 \frac{1}{100000} 1000001(长短半轴相差约60m)。
    \qquad 三轴椭球体虽然描述更精确,但相关计算更加复杂,其赤道的椭圆度并不大,可近似为圆形,因此旋转椭球体的应用最为广泛。
    \qquad 包含南北极点的平面称为子午面,子午面与旋转椭球面的交线成为子午圈(经圈),经过英国格林尼治的经线称为本初子午线(零度经线)。任一经线所在子午面与本初子午面之间的夹角定义为经度 ( ϕ ) (\phi) (ϕ),夹角方向与地球自转方向相同,取值范围 − 180 ° ∼ 180 ° -180°\thicksim180° 180°180°。包含旋转椭球中心且垂直自转轴的平面称为赤道面,赤道面与旋转椭球面的交线称为赤道。平行于赤道的平面与旋转椭球面的交线成为纬圈在这里插入图片描述
    \qquad 对于地球旋转椭球体而言,确定其三维形状参数的关键在于确定二维子午圈椭圆。

    1.地球几何形状

    \qquad 建立地心右手坐标系,即地心地固坐标系 ( E C E F ) (ECEF) (ECEF),坐标原点在地心,与地球固联,跟随地球自转一起相对于惯性空间转动。取子午圈椭圆上任意一点 P P P与地心 O O O的连线 P O PO PO O X e OX_e OXe的夹角成为地心纬度,记为 ϕ \phi ϕ,取值范围 − 90 ° ∼ 90 ° -90°\sim90° 90°90°。南纬为负,北纬为正。过 P P P点的法线 P Q PQ PQ O X e OX_e OXe的夹角称为地理纬度。记为 L L L,取值范围 − 90 ° ∼ 90 ° -90°\sim90° 90°90° P O PO PO称为地心垂线, P Q PQ PQ称为地理垂线。
    \qquad 椭圆方程为:
    x 2 R e 2 + z 2 R p 2 = 1 (1.1) \frac{x^2}{R_e^2}+\frac{z^2}{R_p^2}=1 \tag{1.1} Re2x2+Rp2z2=1(1.1) \qquad 其中 R e , R p R_e,R_p Re,Rp分别为椭圆的长半轴和短半轴。
    \qquad 椭圆扁率(或称椭圆度,flattening)定义为:
    f = R e − R p R e (1.2) f=\frac{R_e-R_p}{R_e} \tag{1.2} f=ReReRp(1.2) \qquad 椭圆偏心率定义为:
    e = R e 2 − R p 2 R e (1.3) e=\frac{\sqrt{R_e^2-R_p^2}}{R_e}\tag{1.3} e=ReRe2Rp2 (1.3) \qquad 第二偏心率定义为:
    e ′ = R e 2 − R p 2 R p (1.4) e'=\frac{\sqrt{R_e^2-R_p^2}}{R_p} \tag{1.4} e=RpRe2Rp2 (1.4) \qquad 且有 e ′ 2 = e 2 1 − e 2 (1.5) e'^2=\frac{e^2}{1-e^2} \tag{1.5} e2=1e2e2(1.5) \qquad 由式(1.2)和式(1.3)分别可得:
    R p = ( 1 − f ) R e (1.6) R_p =(1-f)R_e \tag{1.6} Rp=(1f)Re(1.6) R p = R e 1 − e 2 (1.7) R_p =R_e\sqrt{1-e^2} \tag{1.7} Rp=Re1e2 (1.7) \qquad 比较式(1.6)和式(1.7)可得:
    f = 1 − 1 − e 2 (1.8) f=1-\sqrt{1-e^2}\tag{1.8} f=11e2 (1.8) e 2 = 2 f − f 2 (1.9) e^2=2f-f^2 \tag{1.9} e2=2ff2(1.9) \qquad 将椭圆方程两端对 x x x求导,并将式(7)带入得: 2 x R e 2 + 2 z ⋅ d z d x R e 2 ( 1 − e 2 ) = 0 (1.10) \frac{2x}{R_e^2}+\frac{2z\cdot\frac{dz}{dx}}{R_e^2(1-e^2)}=0 \tag{1.10} Re22x+Re2(1e2)2zdxdz=0(1.10) \qquad 整理得: d z d x = − ( 1 − e 2 ) x z (1.11) \frac{dz}{dx}=-(1-e^2)\frac{x}{z} \tag{1.11} dxdz=(1e2)zx(1.11) \qquad 式中 d z d x \frac{dz}{dx} dxdz表示椭圆在 P P P点的切线 P B PB PB的斜率。切线PB与法线PQ相互垂直,法线PQ的斜率为 t a n L tanL tanL,则有 d z d x t a n L = − ( 1 − e 2 ) x z t a n L = − 1 (1.12) \frac{dz}{dx}tanL=-(1-e^2)\frac{x}{z} tanL=-1 \tag{1.12} dxdztanL=(1e2)zxtanL=1(1.12) 由式(1.12)可得: z = x ( 1 − e 2 ) t a n L (1.13) z=x(1-e^2)tanL \tag{1.13} z=x(1e2)tanL(1.13) \qquad 将式(1.7)和式(1.13)带入椭圆方程,可求得以地理纬度为参数的椭圆方程: { x = R e 1 − e 2 s i n 2 L c o s L z = R e ( 1 − e 2 ) 1 − e 2 s i n 2 L s i n L (1.14) \begin{cases} x&=\frac{R_e}{\sqrt{1-e^2sin^2L}}cosL\\ z &=\frac{R_e(1-e^2)}{\sqrt{1-e^2sin^2L}}sinL\end{cases} \tag{1.14} {xz=1e2sin2L RecosL=1e2sin2L Re(1e2)sinL(1.14) \qquad 记线段长度 P Q ‾ = R N \overline{PQ}=R_N PQ=RN,则有: x = R N c o s L (1.15) x=R_NcosL \tag{1.15} x=RNcosL(1.15) \qquad 比较式(1.14)和式(1.15)可得: R N = R e 1 − e 2 s i n 2 L (1.16) R_N=\frac{R_e}{\sqrt{1-e^2sin^2L}} \tag{1.16} RN=1e2sin2L Re(1.16) \qquad 式(1.14)可化简为: { x = R N c o s L z = R N ( 1 − e 2 ) s i n L (1.17) \begin{cases} x&=R_NcosL\\ z &=R_N(1-e^2)sinL\end{cases} \tag{1.17} {xz=RNcosL=RN(1e2)sinL(1.17) \qquad 对于同一点的地理纬度和地心纬度是存在偏差的,或者说地理垂线和地心垂线之间是存在偏差的。对于地心纬度有 tan ⁡ φ = z x \tan\varphi =\frac{z}{x} tanφ=xz由式(1.13)可得: tan ⁡ φ = ( 1 − e 2 ) t a n L (1.18) \tan\varphi=(1-e^2)tanL \tag{1.18} tanφ=(1e2)tanL(1.18) 地理纬度与地心纬度间的偏差记作 Δ L = L − φ \Delta L=L-\varphi ΔL=Lφ,则有 t a n Δ L = t a n ( L − φ ) = t a n L − t a n φ 1 + t a n L t a n φ = t a n L − ( 1 − e 2 ) t a n L 1 + t a n L ( 1 − e 2 ) t a n L = e 2 s i n L c o s L 1 − e 2 s i n 2 L (1.19) \begin{aligned} tan\Delta L&=tan(L-\varphi)=\frac{tanL-tan\varphi}{1+tanLtan\varphi} \\ &= \frac{tanL-(1-e^2)tanL}{1+tanL(1-e^2)tanL} =\frac{e^2sinLcosL}{1-e^2sin^2L} \end{aligned} \tag{1.19} tanΔL=tan(Lφ)=1+tanLtanφtanLtanφ=1+tanL(1e2)tanLtanL(1e2)tanL=1e2sin2Le2sinLcosL(1.19) \qquad 若将 Δ L \Delta L ΔL e e e视为无穷小量,则式(1.19)可近似为: Δ L ≈ e 2 s i n L c o s L ( 1 + e 2 s i n 2 L ) = e 2 ( 1 + e 2 s i n 2 L ) 2 s i n 2 L (1.20) \Delta L\approx e^2sinLcosL(1+e^2sin^2L)=\frac{e^2(1+e^2sin^2L)}{2}sin2L \tag{1.20} ΔLe2sinLcosL(1+e2sin2L)=2e2(1+e2sin2L)sin2L(1.20) \qquad 进一步有: Δ L ≈ e 2 2 s i n 2 L ≈ f s i n 2 L (1.21) \Delta L \approx \frac{e^2}{2}sin2L\approx fsin2L \tag{1.21} ΔL2e2sin2Lfsin2L(1.21) \qquad 若取椭球扁率 f = 1 298.257 f=\frac{1}{298.257} f=298.2571,则在纬度 L = 45 ° L=45° L=45°处的地理纬度与地心纬度的偏差 Δ L = 11. 5 ′ \Delta L=11.5' ΔL=11.5

    2.曲率半径

    \qquad 在微分几何中,我们把曲率的倒数称为曲率半径,即 r = 1 k r=\frac{1}{k} r=k1曲率是指平面上一条曲线在特定点的弯曲程度。直线不弯曲,其各点的曲率均为0,圆上各点的曲率都一样,但小圆比大圆弯曲的厉害。
    \qquad 一条曲线 A B AB AB的平均曲率 K ‾ A B \overline{K}_{AB} KAB定义为动点沿AB从A移动到B时的切线转过的角度 ∣ Δ α ∣ |\Delta \alpha| Δα与AB的长度 ∣ Δ s ∣ 之 比 |\Delta s|之比 Δs,即 K ‾ A B = ∣ Δ α ∣ ∣ Δ s ∣ (2.1) \overline{K}_{AB}=\frac{|\Delta \alpha|}{|\Delta s|} \tag{2.1} KAB=ΔsΔα(2.1) \qquad 曲线AB上A点处的曲率 K A K_A KA定义为: K A = lim ⁡ P 沿 A B ‾ → A K ‾ A B = lim ⁡ Δ s → 0 ∣ Δ α Δ s ∣ (2.2) K_A=\lim_{P沿\overline{AB}\rarr A}\overline{K}_{AB}=\lim_{\Delta s\rarr0}|\frac{\Delta \alpha}{\Delta s}| \tag{2.2} KA=P沿ABAlimKAB=Δs0limΔsΔα(2.2) \qquad 假设函数 y ( x ) y(x) y(x)二阶可导,曲线上一点 ( x , y ( x ) ) (x,y(x)) (x,y(x))处的曲率 K = K ( x ) K=K(x) K=K(x)

    在这里插入图片描述
    \qquad 如上图所示,曲线C的方程为 y = y ( x ) y=y(x) y=y(x),点 M ( x , y ( x ) ) , M ′ ( x + Δ x , y ( x + Δ x ) ) M(x,y(x)),M'(x+\Delta x,y(x+\Delta x)) M(x,y(x))M(x+Δx,y(x+Δx)),点M处的切线倾角 α = α ( x ) \alpha=\alpha(x) α=α(x),则 t a n ( α ) = y ( x ) ′ tan(\alpha)=y(x)' tan(α)=y(x)
    \qquad 在曲线C上取一点 M 0 ( x 0 , y ( x 0 ) ) M_0(x_0,y(x_0)) M0(x0,y(x0))作为度量弧段的基点,弧 M 0 M M_0M M0M的长度为 s ( x ) s(x) s(x),则曲率定义为: K = K ( x ) = lim ⁡ M ′ 沿 C → M K ‾ M M ′ = lim ⁡ M ′ 沿 C → M ∣ α ( x + Δ x ) − α ( x ) s ( x + Δ x ) − s ( x ) ∣ = lim ⁡ Δ x → 0 ∣ [ α ( x + Δ x ) − α ( x ) ] / Δ ( x ) [ s ( x + Δ x ) − s ( x ) ] / Δ ( x ) ∣ = ∣ α ′ ( x ) ∣ ∣ s ′ ( x ) ∣ (2.3) \begin{aligned} K=K(x)&=\lim_{M'沿C\rarr M}\overline{K}_{MM'}=\lim_{M'沿C\rarr M}|\frac{\alpha(x+\Delta x)-\alpha(x)}{s(x+\Delta x)-s(x)}|\\ &=\lim_{\Delta x\rarr0}|\frac{[\alpha(x+\Delta x)-\alpha(x)]/\Delta(x)}{[s(x+\Delta x)-s(x)]/\Delta(x)}|\\ &=\frac{|\alpha'(x)|}{|s'(x)|} \end{aligned} \tag{2.3} K=K(x)=M沿CMlimKMM=M沿CMlims(x+Δx)s(x)α(x+Δx)α(x)=Δx0lim[s(x+Δx)s(x)]/Δ(x)[α(x+Δx)α(x)]/Δ(x)=s(x)α(x)(2.3) \qquad 因为 t a n ( x ) = y ′ ( x )    ⟹    α ( x ) = a r c t a n ( y ′ ( x ) )    ⟹    α ′ ( x ) = y ′ ′ ( x ) 1 + y ′ ( x ) (2.4) tan(x)=y'(x)\implies \alpha(x)=arctan(y'(x)) \implies \alpha'(x)=\frac{y''(x)}{1+y'(x)} \tag{2.4} tan(x)=y(x)α(x)=arctan(y(x))α(x)=1+y(x)y(x)(2.4) \qquad s ′ ( x ) = 1 + y ′ ( x ) 2 s'(x)=\sqrt{1+y'(x)^2} s(x)=1+y(x)2 及式(2.4)带入式(2.3)得 y = y ( x ) y=y(x) y=y(x) ( x , y ( x ) ) (x,y(x)) (x,y(x))处的曲率为: K ( x ) = ∣ α ′ ( x ) ∣ ∣ s ′ ( x ) ∣ = ∣ y ′ ′ ( x ) ∣ ( 1 + y ′ ( x ) 2 ) 3 2 (2.5) K(x)=\frac{|\alpha'(x)|}{|s'(x)|}=\frac{|y''(x)|}{(1+y'(x)^2)^{\frac{3}{2}}} \tag{2.5} K(x)=s(x)α(x)=(1+y(x)2)23y(x)(2.5) \qquad K = ∣ y ′ ′ ∣ ( 1 + y ′ 2 ) 3 2 (2.6) K=\frac{|y''|}{(1+y'^2)^{\frac{3}{2}}} \tag{2.6} K=(1+y2)23y(2.6) \qquad 由于地球近似为一个参考椭球,所以地球表面的曲率半径不是一个常数,而是随位置和方位的不同而发生变化。因此在惯导系统中进行经纬度解算时,必须考虑地球半径变化的情况。在求取地球表面某一点M的纬度时,首先根据北向速度 V N V_N VN求出纬度的变化率 L ˙ = V N R N \dot{L}=\frac{V_N}{R_N} L˙=RNVN,对其积分得到纬度 L L L
    \qquad 根据曲率半径的定义,可得地球子午面内主曲率半径 R N R_N RN为: R N = 1 K = [ 1 + ( d z e d x e ) 2 ] 3 2 d 2 z e d 2 x e (2.7) R_N=\frac{1}{K}=\frac{[1+(\frac{dz_e}{dx_e})^2]^{\frac{3}{2}}}{\frac{d^2z_e}{d^2x_e}}\tag{2.7} RN=K1=d2xed2ze[1+(dxedze)2]23(2.7) \qquad 其中 d z e d x e = − c o t L (2.8) \frac{dz_e}{dx_e}=-cotL \tag{2.8} dxedze=cotL(2.8) \qquad 由(2.7)得 d 2 z e d 2 x e = 1 s i n 2 L d L d x e (2.9) \frac{d^2z_e}{d^2x_e}=\frac{1}{sin^2L}\frac{dL}{dx_e} \tag{2.9} d2xed2ze=sin2L1dxedL(2.9) \qquad 由(1.14)得 d x e d L = ( 1 − e ) 2 R e t a n L s e c 2 L [ ( 1 − e ) 2 R e t a n 2 L + 1 ] 3 2 (2.10) \frac{dx_e}{dL}=\frac{(1-e)^2R_etanL sec^2L}{[(1-e)^2R_etan^2L+1]^{\frac{3}{2}}} \tag{2.10} dLdxe=[(1e)2Retan2L+1]23(1e)2RetanLsec2L(2.10) \qquad 将式(2.9)带入式(2.8)的: d 2 x e d 2 L = 1 s i n 2 L − [ ( 1 − e ) 2 R e t a n 2 L + 1 ] 3 2 ( 1 − e ) 2 R e t a n 2 L + 1 = − [ ( 1 − e ) 2 R e t a n 2 L + 1 ] 3 2 ( 1 − e ) 2 R e t a n 3 L (2.11) \frac{d^2x_e}{d^2L}=\frac{1}{sin^2L} \frac{-[(1-e)^2R_etan^2L+1]^{\frac{3}{2}}}{(1-e)^2R_etan^2L+1} =-\frac{[(1-e)^2R_etan^2L+1]^{\frac{3}{2}}}{(1-e)^2R_etan^3L} \tag{2.11} d2Ld2xe=sin2L1(1e)2Retan2L+1[(1e)2Retan2L+1]23=(1e)2Retan3L[(1e)2Retan2L+1]23(2.11) \qquad 展开 ( 1 − e ) 2 (1-e)^2 (1e)2项并略去 e 2 e^2 e2项得: R N ≈ ( 1 − 2 e ) R e ( 1 − 2 e s i n 2 L ) 3 2 ≈ R e ( 1 − 2 e + 3 e s i n 2 L ) (2.12) R_N\approx\frac{(1-2e)R_e}{(1-2esin^2L)^{\frac{3}{2}}} \approx R_e(1-2e+3esin^2L)\tag{2.12} RN(12esin2L)23(12e)ReRe(12e+3esin2L)(2.12) \qquad 卯酉圈的主曲率半径 R M R_M RM为: R M = x e c o s L = R e [ ( 1 − e ) 2 s i n 2 L + c o s 2 L ] 1 2 (2.13) R_M=\frac{x_e}{cosL}=\frac{R_e}{[(1-e)^2sin^2L+cos^2L]^{\frac{1}{2}}} \tag{2.13} RM=cosLxe=[(1e)2sin2L+cos2L]21Re(2.13) \qquad 展开 ( 1 − e ) 2 (1-e)^2 (1e)2项并略去 e 2 e^2 e2项得: R M ≈ R e ( 1 − 2 e s i n 2 L ) 1 2 ≈ R e ( 1 + e s i n 2 L ) (2.14) R_M\approx\frac{R_e}{(1-2esin^2L)^{\frac{1}{2}}}\approx R_e(1+esin^2L) \tag{2.14} RM(12esin2L)21ReRe(1+esin2L)(2.14)

    参考文献
    [1] 同济大学数学系. 高等数学.第7版[M]. 高等教育出版社, 2014.
    [2] 严恭敏, 翁浚. 捷联惯导算法与组合导航原理[M]. 西安:西北工业大学出版社, 2020.
    [3] 朱家海. 惯性导航[M]. 北京:国防工业出版社,2008.

    展开全文
  • 直动从动件凸轮廓线曲率半径的CAD动画计算.pdf
  • 手机和Matlab相结合测量牛顿环的曲率半径.pdf
  • 曲率、曲率圆和曲率半径

    千次阅读 2021-04-21 16:20:53
    一个圆半径越小,看起来就越弯曲;半径越大,看起来就越平;半径无限大,看起来就像一条直线,几乎不弯曲了。 —— 引言 曲率是描述曲线弯曲程度的量。...D为曲率中心,为曲率半径。 根据以上分析...

    一个圆半径越小,看起来就越弯曲;半径越大,看起来就越平;半径无限大,看起来就像一条直线,几乎不弯曲了。      —— 引言

    曲率是描述曲线弯曲程度的量。从上面左图可以看出,弯曲程度越大、转角越大;从右图可以看出,转角相同时,弧越短、弯曲程度越大

    如上图所示,曲线C在M点的曲率定义为:

    如上图所示,设曲线y=f(x)在点M处的曲率为,在点M处的曲线的法线上,在凹的一侧取一点D,使得:


    以D为圆心、为半径画圆,称此圆为曲线在点M处的曲率圆。D为曲率中心,曲率半径

    根据以上分析,圆上所有点的曲率相等都等于其半径的倒数

    展开全文
  • #资源达人分享计划#
  • #资源达人分享计划#
  • 波前差法测量焦距和曲率半径
  • %%%求最小曲率半径matlab源程序整个过程以下是对一组数据拟合后,并求其最小曲率半径!%%%%%%%源程序%第一组数据程序和结果%曲线拟合x=[0.001466076582203120,0.001466076582203120,0.007867945866887910,0....
  • matlab求曲率半径

    千次阅读 2021-04-18 10:29:54
    图 1 为 90 磁过滤管, 线圈缠绕的方法是内密外松, 计算中视内侧为紧密排列 , 线圈砸数为 N , 弯管管径 为 r , 曲率半径为 R; 图 2 为喇叭形螺线管,......干涉法测透视镜的曲率半径_理学_高等教育_教育专区。大学物理...
  • 通过输送机平面转弯段的受力分析,改进了平面转弯带式输送机转变处的布置结构,经使用,输送带在运行至弯曲段时能自动进行居中调节,效果较好。
  • 我在二维坐标系中把一段样条曲线等间隔分成360个点之后,可以输出一个点的坐标表格,我想问一下要怎么在matlab中读取表格中这些点并且依次计算每个点的曲率半径然后输出成一个曲率半径的数据文件呢(小白一枚,好多...
  • python_多点拟合曲线并计算曲率半径

    千次阅读 2021-08-03 22:59:46
    本文主要讲解:python_多点拟合二次曲线,再选取两边和中点拟合圆,计算曲率半径 主要思路: 使用numpy中的np.polyfit(x, y, 2)和 np.poly1d(f1)两个函数拟合二次函数 选取拟合的二次函数中的两边和中点拟合圆 利用...
  • 用matlab求出最小曲率半径曲线方程实例:y=(-5.324e-015*x.^7)+(-8.452e-012*x^6)+(-5.204e-009*x^5)+(-1.558 e-006*x^4)+(-0.0002253*x^3)+(-0.01515*x^2)+(-1.464*x)+222.9求在区间(0,2)之间的最小曲率半径。...
  • 曲率及其曲率半径及计算一、弧微分 二、曲率及其计算公式 例1 计算等双曲线x y ?1在点(1,1)处的曲率. 例2 抛物线y?ax2?bx?c 上哪一点处的曲率最大? 三、曲率圆与曲率半径 例3 设工件表面的截线为抛物线y?0.4x 2...

空空如也

空空如也

1 2 3 4 5 ... 20
收藏数 5,578
精华内容 2,231
关键字:

曲率半径