[ETCartman]

Данная тема назревала довольно давно, создана отчасти как самостоятельная, отчасти как ветвь данной темы.

1) Что такое Code-Aster? Это открытая расчетная программа индустриального класса (уровня ANSYS и пр.), сертифицированная по ISO, изначально созданная и развиваемая как основной расчетный инструмент французской атомной и энегретической отрасли (EDF)
Как и ANSYS она не является собственно строительной, это МКЭ программа общего назначения (машиностроение, строительство и пр.), но к слову говоря в определенном смысле она более строительная чем ANSYS (специальные модели для грунтов, бетона, вантовые элементы с которыми в ансисе большие траблы) и проч. В силу многообразия функций и обширной поддержки она естественно лучше чем другие открытые программы приспособлена для адаптации под разные специальные нужды. То есть таким образом это не продвинутый калькулятор ++ вроде скад-лира, а тяжелый пакет к которому применим глагол "моделировать". То есть воспроизводить реальное поведение конструкций с учетом нелинейностей, а не просто считать внутренние усилия по линейной теории и затем прогонять их через арифметические алгоритмы норм. К слову усилия в Code-Aster тоже естественно можно получать и потом обсчитывать по ним с учетом требований снипов и сп. Хотя как калькулятор она едва ли в конечном итоге окажется удобной.


2) имеющиеся на сайте материалы выложенные в dnl находятся поиском по ключевым словам. Можете просто нажать ссылку ПО для поиска материалов по открытому и бесплатному софту (до кучи).
Официальный сайт и сайты поддерживающих компаний:
http://www.code-aster.org (английский и французский переключается - там же есть активный форум на двух языках).
http://www.salome-platform.org/ (это основной пре-построцессор, тоже открытый)
http://code-aster.de/ ,
http://www.alneos.com
Конечно же в обязательном порядке
http://caelinux.com/CMS/
и его вики
http://caelinux.org/wiki/index.php/Main_Page
(форум там есть). Я еще напишу про этот сайт ниже.
http://necs.fr/ (c виндоус версией которая есть также тут http://sourceforge.net/projects/asterwin/
Русскоязычный форум (не особо активный и продуктивный)
http://freecae.ru/
Проект с попыткой перевода на английский
http://codeaster-intl.sourceforge.net/ (http://sourceforge.net/projects/codeaster-intl/)
там же можете взять текстовый файл интернационализации и переводить на русский, башкирский, татарский и другие языки народов СССР
Интерфейс для Python
http://www.libaster.org/
Ссылок и сайтов еще очень много и я просто оставлю тут место для них (и также, надеюсь, они будут появляться в процессе обсуждения)

3) Как начать изучать и пользоваться? Хотя виндоус версия есть, она прошла большую часть верификационных тестов, но она старая и для открытой программы это критично. Открытая программа - это на самом деле не одна а куча программ, куча форматов файлов которые меняются от версии к версии иногда. Все в совокупности это приведет к тому что вы будете не понимать что и от чего глючит. Кроме того большая часть англоязычных примеров относится к более новым версиям.
К виндоус версии лучше обращаться после того как будут получены уверенные навыки владения в более подходящей системе.
Как установить линукс если вы никогда не имели с ним дело? И какой лучше дистрибутив взять.
дистрибутив однозначно лучше этот
заполняете форму (для статистики), выбираете 11 версию (последняя на сегодняшний день) и кнопку download
далее вариантов много
3.1) прожечь iso образ на dvd и запускать как живой диск (посмотреть можно но работать практически нет).
3.2) установить второй системой на компьютер. это лучший вариант но не всегда. во первых до этого (лучше всего) заранее приготовить места на диске (три раздела для линукс - swap, /home и / . во вторых многие виндоус программы постоянно затирают загрузчик и потом могут быть проблемы с запусканием в том числе и виндоус. это не проблема для оптыных юзверей но не рекомендуется для начинающих.
3.3) Установить на отдельную флешку и загружать ее с usb - неплохой вариант, но тоже требует гугления и некоторых усилий. Тоже опять же загрузчик.
3.4) Скачать виртуальный плеер (виндоус программа) http://www.vmware.com/products/player/ , установить его, подцепить загруженный iso при запуске диалога. Получится линукс работающий из windows. То есть конечно это будет не нормальный линукс (быстрая и удобная система), но что то внешне очень похоже. Проблемы с доступом к дискам, видеокартой и пр. Но так или иначе для начального обучения этого может быть достаточно. Я сам лично давно не пользовался, ввиду постоянного использования линукс как основной системы
К слову такой вариант лучше чем использование виндоус версии еще и потому что на caelinux-е все предустановлено - то есть не только code-aster но и куча других программ которые могут быть полезны с ним

4) Литературы на английском по code-aster к настоящему моменту не так уж и мало. Тут выкладываю файл с примерами по расчету стержневых систем и кабелей. Не смотря на кажущуюся сложность логика работы программы очень простая. даже проще чем у обычных наверно. основана на стандартных форматах сеток (unv, med, mail) и именованных группах узлов и элементов этой сетки. К этим группам прикладываются нагрузки, закрепления, назначаются материалы, сечения и пр. с помощью текстового командного файла который формируется в блокноте или в интеративной базе команд Eficas. Понимание команд и перевод Eficas-а я является ключом к освоению всей программы. К слову помомо стандартных команд язык code-aster предлагает и стандартные функции управления моделью и результатами (как ANSYS APDL) и плюс интерфейс к пайтону вообще делает эту программу клондайком для тех кто имеет соответствующие навыки или не против их заиметь.

5) Соучастниками создания темы являются члены группы разработчиков бесплатного софта (см. самую первую ссылку) и в частности DEM
То что открыл именно я - связано с сдвижкой по времени, потому что когда форум будет активен я буду спать. А потом поеду на работу и приеду когда все участники будут спать.
Актуальным вопросом является трансляция модели из одной из общераспространенных программ (например скад-лира) в достаточную по функциональности открытую (из которых Code-Aster по большому счету самый привлекательный претендент). Однако с моей точки зрения (которая основана на некотором опыте использования программ обеих типов - то есть строительных калькуляторов и тяжелых классических) - путь это почти что тупиковый. Можно будет сравнительно легко транслировать какие то простые вещи. Например сетку. Или если постараться то сетку с группами. Но практически любая реальная модель после трансляции будет требовать такой глубокой доработки что смысл практический это действие потеряет. Это чисто в силу разницы и несовместимости идеологий, математики конечных элементов, типов нагрузок и видов расчетов и так далее. На мой взгляд более простым (и с большим количеством последующих полезных применений) будет являться способ - от абстрактной модели конструкции к любой из программ.
То есть легче будет создать собственный описательный язык конструкции и уже из создаваемой на его основе базы данных транслировать модель и в скад и в code-aster и при надобности куда угодно (в calculix или ANSYS и тд). При таком подходе надо будет записывать текстовый файл в определенном формате - что достигается легко даже с самыми элементарными навыками программирования. Да и идеологически этот подход более верный. Допустим завтра скад вообще аннулирует текстовые файлы. Или изменит их сильно. Практически придется все начинать с нуля. Да и почему именно скад? Цель то инженера - свободно использовать любые доступные инструменты эффективно. Для создания такого описательного языка в принципе не надо писать какую то отдельную программу. Нужно просто взять за основу какой то сеточный формат с группами который легко читается кодом на любом языке и наиболее распростраен в качестве формата для обмена. претендентов на такой формат не много. если брать коммерческие программы то наверно самым удобным будет формат abacus (и он же calculix) или nastran. Другой вариант - MSH, UNV или MED. Я бы выбрал если честно UNV как самый универсальный и несложный. К тому же с большим количеством готовых конвертеров изо всего и во все. Но этот выбор не принципиален и можно даже выбрать несколько таких форматов. Это всего лишь стандартная сетка для численных задач и все. Сетка то что все конечно элементные программы объединяет (хотя может быть и разной). В конечном итоге если выбирается тот же unv то и сетка не важна. Разбивается она в том же GMSH на раз а группы привязаны к геометрическим элементам (пластинам целиком, стержням и так далее).
Собственно база данных состоит из массивов которые могут заполняться хоть где и хоть как. Например в электронных таблицах (что естественно коль речь идет о массивах данных в виде таблиц). таким образом вы можете строить модель в скаде не открывая его совсем. или с небольшими правками. задавать нагрузку сразу на все плиты просто поменяв одну цифру в таблице. Сгустить сетку не трогая нагрузок вообще, сделать параметрическую модель здания и только менять размеры. И сделать такую передачу будет на порядок проще чем выуживать что то откуда то и пытаться приспособить голову бобра к телу осла.
Я потом еще продолжу на эту тему, пока все.

[DEM]

Огромное спасибо за начатую тему, вопросы есть т.к. при первоначальном опыте работы с программой было довольно много вопросов....
По части универсального формата описания модели вы правы, это сейчас очень требуется...
Т.к. начаты проработки именно в этом направлении...

[СергейД]

поскольку меня волнуют только нелинейные расчеты=
пожалуйста укажите ссылку на верификационные задачи и расчеты жб и грунтов.
я попробую порешать для сравнения.

[MrWhite]

Очень любопытно. А что за модели грунтов и бетона? И собственно строителей больше интересует железобетон. Сам бетон не так часто нужен и легко моделируется без специальных моделей.
Просто фраза смутила про то что Code-Aster более строительный чем Ansys, хотя уж в чем в чем, а в модели ж/бетона и грунта у меня к нему претензий нет. Про ванты не компетентен, хотя разница между "большие траблы" и "нет возможности" большая

[100k]

А что там с API? Оно таки есть, это понятно, а подробнее можно...

[DEM]

Сообщение от 100k А что там с API? Оно таки есть, это понятно, а подробнее можно...

НАсколько я знаю есть интерфейс работы с Питоном.....

[MrWhite]

Сообщение от DEM НАсколько я знаю есть интерфейс работы с Питоном.....

Под API обычно понимают несколько иное

[DEM]

MrWhite
Я в курсе, но даю инфу котрая есть, т.к. скорее всего он имел ввиду полноценое API как у Robota. Такого не будет скорее всего....

[ANKR]

Правильно ли я понимаю - версии Code-Aster под видновоз нету ?

[viewpoint]

СергейД, вот ссылка на верификационные задачи с официального сайта...
http://www.code-aster.org/V2/doc/def...ndex.php?man=V
На странице можно переключиться на английский язык, но там машинный перевод со всеми вытекающими...

По железобетону запускаем поиск по ключевым словам "beton arme". Нужны разделы [V5] (динамическая нелинейность) и [V6] (статическая нелинейность). Как пример,
http://www.code-aster.org/V2/doc/def.../v5.06.106.pdf
http://www.code-aster.org/V2/doc/def.../v5.06.107.pdf
http://www.code-aster.org/V2/doc/def.../v5.06.108.pdf
http://www.code-aster.org/V2/doc/def.../v6.02.119.pdf

