首页/文章/ 详情

圆锥曲线和天体运动,以及微分方程

10月前浏览427

圆锥曲线和天体运动,以及二阶线性偏微分方程

什么是圆锥曲线

如图所示,圆锥曲线,是由一个圆锥被一个平面切割所产生的各种曲线。 圆锥是由一条与轴线保持固定夹角 α 的直线绕该轴线旋转一周而形成的立体图形。 当用不同角度的平面去切割这个圆锥,横截面与圆锥轴线之间的夹角 β,直接决定了我们得到哪种曲:

  • 当 0° ≤ β < α 时:平面与圆锥的两叶都相交,其交线是一条开放且有两支的双曲线(如图中橙色面所示)。

  • 当 β = α 时:平面恰好与圆锥的一条母线平行,此时得到的是一条开放的抛物线(如图中蓝色面所示)。

  • 当 α < β ≤ 90° 时:平面只与圆锥的一叶相交,交线是一个封闭的椭圆(如图中绿色面所示)。

  • 当 β = 90° 时:这是一个非常特殊的情况,截面与圆锥的轴线垂直,此时得到的便是一个完美的圆(如图中红色面所示)。我们可以将圆视为椭圆的一个特例。

因此,双曲线、抛物线、椭圆和圆,这四种看似形态各异的曲线,实际上本质是统一的,它们都是平面与圆锥相交的所产生的曲线,统称为“圆锥曲线”。

关于x和y的一般二元二次方程

一般形式可以写作:

 

将该方程写成矩阵形式:

 

即:

 

其中    是实对称矩阵,他的行列式为:

 

将    作为判别式

  • 当      时,方程变为椭圆形(退化情况有点椭圆虚椭圆)
  • 当      时,方程变为抛物线(退化情况有两条平行直线)
  • 时,方程变为双曲线(退化情况有两条相交直线)

天体运行物理规律

  • 当天体的动能不足以克服引力源的引力势能时(总机械能<0)
    • 天体在一个封闭的,周期性的轨道上做椭圆运动(条件极其苛刻的情况下做圆周运动几乎不存在)
    • 例如地球围绕太阳,月球围绕地球,哈雷彗星光顾地球
  • 当天体的动能恰好够克服引力源的引力势能时(总机械能=0)
    • 由于条件极其苛刻因此几乎不存在,但是在地球上理想情况下的抛物运动总是抛物线轨迹
  • 当天体的动能大于引力源的引力势能时(总机械能>0)
    • 天体从某个方向飞来,然后在中心天体的引力作用下发生偏转(引力弹弓效应)然后朝另一个方向飞去做双曲线运动
    • 例如星际彗星‘奥陌陌’,最近的‘阿特拉斯’以及人类发射的星际探测器

二阶线性偏微分方程

一般形式可以写作:

 

类似于一般的二元二次方程,其判别式:

 
  • 当      时,对应的是一个椭圆方程,他是封闭、有限的图形。椭圆上没有特定的“方向性”,它是匀称的。
       最典型的例子是拉普拉斯方程    
 
    • 椭圆型方程没有实的特征线。这意味着信息传播速度是无穷大的。
    • 其描述的是平衡态、稳态问题,如稳定的温度场、静电场。某一点的值受区域内所有点的影响,同时也影响所有点,没有时间维度和因果关系。这类似于椭圆的“封闭性”和“整体性”。
  • 当      时,对应的是一个抛物线方程,他是开放、无限的,但它只有一个开放方向,具有独特的对称轴,指明了主导方向。
       最典型的例子是热传导方程    

       
     
    • 有一族实的特征线
    • 描述的是扩散过程,如热传导、粘性流体的扩散。它描述的是随时间单向演化的系统。信息理论上可以瞬间传遍全域,但强度急剧衰减(比如在远处加热,近处只会感受到微弱的温度变化)。这类似于抛物线的“单向开放性”和“主导方向”。
  • 时,对应的是一个双曲型方程,有两支,是开放、无限的,并且有两个明确的渐近方向。
       最典型的例子是波动方程    

       
     
    • 有两族实的特征线
    • 描述的是波动现象,如声波、振动。信息以有限的速度 c沿特征线传播,具有清晰的因果律。某一点的扰动只影响其未来的一个锥形区域(影响域)。这完美地类比了双曲线的两支,代表了波前向前和向后传播的方向。

总结

圆锥曲线基于几何形状(不变量)分类,偏微分方程基于解的特性(不变量)分类,却具有可以类比的性质,相同的判别式,说明截然不同的物理现象,完全可能共享相同的数学本质

后续将分享神经网络解微分方程的文章,本账号做分享的全部文章代码后续都会放在知识星球,如有需要请关注我吧~~~🤞🤞🤞


来源:我的博士日记
振动电场理论
著作权归作者所有,欢迎分享,未经许可,不得转载
首次发布时间:2025-10-26
最近编辑:10月前
此生君子意逍遥
博士 签名征集中
获赞 62粉丝 130文章 153课程 0
点赞
收藏
作者推荐

利用奇异值分解求解率无关晶体塑性本构

参考文献:《A computational procedure for rate-independent crystal plasticity》文章doi:10.1016/j.tws.2024.112610在速率无关晶体塑性中,长期悬而未决的三件事:(i) 如何判定哪些滑移系“活动”;(ii) 活动滑移系上的剪切增量是多少;(iii) 多滑移导致的非唯一性如何稳定、唯一地求解。作者提出一个稳健、唯一的计算流程来同时解决三点,并证明它与低速率敏感度的速率相关模型( )在预测上基本一致。详细的算法总结如下:作者使用该方案,对 fcc 单晶三种典型取向([001]、[111]、[236]),无硬化和强潜在硬化( q=1.4q=1.4)两类情形进行了数值算法的对比,其数值结果表明速率无关求解与速率相关( )预测几乎不可区分;同时多晶(343 晶粒)压缩及 Taylor 模型亦吻合良好。同时率相关和率无关模型模拟的多晶的取向演化特征也基本无差异。此外作者的通过SVD奇异值分解方案的率无关数值框架的计算效率相对率无关高一些,而作者对于umat隐式实现方案的一致性雅可比通过引入小的率相关系数实现的近似雅可比方案。显示vumat则不需要雅可比。基于作者数值方案,这里使用SVD奇异值把率相关模型改为对应的率无关模型(采用显示vumat进行实现)。实现效果如下:应力分布结果(变形5%时):活跃的滑移系同个数(变形5%时):迭代过程中的最大残差(变形5%时): 在这个数值实现过程中,vumat的率无关是略快于经典的率相关求解方案的(大概1.5倍左右),同时计算结果与率相关方案(低率相关系数)的应力分布和整体的流动应力结果几乎没有区别。对率无关模型感兴趣的可以阅读这篇经典的率无关算法实现的方案。来源:我的博士日记

未登录
还没有评论
课程
培训
服务
行家
VIP会员 学习计划 福利任务
下载APP
联系我们
帮助与反馈