Перейти к публикации

FFT, Ceйсмика, PSD, Transient


Рекомендованные сообщения

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

Это все пустое. Я подобное пропускаю мимо ушей.

Вы бы лучше обдумывали свои ответы.

У вас получается, чем больше возбуждающая частота, тем больший шаг по времени нужно брать. :unsure:

Вопрос остался без ответа.

Расчет проводили с полной матрицей жесткости?

Ссылка на сообщение
Поделиться на других сайтах


По теме.

FFT - очень быстрое преобразование Фурье. Оно так и называется (Fast Fourier transform).

Быстрое оно за счет того, что в качестве исходных данных берется 2n временных отчетов, например, 2048 или 32768.

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

Здесь лучше использовать cfft, которое намного медленней, но требует 2*n временных отчетов сигнала.

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

Спасибо вам, сейчас попробую. О результатах отпишусь.

Господа, развели дискуссию :wallbash: очень сложно найти действительно что-то ценное.

Дмитрий, расскажите подробнее как вы делали этот расчет, сделайте полную выкладку. Это будет гораздо полезнее статьи. Говорить великую фразу -"Read the fuc...g manual" - я тоже могу. Но если вы все разложите по полочкам, скажете на какую кнопку нажать, то вы сдвинете большой камень непонимания. Да, возможно, ваш метод расчета будет не самым лучшим и 100% верным, но ошибки которые были сделаны там, помогут не делать их в будущем тем кто только изучает МКЭ и данное ПО. С уважением.

Ссылка на сообщение
Поделиться на других сайтах

"скажете на какую кнопку нажать" - в расчетах важно понимание, а не число - учат классики. Вы как тетеньки - конструкторши из Промстройпроекта, когда-то, рассуждаете, которых зачем-то за Скад засадили, а они бедняги уже и механику давным-давно позабыли. С таким же вопросом обычно приходили :rolleyes:

Ссылка на сообщение
Поделиться на других сайтах

"скажете на какую кнопку нажать" - в расчетах важно понимание, а не число - учат классики. Вы как тетеньки - конструкторши из Промстройпроекта, когда-то, рассуждаете, которых зачем-то за Скад засадили, а они бедняги уже и механику давным-давно позабыли. С таким же вопросом обычно приходили :rolleyes:

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

Ссылка на сообщение
Поделиться на других сайтах

Да какая тут особая физика. Поставьте бутылку пива на стол да стучите сбоку по нему, вот и вся физика. Если поставите на мягкую подставку - упадет. Решите эту задачу и совпадет с экспериментом, можете с чистой совестью и домики считать :rolleyes:

Ссылка на сообщение
Поделиться на других сайтах

Да какая тут особая физика. Поставьте бутылку пива на стол да стучите сбоку по нему, вот и вся физика. Если поставите на мягкую подставку - упадет. Решите эту задачу и совпадет с экспериментом, можете с чистой совестью и домики считать :rolleyes:

:rolleyes: А теперь загоните это в Ansys!) И скажите - Жми туда, и будет счастье!))

А серьезно, скиньте макрос применения команды RESP в Wb с выводом результатов в нем же. Или он выводит это дело в Outpute? или я подгружаю картинку из APDL?

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

Ссылка на сообщение
Поделиться на других сайтах

В борьбе обретешь ты счастье свое - был лозунг у социал-революционеров. Актуален и сейчас. :rolleyes:

Кстати на www.fea.ru было описание как трясли кирпич стоящий на меньшей грани на вибрирующем столе. Модель, зависимость от частоты когда упадет.

Ссылка на сообщение
Поделиться на других сайтах

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

Если это можно сбросить в Word, то там уже и до Excel'я недалеко.
Ссылка на сообщение
Поделиться на других сайтах
Дмитрий Фролов

Дмитрий, расскажите подробнее как вы делали этот расчет, сделайте полную выкладку. Это будет гораздо полезнее статьи. Говорить великую фразу -"Read the fuc...g manual" - я тоже могу. Но если вы все разложите по полочкам, скажете на какую кнопку нажать, то вы сдвинете большой камень непонимания. Да, возможно, ваш метод расчета будет не самым лучшим и 100% верным, но ошибки которые были сделаны там, помогут не делать их в будущем тем кто только изучает МКЭ и данное ПО. С уважением.

