Ухин В.

Delta Design для большого количества инженеров становится основным инструментом разработки печатных плат. С появлением в его составе модуля SimPCB стало проще и удобней проектировать устройства с контролем импеданса. Наша команда занимается разработкой этого модуля и на данный момент добилась, как мы считаем, неплохих результатов. Ознакомиться с ними Вы можете в статьях (ссылка), либо непосредственно в программе Delta Design (ссылка). 

Стоп!!! А где же здесь Polar SI9000? - спросит нетерпеливый читатель. Не торопись, друг, история только начинается.

Не так давно один из пользователей обратился к нам со следующей проблемой: необходимо спроектировать копланарную линию передачи с волновым сопротивлением 75 Ом на материале RO4003С с параметрами: толщина диэлектрика 0.203 мм, диэлектрическая проницаемость 3.38. На рисунке 1 представлен фрагмент спецификации материала.

Рис. 1. Фрагмент спецификации на материал RO4003С
Рис. 1. Фрагмент спецификации на материал RO4003С

В SimPCB получены следующие параметры копланарной линии (КМ-1Н) (рис. 2).

Рис. 2. Результат расчета линии КМ-1Н
Рис. 2. Результат расчета линии КМ-1Н

Ширина проводника составляет 3.6 мм, расстояние от полигона до сигнального проводника 0.58 мм.    

SimPCB очень молодое приложение, поэтому специалисты часто перепроверяет полученные в нем значения с помощью других программных средств. В данном случае инженер пересчитал структуру в Polar SI9000 V22.04.00 (линия Surface Coplanar Waveguide 1B). Результат показан на рисунке 3.

Рис. 3. Результат расчета линии Surface Coplanar Waveguide 1B
Рис. 3. Результат расчета линии Surface Coplanar Waveguide 1B

Ширина проводника 4.5 мм, расстояние от полигона до сигнального проводника  0.58 мм. Отличия в геометрии существенные. Разница в ширине проводника составляет  почти 1 мм! Кому верить? У Polar SI9000 есть авторитет, а SimPCB только появился на рынке.

Конечно, данное обращение нас немного расстроило. Неужели мы где-то допустили ошибку? 

Если параметр W1 и W2 из Polar SI900 подставить в SimPCB, то получим следующее значение импеданса (рис. 4).

Рис. 4. Импеданс линии КМ-1Н
Рис. 4. Импеданс линии КМ-1Н

Волновые сопротивления в программах отличаются друг от друга на 3.3 Ом SimPCB импеданс составляет 71.8341 Ом, в Polar SI9000 - 75.13 Ом) или почти на 4.5 процента. Кто-то из читателей может сказать, что здесь нет ничего страшного, так как разница не превышает стандартные 10 процентов. Однако это был принципиальный момент, так как в процессе разработки SimPCB наша команда сравнивала полученные значения, в том числе, и с Polar SI9000. Погрешность не превышала 2 процентов. На рисунке 5 показан фрагмент наших тестов.

Рис.5. Результат тестов
Рис.5. Результат тестов

Polar SI9000 серьезная программа, которой доверяют большое количество специалистов со всего мира, поэтому мы приступили к поиску ошибки в нашем солвере. Проверив всё несколько раз и ничего не обнаружив, было принято решение привлечь независимого эксперта в виде программы Ansys. На рисунке 6 показан фрагмент модели копланарной линии передачи и импеданс, рассчитанный в ней.

Рис.6. Модель линии передачи и значение импеданса
Рис.6. Модель линии передачи и значение импеданса

Итак, получилось три значения волнового сопротивления:

Polar SI9000 - 75.13 Ом;

SimPCB - 71.83 Ом;

Ansys - 71.15 Ом.

Из полученных результатов наша команда сделала вывод, что с большой долей вероятности в Polar SI9000 расчет данной структуры осуществляется неверно. 

Этого заключения нам было недостаточно и мы приняли решение продолжить исследование, а именно, разобраться, где именно ошибается Polar SI9000

Вычислительный модуль Polar SI9000, как и в SimPCB, основан на методе граничных элементов. Эта информация присутствует на сайтеPolar Instruments. При данном методе задается (автоматически или вручную) область, в которой выполняется расчет. Её габаритные размеры должны превышать размер структуры (ссылка). Например, для выше указанной копланарной линии, область может выглядеть так (рис. 7).

Рис. 7. Линия передачи и расчетная область
Рис. 7. Линия передачи и расчетная область

Вокруг структуры находится воздух с диэлектрической проницаемостью равной единице.

Так как же можно проверить данную структуру? Представим, что материал с толщиной H1 тоже воздух, тогда при любом H1 и постоянстве оставшихся параметров линии передачи программа должна выдавать одинаковые значения волнового сопротивления. Проверим это в SimPCB и Polar SI9000. На рисунке 8 показан результат вычислений в SimPCB.

Рис. 8 Результат вычислений в SimPCB
Рис. 8 Результат вычислений в SimPCB

Видно, что Z0 не зависит от H1, что соответствует теории.

