设计准备
所需工具
MATLAB 2021b ; solidworks 2021 ; TracePro;
设计要求
设计一个焦距 f = 200mm, n = 1.5, r = 50mm 的菲涅尔透镜。其中 n
是菲涅尔透镜的 折射率，r 是透镜的球冠半径。
设计原理
边光原理
原理：若一个具有一定光展的光源发出的光线经过一个或几个序列化的单调正则光学表面的光学作用后投向目标面形成光斑，则整个光学过程具有以下两个性质：一是光源的边光成为光斑的边光；二是光源的内光映射为光斑的内光，并保持光线的拓扑结构不变。
根据边光原理，在设计集光器的外形尺寸时，只需要考虑边光，进而可以简化设计的难度。
近轴光线在球面光学系统的光路计算
光线在光轴附近很小的区域内，这个区域称为近轴区，近轴区内的光线称为近轴光线。根据基本的几何光学原理，光线在界面的传播可以用如下的公式表示。
i = l − r r u i ′ = n n ′ i u ′ = u + i − i ′ l ′ = r ( 1 + i ′ u ′ ) \begin{array}{c} i=\frac{l-r}{r} u \\ i^{\prime}=\frac{n}{n^{\prime}} i \\ u^{\prime}=u+i-i^{\prime} \\ l^{\prime}=r\left(1+\frac{i^{\prime}}{u^{\prime}}\right) \end{array} i=rl−rui′=n′niu′=u+i−i′l′=r(1+u′i′) 整理该公式，可以得到光线在物空间与像空间的成像关系。
n ′ ( 1 r − 1 l ′ ) = n ( 1 r − 1 l ) = Q n ′ u ′ − n u = ( n ′ − n ) h r n ′ l ′ − n l = n ′ − n r \begin{array}{c} n^{\prime}\left(\frac{1}{r}-\frac{1}{l^{\prime}}\right)=n\left(\frac{1}{r}-\frac{1}{l}\right)=Q \\ n^{\prime} u^{\prime}-n u=\left(n^{\prime}-n\right) \frac{h}{r} \\ \frac{n^{\prime}}{l^{\prime}}-\frac{n}{l}=\frac{n^{\prime}-n}{r} \end{array} n′(r1−l′1)=n(r1−l1)=Qn′u′−nu=(n′−n)rhl′n′−ln=rn′−n
光在透镜中的传播
根据焦点的物会成像在无限远的性质，可以得到焦距与球面曲率半径的关系()。
f
1
=
−
r
1
n
−
1
,
f
1
′
=
n
r
1
n
−
1
f
2
=
n
r
2
n
−
1
,
f
2
′
=
−
r
2
n
−
1
\begin{array}{c} f_{1}=-\frac{r_{1}}{n-1}, \quad f_{1}^{\prime}=\frac{n r_{1}}{n-1} \\ f_{2}=\frac{n r_{2}}{n-1}, \quad f_{2}^{\prime}=-\frac{r_{2}}{n-1} \\ \end{array}
f1=−n−1r1,f1′=n−1nr1f2=n−1nr2,f2′=−n−1r2
Δ
=
d
−
f
1
′
+
f
2
\Delta=d-f_{1}^{\prime}+f_{2}
Δ=d−f1′+f2
又由于透镜由两个球面组成，故根据组合光学系统的转面公式[转面公式]{reference-type=“ref”
reference=“转面公式”}，可以得到两组合光学系统的等效焦距、等下焦点。当透镜两边都为空气时，可以得到如下的公式
f
′
=
−
f
=
−
f
1
′
f
2
′
Δ
=
n
r
1
r
2
(
n
−
1
)
[
n
(
r
2
−
r
1
)
+
(
n
−
1
)
d
]
f^{\prime}=-f=-\frac{f_{1}^{\prime} f_{2}^{\prime}}{\Delta}=\frac{n r_{1} r_{2}}{(n-1)\left[n\left(r_{2}-r_{1}\right)+(n-1) d\right]}
f′=−f=−Δf1′f2′=(n−1)[n(r2−r1)+(n−1)d]nr1r2
当透镜的厚度d很小，可以忽略时，透镜的焦距的倒数（即光焦距）可以如下表示
Φ
=
1
f
′
=
(
n
−
1
)
(
ρ
1
=
ρ
2
)
\Phi=\frac{1}{f^{'}}=(n-1)(\rho_1=\rho_2)
Φ=f′1=(n−1)(ρ1=ρ2)
ρ
1
,
ρ
2
\rho_1,\rho_2
ρ1,ρ2
为透镜两个折射面的曲率半径的倒数。一般的透镜的折射率为1.5，故可以近似的认为焦距与半径的关系为
f
=
2
R
f=2R
f=2R
菲涅尔透镜原理
透镜在成像时对光线的偏折起决定作用的是其表面曲率，而透镜本身的厚度对成像贡献较少。一般的透镜的厚度大，主要是为了满足表面曲率或大孔径的要求。为此可以减少轴向厚度，同时保持表面曲率的不变。
保持曲率不变，下沉每一段曲线，如图所示
常见的切割方法主要有两种，分别是纵切法和横向法。切割方法如下图所示。
在实际加工时，透镜厚度的减少，使得透镜曲面位置提前，使得某一环带的出射点也会相应的提前。因此，实际上不同的环带的出射光线汇聚点会有些许的差别。又由于差别较小，故本实验忽略该微小差别。
透镜表面的圆弧可以由一下的参数方程表示
{ x = x 0 + R cos ( θ ) y = y 0 + R sin ( θ ) \left\{\begin{array}{l} x=x_{0}+R \cos (\theta) \\ y=y_{0}+R \sin (\theta) \end{array}\right. {x=x0+Rcos(θ)y=y0+Rsin(θ)
操作步骤
本实验的操作步骤如下:
1 、根据焦点与半径的关系公式，求解平凸透镜的表面曲率
2 、利用MATLAB编程求解菲涅尔透镜的光学母线数据。
3、将上一步所得到的光学母线导入到solidworks中，进行三维建模得到菲涅尔透镜。将三维模型保存为sat文件。
4 、将sat文件导入到TracePro，添加
光源，接收面。设置菲涅尔透镜的材料，菲涅尔透镜的透光性，光源的属性(是否反射，是否吸收)。
5、进行光学仿真
6、优化分析
7、优化设计
步骤1、2:在焦距为200
mm，透镜的口径为50mm，划分的总段数m为9，每段内共有1000个点。采用纵切法，得到如下的光学母线(程序见附录)
步骤3：将得到的菲涅尔透镜的光学母线导入到solidworks中。由于导入的点是离散的，故solidworks会进行插值，进而导致透镜部分区域为一个点，厚度为0.故对曲线进行平移，使得所得到的菲涅尔透镜具有一定的厚度。
步骤4：因为材料库并没有折射率为1.5
的材料，只能选取折射率最接近1.5的材料。本实验所选取的菲涅尔透镜的材料为塑料pmma，折射率为1.4935.光源为圆形的光源，光线垂直于该平面，为平行光。光通量为100流明。
如图
,发现光线并未聚焦在200毫米处，与理想的焦点的距离差距较大。且光束的最小光斑比较大。环带的数量，所用材料的折射率并非完全等于1.5都会造成焦距的偏离。
材料折射率为1.5的材料难以找到，材料的属性难以修改。所以需要修改环带数，进行优化。
当环带数目为7时，得到如下的结果。
对比两个结果，发现7个环带的菲涅尔透镜的焦距更接近200，且光束的最小口径也相对较小。但受限制与材料，以及solid
work将导入数据进行差值而带来的误差，焦距并不能很好地达到理论值。
参考文献： 张航，严金华. 非成像光学设计. pages 76–84, 2016.
抛物线函数
主程序
%菲涅尔透镜母线主程序
clear;
clc;
close;
%f:焦距
f = 200;
rho = 50;
%透镜半径
R = f/2;
%m: 划分的段数
m = 5;
%num ：每段内 ，点的总数量
num=100;
site3 = split_sections(R,rho,m,num);
% 在solidworks中，竖直线段会被拟合成曲线，选择分段导出坐标
for i=1:m
filename = 'z'+string(i)+'.txt';
z_section = site3(1+num*(i-1):num*(i),1:3);
save(filename,'z_section','-ascii');
end
% 导出坐标到z.txt
save('site_fresnel_5.txt','site3','-ascii');
母线划分函数
%菲涅尔透镜的母线划分函数
% R:平凸透镜曲率半径
% r:球冠半径
% n:分成n段
function site=split_sections(R,rho,m,Point_per_sec)
Delta_r = rho/m;
% 圆心坐标
Cir_cen = [0;-sqrt(R^2-rho^2)];
%求解各分段的边界的角度
theta=zeros(1,m+1);
for i=1:m+1
theta(i) = acos((rho-(i-1)*Delta_r)/R);
end
%% 初始化分配内存
x = zeros(m*Point_per_sec,1);
y = x;
y_theta = Cir_cen(2)+R*sin(theta);
%% 每个section的坐标
for j=1:m
range = linspace(theta(j),theta(j+1),Point_per_sec);
x(1+Point_per_sec*(j-1):Point_per_sec*j) = Cir_cen(1)+R*cos(range);
y(1+Point_per_sec*(j-1):Point_per_sec*j) = Cir_cen(2)+R*sin(range)-y_theta(j);
end
%% 存入矩阵z，画图
site = [x,zeros(m*Point_per_sec,1),y];
hold on
plot(x,y)
yline(0,'--')
axis equal
exportgraphics(gca,"母线_5.jpg")
return
:::
j-1):Point_per_secj) = Cir_cen(2)+Rsin(range)-y_theta(j);
end
%% 存入矩阵z，画图
site = [x,zeros(m*Point_per_sec,1),y];
hold on
plot(x,y)
yline(0,‘–’)
axis equal
exportgraphics(gca,“母线_5.jpg”)文章来源：https://www.toymoban.com/news/detail-678771.html
return文章来源地址https://www.toymoban.com/news/detail-678771.html
未经同意，不能转载。转载相关事宜请联系QQ：2251856840，或者发邮件（2251856840@qq.com）申请。
到了这里，关于菲涅尔透镜的文章就介绍完了。如果您还想了解更多内容，请在右上角搜索TOY模板网以前的文章或继续浏览下面的相关文章，希望大家以后多多支持TOY模板网！
