CN105842834B - 一种由多个非球面反射镜组成的太阳能聚光系统设计方法 - Google Patents

一种由多个非球面反射镜组成的太阳能聚光系统设计方法 Download PDF

Info

Publication number
CN105842834B
CN105842834B CN201610278112.2A CN201610278112A CN105842834B CN 105842834 B CN105842834 B CN 105842834B CN 201610278112 A CN201610278112 A CN 201610278112A CN 105842834 B CN105842834 B CN 105842834B
Authority
CN
China
Prior art keywords
point
coordinate system
vector
aspheric
vertex
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Expired - Fee Related
Application number
CN201610278112.2A
Other languages
English (en)
Other versions
CN105842834A (zh
Inventor
秦华
张爽
李文瑞
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Shandong University of Technology
Original Assignee
Shandong University of Technology
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Shandong University of Technology filed Critical Shandong University of Technology
Priority to CN201610278112.2A priority Critical patent/CN105842834B/zh
Publication of CN105842834A publication Critical patent/CN105842834A/zh
Application granted granted Critical
Publication of CN105842834B publication Critical patent/CN105842834B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G02—OPTICS
    • G02B—OPTICAL ELEMENTS, SYSTEMS OR APPARATUS
    • G02B19/00—Condensers, e.g. light collectors or similar non-imaging optics
    • G02B19/0004—Condensers, e.g. light collectors or similar non-imaging optics characterised by the optical means employed
    • G02B19/0019—Condensers, e.g. light collectors or similar non-imaging optics characterised by the optical means employed having reflective surfaces only (e.g. louvre systems, systems with multiple planar reflectors)
    • G—PHYSICS
    • G02—OPTICS
    • G02B—OPTICAL ELEMENTS, SYSTEMS OR APPARATUS
    • G02B19/00—Condensers, e.g. light collectors or similar non-imaging optics
    • G02B19/0004—Condensers, e.g. light collectors or similar non-imaging optics characterised by the optical means employed
    • G—PHYSICS
    • G02—OPTICS
    • G02B—OPTICAL ELEMENTS, SYSTEMS OR APPARATUS
    • G02B19/00—Condensers, e.g. light collectors or similar non-imaging optics
    • G02B19/0033—Condensers, e.g. light collectors or similar non-imaging optics characterised by the use
    • G02B19/0038—Condensers, e.g. light collectors or similar non-imaging optics characterised by the use for use with ambient light
    • G02B19/0042—Condensers, e.g. light collectors or similar non-imaging optics characterised by the use for use with ambient light for use with direct solar radiation
    • G—PHYSICS
    • G02—OPTICS
    • G02B—OPTICAL ELEMENTS, SYSTEMS OR APPARATUS
    • G02B27/00—Optical systems or apparatus not provided for by any of the groups G02B1/00 - G02B26/00, G02B30/00
    • G02B27/0012—Optical design, e.g. procedures, algorithms, optimisation routines

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Optics & Photonics (AREA)
  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Sustainable Development (AREA)
  • Toxicology (AREA)
  • Lenses (AREA)

Abstract

本发明公开了一种由多个非球面反射镜组成的太阳能聚光系统设计方法。确定组成聚光系统的反射镜数目，把所有反射镜顶点规则的放在YOZ平面内；令入射太阳光束平行于XOY面，其入射单位方向矢量
Figure 735928DEST_PATH_IMAGE001
，取聚焦点P f (x f ,‑x f tanθ,0)在过O点的反射光线上；对于顶点位于O点的反射镜，取不过O点的任一反射光线上一点P′′(x′,y′′,z′′)，P′′与P f 之间的距离R，n条光线的
Figure 325172DEST_PATH_IMAGE002
作为粒子群算法优化程序中的适应度函数；对于顶点在O′(a,b,c)点的非球面，a=0,则以O′点为原点建立新坐标系，新坐标系中
Figure 104909DEST_PATH_IMAGE003
变为
Figure 906643DEST_PATH_IMAGE004
，P f 变为P f ′，用
Figure 517753DEST_PATH_IMAGE004
代替
Figure 340216DEST_PATH_IMAGE003
，P f ′代替P f ，优化得到顶点在O′处的非球面参数
Figure 279353DEST_PATH_IMAGE005
。

Description

一种由多个非球面反射镜组成的太阳能聚光系统设计方法
技术领域
本发明属于太阳能利用领域，特别涉及一种由多个非球面反射镜组成的高聚光比反射型太阳能聚光系统的设计方法。
背景技术
用太阳能反射系统直接把太阳光汇聚在一小的区域内以形成高温热源，一直是太阳能利用领域重要的研究方向，在太阳灶以及利用太阳能点火装置中广泛应用，在太阳能光热发电中也得到了越来越多的实践。顾名思义, 太阳能聚光系统就是把大面积的太阳辐照光收集或者聚焦到一小的接收器或输出装置上的工具。
在212 BC，古希腊物理学家阿基米德就设计了一种太阳能聚光装置，烧毁了罗马舰队。据说这种装置是一种凹面的抛光金属镜面，成百上千的这样的金属镜面都反光在同一条船只上，从而使其起火。
2005年10月，来自麻省理工学院的一群学生验证了这个实验，使用了127块1英尺方镜（30cm），聚焦在位于100英尺（30m）远处的实物模型船上。火光烧毁了船上的一块，但只有在天空晴朗且船静止不动被照射10分钟的情况下才能实现。由此推断，这种装置只有在一定条件下才可行。麻省理工学院研究小组在电视秀“流言终结者”中重复了这个实验，他们在旧金山用一木制渔船作为目标，渔船有被烧焦的痕迹，且还引起了小的火焰。木板的自燃温度大约300 °C。
为了提高聚焦温度，就要增加光学聚光比。槽式太阳能聚光器的理论设计聚光比为285，碟式太阳能聚光器的聚光比已达14000。目前，国内企业生产的大型槽式太阳能聚光器的聚光比一般不会超过80，大型碟式太阳能聚光器的聚光比一般在1000左右。为了提高聚光比，文献Improving the concentration ratio of parabolic troughs using asecond-stage flat mirror采用次级装置和通过优化焦点位置和太阳入射角度来提高聚光比；文献High-concentration solar dishes based on pneumatic reflectingmembranes用充气的弹性塑料薄膜提高蝶式太阳能聚光镜的聚光比；文献Highconcentration thin profile solar concentrator utilizing toroidal confocalrelays用环形共焦中继装置来提高太阳能聚光镜的聚光比；文献Optimized solarthermal concentrator system based on free-form trough reflector和文献Improvedirradiance distribution on high concentration solar cell using free-formconcentrator 采用自由曲面增加提高太阳能聚光器的聚光性能；文献Optical design ofan aspherical cylinder-type reflecting solar concentrator采用非球面面型增加槽式太阳能聚光器的聚光比至285。
非球面面型由多个非球面参数决定，因此设计自由度增多，聚焦能力特别是偏轴聚焦能力增强。但是自由度增多，也意味着设计难度增大。
到目前为止还没有一种完整有效的方法用来设计由多个非球面反射镜组成的太阳能聚光系统，以增加太阳能聚光系统的聚光比。
发明内容
为了提高太阳能聚光系统的聚光比，本发明采用由多个非球面反射镜组成的太阳能聚光系统以增加太阳光的聚光比，提出了这种聚光系统一种完整的设计方法。使用这种方法可以设计出只通过一次反射就能达到或超过聚光比14000的太阳能聚光系统。
本发明解决其技术问题所采用的技术方案是采用如下方法步骤：
1. 在直角坐标系XYZ中，对于8次偶次非球面方程：
Figure 861698DEST_PATH_IMAGE001
式中
Figure 821564DEST_PATH_IMAGE002
，方程式第一项相当于一个二次曲面, C为二次曲面的顶点曲率, 1+a 2为二次曲面系数, a 2与二次曲面的离心率有关; a 4, a 6, a 8为非球面方程多项式中各高次项的系数；
太阳能聚光系统中有m个8次非球面反射镜，2≤m≤200，这些非球面反射镜顶点都在YOZ平面内；即把(0,0,20hmax),(0,0,-20hmax), (0,200,20hmax), (0,200,-20hmax)四点之间的平面划分为Y、Z方向都间隔2hmax cm的网格，非球面反射镜顶点都位于网格的节点上；网格节点的坐标值就是聚光系统中各非球面反射镜顶点(a, b, c)的坐标值。
2.令入射太阳光束平行于XOY面，且与X轴负方向成θ角，在XYZ坐标系中，入射光线单位方向矢量
Figure 738704DEST_PATH_IMAGE003
，α 1= cosθ，β 1=sinθ，γ 1=0；对于顶点在O点的非球面，过O点的反射光线与X轴负方向亦成θ角，取聚焦点在这条光线上，聚焦点坐标为P f (x f , -x f tanθ,0)。
3. 在XYZ坐标系中，对于顶点在O点的非球面反射镜，不过顶点O的反射光线单位方向矢量
Figure 49600DEST_PATH_IMAGE004
，其上任一点坐标P′′(x′′,y′′,z′′)，x′′、z′′与y′′的关系由下式决定
Figure 710388DEST_PATH_IMAGE005
其中，x、y和z为光线与非球面交点P的三个坐标值，令y′′=-x f tanθ，则P′′点与聚焦点P f 之间的距离为
Figure 774028DEST_PATH_IMAGE006
。
4. 取n条光线均匀照射在非球面上，n≥10，则n条光线所有R值的和为
Figure 127649DEST_PATH_IMAGE007
，i表示第i条光线。
5. 取5 cm≤ hmax≤ 20 cm，hmax为h的最大值，以
Figure 660262DEST_PATH_IMAGE008
作为粒子群算法优化程序中的适应度函数，优化非球面参数
Figure 859162DEST_PATH_IMAGE009
，得到使
Figure 528041DEST_PATH_IMAGE008
极小的一组数
Figure 52563DEST_PATH_IMAGE009
，以这组数为参数的非球面就是顶点在O点非球面聚光镜面。
6. 对于顶点不在O点，而在O′(a, b, c)点的非球面反射镜，则以O′点为新坐标系原点，以过O′点的入射光线与O′P f 夹角平分线矢量作为新坐标系的X′轴负方向，矢量
Figure 72471DEST_PATH_IMAGE010
与X′轴负方向单位矢的矢积作为Z′轴正方向单位矢建立新坐标系，由此新坐标系的Y′轴正方向也就确定了。
7. 新坐标系中，过O′点的入射光线与X′轴负方向的夹角为θ 1，由θ 1得到以新坐标系为参考系的入射光线单位矢量
Figure 75063DEST_PATH_IMAGE011
，其中α 1′=cosθ 1，β 1′=sinθ 1，γ 1′=0，新坐标系中聚焦点的坐标P f ′(x f ′, -x f ′tanθ 1, 0)。
8. 用
Figure 864027DEST_PATH_IMAGE011
代替
Figure 310183DEST_PATH_IMAGE003
，用P f ′(x f ′, -x f ′tanθ 1, 0)代替P f (x f , -x f tanθ, 0)，在新坐标系X′Y′Z′中，重复以上步骤3、4、5，优化得到顶点在O′处的非球面反射镜参数
Figure 551808DEST_PATH_IMAGE009
。
所述步骤3包括下述步骤：
（a）已知不过顶点O的入射光线任一点、入射单位矢量
Figure 358090DEST_PATH_IMAGE010
与非球面方程，求得入射光线与非球面交点P (x, y, z)；
（b）由交点P (x, y, z) 与非球面方程得到P处非球面法矢量
Figure 1561DEST_PATH_IMAGE012
；
（c）由
Figure 602307DEST_PATH_IMAGE010
和法线矢量
Figure 596808DEST_PATH_IMAGE012
，用反射定律的矢量公式得到反射光线矢量
Figure 941201DEST_PATH_IMAGE013
：
Figure 704758DEST_PATH_IMAGE014
其中
Figure 476405DEST_PATH_IMAGE015
所述步骤6包括下述步骤：
（a） 在坐标系XYZ中，矢量
Figure 692623DEST_PATH_IMAGE016
与-
Figure 89975DEST_PATH_IMAGE010
夹角2θ 1为：
Figure 708038DEST_PATH_IMAGE017
其夹角的平分线矢量为：
Figure 916165DEST_PATH_IMAGE018
Figure 619679DEST_PATH_IMAGE019
令:
Figure 305876DEST_PATH_IMAGE020
Figure 778445DEST_PATH_IMAGE021
Figure 891895DEST_PATH_IMAGE022
则此角平分线矢量的方向余弦为：
Figure 82705DEST_PATH_IMAGE023
以此夹角平分线矢量作为新坐标系的X′轴负方向，则X′轴正方向的方向余弦为：
Figure 572592DEST_PATH_IMAGE024
（b）矢量
Figure 915980DEST_PATH_IMAGE025
与X′轴构成的平面法矢量就是
Figure 465910DEST_PATH_IMAGE026
：
Figure 144016DEST_PATH_IMAGE027
Figure 906435DEST_PATH_IMAGE028
Figure 353597DEST_PATH_IMAGE026
的方向余弦：
Figure 74428DEST_PATH_IMAGE029
Figure 239831DEST_PATH_IMAGE030
Figure 71520DEST_PATH_IMAGE031
新坐标系的Y′轴的方向余弦为
Figure 373189DEST_PATH_IMAGE032
Figure 999342DEST_PATH_IMAGE033
令：
Figure 901308DEST_PATH_IMAGE034
所述步骤7包含下述步骤：
（a）
Figure 271109DEST_PATH_IMAGE035
（b）坐标系XYZ中的坐标（x, y, z）与坐标系X′Y′Z′中的坐标（x′, y′, z′）满足下列关系:
x= x′cosα 1′+ y′cosα 2′+ z′cosα 3′+a
y= x′cosβ 1′+ y′cosβ 2′+ z′cosβ 3′+b
z= x′cosγ 1′+ y′cosγ 2′+ z′cosγ 3′+c
由这三个表达式的逆变换就可以得到XYZ坐标系中聚焦点坐标P f （x f , -x f tanθ,0）到X′Y′Z′坐标系中聚焦点坐标P f ′的变换。
本发明的有益效果
本发明采用非球面反射镜及多个非球面反射镜组合的方式增加太阳能点聚焦的聚光比，利用粒子群算法设计出只通过一次反射就能实现高聚光比的太阳能聚光系统，以方便实现太阳能的高温热利用。
用非球面代替二次曲面使反射镜面的设计有更多个自由度，使聚光镜镜面面型有更多选择，与旋转抛物面相比，利用非球面反射镜可以更好的实现偏轴聚焦。在确定偏轴聚焦点的情况下，利用本发明方法可以方便地确定非球面反射镜方程的参数C, a 2, a 4, a 6, a 8。
对于使用本发明方法的设计结果可以根据具体的非球面方程通过数控机床加工成模具，实现流水线生产，加工工艺简单。此外，本发明中的方法是通过计算机编程实现的，完全可以自动寻找到最佳组合的非球面参数C, a 2, a 4, a 6, a 8, 具有快时、方便等优点，在太阳能应用领域特别是在非球面太阳能聚光镜设计领域具有一定的应用前景。
下面结合附图和实施例对本发明进一步说明
附图说明：
图1：太阳能聚光系统中各非球面反射镜顶点在YOZ平面内的位置分布示意图；
图2：实施例1中由21个非球面反射镜组成的太阳能聚光系统，各非球面反射镜顶点位置示意图；
图3：非球面反射镜顶点在O处时，其过O点的反射光线与不过O点的反射光线的相对位置关系示意图；
图4：非球面反射镜顶点不在O点而是在O′点时，以O′点为原点的新坐标系X′Y′Z′建立示意图；
图5：实施例1中，一非球面反射镜顶点在O′（0,200, -20）点时， 对太阳光的聚光效果图；
图6：单个非球面反射镜聚光效果示意图；
图7：实施例1中所设计的由21个非球面反射镜组成的太阳能聚光系统对太阳光束反射聚光模拟图，其中图7（a）和（b）分别表示从不同角度观察其聚光效果的模拟图示；
图8：6个非球面反射镜的聚光系统对太阳光聚光效果模拟图。
图中：1.非球面反射镜 2.入射光束 3.反射光束 4.聚焦光斑
具体实施方式：
实施例1：
所述一种由多个非球面组成的太阳能聚光系统设计方法包含步骤如下：
1. 在直角坐标系XYZ中，对于8次偶次非球面方程，
Figure 692863DEST_PATH_IMAGE001
式中，
Figure 755497DEST_PATH_IMAGE002
，方程式第一项相当于一个二次曲面, 公式中C为二次曲面的顶点曲率, 1+a 2为二次曲面系数, a 2与二次曲面的离心率有关; a 4, a 6, a 8为非球面方程多项式中各高次项的系数；
请参阅图2,太阳能聚光系统中有m个8次非球面反射镜，2≤m≤200，在此取m=21，这些非球面顶点都位于图2网格中的粗线上，即21个非球面反射镜顶点坐标为（0, 200,200）、（0, 200,180）、（0, 200,160）、（0, 200,140）、（0, 200,120）、（0, 200,100）、（0,200,80）、（0, 200,60）、（0, 200,40）、（0, 200,20）、（0, 200,0）、（0, 200,-20）、（0, 200,-40）、（0, 200,-60）、（0, 200,-80）、（0, 200,-100）、（0, 200,-120）、（0, 200,-140）、（0,200,-160）、（0, 200,-180）、（0, 200,-200）。
2. 请参阅图3,令入射太阳光束平行于XOY面，且与X轴负方向成θ角，θ=250，在XYZ坐标系中，入射光线单位方向矢量
Figure 364333DEST_PATH_IMAGE003
，α 1= cos250，β 1=sin250，γ 1=0；对于顶点在O点的非球面，过O点的反射光线与X轴负方向亦成θ=250角，取聚焦点在这条光线上，聚焦点坐标为P f (x f , -x f tanθ, 0)，取x f =-400，则P f (-400, 186.5231, 0) 。
3. 请参阅图3, 在XYZ坐标系中，对于顶点在O点的非球面，不过顶点O的反射光线单位方向矢量
Figure 803405DEST_PATH_IMAGE004
，其上任一点坐标P′′(x′′,y′′,z′′)，x′′、z′′与y′′的关系由下式决定
Figure 79665DEST_PATH_IMAGE005
其中，x、y和z为光线与非球面交点P的三个坐标值，令y′′=-x f tanθ=-400×tan250=186.5231，则P′′点与聚焦点P f 之间的距离为
Figure 313201DEST_PATH_IMAGE036
4.取n条光线均匀照射在非球面上，n≥10，令n=20，则20条光线所有R值的和为
Figure 674912DEST_PATH_IMAGE037
。
5．取hmax=10 cm，以
Figure 137248DEST_PATH_IMAGE008
作为粒子群算法优化程序中的适应度函数，优化非球面参数
Figure 533595DEST_PATH_IMAGE009
，得到使
Figure 938031DEST_PATH_IMAGE008
极小的一组非球面参数C=-1/910.0, a 2=91.0579, a 4=-1.7817×10-8, a 6=1.7726×10-11, a 8=-3.4078×10-14,以这组数为参数的非球面就是顶点在O点非球面反射镜面；
6. 请参阅图4和图5,对于顶点不在O点，而在O′(a, b, c)点的非球面，如O′(0,200, -20), 则以O′点为新坐标系原点，以过O′点的入射光线与O′P f 夹角平分线矢量作为新坐标系的X′轴负方向，矢量
Figure 521459DEST_PATH_IMAGE010
与X′轴负方向单位矢的矢积就是Z′轴正方向单位矢，因此新坐标系的Y′轴正方向也就确定了；
7. 新坐标系中，过O′(0, 200, -20)点的入射光线与X′轴负方向的夹角为θ 1=11.6187o，由θ 1得到以新坐标系为参考系的入射光线单位矢量
Figure 36754DEST_PATH_IMAGE011
，其中α 1′=cos11.6187o，β 1′=sin11.6187o，γ 1′=0，新坐标系中聚焦点的坐标P f ′(x f ′, -x f ′tanθ 1,0)，x f ′=-392.5153，-x f ′tanθ 1=80.7053；
8.用
Figure 287607DEST_PATH_IMAGE011
代替
Figure 862945DEST_PATH_IMAGE003
，用P f ′(x f ′, -x f ′tanθ 1, 0)代替P f (x f , -x f tanθ, 0)，在新坐标系X′Y′Z′中，重复步骤3、4、5，优化得到顶点在O′(0, 200, -20)处的非球面反射镜参数C=-1/802.2545, a 2= 81.2159, a 4=1.1144×10-8, a 6=-1.4847×10-11, a 8=3.7373×10-14, 用此数据得到的非球面反射镜聚光效果如图5所示，图6示出了非球面反射镜较为清晰的聚光效果示意图。
所述步骤3包括下述步骤：
（a）已知不过顶点O的入射光线任一点、入射单位矢量
Figure 199248DEST_PATH_IMAGE010
与非球面方程，求得入射光线与非球面交点P (x, y, z)；
（b）由交点P (x, y, z) 与非球面方程得到P处非球面法矢量
Figure 518234DEST_PATH_IMAGE012
；
（c）由
Figure 358014DEST_PATH_IMAGE010
和法线矢量
Figure 365240DEST_PATH_IMAGE012
，用反射定律的矢量公式得到反射光线矢量
Figure 923260DEST_PATH_IMAGE013
：
Figure 45937DEST_PATH_IMAGE014
其中
Figure 5802DEST_PATH_IMAGE015
。
所述步骤6包括下述步骤：
（a） 在坐标系XYZ中，矢量
Figure 922943DEST_PATH_IMAGE016
与-
Figure 233838DEST_PATH_IMAGE010
夹角2θ 1为：
Figure 894627DEST_PATH_IMAGE038
其夹角的平分线矢量为：
Figure 443420DEST_PATH_IMAGE018
Figure 797041DEST_PATH_IMAGE019
令:
Figure 80386DEST_PATH_IMAGE039
Figure 279286DEST_PATH_IMAGE040
Figure 213744DEST_PATH_IMAGE041
则此角平分线矢量的方向余弦为：
Figure 472687DEST_PATH_IMAGE042
Figure 227016DEST_PATH_IMAGE043
Figure 495187DEST_PATH_IMAGE044
以此夹角平分线矢量作为新坐标系的X′轴负方向，则X′轴正方向的方向余弦为：
Figure 284151DEST_PATH_IMAGE045
（b）矢量
Figure 713995DEST_PATH_IMAGE025
与X′轴构成的平面法矢量就是
Figure 221200DEST_PATH_IMAGE026
：
Figure 276750DEST_PATH_IMAGE027
Figure 920221DEST_PATH_IMAGE028
Figure 520966DEST_PATH_IMAGE026
的方向余弦：
Figure 249888DEST_PATH_IMAGE046
Figure 859861DEST_PATH_IMAGE047
Figure 623417DEST_PATH_IMAGE048
新坐标系的Y′轴的方向余弦为
Figure 395064DEST_PATH_IMAGE049
Figure 876861DEST_PATH_IMAGE050
令：
Figure 24946DEST_PATH_IMAGE051
Figure 377430DEST_PATH_IMAGE052
Figure 70710DEST_PATH_IMAGE053
所述步骤7包下述步骤：
（a）
Figure 774224DEST_PATH_IMAGE054
（b）坐标系XYZ中的坐标（x, y, z）与坐标系X′Y′Z′中的坐标（x′, y′, z′）满足下列关系:
x= x′cosα 1′+ y′cosα 2′+ z′cosα 3′+a
y= x′cosβ 1′+ y′cosβ 2′+ z′cosβ 3′+b
z= x′cosγ 1′+ y′cosγ 2′+ z′cosγ 3′+c
由这三个表达式的逆变换就可以得到XYZ坐标系中聚焦点坐标P f （x f , -x f tanθ,0）到X′Y′Z′坐标系中聚焦点坐标P f ′的变换。
如：O′(0, 200, -20), 即a=0，b=200，c=-20，x= x f =-400，y=-x f tanθ=186.5231, z=0, 由上面三个表达式的逆变换得到x f ′=x′=-392.5153,-x f ′tanθ 1= y′=80.7053, z′=0，在新坐标系X′Y′Z′中聚焦点坐标P f ′(-392.5153, 80.7053, 0)。
对于非球面顶点分别在XYZ坐标系中的坐标点（0, 200,200）、（0, 200,180）、（0,200,160）、（0, 200,140）、（0, 200,120）、（0, 200,100）、（0, 200,80）、（0, 200,60）、（0, 200,40）、（0, 200,20）、（0, 200,0）、（0, 200,-40）、（0, 200,-60）、（0, 200,-80）、（0,200,-100）、（0, 200,-120）、（0, 200,-140）、（0, 200,-160）、（0, 200,-180）、（0, 200,-200）时，重复步骤6、7、8，得到表1的非球面参数
Figure 726000DEST_PATH_IMAGE009
数据，表1中的最后一列给出了此非球面反射镜在聚光点处的最大光斑半径max（R）。图7是用表1中的数据模拟实施例1所设计的太阳能聚光系统的聚焦效果，图7（a）和（b）是从不同角度观察实施例1所设计的聚光系统对太阳光束的聚焦效果示意图，为了更清楚的观察聚光系统中各非球面反射镜、聚焦点及光束之间的相对位置，图8示出了有6个非球面反射镜的聚光系统。
Figure DEST_PATH_IMAGE056AAAAAA

Claims (5)

1.一种由多个非球面反射镜组成的太阳能聚光系统设计方法，在直角坐标系XYZ中，对于8次偶次非球面方程，
Figure FDA0002598487700000011
式中
Figure FDA0002598487700000012
方程式第一项相当于一个二次曲面，C为二次曲面的顶点曲率，1+a2为二次曲面系数，a2与二次曲面的离心率有关；a4，a6，a8为非球面方程多项式中各高次项的系数，其特征在于包含下述步骤：
(1)太阳能聚光系统中有m个8次非球面反射镜，2≤m≤200，这些非球面顶点都在YOZ平面内；
(2)入射太阳光束平行于XOY面，且与X轴负方向成θ角，在XYZ坐标系中，入射光线单位方向矢量
Figure FDA0002598487700000013
α1＝cosθ，β1＝sinθ，γ1＝0；对于顶点在O点的非球面反射镜，过O点的反射光线与X轴负方向亦成θ角，取聚焦点在这条光线上，聚焦点坐标为Pf(xf，-xftanθ，0)；
(3)在XYZ坐标系中，对于顶点在O点的非球面，不过顶点O的反射光线单位方向矢量
Figure FDA0002598487700000014
其上任一点坐标P″(x″，y″，z″)，x″、z″与y″的关系由下式决定
Figure FDA0002598487700000015
其中，x、y和z为光线与非球面交点P的三个坐标值，令y″＝-xftanθ，则P″点与聚焦点Pf之间的距离为
Figure FDA0002598487700000016
(4)取n条光线均匀照射在此非球面反射镜上，n≥10，则n条光线所有R值的和为
Figure FDA0002598487700000017
i表示第i条光线；
(5)取5cm≤hmax≤20cm，hmax为h的最大值，以
Figure FDA0002598487700000018
作为粒子群算法优化程序中的适应度函数，优化非球面参数C,a2,a4,a6,a8，得到使
Figure FDA0002598487700000019
极小的一组数C,a2,a4,a6,a8，以这组数为参数的非球面就是顶点在O点非球面反射聚光镜面；
(6)对于顶点不在O点，而在O′(a，b，c)点的非球面，则以O′点为新坐标系原点建立新坐标系，以过O′点的入射光线与O′Pf夹角平分线矢量作为新坐标系的X′轴负方向，矢量
Figure FDA00025984877000000210
与X′轴负方向单位矢的矢积作为Z′轴正方向单位矢建立新坐标系，由此新坐标系的Y′轴正方向也就确定了；
(7)新坐标系X'Y'Z'中，过O′点的入射光线与X′轴负方向的夹角为θ1，由θ1得到以新坐标系为参考系的入射光线单位矢量
Figure FDA0002598487700000021
其中α1′＝cosθ1，β1′＝sinθ1，γ1′＝0，新坐标系中聚焦点的坐标Pf′(xf′，-xf′tanθ1，0)；
(8)用
Figure FDA0002598487700000022
代替
Figure FDA0002598487700000023
用Pf′(xf′，-xf′tanθ1，0)代替Pf(xf，-xftanθ，0)，在新坐标系X′Y′Z′中，重复以上步骤(3)、(4)、(5)，优化得到顶点在O′处的非球面反射镜参数C,a2,a4,a6,a8。
2.根据权利要求1的一种由多个非球面反射镜组成的太阳能聚光系统设计方法，其特征在于所述步骤(1)包括下述步骤：
(a)把(0，0，20hmax)，(0，0，-20hmax)，(0，200，20hmax)，(0，200，-20hmax)四点之间的平面划分为Y、Z方向都间隔2hmax的网格，非球面反射镜顶点都位于网格的节点上；
(b)网格节点的坐标值就是聚光系统中各非球面顶点O′(a，b，c)的坐标值。
3.根据权利要求1的一种由多个非球面反射镜组成的太阳能聚光系统设计方法，其特征在于所述步骤(3)包括下述步骤：
(a)已知不过顶点O的入射光线任一点、入射光线单位方向矢量
Figure FDA0002598487700000024
与非球面方程，求得入射光线与非球面交点P(x，y，z)；
(b)由交点P(x，y，z)与非球面方程得到P处非球面法矢量
Figure FDA0002598487700000025
(c)由
Figure FDA0002598487700000026
和法线矢量
Figure FDA0002598487700000027
用反射定律的矢量公式得到反射光线单位方向矢量
Figure FDA0002598487700000028
Figure FDA0002598487700000029
其中g＝-2|α1αN+β1βN+γ1γN|。
4.根据权利要求1的一种由多个非球面反射镜组成的太阳能聚光系统设计方法，其特征在于所述步骤(6)包括下述步骤：
(a)在坐标系XYZ中，矢量
Figure FDA0002598487700000031
与
Figure FDA0002598487700000032
夹角2θ1为：
Figure FDA0002598487700000033
其夹角的平分线矢量为：
Figure FDA0002598487700000034
令：
Figure FDA0002598487700000035
Figure FDA0002598487700000036
则此角平分线矢量的方向余弦为：
Figure FDA0002598487700000037
以此夹角平分线矢量作为新坐标系的X′轴负方向，则X′轴正方向的方向余弦为：
cosα′1＝-Dirx，cosβ′1＝-Diry，cosγ′1＝-Dirz
(b)矢量
Figure FDA0002598487700000038
与X′轴构成的平面法矢量就是
Figure FDA0002598487700000039
Figure FDA00025984877000000310
Figure FDA00025984877000000311
的方向余弦：
Figure FDA00025984877000000312
Figure FDA00025984877000000313
Figure FDA00025984877000000314
新坐标系的Y′轴的方向余弦为
Figure FDA0002598487700000041
令：
Figure FDA0002598487700000042
5.根据权利要求1的一种由多个非球面反射镜组成的太阳能聚光系统设计方法，其特征在于所述步骤(7)包下述步骤：
(a)
Figure FDA0002598487700000043
(b)坐标系XYZ中的坐标(x，y，z)与坐标系X′Y′Z′中的坐标(x′，y′，z′)满足下列关系：
x＝x′cosα1′+y′cosα2′+z′cosα3′+a
y＝x′cosβ1′+y′cosβ2′+z′cosβ3′+b
z＝x′cosγ1′+y′cosγ2′+z′cosγ3′+c
由这三个表达式的逆变换就可以得到XYZ坐标系中聚焦点坐标Pf(xf，-xftanθ，0)到X′Y′Z′坐标系中聚焦点坐标Pf′的变换。
CN201610278112.2A 2016-05-03 2016-05-03 一种由多个非球面反射镜组成的太阳能聚光系统设计方法 Expired - Fee Related CN105842834B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610278112.2A CN105842834B (zh) 2016-05-03 2016-05-03 一种由多个非球面反射镜组成的太阳能聚光系统设计方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610278112.2A CN105842834B (zh) 2016-05-03 2016-05-03 一种由多个非球面反射镜组成的太阳能聚光系统设计方法

Publications (2)

Publication Number Publication Date
CN105842834A CN105842834A (zh) 2016-08-10
CN105842834B true CN105842834B (zh) 2020-11-06

Family

ID=56590003

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610278112.2A Expired - Fee Related CN105842834B (zh) 2016-05-03 2016-05-03 一种由多个非球面反射镜组成的太阳能聚光系统设计方法

Country Status (1)

Country Link
CN (1) CN105842834B (zh)

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102866491A (zh) * 2012-10-17 2013-01-09 山东理工大学 高次柱面反射型太阳能聚光镜
CN103487925A (zh) * 2013-09-10 2014-01-01 山东理工大学 一种组合非球面反射型太阳能聚光镜镜面的设计方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20080049342A1 (en) * 2006-08-24 2008-02-28 David Timothy Dobney Balanced Radiation Concentrating Device

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102866491A (zh) * 2012-10-17 2013-01-09 山东理工大学 高次柱面反射型太阳能聚光镜
CN103487925A (zh) * 2013-09-10 2014-01-01 山东理工大学 一种组合非球面反射型太阳能聚光镜镜面的设计方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
《组合非球面太阳能聚光镜的光学设计》;秦华，冯东太，葛硕硕，王勇;《中国光学》;20141031;第7卷(第5期);844-854 *
《高次柱面反射型太阳能聚光镜的光学设计》;秦华 类成新 刘汉法 葛硕硕;《物理学报》;20130531;第62卷;104215-1-104215-8 *

Also Published As

Publication number Publication date
CN105842834A (zh) 2016-08-10

Similar Documents

Publication Publication Date Title
Dähler et al. Optical design and experimental characterization of a solar concentrating dish system for fuel production via thermochemical redox cycles
US6620995B2 (en) Non-imaging system for radiant energy flux transformation
Coughenour et al. Dish-based high concentration PV system with Köhler optics
Meng et al. A novel free-form Cassegrain concentrator for PV/T combining utilization
Cooper et al. Theory and design of line-to-point focus solar concentrators with tracking secondary optics
CN103487925B (zh) 一种组合非球面反射型太阳能聚光镜镜面的设计方法
Tsai Improved irradiance distribution on high concentration solar cell using free-form concentrator
Tsai Optimized solar thermal concentrator system based on free-form trough reflector
Shanks et al. High-concentration optics for photovoltaic applications
Hernández et al. High-performance Köhler concentrators with uniform irradiance on solar cell
Wardhana et al. Designing and modeling a novel dual parabolic concentrator with three degree of freedom (DOF) robotic arm
Lee et al. Calculation of optical efficiency for the first central-receiver solar concentrator system in Korea
Tavakol-Moghaddam et al. Optimal design of solar concentrator in multi-energy hybrid systems based on minimum exergy destruction
CN105842834B (zh) 一种由多个非球面反射镜组成的太阳能聚光系统设计方法
US9244262B2 (en) Low cost focussing system giving high concentrations
US20150207455A1 (en) Dense-Array Concentrator Photovoltaic System Utilising Non-Imaging Dish Concentrator And Array Of Crossed Compound Parabolic Concentrators
Chong et al. Solar flux distribution analysis of Non-Imaging Planar Concentrator for the application in concentrator photovoltaic system
Chiromawa et al. Concept of Bee-eyes array of Fresnel lenses as a solar photovoltaic concentrator system
Berryman et al. Novel solar concentrator geometry for point focus, low-cost applications
Ferrer-Rodríguez et al. Optimization of an ultra-high CPV Cassegrain-Koehler unit with 2000× concentration ratio
Pavlović et al. Optical design of a solar parabolic thermal concentrator based on trapezoidal reflective petals
Georgiou et al. Caustics as an alternate of ray tracing to evaluate heliostat mirrors
Guo et al. Two new methods used to simulate the circumferential solar flux density concentrated on the absorber of a parabolic trough solar collector
El Himer et al. Improvement of optical performances using the hybrid CPV
CN105759412B (zh) 一种由190个非球面反射镜组成的太阳能聚光系统

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20201106

Termination date: 20210503

CF01 Termination of patent right due to non-payment of annual fee