[fasol]

Нужно построить спектр ответа в заданной точке конструкции при действии акселерограммы (полноценный нестационарный расчет).
Спектры ответа строю функцией RESP, с ней все понятно:
RESP, IR, номер_таблицы_частот, LDTAB, 3, 0.05, 0.02, 0, 15

На вход она, помимо всего прочего, принимает таблицу частот, заполненной FILLDATA или DATA. С этим сложнее:
FILLDATA, номер_таблицы_частот, LSTRT, LSTOP, LINC, VALUE, DVAL
Пробовал оставлять последние 5 параметров пустые, пробовал заполнять - разницы особой нет.
Построенный спектр очень сильно отличается от результатов работы скадовского Редактора Акселерограмм, если скормить тому поэтажную акселерограмму.
Что я делаю не так?

[СергейД]

трудно без контекста советовать..

ktm=1000! точек на графике (можно запросить *get)
fmax=20! от 0.1 до 20 гц частоты учитываем
dfr=fmax/ktm! шаг по частотам
plvar,2! график ускорений
filldata,10,,,,0.1,dfr
plvar,10!проверьте таблицу

dt=0.02! какой шаг у вас ?
xvar,10! по абсциссе дб частоты

demp=0.3/6.28! при декременте 0.30

resp,3,10,2,3,demp,dt,,,1! спектр ускорений
plvar,3

[Бахил]

Эта задача решается элементарно в екселе. Строчек на бейсике будет не больше, чем в ансисе.

[СергейД]

тут несколько строк просто для читабельности
в принципе всего- три строки.

filldata,10,,,,0.1,1000
xvar,10! по абсциссе дб частоты
resp,3,10,2,3,0.05,0.02,,,1! спектр ускорений

все троллишь, Бахил?

[Бахил]

Сообщение от СергейД все троллишь, Бахил?

Да. На бейсике побольше будет.

зато более понятно

[СергейД]

кому как...
ты уверен, что это то что нужно???

у человека макрос под ансис со всей процедурой расчета. модель-расчет-обработка
речь идет о том, чтобы в нее добавить несколько строк. точнее, поправить.

зачем тут что-то передавать в эксель. там программировать, потом обратно в ансис загружать???
я понимаю еще питон применить.

[Бахил]

Offtop: Да шучу я. Не бери в голову.

[СергейД]

кстати, планируем организовать обучение ansys на курсах.
но я только на базе apdl смогу научить..
WB потом самому можно изучить.

[fasol]

Сообщение от СергейД filldata,10,,,,0.1,1000 xvar,10! по абсциссе дб частоты resp,3,10,2,3,0.05,0.02,,,1! спектр ускорений

Спасибо, стало намного лучше.
Как оказалось, на спектр сильно влияет практически все, а особенно максимальная частота и количество точек в графике, см. вложение.
Причем при увеличении количество точек (больше 1000), становится еще хуже - пик переезжает в другую половину графика.
В итоге наиболее точно ансис со скадом сходится при максимальной частоте 25 Гц и количестве точек как в акселерограмме (756).
Там же просто разложение Фурье, откуда такая чувствительность ко входным данным?

[Бахил]

Разный масштаб абсциссы.

[fasol]

Сообщение от Бахил Разный масштаб абсциссы.

Лишь бы ляпнуть, да?

[Бахил]

По твоим графикам только одно можно сказать: СКАД считает правильно.
А с ансисом ты не разобрался.

----- добавлено через ~2 мин. -----
Да, и судя по графикам, по абсциссе отложена не частота, а период.
Можешь проверить по элементарной процедуре. Ну если разберёшься.

[СергейД]

число точек = число шагов решения=756
по горизонтали должны быть частоты.

вышли в личку акселерограмму.

[fasol]

Выслал.

[Бахил]

Offtop: ...

[СергейД]

все очень хорошо совпало со скадовскими графиками.

[fasol]

Методом тыка было установлено, что максимальная частота никак не влияет на график спектра.
Заметно влияет шаг RESP, IR, номер_таблицы_частот, LDTAB, 3, 0.05, 0.02.
Причем при совпадении шага с дискретностью акселерограммы (0.02) сходимость со СКАДом наименьшая.

[СергейД]

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

[fasol]

Сообщение от СергейД при уменьшении шага до 0.005с и менее точность разложения в ансис возрастает до необходимой

Я как раз это на графике и привел.
0.005 почти не отличается от 0.01, фаза практически совпадает со скадовской, а вот амплитуда на 20% меньше.
Дальнейшее уменьшение шага ухудшает картину.

