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

Прочносной анализ Анизотропного тела


mgsu_engineer-pgs

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

Перед началом серьезной диссертационной работы пытаюсь численно смоделировать анизотропный материал. А именно для сосны. Для начала моделирую ситуацию по ГОСТ 16483.3-84 статический изгиб балки 20х20х300мм ,а далее буду переходить к более сложной постановке задачи с топологической и параметрической оптимизацией.

Использую solid 186,цилиндрическую координатную систему (для учета годовых колец). При расчете кремния и silicon anisotropic проблем не возникает. При задании данной матрици выдает ошибку. Подскажите пожалуйста  чем это может быть связано?

Screenshot_3.jpgМатрица податливости для сосны.

Screenshot_1.jpgошибка.

Screenshot_5.jpgграничные условия и нагрузка.Screenshot_6.jpgмесная координатная система.

Screenshot_7.jpgСетка КЭ. Solid 186.

 

 

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


Может я неправильно рассчитал коэффициенты?:wallbash:

Просто усомнился т.к. для силикона они немного другие

 

Screenshot_1.jpg

Screenshot_8.jpg

Screenshot_9.jpg

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

Screenshot_3.jpg

 

то получается, что

 

9261c17b52868ae5cb75269844c9f909.png

 

Так что таки да)

8 hours ago, mgsu_engineer-pgs said:

 

Может я неправильно рассчитал коэффициенты?

 

 

А из каких книг данные по древесине и теория, подскажите?

Изменено пользователем piden
Ссылка на сообщение
Поделиться на других сайтах
9 hours ago, mgsu_engineer-pgs said:

Может я неправильно рассчитал коэффициенты?

Вы вычислили матрицу соответствия S (compliance matrix), а требуется матрица податливости (elastic coefficient matrix). Она обратна S. 

D = S^-1

 

a091026fefded7de75700c6aa9801cda.png

 

но лучше проверьте.

1 hour ago, piden said:

и теория, подскажите?

Уже нашел)

https://engineering.ucsb.edu/~sumita/courses/Courses/ME141B/HopcroftJMEMS.pdf

http://www.kaajakari.net/~ville/research/tutorials/elasticity_tutorial.pdf

Изменено пользователем piden
Ссылка на сообщение
Поделиться на других сайтах

Спасибо огромное piden! Очень помогли.

Теория книга Ашкенези Е.А. Анизотропия древесины. стр.29,38.

Данные по древесине книга Уголев Б.Н. древесиноведение стр. 181-182.

А откуда вы скрин взяли матрици then D=      ? Очень хочется почитать))))

Прикладываю файл расчета матрици податливости, по-моему все вышло, скорее всего третий наиверный.

 

 

сосна анизотропия.xlsx

Изменено пользователем mgsu_engineer-pgs
Ссылка на сообщение
Поделиться на других сайтах
13 minutes ago, mgsu_engineer-pgs said:

А откуда вы скрин взяли матрици then D=      ? Очень хочется почитать))))

Вычислил обратную к S. Но сейчас занялся этой темой плотнее, и сразу появились сомнения относительно правильности моих вычислений) Поэтому и не привожу выкладки.

 

В литературе есть описания, как найти коэффициенты матрицы D. Но, подставляя их значения и находя обратную матрицу, у меня не получилось прийти к матрице S, которая описана в той же литературе. К, примеру, пока не понятно как получена матрица коэффициентов (5) из (4) и (6).