выслал вам в личку свой макрос для получения спектра, пользуйтесь)

Ссылка на сообщение
Поделиться на других сайтах

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

Попробуйте через экспорт

Экспортирует всю матрицу в виде форматированного текста в указываемый далее файл.

Ссылка на сообщение
Поделиться на других сайтах

начальный 0.0025, минимальный 0.00002. если вы эту информацию не найдете в хелпе ансиса, то что вы вообще тут делаете?

Резкое высказывание.

Такой информации в хелпе нет и быть не может. Это же входные данные проведенного Вами расчета. :unsure:

А на форуме я умножаю свои знания.

Ссылка на сообщение
Поделиться на других сайтах

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

Ссылка на сообщение
Поделиться на других сайтах

Все зависит от частоты внешней назрузки.

И соответсвенно это тесно связано с количеством собственных частот, участвующих в расчете.

Это если есть желание уменьшить время решения и не возиться с полными матрицами масс и жесткости. :unsure:

Ссылка на сообщение
Поделиться на других сайтах

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

Ссылка на сообщение
Поделиться на других сайтах

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

Я бы тогда говорил о минимальном периоде колебаний.

Насколько помню аналитическое решение есть и для импульса... :unsure:

Ссылка на сообщение
Поделиться на других сайтах

Для функции Дирака есть - называется фундаментальным. Очень удобно для бесконечных и полубесконечных областей типа грунта использовать. С помощью логарифмической линейки элементарно осадки зданий, кварталов, городов считать <noindex>http://www.pinega3.narod.ru/mg/mg.htm</noindex> :rolleyes:

Ссылка на сообщение
Поделиться на других сайтах

С помощью логарифмической линейки элементарно осадки зданий, кварталов, городов считать

Аналоговый компьютер всегда был дешевле и доступнее цифрового. :unsure:

Россия - матушка.

Ссылка на сообщение
Поделиться на других сайтах

Господа, использовал ваши рекомендации и провел сравнение двух методов для сейсмики. Transiet и RS. Делал это через WB, просто так удобнее. Спешу поделится результатами и интересно узнать ваше мнение.

Брал произвольную функцию ( в табличном виде решил не задавать, на данном этапе)

post-35250-1354777914.jpg

использовал fft функцию. ( быстрое преобразование фурье). cfft пока не использовал, было лень возится с табличными данными.

post-35250-1354778065_thumb.jpg

Считал в Transient анализе приложил ускорение в табличном виде ( было много мата ,по поводу вставки таблиц и работы со временем)

post-35250-1354778148_thumb.jpg

post-35250-1354778184_thumb.jpg

post-35250-1354778425.jpg

Полученные результаты 14 МПа.

Потом данную модель, с такой же сеткой, рассчитал в модальном анализе, определил собственные частоты

post-35250-1354778215.jpg

Приложил линейно-спектральную диаграмму. WB больше 10 точек по какой то причине отказывался прикладывать. Главное нужно было установить Scale фактор 9,81 ( понимаю что связано с гравитацией, но почему?)

post-35250-1354778197_thumb.jpg

И результаты практически полностью совпадают

post-35250-1354778208_thumb.jpg

Время расчета Transient

CP Time (sec) = 55812.195

Elapsed Time (sec) = 14373.000

Модальный анализ

CP Time (sec) = 29.859

Elapsed Time (sec) = 11.000

RS анализ

CP Time (sec) = 2.808

Elapsed Time (sec) = 4.000

Жду ваших комментариев.

Ссылка на сообщение
Поделиться на других сайтах

Присоединяйтесь к обсуждению

Вы можете опубликовать сообщение сейчас, а зарегистрироваться позже. Если у вас есть аккаунт, войдите в него для написания от своего имени.
Примечание: вашему сообщению потребуется утверждение модератора, прежде чем оно станет доступным.

Гость
Ответить в тему...

×   Вставлено в виде отформатированного текста.   Вставить в виде обычного текста

  Разрешено не более 75 эмодзи.

×   Ваша ссылка была автоматически встроена.   Отобразить как ссылку

×   Ваш предыдущий контент был восстановлен.   Очистить редактор

