(U4.53.01) Оператор DYNA NON LINE

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


Содержание

1 Краткое описание

Вычисляет динамические изменения структуры, которая имеет форму или материальный состав с нелинейным поведением. Это могут быть, например, структуры с нелинейностью материала (пластичность или большие геометрические смещения) (R5.05.05). Синтаксис этой команды подобен синтаксису оператора STAT_NON_LINE (U4.51.03).

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

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

Продукт понятие типа evol_noli.



2 Синтаксис

 dynanl [evol_noli] = DYNA_NON_LINE
 ( ◊ reuse = dynanl , [evol_noli]    
   ♦ MODELE = mo , [model]     
   ♦ CHAM_MATER = chmat , [cham_mater]    
   ◊ MODE_STAT = modestat , [mode_meca]    
   ◊ CARA_ELEM = carac , [cara_elem]     
   ◊ MASS_DIAG = / “OUI” , 
                 / “NON” ,
   ◊ EXCIT = _F ( ◊TYPE_CHARGE = / “FIXE_CSTE” , [DEFECT]
                                 / “SUIV” ,
                                 / “DIDI” ,       
                  ◊CHARGE = chi , [char_meca]
                                  [char_cine_meca]       
                  ◊ / FONC_MULT = fi , [function]
                    / DEPL = depl , [function]       
                      VITE = vite , [function]        
                      ACCE = acce , [function]      
                  ◊ MULT_APPUI = / “NON” , [DEFECT]
                                 / “OUI”,       
                  ◊ DIRECTION = (d1, d2, d3) , [l_R]      
                  ◊ NOEUD = lno , [l_noeud]       
                  ◊ GROUP_NO = lgrno , [l_gr_noeud]      
                ) ,    
   ◊ SOUS_STRUC = _F ( смотрите документ [U4.51.03] ) ,    
   ◊ AMOR_RAYL_RIGI = / “TANGENTE” , [DEFECT]                        
                      / “ELASTIQUE”,         
   ◊ AMOR_MODAL = _F (
                       ♦ MODE_MECA = mode , [mode_meca]       
                       ♦ / AMOR_REDUIT = l_amor , [l_R]         
                         / LIST_AMOR = lisamor [listr8]       
                       ◊ NB_MODE = / nbmode , [I]                  
                                   / 9999 , [DEFAUT]       
                       ◊ REAC_VITE = / “OUI” , [DEFECT]                    
                                     / “NON” ,
                     ) ,    
   ♦ | COMP_INCR = _F ( смотрите документ [U4.51.11] ) ,     
     | COMP_ELAS = _F ( смотрите документ [U4.51.11] ) ,    
   ◊ ETAT_INIT = _F ( смотрите документ [U4.51.03]       
                    | VITE = quickly , [cham_no_DEPL_R]       
                    | ACCE = acce , [cham_no_DEPL_R]       
                  ) ,    
   ♦ INCREMENT = _F ( смотрите документ [U4.51.03] )    
   ◊ EXCIT_GENE = _F (
                      ◊FONC_MULT = fomult , [fonction_sdaster]
                      ♦ VECT_GENE = vecgen , [vect_asse_gene]       
                   ),    
   ◊ NEWTON = _F ( see the document [U4.51.03] ) ,   
   ◊ RECH_LINEAIRE = _F ( see the document [U4.51.03] ) ,     
   ◊ SOLVEUR = _F ( смотрите документ [U4.50.01] ) ,
   ◊ CONVERGENCE = _F ( смотрите документ [U4.51.03] ) ,    
   ◊ MODE_VIBR = _F (
                    ◊ NB_FREQ = / 3 , [DEFAUT]                   
                                / nbfreq , [I]       
                    ◊ MATR_RIGI = / “ELASTIQUE” , [DEFECT]                     
                                  / “TANGENTE” ,                     
                                  / “SECANTE” ,       
                    ◊ BANDE = intba , [listr8]       
                    ◊ / LIST_INST = liste_r8 , [listr8]          
                      / INST = l_r8 , [R]        
                      / PAS_CALC = npas , [I]       
                    ◊ PRECISION = / 1.E-6 , [DEFAUT]                     
                                  / prec , [R]       
                    ◊ CRITERE = / “RELATIF” , [DEFECT]                    
                                / “ABSOLU” ,              
                   ) ,    
   ◊ CRIT_STAB = _F ( смотрите документ [U4.51.03]) ,    
   ◊ ENERGIE = _F ( )    
   ◊ ARCHIVAGE = _F ( смотрите документ [U4.51.03]) ,    
   ♦ SCHEMA_TEMPS =_F (
                       ♦ SCHEMA = / “NEWMARK” ,
                                  / “HHT” ,
                                  / ”THETA_METHODE“ ,                   
                                  / “KRENK” ,                   
                                  / “DIFF_CENT” ,                   
                                  / “TCHAMWA” ,       
                       ◊ COEF_MASS_SHIFT = / 0. , [DEFAUT]                            
                                           / coef , [R]       
                       { If SCHEMA = “NEWMARK”          
                         ◊ BETA = / 0.25 , [DEFAUT]                    
                                  / beta , [R]          
                         ◊ GAMMA = / 0.5 , [DEFAUT]                  
                                   / gamm , [R]       
                       }       
                       { If SCHEMA = “HHT”
                         ◊ ALPHA = / -0.3 , [DEFAUT]                     
                                   /alph , [R]          
                       ◊ MODI_EQUI = / “OUI” ,
                                     / “NON” , [DEFECT]        
                       }       
                       { If SCHEMA = 'THETA_METHODE'          
                         ◊ THETA = / 1. , [DEFAUT]                     
                                   /theta , [R]       
                       }       
                       { If SCHEMA = 'KRENK'          
                         ◊ KAPPA = / 1. , [DEFAUT]                     
                                   / kappa , [R]       
                       }       
                       { If SCHEMA = 'TCHAMWA'          
                         ◊ PHI = / 1.05 [DEFAUT]                   
                                 / phi , [R]       
                       }       
                       { If SCHEMA = / “NEWMARK”,  
                                     / “HHT” ,                     
                                     / “THETA_METHODE” ,                     
                                     / “KRENK” ,
                                     ♦ FORMULATION = / “DEPLACEMENT” ,                                     
                                                     / “VELOCITY” ,                                     
                                                    / “ACCELERATION” ,       
                       }       
                       { If SCHEMA = / “DIFF_CENT”,                     
                                     / “TCHAMWA”,                      
                                     ♦ FORMULATION = 'ACCELERATION ' ,                     
                                     ◊ STOP_CFL = / “OUI” , [DEFECT]                                   
                                                  / “NON”,       
                       } ,
                       ) ,
   ◊ OBSERVATION = _F ( смотрите документ [U4.51.03]) ,    
   ◊ AFFICHAGE   = _F ( смотрите документ [U4.51.03]) ,    
   ◊ PROJ_MODAL = _F (      
                       ♦ MODE_MECA = mode , [mode_meca]       
                       ◊ NB_MODE = / nbmode , [I]                   
                                   / 9999 , [DEFAUT]       
                       ◊ / MASS_GENE = massgen [matr_asse_gene_R]            
                           RIGI_GENE = rigigen [matr_asse_gene_R]          
                         / AMOR_GENE = amorgen [matr_asse_gene_R]        
                       ◊ DEPL_INIT_GENE = deplgen , [vect_asse_gene]        
                       ◊ VITE_INIT_GENE = vitegen , [vect_asse_gene]        
                       ◊ ACCE_INIT_GENE = accegen , [vect_asse_gene]        
                     ) ,     
   ◊ INFO = / 1 , [DEFAUT]             
            / 2 ,    
   ◊ TITRE = tx , [kN]
 )