(Взято из https://engineering.ucsb.edu/~sumita/courses/Courses/ME141B/HopcroftJMEMS.pdf )

 

beabdc70049cb768f1b4e6e019173789.png

 

 

 

 

10 minutes ago, mgsu_engineer-pgs said:

Теория книга Ашкенези Е.А. Анизотропия древесины. стр.29,38.

Спасибо!

 

А вот метод вычисления коэфф. D из курсов абакуса.



 

5421e8c834242fd7382c1fae63707dcc.png


1ed544a0ee16b6bc1724df1b10b77c3e.png


2ab6e00ae59bb64cda15913553bacc27.png

 

 

 

 

И снова у меня не получается прийти к матрице S, которая должна иметь вид

 

2b58d762202ac7929522a3eb1e832c69.png

 

(

 

 

Изменено пользователем piden
Ссылка на сообщение
Поделиться на других сайтах

Я прикрепил файл эксель , по-

10 часов назад, piden сказал:

Я имел ввиду вот этот скрин.

 

a091026fefded7de75700c6aa9801cda.png

 

 

моему последняя верная?

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

Скажем так, она похожа на верную, но...

Лучше взять какой-нибудь известный пример и проверить.

Сам планирую взять пример из NAFEMS Benchmarks из справки SolidWorks и проверить в ансисе и абакусе...

 



 

f5cd25c3bb72850025fcae1465ec3a50.png

Ссылка на сообщение
Поделиться на других сайтах
19 minutes ago, mgsu_engineer-pgs said:

по-моему последняя верная?

Возможно. Здесь описан такой же подход. Но, повторюсь, лучше проверить)

Ссылка на сообщение
Поделиться на других сайтах
1 минуту назад, piden сказал:

Возможно. Здесь описан такой же подход. Но, повторюсь, лучше проверить)

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

Screenshot_12 — копия.jpg

Screenshot_10 — копия.jpg

Screenshot_11 — копия.jpg

Ссылка на сообщение
Поделиться на других сайтах
4 часа назад, piden сказал:

К, примеру, пока не понятно как получена матрица коэффициентов (5) из (4) и (6).

С моей стороны глупо предположить, что (6) просто решили относительно напряжений ?

(4) - это определитель, сигмы получаются по правилу Крамера. Первое вроде сходится.

 

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

http://old.exponenta.ru/soft/Mathemat/pinega/a11/a11.asp   когда-то экспериментировал с похожими матрицами , может пригодится :)

Ссылка на сообщение
Поделиться на других сайтах
53 minutes ago, Борман said:

Первое вроде сходится.

Толи я уже совсем...

Можешь показать свои выкладки?

 

У меня так:



db17f35b116796cbf5b74fcbabad1a0c.png

 

 

@Fedor , это вы все в математике набирали? КАК?!

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

http://old.exponenta.ru/soft/Mathemat/pinega/a11/a11.asp   когда-то экспериментировал с похожими матрицами , может пригодится :)

 

Пальчиками стучал по клавишам :)

Ссылка на сообщение
Поделиться на других сайтах
10 minutes ago, Fedor said:

Пальчиками стучал по клавишам :)

А как, к примеру, вставить формулу в текст и вычислить ее тут же? Или отдельно вычисляли и ответ копировали? :blink:

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

Отдельно.  И вставлял в html  frontpage насколько помню, а сейчас его вроде и не поддерживают  :)

 

https://ru.wikipedia.org/wiki/Microsoft_Expression_Web  сейчас вроде этот годится, но не пробовал ... Сейчас наверное проще просто печатать в pdf из Math , тогда не знал такой технологии...

Изменено пользователем Fedor
Ссылка на сообщение
Поделиться на других сайтах

:worthy::worthy::worthy:

 

Я сначала подумал, что и текст, и формулы прямо во front end'e набирались, потом все дело экспортировалось каким-то образом.

Чтобы даже форматировать не пришлось.

Изменено пользователем piden
Ссылка на сообщение
Поделиться на других сайтах
59 минут назад, piden сказал:

Толи я уже совсем...

Можешь показать свои выкладки?

Одна компонента косячная... наверное в книжке.

Безымянный.png

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

Из math все забирал картинками. Хотел еще и nb разместить, да искать надо, да и narod умер что-то сохранилось, а доступа по ftp не знаю

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

DET расходится...

 

Была бы статья не 2010 года, а 1910 - было бы все правильно.

 

