оставаться в /solve или использовать Restart?

[mazam 0]

Уважаемые господа,

 

Пожалуйста помогите с нелинейной задачей, так как уже довольно долго не могу решить.

Еще прошу снисходительности, так как опыта не так уж много у меня в модлировании.

 

Моя проблема состоит в том, что я хочу разделить решение на 2 этапа.

 

1ый этап - задание нагрузки и constrains и получение решения о нелинейном перемещении.

2ой этап - на уже деформированное тело, задать нелинейные перемещения, зависимые от предыдущего решения, и изменить constrains

 

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

 

все время оставаясь в /solve

или через Рестрат

 

В чем проблема, не могу понять.

 

Прмерный код с restart:

!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

finish
/clear,start
/title,experiment_test1

/units,si !Use SI unites
/filnam,experiment_test1

/prep7

 

. . . . !Создание геометрии, meshing, задание нагрузки и DOF 

finish

 

!!!!!!1st Solution  

/solu
antype,static
deltim,0.01,1E-3,.5
 outpr,basic,last        
outres,all,all
rescont,define,all,all,1
cnvtol,m,5e-3
nlgeom,on
solve
/eshape,1
/post1                                                      !получение перемещений после первого этапа
plnsol,u,sum,2,1
nsel,s,loc,x,0
*get,node_numb,node,0,num,max               
*get,now_x,node,node_numb,u,x
*get,now_y,node,node_numb,u,y

  finish

parsav,all

 

 

!!!!2nd Solution with RESTART

 

/clear,nostart
      /filnam,experiment_test1

      /solu
      nlgeom,on
      antype,,rest
      ddele,all

      ...   !Новая нагрузка и новые DOF (Основываясь на данных *get о перемещениях после первого этапа )

   ! Вопрос:  Если выбирать узлы по xyz location положению, то  считать, что тело уже деформированно и переместилось или как если бы первого решения не существовало?

 

       allsel,all
!       set,2
      solve
      finish
      /post1
      set,2
      plnsol,u,sum,2,1
      /eshape,1

/eof

!!!!!!!!!!!!!!!!!!!

 

 

Прмерный код, не выходя из /solve:

!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

finish
/clear,start
/title,experiment_test1

/units,si !Use SI unites
/filnam,experiment_test1

/prep7

 

. . . . !Создание геометрии, meshing, задание нагрузки и DOF 

finish

 

!!!!!!1st Solution            !Производится только для получения данных в постпроцессоре

/solu
antype,static
deltim,0.01,1E-3,.5
 outpr,basic,last        
outres,all,all
rescont,define,all,all,1
cnvtol,m,5e-3
nlgeom,on
solve
/eshape,1
/post1                                                      !получение перемещений после первого этапа
plnsol,u,sum,2,1
nsel,s,loc,x,0
*get,node_numb,node,0,num,max               
*get,now_x,node,node_numb,u,x
*get,now_y,node,node_numb,u,y

  finish

 

!!!!2nd Solution  проделать то же решение, но уже не выходя в постпроцессор

 

      allsel,all
      /solu
      nlgeom,on
!      time,1
!      nsubst,50
      solve

 

!!!!!3rd Solution    - Задание новых нагрузок, перемещений  и граничных условий на основе *get из первого решения
!      time,2
!      nsubst,100
      nlgeom,on
      ddele,all
 
      ...   !Новая нагрузка и новые DOF (Основываясь на данных *get о перемещениях после первого этапа )

      

      allsel,all
      solve
      /post1
      set,2
      plnsol,u,sum,2,1 ! 
      allsel,all

 

!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

 

В обоих случаях АНСИС не видит результатов после первого решения. В чем может быть проблема?

 

 

[soklakov 1 476]

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

А как должно быть по-Вашему?

[mazam 0]

По-моему,

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

Изменено 22 января 2014 пользователем mazam

[soklakov 1 476]

Это понятно. Но если Вы убрали граничные условия, то исчезнут и напряжения ими вызванные. Вы уверены, что не запутались?

[mazam 0]

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

 

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

 

Простите, если немного непонятно изъясняюсь.

Так как же поступить тогда, если для второго решения мне constrains первого решения не нужны?

Изменено 22 января 2014 пользователем mazam

[Makarrr 15]

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

[soklakov 1 476]

сохраняя напряжения и деформации из решения А, я создаю новые граничные условия и нагрузки Б

И с чего вдруг напряжения и деформации полностью сохраняются?

Как заметил @@Makarrr, не хватает ясности физики процесса.