Ниже (рис. 9 ) представлен результат, полученный в Polar SI9000.

Рис. 9. Результат вычислений в Polar SI9000
Рис. 9. Результат вычислений в Polar SI9000

Теперь уже нет сомнений в том, что в Polar SI9000 при расчете структуры Surface Coplanar Waveguide 1B присутствует ошибка, так как Z0 зависит от H1. При H1=0.103 мм импеданс составляет 95.8 Ом. При увеличении толщины диэлектрика волновое сопротивление снижается до 81.6 Ом. Примерно такое значение импеданса (81.6 Ом) должно быть при любом значении H1.   

Помним, что волновое сопротивление вычисляется по формуле:

где L - индуктивность линии передачи на единицу длины, С - емкость линии передачи на единицу длины.

Высокий импеданс при H1=0.103 мм говорит о том, что у структуры первичный параметр C меньше, чем, например, при H1=1.003 мм. Ёмкость зависит в том числе и от габаритных размеров опорного слоя. То есть можно предположить, что горизонтальный размер расчетной области структуры Surface Coplanar Waveguide 1B меняется в зависимости от Н1, и при малом H1 в область попадает лишь небольшая часть проводников с нулевым потенциалом. Ниже представлены поясняющие рисунки 10 и 11.

Рис. 10. Расчетная область при малом H1
Рис. 10. Расчетная область при малом H1
Рис. 11. Расчетная область при большом H1
Рис. 11. Расчетная область при большом H1

С ростом H1 увеличивается площадь полигона в расчетной области, что приводит к повышению емкости, снижению индуктивности и импеданса.  Если теория верна, то получается, что при малых значениях H1 и при больших значениях W1, W2 или D1 опорная плоскость для структуры Surface Coplanar Waveguide 1B  представляется в виде проводника, ширина которого может быть много меньше ширины сигнального.  

Проверить это достаточно просто, используя только Polar SI9000, так как в нем присутствует следующая модель линии передачи (рис. 12).

Рис. 12 Модель линии передачи Surface Coplanar Strips 1B
Рис. 12 Модель линии передачи Surface Coplanar Strips 1B

Подберем размер G1=G2 такой, чтобы импедансы структур Surface Coplanar Waveguide 1B и Surface Coplanar Strips 1B совпадали. При этом все остальные параметры одинаковые в обеих линиях передачи. 

Итак, пусть H1=0.203 мм, Er1=1, W1=W2=4.5 мм, D1=0.58 мм, T1=0.018 мм. Результаты расчета представлены на рисунках 13 и 14.

Рис. 13 Результат расчета структуры Surface Coplanar Waveguide 1B
Рис. 13 Результат расчета структуры Surface Coplanar Waveguide 1B
Рис. 14 Результат расчета структуры Surface Coplanar Strips 1B
Рис. 14 Результат расчета структуры Surface Coplanar Strips 1B

Параметр G1=G2=2.85 мм. Получается, что горизонтальный размер вычисляемой области для линии Surface Coplanar Waveguide 1B составляет примерно 4.5 мм + 2*0.58 мм + 2*2.85 мм = 11.36 мм, а полигон представляется в виде двух проводников шириной 2.85 мм. Уменьшим H1 до 0.103 мм и найдем значение G1=G2. На рисунках 15 и 16 представлен результат.

Рис. 15. Расчет структуры Surface Coplanar Waveguide 1B
Рис. 15. Расчет структуры Surface Coplanar Waveguide 1B
Рис. 16. Расчет структуры Surface Coplanar Strips 1B
Рис. 16. Расчет структуры Surface Coplanar Strips 1B

Ширина полигона стала ещё меньше, так как G1=G2=1.1 мм. При этом размер области около 7.86 мм. 

Несложно заметить, что увеличение толщины диэлектрика на 0.1 мм приводит к увеличению расчетной области примерно на 1.75 мм. То есть, если диэлектрик будет толщиной 0.303 мм, то размер полигона составит 4.6 мм с каждой стороны. Проверим это выполнив расчет для двух структур. Результат показан на рисунках 17 и 18

Рис. 17. Расчет структуры Surface Coplanar Waveguide 1B
Рис. 17. Расчет структуры Surface Coplanar Waveguide 1B
Рис. 18. Расчет структуры Surface Coplanar Strips 1B
Рис. 18. Расчет структуры Surface Coplanar Strips 1B

Наше предположение оказалось верным. Графическое представление линии Surface Coplanar Waveguide 1B при определенных значениях H1, D1, W1 и W2 не соответствует структуре, для которой выполняется расчет. Рисунок предполагает, что размер полигона должен превышать ширину сигнальной линии. На практике же оказывается, что это не всегда так и зависит от геометрических размеров ячейки. 

Результаты исследований были отправлены пользователю с заключением, что при решении данной задачи он смело может довериться модулю SimPCB

Какие же напрашиваются выводы? 

  1. Разработчикам вычислительных программ необходимо хотя бы в первом приближении знакомить пользователей с методами и принципами расчета, что повысит доверие к продукту и позволит специалистам избежать ошибок. 

  2. В Polar SI9000 необходимо с осторожностью рассчитывать копланарные структуры.

