(U4.81.04) Оператор CALC CHAMP
Содержание
- 1 1 Цель
- 2 2 Синтаксис
- 2.1 2.1 Операнды RESULTAT/MODELE/CHAM_MATER/CARA_ELEM/EXCIT
- 2.2 2.2 Выбор сеток связанных с вычислениями
- 2.3 2.3 Выбор последовательных номеров
- 2.4 2.4 Местоположение Dans
- 2.5 2.5 Зависимость областей
- 2.6 2.6 Операнды для опций механики
- 2.6.1 2.6.1 Вычисление опций напряжений (операнд FORCED)
- 2.6.2 2.6.2 Вычисление опций деформаций (Операнд DEFORMATION)
- 2.6.3 2.6.3 Опции извлечения внутренних переменных (операнд VARI_INTERNE)
- 2.6.4 2.6.4 Вычисление опций энергии (операнд ENERGIE)
- 2.6.5 2.6.5 Вычисление опций критериев (операнд CRITERES)
- 2.6.6 2.6.6 Вычисление опции гидравлических течений (операнд HYDRAULIC)
- 2.7 2.7 Операнды термических опций
- 2.8 2.8 Операнды для акустических опций
- 2.9 2.9 Операнд для узловых сил и реакций
- 2.10 2.10 Вычисления в области пользователя
- 2.11 2.11 Операнд TITRE
1 Цель
Создать или добавить результат вычисляя области поэлементно или при помощи узлов ( вынужденных, напряжений,…).
Полученный результат либо создается, либо модифицируется, то есть вызов CALC_CHAMP выполняется следующим образом:
resu = CALC_CHAMP ( RESULTAT = resu…, reuse = resu,…) или resu1 = CALC_CHAMP ( RESULTAT = resu,…)
2 Синтаксис
resu [*] = CALC_CHAMP ( ◊reuse = resu, ◊MODELE =mo , [model] ◊CHAM_MATER =chmater , [cham_mater] ◊CARA_ELEM =carac , [cara_elem] ◊EXCIT = _F ( ♦CHARGE = l_charge, [l_char_meca] ◊ / COEF_MULT = cm, [R] / FONC_MULT = Fm, [function] ◊TYPE_CHARGE = “FIXE', ) ◊#Sélection сеток связанных с вычислениями /TOUT = 'OUI”, [DEFECT] / | GROUP_MA =l_grma , [l_gr_maille] | NET =l_mail , [l_maille] ◊#Sélection последовательных номеров /TOUT_ORDRE = 'OUI', /NUME_ORDRE =l_nuor , [l_I] /LIST_ORDRE =l_nuor , [listis] /NUME_MODE =l_numo , [l_I] /NOEUD_CMP =l_nomo , [l_K16] /NOM_CAS =nocas , [K16] /INST =l_inst , [l_R] /FREQ =l_freq , [l_R] /LIST_INST =l_inst , [listr8] /LIST_FREQ =l_freq , [listr8] ◊ | CRITERE = / “RELATIF”, [DEFECT] / “ABSOLU”, | PRECISION = / prec, / 1.0E-6, [DEFAUT] #options для линейных механических результатов ♦RESULTAT =resu , #options вычисления напряжений (элементы непрерывных сред 2D и 3D) ◊CONTRAINTE = | 'EFGE_ELGA' | 'EFGE_ELNO' | 'EFGE_NOEU' | 'SIEF_ELGA' | 'SIEF_ELNO' | 'SIEF_NOEU' | 'SIGM_ELGA' | 'SIGM_ELNO' | 'SIGM_NOEU' | 'SIPM_ELNO ' | 'SIPO_ELNO' | 'SIPO_NOEU' | 'SIRO_ELEM' #options вычисления деформаций ◊DEFORMATION = | 'DEGE_ELGA' | 'DEGE_ELNO' | 'DEGE_NOEU' | 'EPME_ELGA' | 'EPME_ELNO' | 'EPME_NOEU' | 'EPSG_ELGA' | 'EPSG_ELNO' | 'EPSG_NOEU' | 'EPSI_ELGA' | 'EPSI_ELNO' | 'EPSI_NOEU' | 'EPVC_ELGA' | 'EPVC_ELNO' | 'EPVC_NOEU' #options вычисления энергии ◊ENERGIE = | 'DISS_ELEM' | 'DISS_ELGA' | 'DISS_ELNO' | 'DISS_NOEU' | 'ECIN_ELEM' | 'ENEL_ELEM' | 'ENEL_ELGA' | 'ENEL_ELNO' | 'ENEL_NOEU' | 'EPOT_ELEM' | 'ETOT_ELEM' | 'ETOT_ELGA' | 'ETOT_ELNO' | 'ETOT_NOEU' #options вычисления ◊CRITERES criteria = | 'EPEQ_ELGA' | 'EPEQ_ELNO' | 'EPEQ_NOEU' | 'EPMQ_ELGA' | 'EPMQ_ELNO' | 'EPMQ_NOEU' | 'SIEQ_ELGA' | 'SIEQ_ELNO' | 'SIEQ_NOEU' #options интерполяции и извлечения внутренних переменных ◊VARI_INTERNE = | 'VARC_ELGA' #options для нелинейных результатов (получаются при использовании STAT_NON_LINE или DYNA_NON_LINE): ♦RESULTAT =resu , / [evol_noli] #options вычисления напряжений (элементы непрерывных сред 2D и 3D) ◊CONTRAINTE = | 'EFGE_ELGA' | 'EFGE_ELNO' | 'EFGE_NOEU' | 'SIEF_ELNO' | 'SIEF_NOEU' | “SIGM_ELGA” | 'SIGM_ELNO' | 'SIGM_NOEU' | 'SIPO_ELNO' | 'SIPO_NOEU' | 'SIRO_ELEM' #options вычисления деформаций ◊DEFORMATION = | 'DEGE_ELGA' | 'DEGE_ELNO' | 'DEGE_NOEU' | 'EPFD_ELGA' | 'EPFD_ELNO' | 'EPFD_NOEU' | 'EPFP_ELGA' | 'EPFP_ELNO' | 'EPFP_NOEU' | 'EPME_ELGA' | 'EPME_ELNO' | 'EPME_NOEU' | 'EPMG_ELGA' | 'EPMG_ELNO' | 'EPMG_NOEU' | 'EPSG_ELGA' | 'EPSG_ELNO' | 'EPSG_NOEU' | 'EPSI_ELGA' | 'EPSI_ELNO' | 'EPSI_NOEU' | 'EPSP_ELGA' | 'EPSP_ELNO' | 'EPSP_NOEU' | 'EPVC_ELGA' | 'EPVC_ELNO' | 'EPVC_NOEU' #options вычисления энергии ◊ENERGIE = | 'DISS_ELEM' | 'DISS_ELGA' | 'DISS_ELNO' | 'DISS_NOEU' | 'ENEL_ELEM' | 'ENEL_ELGA' | 'ENEL_ELNO' | 'ENEL_NOEU' | 'ETOT_ELEM' | 'ETOT_ELGA' | 'ETOT_ELNO' | 'ETOT_NOEU' #options вычисления of ◊CRITERES criteria = | 'DERA_ELGA' | 'DERA_ELNO' | 'DERA_NOEU' | 'ENDO_ELGA' | 'ENDO_ELNO' | 'ENDO_NOEU' | 'EPEQ_ELGA' | 'EPEQ_ELNO' | 'EPEQ_NOEU' | 'EPMQ_ELGA' | 'EPMQ_ELNO' | 'EPMQ_NOEU' | 'INDL_ELGA' | 'PDIL_ELGA' | 'SIEQ_ELGA' | 'SIEQ_ELNO' | 'SIEQ_NOEU' #options интерполяции и извлечения внутренних переменных ◊VARI_INTERNE = | 'VAEX_ELGA' ♦NOM_VARI = (cf [#2.4.3.]) | 'VAEX_ELNO' ♦NOM_VARI = (cf [#2.4.3.]) | 'VAEX_NOEU' ♦NOM_VARI = (cf [ #2.4.3. ]) | 'VARC_ELGA' | 'VARI_ELNO' | 'VARI_NOEU' #options вычисления гидравлических течений (элементы THM) ◊HYDRAULIQUE = | 'FLHN_ELGA' термические #options ♦ RESULTAT =resu , / [evol_ther] ◊THERMIQUE= | 'FLUX_ELGA' | 'FLUX_ELNO' | 'FLUX_NOEU' | 'HYDR_NOEU' | 'SOUR_ELGA' | 'ETHE_ELEM' акустические #options ♦RESULTAT =resu , / [acou_harmo] ◊ACOUSTIQUE = | 'PRAC_ELNO' | 'PRAC_NOEU' | 'PRME_ELNO' | 'INTE_ELNO' | 'INTE_NOEU' #options для обобщенных сил узловых реакций ♦RESULTAT =resu , ◊FORCE = | 'FORC_NODA' | 'REAC_NODA' #calcul области пользователя ◊CHAM_UTIL = _F ( ♦NOM_CHAM =ncham , ♦ / CRITERE= | 'VMIS', | 'INVA_2' | 'TRACE', /FORMULE =l_forme , [formula] ♦NUMERIQUE_CHAM_RESU=numeric , [I] ), ◊TITRE = title , [l_Kn] ◊INFO = / 1, [DEFAUT] /2 , )
2.1 Операнды RESULTAT/MODELE/CHAM_MATER/CARA_ELEM/EXCIT
2.1.1 Операнд RESULTAT
♦RESULTAT =resu
Имя структуры данных дополняющее результат. Этот аргумент может быть тем же, что использовался для понятия, дополненный оператором или другим именем, которое создаст новый результат структуры данных.
Примечание:
- В большинстве ситуаций, структура данных resu содержит всю необходимую информацию с вычислением опций: элементарная модель, материальную область, характеристики, нагружения. Ключевые слова MODELE , CHAM_MATER , CARA_ELEM и EXCIT при этом бесполезны.
2.1.2 Операнды MODELE / CHAM_MATER / CARA_ELEM
◊MODELE = Mo
Имя модели на которой вычисляются силы, напряжения, деформаций, и т.д.
Необязательно, потому что может быть извлечено из результата.
◊CHAM_MATER = chmater
Материальная область связанная с моделью Mo. Это ключевое слово опционально и должно предоставляться только в исключительных случаях (например добровольная поправка материала).
◊CARA_ELEM = carac
Элементарные характеристики связанные с моделью Mo, если она содержит структурные элементы или если на изопараметрические элементы влияет локальная система координат анизотропии.
Это ключевое слово опционально и может быть извлечено из результата.
2.1.3 Ключевое слово EXCIT
Данное ключевое слово (опционально) позволяет определить термические или механические нагружения чтобы использовать их для вычисления опций, вместо тех, которые были полезны в вычислении структур данных определенных под ключевым словом RESULTAT.
Определение этого ключевого слова идентично тем, командам, которые строят структуры данных resu: смотри команды MECA_STATIQUE (U4.51.01), STAT_NON_LINE (U4.51.03), DYNA_LINE_HARM (U4.53.11), и DYNA_LINE_TRAN (U4.53.02).
2.2 Выбор сеток связанных с вычислениями
Ключевые слова TOUT, GROUP_MA и MESH позволяют пользователю выбирать сетки, на которых он желает выполнять элементарные вычисления постобработки.
/TOUT = “OUI”
Все сетки (проведенные по конечным элементам) будут обработаны. Это значение по умолчанию.
| GROUP_MA=l_grma | MAILLE=l_maille
Только сетки включенные в l_grma и/или l_maille будут обработаны.
2.3 Выбор последовательных номеров
Использование ключевых слов TOUT_ORDRE, NUM_ORDRE, INST, FREQ описано в документе(U4.71.00).
2.4 Местоположение Dans
Одни только области продолжения документа точно не определяют расположение полей.
На самом деле расположение дано в имени области ( как и опции):
- Область по элементам: *_ELEM
- Область на гауссовских точках по элементам: *_ELGA
- Область с узлами по элементам: *_ELNO
- Область с узлами: *_NOEU
Области, в большинстве случаев, вычисляются на гауссовских точках (*_ELGA).
Области с узлами по элементам (*_ELNO) получаются экстраполяцией начиная с области на гауссовских точках ( детальный метод в (R3.06.03)).
Области с узлами (*_NOEU) получаются начиная с областей с узлами по элементам при помощи вычисления простого арифметического среднего (не сбалансированного размером сетки) значений записанных на элементах в данном узле.
Примечание 1:
- Для вычисления эквивалентов, области с узлами по элементам (*_ELNO) не получают экстраполяцией начиная с области на гауссовских точках. Экстраполяцию выполняют на области напряжения или растяжения и затем вычисляют область эквивалента.
Примечание 2:
- Средние с узлами вычисленных областей в локальной системе координат уместны только если углы между ними слабы. В противном случае, они не имеют смысла.
Примечание 3:
- Когда ключевые слова GROUP_MA или MESH обозначены для вычисления опции *_NOEU, арифметическое среднее берется на выбранных сетках. Таким образом, этот результат отличен от того, который получен при помощи полного вычисления.
Например:
- В одном случае структура под нагрузкой σxy аналогична:
- В полном вычислении σxy обращается в ноль на M1∩M2 как среднее двух противоположных величин. Эти величины далеки от ноля, как показывает вычисления только на M1. Величины на границах требуемых областей таким образом интерпретируются с предупреждением.
Для структурных элементов, которые имеют субточки ( мультиволоконные плиты, оболочки, балки, трубы), области типа *_ELGA и *_ELNO вычисляются на всех субточках. Для нахождения области с одной субточкой (например слой или уровень), необходимо сделать извлечение через оператор POST_CHAMP (опции EXTR_COQUE, EXTR_PMF и EXTR_TUYAU).
Более того эта промежуточная операция необходима для вычисления области типа *_NOEU type для этих структурных элементов: области типа *_NOEU в действительности никогда не имеют субточек.
Наконец, вычисление опций энергии никогда не производит области на субточках. В действительности для структурных элементов , области интегрируются по толщине ( интегрирование выполняется на субточках).
2.5 Зависимость областей
Вычисление области может потребовать предварительного вычисления одной или более других областей. Таким образом, например, для вычисления области *_NOEU необходимо иметь ту же самую область по элементам для узлов *_ELNO и гауссовских точек *_ELGA.
Зависимость решается оператором CALC_CHAMP, который автоматически проводит вычисления в промежуточной области. Таким образом, пользователю необязательно знать порядок зависимости опций.
Только области точно требуемые пользователем сохраняются в результате структуры данных.
2.6 Операнды для опций механики
2.6.1 Вычисление опций напряжений (операнд FORCED)
Компоненты областей напряжений и сил детально представлены в документе (U2.01.05).
| 'EFGE_ELGA' | 'EFGE_ELNO' | 'EFGE_NOEU'
Вычисление обобщенных сил (структурных элементов).
Его действие состоят из извлечения сил находящихся в области SIEF_ELGA/STRX_ELGA ( в случае балочных элементов/труб или дискретных), из вычисления при помощи интегрирования напряжений ( в случае мульти волоконных балочных элементов или плит и оболочек).
Примечание 1:
- Область EFGE_ELNO не всегда является экстраполяцией области EFGE_ELGA; в частности для линейного вычисления где эта область вычисляется прямо начиная с замещения. Это так потому, что некоторые компоненты не вычисляются в нелинейности.
Примечание 2:
- Для смещенных плит, силы вычисляются в плоскости сетки. При желании получить эти силы в средней плоскости плиты, команда POST_CHAMP должна использоваться / COQUE_EXCENT.
| 'SIEF_ELGA' | 'SIEF_ELNO' | 'SIEF_NOEU'
Вычисление напряженных состояний (силовых или сил обобщенных относительно моделизации) начиная со смещений (линейная теория упругости), смотри (U2.01.05).
Заметка:
- Область “SIEF_ELGA” подсчитывается при помощи нелинейных операторов разложения. Результат структуры данных всегда представляется типом evol_noli.
| 'SIGM_ELGA' | 'SIGM_ELNO' | 'SIGM_NOEU'
Вычисление напряженных состояний.
Это фактически об извлечении напряжений содержащихся в области SIEF_ELGA, смотри (U2.01.05).
| 'SIPO_ELNO' | 'SIPO_NOEU'
Вычисление напряжений в сечении балки(бруска) разбитого на составляющие обобщенной силы.
Лист компонентов области:
- Contribution of the normal force - Составляющая нормальной силы
- Contribution of the bending moment - Составляющая сгибающего момента
- Contribution of the shearing force - Составляющая сдвигающего момента
- shear coefficient in direction – коэф. сдвига в направлении
- Contribution of the twisting moment - Составляющая скручивающего момента
Силы описанные выше представлены в локальной системе координат, то есть главная база инерции поперечного сечения (R3.08.01).
Значения σxx относительно двух сгибающих моментов являются максимальными значениями вычисленными в Ymin, Ymax с одной стороны и в Zmin, Zmax с другой стороны (за исключением основного сечения, в котором пользователь предоставляет локальные экстремумы при помощи ключевых слов, RZ и RT смотри AFFE_CARA_ELEM (U4.42.01)).
Для прямоугольного сечения:
- подсчитывается значение SMFY в z=HZ /2,
- подсчитывается значение SMFZ в y=HY /2.
Для круглого сечения, подсчитывается значения SMFY и SMFZ для y и z меньше R.
| 'SIPM_ELNO'
Вычисление напряжений в максимуме и минимуме сечений балки начинают с обобщенных сил (линейная теория упругости).
То же замечание, что и для SIPO_ELNO применимо в случае основного сечения.
| 'SIRO_ELEM'
Вычисление напряжений спроецированных на поверхность объема (например на облицовки гидравлических устройств).
Список компонентов области:
- In total reference – в базисе
- Eigenvalue – собственное значение
Области извлекаются начиная с областей напряжений посчитанных на объемных сетках (MODELISATION= ` 3D' или “3D SI'):
- Идентификация объемных сеток соответствующих обрывам поверхностей сеточной группы;
- Восстановление напряжений 3D чтобы присвоить им размеры узлов;
- Среднее каждого компонента тензора напряжений в центре сторон элементов;
- Берем базис образованный нормальным вектором n на обрыве и плоскостью обрыва. Получаем ранее отмеченный тензор [σ].
- Вычисляется
который является вектором коллинеарным с
. Тогда
это вектор представляющий сдвиги которые ничтожно малы в случае граней сонаправленных/противоположнонаправленных торможению. Следует заметить, что
и SIGN показывают наличие растяжение, если они положительно и сжатия, если отрицательны.
- Перемещаемся в место предполагаемых малых сдвигов
- Ищем главные векторы относящиеся к σ2D : так получаем векторы
, которые находятся в области разрыва и собственных значений SIGT1 и SIGT2.
Примечание 1:
- В случае разрывов погруженных в объеме, выбирается объемная сетка, сторона “-” которой соотнесена с нормалью разрыва. Благодаря командам MODI_MAILLAGE / ORIE_PEAU_3D / GROUP_MA_VOLU у пользователя есть возможность переориентировать эту нормаль по своему желанию. Эта условность соответствует тому, что получается на облицовках внешних к упору, если нормаль с обрывами «вне границ».
Примечание 2:
- Если TOUT= is informed'OUI”, список сеток фильтруется , оставляя только сетки поверхности.
2.6.2 Вычисление опций деформаций (Операнд DEFORMATION)
Компоненты областей деформаций детально описаны в (U2.01.05).
| 'DEGE_ELGA' | 'DEGE_ELNO' | 'DEGE_NOEU'
Вычисление обобщенных деформаций начинают со смещений.
Эта опция имеет значение только для структурных элементов (балка, пита, труба).
Обобщенные деформации получают в локальной системе координат элемента.
| 'EPFD_ELGA' | 'EPFD_ELNO' | 'EPFD_NOEU'
Вычисление деформаций, для моделей BETON_UMLV_FP и BETON_BURGER_FP.
| 'EPFP_ELGA' | 'EPFP_ELNO' | 'EPFP_NOEU'
Вычисление основной деформации ползучести, связанной с моделью GRANGER_FP, моделью BETON_UMLV_FP или моделью BETON_BURGER_FP.
| 'EPME_ELGA' | 'EPME_ELNO' | 'EPME_NOEU'
Вычисление механических деформаций начинают со смещений.
Это вычисление выполняется в теории «малых смещений». Вычисленные напряжения равны общим прогибам минус термические деформации.
| 'EPMG_ELGA' | 'EPMG_ELNO' | 'EPMG_NOEU'
Вычисление механических деформаций начинают со смещений.
Это вычисление выполняется в теории «больших смещений». Вычисленные напряжения равны общим прогибам минус термические деформации.
| 'EPSG_ELGA' | 'EPSG_ELNO' | 'EPSG_NOEU'
Вычисление деформаций Грина-Лагранжа.
| 'EPSI_ELGA' | 'EPSI_ELNO' | 'EPSI_NOEU'
Вычисление деформаций начиная со смещений.
Для структурных элементов это деформации получают в локальной системе координат элемента.
| 'EPSP_ELGA' | 'EPSP_ELNO' | 'EPSP_NOEU'
Вычисление unelastic деформаций начиная с области смещения u, напряжения σ, температуры T, возможных unelastic деформаций εa, и внутренних переменных,
- где εfl это чистая деформация «усушки???» Гренгера.
| 'EPVC_ELGA' | 'EPVC_ELNO' | 'EPVC_NOEU'
Вычисление деформаций относящихся к командным переменным. Для времени определяются только следующие компоненты:
- термические деформации: EPTHER_L, EPTHER_T, EPTHER_N такие как:
(если материал изотропный, то три компоненты), T это температура, αi термический коэффициент расширения;
- усадка усушки EPSECH (используется для моделей описывающих поведение бетона???)
, S это командная переменная усушки и Kdessic коэффициент усадки высыхания;
- усадка гидратации EPHYDR (используется для моделей описывающих поведение бетона)
, h это командная переменная гидратации и Bendog эндогенный коэффициент усадки.
- Деформации относящиеся к жидкостному давлению (для термо-гидро-механики с разложением по последовательности): EPPTOT такие как:
, ptot это командная переменная неподвижного давления жидкости, b это коэффициент Biot, K модуль эластичности.
2.6.3 Опции извлечения внутренних переменных (операнд VARI_INTERNE)
| 'VAEX_ELGA' | 'VAEX_ELNO' | 'VAEX_NOEU'
Извлечение внутренних переменных только в THM.
Цель этой опции быть способной проводить пост-обработку внутренних переменных в THM в удобном для пользователя виде. Принципом этих областей является извлечение из области VARI_ELGA (или VARI_ELNO) ее (одну и только одну) внутреннюю переменную, которая нас интересует, через ключевое слово без надобности узнавать ее имя в области VARI_*.
Список возможных компонентов области (область имеет только один компонент, который выбирает пользователь через NOM_VARI):
DPORO изменение пористости материала DRHOLQ изменение плотности материала DPVP изменение парового давления SATLIQ объемное насыщение liquid EVP накопленная пластическая деформация IND_ETA индикатор механического состояния D Величина ущерба IND_END Индикатор ущерба TEMP_MAX максимальная температура GAMP критический девиатор пластического напряжения cumulated PCR давление SEUIL_HYD водный порог IND_HYD индикатор водной необратимости PCOHE давление сцепления COMP_ROC изотропное поведение rock SEUIL_ISO порог ANG_DEV пороговый угол девиатора X11 Компоненты тензора кинематического затвердевания X22 Компоненты тензора кинематического затвердевания X33 Компоненты тензора кинематического затвердевания X12 Компоненты тензора кинематического затвердевания X13 Компоненты тензора кинематического затвердевания X23 Компоненты тензора кинематического затвердевания DIST_DEV расстояние обобщенное с пороговым девиатором DEV_SUR_CRIT Сообщение открывающее пороговый девиатор и пороговое девиаторное критическое DIST_ISO Дистанция приведенная к изотропному порогу NB_ITER Число внутренних итераций ARRET Значение локального теста остановки итеративного процесса NB_REDE Число местных пересокращений временного шага SIGNE Знак стянутого произведения девиаторного напряжения на девиаторную пластическую деформацию
◊NOM_VARI = / nom_vari, [TXM] имя внутренней переменной.
Примечание 1:
- Когда извлекаемая переменная не формирует часть внутренних переменных модели, появляется предупреждение, но область ведет себя также как и с R8VIDE () (очень большое число около 1.0E+308).
Примечание 2:
- Область VAEX_NOEU подсчитыватся начиная с VAEX_ELNO и не извлечением области VARI_NOEU.
| 'VARC_ELGA' Вычисление командных переменных используется для механических вычислений.
Список компонентов области: TEMP HYDR SECH CORR Смотри документацию IRRA команды AFFE_MATERIAU (U4.43.03) PTOT для определения каждого компонента. DIVU NEUT1 NEUT2
Заметка:
- Все компоненты систематически вычисляются. Переменные, которые не определены инициализируются со значением R8VIDE () (очень большое число около 1.0E+308).
| 'VARI_ELNO' | 'VARI_NOEU' Вычисление внутренних переменных.
Список компонентов области: V1 внутренняя переменная 1 … VI внутренняя переменная I … Vn внутренняя переменная N
Число и тип этих внутренних переменных специфичен для каждой модели поведения (см (U4.51.11)).
Заметка:
- Область “VARI_ELGA” подсчитывается нелинейными операторами разложения. Результат структуры данных всегда представляется типом evol_noli.
2.6.4 Вычисление опций энергии (операнд ENERGIE)
| 'DISS_ELEM'
Вычисление энергии выделяемой при разрушении. Полученная область имеет только один компонент по имени “ENDO”.
Список компонентов области: ENDO энергия выделяемая при разрушении
Примечание:
- Действительно только для элементов DKTG и модели GLRC_DM. Их формулировка дана в (R7.01.32).
| 'DISS_ELGA' | 'DISS_ELNO' | 'DISS_NOEU'
Вычисление плотности энергии выделяемой при разрушении. Полученная область имеет только один компонент по имени “ENDO”.
Список компонентов области: ENDO энергия выделяемая при разрушении
Примечание:
- Действительно только для элементов DKTG и модели GLRC_DM. Их формулировка дана в (R7.01.32).
| 'ECIN_ELEM'
Вычисление кинетической энергии.
Список компонентов области: --------------------------------------------------- TOTALE кинетическая энергия --------------------------------------------------- Дополнительные компоненты для плит и оболочек: MEMBRANE составляющие кинетической энергии (см (R3.07.03)) Contributions --------------------------------------------------- Дополнительные компоненты для скрученных балок: PLAN_XY составляющие кинетической энергии (см (R3.08.01)) PLAN_XZ --------------------------------------------------- Дополнительные компоненты для дискретных элементов: DX DY DZ составляющие кинетической энергии DRX DRY DRZ
| 'ENEL_ELEM'
Вычисление энергии упругой деформации.
Список компонентов области: ---------------------------------------------- TOTALE энергия упругой деформации ---------------------------------------------- Дополнительные компоненты для плит и оболочек: MEMBRANE Contributions BENDING with elastic strain energy (см (R3.07.03))
| 'ENEL_ELGA' | 'ENEL_ELNO' | 'ENEL_NOEU'
Вычисление плотности энергии упругой деформации.
Список компонентов области: ---------------------------------------------- TOTALE энергия упругой деформации ---------------------------------------------- Дополнительные компоненты для плит и оболочек: MEMBRANE Contributions BENDING with elastic strain energy (см (R3.07.03))
| 'EPOT_ELEM'
Вычисление потенциальной энергии деформации, начиная с области смещения u и областей температуры T:
Список компонентов области: ---------------------------------------------- TOTALE потенциальная энергия ---------------------------------------------- Дополнительные компоненты для плит и оболочек: MEMBRANE Contributions BENDING with potential energy (см (R3.07.03)) ---------------------------------------------- Дополнительные компоненты для прямых балок: TRAC_COM TORSION Составляющие потенциальной FLEX_Y энергии (см (R3.08.01)) FLEX_Z ---------------------------------------------- Дополнительные компоненты для скрученных балок: PLAN_XY Составляющие потенциальной PLAN_XZ энергии (см (R3.08.01)) ---------------------------------------------- Дополнительные компоненты для дискретных элементов: DX DY DZ Составляющие потенциальной энергии DRX DRY DRZ
- для элементов непрерывных сред 2D и 3D:
- для балочных элементов:
- для оболочных элементов и оболочек:
| 'ETOT_ELEM '
Вычисление прироста полной энергии деформации запускает течение времени и предыдущее время.
Список компонентов области: TOTALE прирост полной энергии деформации
| 'ETOT_ELGA” | 'ETOT_ELNO' | 'ETOT _NOEU'
Вычисление прироста плотности полной энергии деформации запускает течение времени и предыдущее время.
Список компонентов области: TOTALE прирост полной энергии деформации
2.6.5 Вычисление опций критериев (операнд CRITERES)
| 'DERA_ELGA' | 'DERA_ELNO' | 'DERA_NOEU'
Вычисление местного индикатора пропуска и индикатора потерь радиальности (R4.20.01).
Список компонентов области: ----------------------------------------------------------------------------- DCHA_V Индикатор разгрузки рассчитанный на локальном тензоре напряжений ----------------------------------------------------------------------------- DCHA_T Индикатор разгрузки рассчитанный на общем тензоре напряжений ----------------------------------------------------------------------------- IND_DCHA Индикатор, позволяющий узнать останется ли разгрузка эластичной или есть риск пластификации при использовании чистого кинематического затвердевания ----------------------------------------------------------------------------- VAL_DCHA Показывает пропорцию выхода критерия в случае жесткой разгрузки ----------------------------------------------------------------------------- X11 X22 Компоненты кинематического тензора используемого X33 для вычисления IND_DCHA X12 X13 X23 ----------------------------------------------------------------------------- RADI_V Индикатор изменения направления напряжений между моментом времени t и t+Δt ----------------------------------------------------------------------------- ERR_RADI Ошибка η вызванная дискретизацией во времени, напрямую связанная с вращением нормали на поверхности нагрузки -----------------------------------------------------------------------------
DCHA_V и DCHA_T показывают, в обоих случаях, относительное изменение нормали напряжений в рамках значений Von Mises:
ND_DCHA может принимать следующие значения:
- 0 : свободное начальное значение;
- 1 : при эластичной нагрузке;
- 2 : при пластичной нагрузке;
- -1 : при незначительной эластичной разгрузке (независимо от типа затвердевания);
- -2 : при значительной разгрузке (может произойти пластификация с кинематическим затвердеванием).
RADI_V задано следующим соотношением:
Количество равно нулю, когда радиальность сохраняется с течением времени.
ERR_RADI это угол включающий n−, нормаль с критерием пластичности в начале временного шага (мгновенное t−), и n+, нормаль с критерием пластичности подсчитанная в конце временного шага (мгновенное t+) следующим образом:
Это позволяет измерить ошибку (также используется для уточнения временного шага (U4.51.11)). Это рабочий критерий для эластопластичного характера изменения Von Mises с изотропным затвердеванием, линейно кинематический и смешанный VMIS_ISOT_LINE,VMIS_ISOT_TRAC, VMIS_ISOT_PUIS, VMIS_CINE_LINE, VMIS_ECMI_LINE, VMIS_ECMI_TRAC, и для упруго-вязко-пластичного характера изменения Chaboche:VMIS_CIN1_CHAB, VMIS_CIN2_CHAB, VMIS_CIN2_MEMO, VISC_CIN1_CHAB, VISC_CIN2_CHAB, VISC_CIN2_MEMO.
Примечание:
- Вычисление этих опций требует сравнения областей напряженности в моменты времени ti и ti+1. Результат собирается в последовательность номеров, связанных с временем ti. Индикатор разгрузки вычисляется таким образом:
.
- По умолчанию вычисление производится для неморов последовательности от 1 до n-1. Но если определять список времен (с наличием «дыр»), вычисление будет связано только с требуемыми временами, но всегда будет производится сравнение между временем ti и ti+1 в списке времен используемых для нелинейных вычислений.
| 'ENDO_ELGA' | 'ENDO_ELNO' | 'ENDO_NOEU'
Вычисление разрушения d начинают с тензора напряжений и накопленной пластической деформации p.
Список компонентов области: ----------------------------------------------------------------------- TRIAX Показатель трехосности ----------------------------------------------------------------------- SI_ENDO Эквивалентной напряжение разрушения Lemaître-Sermage ----------------------------------------------------------------------- COENDO Стандартизированное ограничение напряжения Lemaître-Sermage ----------------------------------------------------------------------- DOM_LEM Разрушение Lemaître-Sermage -----------------------------------------------------------------------
Показатель трехосности α задан следующим выражением:
И эквивалентное напряжение разрушения σ*:
Кинетика разрушения задается моделью Lemaître-Sermage:
- где S и s это коэффициенты характеристик материала и pseuil это порог разрушения относительно энергии запасенной в материале (если s=1 получаем обычную модель Lemaître).
| 'EPEQ_ELGA' | 'EPEQ_ELNO' | 'EPEQ_NOEU' | 'EPMQ_ELGA' | 'EPMQ_ELNO' | 'EPMQ_NOEU'
Вычисление подсчитываемых эквивалентов деформаций начинают с областей EPSI_*, или EPME_*).
Список компонентов области: ------------------------------------------------ INVA_2 эквивалентная деформация Von Mises ------------------------------------------------ PRIN_1 Основные деформации, PRIN_2 упорядоченные в порядке PRIN_3 возрастания ------------------------------------------------ INVA_2SG эквивалентная деформация Von Mises со знаком следа ε ------------------------------------------------ VECT_1_X VECT_1_Y VECT_1_Z VECT_2_X Основные направления VECT_2_Y VECT_2_Z VECT_3_X VECT_3_Y VECT_3_Z ------------------------------------------------
Эквивалентные деформаций Von Mises заданы следующим выражением:
Примечание:
- Стоит отметить, что эквивалентные деформации полученные начиная с EPSI_* и EPME_* идентичны. На самом деле, различием между двумя тензорами является сферический тензор (термическая деформация). В то время как эквивалентную деформацию получают начиная с второго инварианта девиатора, сферический тензор «исчезает» когда берется девиатор.
| 'INDL_ELGA'
Вычисление индикатора местоположения, основанного на акустическом тензоре (критерий RICE).
Список компонентов области: Индикатор местоположения INDICE 0 если det(N.H.N) > 0 1 если нет, что соответствует основанию нормали DIR1 Первая зона нормали DIR2 Вторая зона нормали DIR3 Третья зона нормали DIR4 Четвертая зона индикации
Такое местоположение определяет состояние, в котором локальная проблема интеграции характера изменения теряет свою уникальность. Это определяется как: det(N.H.N) ≤ 0 , где H назначает тангенциальный оператор, а N нормаль с направлениями местоположения.
Примечание:
- Этот метод работает только в случае 2D и для составных законов типа DRUCK_PRAGER и HUJEUX.
| 'PDIL_ELGA'
Вычисление модуля жесткости микрокомпьютерного расширения.
Список компонентов области: A1_LC2 модуль жесткости микрокомпьютерного расширения
Опция PDIL_ELGA создает в ?? сред второго градиента расширения значение модуля A1_LC2, позволяя контролировать интервал нетривиального решения первоначально однородной проблемы (R5.04.03).
Вычисление A1_LC2 выполняется через оценку функции зависящей от геометрического направленного значения подразумеваемого материала. Текущая назначенная угловая дискретизация равна 0.1°.
Примечание:
- Этот метод работает только для составных законов типа DRUCK_PRAGER и HUJEUX.
| 'SIEQ_ELGA' | 'SIEQ_ELNO' | 'SIEQ_NOEU'
Вычисление эквивалентных напряжений начиная с областей напряжений.
Список компонентов области: ---------------------------------------------- VMIS эквивалентное напряжение Von Mises ---------------------------------------------- TRESCA сила Tresca ---------------------------------------------- PRIN_1 Основные напряжения, PRIN_2 упорядоченные в порядке возрастания PRIN_3 ---------------------------------------------- VMIS_SG эквивалентное напряжение Von Mises со знаком следа σ ---------------------------------------------- VECT_1_X VECT_1_Y VECT_1_Z VECT_2_X VECT_2_Y Основные направления VECT_2_Z VECT_3_X VECT_3_Y VECT_3_Z ---------------------------------------------- TRSIG След (черта, отпечаток) σ ---------------------------------------------- TRIAX Показатель трехосности ----------------------------------------------
Эквивалентное напряжение Von Mises задается следующим выражением:
Уровень трехосности задается следующим выражением:
2.6.6 Вычисление опции гидравлических течений (операнд HYDRAULIC)
| 'FLHN_ELGA'
Вычисление опции гидравлических течений в THM Фij = Mij∙ν на гранях элементов (2D или 3D) начинают с вектора течения с узлами. Mij это гидравлический вектор течения компоненты ij.
Список компонентов области: FH11 FH22 FH12 FH21
2.7 Операнды термических опций
2.7.1 Операнд THERMAL
| 'FLUX_ELGA' | 'FLUX_ELNO' | 'FLUX_NOEU'
Вычисление потока тепла начинают с температуры.
Список компонентов области: ------------------------------------------------ FLUX Потоки тепла в трех направлениях FLUY пространства (в среднем слое оболочек) FLUZ ------------------------------------------------ Дополнительные компоненты для оболочек: ------------------------------------------------ FLUX_INF Потоки тепла в трех направлениях FLUY_INF пространства в нижнем слое FLUZ_INF ------------------------------------------------ FLUX_INF Потоки тепла в трех направлениях FLUY_INF пространства в верхнем слое FLUZ_INF ------------------------------------------------
| 'HYDR_NOEU'
Вычисление гидратации.
Список компонентов области: HYDR Гидратация
Примечание:
- Область “HYDR_ELNO” вычисляется первоначально нелинейным термическим оператором THER_NON_LINE для конкретной моделизации (R7.01.12).
| 'SOUR_ELGA'
Вычисление источника тепла.
Список компонентов области: SOUR источник тепла
Этот источник рассчитывается начиная с электрического потенциала через модель Ohm. Этот электрический потенциал должен вычисляться оператором THER_LINEAIRE (U4.54.01) при помощи нахождения необходимых аналогий. Этот источник может затем быть использован в термических вычислениях через ключевое слово SOUR_CALCULEE команды AFFE_CHAR_THER (U4.44.02).
| 'ETHE_ELEM'
Вычисление термической энергии до равновесия начинают с области температуры T .
Список компонентов области: Thermal Полная энергия
2.8 Операнды для акустических опций
2.8.1 Операнд ACOUSTIC
| 'INTE_ELNO' | 'INTE_NOEU'
Вычисление акустической интенсивности. Определения в (R4.02.01).
Список компонентов области:
INTX_R | Акустическая интенсивность, действительная часть относительно оси X |
INTY_R | Акустическая интенсивность, действительная часть относительно акустической оси Y |
INTZ_R | Акустическая интенсивность, действительная часть относительно акустической оси Z |
INTX_I | Акустическая интенсивность, мнимая часть относительно оси X |
INTY_I | Акустическая интенсивность, мнимая часть относительно акустической оси Y |
INTZ_I | Акустическая интенсивность, мнимая часть относительно оси Z |
| 'PRAC_ELNO' | 'PRAC_NOEU'
Вычисление давления на узлы частично действительное, мнимая часть и децибелы.
Список компонентов области: PRES_I Акустическое давление, действительная часть PRES_R Акустическое давление, мнимая часть DB Акустическое давление в децибелах
| 'PRME_ELNO'
Вычисление давления на узлы для элементов FLUID.
Список компонентов области: DB Акустическое давление в децибелах
2.9 Операнд для узловых сил и реакций
2.9.1 Операнд FORCE
| 'FORC_NODA'
Вычисление узловых обобщенных сил начинают с напряжений в Gauss точках.
Список компонентов области: ---------------------------------------------------- DX DY Узловые силы DZ ---------------------------------------------------- Дополнительные компоненты для структурных элементов: DRX DRY Узловые силы DRZ ----------------------------------------------------
Вычисление производится следующим образом:
- где σK напряжениями в Gauss точках элемента K;
- B оператор конечных элементов обобщенных деформаций;
- uK обобщенное элементарное смещение
- где B это матрица соединяющая деформации первого порядка с смещениями.
Размерность узловых сил двойственно для uK чтобы дать работу (в Джоулях).
Для балочных элементов и дискретных элементов, напряжения в Gauss точках фактически являются узловыми силами обобщенными по отношению к элементу (находится при помощи произведения жесткости матрицы элемента, смещения и подсчетом сил термического происхождения и распределенных сил. Вычисление узловых сил выполняется при помощи проецирования узловых сил содержащихся в области с именем “SIEF_ELGA” в общем отношении. Затем суммируют элементы указанные выше. Компоненты DX, DY и DZ дают силы и DRX, DRY DRZ моменты.
Для осесимметричных элементов, интеграция в theta выполняется в секторе в 1 радиан. При желании вычислить интеграл поверхностной силы на всем диске необходимо умножить на 2π. Для элементов плоских деформаций, вычисление выполняется на ленте на единицу ширины. Вычисленные узловые силы в таком случае как силы на единицу длины. При желании подсчитать узловые силы оказываемые на структуру шириной l, необходимо умножить результат в D_PLAN на l, при таком близком предположении плоская деформация отлична для двух сторон. Поэтому получаем приблизительный результат.
Для твердых элементов (3D, 2D and брусья), FORC_NODA в общем имеет Размерность силы. Именно в области узлов сетки значение в узле получается начиная с напряжений вычисленных в конвергентных элементах с этим узлом, поэтому их значения меняются когда меняется сетка. В отсутствие распределенных нагрузок, равновесие налагает их недействительность на внутренний узел, в то время как они отвечают на реакцию на опору где налагается кинематическое отношение (случай наложенного смещения).
В случае оболочек, компоненты DX, DY и DZ дают FORC_NODA (Размерность силы) в общем отношении к сетке. Эти компоненты построены с нормальной силой и режущей кромкой оболочки. Компоненты DRX, DRY и DRZ дают FORC_NODA (размерность одного момента) в общем отношении к сетке, построенной со сгибающими моментами оболочки.
В гидравлике, обобщенные узловые силы связанные с каждым компонентом относящемуся к потоку. Если отметить QTσ0 результатом FORC_NODA, для гидравлических уравнений, тогда для временного шага Δt, то:
В FORC_NODA :
- со степенью свободы PRE1 связан поток воды
- со степенью свободы PRE2 связан поток газа
- со степенью свободы TEMP связан поток тепла
с q поток тепла и Mw, Mvp , Mas и Mad гидравлические потоки жидкой воды, испарений, сухого воздуха и воздуха растворенного в жидкости. Эти данные относятся к обобщенным напряжениям M11 ,M12 ,M21 ,M22.
| 'REAC_NODA'
Вычисление узловых сил обобщенных реакций с узлами, напряжений в Gauss точках.
Список компонентов области: --------------------------------------------------- DX DY узловые силы DZ --------------------------------------------------- Дополнительные компоненты для структурных элементов: DRX DRY узловые силы DRZ ---------------------------------------------------
Для концепций результата типа evol_elas, mult_elas, fourier_elas или evol_noli, это вычисление выполняется так:
- где f - это объемные силы
- Fs - обобщенные поверхностные силы
- Fi - специфические силы с узлом i
Если отметить RK вектор узловых реакций элемента K, начинает с обобщенных узловых сил:
Другими словами отсекаем узловые силы применяя внешние силы к элементу K .
Стоит отметить, что нагрузочная температура не появляется во внешних силах.
В динамике, для выполнения узловых реакций, советуют дополнительно убрать силы инерции (ускорения) и демпфирования (скорости). На данный момент в Code_Aster эффектами демпфирования на узловых реакциях пренебрегают.
Для концепций результата типа mode_meca (полученных при помощи модальных вычислений) формула такова:
- где M это матрица массы
- ω собственная пульсация
- u область смещения
Для концепций результата типа dyna_trans при помощи динамических вычислений линейного переходного процесса (DYNA_LINE_TRAN, или DYNA_TRAN_MODAL при помощи REST_GENE_PHYS), типа dyna_harmo полученных при помощи гармонических вычислений (DYNA_LINE_HARM) или типа evol_noli полученных при помощи вычислений нелинейной динамики переходного процесса (DYNA_NON_LINE) выражение таково:
- где M матрица массы;
- ü область ускорения;
- L вектор применения внешних сил
Примечание 1:
- Узловые реакции нулевые в любой внутренней точке модели и не нулевые в точках границы приведенной в граничное кинематическое состояние (или с силами контакта).
- Однако, факт пренебрежения вкладом демпфирования в динамику может создать легкую вариацию сравнимую с точным результатом.
Примечание 2:
- Если обозначено ключевое слово GROUP_MA, опции “FORC_NODA” и “REAC_NODA” вычисляются следующим образом:
FKвычисляется только если запрошенные элементы потом собираются. Результат отличен от общего вычисления на всей области затем уменьшается до требуемых элементов. - Признанный метод позволяет измерить реакцию одной части модели на другую.
2.9.2 Пример 1: Структура нагруженная узловой силой (2 элемента QUAD4)
Согласно этому примеру, реакции на узлах (2) почти равны узловым силам (1) минус нагрузка. Они показывают реакции на опору структуры.
Если ограничено вычисление с сеткой M1, силы (3) с узлами принадлежащими к границе вводят M1 и M2 по разному. Они показывают реакцию модели сформированной M1 на модель сформированную M2. Стоит отметить, что узловые нагрузки поделены пополам, потому что две сетки вносят свой вклад в него. Узловые реакции (4) также равны узловым силам минус нагрузка.
На вычислениях ограниченных с сеткой M2, узловые сил (5) относительно OX противоположного знака к вычислениям ограниченным с сеткой M1, иллюстрируя принцип действия и противодействия.
2.9.3 Пример 2: Структура с нагружаемой температурой
В этом примере, узловые силы и узловые реакции совпадают, потому что единственной нагрузкой является нагружаемая температура.
Если ограниченно вычисление с сеткой M2, силы относительно OY остаются теми же, но меняются относительно OX.
2.10 Вычисления в области пользователя
Ключевое слово CHAM_UTIL позволяет вычислить области как неопределенные, известные как «пользовательские», потому что имя будет влиять на концепцию их результата.
Может быть несколько вхождений в CHAM_UTIL с целью подключения вычисления нескольких областей.
Обработка проводится в конце команды CALC_CHAMP, здесь доступны вычисленные области с предшествующими ключевыми словами (STRESS, DEFORMATION…).
Также запрашивается вычисление заданного критерия, или применение одной или более формул для вычисления других областей.
2.10.1 Операнд NOM_CHAM
Действует область, в которой выполняются вычисления. Полученная область будет того же типа: ELGA, ELNO или NOEU.
2.10.2 Операнд CRITERE
Запрашивает вычисление заданного критерия. Критериями являются (параграф 2.6.5 предоставляет определения каждого критерия):
- VMIS (для областей напряжения),
- INVA_2 (для областей деформаций),
- TRACE (для областей напряжения или деформаций).
Каждый из этих критериев создает компонент (называемый X1). Одним из интересов является возможность вычислить INVA_2 в каждой области деформации.
2.10.3 Операнд FORMULATES
Позволяет подсчитать любое выражение функции компонентов области предоставленной для NOM_CHAM.
Полученная область будет содержать столь же много компонентов сколько в приведенных формулах: с первой формулой совпадает компонент X1, со второй X2, и т.д. Вплоть до 30 компонентов могут быть созданы таким образом. Примером формул позволяющих найти критерий VMIS и INVA_2 могут быть найдены во второй части теста sslv104a.
2.10.4 Операнд NUME_CHAM_RESU
Полученная область должна быть упорядочена, одним образом, в концепции результата. Области «пользователь» при этом нумеруются при помощи NUME_CHAM_RESU и типа области.
Имя области при этом будет типа UT01_ELGA, UT22_NOEU, и т.д.
2.10.5 Пример вычислений в пользовательской области
Результат области UT02_ELGA с двумя компонентами. X1 это след SIGM_ELGA (сравнимо с компонентом TRSIG of SIEQ_ELGA) и X2 это эквивалентное напряжение Von Mises (компонент VMIS из SIEQ_ELGA).
fTrace = FORMULA (NOM_PARA = (“SIXX”, “SIYY”, “SIZZ”), VALE = “““SIXX + SIYY + SIZZ”””) fVonMis = FORMULA (NOM_PARA = (“SIXX”, “SIYY”, “SIZZ”, “SIXY”, “SIXZ”, “SIYZ”), VALE = “““ sqrt (3. /2. * (SIXX ** 2 + SIYY ** 2 + SIZZ ** 2 + 2*SIXY ** 2 + 2*SIXZ ** 2 +2*SIYZ ** 2) - 1. /2. * fTrace (SIXX, SIYY, SIZZ) ** 2) ”””) LMBO = CALC_CHAMP (reuse=RES, RESULT=RES, CHAM_UTIL=_F (NOM_CHAM=`SIGM_ELGA`, FORMULE= (fTrace, fVonMis), NUME_CHAM_RESU=2, ) , )
2.11 Операнд TITRE
◊TITRE = Titre
заголовок, который желаете дать результату выполнения команды (U4.02.01).