[mazam 0]

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

Имеется гибкая полоска неоднородного материала, по свойствам подобная резине. Ее надо свернуть в трубу (без нахлеста).

Данные соответствуют реальному образцу.

Если задать перемещения точек (плоская полоска - труба), то АНСИС будет отправлять узлы по наименьшему расстоянию. И поэтому, вся задача разваливается.

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

  

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

 

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

 

Если у вас есть другие идеи, как решить задачу, я буду очень рада.

 

 

[soklakov 1 476]

Освободить оба конца и удалить g, зафиксировать середину, и дать перемещения узлам

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

Но отпуская модель, Вы ее разгружаете, и она, естественно, выпрямляется.

[dolivanov 7]

 

Освободить оба конца и удалить g, зафиксировать середину, и дать перемещения узлам

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

Но отпуская модель, Вы ее разгружаете, и она, естественно, выпрямляется.

 

А как же галочка в prep7 в настройках решателя "calculate prestress effects"? Включение этой опции позволяет передавать НДС нагруженной конструкции для решения задачи устойчивости. Разве здесь это не будет работать?

[soklakov 1 476]

Разве здесь это не будет работать?

Думаю, нет. Но лучше попробовать и убедиться, чем не пробовать и сомневаться.

[mazam 0]

Вы имеет ввиду, вот так?

/prep7

...

/solu

solve

/post1

*get   !Get displacements

 allsel,all
      /solu          !Solve again
      nlgeom,on
      pstres,on

      solve

 nlgeom,on
      ddele,all
      acel,,0
     nsel,s,....

     d, ...   *на основе получ *get

...

     solve
      /post1

 

Я пробовала, все равно не видит первое НДС в ходе первой  задачи

 

Попробую подумать, как обойтись без ddele.

[Makarrr 15]

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

Я так попробовал сделать для вашей задачи. Вроде получается, но очень долго...

И исчё: задачка симметричная, берите половинку.

[mazam 0]

Спасибо огромное за помощь и советы. На самом деле свернуть полоску из плоского состояния в круглое - получается нет проблем. Без всяких промужуточных решений. Я использовала beam188, keyopt,1,3,2. Для такого элемента ROTZ пропорционален приложенному моменту. То есть я просто фиксирую центральную часть и прикладываю момент +/- EI/R , где R - это радиус к концам полоски. nlgeom, on. Делю нагрузку на 100 шагов и задача решается. 

 

Сейчас у меня следующая проблема. Вернее 2.

1) Я не могу запустить рестарт. После того, как полоска свернута, сохраняя НДС, я дополнительно фиксирую 6 узлов и прикладываю на концы - моменты, равные нулю. Тогда полоска должна распрямиться и зафиксированные узлы показать реакции опоры. Так вот, что-то не так с кодом, так как НДС почему-то не сохраняется после первого решенияб и задача не решается.

Вот примерный код:

 

finish
/clear,start
/title,internal stiffness
/units,si !Use SI unites
/filnam,internal_stiffness

/prep7

 ...!Строю геометрию, делю на конечные элементы, фиксирую центральный узел и прикладываю моменты на концах

/solu
antype,static
nlgeom,on
deltim,0.01,1E-3,.5          
outres,all,all
rescont,define,all,all
cnvtol,m,5e-3
rescont,define,all,all
nsubst,100
autots,on
solve
finish
/post1
plnsol,u,sum,2,1
/eshape,1 

!Second move

/solu
time,2
antype,,rest,1,15         !Открываю НДС полоски после сворачивания из плоского состоянив трубу
nlgeom,on

...! Фиксирую 6 точек по периметру и прикладываю моменты, равные нулю - освобождаю полоску

 autots,on
outres,all,all
solve
finish
/post1
set,last
plnsol,u,sum,2,1
/eshape,1

 

Задача не решается, в чем проблема? Что-то не так с DOF?

 

2) Если я хочу смоделировать вместо 6 точек цилиндрическую камеру, как показано на рисунке выше, что я должна делать? Я думала, опять использовать Рестарт и смоделировать контакную задачу, и использовать Birth and Death Elements. Возможно ли это или лучше решать задачу, как вы посоветовали с прессом?

Если нужно, могу приложить полный АПДЛ код.

[Mausing 0]

надо свернуть в трубу (без нахлеста).

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

 

 

Если нужно, могу приложить полный АПДЛ код.

 

Будьте любезны.

Изменено 14 апреля 2014 пользователем Mausing