MCNP快速入门1分析.pptVIP

问题类型（MODE）卡 如果不给出MODE卡，则缺省形式是MODE N，即缺省值是中子输运问题。 格式： MODE x1 … xi xi ＝ N，中子输运。 P，光子输运。 E，电子输运。 几何卡 几何卡有以下几类： 助记符 卡片类型 VOL 栅元体积 AREA 曲面面积 U Universe TRCL 栅元变换 LAT 栅格 FILL 填充卡 TR 坐标变换 VOL 体积卡（可选） 体积卡用于输入各个栅元的体积。 格式： VOL x1 x2 … xi 或 VOL NO x1 x2 … xi xi ＝ 栅元 i 的体积，i ＝1，2，…，栅元总数 NO ＝ 不计算栅元的体积和面积 缺省：MCNP 将计算所有栅元的体积，除非在体积卡上出现“NO”。如果在体积卡上没有填写某一栅元的体积，则将使用计算的体积值。 在体积卡上，如果输入的项数不等于栅元总数，将产生致命错误。可以用 nJ 格式跳过不想输入体积的栅元。 MCNP 在计算栅元体积的同时，也计算栅元的质量和面积。 AREA 面积卡（可选） 面积卡用于输入各个曲面的面积。 格式： AREA x1 x2 … xi xi ＝ 曲面 i 的面积，i ＝1，2，…，曲面总数 缺省：MCNP 将计算所有曲面的面积。在面积卡上没有填写的面积，将使用计算的面积值。 记录计数所需的面积也可以由 SDn 卡提供。如果 AREA 卡和 SDn 卡都未提供记录计数所需的面积，而 MCNP 又无法计算该面积，将产生致命错误。 记录计数时，通常要用 SDn 卡提供的面积。当记录面积是整个曲面时，也可以用 AREA 卡提供的面积；否则，如果只记录一段曲面的计数，则只能用 SDn 卡提供的面积。 例子：如图结构，三个栅元 （一个外部的立方体内 包含一个内部的立方体， 中心是一个球） 栅元3描述如下： 3 0 8 -9 -10 11 -12 13 #2 #1 则MCNP不能计算栅元3的体积。 如果用12个面来描述，则能够 计算其体积。 TRn 坐标变换卡 TRn O1 O2 O3 B1 B2 B3 B4 B5 B6 B7 B8 B9 M n ＝变换号， 1≤ n ≤999 O1，O2，O3 ＝坐标变换向量的位移。 B1至B9 ＝坐标变换的坐标旋转矩阵。 M ＝1， 表示位移是辅助坐标系原点 相对于基本坐标系的位移。 ＝-1，表示位移是基本坐标系原点 相对于辅助坐标系的位移。 缺省值为：TRn 0 0 0 1 0 0 0 1 0 0 0 1 1 *TRn表示 Bi是角度而非角度的余弦，Bi的角度范围从0°至180° ，其排列顺序如下： 元素 B1 B2 B3 B4 B5 B6 B7 B8 B9 轴 x,x’ y,x’ z,x’ x,y’ y,y’ z,y’ x,z’ y,z’ z,z’ Bi的意义与M无关。用户不必给出全部的Bi值，使用nJ 格式可以跳过不用的项。Bi可按以下五种方式中的任意一种方式指定： 给出全部 9 个元素。 给出矩阵中列或行方向上的两个向量（6个值），MCNP将用矢积产生矩阵的第三个向量。 给出矩阵中列和行方向上各一个向量（5个值），其中的公共分量必须小于1，MCNP将用欧拉角方法完成该矩阵。 给出矩阵中列或行方向上的一个向量（3个值），MCNP将用某种任意的方法生成其它两个向量。 不给出Bi的值， MCNP将生成单位矩阵。 一个问题中最多可有999个坐标变换卡。 MCNP所有类型的曲面都具有某种对称性， 当对称轴平行于某坐标轴时，即使是较复杂的曲面都不难给出其描述的参数。 但是，当曲面对称轴不平行于坐标轴时，给出描述曲面的参数难度就大多了。 在这种情况下，坐标变换功能就显得极为有用。用户可以在一个以曲面对称轴为坐标轴的辅助坐标系中描述该曲面，然后建立该辅助坐标系与基本坐标系之间的坐标变换关系，通过TRn卡将该曲面的描述转换到基本坐标系中。 例子：高14cm、半径4cm的圆柱体，中心在(0,10,15)，对称轴在YZ平面，与Y轴成30度角，如图所示。 则用(x,y,z)坐标系来描述曲面就比较简单明了。 1 1 CY 4 2 1 PY -7 3 1 PY 7 TR1 0 10 15 1 0 0 0 .866 .5 0 -.5 .866 或 *TR1