[Бахил]

Интересно взглянуть на акселерограмму.
Шаг вообще-то не должен сильно влиять.

[СергейД]

там частота 10 гц. период 0.1 с. при 0.02 с пять точек на две полуволны. естественно, график с изломами. пики пропускаются.

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


при 0.005 20 точек на период- график более менее гладкий.
я просто для наглядности сгенерировал синусоиду 10 гц и сохранил с различными шагами по времени.
все очевидно.

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

[fasol]

Сообщение от СергейД я просто для наглядности сгенерировал синусоиду 10 гц и сохранил с различными шагами по времени.

Покажите, если не трудно, как у вас сходятся спектры акселерограммы, что высылал вам.
Аппроксимации в скаде точно нет.

[СергейД]

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

[fasol]

Сообщение от СергейД про скад откуда вы знаете что у них внутри?

Потому что точек на графике там заметно меньше.

Синусоида у меня тоже более-менее сходится (шаг акселерограммы 0.02).

[fasol]

Есть у кого-нибудь возможность прогнать акселерограмму в Одисее или схожем ПО?
Появилось подозрение, что Скад временами привирает.
Акселерограмма - https://pastebin.com/DzeVrphU
Результирующий скадовский спектр 5% - https://pastebin.com/TKMX3nXP

[СергейД]

опять 0.02...
при таком грубом шаге будут все те же коллизии

затребуйте с меньшим шагом или сами ее интерполируйте сплайном

[fasol]

Сообщение от СергейД опять 0.02...

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

[Бахил]

11,6 Гц резонанс 302

[СергейД]

в ансис обрабатывается график с линейным изменением между точками.
если шаг 0.02с (50 гц) а функция изменяется в основном 10-12 гц,
то как вы две полуволны правильно опишете пятью точками? будут изломы графика.

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

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


2 бахил. вы по что сейчас?. своей программой считали или какой?
создайте синусоиду10 гц, запишите ее с шагом 0.02с и определите ее спектр при нулевой и ненулевом демпф. что получится?

[fasol]

Сообщение от СергейД получите той или иной интерполяцией перед ансис дополнительные точки на графике и считайте себе спектр в ансис.

Попробую.

Сообщение от СергейД в ноль остаточные перемещения приходят?

Нет, акселерограмма не сбалансированная, как и любая синтезированная. Но это никак не влияет на спектр, я проверял.

[Бахил]

Offtop: Кнопкотыки...
0,02 - шаг периода для которого считается максимальная реакция линейного осциллятора.
Для данной акселерограммы максимальная реакция осциллятора= 302 при его частоте 11,6 Гц. Для периода 0,1с примерно 297.
Странная акселерограмма. Бункер считаешь?

[СергейД]

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

вопрос и в использовании этой акс-ммы для расчетов отклика.
какие частоты у нагружаемого оборудования?
волновать должны пики спектра нагрузки в зоне низших собственных частот.
если эти частоты 30-50 гц, то разницы при пиках спектров 9 или 11 гц особой не заметишь.

[Бахил]

Да.

[fasol]

Сообщение от СергейД получите той или иной интерполяцией перед ансис дополнительные точки на графике и считайте себе спектр в ансис.

Сделал линейную интерполяцию результирующей сейсмограммы, увеличив количество точек в два раза.
Потом через псевдозадачу как у вас снял NSOL и опять прогнал RESP - результат еще хуже, фаза на месте, а вот амплитуда просела на 20%.
Слишком трудоемкая задача оказалась, возвращаюсь к инструментам скада.

Сообщение от СергейД какие частоты у нагружаемого оборудования? волновать должны пики спектра нагрузки в зоне низших собственных частот. если эти частоты 30-50 гц, то разницы при пиках спектров 9 или 11 гц особой не заметишь.

Такой подход меня не устраивает, это уже шаг назад.

[СергейД]

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

[fasol]

Сообщение от СергейД гладким сплайном.

А функционал APDL позволяет это сделать?
Этож нужно ловить несколько точек до и несколько после рассматриваемого интервала и аппроксимировать полиномом?

[СергейД]

сделайте в матлабе или еще где-то.

в ансис можно изловчиться-извратиться= вариантов много.
например, можно цифры в координаты KP превратить и провести сплайн (можно до 200 точек в одной порции). разбить линии ldiv и сохранить координаты всех точек.

[Бахил]

Интерполировать исходную акселерограмму бесполезно