Многие вопросы возможно решить используя обратную связь. В представленном случае специалист получил помощь от разработчиков SimPCB и подтверждение качества приобретенного продукта.

Комментарии (11)


  1. Flammmable
    30.08.2024 10:34

    Попробовал только что на своём Si8000.
    Действительно, выставив Er в единицу и меняя толщину H1 импеданс начинает стремительно уплывать.

    Позорище-то какое... Не ожидал такого от Polar.

    А вы - красавчики!


    1. Viktor_Uhin
      30.08.2024 10:34

      Спасибо :)


      1. Flammmable
        30.08.2024 10:34

        У меня есть мечта.
        Взять медную пластину толщиной в 5-10мм.
        Вырезать из неё модель дорожки трапециевидного сечения. А также сделать модели боковых полигонов и нижнего полигона.
        Полирнуть всё пастой ГОИ, чтоб блестело на фотках.
        Выфрезеровать из толстенного оргстекла модель конформной паяльной маски - толщиной 3-4мм.
        И ещё вырезать несколько листов из разных диэлектриков - миллиметров по 30-40.

        Сделать это всё пересобираемым на шестигранных винтах и добавить туда коаксиальных разъёмов.

        Взять surf-van, вроде такого:

        ...и отправится в Новосибирск, где хранится ГЭТ 75-2023, государственный первичный эталон единицы волнового сопротивления.

        Потому, что как на конференции кого ни спросишь (из разработчиков абсолютно любого симулятора любых физических процессов):
        - А вы проверяли ваш симулятор на реальных объектах?, звучит ответ
        - Вы знаете, нет, но мы сравнивали с ANSYS, а на него не жаловались. Ну вроде бы.

        Есть что-то глубоко порочное в такой практике.


        1. Viktor_Uhin
          30.08.2024 10:34
          +1

          Тачка огонь!!!

          Сравнение с реальными объектами конечно лучше. Я бы даже сказал: эксперимент всему голова!!! Мы, кстати, делали такое для нескольких типов линий. Измерения проводил РЕЗОНИТ, и мы сами через КСВ. Результат хороший, в 10 процентов укладывались. В одной из статей есть эта информация.


          1. Flammmable
            30.08.2024 10:34

            Проблема со, скажем так, рабочими печатными платами в том, что там начинается влияние технологии производства.
            И трапециевидность может быть не такой уж трапециевидной.
            И шаг плетения стекловолокна уже даёт о себе знать.

            "Пристреливать" симулятор по дорожкам на реальной печатной плате, импеданс которых и так может гулять на ±5-10% - неблагодарное дело.

            А вот если бы механообработкой сделать большую референсную модель (модели), на которой куда стабильней и предсказуемей были бы и геометрические и электрические параметры - вот это было бы круто.

            Мне кажется, на выставке от такого стенда потенциальные заказчики бы пис... ээ... были бы в полном восторге. Забронить векторный анализатор. И повесить Г-образный ключ на цепочке (или сделать вообще на магнитных клипсах), чтобы посетители стенда могли бы сами собрать из такого конструктора модель, просимулировать её и тут же промерить.

            P.S. Стенд оформить в стиле сёрф-ретро, конечно же )))))


  1. Karlson_rwa
    30.08.2024 10:34

    Вы, конечно, молодцы, но немного лукавите, пользуясь, вероятно, крякнутой версией полара. В актуальной 24-й версии вот такой набор dll: CalcEngineBEMDll.dll v23.05.10.0 / PSec.dll v22.04.06.0 / PolarGraphDll.dll v20.10.0.0 / PEGRP32F.dll v 8.0.0.62 / TRCPlus.exe v24.01.01 , что, очевидно, слегка отличается от имеющегося в 22-й версии на пару мажорных версий. Вот повторить бы ваш эксперимент с последней версией полара, вот тогда действительно молодцы без всяких оговорок.


    1. Skat-Pro Автор
      30.08.2024 10:34

      Не какого лукавства, работаем с тем что есть. Если дадите нам доступ к более свежей версии Polar, будем вам признательны. Проверим расчеты в более свежей версии.


      1. Karlson_rwa
        30.08.2024 10:34

        Если б у меня была 24 версия, я бы с радостью показал бы скриншот расчета из неё. Но, увы :)


  1. sdy
    30.08.2024 10:34

    На моей 22 версии Si9000, импеданс с выставленным Eг = 1 не зависит от H1

    Не знаю, что там и как считает, но если Polar считает плохо, то надо бы все таки для начала в их техподдержку обратиться, а потом уже выводы делать


    1. Viktor_Uhin
      30.08.2024 10:34

      Здорово, если они исправили ошибку. А какая версия Polar у Вас? Можете написать? Мы использовали вот эту: Polar SI9000 V22.04.00. И еще, если не сложно, пришлите пожалуйста результат вашего расчета картинкой из Polar, случай, который представлен на рисунке 9 данной статьи.


    1. Flammmable
      30.08.2024 10:34

      Скриншот можете приложить?