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

Оптимизация в ансис


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

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

Это код прямой задачи. Файл u2.txt

/PREP7  

!*

R=0.3

!*  

ET,1,SOLID185   

!*   

MPTEMP,,,,,,,,  

MPTEMP,1,0  

MPDATA,EX,1,,2e11   

MPDATA,PRXY,1,,0.3  





BLC4,0,0,5,1,1  

CYL4,2.5,0.5,R, , , ,1

VSBV,1,2



ESIZE,0.1

VSWEEP,ALL



VSUM,ALL

/SOL



DA,5,ALL,ALL

SFA,6,1,PRES,100 



SOLVE

FINISH  

/POST1  



*GET,VV,VOLU,3,VOLU 

FINISH

Это код оптимизации

/INPUT,'u2','txt',,, 0

/OPT

OPANL,'u2','txt',' '

OPVAR,R,DV,0.1,0.4, ,

OPVAR,VV,OBJ, , , ,

OPTYPE,RUN  

OPFRST,10,2,1,

OPEXE

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


Еще раз всем добрый день. Немного разобрался с оптимизацией и вот возник такой момент:

У меня задача состоит из нескольких анализов и чтобы передать скалярные параметры и массивы из одного в другой я их сохраняю в файл и считываю в следующем анализе. Ну в общем-то как это и в хэлпе описывается через команды PARSAV и PARRES. И вот я какую штуку заметил. Если к такой задаче применять оптимизацию, то переменная задачи на втором шаге оптимизации обнуляется ну а точнее теряется. Я проверял на простой задаче, когда шаг нагружения один и один анализ - оптимизация работает, далее добавил еще шаг и сменил анализ (внутри задачи обязательно надо сохранить и считать параметры) так вот оптимизация не пошла, т.к. обнулила параметр. Кто с таким сталкивался? Спасибо

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

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

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

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

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

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

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

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

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

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

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




  • Сообщения

    • ДОБРЯК
      Вы не слышите, то что вам говорят. :=) И поэтому у вас проблемы. И Ансис и СВ и Настран и ИСПА и ... делают численную факторизацию матрицы при определении первых собственных частот незакрепленной конструкции.  Подождите немножко, я подготовил пример на две степени свободы для незакрепленной конструкции и будем его решать.
    • Orchestra2603
      Нулевая собственная частота означает нулевое собственное значение матрицы M^-1 * K, а значит нулевок СЗ у К (или бесконечно большое СЗ для M, но такого я не встречал никогда в реальных задачах). Т.е. матрица К - вырожденная. Это равнозначные утверждения про вырожденность матрицы жесткости и нулевую частоту. Почитайте работу по ссылке. Там рассматривается случай вырожденной матрицы M и алгоритм получения вектров. Случай с вырожденной К также приводится к этому случаю простыми манипуляциями.   Работа опубликована в 1972 году! Модальник Ансис, например, работает с вырожденной матрицей жесткостью вполне нормально. Он показывает 6 первых нулевых СЗ (для свободно "парящей" системы) и показывает формы соответствующие перемещениям, как жесткого целого. Все остальные СЗ (в подавляющем большинстве случаев, за исключением случаев некорректной постановки расчета) соотвтетсуют ожидаемым частотам. Не работал в Simulation, но я не буду удивлен, если там это работавет так же. Матлаб и Маткад тоже вполне умеет находить СЗ и СВ для случая с одной вырожденной матрицей. У вас не возникает ощущения, что вы ищете проблему, где ее нет?   Да, где интересуют очень высокие формы, там возникают вопросы, связанные с корректностью выбора функции формы, размеров элемента и т.п.  В моей практике с этим не сталкивался, но, наверное, если кто-то стоячие волны, например, исследует в ультразвуковом диапазоне - там, наверное, действиетолно возникают проблемы. Но это проблема самого МКЭ, а не каких-то алгоритмов или методов отыскания СЗ и СВ.   Так, по-вашем, что такое я  прислал? Это какая-то не такая факторизация?  Я же показал, есть такие методы факторизации. У меня нет оснований полагать, что все вокруг по всему миру идиоты, этого не знают, и только я молодец. У вас есть?  
    • yunoleg
      https://drive.google.com/file/d/1o5b62TuvZQalbDJ9_aG81fc6I_Ojsg7d/view?usp=drive_link, https://drive.google.com/file/d/17bjDQq5UQ9HaZRu9QIGNjD5Iy9rEQ9XB/view?usp=drive_link, https://drive.google.com/file/d/1v5rY8NXGV0dsvizuyDwHdnQE53M0QlRi/view?usp=drive_link, https://drive.google.com/file/d/1xCjwwPhNUiLshQmCz5GVwO6nwrL1teew/view?usp=drive_link, https://drive.google.com/file/d/1biXChIeRXUPkyRq-vUAX7QTYDSAFajwG/view?usp=drive_link, https://drive.google.com/file/d/1hmQcWIwtIMVSaAKob7ET6JHxK9BiN1XO/view?usp=drive_link, https://drive.google.com/file/d/1wyzqJB9iwENMyWHD1EQgyUCw7iED8Blm/view?usp=drive_link, https://drive.google.com/file/d/1WCtff2C-_jPQniOz8QPRVTEZ9lNIaXIZ/view?usp=drive_link, https://drive.google.com/file/d/1PNeqoEomozDWPtBuB7_BMwNQLdkd7F13/view?usp=drive_link, https://drive.google.com/file/d/1zWW6U-gzl9ypvOWv1AQcW18Xb4AVRJXB/view?usp=drive_link, https://drive.google.com/file/d/1mtitt98SGu63Fr_db5HSVfE43-j6utgG/view?usp=drive_link, https://drive.google.com/file/d/1TeWtLQYDtWxOnOphgJ82brZWgapoftxT/view?usp=drive_link, https://drive.google.com/file/d/1zWk7w8BxpU_uwj-LP0H23eoR8nphLeNQ/view?usp=drive_link, https://drive.google.com/file/d/1YCZDInHYxQyLo_xQQgc0sUOcuUNrv1c9/view?usp=drive_link, https://drive.google.com/file/d/16bIqWYg8ryKRFa8LTfSpaIySo1_4VdGs/view?usp=drive_link, https://drive.google.com/file/d/1x2xqlaFCfXRIwoTxvSjxciKuJ04Ao7-0/view?usp=drive_link, https://drive.google.com/file/d/18YxxWmy2zgUzbMai9lD4XTgQeJUbGVrd/view?usp=drive_link, https://drive.google.com/file/d/16ZBXioHIXFtlYMhBQFopycwDPRdDYC_x/view?usp=drive_link, https://drive.google.com/file/d/112_xqezNot4ecDOMbmvtLN_nOnEQR_tR/view?usp=drive_link, https://drive.google.com/file/d/1Qx6eOwQqD9dcxrpeudbXJL8HQTZzQe2l/view?usp=drive_link, https://drive.google.com/file/d/1tUJ14IpIHF49EXBQ85tdrEzrWDQ0zZSK/view?usp=drive_link, https://drive.google.com/file/d/1fJarfkIcPR1B2_pMod10R8iO09HqqpEn/view?usp=drive_link, https://drive.google.com/file/d/1LsWkJXPuwBL3HlhbzehJpCOjZ0Gf6Eaq/view?usp=drive_link, https://drive.google.com/file/d/1Nm2j1JAn3tTXxznxdY9Y6bpD_ZNhjy0_/view?usp=drive_link c воздухом все в порядке. Все концевики проверены. тоже все в норме
    • Леонид124
    • karlf
      В общем, сделал образ с родного диска, накатил его на новый диск и докинул туда недостающий файл и всё поехало
    • vo.aand
      Всем привет, я провожу простой расчет напряжений в SolidWorks, решила попробовать так, чтобы получить достоверную картинку, а не только значения по формуле, так вот результирующие в программе получаются в 10^4 раз больше.  Размеры листовой детали L x B x H=6.5 x 1.4 x 0.016 м, размер выреза в детали L x B x H=2.5 x 0.7 x 0.016 м. Размеры модели в м, но при 500 000 Н напряжения получаются очень большими равными 89 670 МПа это максимум, в нейтральной зоне 67 680 МПа. Если считать через формулу то будет так σ=N/F= 500 000 Н/ (0.016 x 1.5) м^2 = 22 321 429 Па= 22,32 МПа Изображение по ссылке: https://disk.yandex.ru/i/jFPT3Yn7W8hlpA.  
    • gudstartup
      я думал вы сами пересчитывали. а про параметры серво они вам ничего не сказали? может спросить отчего медленее едет....
    • ope@gusarm.ru
      производители станка сказали какие поставить в зависимости от шага винта и редуктора. Линейно то правильно считает, значит тут все так должно быть. Если б тут была ошибка то при передвижении на метр ехало больше или меньше.
    • gudstartup
      как считали flexible feed gear а именно 2084\2085 программой?
    • ope@gusarm.ru
×
×
  • Создать...