Ошибки мешера в gmesh при обработке сверхтонкостенной геометрии

Доброго времени суток!
Возникла идея написать генератор скриптов .geo для программы gmesh для выгрузки геометрии и сетки конечных элементов из программы для расчета трансформаторов с последующей загрузкой этих сеток в различные решатели и анализа трансформатора на реакции на различные условия. Какие решатели будут использоваться - заранее не известно, но разные, так что к конкретным решателям привязываться не целесообразно. Поэтому, было принято решение выгружать геометрию обмоток послойно - слой обмотки, изоляции, обмотки... Скрипты пишу при помощи нейросети, и она сказала, что в случае, когда заранее не известно, какие решатели будут использоваться и сколько - такой вариант является компромиссным. Но, вероятнее всего, будут использоваться решатели, рассчитывающие нагрев и поля рассеяния. Структурно это должно выглядеть так: производится расчет трансформатора, после запускается скрипт, который извлекает из программы основные параметры трансформатора и на их основе рассчитывает геометрию (расположение в пространстве различных элементов), после чего генерирует автоматически .geo скрипт, который затем открывается в программе gmesh, которая строит сетку конечных элементов.
Во вложении привожу пример того, что скрипт генерирует на данном этапе. Там сердечник, каркас и первая (внутренняя) обмотка. В процессе построения сетки для обмотки возникают ошибки. Как я понимаю, мешеру не нравится сверхтонкая геометрия слоев обмотки (в данном случае 0,4мм обмотка и 0,2мм изоляция), а так же то, что профиля на 0.000001 входят в обмотку. Если между профилями и слоями задать выше указанный (положительный) зазор, ошибка исчезает. Но нейросеть сказала, что такие микроскопические зазоры могут сломать решатель из-за наличия микроскопических объемов с вакуумом. К тому же, для профилей и обмоток строятся разные сетки, а по правильному их нужно соединить в одну, что бы решатель понимал, что эти элементы соприкасаются между собой. Поэтому, я задал отрицательный зазор между профилями и обмоткой, что бы они врезались в обмотку. Мешер не понимает, как строить сетку и ломается. А я не понимаю, как решить эту проблему.
Если кто-то в этом разбирается, посмотрите, пожалуйста, код и модельку. Думаю, если там сама логика не правильная, вряд ли осилю полностью переписать генератор. Но если можно решить эту проблему какими-то настройками, подскажите, буду благодарен. Открывать нужно файл trans_assembled.geo - это "точка входа" в программу. И файл с обмоткой должен называться winding1.geo. Приписки в скобках нужно убрать - я добавил их для ясности.
Так же хочу отметить, что опыта ни в .geo скриптах, ни в использовании решателей у меня нет. Так что я даже не знаю, какая сетка им нужна. И до какой степени детализации нужна модель - тоже. Так что не судите строго. Если это заранее не будет работать - так и скажите.
Заранее благодарю.

Вложения

trans.zip (24.1 Кб, 0 просмотров)

0