计算机模拟物理作业1.docxVIP

计算机模拟物理第一次作业工物02黄丽达 2010011686工物02 王岩 2010011707工物02 胡丕丕 2010011681行星运动方程的描述在行星运动过程中受到太阳对其的引力作用，在直角坐标系下，其受力如（1-1）式：（1-1）其中，外有引力常量，M为太阳质量，，为指数偏离度。1.1 地球绕太阳运动在天体的实际运动中，空间尺度和时间尺度都非常大，因此，取地球到太阳的平均距离为单位长度（1AU），地球绕太阳运动的周期时间为单位时间（1year），记：A=1AU=1.496x1011mB0=1year=365x24x3600s=3.154x107sB0=1year是平方反比的中心力场情况下地球做圆周运动的周期，如果引力不满足平方反比规律，有一个偏离，此时的周期为B，单位为年（year）。当使用A和B作单位时，记x=XA，y=YA，t=TB，此时有（1-2）将（1-2）代入（1-1）得到：（1-3）记，采用新的单位制后，运动方程为：（1-4）二行星运动记：A=1AU=1.496x1011m，B=1year=365x24x3600s=3.154x107s，不考虑两行星之间的相互作用时，每个行星的运动单独考虑。（1-5）考虑两行星之间的相互作用时，对行星1，加速度分量如下：（1-6）记（1-7）同理，对行星2：（1-8）记（1-9）作业采用Visual Studio 2012编程生成轨道的数据，用matlab作图。地球绕太阳运动数值模拟实验太阳坐标取为（0,0）。选择合适的初速度，使行星运动为圆周运动。用Euler-Cromer方法模拟。用Euler方法结果如何。选取初速度时，我们假定地球绕太阳作圆周运动，则其速度满足使用Euler-Cromer方法，时间间隔设为0.1，模拟10个周期，得到的轨道如图2-1：时间间隔设为0.01，模拟10个周期，得到的轨道如图2-2：时间间隔设为0.001，模拟10个周期，得到的轨道如图2-3：使用Euler方法，时间间隔设为0.1，模拟10个周期，得到的轨道如图2-4：时间间隔设为0.001，模拟10个周期，得到的轨道如图2-5：时间间隔设为0.00001，模拟10个周期，得到的轨道如图2-6：综上所述，行星绕地球运动的轨道模拟结果为圆形。模拟时时间间隔越小，轨道越接近圆形；Euler-Cromer方法与Euler方法比较，相同时间间隔的情况下，前者的模拟轨道更接近圆形，即结果更接近真实值。如果行星间的吸引力不是而是，其中是小量，取为0.05.取行星在t=0时的位置及速度为，什么会发生？轨道会不会重叠？说明这不是由Δt的选取所致。行星会飞离或落入恒星吗？根据上问的结果，我们选用Euler-Cromer方法进行模拟。时间间隔设为0.001，模拟10个周期，得到的轨道如图2-7：时间间隔设为0.001，模拟1000个周期，得到的轨道如图2-8：时间间隔设为0.0001，模拟10个周期，得到的轨道如图2-9：行星运动的轨道不再是圆，而是在一个圆环形的区域内运动，其轨道会发生重叠。由2-7和2-9可知，这不是Δt的选取导致的，行星不会飞离或落入恒星。如果吸引力，先分析当，u=0时的v值。Δt必须是多少才能使计算的圆轨道在数圈内维持圆周？改变前面所确定的v大小2%，所得的轨道有何性质？此时行星受力满足，，即初速度时间间隔设为0.01，模拟10个周期，得到的轨道如图2-10：时间间隔设为0.001，模拟10个周期，得到的轨道如图2-11：时间间隔设为0.0001，模拟10个周期，得到的轨道如图2-12：时间间隔设为0.00001，模拟10个周期，得到的轨道如图2-13:改变v的大小，使v增大2%，v=6.4056时间间隔设为0.00001，模拟10个周期，得到的轨道如图2-14:改变v的大小，使v减小2%，v=6.1544时间间隔设为0.00001，模拟10个周期，得到的轨道如图2-15:综上，如果吸引力，初速度仍为6.28（天文学量纲），Δt至少小于0.001时，轨道才能在数圈内维持圆周。v增大2%，行星将向远处飘移逐渐远离太阳，v减小2%，行星将逐渐接近太阳。二行星运动数值模拟实验二行星m1、m2绕太阳M运行，m1/M=0.001，m2/M=0.01。取Δt=0.001，太阳坐标为（0,0）；（x1,y1）=（1,0），（u1,v1）=（0，）；（x2,y2）=（41/3,0），（u2,v2）=（0，）化为天文单位，v1=6.2809，v2=4.9581不考虑行星相互作用时，描述行星的轨道形状和周期。时间间隔0.001，模拟10个周期，得到的轨道如图2-14，两个行星的轨道都是以太阳为圆心的圆。其中红色为太阳的位置，蓝色的圆是行星1的轨道，周期为1B，黑色的圆是行星2的轨道，周期为2B