Интерполяционный полином

[Kikurato 0]

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

 

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

​

[Борман 2 375]

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

Можно чуть подробнее, чему вы хотите научится ? Писать собственные типы конечных элементов, или создавать КЭ-модели на встроенных элементах ?

 

 

 

простроение на основе файлов

Вы имеете ввиду командных файлов ?

 

PS. Правильно писать "NASTRAN". 

[Kikurato 0]

Создавать новый тип конечного элемента, или вносить изменение в полином одного из предоставленных (читал хелп по 12 версии, там были пользовательские типы конечных элементов, в которые можно было вносить изменения в саму структуру)

 

 

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

Изменено 27 августа 2015 пользователем Kikurato

[dbarlam 184]

Процесс следующий - садимся и пишем на бумажке полиномы / либо любые другие интерполяционные функции на Ваше усмотрение

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

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

для удобства проведения численного интегрирования и беря производные записываем матрицу жесткости элемента

А ее программируем - все в литературе описано , но первый этап - на бумаге

[Kikurato 0]

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

[Fedor 1 597]

Зачем на бумаге то. Еще лет 15 назад получал в Mathematica.  Там можно сишные или фортрановские коды сразу и заказать :)

К сожалению ссылки поудаляли мои, ну да найдете на экспоненте ру  :)

[Andron 5]

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

Programmer's Manual for Mechanical APDL уже читали?

[Kikurato 0]

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

[Andron 5]

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

Может у Вас другой документ( а может враги его испортили  ), но там вообще-то идет речь о коде на языке ФОРТРАН( и/или С), изменение текста которого самого по себе ничего не даст. Как любой программный код, он должен быть откомпилирован и слинкован (см. пп. 5.9 Compiling and Linking UPFs on Unix/Linux Systems или 5.10 Compiling and Linking UPFs on Windows Systems - название может отличаться в зависимости от используемой версии Ансис ). Таким образом, Вы получите на выходе новый исполняемый файл Ансиса.

 

P.S. Не знаю как в последних версиях Ансис - а на версиях 11.0 - 13.0 нужно было дополнительно компилятор фортрановский ставить определенных версий. Рекомендую Вам внимательно читать документацию - там описано все довольно подробно.

[Fedor 1 597]

Для С++ нужен VS 2010  как написано в документации к 16.

[Chardash 306]

Всем доброго времени суток! краем уха услышал речь про Mathematica. Может подскажете, как отобразить в фазовом портрете кривые, получились траектории в моем седле. Извиняюсь, что не в тему. ) На спец форумах все спят

ура) тут все так же и Fedor здесь 

 

+ тут нужно еще сепаратистов сепаратрисы нарисовать.

 

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

 

{-2 x + (6/7) y, 14 x - 3 y}

[Fedor 1 597]

Нарисуйте отдельные кривые и через Show[... ] изобразите их на одном графике...
Пример :
Show[Plot[x^2, {x, 0, 3.5}], ListPlot[{1, 4, 9}]] и вроде можно вкладывать Show друг в друга типа рекурсивности...
Я так обычно парамертически строю M-N огибающую кривую прочности армированной колонны и из ансиса беру кучу точек N-M слегка обрабатываю в ультраэдите и все рисую на одном листочке через Show. Сразу гарантия что все колонны будут гарантированно прочными. При этом иногда и еще множество различных нагружений в одной куче использую. :)
Иногда разбиваю на несколько подмножеств кривых армирование, чтобы не все колонны одинаково армировать ... И все на одном листочке помещается. Очень удобно и наглядно

https://en.wikipedia.org/wiki/Phase_planeтут гляньте, как-то рисуют...

https://en.wikipedia.org/wiki/Phase_portrait

http://onlinesciencetools.com/tools/phaseportrait  можете попробовать тут нарисовать, а потом ножницами вырезать и вставить куда надо. Например в ту же математику  :)

[Chardash 306]

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

Пример :

Show[Plot[x^2, {x, 0, 3.5}], ListPlot[{1, 4, 9}]] и вроде можно вкладывать Show друг в друга типа рекурсивности...

Я так обычно парамертически строю M-N огибающую кривую прочности армированной колонны и из ансиса беру кучу точек N-M слегка обрабатываю в ультраэдите и все рисую на одном листочке через Show. Сразу гарантия что все колонны будут гарантированно прочными. При этом иногда и еще множество различных нагружений в одной куче использую. :)

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

https://en.wikipedia.org/wiki/Phase_planeтут гляньте, как-то рисуют...