С грунтами оказалось все сложнее, многие верификационные задачи находятся в закрытом доступе (
Искать надо по ключевому слову "sol".
http://www.code-aster.org/V2/doc/def.../v7.30.100.pdf

MrWhite , о моделях грунтов и бетона, используемых в Code Aster, можно почитать в документации на официальном сайте...
http://www.code-aster.org/V2/doc/def...ndex.php?man=R
[R7.01.11] [R7.01.13] [R7.01.14] [R7.01.15] [R7.01.16] [R7.01.17] [R7.01.18] [R7.01.22] [R7.01.23] - это документация по грунтам;
[R7.01.24] [R7.01.01] [R7.01.03] [R7.01.04] [R7.01.06] [R7.01.07] [R7.01.08] - это документация по бетону;

ETCartman, Так как многие совсем не знакомы с Code Aster нужно акцентировать внимание, что Code Aster - это решатель, к нему необходимы пре-постпроцессоры такие, как Salome или GMSH, ссылки на которые вы уже дали.
До кучи ссылка на русский сайт с небольшими примерами Salome и Code Aster...
http://www.laduga.ru/salome/salome.shtml

ANKR , версия для винды есть, но она всегда отстает на одно поколение от последней версии для Linux. Как написал ETCartman, для опенсорс программы это критично.

[ANKR]

для начала освоения и простых задач я думаю не критично
aster - скачал

[DEM]

Code-Aster скачал и запустил, создание схемы пока не получается, окно создания элементов постоянно скрывается за окном просмотра....
Salome for Win 32-ни под каким соусом не хотит на 64 битной машине запускаться, но это и не мудрено....
Скорее всего буду пробовать Линь на флешку установить и пробовать, это как первоначальный вариант, если все пойдет, то буду ставить вторую систему....

[ANKR]

DEM, Salome for Win 32 у меня на WinXP 64 bit запустилась без проблем. Причем несколько разных версий. Последнюю качал из dnl

[DEM]

ANKR
Попробую с dnl запустить....

[Boxa]

Самый лучший принцип обучения это "делай как я" и после освоения базовых принципов переход к "творчеству".
Не мог бы автор темы наглядно показать (желательно видео) расчет ЖБ здания в упомянутых выше программах.
т.е. на примере 2х этажного Жб домика, от создания расчетной схемы, до получения результатов.
не особо наглея попрошу в линейной и физ. нелинейной постановке.

Заранее спасибо.

[ETCartman]

Версии под виндоус работают плохо, для обучения как раз наименее годится. Сборка salome под виндоус экспериментальная, виснет и вылетает постоянно. Для обучения как раз это плохо. Если под виндоус то лучше однозначно с виртуальным плеером (как я описал в пункте 3) - тогда можно будет спокойно разбирать примеры (и кстати использовать визарды для code-aster в Salome-MECA как в видеопримерах из dnl )

[DEM]

У придется на ноут, ставить Линь на комп....
Тут я вообще полный профан, какую сборку то лучше ставить???

[ETCartman]

лучше сборку caelinux2011 с виртуальным плеером прямо из виндоус. то есть это не линукс но вариант самый простой (я описал в п.3 в топике)

Скачать виртуальный плеер (виндоус программа) http://www.vmware.com/products/player/ , установить его, подцепить загруженный iso при запуске диалога.

Касаемо единого формата данных для всех программ.
Я собственно писал что логика строительных калькуляторов и классических программ очень разная и перегнать данные из строительной в общую сложно по причине разных подходов. Да собственно даже перегонка из скада в лиру дело тоже не простое, потому что в лире есть такие вещи которых в скаде в принципе нет. или из скада в wolsink framework. Последний считает с учетом стесненного кручения, нужны секториальные характеристики и данные о стеснении деформаций по концам. Но в данном случае противоречия не столь обширные и с допущениями экспорт возможен.
Программы типа ANSYS и Code-Aster требуют отдельного блока "материалы". В скад и пр. материал задается при определении сечения (EA, EJ)
также и плотность не отдельным блоком а в составе массива сечений.
Стержневые элементы для корректной работы в двух программах придется даже разбивать по разному. Чтобы иметь возможность правильно приложить нагрузки и тд чтобы потом можно было сравнивать. Это слабое место потому что если использовать форматы типа UNV или MED - там логика разбивки машиностроительная. Группы: в скаде группы есть но не позволяют управлять ими поименно. А номера элементов переданные в разные программы скорее всего будут не совпадать.
Проектировать подобный формат дело интересное. Вот например разработчики CalculiX не парились и сразу просто взяли за основу формат abacus (они его постоянные пользователи). Что очень упростило им жизнь. Если придумывать свой универсальный формат то например только одних разновидностей карт со свойствами сечений будет около десятка (через стандартные характеристики, жб сечение с арматурой, стальное сечение стандартного типа или произвольное - для нелинейного анализа с пластикой). Иначе весь смысл теряется напрочь - никакими преимуществами code-aster воспользоваться будет нельзя если вы упрощаете ввод до скадовского.
Я думаю на данном этапе достаточно просто изучать и переводить aster - потом уже выработка общего формата станет формальностью.

[Кулик Алексей aka kpblc]

ETCartman, отредактируй первый пост (используй тэги [more][/more]), чтобы текст не был на 2 экрана.

[ETCartman]

тут я добавил в архиве программку для виндоус на французском которая что то там делает для code-aster (загрузить сетку, создать командный файл для расчета и группы и выгрузить сетку - по моему так). в общем можете ее добавить в архив материалов. таких вообще говоря тьма тьмущая. можете глянуть кто хочет (там ссылка на сайт есть тоже)

также здесь я проиллюстрирую схему для общего формата данных.
и поясню как она работает.
1) за основу берется любой формат (текстовый или бинарный) который допускает различное представление сеток с узлами и элементами которые могут входить в именованные группы. Ну то есть кто начал читать руководство которое я подцепил в топике тот уже в курсе что я имею в виду.
2) используется общий формат данных в достаточно произвольной форме (на самом деле из анализа данных которые могут быть введены в code-aster, scad, ansys и тд для той или иной задачи. То есть по меньшей мере это какие то таблицы с цифрами. Таблица "линейные свойства материалов" и в колонках: номер типа, имя, модуль юнга, коэффициент пуассона. Другая таблица - "свойства стержневых элементов такого то типа (простая балка тип 5 в скад или универсальный нелинейный элемент Beam189 в ANSYS" Не надо придумывать все и сразу. Важно классифицировать возможные случаи и закодировать таблицы ключевыми словами, чтобы те же характеристики потом не повторялись.
Кстати говоря - может быть таблица исходных данных для сейсмики. Они почти во всех нормах мира одинаковы и сам расчет в сущности арифметическое преобразование результатов модального анализа (который делают все программы). То же относится и к ветровым пульсациям. Но весь цимес в том что можно также будет предусмотреть виды расчетов которых в обычных программах нет. Например данные для расчета на усталость. И так далее
При написании программы эти таблицы будут двумерными массивами. так что тут понятно для чего.
3) Пишется конвертер с группами для любого решателя (ну например для скад или лиры). Для этого естественно нужно будет считать исходную сетку и переписать ее в формате скад (или чего угодно). Причем не важно что скад именованных групп не поддерживает в том же смысле что и code-aster (то есть нельзя к ним в текстовом файле напрямую обращаться). в скаде есть только номера. Но номера есть в том формате (они ассоциированы с именами) который будет выбран в качестве исходного. И например если использовать тот же пайтон - то есть готовые скрипты (открытые) которые эти форматы читают. останется только написать небольшой фрагмент кода для записи в скад.
ЧТо это дает. Почему такой подход лучше чем просто писать конвертер из чего то во что то? Потому что во первых так есть возможность сразу писать код для ковертирования данных в любую программу. И я например могу пособить по части того же calculiX. То есть любой труд в данном случае может быть использован сразу во многих целях. и проще будет отследить что какой конвертер конвертирует а что нет, и как он это делает. Модель введенная и заданная один раз может быть автоматом перегнана сразу во много решателей.
В чем реализовать данный (новый) формат модели? Ну конечно можно написать какой то интерфейс или сразу несколько. можно взять интерфейс электронных таблиц. Если взять тот же опен офис то там пайтон есть как встроенный язык (правда версии теперь разные с выходом Libre4) Можно взять эксель и использовать VBA. Можно сохранить таблицы в txt и читать их чем угодно. В данном случае это не важно. Какое ваше мнение?

[DEM]

ETCartman
Ну распарсить, я смогу почти любой файл(вот пример программы), самое главное чтобы структура была регулярная....
Все что считано из файла, находится внутри программы в виде массива Array[1..20] of TStringList....
Так что переписать все это могу в любом виде....
Нужно только найти подходящий формат и с 10 ток примеров к нему, чтобы можно было составить алгоритм.....
PS. прошу не обращать внимания на ошибки возникающие в работе программы, т.к. делал её для тестирования считывания массива из txt файла.....

Пример файла Скада, тоже прилагаю....

[ETCartman]

Сообщение от DEM ETCartman Ну распарсить, я смогу почти любой файл(вот пример программы), .

Спасибо, очень полезно кстати.
Насчет code-Aster наверно какое то время займет изучение примеров и команд
но что касается формата для универсального конвертера то
то вот http://dip28p.web.fc2.com/calculix/g...vert/index.htm (по ссылке http://www.caelinux.org/wiki/index.php/Proj:UNVConvert )
можно скачать пайтон скрипт который читает формат UNV (I-Deas, поддерживается почти везде)
и вот для питонистов ссылка http://code.google.com/p/meshconvert...=26&spec=svn26
то есть кто знает при этом формат скада - дело за малым - только запись в этот формат
вот статья про меширование

[DEM]

ETCartman
По позже может скачаю, просмотрю алгоритм....
Мне нужны будут примеры, от простых к сложным, чтобы можно было функции отрабатывать(думаю тут вы сможете помочь)...
Правда потом может фигня одна получится, в СКАДе и Лире, номера узлов объемных элементов идут не абы как...
Кстати прога сделана на Лазаре, так что потом надо будет её переделать, и она будет кросплатформенной....

[ETCartman]

Сообщение от DEM ETCartman Правда потом может фигня одна получится, в СКАДе и Лире, номера узлов объемных элементов идут не абы как...

так они и везде не абы как. даже если какая то разница есть - можно соответственно поменять в записи в заданном порядке. не могу пока найти спеку для UNV с рисунками. там все должно быть и в том числе номера. думаю что как раз это все стандартно. кстати в лире есть поддержка сетки GMSH вроде бы с какой то версии они сделали (но по моему без групп).
Примеры то есть только разбирать их надо. Проще всего начать с простого объемного с тетраэдрами тела где минимум характеристик и узловая нагрузка. такое считается в caelinux интерактивно - даже кодастер не надо учить чтобы посчитать (примеры видео есть в dnl)
и с чтения записи UNV формата с группами. Потом оболочки, и балки.

[DEM]

Да я уже смотрел некоторые примеры из dnl там довольно просто все, только разбито на 3-4 файла описание....

[ETCartman]