3 Операнды

3.1 Операнды MODELE/CHAM_MATER/CARA_ELEM

Данные операнды имеют то же значение, что и в документе (U4.51.03).

3.2 Операнды MODE_STAT/MASS_DIAG

 ◊ MODE_STAT = modestat 

Имя статического параметра, необходимого в случае сейсмического расчета с возбуждениями мульти-опоры (R4.05.01).

 ◊ MASS_DIAG = / "OUI" , 
               / "NON" , 

Опция использует явную схему [Bib2] и решает, используя диагональную смешанную матрицу масс. Эта опция доступна не для всех типов элементов, в особенности дискретных (в этом случае, надо решать с последовательной матрицей масс). С неявной схемой не возможно использовать смешанную матрицу (значительный рост ошибки).

3.3 Ключевое слово EXCIT

 ◊EXCIT = _F

Это ключевое слово позволяет описать при каждом возникновении нагрузку (напряжения и граничные условия), а также множитель и/или тип нагрузки Операнды имеют похожее значение, что и в документе (U4.51.03), но есть несколько особенностей касающихся динамики.

Важное замечание для временных схем:

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

3.3.1 Операнды CHARGES/FONC_MULT

 ♦ LOAD = chi
 ◊ FONC_MULT = fi 

Данные операнды имеют то же значение, что и в документе (U4.51.03).

