Аналитический расчет конструктивных параметров тонкого триплета-апохромата
Разработка метода расчета конструктивных параметров апохромата-триплета в области приближения бесконечно-тонких линз. Определение численных критериев для расчета степени апохроматизации триплета. Рассмотрение критерия исправления вторичного спектра.
Рубрика | Физика и энергетика |
Вид | статья |
Язык | русский |
Дата добавления | 15.01.2021 |
Размер файла | 3,0 M |
Отправить свою хорошую работу в базу знаний просто. Используйте форму, расположенную ниже
Студенты, аспиранты, молодые ученые, использующие базу знаний в своей учебе и работе, будут вам очень благодарны.
Размещено на http://www.allbest.ru/
Московское высшее общевойсковое командное училище
Аналитический расчет конструктивных параметров тонкого триплета-апохромата
Дмитерко Руслан Алексеевич, кандидат физико-математических наук, доцент кафедры ОПД
Страхов Андрей Александрович, преподаватель кафедры ОПД
Россия, Москва
Аннотация
Разработан аналитический метод расчета конструктивных параметров апохромата-триплета в области приближения бесконечно-тонких линз. Для расчета условия апохроматизации применена дисперсионная формула Бучдала (Н.А. БисИБаМ) 3-й степени. Выводятся расчетные формулы для вычисления оптических сил тонких линз триплета из условия минимизации некоторой функции ошибки ББ - условие апохроматизации объектива. Устанавливаются численные критерии для определения степени апохроматизации триплета. Приведены конструктивные параметры оптической схемы триплета-апохромата с использованием стекол с нормальным ходом дисперсии.
Ключевые слова: метод расчета, апохромат, вторичный спектр, тонкий триплет, дисперсионная формула.
Annotation
Dmiterko Ruslan Alekseevich Doctor of Philosophy (Ph.D) in phisics and mathematics associate Professor of OPD Department Moscow higher arms command school
Russia, Moscow Strakhov Andrei Alexandrovich teacher of the OPD Department Moscow higher arms command school
Russia, Moscow
ANALYTICAL CALCULATION OF DESIGN PARAMETERS OF A THIN APO-TRIPLET
An analytical method has been developed for calculating the design parameters of an apo-triplet in the region of infinitely thin lenses approximation. To calculate the achromatization condition, the 3rd degree dispersion formula of H. A. Buchdahl was applied. Calculation formulas for calculating the optical forces of thin triplet lenses are derivedfrom the condition for minimizing a certain SS error function - the lens achromatization condition. Numerical criteria for determining the degree of triplet achromatization are established. The design parameters of the optical scheme of apo-triplet using glasses with a normal course of dispersion are given.
Keywords: calculation method, apo-lens, secondary spectrum, thin triplet, dispersion formula.
Введение
Проблема апохроматизации оптических систем различного назначения остается актуальной и по сей день, несмотря на заметные достижения в этой области. Разработка новых оптических материалов с особым ходом дисперсии позволила облегчить задачу коррекции вторичного спектра. Вместе с тем, некоторые физические свойства и высокие цены особых стекол ограничивают их применение в оптических системах специального назначения и в системах, работающих в экстремальных условиях. Классические материалы (стекла с нормальным ходом дисперсии) в этом смысле являются, в большинстве своем, более подходящими для построения систем с исправленным вторичным спектром. Однако, как показано в ряде работ [1, 3, 4, 10, 11], коррекция вторичного спектра стеклами с нормальным ходом дисперсии затруднительна и связана с необходимостью применения усложненных схемных и конструктивных решений.
Целью данной работы является поиск новых технических решений для создания компактных оптических систем с исправленным вторичным спектром, в том числе, на основе стекол с нормальным ходом дисперсии. Разработанный алгоритм рассматривается на примере тонкого объектива- триплета.
Алгоритм апохроматизации тонкого триплета
В периодической литературе опубликован ряд работ [4, 6, 7], в которых авторы используют представление зависимости показателя преломления оптической среды в виде дисперсионной формулы Бучдала (1). Такое представление дает существенные преимущества по сравнению с известными дисперсионными формулами Шотта, Зельмайера и др. при проведении численных расчетов в плане их упрощения.
где - цветовая координата
X - длина волны, мкм,
Х0 -базовая длина волны (принимается равной 0,574 мкм [6]), s - порядок дисперсионной формулы,
Vi - коэффициент дисперсионной формулы.
Для определения коэффициентов Vi дисперсионной формулы (1) авторами был разработан практический алгоритм расчета. В качестве исходных данных были использованы каталожные значения коэффициентов дисперсионных формул. Исходя из этих данных, были вычислены коэффициенты дисперсионной формулы Бучдала 3-й, 4-ой и 5-й степеней для 500 стекол из каталогов фирм Schott, Ohara и ЛЗОС - основных изготовителей оптического стекла, что позволило при разработке алгоритма апохроматизации тонкого триплета использовать дисперсионную формулу Бучдала.
Алгоритм основан на методе квадратов и, в отличие от работы [9], где приведены соотношения для расчета коэффициентов Vi дисперсионной формулы (1) только исходя из коэффициентов дисперсионной формулы Зельмайера, позволяет рассчитывать указанные коэффициенты для тех оптических сред, которые имеют запись дисперсионных соотношений в виде формул Шотта, Г ерцбергера и Резника.
Запишем выражения для оптической силы тонких линз триплета от цветовой координаты т в следующем виде (здесь используется дисперсионная формула Бучдала 3-ей степени):
Где
В выражениях (2) применяются нормированные значения коэффициентов дисперсионной формулы Бучдала [7], [8]. Для 1 -й линзы они определяются следующим образом (для 2-й и 3-й линз - аналогично):
Оптическую силу тонкого триплета в зависимости от цветовой координаты ю определим как сумму оптических сил составляющих линз:
Условием апохроматизации триплета является независимость его оптической силы от цветовой координаты т. Используя нормировку Ф1+Ф2+Ф3 = 1, составим целевую функцию SS (функцию ошибки), подлежащую минимизации для выполнения этого условия:
Подставим в (5) выражение для функции F(ro) из (4) и после очевидных преобразований запишем функцию SS как степенную функцию от независимой переменной њ:
Учитывая, что функция ошибки SS является сугубо положительной функцией (как квадрат некоторого выражения) одной переменной ю, поставим в соответствие ей некоторое положительное число INT, являющееся площадью фигуры, заключенной между осью абсцисс и функцией SS. Для этого найдем аналитически определенный интеграл от функции SS в пределах от Ю1 = - 0,112602 (голубая линия водорода F) до Ю2 = 0,068237 (красная линия водорода С):
где соответствующие постоянные после проведения интегрирования в указанных пределах равны:
Функция INT(Zi,Z2,Z3) является функцией трех независимых переменных Zi, Z2 и Z3. Точка локального минимума этой функции определяет условия, при которых триплет является апохроматом, так как при этом его оптическая сила в спектральном диапазоне длин волн F - C мало отличается от единицы. Необходимо также отметить, что значение функции INT указывает на степень апохроматизации объектива.
Для определения точки локального экстремума функции INT продифференцируем ее по трем переменным Zi,Z2 и Z3. Для упрощения вычислений предварительно запишем функцию INT в следующем виде: INT(Zi,Z2,Z3) = eiZi2+(e2Z2+e3Z3+e4)Zi+e5Z22+(e6Z3+e7)Z2+e8Z32+e9Z3+ki, (6) где вспомогательные величины еь..е9 определяются выражениями:
Проведя дифференцирование функции ЮТ, записанной в виде (6), получаем частные производные. Для нахождения точки локального экстремума составим систему уравнений:
Решением системы уравнений (7) являются следующие значения:
где вспомогательная величина
Необходимо отметить, что для получения значимых значений величин Zl, Z2 и Zз, все промежуточные вспомогательные постоянные должны вычисляться с не менее чем 30-ю значащими цифрами, так как порядок величины Zo составляет 10-16.
Подставляя вычисленные значения величин Zl*, Z2* и Zз* в выражение (6), можно определить значение функции ЮТ и возможную степень апохроматизации объектива. Предварительно можно определить, что апохроматы имеют значения INT < 1-10"10, а ахроматы соответственно INT > Ы0"8.
Полученные значения величин Z1*, Z2* и Z3* являются исходными данными для определения некоторых начальных значений радиусов кривизны шести оптических поверхностей R1...R6. Дальнейшая оптимизация оптической системы при конечных значениях апертуры и толщин воздушных промежутков (линз) проводится при помощи специализированных программных комплексов, например ZEMAX [5].
Расчет объективов
Приведем результаты расчета некоторых оптических систем с исправленным вторичным спектром в соответствие с разработанным алгоритмом.
Рисунок 1 Результаты оптимизации системы К8-ОК4-ТК12
В качестве пробного теста работы алгоритма рассчитаем известную апохроматическую систему К8-ОК4-ТК12. Вычислим указанные в разделе 2 величины Zl*, Z2* и Zз* (результаты промежуточных вычислений вспомогательных величин опускаем) и определим значение параметра ЮТ, отвечающее апохроматическому решению:
Zl* = - 0.0105 дптр Z2* = 3.1680 дптр Zз* = - 2.1574 дптр Параметр ШТ = 0.589540"12.
Результаты оптимизации системы с апертурой 100 мм и относительным отверстием 1:10 приведены на рисунке 1. Оптимизация проводилась в среде программного комплекса оптимизации оптических систем ZEMAX.
Из представленного расчета видно, что с точки зрения устранения хроматизма применение компонента из стекла К8 не является оптимальным ввиду малой оптической силы данного компонента. Это связано с тем, что частная дисперсия PFe стекла К8 меньше соответствующей частной дисперсии стекла ОК4: для К8: PFe = 0.453754, для ТК12 PFe = 0.454709, для ОК4 PFe = 0,454418.
Для выравнивания оптических сил отрицательных компонентов взамен стекла К8 подберем такое стекло, чтобы его частная дисперсия находилась примерно в середине диапазона ОК4-ТК12. Например, для стекла N-PSK3 частная дисперсия PFe= 0.454503. Произведем соответствующий расчет системы N - PSK3 - ОК4 - ТК12.
Аналогично предыдущему расчету определим значение параметра INT для системы N - PSK3 - OK4 - TK12:
Z1 = - 0.1993 дптр,
Z2 = 3.1724 дптр,
Z3 = - 1.9731 дптр.
Параметр INT = 0.583040"12.
Оптические силы отрицательных компонент выровнялись. Параметр INT уменьшился по величине.
Рассчитаем теперь систему на стеклах с нормальным ходом дисперсии.
Из литературы [2] известно, что триплеты на сверхтяжелых лантановых кронах имеют удовлетворительное исправление вторичного спектра в видимой области. Одной из таких систем является система ТК21-СТК9-ЛК3, найденная с помощью разработанного алгоритма. Расчетные параметры имеют следующие значения:
Z1 = 5.6698 дптр,
Z2 = - 7.8583 дптр,
Z3 = 3.1885 дптр.
Параметр INT = 0.342240"12
Результаты оптимизации системы с параметрами D = 100 мм и F/D = 1:10 приведены на рисунке 2.
Рисунок 2 Результаты оптимизации системы ТК21-СТК19-ЛК3
Из приведенных данных видно, что исправление вторичного спектра в этой системе достаточно для того, чтобы классифицировать ее как супер апохромат. Таким образом, параметр INT адекватно отражает степень апохроматизации (исправления вторичного спектра) оптической системы.
Обсуждение результатов
Предварительные результаты проведенных расчетов показывают, что возможно аналитически вычислить некоторый параметр INT, позволяющий выделить оптические системы с исправленным вторичным спектром (апохроматы). Одновременно рассчитываются оптические силы всех трех компонентов. Необходимо указать, что с помощью приведенных соотношений рассчитываются параксиальные конструктивные параметры и обеспечивается исправление вторичного спектра (при достаточно малом значении параметра INT) только в параксиальной области.
Для исправления хроматизма на широком зрачке приходится учитывать влияние сферохроматической аберрации. Именно эта аберрация, т.е. зависимость сферической аберрации от длины волны, не позволяет в подавляющем количестве случаев исправить вторичный спектр на широком зрачке. Единственным пока обнаруженным исключением (не рассматриваются системы на особых стеклах типа фтор фосфатный крон ОК4) являются системы на стеклах СТК9 (СТК19), представленные в статье и позволяющие конструктору рассчитывать апохроматические системы с широким зрачком.
Приведенный в статье алгоритм может быть распространен на апохроматы, состоящие из четырех и пяти линз. При использовании четырех и более стекол в оптической схеме сильно возрастает как количество возможных вариантов, так и вероятность получения систем с хорошим исправлением вторичного спектра на широком зрачке. В отличие от параксиальной области, где безразлично взаимное положение стекол в системе, в оптических системах с широким зрачком данный фактор может играть решающую роль в решении задачи исправления вторичного спектра.
Дальнейшее развитие алгоритма, на взгляд авторов, должно заключаться в учете в функции ошибки SS некоторых параметров, зависящих от коэффициентов Зейделя S1 (сферическая аберрация) и S2 (кома). Такая запись функции ошибки позволит сразу рассчитывать конструктивные параксиальные параметры тонкого апохромата с исправленной комой, т. е. апланата. апохромат триплет линза
В связи с большим количеством возможных вариантов сочетаний оптических сред представляет интерес дальнейшее развитие алгоритма с чисто математической точки зрения в направлении возможности нахождения аналитического решения некоторого неравенства при задании конкретного значения параметра INT, соответствующего варианту апохроматического решения. Тогда при задании двух конкретных оптических сред возможно определение области, где находится 3-я оптическая среда, отвечающая за существование апохроматического решения.
Заключение
В данной работе представлены аналитические соотношения для вычисления параксиальных конструктивных параметров тонкого триплета-апохромата. Алгоритм с математической точки зрения является по своей сути универсальным и позволяет расширить область своего действия на четырех и пяти линзовые объективы, а также на системы с конечным расстоянием между двумя сложными оптическими компонентами. Для проведения расчета требуется представление показателей преломления оптических сред в виде дисперсионной формулы Бучдала. Критерием исправления вторичного спектра является параметр INT. Представлена оптическая система триплета с исправленным вторичным спектром на основе найденной триады стекол с нормальным ходом дисперсии.
Литература
1. Чуриловский В.Н. Теория хроматизма и аберраций третьего порядка. «Машиностроение», Ленинград, 1968. - 312 с.
2. Попов Г.М. Современная астрономическая оптика. М.; Наука, 1988. - 192с.
3. P.N. Robb, “Selection of optical glasses. 1: two materials”, Appl. Opt. 24 (12), 1864-1877 (1985).
4. R.D. Sigler, “Glass selection for airspaced apochromats using the Buchdahl dispersion equation”, Appl. Opt. 25 (23), 4311-4320 (1986).
5. ZEMAX 12 R2 SP2 IE, S/N 34581. Optical Design Program. 2012 г.
6. M. Bolser. Mercado/Robb/Buchdahl Coefficients - An Update of 243 Common Glasses. Proc. of SPIE, v.4832 (2002), p.525.
7. R.I. Mercado and P.N. Robb, “Calculation of refractive Indices using Buchdahl's Chromatic Coordinate”, Appl. Opt. 22, p. 1198 (1983).
8. Chia-Ling Li and Jose Sasian, ”Adaptive dispersion formula for index interpolation and chromatic aberration correction”, Opt. Expr. 22(1), 1193-1202 (2014).
9. P.J. Reardon and R.A. Chipman, ”Buchdahl's glass dispersion coefficients calculated from Schott equation constant”, Appl.Opt. 28(15), 3520-3523 (1989).
10. J.L. Rayces and M. Rosete-Aguilar, ''Selection of glasses for achromatic doublets with reduced secondary spectrum. I. Tolerance conditions for secondary spectrum, spherochromatism, and fifth-order spherical aberration”, Appl. Opt. 40(31), 5663-5676 (2001).
11. S. Chatterjee and L. Hazra, ''Structural design of cemented triplets by genetic algorithm”, Opt. Eng. 43(2), 432-440 (2004).
Размещено на Allbest.ru
...Подобные документы
Определение мощности электрокалорифера. Осуществление теплового расчета нагревательных элементов. Выбор вентилятора и определение мощности электродвигателя для его привода. Расчет конструктивных параметров нагревательного устройства и сети подключения.
курсовая работа [597,3 K], добавлен 17.01.2012Расчет параметров теплообменивающихся сред по участкам. Обзор основных параметров змеевиковой поверхности. Выбор материалов, конструктивных размеров. Распределение трубок по слоям навивки. Определение параметров кипящей среды и коэффициентов теплоотдачи.
курсовая работа [1,3 M], добавлен 16.08.2012- Определение эффективности действия ударника по преграде и его рациональных конструктивных параметров
Определение рациональных конструктивных параметров кумулятивного снаряда. Расчет изменения скорости схлопывания кумулятивной облицовки, изменения угла схлопывания, скорости элементов кумулятивной струи, изменения диаметра и глубины отверстия в преграде.
контрольная работа [538,8 K], добавлен 04.09.2010 Методы расчета простых и сложных заземлителей в однородной и неоднородной среде. Обоснование необходимости определения показателей надежности при проектировании заземляющих устройств. Выбор метода контроля основных параметров заземляющих устройств.
дипломная работа [2,9 M], добавлен 13.06.2012Задача расчета режима как определение характерных параметров режима, необходимые исходные данные и основные этапы. Особенности метода расчета режима при заданном напряжении в конце и в начале линии электропередач, их отличия, интерпретация результатов.
презентация [470,5 K], добавлен 20.10.2013Расчет конструктивных и технологических параметров поперечно-струйной турбины, водоводов и водоприемника. Количество вырабатываемой электроэнергии за год и объем плотины для гидроэлектростанции, работающей при расходе воды Qн=0,8м/c2 и напоре сопла Нс=6м.
контрольная работа [1,4 M], добавлен 24.09.2013Нахождение главных и конструктивных размеров магнитопровода и обмоток. Проведение электромагнитного и теплового расчета трансформатора. Вычисление параметров трансформатора для определения токов однофазного, двухфазного и трехфазного короткого замыкания.
курсовая работа [566,5 K], добавлен 22.09.2021Особенности конструкций газографитовых ядерных реакторов. Выбор и обоснование основных элементов активной зоны. Расчет бесконечного коэффициента размножения, спектра и ценностей нейтронов в активной зоне. Определение параметров двухгруппового расчета.
курсовая работа [1,5 M], добавлен 14.05.2015Метод расчета параметров измерительного механизма магнитоэлектрической системы, включенного в цепь посредством шунта. Определение мощности вольтметра и амперметра. Измерение активной мощности в цепях трехфазного тока. Выбор измерительной аппаратуры.
курсовая работа [647,1 K], добавлен 26.04.2014Особенности расчета параметров схемы замещения ЛЭП. Специфика выполнения расчета рабочего режима сети с учетом конденсаторной батареи. Определение параметров рабочего режима электрической сети итерационным методом (методом последовательных приближений).
курсовая работа [890,7 K], добавлен 02.02.2011Применение метода контурных токов для расчета электрических схем. Алгоритм составления уравнений, порядок расчета. Метод узловых потенциалов. Определение тока только в одной ветви с помощью метода эквивалентного генератора. Разделение схемы на подсхемы.
презентация [756,4 K], добавлен 16.10.2013Описание конструкций теплообменников и обоснование выбора проектируемого теплообменника. Проведение технологического и гидравлического расчета. Элементы механического расчета: определение параметров обечайки, решетки, выбор и обоснование крышки и опор.
курсовая работа [1,7 M], добавлен 16.11.2012Определение напора и расхода воды для гидроэлектростанции, диаметра рабочего колеса, частоты вращения турбины, высоты всасывания и подбор генератора. Расчет энергетических и конструктивных параметров комбинированной ветроэлектрической энергоустановки.
курсовая работа [166,2 K], добавлен 26.12.2015Методика расчета гидравлической системы с параллельными и последовательными линиями. Определение характеристик простых трубопроводов. Упрощение гидравлической системы. Построение характеристики насоса. Определение параметров рабочих циклов гидросистемы.
учебное пособие [429,5 K], добавлен 06.12.2011Расчет объемов и энтальпий воздуха, а также продуктов сгорания топлива. Тепловой баланс котлоагрегата. Определение параметров теплообмена в топке. Порядок и методика расчета водяного экономайзера, аэродинамических параметров. Невязка теплового баланса.
курсовая работа [220,1 K], добавлен 04.06.2014Расчет площади сечения и формы токоведущей жилы. Оценка зависимости напряженности электрического поля в толще изоляционного слоя. Определение электрических параметров кабеля. Расчет тепловых сопротивлений конструктивных элементов и окружающей среды.
курсовая работа [218,5 K], добавлен 10.01.2015Порядок расчета теоретически необходимого количества воздуха для сгорания топлива. Определение параметров процессов впуска. Вычисление основных параметров процесса сгорания, индикаторных и эффективных показателей двигателя. Основные показатели цикла.
контрольная работа [530,4 K], добавлен 14.11.2010Определение параметров Т-образной схемы замещения трансформатора. Составление полных векторных диаграмм преобразователя для активно-индуктивной и активно-емкостной нагрузок. Расчет изменения вторичного напряжения аналитическим и графическим методами.
задача [229,0 K], добавлен 04.12.2010Расчет выпрямителей с емкостной реакцией нагрузки. Методика расчета ключевых стабилизаторов напряжения. Программные средства моделирования схем источников вторичного электропитания. Алгоритмы счета и программная реализация стабилизаторов напряжения.
дипломная работа [704,4 K], добавлен 24.02.2012Определение фокусных расстояний линз и зеркал, наблюдение и оценка их аберраций. Свойства линз и сферических зеркал превращать расходящиеся гомоцентрические пучки лучей в гомоцентрические сходящиеся пучки, виды аберрации. Формула сферического зеркала.
лабораторная работа [59,3 K], добавлен 20.02.2010