Да GMSH довольно специфическая программа, но там много функций на самом деле. Есть еще netgen который легче для новичка (загрузить модель и нажать одну кнопку по простомy, в крайнем случае пометить номера нужных граней). И из netgen легко получить UNV с группами который потом хоть куда - править или конвертировать. тут пример такой сетки (посчитанный в CalculiX - 2000 фунтов нагрузка сверху, единицы имперские). Можете открыть этот файл в gmsh например. посмотреть там группы. Попробовать перекинуть в скад и приложить тонну сверху
Я попробую подумать какие могут быть примеры сделать именно для GMSH. Я сам его не сильно хорошо знаю. Но я абсолютно уверен что этот путь самый благодарный. и в том числе (поскольку мне вопросы СНиПов-СП не интересны сейчас, как часть вашего проекта для меня бесполезная) это (общий формат и свободная манипуляция данными) будет полезно и мне - в том числе и со скадом который у меня есть тоже.
Это вообще мог бы быть самый популярный у строителей открытый проект в принципе, не только для русскоязычных пользователей

ps http://www.sdrl.uc.edu/PDF/test_univ...le_formats.pdf все что нашел по unv - там не только сетка а все до кучи на самом деле.
но с другой стороны выбор формата не так и важен потому что почти все они взаимозаменяемы.
вот спек для msh2 (поддерживается тоже в netgen) http://geuz.org/gmsh/doc/texinfo/gms...II-file-format
если можно конвертировать из одного - то можно и из другого

[DEM]

Offtop: ETCartman
Еще бы времени у меня по больше было....
И денег чтобы этим заниматься....
А то на работе грузят всякой фигней.....

[ETCartman]

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

[DEM]

Пара глупых вопросов....
1. Как изменить разрешение экрана в VMware Player, чтобы виртуальная машина, была не середине экрана в маленьком окошке...
2. Как изменить язык в caelinux2011,хотя бы интерфейса...
А то блина что в Убунте ни в зуб ногой, что в VMware Player ни фига не разобрался пока...

[andsan]

Сообщение от DEM Как изменить разрешение экрана в VMware Player

В загруженной caelinux->меню вверху экрана->System->Preferences->Monitors->выбираете разрешение побольше.

[DEM]

andsan
Неее...
Через плеййер так не получалось...
Другой плейер растягивал изображение, тоже выглядело не очень...
В итоге сделал загрузочную флешку и все ОК...
Сейчас разбираюсь с прогами....
Интерфейс конечно а ля СКАД...

[ANKR]

Сообщение от DEM andsan Неее... Через плеййер так не получалось... Другой плейер растягивал изображение, тоже выглядело не очень... В итоге сделал загрузочную флешку и все ОК... Сейчас разбираюсь с прогами.... Интерфейс конечно а ля СКАД...

Я тож хочу так. DEM, какой прогой делал ?

[DEM]

Universal-USB-Installer-1.9.2.5 и unetbootin-windows-583
Любая подойдет...

[andsan]

Сообщение от DEM Неее...

Ну, не знаю. У меня нормально получилось в VMware Player.

[DEM]

andsan
Та не с помощью unetbootin-windows-583 по удобнее будет....

[ANKR]

DEM, на 100 % поудобнее. а главное вообще никакого напряга с созданием флэшки

[DEM]

В общем то потихоньку разбираюсь с созданием схем....
А как нагрузки прикладывать и создавать комбинации нагрузок не фига не пойму...
ETCartman
Вы не могли бы сделать какое то видео с примером расчета простой плиты....

[Yaroshutin]

Очень волнует один вопрос... в документации (R3.08.03) нашел сведения о расчете свойств сечений произвольной формы. Я правильно понимаю, что свойства этих сечений можно назначать стержневым элементам?
И есть ли возможность использовать тонкостенные сечения (набранные из листов, в пределах которых напряжения постоянны по толщине)? Это было бы очень удобно для формирования сечений пролетных строений мостов.

[DEM]

Урок работы с GMSH от ETCartman-а
http://dwg.ru/dnl/11890
http://dwg.ru/dnl/11892

[Chardash]

На сайте ibm.com , есть неплохая информация по Salome и Code Aster.

[iliorik]

Вышла бесплатная книга по Code-Aster
http://framabook.org/14-beginning-with-code_aster/
Очень любопытно.

[Askemann]

Файлы расчета:
https://www.dropbox.com/sh/unb4pmxki...nbPfh0Exa?dl=0

Кто из знающих, может посмотрите на адекватность?
Не знаю как посмотреть деформированную форму в Паравис. Фильтры неактивны.
И хотелось бы посмотреть не всю конструкцию, а отдельные группы из геометрии, возможно так сделать?

Вопрос по контакту. На металлоконструкцию в файле надо сверху поставить такую же. Поставить просто, без креплений (в реальности между ними только слой резиновых уплотнений). В СолидВоркс Симулейшн это вроде задается контактом без проникновения. Подскажите, как задать контакт в этом случае.

Вопрос по моделированию.
Как часть модели из одного файла перенести в другой файл? Надо импортировать (например в .brep), а потом экспортировать в новую модель?

----- добавлено через ~2 мин. -----

Сообщение от ETCartman Русскоязычный форум (не особо активный и продуктивный) http://freecae.ru/

Уже давно не открывается. наверное уже нет его.

[ETCartman]

Нет времени пока посмотреть. последний раз когда я брался за изучение - возникла тоже проблема та же -- с заданием выдачи нужных результатов.Там еще насколько я понял, команды поменялись в последних версиях, поэтому старые примеры с caelinux вики не идут
Предположу (поскольку там опция задания всего и вся очень гибкие и подробные) что надо указывать группы для задания выдачи отдельно.
Я кстати добавил в свой бложек новые материалы по code-aster
И в принципе задача с контактом решается (относительно просто) в calculix (я решал когда то, в принципе оно есть в тестовых примерах к программе) - полегче все чем в code-aster
PS "легкость" в данном случае в смысле освоения. У calculix порог вхождения намного ниже - по сути, те кто знаком с принципом работы консольных приложений вообще проблем не испытывают. В code-aster философия помимо всего прочего построена на интерфейс из пайтон, да и вообще функционально эта программа намного сложней.

[AlexKaz]

Посчитал 18-ти слойный композит в CA, элемент - шелл 4 узла вроде. Первая собственная частота определилась с погрешность 25 % по сравнению с экспериментом (нет в таблице) и твердотелкой в ANSYS. По правде, реальный образец и образец в ANSYS имели дополнительные накладки из текстолита, так что результат неплохой.

----- добавлено через ~14 мин. -----

Сообщение от iliorik Вышла бесплатная книга по Code-Aster http://framabook.org/14-beginning-with-code_aster/ Очень любопытно.

Я писал Жану Пьеру Абри (автору) на одном из интернет форумов, где он консультирует инженеров, на счёт разрешения перевода книжки на русский. Он не ответил, и я не стал интересоваться, позволяет ли лицензия книги перевести её на русский без нарушения авторских прав. Книжка теперь немного устарела, поскольку примеры в ней довольно простые на мой взгляд. Сегодня есть возможность заняться переводом руководств поинтереснее и составить сборник с задачками посложнее. Но только в случае если инженерному сообществу это будет реально необходимо, и за финансирование, которое отобьёт потерянное время. Поколение мальчиков и девочек, выросшее на SW и WB, и не открывавшее ни разу сопромат, как бы не особо тянется мозгами в текстовые интерфейсы - нет мотивации.

Сам CA - дело либо частных микрофирм, которые не могут позволить себе выделить 50 килобаксов на софт (хоррошая такая зп на Западе за 4-6 месяцев, а про нас вообщу молчу), либо узких спецов - а они без труда читают и в оригинале. Ладуга частично перевела на русский руководство, но стиль их перевода - лопатить весь текст целиком. В общем, все надежды только на частные фирмы и их руководство - будет спрос на этот софт, будут и переводы, навески в виде библиотек макросов, и нормальный графический интерфейс взамен мозголомной Eficas. Как раз отсутствие привычного GUI как в других пакетах и делает CA ужастиком для новичков.


По моим ощущениям Code_Aster больше для научников и ориентирован на проведение инновационных исследований, медленнее считает, больше (много больше с непривычки) времени уходит на подготовку задачи, макро-язык вроде APDL нестандартизирован (пользователь вводит в *.comm файл на мой взгляд много лишнего) - пиши макросы сам, для производства толком не предназначен - есть ли в РФ такие инженеры, которые могут использовать потенциал CA?

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

В общем, CA - это дело не фирмачей и рядовых инженеров, это современная большая наука, для которой в стандартных программах нет нужного функционала. А всем остальным дано право пользоваться готовыми вылизаными специализированными пакетами для расчётов конкретных хорошо изученных задач.
AlexKaz вне форума вставить имя Обратить внимание модератора на это сообщение

[AlexKaz]

Пример расчёта композита в Code_Aster со всеми файлами отправил в тему на reaa http://www.reaa.ru/cgi-bin/yabb/YaBB.pl?num=1496578938

[ETCartman]

Я не думаю что микрофирмам нужен Code_aster - для этого есть масса других решений. От бесплатного CalculiX и его дешевых производных типа Mecway до программ типа Solidworks Simulation (в которых есть до 90% функций тяжелых программ), и те что есть сделано удобней и лучше. SWS стоит 5 или 8 тыс, но точно не 50. Есть урезанные ANSYSы - тоже дешево относительно
Это вообще то даже не программа (в привычном понимании), а набор классов для Python - то есть полуфабрикат. На его основе пишутся программы или сервисы (как simscale) https://www.simscale.com/
Освоить CA для обычного инжиниринга очень сложно, если это не что то совсем простое - http://engineering.moonish.biz/efficient/ .

[AlexKaz]

Сообщение от ETCartman Освоить CA для обычного инжиниринга очень сложно, если это не что то совсем простое - http://engineering.moonish.biz/efficient/ .

Работает, но немного странно. Во всяком случае запускается.

По примерам использования - есть товарищ из Японии, читать лучше всего в гугл-транслейт с транслитераций в английский https://translate.google.com/transla...FFE_CHAR_MECA-

На сайте Code_Aster появился специальный раздел с дотошными туториалами, с каждым уровнем сложность растёт http://code-aster.org/spip.php?rubrique68

Сообщение от ETCartman От бесплатного CalculiX и его дешевых производных типа Mecway

CalculiX приживается в среде FreeCAD (кстати, к FreeCAD пытаются прикрутить OpenFOAM). Мы в лаборатории на этой неделе показываем школьниками простейшие опыты по колебаниям и анализу напряжений (оптико-поляризационный метод - гоняем поляризованный свет через полимерные пластинки и получаем разноцветные бабочки). После опытов на оборудовании показываю им во FreeCAD как провести простейший анализ на прочность и собственные частоты. Ребята в итоге могут сравнить картинки напряжений в пластинках в лабе и на компе, увидеть формы колебаний незакреплённой квадратной тонкой пластины (фигуры Хладни по сути). Пластинки и там, и там с острыми вырезами, трещинами, отверстиями.