3.3.2 Оператор TYPE_CHARGE

 ◊ TYPE_CHARGE = / “FIXE_CSTE” , [DEFECT]
                 / “SUIV” ,
                 / “DIDI” ,

Данные операнды имеют то же значение, что и в документе (U4.51.03), но нагрузка не может контролироваться в динамике, и поэтому tchi не может быть FIXE_PILO.

3.3.3 Операнды MULT_APPUI/ACCE/VITE/DEPL/DIRECTION/NODE/GROUP_NO

 ◊ MULT_APPUI = / “NON” , [DEFECT]
                / “OUI” ,
 ♦ ACCE = ac , [function]
 ♦ VITE = vi , [function]
 ♦ DEPL = dp , [function]
 ◊ DIRECTION = (dx , Dy , dz , drx , dry , drz) , [l_R]
 ◊ / NOEUD =lno , [l_noeud]
   / GROUP _ NO = lgrno , [l_groupe_no]

В случае возбуждения мульти-опоры (MULT_APPUI = "OUI"), у других операндов есть точно то же самое значение как в факторе ключевого слова EXCIT оператора DYNA_TRAN_MODAL (U4.53.21). В этом случае, области "DEPL", "VITE", "ACCE" соответствуют соответственно смещениям, скоростям и ускорениям относительного движения по сравнению с движением возбуждения двигателя. Новые поля "DEPL_ABSOLU', "VITE_ABSOLU', "ACCE_ABSOLU" тогда создаются и, соответственно, переписываются в смещениях, скоростях и ускорениях абсолютного движения, добавляют движение мульти-опоры двигателя и относительное движение по сравнению с этим движением мульти-опоры двигателя.

3.4 Ключевое слово CONTACT

 ♦ CONTACT = contact

Это ключевое слово позволяет активировать решение трения контакта или принятие в расчет односторонней связи. contact является понятием, следующим из оператора DEFI_CONTACT (U4.44.11). У операнда есть то же самое значение как в документе (U4.51.03).

3.5 Ключевое слово SOUS_STRUC

 ◊ SOUS_STRUC = _F

Это ключевое слово позволяет уточнить, какими являются нагрузки, которые надо использовать для статических подструктур, которые являются тогда обязательно частью модели. Данный операнд имеет то же значение, что и в документе (U4.51.03).

3.6 Ключевые слова COMP_INCR и COMP_ELAS

Синтаксис этих ключевых слов схож с несколькими командами, описанными в документе (U4.51.11). Все модели поведения, поддерживаемые в STAT_NON_LINE, доступны также и в DYNA_NON_LINE, при условии, что предусмотрен расчет матрицы масс затронутых элементов.

3.7 Ключевое слово ETAT_INIT

 ◊ ETAT_INIT = _F

Это ключевое слово определяет начальные условия задачи. Операнды данного ключевого слова имеют то же значение, что и в документе (U4.51.03).

