(U4.81.04) Оператор CALC CHAMP

Материал из ru.wiki.laduga.ru
Перейти к: навигация, поиск


Содержание

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 аналогична:
U4.81.04 2.4.png
В полном вычислении σ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'

Вычисление напряжений в сечении балки(бруска) разбитого на составляющие обобщенной силы.

Лист компонентов области:

U4.81.04 2.6.1 1.png
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'

Вычисление напряжений спроецированных на поверхность объема (например на облицовки гидравлических устройств).

Список компонентов области:

U4.81.04 2.6.1 2.png
In total reference – в базисе
Eigenvalue – собственное значение

Области извлекаются начиная с областей напряжений посчитанных на объемных сетках (MODELISATION= ` 3D' или “3D SI'):

  • Идентификация объемных сеток соответствующих обрывам поверхностей сеточной группы;
  • Восстановление напряжений 3D чтобы присвоить им размеры узлов;
  • Среднее каждого компонента тензора напряжений в центре сторон элементов;
  • Берем базис образованный нормальным вектором n на обрыве и плоскостью обрыва. Получаем ранее отмеченный тензор [σ].
  • Вычисляется U4.81.04 2.6.1 3.png который является вектором коллинеарным с U4.81.04 2.6.1 4.png. Тогда U4.81.04 2.6.1 5.png это вектор представляющий сдвиги которые ничтожно малы в случае граней сонаправленных/противоположнонаправленных торможению. Следует заметить, что U4.81.04 2.6.1 6.png и SIGN показывают наличие растяжение, если они положительно и сжатия, если отрицательны.
  • Перемещаемся в место предполагаемых малых сдвигов U4.81.04 2.6.1 7.png
Ищем главные векторы относящиеся к σ2D : так получаем векторы U4.81.04 2.6.1 8.png, которые находятся в области разрыва и собственных значений 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' 

Вычисление механических деформаций начинают со смещений.

Это вычисление выполняется в теории «малых смещений». Вычисленные напряжения равны общим прогибам минус термические деформации.

U4.81.04 2.6.2 1.png
 | 'EPMG_ELGA'
 | 'EPMG_ELNO'
 | 'EPMG_NOEU'

Вычисление механических деформаций начинают со смещений.

Это вычисление выполняется в теории «больших смещений». Вычисленные напряжения равны общим прогибам минус термические деформации.

U4.81.04 2.6.2 2.png
 | 'EPSG_ELGA'
 | 'EPSG_ELNO'
 | 'EPSG_NOEU'

Вычисление деформаций Грина-Лагранжа.

U4.81.04 2.6.2 3.png
 | 'EPSI_ELGA'
 | 'EPSI_ELNO'
 | 'EPSI_NOEU'

Вычисление деформаций начиная со смещений.

U4.81.04 2.6.2 4.png

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

 | 'EPSP_ELGA'
 | 'EPSP_ELNO'
 | 'EPSP_NOEU'

Вычисление unelastic деформаций начиная с области смещения u, напряжения σ, температуры T, возможных unelastic деформаций εa, и внутренних переменных,

U4.81.04 2.6.2 5.png
где εfl это чистая деформация «усушки???» Гренгера.
 | 'EPVC_ELGA'
 | 'EPVC_ELNO'
 | 'EPVC_NOEU'

Вычисление деформаций относящихся к командным переменным. Для времени определяются только следующие компоненты:

  • термические деформации: EPTHER_L, EPTHER_T, EPTHER_N такие как:
U4.81.04 2.6.2 6.png (если материал изотропный, то три компоненты), T это температура, αi термический коэффициент расширения;
  • усадка усушки EPSECH (используется для моделей описывающих поведение бетона???) U4.81.04 2.6.2 7.png, S это командная переменная усушки и Kdessic коэффициент усадки высыхания;
  • усадка гидратации EPHYDR (используется для моделей описывающих поведение бетона) U4.81.04 2.6.2 8.png, h это командная переменная гидратации и Bendog эндогенный коэффициент усадки.
  • Деформации относящиеся к жидкостному давлению (для термо-гидро-механики с разложением по последовательности): EPPTOT такие как: U4.81.04 2.6.2 9.png, 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'

Вычисление кинетической энергии.

U4.81.04 2.6.4 1.png
 Список компонентов области:
 ---------------------------------------------------
 TOTALE           кинетическая энергия 
 ---------------------------------------------------
 Дополнительные компоненты для плит и оболочек:
 MEMBRANE         составляющие кинетической энергии (см (R3.07.03))
 Contributions 
 ---------------------------------------------------
 Дополнительные компоненты для скрученных балок:
 PLAN_XY          составляющие кинетической энергии (см (R3.08.01)) 
 PLAN_XZ 
 ---------------------------------------------------
 Дополнительные компоненты для дискретных элементов:
 DX
 DY
 DZ               составляющие кинетической энергии 
 DRX
 DRY
 DRZ
 | 'ENEL_ELEM'

Вычисление энергии упругой деформации.

U4.81.04 2.6.4 2.png
 Список компонентов области:
   ----------------------------------------------
   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:
U4.81.04 2.6.4 3.png
  • для балочных элементов:
U4.81.04 2.6.4 4.png
  • для оболочных элементов и оболочек:
U4.81.04 2.6.4 5.png
 |   '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:

U4.81.04 2.6.5 1.png , нормаль является функцией поведения (изотропное затвердевание или линейная кинематика)

ND_DCHA может принимать следующие значения:

  • 0 : свободное начальное значение;
  • 1 : при эластичной нагрузке;
  • 2 : при пластичной нагрузке;
  • -1 : при незначительной эластичной разгрузке (независимо от типа затвердевания);
  • -2 : при значительной разгрузке (может произойти пластификация с кинематическим затвердеванием).

RADI_V задано следующим соотношением:

U4.81.04 2.6.5 2.png

Количество равно нулю, когда радиальность сохраняется с течением времени.

ERR_RADI это угол включающий n, нормаль с критерием пластичности в начале временного шага (мгновенное t), и n+, нормаль с критерием пластичности подсчитанная в конце временного шага (мгновенное t+) следующим образом:

U4.81.04 2.6.5 3.png

Это позволяет измерить ошибку (также используется для уточнения временного шага (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. Индикатор разгрузки вычисляется таким образом: U4.81.04 2.6.5 4.png.
По умолчанию вычисление производится для неморов последовательности от 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
   -----------------------------------------------------------------------

Показатель трехосности α задан следующим выражением:

U4.81.04 2.6.5 5.png

И эквивалентное напряжение разрушения σ*:

U4.81.04 2.6.5 6.png

Кинетика разрушения задается моделью Lemaître-Sermage:

U4.81.04 2.6.5 7.png
где 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 заданы следующим выражением:

U4.81.04 2.6.5 8.png

Примечание:

Стоит отметить, что эквивалентные деформации полученные начиная с 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 задается следующим выражением:

U4.81.04 2.6.5 9.png

Уровень трехосности задается следующим выражением:

U4.81.04 2.6.5 10.png

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 
 ----------------------------------------------------

Вычисление производится следующим образом:

U4.81.04 2.9.1 1.png
где σK напряжениями в Gauss точках элемента K;
B оператор конечных элементов обобщенных деформаций;
uK обобщенное элементарное смещение
U4.81.04 2.9.1 2.png обобщенные узловые силы
где 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, то:

U4.81.04 2.9.1 3.png

В FORC_NODA :

  • со степенью свободы PRE1 связан поток воды U4.81.04 2.9.1 4.png
  • со степенью свободы PRE2 связан поток газа U4.81.04 2.9.1 5.png
  • со степенью свободы TEMP связан поток тепла U4.81.04 2.9.1 6.png

с 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, это вычисление выполняется так:

U4.81.04 2.9.1 7.png
где f - это объемные силы
Fs - обобщенные поверхностные силы
Fi - специфические силы с узлом i

Если отметить RK вектор узловых реакций элемента K, начинает с обобщенных узловых сил:

U4.81.04 2.9.1 8.png

Другими словами отсекаем узловые силы применяя внешние силы к элементу K .

Стоит отметить, что нагрузочная температура не появляется во внешних силах.

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

Для концепций результата типа mode_meca (полученных при помощи модальных вычислений) формула такова:

U4.81.04 2.9.1 9.png
где M это матрица массы
ω собственная пульсация
u область смещения

Для концепций результата типа dyna_trans при помощи динамических вычислений линейного переходного процесса (DYNA_LINE_TRAN, или DYNA_TRAN_MODAL при помощи REST_GENE_PHYS), типа dyna_harmo полученных при помощи гармонических вычислений (DYNA_LINE_HARM) или типа evol_noli полученных при помощи вычислений нелинейной динамики переходного процесса (DYNA_NON_LINE) выражение таково:

U4.81.04 2.9.1 10.png
где M матрица массы;
ü область ускорения;
L вектор применения внешних сил

Примечание 1:

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

Примечание 2:

Если обозначено ключевое слово GROUP_MA, опции “FORC_NODA” и “REAC_NODA” вычисляются следующим образом:
FKвычисляется только если запрошенные элементы потом собираются. Результат отличен от общего вычисления на всей области затем уменьшается до требуемых элементов.
Признанный метод позволяет измерить реакцию одной части модели на другую.

2.9.2 Пример 1: Структура нагруженная узловой силой (2 элемента QUAD4)

U4.81.04 2.9.2 1.png
U4.81.04 2.9.2 2.png

Согласно этому примеру, реакции на узлах (2) почти равны узловым силам (1) минус нагрузка. Они показывают реакции на опору структуры.

Если ограничено вычисление с сеткой M1, силы (3) с узлами принадлежащими к границе вводят M1 и M2 по разному. Они показывают реакцию модели сформированной M1 на модель сформированную M2. Стоит отметить, что узловые нагрузки поделены пополам, потому что две сетки вносят свой вклад в него. Узловые реакции (4) также равны узловым силам минус нагрузка.

На вычислениях ограниченных с сеткой M2, узловые сил (5) относительно OX противоположного знака к вычислениям ограниченным с сеткой M1, иллюстрируя принцип действия и противодействия.

2.9.3 Пример 2: Структура с нагружаемой температурой

U4.81.04 2.9.3.png

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

Если ограниченно вычисление с сеткой 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).