[AlexKaz]

Видео и статьи в дополнение к ссылкам выше.

https://www.youtube.com/user/t09aki/...f_id=0&sort=dd
https://translate.googleusercontent....8Q8W5jJDhXRrbg - то же, что и http://salome-meca.cocolog-nifty.com/blog/ , только в переводе.

[ETCartman]

Code_Aster можно начать осваивать и даже некоторые задачи сравнительно просто решать, но много минусов и продвигаться очень сложно.
1) под виндоус нет. вернее есть но неофициально и только по запросу. Для полноценного использования нужно устанавливать Salome-MECA, которого под виндоус вообще нет 2) Не программа а набор классов к пайтон. Нужно знать пайтон. Пайтон такой специфический язык (отступы и прочее), что любой найденный где то скрипт с большой вероятностью не запустится без коррекции. 3) CA Меняют синтаксис от версии к версии, старые примеры не работают. Сам он очень очень громоздкий, многоходовый (в отличие от CalculiX) - простое действие в несколько операторов. Даже если разберешься - быстро забудешь (потому народ и публикует примеры).

[MrWhite]

Сообщение от ETCartman Code_Aster можно начать осваивать и даже некоторые задачи сравнительно просто решать, но много минусов и продвигаться очень сложно. 1) под виндоус нет. вернее есть но неофициально и только по запросу. Для полноценного использования нужно устанавливать Salome-MECA, которого под виндоус вообще нет 2) Не программа а набор классов к пайтон. Нужно знать пайтон. Пайтон такой специфический язык (отступы и прочее), что любой найденный где то скрипт с большой вероятностью не запустится без коррекции. 3) CA Меняют синтаксис от версии к версии, старые примеры не работают. Сам он очень очень громоздкий, многоходовый (в отличие от CalculiX) - простое действие в несколько операторов. Даже если разберешься - быстро забудешь (потому народ и публикует примеры).

Не знаком с Code_Aster, но Python не минус - это плюс. Его "специфичность" может осложнить работу только совсем уж новичкам, не втаскивающим в структурное программирование.

[Meknotek]

Сообщение от ETCartman Code_Aster можно начать осваивать и даже некоторые задачи сравнительно просто решать, но много минусов и продвигаться очень сложно. 1) под виндоус нет. вернее есть но неофициально и только по запросу. Для полноценного использования нужно устанавливать Salome-MECA, которого под виндоус вообще нет 2) Не программа а набор классов к пайтон. Нужно знать пайтон. Пайтон такой специфический язык (отступы и прочее), что любой найденный где то скрипт с большой вероятностью не запустится без коррекции. 3) CA Меняют синтаксис от версии к версии, старые примеры не работают. Сам он очень очень громоздкий, многоходовый (в отличие от CalculiX) - простое действие в несколько операторов. Даже если разберешься - быстро забудешь (потому народ и публикует примеры).

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

Про особенности в виде "отступов" можно говорить, если иметь слабое представление о языках программирования и о программировании вообще. (Если бы Вы поработали программистом-джуном хотя бы пару дней, Вам бы сразу объяснили, что и в каком-нибудь С++ КРАЙНЕ ЖЕЛАТЕЛЬНО делать отступы, если хотите вообще работать там). "Отступы" - неотъемлемая часть исходного кода на любом языке программирования, если Вы работаете в команде или хотя бы просто читали "стили оформления" для конкретного языка.
https://en.wikipedia.org/wiki/Indent_style
C style guide
(тысячи их, гуглятся легко, приводить не буду)

[ETCartman]

А там уже достаточно материалов на английском, кстати говоря. Python конечно хороший язык, но при практических расчетах желательно вообще быстро строить модели и задавать граничные условия избегая всякого программирования и запоминания специфических вещей о том или ином языке.
Что кстати там тоже возможно отчасти для простых задач (Eficas), но по совокупности разных вещей программа не проста для начинающего, для распространения и для практической работы. Даже если не считать сложностью наличие ее только под Linux.
В отличие от Calculix например - он сравнительно прост и для того и для другого и тоже может быть использован вместе с Salome-MECA (при использовании конвертеров сеток).

[AlexKaz]

ETCartman, подскажите пожалуйста, Вам попадались free-коды для моделирования процессов обработки металлов давлением? Сейчас осматриваю Free Tochnog http://tochnog.sourceforge.net/ . Code_Aster применительно к Metalforming в сущности не приспособлен: медленный, необходимо извращаться чтобы уничтожать элементы, возможно будут проблемы c перестроением сетки
https://www.code-aster.org/forum2/vi...c.php?id=16068
https://code-aster.org/forum2/viewtopic.php?id=16260

Про-версия Tochnog по цене может потягаться с коммерческими кодами, к слову.

[ETCartman]

В CalculiX решаются такие задачи (то есть учет пластических деформаций, контакт там есть, а больше ничего не надо - сетку адекватно задаче и жесткое нагружение)
Нагружение-разгрузка задается командой *AMPLITUDE
деформированная форма в stl может быть тоже получена в CGX
Процесс как правило сходится до очень значительных деформаций при надлежащей (желательно свипленной) сетке.
Также можно в украинской программе Impact-fem (спонсирована минобразования) - оболочки со слоями (точками интегрирования по толщине и контакт там есть)

[AlexKaz]

Спасибо!

Хочется не упустить нескольких зайцев: быстрый код, разнообразные модели поведения, возможность разрушения. И всё это предназначено для листовых изделий.
Вы ранее выкладывали пример пользовательского материала. Попробую поупражняться со скоростью деформирования. Простейшая модель разрушения (через близкий к нулю модуль Юнга) видимо под вопросом.

Input arguments: ! ! STRAN(NTENS) - total strain at beginning of increment ! DSTRAN(NTENS) - strain increments ! TIME(1) - step time at beginning of increment ! TIME(2) - total time at beginning of increment ! DTIME - time increment ! TEMP - temperature at beginning of increment ! DTEMP - increment of temperature ! PREDEF - interp. values of predef. field variables ! DPRED - increment in predef. field variables ! CMNAME - name given on *material option. ! NDI - no. of direct stress/strain components ! NSHR - no. of shear stress/strain components ! NTENS - total no. of stress/strain components ! NSTATV - no. of state variables ! PROPS(NPROPS) - array with material properties ! NPROPS - no. of material properties ! COORDS(3) - coordinates at point ! DROT(3,3) - rotation increment vector ! CELENT - characteristic element length ! DFGRD0(3,3) - deformation increment at start of step ! DFGRD1(3,3) - deformation increment at end of step ! NOEL - element number ! NPT - integration point number ! LAYER - layer number (not used) ! KSPT - section point number (not used) ! KSTEP - step number ! KINC - increment number Output arguments: ! ! DDSDDE(NTENS,NTENS) - the elastoplastic constitutive tensor ! STRESS(NTENS) - array with stresses ! STATEV(NSTATEV) - state variables** ! SSE - specific elastic strain energy ! SPD - specific plastic dissipation ! SCD - specific "creep" dissipation ! RPL - volumetric heat generation* ! DDSDDT(NTENS) - stress variation due to temperature* ! DRPLDE(NTENS) - variation of RPL due to strain* ! DRPLDT - variation of RPL due to temperature*

В добавок к Tochnog нашлись чешский OOFEM http://www.oofem.org и getfem++ http://getfem.org .

Обновляется CalculiX for CUDA, новая ссылка https://www.openaircraft.com/calculix-extras/

[AlexKaz]

Скомпилировал CalculiX-Extras с поддержкой CUDA. Решатель CUDA работает ну ооочень медленно, процедура библиотеки cusp как-то слишком неторопливо выгружает матрицы из RAM в GPU-RAM. Тестовая задача с помощью Cholmod решается за доли секунды, на CUDA - 36 секунд, фактически 36 секунд идёт выгрузка в GPU RAM и доли секунд - решение. Может автор писал поддержку CUDA только для CFD - не проверял.

# ~/CalculiX/ccx/test_simplebeam $ ./ccx_2.12 -i static ************************************************************ CalculiX Version 2.12, Copyright(C) 1998-2017 Guido Dhondt CalculiX Extras version 2.12, Copyright(C) 2013-2017 Peter Gustafson CalculiX comes with ABSOLUTELY NO WARRANTY. This is free software, and you are welcome to redistribute it under certain conditions, see gpl.htm ************************************************************ You are using an executable made on Вс июл 9 23:33:08 YEKT 2017 The numbers below are estimated upper bounds number of: nodes: 12133 elements: 5583 one-dimensional elements: 0 two-dimensional elements: 0 integration points per element: 4 degrees of freedom per node: 3 layers per element: 1 distributed facial loads: 0 distributed volumetric loads: 0 concentrated loads: 75 single point constraints: 450 multiple point constraints: 1 terms in all multiple point constraints: 1 tie constraints: 0 dependent nodes tied by cyclic constraints: 0 dependent nodes in pre-tension constraints: 0 sets: 4 terms in all sets: 30074 materials: 1 constants per material and temperature: 2 temperature points per material: 1 plastic data points per material: 0 orientations: 0 amplitudes: 2 data points in all amplitudes: 2 print requests: 0 transformations: 0 property cards: 0 STEP 1 Static analysis was selected Decascading the MPC's Determining the structure of the matrix: number of equations 35949 number of nonzero lower triangular matrix elements 1158492 Using up to 1 cpu(s) for the stress calculation. Using up to 1 cpu(s) for the symmetric stiffness/mass contributions. Using CUDA based on CUSP CG SOLVER CUDA v8.0 Thrust v1.8 Cusp v0.4 Assembled stiffness matrix on CUDA device in = 0.469202 seconds Diagnonal preconditioner Preconditioning time = 0.000263 seconds CUDA iterative solver time = 36.0767 seconds Solver converged to 1e-06 relative tolerance after 19364 iterations Using up to 1 cpu(s) for the stress calculation. Job finished # ~/CalculiX/ccx/test_simplebeam $

Написал здесь т.к. тему с CalculiX закрыли.

[AlexKaz]

EDF выложила демо AsterStudy - нового модуля для Salome-Meca, который облегчит набор comm-файла (от EFICAS отказались, в свежем билде этого редактора нет). По планам народный дистр с новинкой соберут в 2018-19.
http://www.code-aster.org/spip.php?article1001
видео на ютубе youtube.com/watch?v=LLghq2NtvJA

[ETCartman]

Я работал с Eficas
Конечно это не просто, не то, чтобы технически, а потому что оно работает в одну сторону. Только с определенной версией Aster и ошибки не распознает и не правит.
Там давно уже назрела необходимость в чем то более продвинутом.

[AlexKaz]

Написал на Гиктаймс заметку о том, как сделать бинарник CalculiX Extras поддержкой CUDA и CHOLMOD.
https://geektimes.ru/post/291475/

[DEM]