×   Вы не можете вставить изображения напрямую. Загрузите или вставьте изображения по ссылке.

  • Сейчас на странице   0 пользователей

    Нет пользователей, просматривающих эту страницу.




  • Сообщения

    • Горыныч
      Только Вы решаете. Я Вам из собственного опыта привел рекомендации. Если страшно, то не берите, только и всего. Сломать можно все что угодно. 
    • Orchestra2603
      Ладно. Я примерно понял, что вы имеете в виду. В теории, даже чтобы выделить n-ный базис нужно зафиксировать какие-то значения n раз. Только надо так их задавать, чтобы каждый раз получались линейно-независимые вектора. Я бы не называл это "закреплением", потому что возникает ложная кажущаяся параллель со статическим решением, где мы ищем единственное решение и для этого фиксируем какие-то компоненты, хотя здесь цель расчета принципиально иная. Это, как по мне, приводит к путанице, но если вам нравится это так называть, то - пожалуйста. Кто я такой, чтобы вам запрещать)   Я не могу сказать, что решать нужно всегда. В моем понимании в определенных итерационных схемах собственные вектора получаются полседовательными приближениями сами по себе. Применяя всякие манипуляции к матрицам (сдвиги и т.д.) можно контролировать сходимость к наибольшему, наименьшему, i-му и т.д. собственному вектору. Я не прав?   Что касается факторизаций, то для любых квадратных матриц существует разложение Шура, где ортогональными преобразованиями матрица приводится к треугольной форме. В частности, для обобщенной задачи, вот выдержка из статьи. Алгоритм получаения такого разложения для обобщенной задачи называется QZ алгоритм. Он вовсю используется в Матлабею Я не вижу никаких препятствий принципиальных, чтобы получить такую факторизацию матриц.     Т.е. похоже, что, да, все же приходится решать СЛАУ, в общем случае для треугольной матрицы A-lambda*B. И покуда lambda будет собственным значением, вся матрица A-lambda*B будет вырождена. В частности, для нулевых СЗ A-lambda*B = A, и А, сама собой, вырождена.   Но.. не могу пока сказать точно, но вроде как если A и B - симметричны, то QAZ и QBZ из треугольных превращаются в диагональные, и тогда столбцы Q становятся автоматически собственными векторами (разложение Шура переходит в спектральное), т.е. и решать ничего не приходится. Но где-то пишут, что нужна положительная определенность, где-то пишут, что только симметрии достаточно... не могу пока точно сказать.  
    • ДОБРЯК
      Я вижу, что вы не поняли. А говорили, что это знали еще в детском саду... https://ru.wikipedia.org/wiki/%D0%9C%D0%B5%D1%82%D0%BE%D0%B4_%D0%93%D0%B0%D1%83%D1%81%D1%81%D0%B0 Читайте внимательно. :=)   Ме́тод Га́усса — классический метод решения системы линейных алгебраических уравнений (СЛАУ). Назван в честь немецкого математика Карла Фридриха Гаусса. Это метод последовательного исключения переменных, когда с помощью элементарных преобразований система уравнений приводится к равносильной системе треугольного вида, из которой последовательно, начиная с последних (по номеру), находятся все переменные системы.  
    • Shura762
      (дисклеймер ) не взлетит... особенно учитывая его отношение, как человека колбасит, и почему в ТФ "узел" а не точка... перепутано местами "ты" , "на" и пропустили слово "пошёл" ... можно сказать я работаю в таком коллективе, и ЧЁ, плевали на ТФ с высокой колокольни, если молодежь можно еще как-то замотивировать, показать как надо работать, то старшее поколение делает на "отвали моя черешня" чертежи выдали тяп-ляп и что тебе надо то еще .... (дальше идиоматические выражения в рамках законодательства) p.s. я догадываюсь почему "узел" и считаю что правильное название дали.
    • lem_on
      Термообработка намного эффективнее, но из за размеров можно использовать локальное воздействие на сварочные швы вибрацией либо "проковкой" 
    • maxx2000
    • lem_on
    • Павлуха
      Да повесть на команду S и не париться.
    • gudstartup
      та и ломают не от ума! так их не ремонтируют а просто новые устанавливают
    • lem_on
      Та хз, мне ехать а не шашечки. Жалоб нет , значить всех все устраивает.    С дуру можно много чего сломать, а по цене ремонта направляющие качения нам ного больше стоят чем скольжения. 
×
×
  • Создать...