https://en.wikipedia.org/wiki/Phase_portrait

http://onlinesciencetools.com/tools/phaseportrait  можете попробовать тут нарисовать, а потом ножницами вырезать и вставить куда надо :)

Fedor, Спасибо Вам за ссылки, особенно последняя. Буду пробовать, но методом тыка этот калькулятор не заработал. Я сделал расчет врукопашку, да, отображаю через Show. Но рисует кривые только в одной чертверти, уравнения сепаратрис вручную задал в математике и так же их отобразил. Хорошо бы, если все автоматом...

 

Если нужно скину код. Нашел примеры в справке, я так понял у меня проблема с  t. 

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

[Fedor 1 597]

Если я правильно понял, Вы привели правую часть. Попробуйте через ParametricPlot рисовать. Я подумал что в векторе правая часть системы ДУ приведена, решил через DSolve и для одной ветки рисовал. Вроде и в другую сторону внизу загибались...

[Chardash 306]

да, так и делаю

eqv = DSolve[{x'[t] == -2 x[t] + (6/7) y[t] && y'[t] == 14 x[t] - 3 y[t], x[0] == c, y[0] == c}, {x[t], y[t]}, t]
r1 = ParametricPlot[Evaluate[Table[{x[t], y[t]} /. eqv[[1]], {c, 10}]], {t, 0, 1}, 
  PlotRange -> All, AspectRatio -> 1]

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

[Fedor 1 597]

x[0] == c, y[0] == c думаю что из - за этого. Выкиньте и получите решение с C[1] и С[2]. Их пробуйте менять чтобы ветки загибались.

DSolve[{x'[t] == -2 x[t] + 6/7 y[t], y'[t] == 14 x[t] - 3 y[t]}, {x[t], y[t]}, t] я так делал. И ручками переопределял C[1] и С[2] в решении из экспонент и потом их менял на числовые значения...

eqv = DSolve[{x'[t] == -2 x[t] + (6/7) y[t] &&

y'[t] == 14 x[t] - 3 y[t], x[0] == c, y[0] == c}, {x[t], y[t]}, t]

r1 = ParametricPlot[

Evaluate[Table[{x[t], y[t]} /. eqv[[1]], {c, -10, 10}]], {t, -0.3,

1}, PlotRange -> All, AspectRatio -> 1] так попробуйте :)

[Chardash 306]

большое спасибо, Fedor, пытался через массив константы задать, не могу пока быстро, если их несколько. С одной проще  C[1] -> a  и все. Отличные графики, вручную тоже получилось) Время с -1 брал, чуть не  подобрал.

 

Кстати, отличный бесплатный редактор для ANSYS Sublime3. 

https://www.meshparts.de/blog/Dem_SublimeText_Editor_die_Ansys_APDL_Hilfe_beibringen

https://www.meshparts.de/blog/Definitely_one_of_the_best_text_editors_for_Ansys_APDL

[Fedor 1 597]

a =.; c =.;

eqv = DSolve[{x'[t] == -2 x[t] + (6/7) y[t] &&

y'[t] == 14 x[t] - 3 y[t], x[0] == a, y[0] == c}, {x[t], y[t]}, t]

r1 = ParametricPlot[

Evaluate[

Table[{x[t], y[t]} /. eqv[[1]], {a, -1, 1}, {c, -4, 4}]

],

{t, -0.2, 1.0}, AspectRatio -> 1]

Попробуйте так, что-то разумное получается :)

Спасибо за ссылки, но я привык к UE, там и для ансиса можно настроить и главное можно со столбиками работать, что удобно и есть язык скриптов с помощью которого забираю и автоматически преобразую данные из автокада. Где-то в железобетоне приводил коды :)

[Chardash 306]

ну вообще круто. теперь как в книге 

 

да, у меня экзамен скоро(пересдаю диффуры, не успел в прошлом семестре), поэтому рад вдвойне.

 

 

 

Где-то в железобетоне приводил коды :) 

да, я помню.

[Fedor 1 597]

Удачи. Спасибо развлекли. Я последний раз дифурами лет 30-35 назад занимался. Уравнение Ван Дер Поля и всякая нелинейная механика со странными аттракторами  :)

https://ru.wikipedia.org/wiki/Осциллятор_Ван_дер_Поля

https://ru.wikipedia.org/wiki/%D0%A2%D0%B5%D0%BE%D1%80%D0%B8%D1%8F_%D1%85%D0%B0%D0%BE%D1%81%D0%B0:)