AlexKaz
Надо попробовать прогнать тестовую задачку, насколько увеличивается производительность.
Вообще конечно не хватает рускоязычных мануалов по Code-Aster.
Я кстати тоже с Фрикадом баловался, даже расчетные схемы из СКАДа туда пытался импортировать.
Но возникла проблема, для нового СКАДа надо бы API переделать под пайтон.
Я за это может быть опять займусь, может через годик, сейчас надо постепенно реализовывать более денежные проекты, которые позволят развивать те которые пока не приносят дохода.

[AlexKaz]

Сообщение от DEM Вообще конечно не хватает рускоязычных мануалов по Code-Aster.

Пока искал, как можно провести оптимизацию в C_A, перевёл для себя несколько туториалов https://cloud.mail.ru/public/CK1M/q7apYpCz8/
Наверное ещё выложу перевод динамики и на этом остановлюсь.
FEM в целом порядком меня подзадалбал, наверное пора валить с этой карусели.

[DEM]

Сообщение от AlexKaz FEM в целом порядком меня подзадалбал, наверное пора валить с этой карусели.

Offtop:
В строительстве тоже глухо...

[AlexKaz]

Тем временем, французы выпускают Salome-Meca 2017 со встроенным менеджером ключевых слов. Слегка не то, что я ожидал, да и Хэлп к модулю ещё не перевели на Eng. Радует ограниченное количество ключевых слов.
It is free for all the humans

http://forum.dwg.ru/attachment.php?a...1&d=1508415108

[AlexKaz]

Появились видео-уроки работы в новом интерфейсе Salome-Meca+Code_Aster на примере простеньких задач https://www.youtube.com/user/dbpatankar/videos

[ETCartman]

Интересное практическое приложение (на базе Code_Aster)
https://www.code-aster.org/spip.php?article1025
Расчет многоэтажок деревянных, с учетом фактического закрепления и в том числе расчет крепежа

[AlexKaz]

Сообщение от ETCartman В CalculiX решаются такие задачи (то есть учет пластических деформаций, контакт там есть, а больше ничего не надо - сетку адекватно задаче и жесткое нагружение) Нагружение-разгрузка задается командой *AMPLITUDE деформированная форма в stl может быть тоже получена в CGX Процесс как правило сходится до очень значительных деформаций при надлежащей (желательно свипленной) сетке. Также можно в украинской программе Impact-fem (спонсирована минобразования) - оболочки со слоями (точками интегрирования по толщине и контакт там есть)

Тяжело решаются. У меня с ходу не пошла листовая штамповка.
http://forum.dwg.ru/attachment.php?a...1&d=1529234657

----- добавлено через ~4 мин. -----
В LS-Dyna удаётся решить похожую задачу неявным солвером минут за 10-20 на паре ядер.

----- добавлено через ~3 ч. -----
И с сеточкой получше нет сходимости, пока сделаю вывод, что CalculiX отпадает.

[AlexKaz]

Сообщение от ETCartman Интересное практическое приложение (на базе Code_Aster) https://www.code-aster.org/spip.php?article1025 Расчет многоэтажок деревянных, с учетом фактического закрепления и в том числе расчет крепежа

Презентация весной 2018 на эту тему, хорошо так попали в эксперимент.
https://www.code-aster.org/V2/UPLOAD...lis-expose.pdf

[DEM]

Сообщение от AlexKaz Презентация весной 2018 на эту тему, хорошо так попали в эксперимент. https://www.code-aster.org/V2/UPLOAD...lis-expose.pdf

Это вы называете хорошо?
Это всего лишь опять же говорит, о том что расчет на сейсмику, довольно условен.
Зададим пару статических нагрузок, по своей схеме, или наоборот зададим не статическую нагрузку а какой нибудь маятник на этаже поставим, и все пойдет не так как планировалось.

[MrWhite]

Сообщение от DEM Это вы называете хорошо? Это всего лишь опять же говорит, о том что расчет на сейсмику, довольно условен. Зададим пару статических нагрузок, по своей схеме, или наоборот зададим не статическую нагрузку а какой нибудь маятник на этаже поставим, и все пойдет не так как планировалось.

Это можно ко всему остальному отнести, а не только к сейсмике.
А так - есть эксперимент и есть по нему проведенное численное моделирование. Вероятностная оценка - другой раздел.

[СергейД]

появилась необходимость переходить с ansys на "бесплатное".
почитал тему.склоняюсь все же к calculix.
насколько удобно там решаются контактные задачи?
есть ли верификационное руководство?

[ETCartman]

по code-aster на сайте laduga появился полный русский перевод справки
по calculix все есть на вебсайте
обе программы разрабатываются и работают нормально в основном под linux
Под виндоус тоже работают но иногда вылазят ошибки
Учиться можно, работать не удобно

[vanAvera]

Сообщение от СергейД появилась необходимость переходить с ansys на "бесплатное".

Импортозамещение?

[СергейД]

у нас планируются на 2019 работы по импортозамещению МКЭ.
с выпуском отчета=
МКэ в железнодорожных задачах
справка по основам и теории
верификация сравнением с теорией с экспериментами и с ансис -абакус
туториалы
примеры решения практических задач

с применением в договорных работах
и с последующим обучением специалистов отрасли и прочих желающих.

посоветуйте. code-aster смотрел несколько лет назад.
пока что склоняюсь к calculix.

[Бахил]

Offtop: Бесплатный сыр бывает только в мышеловках.

Сообщение от СергейД посоветуйте.

Даже не стоит пробовать ни то ни другое.

[MrWhite]

Сообщение от СергейД у нас планируются на 2019 работы по импортозамещению МКЭ.

Если импортозамещение - тогда почему не Fidesys?

[СергейД]

для второй мыши сыр бесплатный.

+надо=финансы выделяют.

----- добавлено через ~3 мин. -----
фидесис денег немалых стоит. а толку от него мало.

[Boxa]

Сообщение от Бахил Даже не стоит пробовать ни то ни другое.

Мне так же интересны перспективы... Вы не могли бы развить и обозначить аргументы, реально очень интересно.

[MrWhite]

Сообщение от СергейД фидесис денег немалых стоит. а толку от него мало.

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

[Бахил]

Сообщение от СергейД финансы выделяют

Ну если только залезть в исходники и посмотреть что там навояли.

----- добавлено через 54 сек. -----

Сообщение от ETCartman Учиться можно, работать не удобно

Это точно.

[ETCartman]

CalculiX это не столько программа сколько коллекция кодов, им трудно чего то заместить. Дело не в том - хороший или плохой, это смотря в чем и как, а в том чтобы использовать его надо быть фактически со-разработчиком, самому тестировать и верифицировать и так далее. Code-Aster то уже все таки более менее подходит за исключением того что нормальной документации нет даже на английском.
Так что все это пустая трата денег и времени, за исключением если только задача будет поставлена по иному, а не замена шила на мыло
CalculiX хорошо для образования и в странах ЕС многие преподаватели его затем и используют

[СергейД]

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

по сути идет речь о препроцессоре, генерирующем текстовый файл.
юзер может править его потом ручками, если в теме. (как в ls-dyna)
и постпроцессоре, обрабатывающем вывод.

сделать аналогичное в ансис смог бы , но это не финансируется.
а в фидесисе и подобном чужом нереально.

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

[DEM]

Сообщение от СергейД а в фидесисе и подобном чужом нереально.

Да он кажется уже загнулся...

[csp]

Сообщение от СергейД пока пришел к мнению, что удобнее будет задействовать calculix

И что, calculix перестал быть импортным?

[СергейД]

фидесис навязывают. за треть или более цены ансис.
калькуликс- бесплатный.
побеседовал с коллегами- видимо все же коде-астер будем применять как ядро.
препроцессор на основе саломе.

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

[csp]

Сообщение от СергейД видимо все же коде-астер будем применять

Сообщение от СергейД Подмосковье

Вы берега не попутали))))))

[СергейД]

2 csp
посмотрел ваши посты. буквально ни одного без подколов и пр.
вы о чем вообще? тему вообще читали?

2 alexKaz
Прошу написать мне в личку. хотел посоветоваться.

[Boxa]

на всякий случай закину сюда ссылку сегодняшней стать по поводу CalculiX : https://habr.com/post/423359/
Интересно в первую очередь тем, что рассматривается запуск и работа с программой под win.

[Бахил]

Сообщение от Boxa по поводу CalculiX : https://habr.com/post/423359/

Попытался. Это мазохизм какой-то.

[MrWhite]

Сообщение от Бахил Попытался. Это мазохизм какой-то.

Да нельзя работать с бесплатным софтом, это не эффективно как минимум. Обычно и гораздо дороже выходит.

[ETCartman]

вот самые годные примеры по такому стилю работы https://github.com/mkraska/CalculiX-Examples
то есть такой стиль практически бывает не удобен - это больше для скриптования и прочего

[СергейД]

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

вот за сколько коллеги (их фирмы) купили бы оболочку к коде-астеру. ?
с набором примеров верификацонных задач по нелинейностям (грунты, бетон, пластика, контакт)
с примерами типовых задач
с импортом моделей из лиры-скада, с экспортом туда результатов.
с плагинами анализа и оформления результатов?
с сертификатом РААСН со временем.
500 тыр? 900?
какие требования? пишите. хочу понять спрос.

почем купили бы оболочку к ветровой аэродинамике от OPENFOAM?
верификацию сделать сравнением ANSYS и аэротрубой. примеры должны быть.
модели для типовых задач. может быть справочник по результатам из литературы,
к примеру, в софистик встроили за 4900 евро dolfyn (бесплатный слабый продукт)

[DEM]

100-150 тыс руб для начала.
Такими заоблачными ценами вы сразу отпугнете всех.
Сделать варианты подписки.
Сделать вариант расчетов на собственном сервере.
Т.е. моделирование на локальном компе, расчёт на внешнем сервере.
Т.е. простенькие задачи до 1000 элементов расчитывать локально, большие задачи расчитывать на сервере.
Нужно не спрос изучать, он при вашем подходе будет нулевой. А изучать варианты выхода на рынок.

[MrWhite]

Сообщение от СергейД есть выход- сделать бесплатный софт платным... приделав какой-то нормлаьный интерфейс и минимальные удобства... и продавать интерфейс за деньги. какие требования? пишите. хочу понять спрос.

Да это выход. Единственно что сделать нормальный интерфейс, сложнее и дольше, чем расчетное ядро. Банально в разы больше кода.
А так очень много от скорости зависит. Вот даже такой платный и симпатичный продукт как MIDAS GTS мне совершенно не понравился из-за общей медленности.

P.S. А кстати, никто не знает, есть в CalculiX аналог команды Ansys PATH ?

[Ал-й]

Сообщение от MrWhite платный и симпатичный продукт как MIDAS GTS мне совершенно не понравился из-за общей медленности.

А что имеется ввиду? И шла ли речь именно о старом gts, а не о новом gts nx? Который на порядок быстрее...

[СергейД]

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

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

улучшать саломе и пытаться не хочется.
сервис- набор исходников CA для типовых задач.
и отдельные плагины для подготовки правильных исходников со всеми параметрами.

мне этим FLAC3d нравится.
копируешь для подходящей задачи из хелп фрагмент кода (похож на Apdl) и правишь.
и для ls-dyna можно найти примеры.