В динамике модно определить больше значений начальных скоростей и ускорений.

 ♦ / | VITE = quickly
   / | ACCE = acce

Если ключевые слова EVOL_NOLI, DEPL и VITE отсутствуют, предполагается, что начальное состояние со смещениями, скоростями и нагрузками – нулевое, и вычисляются ускорения, соответствующее нагрузке во время instini, определенное операндом INST.

3.8 Ключевое слово INCREMENT

 ♦ INCREMENT = _F

Определяет список времен расчета. Операнды данного ключевого слова имеют то же значение, что и документе (U4.51.03).

3.9 Ключевое слово NEWTON

 ◊ NEWTON = _F

Точные характеристики итерационного метода решения нелинейных задач (метод Ньютона-Рафсона). Операнды данного ключевого слова имеют то же значение, что и документе (U4.51.03).

3.10 Ключевое слово RECH_LINEARIE

 ◊ RECHERCHE_LINEAIRE = _F (
 ◊ METHODE = / “ROPE” [DEFECT]
             / “MIXTE” )

Позволяет активировать линейный поиск. Операнды данного ключевого слова имеют то же значение, что и документе (U4.51.03), но отсутствует метод PILOTAGE.

3.11 Ключевое слово SOLVER

 ◊SOLVEUR = _F

Синтаксис этого ключевого слова такой же как и у команд, описанных в документе (U4.50.01).

3.12 Ключевое слово CONVERGENCE

 ◊CONVERGENCE = _F

Ключевое слово описывает параметры, позволяющие оценить сходимость метода Ньютона, используемого для решения нелинейных задач механики. Операнды данного ключевого слова имеют то же значение, что и документе (U4.51.03).

3.13 Ключевое слово ARCHIVAGE

 ◊ARCHIVAGE = _F

Записывает в файл все результаты или время вычислений. В отсутствие этого ключевого слова все шаги по времени архивируются, включая время вычислений с автоматически созданными далее сокращенными шагами по времени. Операнды данного ключевого слова имеют то же значение, что и документе (U4.51.03), кроме CHAM_EXCLU.

3.13.1 Операнд CHAM_EXCLU

 ◊ CHAM_EXCLU = | “DEPL” 
                | “VITE” 
                | “ACCE” 
                | “SIEF_ELGA”
                | “VARI_ELGA”

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

3.14 Ключевое слово AMOR_RAYL_RIGI

 ◊ AMOR_RAYL_RIGI = / “TANGENTE” , [DEFECT] 
                    / “ELASTIQUE”

Позволяет определить матрицу жесткости K, которая используется для построения уравнения затухания Рэлея C = α*K+β*K.

С значением по умолчанию (“TANGENTE”) матрица K будет такой же, как и при вычислении внутренних сил. Выбор значения “ELASTIWUE” запускает расчет затухания Рэлея с упругой матрицей жесткости. Для смягчения модели или типа GLRC рекомендуем использовать упругую матрицу.

3.15 Ключевое слово AMOR_MODAL

 ◊ AMOR_MODAL = _F

Данное ключевое слово позволяет принять во внимание затухание, эквивалентное модальному затуханию, разложенному на основе способов, рассчитанных через mode_meca type. Это затухание учтено в уравнении динамического баланса как поправка с вторым членом –C*X.

3.15.1 Операнды MODE_MECA/AMOR_REDUIT/LIST_AMOR/NB_MODE

 ♦ MODE_MECA = mode
 ♦ / AMOR_REDUIT = l_amor ,
   / LIST_AMOR = lisamor
 ◊ NB_MODE = nbmode

Метод mode_meca type (введенный операндом MODE_MECA) представляет просчитанную базу, на которой раскладывается модальное затухание. Эта база должна обязательно иметь тот же профиль классификации, что и в динамической системе, определенной параметрами ключевого слова SOLVEUR [3.113.11]. Можно усечь модальную базу методами, определенными с помощью NB_MODE. За неимением, берутся все способы модальной базы. Модальное затухание в сокращенной форме представлено списком условий, число которых меньше или равно числу использованных способов. Если число условий списка строго ниже, то этот список расширяется , пока его размер не достигает числа вычисленных способов.

