Physics Note 1
自己先读他一遍
力学
坐标系转化
例: 球坐标系
科里奥利力
[think] 这里的关键是,其实向量并不擅长表达旋转,所以表达出旋转坐标系下的某位置在真实系下的表达然后再求两次导是很麻烦的.但 一个量在旋转系和静止系下分别的导数的关系是简单的(指的是,同一个向量,选定旋转的基向量与固定的基向量的两种表示后得到的坐标求导的关系)
dtdvI=dtdvR+ω×v
其中I表示静止系下看,R表示在旋转系下看.
直接写成vI=(rcost,rsint)vR(0)+(−rsint,rcost)vR(1),求导得:
v˙I=(rcosωt,rsinωt)vR(0)˙+(−rsinωt,rcosωt)vR(1)˙+(rcosωt,rsinωt)˙vR(0)+(−rsinωt,rcosωt)˙vR(1)
整理一下就是v˙R+ω×v了.ω顺时针向内逆时针向外,遵循右手螺旋定则.
于是就可以推了:
r˙I=r˙R+ω×rRr¨I=(r¨R+ω˙×rR+ω×r˙R)+ω×(r˙R+ω×rR)⇒aI=(aR)+(ω˙×rR)+(2ω×vR)+(ω×(ω×rR))
其中第一项是原来的加速度,第二项是如果非匀速旋转导致的,第三项是科里奥利力导致的,第四项是离心力导致的(ω×(ω×rR)方向朝内,大小是ω2∣rR∣).
引力的球壳公式
怎么喵的推起来这么麻烦?分出去开篇新文章了.
狭义相对论
首先,他是要处理在两个有匀速相对运动的参考系下的变换,且由于两个参考系等价,这个变换应该是对称的,还要做到不改变物理规律.它解决的变换是B系是怎么样的,那么你在A系看B系是什么样的.
同时性原理:就是同一点同时发生的两件事在另一个系里也应该同时发生.(其实感觉这个也可以是保证因果链,避免果在因之前发生,看起来和保时序是等价的).我们很快发现它保障了可以定义一个点的时间.
光速不变原理:光速在每个系看起来是一样的.
然后一个推导过程是这样的:考虑B系相对A系有一个(v,0)的速度.
那么我们在B系(0,0),(0,1)分别放一个镜子,并让里面有一束光在这两个点之间来回,这相当于一个运行在(0,0)处的钟,我们叫他光钟.
那么根据同时性原理我们可以在0处进行计时:如果你在t1时刻和t2时刻接收到这个光,那么一件发生在[t1,t2]的事就有了时长属性t2−t1,这也是为什么它真的是一个钟.或者说,当我们想变换时间的时候我们需要先知道时间是什么,那么不妨定义成光钟计算出的时长.
那么我们发现,在B系下,显然每隔t1=c2l接到一次;而A系下,光速不变,但光速现在是斜着走的了,所以光速就成了(v,c2−v2),竖直分量笑了,就成了每隔t2=c2−v22l,于是你就算出了时间的变换系数是γ=t1t2=1−(cv)21.
注意方向:t1t2是说同一件事,在里面经过了t1的时间,则在外面经过了t2时间,因为t2>t1,所以我们说钟慢了.
然后你也可以再去考虑距离:因为光速相同,由速度的定义,相同时间里走过的路程也相同,所以你把光钟横着放,考虑光钟的长度即可.
然后你还是考虑长度为x的光钟,在B系下显然是要用t1=c2x1,在A系下由于光速不变,系在跑,就成了t2=v−cx2−v+cx2,令这两个时间相等得到t2=γt1,化简得x2=γ1x1
最后你再考虑不同位置的两个点的时间怎么算?即:B系下相隔x处的两个点的时间差:考虑在2x处同时向两侧发两束光,那么它俩在B系下同时到达,这是一个B系下的同时的事件(但不同地点),而到A系看的话:你先把整个长度2x变换过去成2γx,然后光在以c的速度跑,车在以v的速度跑,你就算出来时间是2γx(c−v1−c+v1)=γc2vx.
你用这三部分就可以推任意变换了:你认为B系的原点的时间是这个系的时间,那么变换(x,t)处的事件,自然就变成:(γx,γt−c2γvx)
如果你x这一维是沿运动方向为正的,那么这个Δt就该减上去了.
最后再简单说几个点:
- 为什么垂直于运动方向的时间和空间都不会变?因为你假设会变,那考虑两个圆环沿轴向彼此接近,它俩互相看对方都比自己大/小,就说不清谁会套到谁了.
- 如果一维空间?那很糟了,不过结论不会变,你可以用A和B之间互相发光推.
感觉把相对论拆成空间,时间,不同空间位置的时间差这三块看比洛伦兹更好记罢了()
功和能
动能定理
F=dtd(mv)∫Fdx=∫dtd(mv)dx=m1∫mvd(mv)=21mv2∣t0t1
保守力与势能与机械能守恒
与路径无关的力可以构造出势能的标量场.使得两点之间运动时该力的功等于势能差.势能场的梯度是力场.
我们看成保守力的势能只能和动能等量转化,一个质点的机械能守恒其实是每个 (保守力+动能守恒)相加 的结果.
稳定平衡点与不稳定平衡点
总势能场求梯度是力场,所以在势能场的极值点处做泰勒展开(物理当然假设无限可微)取一阶项,就知道平衡点处微扰后是类似简谐的稳定在平衡点附近还是类似正反馈偏离平衡.
质点系功能原理
质点系总能量(动能+保守力势能)变化量等于非保守力和外力的做功.
刚才你把保守力视为力场的时候,你是固定了力的一端看另一端,才能产生静态的力场.
现在要考虑两端都在动的情况:不过你发现没有关系,由作用力反作用力的关系,你能知道:
∫Fdx1−∫Fdx2=∫F(dx1−dx2)
对质点来说,你的力场不会旋转,那么其中dx1−dx2恰好是你固定这个力场在一个点后另一个点的相对位移.
但这么想不是很优雅:你需要破坏质点之间的对称.有没有更好的方法呢?
AI大人说,你应该把势能看成关于两者坐标的函数E(x1,x2),然后从空间的平移对称性可以说明E必须只依赖x1−x2,从空间的旋转对称性,可以说明只依赖x1−x2的长度,从而任意两个点的势能都只和相对距离相关.
在这里对E求梯度还可以得到保守力一定是沿两质点方向的结论.(保守力一定是有心力.)
此时你的动能定理变成了:
∫F1dx1−F1dx2=E(x1,x2)−E(x1+Δx1,x2+Δx2)=E(∥x1−x2∥)−E(∥x1+Δx1−x2−Δx2∥)
现在同时对所有点对累加,一侧累加了所有内力的做工,一侧是系统总势能的变化.
这时候你又有了些疑问:无限个点的时候我直接用质量密度积起来凭啥是对的?然后你想了一下发现我们在学物理,所以E和F都关于点的位置连续,而密度必须连续且有界,又因为我们在高斯定理那里考虑过弱有界定理的事,所以对于空间中二次方的力都是可积的(因为dV关于半径是三次方),这是不是甚至保证了绝对可积可以交换积分号啊.
这些性质是普遍满足的,所以我们看起来总可以在次数不超过2的时候乱换积分号.
而这里你只是把做功加起来,势能加起来,不需要交换积分,所以就合法了.
动量
动量定理
F=dtd(mv)⇒Fdt=d(mv)⇒∫Fdt=∫mdv=mΔv
动量守恒
你在学质点系功能定理的时候意识到力还是对点对考虑比较好.而对每个点对显然m1dv1+m2dv2=Fdt−Fdt=0,然后你仍然是都累加起来得到动量守恒.
质点系动量定理
内力已经动量守恒了,然后因为你Fdt=dtd(mx)是线性的,所以你可以全累加起来,然后要交换求和或积分号,得到:
∫∑Fdt=∫∑midvi∫∑mvdt=∫∑mΔxi,1−xi,0=M[(∑MmiΔxi,1)−(∑MmiΔxi,0)]where M=i∑mi
第一个式子是合外力冲量等于总动量变化量
第二个式子告诉你对动量积分得到于质心的位移
恢复系数
定义为碰撞前后的相对速度比:(v1−v2)=e(v1′−v2′)
{mv1+Mv2=mv1′+Mv2′v1−v2=e(v1′−v2′)⇒⎩⎨⎧v1′=v1−m1+m2m2(1+e)(v1−v2)v2′=v2−m1+m2m1(1+e)(v1−v2)
角动量
角动量定理
定义力矩M=r×F
角动量L=r×p=mr×v.
感觉拿面积建立直觉还不如速度乘旋转半径好理解呢
dtdL=dtdr×p+r×dtdp=0+r×F=M
(发现乘法法则对任意双线性的东西都是类似的).
自然你会得到当没有力矩的时候角动量守恒,也就是开普勒第二定律了.
天体椭圆轨道
选定质心为原点,则无外力的情况下它是惯性系.
只有引力,过质心,所以角动量守恒,注意角动量向量有方向的垂直于旋转面,从而是平面运动,且天体满足r×r˙=mL,L为定值.
同时你可以拿机械能守恒再列一个−∥r∥GMm+21mr˙2=E为定值.
然后你把r的分量写开肯定是两个未知数两个方程,然后把一个变量消掉就开始积分吧()由于太麻烦由于我们在学物理所以省略这些神秘的积分步骤,你会得到它是椭圆.
刚体定轴旋转
三维刚体在定轴旋转时等价于压扁了的二维刚体.所以下面都研究二维刚体.
L=∫x∈Vx×(ρ(x)dV)(ω×x)x×(ω×x)=x2ω⇒L=∫x2ωρ(x)dV=ω∫x2ρ(x)dV=ωI
这里I就是转动惯量.单位是m2⋅kg.
因为是二维的,所以x⊥ω可以直接乘.
E=∫x∈V21(ρ(x)dV)(ω×x)2=21(∫x∈Vρ(x)dVx2)ω2=21ω2I
转动惯量之间有关系,方便计算
平行轴定理:已知过质心的l的转轴对应的转动惯量为Il,对与他平行的距离为d的轴m有:
Im=Il+Md2
Im=∫x∈V((x×ul)−d)2ρ(x)dV=∫x∈V(x×ul)2+∫x∈Vd2ρ(x)dV−2∫x∈V(x×ul)d=Il+Md2−2∫x∈V(x×ul)d
其中ul是l的方向单位向量.注意到因为原点取质心,所以最后一项是0,得证.
对二维刚体来说,垂直于平面的轴转动惯量等于平面上两条相交于它的垂直的轴的转动惯量之和:
Iz=Ix+Iy
振动与机械波
简谐运动的合成
物理好像莫名其妙喜欢考三角函数等差数列求和.也就是:
i=L∑Rcos(aix+b)=i=L∑R2sin(2ax)sin(ax(i+21)+b)−sin(ax(i−21)+b)=2sin(2ax)1(sin(ax(R+21)+b)−sin(ax(L−21)+b))=sin(2ax)sin(2(R−L+1)ax)cos(2(L+R)ax+b)
阻尼/受迫振动
微分方程求解的一集.
阻尼振动:
mx¨=−kx−fx˙
常系数齐次线性微分方程.解法是设解erx代入然后比较系数把r解出来:
mr2+fr+k=0Δ=f2−4mkr=2m−f±f2−4mk
从而Δ>0是过阻尼,解出来两个指数衰减的.Δ=0是临界阻尼.此时xerx也是一个特解.Δ<0是欠阻尼,此时不一定两个系数相同,因为是共轭的怎么取都是实数,所以是会产生相位的,且A是吸收了一个常数2,辅助角公式的系数等的,用欧拉公式展开即得公式:
ω=2m4mk−f2x=Ae2m−fxcos(ωx+φ)
受迫振动
变成:
mx¨+fx˙+kx=Fcosωt
改成解mx¨+fx˙+kx=0的通解和原方程的任意一个特解.而特解通过复数:
mx¨+fx˙+kx=Feiωtx=Aeiωt⟹A=−ω2m+iωf+kF
此时可以方便的通过∣A∣得到稳定时的振幅,通过k−ω2m+ωfi的幅角的相反数得到相位.就是稳定时的解.
绳子上的波
这个是讲如何建模一个振动的绳子.
设绳子的力是T,绳子的线密度(m/l)为μ,则波速是什么?
设波的方程是f(x,t),则对一个x0处的小段dl,其只受左右两侧的质点对他的拉力,那么比如右侧对他的拉力的竖直分量是Tsinθ,但我们难以知道sinθ,但可以方便的知道sinθ∼tanθ=∂x∂f(x+dl,t),而左侧是∂x∂f(x,t),所以你得到纵向的力是:
T(∂x∂f(x+dl,t)−∂x∂f(x,t))=Tdl(∂x)2∂2f=dlμ(∂t)2∂2f
解出来就得到波速是:
v=μT
波的能量,干涉
多普勒效应