[ETCartman]

платный софт на основе открытого не новое изобретение, если что
https://www.simscale.com/
https://mecway.com/

[DEM]

СергейД
В фонд Бортника поди заявку собираешься подать
У меня руки не дошли..
Да и не такие там большие деньги а вовлеченность должна быть 100%.
А по части интерфейса, тут ты не прав.
Можно использовать какой нибудь FreeCad и все станет значительно проще.
Главное примеры и возможность работать с коробки.

[СергейД]

разумеется, не я первый это придумал.
специфика=в наборе примеров типовых задач по нелину и в верификации.
причем тут freecad. другие совсем цели. библиотека желдор моделей будет.
это тема РЖД. (выше уже писал)
посмотрим что получится. будет грант- будем делать.

[DEM]

Сообщение от СергейД разумеется, не я первый это придумал. специфика=в наборе примеров типовых задач по нелину и в верификации. причем тут freecad. другие совсем цели. библиотека желдор моделей будет. это тема РЖД. (выше уже писал) посмотрим что получится. будет грант- будем делать.

Фрикад чтобы интерфейс ввода самому не ваять.

[СергейД]

dem!
ты не читал то, что я выше писал.
задача создать набор исходников для коде-астера.
под типовые нелинейные желдор и строит задачи.
с правильными опциями и верфикацией.
вплоть до РААСН, чтобы можно было идти в экспертизу.

для нелина разумно делать небольшие задачи с несложной геометрией, фрагменты
любыми средствами
хоть макросом ansys (нашел такой в интернете), хоть с++

задача рисования моделей не ставится. по умолчанию рисуются в саломе.
с использованием питона.
или читаются из созданной библиотеки.

[DEM]

Сообщение от СергейД dem! ты не читал то, что я выше писал. задача создать набор исходников для коде-астера. под типовые нелинейные желдор и строит задачи. с правильными опциями и верфикацией. вплоть до РААСН, чтобы можно было идти в экспертизу. для нелина разумно делать небольшие задачи с несложной геометрией, фрагменты любыми средствами хоть макросом ansys (нашел такой в интернете), хоть с++ задача рисования моделей не ставится. по умолчанию рисуются в саломе. с использованием питона. или читаются из созданной библиотеки.

Тогда вообще про стоимость можешь не заикаться.
Это будет очень специфическая программа.
И даже ге программа а шаблоны.

----- добавлено через 52 сек. -----
Ну и ветка называется не так и обсуждение в основном другого тут.

[СергейД]

вариант желдор для гранта.
вариант для строительства- для окупаемости.

[MrWhite]

Сообщение от Ал-й А что имеется ввиду? И шла ли речь именно о старом gts, а не о новом gts nx? Который на порядок быстрее...

Не помню. И на порядок быстрее старого - не значит быстро. Мое знакомство с MIDAS закончилось на понимании его скорости и удобстве пользовательского программирования.

[тоннельщик]

Offtop: стоит отметить, что GTS NX действительно очень медленный на задачах очень большой размерности

----- добавлено через ~5 мин. -----
Вообще, коммерциализация наработок под CA + верифицирование РААСН видится вполне себе перспективной. Вполне можно и про стоимость...

[СергейД]

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

лично я сам планировал заняться собственно тестированием и практическими расчетами.
мне для этот красивый интерфейс вообще не нужен. он вторичен и третичен.
этим могут заняться параллельно и даже после верификации другие
мне хватит экспорта из ansys. который уже есть. его нужно проверить и дополнить.

начинать нужно с активных юзеров, котороые бы обосновали начальству полезность.
несколько фирм бы скоординировались и привлекло бы энтузиастов, которые сообща подготовили
мы материалы для начала беседы в раасн.

далее уже вопрос финансовый и организационный. небыстрый. год-два до вериф св-ва РААСН.

но у фрилансера или маленькой фирмы тогда будет недорогой инструмент для серьезных нелинейных расчетов.
на уровне софистика или мидаса.

[gad]

Сообщение от СергейД но у фрилансера или маленькой фирмы тогда будет недорогой инструмент для серьезных нелинейных расчетов. на уровне софистика или мидаса.

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

[СергейД]

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

[DEM]

Сообщение от СергейД я и опрашиваю... хочу узнать, кто уже что-то сделал, чтобы скоординировать действия.

А ты вот у энтузиаста узнай, он вроде тоже один ваяет..
http://selenasys.com/Default.aspx

[AlexKaz]

Речь о сертификации по российским нормам? Они вообще существуют? Code_Aster и CalculiX и так проходят тесты NAFEMS, по крайней мере те из них, которые в этих кодах возможно прорешать.

[MrWhite]

Сообщение от AlexKaz Речь о сертификации по российским нормам? Они вообще существуют? Code_Aster и CalculiX и так проходят тесты NAFEMS, по крайней мере те из них, которые в этих кодах возможно прорешать.

Сертификация по российским нормам это вообще уже к постпроцессорным операциям имеет отношение по сути. Верификация ПО для РААСН это некое дублирование тестов NAFEMS. Не знаю насколько совпадают верификационные примеры.

[AlexKaz]

"но у фрилансера или маленькой фирмы тогда будет недорогой инструмент для серьезных нелинейных расчетов. на уровне софистика или мидаса."

Рано. Для начала софту нужна техподдержка, и даже меньше - тупо обучение. Это литература и уроки по освоению ПО, а они стоят денег (даже за бумагу и издание). А уроков на самом деле на русском языке нет, и литературы нет (не надо строить иллюзий, что всё будет легко и просто "кнопка здесь, кнопка там, и понеслась..."). Если кто впишется - я могу написать методичку. Но никто не впишется.

[gad]

Offtop:

Сообщение от AlexKaz Для начала софту нужна техподдержка, и даже меньше - тупо обучение. Это литература и уроки по освоению ПО, а они стоят денег (даже за бумагу и издание).

Вот за это точно можно и нужно деньги брать.

[СергейД]

речь у меня, разумеется, шла о верификации в РААСН.
получении документа для code-aster., который имеет вес для экспертизы

не исключена ситуация, что всерьез возьмусь за CA.
и хотел бы знать, кто также будет участвовать в процессе. одному это невозможно.
надо 5-7 участников.
пишите в личку, кому интересно.

[AlexKaz]

Сообщение от Boxa на всякий случай закину сюда ссылку сегодняшней стать по поводу CalculiX : https://habr.com/post/423359/ Интересно в первую очередь тем, что рассматривается запуск и работа с программой под win.

Там в статье всё очень тяжело для восприятия. У CalculiX есть нормальный внешний графический лаунчер, где кнопками можно задать нагрузки, материалы, выбрать тип анализа, и сразу запустить на расчёт без ковыряния в коде. Фактически то же самое, только несколько удобнее, реализовано во FreeCAD (и функционал препроцессора продолжает дописываться).

----- добавлено через ~8 мин. -----

Сообщение от СергейД спасибо откликнувшимся. софт- это не более чем ядро. поверх сделать много чего нужно железнодорожного (библиотеки,базы данных, алгоритмы), сделать верификацию на отраслевых задачах, увязать со спецализированными пакетами (не мкэ) нужны сборки и интерфейс, понятные практикам, не владеющим линуксами и тд по сути идет речь о препроцессоре, генерирующем текстовый файл. юзер может править его потом ручками, если в теме. (как в ls-dyna) и постпроцессоре, обрабатывающем вывод. сделать аналогичное в ансис смог бы , но это не финансируется. а в фидесисе и подобном чужом нереально. пока пришел к мнению, что удобнее будет задействовать calculix. для начала нужны верификационные примеры calculix c исходниками. контакт важнее всего.

Разве фидесис не использует CA в качестве решателя? Там какая-то мутная ситуация с софтом. Постпроцессор у них точно забугорный, причём интерфейс местами тупо переведён на русский. Препроцессор, я предполагаю, написан какими-то дизайнерами на пару с несколькими программерами. В любом случае, софт вообще нельзя использовать в строительстве. Это больше напоминает поделку студентов, только красиво обёрнутую.

Новый препроцессор CA в SALOME уже написан самой фирмой EDF. Это лучшее, что могут сейчас написать программеры. Переписывать препост под ваши желания (кнопочный интерфейс) очень муторно из-за специфики формата, и это будет дорого стоить. В РФ Вы не сможете написать препост. Обращайтесь в EDF и заключайте с ними контракт. Если не получается - придётся материться и терпеть (если есть желание юзать CA).

----- добавлено через ~15 мин. -----

Сообщение от СергейД dem! ты не читал то, что я выше писал. задача создать набор исходников для коде-астера. под типовые нелинейные желдор и строит задачи. с правильными опциями и верфикацией. вплоть до РААСН, чтобы можно было идти в экспертизу. для нелина разумно делать небольшие задачи с несложной геометрией, фрагменты любыми средствами хоть макросом ansys (нашел такой в интернете), хоть с++ задача рисования моделей не ставится. по умолчанию рисуются в саломе. с использованием питона. или читаются из созданной библиотеки.

Гранта не хватит, это очень дорого, т.к. придётся заюзать много времени и много людей, которые вообще не будут поначалу вдуплять, что за софт перед ними. Не проще взять какую-нибудь Лиру или что подешевле из русскоязычного, и изваять набор макросов там парой-тройкой сишных программеров?

----- добавлено через ~21 мин. -----

Сообщение от MrWhite Да нельзя работать с бесплатным софтом, это не эффективно как минимум. Обычно и гораздо дороже выходит.

Подтвержаю относительно CA (не весь бесплатный софт неэффективный). CA очень медленный по сравнению с Лирами и прочим вашим строительным специализированным софтом. В разы медленнее ANSYS. Речь не про время решения задачи, а про время подготовки в препроцессоре. Если Вы готовы тратить небесплатное время дорогих спецов - добро пожаловать в OpenSource. CA идеален для научников, но для производственников есть платные более выгодные решения.

[ETCartman]

CalculiX считает по моему быстрее SCAD
Во всяком случае оболочки. Ну или не медленней. Code_Aster не медленней CalculiX из того что я считал
ANSYS, Solidworks Simulation быстрее но не в разы точно. CalculiX с решателем по умолчанию зависит от объема оперативки сильно

[AlexKaz]

В августе я считал автомобильный кофр оболочками, до 100 тыс. элементов. Если Вы захотите формировать не одно поле (например, перемещения), а штук 5 полей результатов, то выгрузка результатов расчёта в *.rmed медленная (иногда доходит до десятков секунд). А решатель работает быстро (для бесплатного проекта), тут у меня нет претензий.
Комп у меня не самый слабый, но RAM не хватает при создании мелкой сетки. Мэшер забивает память до отказа. Процесс не падает, просто (видимо) уходит в swap-файл.
Также Salome периодически падает при большом количестве объектов модели - плоскостей, поверхностей и т.д., в районе нескольких сотен и выше. Похоже на утечку памяти, или повальное сохранение предыдущих действий в препроцессоре в RAM. Это не такой серьёзный глюк, но сохраняться периодически необходимо.