3.15.2 Операнд REAC_VITE

 ◊ REAC_VITE = / “OUI” , [DEFECT]
               / “NON” ,

При значении “OUI” изменяет поправку силы модального затухания на каждой итерации метода Ньютона. При значении “NON” обновляет условие только в начале каждого временного шага.

3.16 Описание схемы интегрирования по времени [bib2] [R5.05.05]

 ♦ SCHEMA_TEMPS =_F ()
 ◊ STOP_CFL = / “OUI” , [DEFECT]
              / “NON” ,
 ◊ FORMULATION = / “DEPLACEMENT” ,
                 / “VELOCITY” ,
                 / “ACCELERATION” ,

Необходимо использовать неявный метод NEWMARK (ключевое слово SCHEME=”NEWMARK”, или модифицированное среднее ускорение: SCHEME=”HHT” с MODI_EQUI="NON"), метод HILBER-HUGHES-TAYLOR (SCHEME=”HHT” с MODI_EQUI="OUI") или THETA_METHODE или, наконец, схему KRENK.

В неявной схеме со смещением доступно текущие скорость или ускорение (ключевое слово FORMULATION = “DEPLACEMENT”, “VELOCITY” или “ACCELERATION”).

Также можно выбрать явный метод центральной разности (ключевое слово SCHEME=“DIFF_CENT”) или диссипативную схему типа TCHAMWA (ключевое слово SCHEME=”TCHAMWA”). Явной схемой можно получить только ускорение (ключевое слово FORMULATION = "ACCELERATION").

Явные схемы, являющиеся условно устойчивым, полезно проверить, если временной шаг, заданный в начале вычисления, оставался в устойчивом состоянии стабильности (условие CFL). Если STOP_CFL="OUI" (дефект), то, если список времен, введенных пользователем, включает один или несколько временных шагов, больших, чем условие стабильности, вычисление останавливается по фатальной ошибке. Если STOP_CFL="NON", то сигнализируется предупреждение и продолжается вычисление.

Во всех случаях критический временной шаг дан в файле сообщений.

Вычисление CFL не запрограммировано для всех элементов (в особенности, дискретные элементы проигнорированы.); CFL, оцененный Code_Aster, может таким образом быть больше (с меньшими штрафами), чем реальный CFL, с рисками грубого расхождения, которое следует из этого. В явной схеме также рекомендуется использовать смешанную массовую матрицу (диагональную): что можно получить с ключевым словом MASS_DIAG = "OUI".

3.16.1 Если SCHEMA = “NEWMARK”

 ◊ BETA = / 0.25 , [DEFAUT]
          / beta , [R]
 ◊ GAMMA = / 0.5 , [DEFAUT]
           / gamm , [R]

Это метод интегрирования по времени NEWMARK, с значениями, задаваемыми параметрами beta и gamm.

Когда не определены ни beta, ни gamm, то это метод трапеций (beta=0.25; gamm=0.5), который, на линейных задачах стабилен и не вносит паразитных потерь (т.н. числовое затухание), но который на нелинейных задачах может быть нестабилен [bib1][bib2].

3.16.2 Если SCHEMA = “HHT”

 ◊ ALPHA = / -0.3 , [DEFAUT]
           / alph , [R]
 ◊ MODI_EQUI = / “OUI” , 
               / “NON” , [DEFECT]