Изменено пользователем Борман
Ссылка на сообщение
Поделиться на других сайтах

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

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

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

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

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

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

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

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

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

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




  • Сообщения

    • Fedor
      "ГОСТ Р 57700.10-2018"  Удивительный пункт 5.3.1.   А если элементы квадратичные или более того кубичные или с более сложной структурой ?    Методы основаны на математике, а на математику тоже госты собираются написать ?   :) Валидация это вам не верификация   https://ru.wikipedia.org/wiki/Валидация    :):):) 
    • Alexey8107
      Вобщем результаты. В обед перед окном еще раз прошелся по своей памяти, сравнил две старые платы(одна от взорванного PSM, вторая от ремонта с рекуперацией. Я точно помню что лет 5 назад ремонтировали SPM), оставшиеся от ремонта. Начал подозревать что они обе слишком подозрительно похожи друг на друга, по сути отличаются только номиналами трансформаторов тока и совсем мелкими деталями... И понимаю что память мне сильно изменяет! Обе платы от PSM, но разной мощности. И действительно в SPM стоит плата хоть и похожа, но отличия значительные. Так что решил сразу замахнуть PSM. Станок сразу запустился. Осталось только загадкой, почему неисправный PSM в соло вчера показывал "-". Сегодня запустил на операционном столе и вижу "1", точно так же в соло. Вобщем пока свяжусь с фануком, пока уладим административные дела, пока проплатим, попробую в нем аккуратно поковыряться.
    • aalex_b
      Добрый день. ЧПУ TNC426. Станок с поворотной головой: есть вертикальное положение и есть горизонтальное. Так же поворотный стол, ось С Вышел из строя HDD. Восстановил данные через акроникс, сделав образ с соседнего аналогичного станка. Станок поехал, все зашевелилось Но теперь не знаю, как настроить ему точки смены инструмента, трансформацию и поворот оси C. Где подсмотреть инструкцию, что измерять и как измерять, куда вносить данные, в какие параметры. Прошу совета и помощи. Никогда это не делал.
    • Jesse
      За прошедшее десятилетие (ухх, как же время летит!) уже были вопросы про локальные пластические деформации в статическом расчете и про локальные формы потери устойчивости. В обоих случаях на все эти "локальности" можно подзабить. В первом случае у меня даже шаблонная фразочка есть, которую я в отчёты вставляю, ссылаясь на ГОСТ Р 57700.10-2018. Это всё хорошо. Но что у нас с цикликой? К примеру, есть вал насоса, который в процессе работы изгибается. В валу есть шпоночные отверствия, всякие переходы вплоть до радиуса R=0.25 мм и прочая мелкая ересь, где развиваются большие напряжения, превышающие предел выносливости материала.. В то время как в основной толще материала напряжения на порядок меньше предела выносливости и там всё ок. Как быть в этой ситуации? Ведь разрушение в процессе многоцикловой усталости происходит в результате зарождения и роста микротрещин. То есть по идее в этом радиусе R=0.25 мм трещинка пойдёт и разрушит весь весь вал за условные 10^12 циклов? Или не разрушит? Если не разрушит, то как обосновать, что не разрушит? С точки зрения физики/МДТТ мы знаем, что в статике у пластичного материала в малой области превысится предел текучести, появится небольшая область течения, и на этом всё. "Пластика" дальше не пойдёт. МОжно даже просчитать и проверить для успокоения души, что пластических шарниров не будет. Но в многоцикловой усталости мы ничего расчетом проверить не можем как все эти "неприятные локальности" себя поведут за 10^12 циклов. Только кривые Веллера и т.п. Ну и результаты упругого расчета..   Какой эмпирико-нормативный опыт есть у человечества, чтоб показать, что разрушения и катастрофического роста трещины НЕ произойдёт за большое кол-во циклов, если в небольшом объеме материала по результатам упругого расчета напряжения превышают предел выносливости?
    • malvi.dp
    • Srgwell
      Библиотеки не все подключены. Там другие какие то были включены. При чем сбиваются иногда
    • Snake 60
      @clavr  Попробовал в 2024-м так и есть, не хочет перепривязывать размеры, настройки не нашел :(
    • Snake 60
    • gudstartup
      @Alexey8107 учитывая то что иногда станок все же включается то похоже на плохой контакт или высохшую емкость у нас на старом svm альфа были случайные ошибки 8 9 именно из за емкости. фанук их умудрился так замаскировать между  двумя  разъемами еле нашли,,,
    • vad0000
      DELEM предлагает такой метод Приходит партия металла, Вы в ЧПУ вводите уголок 100*100., и получаете развёртку. Изготавливаете ее, и далее гнете. Одна сторона у Вас получится 100, если конечно правильно настроен упор, а другая никогда не получится 100, либо больше, либо меньше. Cybelec для каждого материала предлагает 10 различных поправок. Т.е. согнул один раз уголок, и получил поправку на всю партию металла. Пришла новая партия, согнули ещё один уголок. Все понятно?
×
×
  • Создать...