[Бахил]

Сообщение от AlexKaz решатель работает быстро

А исходник открытый есть?

[AlexKaz]

Так прямо на сайте CA и качните. Всё лежит в открытом доступе. https://www.code-aster.org/spip.php?rubrique21

[ETCartman]

для виндовс вроде есть (не пробовал)
https://code-aster-windows.com/download/

[miko2009]

Сообщение от ETCartman ANSYS, Solidworks Simulation быстрее но не в разы точно.

в 6 раз, даже моя древняя K5000 увеличивала на 30% задачи от 40 000 000 неизвестных. При желании на Китайском сервере можно арендовать виртуальную машину (по памяти 2000$ - 168 часов = 21х8 часов) на которую раньше можно было устанавливать "пропатченный" софт и пользоваться вычислительной мощностью сервера, в некоторых задачах прирост мог быть и в 100 раз и более, смотря как сервер загружен, в основном сервера Китайские используют киношники которые обрабатывают компьютерную графику (VFX)

[master_luc]

Сообщение от ETCartman Во всяком случае оболочки.

Прошу прощение, а где в Calculix вы нашли оболочки:
Manual on http://www.dhondt.de/ccx_2.12.pdf page 84 to 112.

[miko2009]

85 стр.

[master_luc]

In CalculiX, S4 and S4R four-node shell elements are expanded into three-
dimensional C3D8I and C3D8R elements, respectively.

[miko2009]

Насколько я понял речь не о полной замене, а о том что можно получить объемный элемент C3D8I за счет " drop the middle nodes" - перемещения узлов.

[master_luc]

Хм-хм...я права даже и не знаю...перевести для вас или оставить в счастливом заблуждение.
В общем и целом: все шелы переходят в Объмники, линейные элементы получают по аналогии с квадратными опусканием узлов по середине ребер.

[miko2009]

Сообщение от master_luc .перевести для вас или оставить в счастливом заблуждение

скорее это в вашем счастливом заблуждении все переходит в объемник. Во время расчета автоматически создаются дополнительные узлы которые задают толщину пластины (и по справке этот тип КЭ рекомендуется применять только при расчете на изгиб) и для вас там даже картинка есть где показывается стык пластин по средней линии см. стр 88 (стык оболочек), дополнительные узлы служат для уточнения нормалей во время разложении матрицы.

[master_luc]

видео , но это ни как не опровергает слова оффициальной документации и код программы, в котором я не нашел пластинок...
"...
6.2.12 Six-node shell element (S6)
This is a general purpose triangular shell element. The node numbering and
the direction of the normal to the surface is shown in Figure 58.
In CalculiX, six-node shell elements are expanded into three-dimensional
wedge elements. The way in which this is done is illustrated in Figure 59. For
more information on shell elements the reader is referred to the eight-node shell
element in the next section.
..."
П.С. Miko2009, вообще-то вопрос был не к вам.

[bigden]

а есть какой-нибудь отечественный открытый решатель с документацией на русском?

[ETCartman]

Сообщение от master_luc Прошу прощение, а где в Calculix вы нашли оболочки:

там несколько типов оболочек. S4,S3,S6,S8 и все то же с R
(Reduced integration)
Могут быть композитные слоистые - с разными материалами по слоям (S6 и S8R)
они же применяются для изгиба с физнелинейностью
Оболочка внутри программы формируется и считается как объемный элемент, что математически эквивалентно "толстой оболочке" по теории - то есть без упрощений. При сопоставлении с другими программами результаты совпадают когда берется толстая оболочка, там где она есть, например Solidworks Simulation
тут еще стоит заметить, что толстая оболочка это просто более точная теория (учет сдвигов) а не физическая толщина. Считать тонкую по теории толстой - Окей, просто больше вычислений.
На русском языке не так много книг по теории соотносимой с компьютерными расчетами, но в принципе кто мешает освоить английский технический?

[Бахил]

Сообщение от bigden а есть какой-нибудь отечественный открытый решатель с документацией на русском?

Улкинсон, Райниш. Линейная алгебра.

[AlexKaz]

В одном из недавних релизов Guido внёс возможность создавать свои элементы. Я ещё не смотрел детали. Знаю лишь, что в прошлых релизах в каждом узле было жёстко задано ровно 3 степени свободы, и код был очень сильно завязан на тройку. Может быть Гвидо всё переписал, и теперь появится шанс добавить привычные балочные и оболочечные элементы.

[СергейД]

уилкинсон, райнш. сборник алгоритмов и программ на языке АЛГОЛ. мир 78?

я перевел на фортран и адаптировал к блочной обработке матриц ВСЕ программы оттуда...
1987-1991 год...
метод RITZIT вошел мне в диплом...

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

но уважаемый БАХИЛ,как обычно, троллит...

[master_luc]

Сообщение от ETCartman Оболочка внутри программы формируется и считается как объемный элемент, что математически эквивалентно "толстой оболочке" по теории - то есть без упрощений.

так и я об этом, что нет там оболочек, а только объемные элементы, что, впрочем, не является крамолой.

Сообщение от ETCartman На русском языке не так много книг по теории соотносимой с компьютерными расчетами, но в принципе кто мешает освоить английский технический?

а кто жаловался на английский или отсутствие литературы: много читаю и много всего есть.

Сообщение от СергейД но уважаемый БАХИЛ,как обычно, троллит...

за Бахилом не заржавеет

[Бахил]

Сообщение от СергейД в каждой программе одна опечатка и одна внесенная ошибка

Что-то ни одной не заметил. Всё очень чётко. Правда не все программы использовал.
Посмотрел исходники астера. Нашёл только разложение матрицы по сингулярным числам в фортране.
Вообще странный продукт. Куча файлов на разных языках. "Я тебя слепила из того что было"(с).

[MrWhite]

Сообщение от Бахил Вообще странный продукт. Куча файлов на разных языках. "Я тебя слепила из того что было"(с).

Ну это как раз не странно, это ожидаемо.

[Vsa125]

Сообщение от Бахил Что-то ни одной не заметил. Всё очень чётко. Правда не все программы использовал. Посмотрел исходники астера. Нашёл только разложение матрицы по сингулярным числам в фортране. Вообще странный продукт. Куча файлов на разных языках. "Я тебя слепила из того что было"(с).

Исходный код здесь https://bitbucket.org/code_aster/cod...44d?at=default

Используется Fortran и Python и небольшие вкрапления C, C++

Вот например матрица жесткости плиты DKQ "codeaster-src / bibfor / elements / dkqrig.F90"

Код:

[Выделить все]

! --------------------------------------------------------------------
! Copyright (C) 1991 - 2017 - EDF R&D - www.code-aster.org
! This file is part of code_aster.
!
! code_aster is free software: you can redistribute it and/or modify
! it under the terms of the GNU General Public License as published by
! the Free Software Foundation, either version 3 of the License, or
! (at your option) any later version.
!
! code_aster is distributed in the hope that it will be useful,
! but WITHOUT ANY WARRANTY; without even the implied warranty of
! MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
! GNU General Public License for more details.
!
! You should have received a copy of the GNU General Public License
! along with code_aster.  If not, see <http://www.gnu.org/licenses/>.
! --------------------------------------------------------------------

subroutine dkqrig(nomte, xyzl, option, pgl, rig,&
                  ener)
    implicit none
#include "asterf_types.h"
#include "jeveux.h"
#include "asterfort/assert.h"
#include "asterc/r8gaem.h"
#include "asterfort/bsthpl.h"
#include "asterfort/dkqbf.h"
#include "asterfort/dkqshp.h"
#include "asterfort/dxqgm.h"
#include "asterfort/dxmate.h"
#include "asterfort/dxqbm.h"
#include "asterfort/dxqloc.h"
#include "asterfort/dxqlocdri1.h"
#include "asterfort/dxqlocdri2.h"
#include "asterfort/dxqlocdri3.h"
#include "asterfort/dxqlocdri4.h"
#include "asterfort/dxqloe.h"
#include "asterfort/dxqloe_NV.h"
#include "asterfort/elrefe_info.h"
#include "asterfort/gquad4.h"
#include "asterfort/jevech.h"
#include "asterfort/jquad4.h"
#include "asterfort/r8inir.h"
#include "asterfort/utbtab.h"
#include "asterfort/utctab.h"
#include "asterfort/utpvgl.h"
#include "blas/dcopy.h"
#include "blas/dscal.h"
    real(kind=8) :: xyzl(3, *), pgl(*), rig(*), ener(*)
    character(len=16) :: option, nomte
!
!     MATRICE DE RIGIDITE DE L'ELEMENT DE PLAQUE DKQ
!     ------------------------------------------------------------------
!     IN  XYZL   : COORDONNEES LOCALES DES QUATRE NOEUDS
!     IN  OPTION : OPTION RIGI_MECA OU EPOT_ELEM
!     IN  PGL    : MATRICE DE PASSAGE GLOBAL/LOCAL
!     OUT RIG    : MATRICE DE RIGIDITE
!     OUT ENER   : TERMES POUR ENER_POT (EPOT_ELEM)
!     ------------------------------------------------------------------
    integer :: ndim, nno, nnos, npg, ipoids, icoopg, ivf, idfdx, idfd2, jgano
    integer :: multic, i, jcoqu, jdepg
    real(kind=8) :: wgt
    real(kind=8) :: df(9)=0.0, dm(9)=0.0, dmf(9)=0.0, dc(4)=0.0, dci(4)=0.0
    real(kind=8) :: df2(9)=0.0, dm2(9)=0.0, dmf2(9)=0.0
    real(kind=8) :: dmc(3, 2)=0.0, dfc(3, 2)=0.0
    real(kind=8) :: bf(3, 12)=0.0, bm(3, 8)=0.0
    real(kind=8) :: xab1(3, 12)=0.0, depl(24)=0.0, caraq4(25)=0.0, jacob(5)=0.0
    real(kind=8) :: qsi=0.0, eta=0.0
    real(kind=8) :: flex(144), memb(64), mefl(96)
    real(kind=8) :: t2iu(4)=0.0, t2ui(4)=0.0, t1ve(9)=0.0
    real(kind=8) :: bsigth(24)=0.0, enerth=0.0, excent=0.0, un, ctor=0.0
    aster_logical :: coupmf=.false., exce=.false., indith=.false.

!
!   LOCAL VARIABLES FOR COEF_RIGI_DRZ

    integer :: j, ii, jj,irot
!    integer :: iishp,jjshp
    integer, parameter :: npgmx=9
    real(kind=8) :: shp(3,4,npgmx), shpr1(3,4,npgmx), shpr2(3,4,npgmx), bb(12,npgmx)
    real(kind=8) :: gshp1(3,4), gshp2(3,4)
    real(kind=8) :: dArea=0.0, gam=0.0, epais=0.0, fact=0.0, gm(3, 4)
    real(kind=8) :: gmemb(4,4), btgmemb(8,4), gmefl(4,12)
