diff --git a/README.md b/README.md index 763e275..94d9b9e 100644 --- a/README.md +++ b/README.md @@ -80,8 +80,32 @@ Dose Estimation by Simulation of China Space Station

AE8/AP8模型结果

AE9/AP9模型结果
+### 辐射环境拟合 +#### 银河宇宙射线 +[^4] + +#### 太阳辐射 +[^5] + +#### 俘获辐射带 +1. 电子 + 根据文献，在轨道较高（$L>2.5$）的内带中，电子能谱中有$\sim64\%$符合指数分布，$\sim4\%$符合幂律分布。[^6] + + 但在本例中，轨道较低（$L\sim1.06$），因此使用指数分布$j = j_0e^{-\frac{E}{E_0}}$和幂律分布$j=j_0E^{-\alpha}$进行拟合，随后使用`Origin Lab`对两者进行拟合比较，最终选用指数分布模型$j_0=1.333680\times10^6,\ E_0=0.0824$。 + +
+ + | Model | $R^2$ | BIC | BIC diff | + |:-----:|:-------:|:---------:|:--------:| + | Ex | 0.98107 | 402.50208 | 0 | + | Pow | 0.94279 | 423.51296 | 21.01088 | + +
+ +2. 质子 + ## 空间站结构 -1. 尺寸与分区[^4] +1. 尺寸与分区[^7]
* 全长：16.6 m @@ -101,8 +125,8 @@ Dose Estimation by Simulation of China Space Station 2. 材料 * 外壳： 3层 2mm 铝 + 10mm 芳纶 + 5mm 铝
- * 5系铝合金[^5] - * 泰普龙[^6] + * 5系铝合金[^8] + * 泰普龙[^9]
* 填充：金属为主，使得总重大致相当 * 内部：空气 @@ -126,6 +150,9 @@ Dose Estimation by Simulation of China Space Station [^1]: [中国空间站轨道参数](http://www.cmse.gov.cn/gfgg/zgkjzgdcs/) [^2]: 程彭超,闵锐.近地空间辐射环境与防护方法概述[J].辐射防护通讯,2017,37(01):14-21. [^3]: Bourdarie S, Xapsos M. The near-earth space radiation environment[J]. IEEE transactions on nuclear science, 2008, 55(4): 1810-1832. -[^4]: [【知识点·航天】“天和”核心舱、“天宫”空间站和新一代载人飞船的最新知识点（干货版）](https://zhuanlan.zhihu.com/p/103709953) -[^5]: [中铝造为“天和”号核心舱披“铠甲”壮“筋骨”](https://m.thepaper.cn/baijiahao_12484370) -[^6]: [泰普龙产品相关知识](https://wenku.baidu.com/view/a2ecf93501d8ce2f0066f5335a8102d276a261cd.html) \ No newline at end of file +[^4]: +[^5]: +[^6]: Zhao H, Johnston W R, Baker D N, et al. Characterization and evolution of radiation belt electron energy spectra based on the Van Allen Probes measurements[J]. Journal of Geophysical Research: Space Physics, 2019, 124(6): 4217-4232. +[^7]: [【知识点·航天】“天和”核心舱、“天宫”空间站和新一代载人飞船的最新知识点（干货版）](https://zhuanlan.zhihu.com/p/103709953) +[^8]: [中铝造为“天和”号核心舱披“铠甲”壮“筋骨”](https://m.thepaper.cn/baijiahao_12484370) +[^9]: [泰普龙产品相关知识](https://wenku.baidu.com/view/a2ecf93501d8ce2f0066f5335a8102d276a261cd.html) diff --git a/assets/trapped-9.csv b/assets/trapped-9.csv new file mode 100644 index 0000000..e7c3ee2 --- /dev/null +++ b/assets/trapped-9.csv @@ -0,0 +1,4 @@ +0.1, 0.2, 0.4, 0.6, 0.8, 1, 2, 4, 6, 8, 10, 15, 20, 30, 50, 60, 80, 100, 150, 200, 300, 400, 700, 1200, 2000 +170.46630985001116, 148.7169777853078, 116.02958215502514, 91.68016573576256, 75.4004563167008, 52.44979368353229, 27.39576993169936, 11.534263356467031, 4.927074386668564, 2.503935512743802, 1.2207035405229754, 0.6406616768835681, 0.32153683341049205, 0.18332837912910016, 0.11496359892828469, 0.08025880550325193, 0.05970056915795244, 0.0384373301355559, 0.022821770094163085, 0.011082462960424989, 0.0047728913059455905, 0.0007325376524683292, 5.947420364867388e-05, 1.918902440238389e-07, 0.0 +0.04, 0.07, 0.1, 0.25, 0.5, 0.75, 1, 1.5, 2, 2.5, 3, 3.5, 4, 4.5, 5, 5.5, 6, 6.5, 7, 8.5, 10 +781425.9929592956, 679139.2146049731, 322379.40275874606, 56156.789327450584, 6157.486909985416, 494.7242655747493, 77.04500465644205, 17.196451096282715, 5.771096902928527, 1.655431341146984, 0.5327835060623752, 0.17776818060047872, 0.05960754201838149, 0.02284581379206805, 0.008191147518478015, 0.0034509863771389026, 0.0010065728536974858, 0.00025474598119156806, 4.308092673195076e-06, 0.0, 0.0 \ No newline at end of file diff --git a/fit.opju b/fit.opju new file mode 100644 index 0000000..04c8708 Binary files /dev/null and b/fit.opju differ diff --git a/utils/ae9ap9.py b/utils/ae9ap9.py index 8763efa..90a0bfa 100644 --- a/utils/ae9ap9.py +++ b/utils/ae9ap9.py @@ -24,8 +24,6 @@ _, ax = plt.subplots(1, 2, dpi=150, figsize=(16, 6)) ax1, ax2 = ax[0], ax[1] -print(proton) - ax1.plot(E_proton, proton) ax1.set_ylim([1e-3, 5 * 1e2]) ax1.set_yscale('log')