本发明属于集热器技术领域，尤其是非抛物线反射镜的蒙特卡洛模拟方法。
背景技术
在光场设备的施工过程中，需要借助软件算法分析影响槽式集热器的光学和能量转化效率的各种因素，用以评估最优的光场设备和施工建造成本，为不同形式集热器或集热回路的可行性提供依据。软件算法的依据为蒙特卡洛法，即当所求解问题是某种随机事件出现的概率，或者是某个随机变量的期望值时，通过某种“实验”的方法，以这种事件出现的频率估计这一随机事件的概率，或者得到这个随机变量的某些数字特征，并将其作为问题的解。本分析方法涉及的计算原理为蒙特卡洛光线追踪法，计算光线自入射到反射镜、反射镜反射后射到集热管的过程，认为轴向无变化。计算涉及到四个随机量：
第一个随机量是入射时太阳光光线锥度角，n条光线将有n个随机的锥度角，在太阳光锥角-0.26°到0.26°范围内平均分布。
第二个随机量是光线入射在反射镜上的位置，认为入射点在水平方向即x方向平均分布。
第三个随机量是入射点位置反射镜的制造误差，该误差服从均值为零的正态分布。
第四个随机量是跟踪系统跟踪误差，服从均值为零的平均分布。
当入射光线数目足够多时，最后的光线反射将接近真实情况。
在现有技术中，存在以下几个：缺少对于非标准抛物线方程的求解，目前无任何论文提出变形后非抛物线反射镜的计算。
技术实现要素：
本发明克服了现有技术中的缺点，提供了非抛物线反射镜的蒙特卡洛模拟方法。
为了解决上述技术问题，本发明是通过以下技术方案实现的：
非抛物线反射镜的蒙特卡洛模拟方法，依据蒙特卡洛光线追踪法，计算涉及以下四个随机量的模拟：入射时太阳光光线锥度角、光线入射在反射镜上的位置、入射点位置反射镜的制造误差，和跟踪系统跟踪误差。
步骤一、根据输入的时刻、经纬度、地基误差计算太阳位置，确定入射角；同时计算出排遮阳系数、iam、ani、末端损失和跟踪角等；
步骤二、确定入射的n条光线的阳光锥角；
步骤三、确定入射光在反射镜上位置；
步骤四、确定入射点反射镜制造误差；
步骤五、联立入射光方程y＝kx+b和反射镜方程y＝(x+tp)^2/4p+c。求出入射点位置(x0,y0)。求导反射镜方程得到该位置的切线的斜率k0，k0+m为考虑反射镜制造误差后的斜率，1/(k0+m)为反射镜在该点的法线斜率。入射光与法线夹角等于反射光与法线夹角，以此计算出反射光斜率。由点线法，可求出反射光的几何方程y＝k’x+b’；
步骤六、联立反射光方程y＝k’x+b’和集热管方程(x+t1)^2/a^2+(y+t2)^2/b^2＝1，求出的解为反射光入射到集热管上的位置点坐标；
步骤七、统计光线携带能量，计算出集热管上热流分布。将反射光打到集热管上的坐标(x,y)转换为极坐标(p，r)，极坐标原点为集热管中心，组成新矩阵(p,r,i’)。以角度p对矩阵(p,r,i’)排序，每1°区间内i相加后得到热量ia，ia除以1°对应集热管面积得到热流密度p。集热效率即为ia之和与入射能量i之和。
在步骤五中，当用非抛物线计算方式时，输入反射镜一组坐标(x,y)值，用matlab拟合功能得到拟合函数，求导得到切线斜率k0，之后步骤与抛物线计算相同。
在步骤六中，反射光方程与集热管方程联立得到的n个解中有部分解是没有意义的，为虚数解，这些虚数解表示光线因为斜率过大或过小会打不到集热管上，即为漏光光线，排除虚数解就能排除漏光光线；同时排除i矩阵中相同位置的值得到i’，同时求出反射光入射到集热管上时的轴向和径向夹角；将排除掉的光线输出到漏光统计图中，将反射镜间隙间未能反射的光线输出到未反射光统计图中。
与现有技术相比，本发明的有益效果是：
本方法适用于任何曲线形式的反射镜，即能按照标准抛物线方程求解，也能计算变形后的反射镜，只需要输入非抛物线反射镜多个点的坐标即可。
附图说明
图1为本方法流程图。
图2为图集热器与太阳位置关系示意图。
图3阳光锥角与能量关系示意图。
图4为集热器形状的示意图。
图5为入射光锥角统计分布图。
图6为反射光倾角的统计分布图。
图7为反射镜漏光在x轴的统计分布图。
图8为径向偏角统计分布图。
图9为总入射偏角统计分布图。
图10为反射光在集热管上的统计分布图。
图11为集热管表面热流密度分布示意图。
图12为集热管上热流密度分布图。
具体实施方式
下面结合附图与具体的实施方式对本发明作进一步详细描述：
如图中所示，非抛物线反射镜的蒙特卡洛模拟方法，按照下列流程进行：
步骤一、根据输入的时刻、经纬度、地基误差计算太阳位置，确定入射角；同时计算出排遮阳系数、iam、ani、末端损失和跟踪角等。计算公式如下：单位为默认国际单位制。
e＝229.18(0.000075+0.001868cosb-0.032077sinb-0.014615cos2b-0.04089sin2b)，e为表达地球椭圆轨道导致全年白天时长不一致的函数。
步骤二、确定入射的n条光线的阳光锥角，平均分布。用matlab的rand命令一次性生成n个-0.26°至0.26°角度(0.26＝0.53/2)，存储在一个矩阵k中,同时按照锥角与能量关系计算出每个锥角的光线携带的能量矩阵i，用rand命令生成n个-e至e之间的随机数，e为跟踪误差，存为矩阵ez。使k＝k+ez+e’,e’为入射偏角，即某个原因造成的集热器与入射光固定偏角。则k为考虑跟踪误差后的入射光角度。输出k的统计分布图(如图2所示)。
步骤三、确定入射光在反射镜上位置，x方向平均分布。用matlab的rand命令一次性生成n个x1-x2之间的x值，存于矩阵x中。反射镜坐标系原点为抛物线顶点，集热管坐标系原点为集热管中心，相当于反射镜坐标系上移一个焦距p。
步骤四、确定入射点反射镜制造误差，正态分布。用matlab的randn命令一次性生成n个均值0标准差为m的数存入矩阵m中。在矩阵k，x，m中，相同位置的数即表示的是同一条入射光的入射角、入射到反射镜上时的x坐标，以及该位置反射镜的光学误差。
步骤五、联立方程：入射光方程y＝kx+b和反射镜方程y＝(x+tp)^2/4p+c。可求出入射点位置(x0,y0)。求导反射镜方程得到该位置的切线的斜率k0，k0+m为考虑反射镜制造误差后的斜率，1/(k0+m)为反射镜在该点的法线斜率。入射光与法线夹角等于反射光与法线夹角，以此计算出反射光斜率。由点线法，可求出反射光的几何方程y＝k’x+b’。其中tp,4p,c为输入值，b由入射光斜率k和经过的抛物线上的点(x0,y0)用点线法求得。
当用非抛物线计算方式时，输入反射镜一组坐标(x,y)值，用matlab拟合功能得到拟合函数，求导得到切线斜率k0，之后步骤与抛物线计算相同。
步骤六、联立方程：反射光方程y＝k’x+b’和集热管方程(x+t1)^2/a^2+(y+t2)^2/b^2＝1，求出的解为反射光入射到集热管上的位置点坐标，解的数目为n个(x,y)值。因为某些光线因为斜率过大或过小会打不到集热管上，所以n个解中有部分解是没有意义的，即虚数解，排除虚数解就能排除漏光光线，同时排除i矩阵中相同位置的值得到i’。同时求出反射光入射到集热管上时的轴向和径向夹角。排除掉的光线输出到漏光统计图中，反射镜间隙间未能反射的光线输出到未反射光统计图中。
步骤七、统计光线携带能量，计算出集热管上热流分布。将反射光打到集热管上的坐标(x,y)转换为极坐标(p，r)，极坐标原点为集热管中心，组成新矩阵(p,r,i’)。以角度p对矩阵(p,r,i’)排序，每1°区间内i相加后得到热量ia，ia除以1°对应集热管面积得到热流密度p。集热效率即为ia之和与入射能量i之和。
本方法中的输入参数与输出参数如下表所示：
表1.输入参数与输出参数
以上对本发明进行了详细说明，但所述内容仅为本发明的较佳实施例，不能被认为用于限定本发明的实施范围。凡依本发明申请范围所作的均等变化与改进等，均应仍归属于本发明的专利涵盖范围之内。