Для MODI_EQUI=”NON” (значение по умолчанию), метод интегрирования по времени (неявная схема интегрирования) является методом измененного среднего ускорения (семейства Newmark) ([bib1], [bib2]), с отрицательным значением alph. Чем больше alph по модулю, тем более важно числовое затухание, привнесенное вычислением. Но это разложение иногда необходимо, в нелинейных случаях, чтобы гарантировать стабильность (кроме присвоения затухания материала в конструкции). Для MODI_EQUI=”OUI”, метод интегрирования по времени (неявная схема интегрирования) является методом HILBER-HUGHES-TAYLOR (HHT или α-метод) [bib2] с отрицательным значением alph. Чем больше alph по модулю, тем более важно числовое затухание, привнесенное вычислением. По сравнению с предыдущей схемой (MODI_EQUI="NON”) измененного среднего ускорения, вызванное числовое затухание является более "выборочным": оно более слабо, с низкой и промежуточной частотой (асимптотически с ненулевой частотой), и растет более быстро с ростом частоты.

Вторая схема основана на первой с одной поправкой уравнения баланса (смещаются во времени внутренние и внешние силы) [bib2].

3.16.3 Если SCHEMA = “THETA_MOTODE”

 ◊ THETA = / 1. , [DEFAUT]
           / theta , [R]

Схема интегрирования по времени является неявной theta-схемой порядка скорости. В случае использования с контактными нагрузками нужно также вызывать непрерывный метод (AFFE_CHAR_MECA/CONTACT/METHODE="CONTINUE") и состав по скорости (FORMULATION="VELOCITY").

theta должна располагаться между 0.5 и 1:0.5 соответствующего минимума числовых потерь, theta=1 соответствует максимуму числовых потерь. theta=1 позволяет найти схему Эйлера.

Эта схема также используется в варианте со смещением.

3.16.4 Если SCHEMA = “KRENK”

 ◊ KAPPA = / 1. , [DEFAUT]
           / kappa , [R]

Это естественная диссипативная неявная схема Кренка. Ее использование таким образом рекомендуется, точно так же, как theta-схема, для неравномерных задач как, таких как ударные задачи. Числовое разложение регулируется параметром kappa, которая должна быть больше или равна 1. Если значение будет близко к 1, то схема не будет привносить потерь. Чем больше kappa, тем больше будут потери.

3.16.5 Если SCHEMA = “DIFF_CENT”

Схема центральных разностей является явной схемой второго порядка семейства Newmark, beta=0 и gamm=0.5. Данная схема не дает числовых потерь.

3.16.6 Если SCHEMA = “TCHAMWA”

 ◊ PHI = / 1.05 , [DEFAUT]
         / phi , [R]

Это альтернатива схеме центральных разностей, разработанная Бертрандом Тчамва и Кристианом Уилгошем. У этой явной схемы есть несколько интересных особенностей. Это не производная Newmark и изменение ее параметра PHI позволяет управлять числовыми потерями на высоких частотах. Когда значение близко к 1, потерь нет. Рекомендуется, исходя из ошибки, выбирать PHI между 1.10 и 1.05, чтобы достичь устойчивости, сравнимой с таковой в методе центральных разностей.

3.16.7 Операнд COEF_MASS_SHIFT

 ◊ COEF_MASS_SHIFT = / 0. [DEFAUT]
                     / coef

Данные в коэффициенте coeff позволяют выполнить смещение матрицы масс M, которая становится:

M ' = M + coef *K

По умолчанию coef имеет нулевое значение, но должно быть ненулевым, чтобы изменить схему в динамике.

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

3.17 Ключевое слово CRIT_STAB

 ◊ CRIT_STAB = _F

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

Однако необходимо отметить, что использование CRIST_STAB на сдвоенных моделях жидких структур обязывает исключать жидкие степени свободы из-за вырожденности их матриц жесткости. Поэтому необходимо комбинировать ключевые слова DDL_EXCLUS и MODI_RIGI=”OUI”. Больше информации можете найти в документах (U4.51.03), (U2.06.11) и (U2.08.04).

3.18 Ключевое слово MODE_VIBR

 ◊ MODE_VIBR = _F (
 ◊ NB_FREQ = / 3 , [DEFAUT]
             / nbfreq , [I]
 ◊ MATR_RIGI = / “ELASTIQUE” , [DEFECT]
               / “TANGENTE” ,
               / “SECANTE” ,
 ◊ BANDE = intba , [listr8]
 ◊ / LIST_INST = liste_r8 , [listr8]
   / INST = l_r8 , [R]
   / PAS_CALC = NPAS , [I]
 ◊ PRECISION = / 1.e-6 [DEFAUT]
               / prec
 ◊ CRITERE = / “RELATIF” , [DEFECT]
             / “ABSOLU” ,
                    ) ,

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

Этот критерий полезен во время переходного вычисления развития вибрирующего ответа нелинейной структуры.

Этот критерий вычисляется следующим образом: в конце временного шага вычисляется det(K-ω2,M)=0. K может или быть упругой матрицей жесткости, или связной матрицей тангенса, или секущей матрицей. M является матрицей масс. Этот модальный анализ разрешен только для симметрических матриц (масс и жесткости).

Ключевое слово NB_FREQ (3 за одну ошибку) указывает на число собственных частот к вычислению.

Собственный режим, соответствующий собственной частотете хранится в S.D. RESULTAT под именем MODE_MECA. Этот собственный режим может извлекаться и визуализироваться (как область смещений или как обычный собственный режим). Это стандартизировано с 1 на самом большом компоненте смещения.

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

Время, в течение которых необходимо произвести расчет колебательного процесса, задается списком времен LIST_INST или INST (liste_r8 или l_r8) или частотой PAS_CALC (весь npas времени). В отсутствие этих ключевых слов вибрационный модальный анализ выполняется со всеми временными шагами. Ключевые слова PRECISION и CRITERE позволяют выбрать времена, см. (U4.71.00).

3.19 Операнд ENERGIE

 ◊ ENERGIE = _F ()

Это ключевое слово позволяет активировать вычисление оценки энергии, ее показа в ходе вычисления и ее хранения в массиве PARA_CALC. Оценка энергии может быть извлечена из этого массива командой RECU_TABLE (U4.71.02).

3.20 Операнд PROK_MODAL

 ◊ PROJ_MODAL = _F

Это ключевое слово позволяет сделать вычисление на модальной основе (или Ритца) вычисленной заранее. Это должно использоваться с явной схемой интегрирования.

3.20.1 Операнды MODE_MECA, NB_MODE

 ♦ MODE_MECA = mode , [mode_meca]
 ◊ NB_MODE = / nbmode , [I]
             / 9999 , [DEFAUT]

Определяют основу, которая будет использоваться (MODE_MECA) и число способов (NB_MODE).

3.20.2 Операнды MASS_GENE, RIGI_GENE, AMOR_GENE

 ◊ / MASS_GENE = massgen , [matr_asse_gene_R]
     RIGI_GENE = rigigen , [matr_asse_gene_R]
     AMOR_GENE = amorgen , [matr_asse_gene_R]

Данные операнды используются вместе, если вы хотите динамически привести часть модели к линейному поведению, точно вычисляя с помощью DYNA_NON_LINE только области с нелинейным поведением. Это позволяет уменьшить размер модели вычисления. В этом случае необходимо вычислять модальную базу Ритца на всех областях: нелинейная область, моделируемая для расчета, вызываемого динамически, уплотняет DYNA_NON_LINE и другие линейные области. Эта база должна быть ортогонализирована по отношению к массе и с линейной жесткости всех областей. Это должно просто быть представлением движений активации всех областей. С другой стороны, MODE_MECA не будет информироваться о том, что режимы получены сокращением базы Ритца к модели вычисления, рассматриваемой DYNA_NON_LINE. Пример вычисления обеспечивается тестом SDNV107A (V5.03.107).

Операнд MASS_GENE позволяет спроецировать матрицу масс всех областей на базу Ритца с диагональным хранением. Операнд RIGI_GENE позволяет спроецировать матрицу жесткости линейных областей, сжатых только по базе Ритца с полным хранением. Операнд AMOR_GENE позволяет спроецировать матрицу затуханий (если она существует) на линейные области, сжатые только на основе Ритца с полным хранением.

3.20.3 Операнды DEPL_INIT_GENE, VITE_INIT_GENE, ACCE_INIT_GENE

 ◊ DEPL_INIT_GENE = deplgen , [vect_asse_gene]
 ◊ VITE_INIT_GENE = vitegen , [vect_asse_gene]
 ◊ ACCE_INIT_GENE = accegen , [vect_asse_gene]

Данные операнды связаны с использованием операндов MASS_GENE, RIGI_GENE и, возможно, AMOR_GENE в ключевом слове PROJ_MODAL. Это используется, чтобы ввести обобщенный вектор, следующий из проекции (с помощью PROJ_VECT_BASE (TYPE=”DEPL”)) смещения области, или скорости или ускорения полной модели (включая область с линейным поведением) на основе этой полной образцовой модальной базы. Этот обобщенный вектор будет действовать как начальное условие развития обобщенных координат вычисления на масштабной модели с нелинейной областью. Тогда будет необходимо также сообщить операндам ключевого слова ETAT_INIT с соответствующим смещением областей, скоростей, ускорений и напряжений, пробную переменную небольшой модели. Пример вычисления обеспечивается испытательным SDNV107C (V5.03.107).

3.21 Ключевое слово EXCIT_GENE

 ◊ EXCIT_GENE = _F (
 ◊ FONC_MULT = fomult , [fonction_sdaster]
 ♦ VECT_GENE = vecgen , [vect_asse_gene]
                   )

Это ключевое слово связано с использованием операндов MASS_GENE, RIGI_GENE и возможно AMOR_GENE в ключевом слове PROJ_MODAL. Это используется, чтобы ввести приложенные силы на областях линейного поведения, сжатого динамически и не смоделированного в расчет, путем вызова явной схемы интегрирования. Эти силы проецируются на базу Ритца, вычисленной на всех областях.

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

3.22 Операнд INFO

 ◊ INFO = inf

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

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

Внимание, файлы .mess могут стать очень важными при INFO = 2.

3.23 Операнд TITRATES

 ◊ TITRE = tx

tx является названием вычисления. Это будет напечатано наверху результатов. См. (U4.03.01).



4 Пример: движение маятника большой амплитуды

 #TITRE simple Pendulum in great oscillation
 #
 #PENDULE CONSTITUE Of an ELEMENT OF CABLE (test SDNL100A).
 #
 RESU=DYNA_NON_LINE (MODELE=MO, CHAM_MATER=CHMAT, CARA_ELEM=CARA, 
                     EXCIT=(_F(LOAD=CHA1),
                            _F(LOAD=CHA2)), 
                     INCREMENT=_F(INST_INIT=0.,LIST_INST=L_INST1), 
                     ARCHIVAGE=_F(LIST_INST=L_INST2), 
                     SCHEMA_TEMPS=_F(SCHEMA='NEWMARK”, 
                                     FORMULATION=”DEPLACEMENT”), 
                     COMP_ELAS=_F(RELATION=“CABLE”, 
                                  DEFORMATION=“GREEN”), 
                     CONVERGENCE=_F(RESI_GLOB_RELA=1.E-6, 
                                    ITER_GLOB_MAXI=100), 
                     NEWTON=_F(REAC_ITER=1)
                    )
 FIN()
  • груз cha1 вынуждает node 1 остаться фиксированным с узлом 2, чтобы переместить в вертикальную плоскость XZ
  • груз cha2 – сила тяжести
  • команда DYNA_NON_LINE определяет:
    • будет использоваться метод интегрирования ”NEWMARK” – Метод трапеций (также называемый средним ускорением), т.к. у “NEWMARK” нет ни одного аргумента
    • Начальное состояние, в момент времени 0, имеет нулевое смещение; смещения будут рассчитывая начиная с исходного положения и нулевой скорости
    • итерационный расчет будет продолжаться до тех пор, пока относительная погрешность будет больше 10-2, но число итераций ограничено числом 100
    • конечная матрица тангенсов линейной системы будет пересчитываться на каждой итерации (ошибка начиная с ключевого слова “NEWTON” отсутствует).