MCNP所有类型的曲面都有某种对称性,当对称轴平行于某些坐标轴时,即使较为复杂的曲面都是不难给出描述参数的。但是,当曲面对称轴不平行于坐标轴时,对曲面给出描述参数难度就大多了。在这种情况下,用户可以在一个较为方便的辅助坐标系来描述该曲面,然后使用如下变换功能,将该曲面的描述转换到问题的基本坐标系中。这个变换是通过坐标变换卡TRn卡实现的。TRn卡指明辅助坐标系与基本坐标系的关系,在辅助坐标系描述曲面是容易的,但这时在该曲面描述卡上要有变换号n,在相应的TRn卡指定这个变换。环面(Torus)不能做坐标变换。
在MCNP的输入文件中,基本上是在辅助坐标系下给出较为简单的曲面描述卡。当辅助坐标系并不相应于基本坐标系时,便需使用TAn卡,用以指明辅助坐标系与基本坐标系之间的关系。
格式:trn O1 O2 O3 B1 B2 B3 B4 B5 B6 B7 B8 B9 M
n:变换号, 1<n<999;
O1 ,O2 ,O3:变换的位移矢量;
B1~B9:变换的旋转矩阵;
M:M=l意味着位移矢量是辅助坐标系的原点在基本坐标系定义的位置。M=-l意味着位移矢量是基本坐标系的原点在辅助坐标系定义的位置。
Trn卡的缺省形式为:
trn 0 0 0 1 0 0 0 1 0 0 0 1 1
一个问题中最多可有10个TRn卡。B矩阵说明两个坐标系轴的方向之间的关系。Bi是在基本坐标系(X Y Z)轴与辅助坐标系(x’y’z’)轴之间的夹角余弦(若在TRn前加了“*”,如*TRn,则由指明余弦值改为指明角度,0o<θ<180o=,其排列规则为:
元素 B1 B2 B3 B4 B5 B6 B7 B8 B9
轴
Bi的意义与M无关。用户不必给出所有Bi值,可用kJ跳过不指定的项,Bi可按以下任一种方式指定:
i)Bi全部给出(9个值);
ii)给出矩阵中列或行方向上三个矢量中的两个(6个值)。程序将由它们的叉乘产生第三个矢量;
iii)给出矩阵中列与行方向上各一个矢量,其中的公共分量必须小于1(共5个值)。程序将用欧拉角方法填完该矩阵;
iv)给出矩阵中列或行方向上的一个矢量(3个值),程序将按某种随机方法生成其他两个矢量;
v)没有输入,程序将建立相同矩阵。
一个矢量由三个元素组成,这些元素可以按矩阵的横向或竖向排列,但不能走对角钱。在所有情况下,MCNP解决任何小的非正交性,并归一这个矩阵。在归一过程,形似(l,0,0)的矢量将保持原样。如果正交性大于0.001弧度时,则给出一个警告信息。
当变换只是一个平移时,使用方式v是合适的。当一些斜置曲面都是关于某个斜轴的旋转面时,使用方式iv)是合适的。曲面不是旋转面(一般情况)时,方式ii)与iii)大体是同等适用的。如果辅助坐标系和基本坐标系中一个是右手坐标系,而另一个是左手坐标系,则使用方式i)是合适的。
使用坐标变换以简化曲面的几何描述,以及建立曲面源问题的坐标系统与写这个曲面源文件中这个问题的坐标系统之间的关系。
通常一个栅元组件将有同一个坐标系统。在这个坐标系统下描述这个组件的曲面,然后用TRn卡指定这个系统与基本坐标系统的关系。
例3-2 17 4 px 5
tr4 7 .9 1.3 0 -1 0 0 0 1 -1 0 0