!    real(kind=8) :: nm1(8), nm2(8), gm1(4), gm2(4)
    real(kind=8) :: bxb(12,12)

    aster_logical :: dri = .false.


!
    call elrefe_info(fami='RIGI', ndim=ndim, nno=nno, nnos=nnos, npg=npg,&
                     jpoids=ipoids, jcoopg=icoopg, jvf=ivf, jdfde=idfdx, jdfd2=idfd2,&
                     jgano=jgano)
!
    un = 1.0d0
    enerth = 0.0d0

    call jevech('PCACOQU', 'L', jcoqu)
    ctor = zr(jcoqu+3)
    excent = zr(jcoqu+4)
    epais = zr(jcoqu)
    exce = .false.
! COEF_RIGI_DRZ ACTIVE = -1 --> dri = true,  dri =  false sinon
    dri = .false.
    if (ctor .lt. 0.0d0 ) dri = .true.
    if (abs(excent) .gt. un/r8gaem()) exce = .true.
!
!     ----- MISE A ZERO DES MATRICES : FLEX ,MEMB ET MEFL :
    call r8inir(144, 0.d0, flex, 1)
    call r8inir(64, 0.d0, memb, 1)
    call r8inir(96, 0.d0, mefl, 1)
!
!     ----- CALCUL DES MATRICES DE RIGIDITE DU MATERIAU EN FLEXION,
!           MEMBRANE ET CISAILLEMENT INVERSEE --------------------------
    call dxmate('RIGI', df, dm, dmf, dc,&
                dci, dmc, dfc, nno, pgl,&
                multic, coupmf, t2iu, t2ui, t1ve)
!     ----- CALCUL DES GRANDEURS GEOMETRIQUES SUR LE QUADRANGLE --------
    call gquad4(xyzl, caraq4)
!
  if(dri) then
     call r8inir(12*npgmx, 0.d0, shp, 1)
     call r8inir(12*npgmx, 0.d0, shpr1, 1)
     call r8inir(12*npgmx, 0.d0, shpr2, 1)
     call r8inir(12*npgmx, 0.d0, bb, 1)

     dArea = 0.0d0

     call r8inir(12, 0.d0, gshp1, 1)
     call r8inir(12, 0.d0, gshp2, 1)
     call r8inir(12,  0.d0,   gm, 1)
     call r8inir(16,  0.d0,   gmemb, 1)
     call r8inir(32,  0.d0, btgmemb, 1)
     call r8inir(48,  0.d0, gmefl, 1)
     call r8inir(144, 0.d0, bxb, 1)

      epais = zr(jcoqu)
     gam   = abs(ctor)*dm(1)
     do ii = 1, npg
!
!        ----- COORDINATES :
        qsi = zr(icoopg-1+ndim*(ii-1)+1)
        eta = zr(icoopg-1+ndim*(ii-1)+2)
!
!        ----- JACOBIAN AND WEIGHT :
        call jquad4(xyzl, qsi, eta, jacob)
        wgt = zr(ipoids+ii-1)*jacob(1)

!
!        ----- LOOP FOR SHP FUNCTIONS :
!
!        -- ELEMENT AREA :
        dArea = dArea + wgt
!
!        -- COMPUTE LINEAR AND ROTATIONAL SHAPE FUNCTIONS AND DERIVATIVES :
        call dkqshp(qsi, eta, caraq4, jacob, &
                    shp(1,1,ii), shpr1(1,1,ii), shpr2(1,1,ii))

          do j = 1,4
            do i = 1,3
              gshp1(i,j) = gshp1(i,j) + shpr1(i,j,ii)*wgt
              gshp2(i,j) = gshp2(i,j) + shpr2(i,j,ii)*wgt
            end do
          end do
     enddo

     do ii = 1, npg


      do j = 1,4
       do i = 1,3
          shpr1(i,j,ii) = shpr1(i,j,ii) - gshp1(i,j)/dArea
          shpr2(i,j,ii) = shpr2(i,j,ii) - gshp2(i,j)/dArea
       end do
      end do


      do i = 1,4
        j = 3*(i-1)
        bb(1+j,ii) = bb(1+j,ii) - shp(2,i,ii)
        bb(2+j,ii) = bb(2+j,ii) + shp(1,i,ii)
        bb(3+j,ii) = bb(3+j,ii) - 2.d0*shp(3,i,ii)
        bb(3+j,ii) = bb(3+j,ii) - shpr1(2,i,ii) + shpr2(1,i,ii)
      end do
   enddo

  endif
!
    do i = 1, npg

        qsi = zr(icoopg-1+ndim*(i-1)+1)
        eta = zr(icoopg-1+ndim*(i-1)+2)
!        ----- CALCUL DU JACOBIEN SUR LE QUADRANGLE --------------------
        call jquad4(xyzl, qsi, eta, jacob)
        wgt = zr(ipoids+i-1)*jacob(1)
!
!        -- FLEXION :
        call dkqbf(qsi, eta, jacob(2), caraq4, bf)
!        ----- CALCUL DU PRODUIT BFT.DF.BF -----------------------------
        call dcopy(9, df, 1, df2, 1)
        call dscal(9, wgt, df2, 1)
        call utbtab('CUMU', 3, 12, df2, bf,&
                    xab1, flex)
!
!        -- MEMBRANE :
        call dxqbm(qsi, eta, jacob(2), bm)
!        ----- CALCUL DU PRODUIT BMT.DM.BM -----------------------------
        call dcopy(9, dm, 1, dm2, 1)
        call dscal(9, wgt, dm2, 1)
        call utbtab('CUMU', 3, 8, dm2, bm,&
                    xab1, memb)
!
!   compute rotational part of membrane B matrix
!
    if(dri) then
!=====================================================================
! ---  CALCUL DE LA PARTIE MEMBRANE DE LA MATRICE DE MASSE =
! ---  LES TERMES SONT EN NK*NP                                      =
!=====================================================================
!
!        call dkqnim(shp(1,1,i), shpr1(1,1,i), shpr2(1,1,i), &
!                          nm1, nm2, gm1, gm2)

!        do i = 1, 8
!            do j = 1, 8
!                memb(i,j) = memb(i,j) + nm1(i) * nm1(j) * wgt
!                memb(i,j) = memb(i,j) + nm2(i) * nm2(j) * wgt
!            end do
!        end do

!        do iishp = 1, 4
!            do jjshp = 1, 4
!                gmemb(i,j) = gmemb(i,j) + gm1(i) * gm1(j) * wgt
!                gmemb(i,j) = gmemb(i,j) + gm2(i) * gm2(j) * wgt
!            end do
!        end do


!        do iishp = 1, 8
!            do jjshp = 1, 4
!                btgmemb(i,j) = btgmemb(i,j) + nm1(i) * gm1(j) * wgt
!                btgmemb(i,j) = btgmemb(i,j) + nm2(i) * gm2(j) * wgt
!            end do
!        end do
!        -- MEMBRANE (DRILLING PART) Gm:
        call dxqgm(shpr1(1,1,i), shpr2(1,1,i), gm)

!        ----- CALCUL DU PRODUIT GMT.DM.GM
        call dcopy(9, dm, 1, dm2, 1)
        call dscal(9, wgt, dm2, 1)
        call utbtab('CUMU', 3, 4, dm2, gm,&
                    xab1, gmemb)
!        ----- CALCUL DU PRODUIT BMT.DM.GM
        call dcopy(9, dm, 1, dm2, 1)
        call dscal(9, wgt, dm2, 1)
        call utctab('CUMU', 3, 4, 8, dm2,&
                    gm, bm, xab1, btgmemb)

!        ----- CALCUL DU PRODUIT gam/Omega*b(x)b

        do irot = 1, 12
          fact = wgt*gam/dArea * bb(irot,i)
          do jj = 1, 12
            bxb(irot,jj) = bxb(irot,jj) + fact * bb(jj,i)
          end do
        end do


  endif

!
!
!        -- COUPLAGE :
        if (coupmf .or. exce) then
          if(dri) then
!           ----- CALCUL DU PRODUIT BMT.DMF.BF -------------------------
            call dcopy(9, dmf, 1, dmf2, 1)
            call dscal(9, wgt, dmf2, 1)
            call utctab('CUMU', 3, 12, 8, dmf2,&
                        bf, bm, xab1, mefl)
!
!   compute product Gmt.Dmf.Bf
!           ----- CALCUL DU PRODUIT GMT.DMF.BF -------------------------
            call dcopy(9, dmf, 1, dmf2, 1)
            call dscal(9, wgt, dmf2, 1)
            call utctab('CUMU', 3, 12, 4, dmf2,&
                        bf, gm, xab1, gmefl)
          else if (.not. dri) then
!           ----- CALCUL DU PRODUIT BMT.DMF.BF -------------------------
            call dcopy(9, dmf, 1, dmf2, 1)
            call dscal(9, wgt, dmf2, 1)
            call utctab('CUMU', 3, 12, 8, dmf2,&
                        bf, bm, xab1, mefl)
          else
             ASSERT(.false.)

          endif
        endif
    end do
!
    if (option .eq. 'RIGI_MECA') then
        if(.not. dri) then
          call dxqloc(flex, memb, mefl, ctor, rig)
        elseif (dri) then
!     Add rotational to stiffness matrix
!
          ctor=0.d0
          call dxqloc(flex, memb, mefl, ctor, rig)
         call dxqlocdri1(gmemb, rig)
         call dxqlocdri2(btgmemb, rig)
         call dxqlocdri3(gmefl, rig)
         call dxqlocdri4(bxb, rig)
         else
          ASSERT(.false.)
        endif
!
!

    else if (option.eq.'EPOT_ELEM') then
        call jevech('PDEPLAR', 'L', jdepg)
        call utpvgl(4, 6, pgl, zr(jdepg), depl)
       if(.not. dri) then
        call dxqloe(flex, memb, mefl, ctor, coupmf,&
                    depl, ener)
       elseif  (dri) then
!        call dxqloe(flex, memb, mefl, abs(ctor), coupmf,&
!                    depl, ener)
  !     Add rotational to stiffness matrix

          ctor=0.d0
          call dxqloc(flex, memb, mefl, ctor, rig)
         call dxqlocdri1(gmemb, rig)
         call dxqlocdri2(btgmemb, rig)
         call dxqlocdri3(gmefl, rig)
         call dxqlocdri4(bxb, rig)
       call dxqloe_NV(coupmf,rig,depl, ener)
       else
         ASSERT(.false.)
       endif
        call bsthpl(nomte, bsigth, indith)
        if (indith) then
            do i = 1, 24
                enerth = enerth + depl(i)*bsigth(i)
            end do
            ener(1) = ener(1) - enerth
        endif
    endif
!
end subroutine

[Бахил]

Сообщение от Vsa125 Вот например матрица жесткости плиты DKQ "codeaster-src

Всё это прекрасно, только внутренние процедуры типа jevech(), dxqloc() и т.д. пустые.

[iliorik]

https://www.freecadweb.org/wiki/Manu...EM_analyses/ru
Кто нибудь пробовал выполнять так расчёты