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

Пьезоэлектрики в Static Structural


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

12 минуты назад, soklakov сказал:
30 минут назад, Pumpov сказал:

Правильно хоть ввел d_21, опуская множитель 10^(-12) ?

а почему опустили?

Когда вводил с такой степенью, соответственно получался результат в 10^12 меньший, совсем абсурдный. Рискнул предположить, что нужно убрать такую степень, т.к. это и есть характерный порядок данной величины. Значит, ошибка в другом?

15 минут назад, soklakov сказал:

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

В первом случае, в Engineering Data я как раз и выбирал ортотропный материал и задавал одинаковые константы, для простоты. И все работало правильно.

А во втором случае, если я через команду "MP" ввожу константы для различных направлений, то это разве не воспринимается программой как неизотропность?

 

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


3 минуты назад, Pumpov сказал:

А во втором случае, если я через команду "MP" ввожу константы для различных направлений, то это разве не воспринимается программой как неизотропность?

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

Вы можете возразить, что 

6 минут назад, Pumpov сказал:

И все работало правильно.

Я лишь отмечаю возможное поле для ошибки.

8 минут назад, Pumpov сказал:

Значит, ошибка в другом?

пока непонятно.

у вас есть проект, где сравниваются два кубика: один через "е", другой через "d"?

Вы довольно много информации выкладываете, но возможно в проекте удастся заметить что-то еще.

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

@Pumpov

Цитата

The [e] matrix is typically associated with the input of the anisotropic elasticity in the form of the stiffness matrix [c], and the permittivity at constant strain [εS]. The [d] matrix is associated with the input of compliance matrix   and permittivity at constant stress [εT] 

 

 

Подумал, что, может быть, стоит обратить внимание.

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

у вас есть проект, где сравниваются два кубика: один через "е", другой через "d"?

Вы довольно много информации выкладываете, но возможно в проекте удастся заметить что-то еще.

Вот этот небольшой проект:

Kubik_Piezo.wbpz

19 минут назад, soklakov сказал:
Цитата

The [e] matrix is typically associated with the input of the anisotropic elasticity in the form of the stiffness matrix [c], and the permittivity at constant strain [εS]. The [d] matrix is associated with the input of compliance matrix   and permittivity at constant stress [εT] 

 

 

Подумал, что, может быть, стоит обратить внимание.

Посмотрю...

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

@Pumpov

Цитата

The program converts a piezoelectric strain matrix [d] matrix to a piezoelectric stress matrix [e] using the elastic matrix at the first defined temperature. To specify the elastic matrix required for this conversion, use the TB,ANEL command (not the MP command).

 

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

The program converts a piezoelectric strain matrix [d] matrix to a piezoelectric stress matrix [e] using the elastic matrix at the first defined temperature. To specify the elastic matrix required for this conversion, use the TB,ANEL command (not the MP command).

Точно, сам об этом вспомнил вчера перед сном...

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

Привет всем!

Я сейчас хочу провести несколько верификационных расчетов по прогибу тонких пьезоэлектрических пластин.

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

Я хочу упростить расчет. Сначала, пришла в голову мысль рассматривать симметричную задачу - 1/4 диска, а потом и вообще подумал - решить осесимметричную задачу в 2D. Создал геометрию в виде двух тонких прямоугольников, сделал Surface Body, в дереве Geometry выбрал для пункта 2D Behavior - Axisymmetric. Но вот тогда какой тип элемента прописывать, SURF154? А у него нет степени свободы для электрического потенциала, не могу найти подходящий тип КЭ. И ещё, в DesignModeler я не объединял два тела в один Part, а почему-то контакт не создался между прямоугольничками. Что-то я запутался в трёх соснах. Помогите, пожалуйста, кто знает.

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

Но вот тогда какой тип элемента прописывать, SURF154?

planeXXX

Простой путь поиска - посмотреть справку по элементам в той задаче, которую считали в 3д. Там будет указана плоская версия таких же элементов.

6 часов назад, Pumpov сказал:

И ещё, в DesignModeler я не объединял два тела в один Part, а почему-то контакт не создался между прямоугольничками.

по дефолту контакты edge-edge не ищутся. можно сменить в наcтроqках поиска контакта.

 

6 часов назад, Pumpov сказал:

Помогите, пожалуйста, кто знает.

welcome

6 часов назад, Pumpov сказал:

SURF154? А у него нет степени свободы для электрического потенциала

он вообще для нагрузок и трехмерный, вроде.

6 часов назад, Pumpov сказал:

Что-то я запутался в трёх соснах.

стоит юзать расширялки

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

Благодарю за развёрнутый ответ по пунктам. )

Буду пробовать.

 

Получилось посчитать! Хочу теперь уточнить, корректно ли я посчитал.

Есть формула, выражающая значение кривизны мономорфного слоя (склеенные пьезодиск и пассивный диск, одного диаметра) при подаче электрического напряжения V:Формула_Кривизны мономорфного слоя.png.

Индекс "1" относится к пьезодиску, индекс "2" - к пассивному диску. У меня получилось отличие от теории в районе 4%. Диаметр дисков брал 50 мм, толщины пьезодиска и пассивного слоя - 0,6 и 0,5 мм соответственно. Возможно, тут геометрическая нелинейность не дала более точного результата.

 

Я использовал элемент plane223. Правильно ли я в командной вставке прописал для материала: et,matid,223,1001,,1 ? KEYOPT(3) поставил согласно осесимметричному случаю.

Ещё в литературе (книга А.В. Наседкина) нашёл, что для 2D-случая пьезоконстанты выглядят вот так:

Пьезомодули_2D.png

И для задания с помощью команды tbdata этих пьезомодулей надо брать номера согласно таблице: Пьезомодули_2D_вызов.png .

Только я не пойму если всего 8 констант, то как, например, появился e15  и оказался десятым?

Правильно ли, нужно вычислить e_ij как для трехмерного случая (по формуле [e] = [c][d]), а потом написать в командной вставке (после объявления типа элемента) для пьезоматериала:

tb,piez,matid

tbdata,2,e_31
tbdata,5,e_33
tbdata,8,e_31

tbdata,10,e_15

?
 

 

Ссылка на сообщение
Поделиться на других сайтах
  • 2 недели спустя...
02.03.2022 в 10:16, Pumpov сказал:

Есть формула

надеюсь, очевидно, что смотреть формулу не хочется.

02.03.2022 в 10:16, Pumpov сказал:

У меня получилось отличие от теории в районе 4%.

велоиколепно ж!

02.03.2022 в 10:16, Pumpov сказал:

что для 2D-случая пьезоконстанты выглядят вот так:

там какая-то шляпа с осями... у Ансиса неочевидное назначение осей. и про это он пишет в справке.

02.03.2022 в 10:16, Pumpov сказал:

Только я не пойму если всего 8 констант, то как, например, появился e15  и оказался десятым?

может это е16?))

вообще, елси праивльно помню, то е33 это е11... а вот что это значит)))

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

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

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

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

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

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

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

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

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

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

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




  • Сообщения

    • gudstartup
      так он так и позиционируется по m19 pmc выдает g70.6 а чпу отвечает f45.7 но ориентацию и смещение в 4077 он отрабатывает нормально шпиндель встает ровно и смена происходит хорошо. вопрос почему после ввода команды управления скоростью он все еще продолжает контролировать число импульсов между нуль метками хотя в принципе уже должен отменить позиционный контроль и просто считать обороты по 0 метке как он это делает без М19?
    • Александр1979
      SP9047 SSPA:47 ILLEGAL SIGNAL OF POSITION CODER "The relationship between the A/B phase and 1-rotation signal is incorrect (Pulse interval mismatch)." "Неправильное значение счетчика импульсов сигнала на энкодере ALPHAi. На фазах A и B энкодера за один оборот шпинделя насчитывается 4096 импульсов обратной связи. Программное обеспечение по управлению шпинделем проверяет количество импульсов на фазах A и B, соответствующее энкодеру, при каждой генерации сигнала одного оборота. Данный аварийный сигнал срабатывает, если регистрируется число импульсов, нарушающее заданный диапазон."
    • gudstartup
      после ввода 789 в edit mod надо перейти в режим ручное упр. нажать mod а потом помощь и откроется wzw.hlp  с возможностью выполнения команд управления atc. если не появляется содержимое wzw.hlp то возможно что у вас другой код хотя 789 'это стандарт для dmg в этом случае надо смотреть ваш plc раздел у меня пароль прописан в файле wzw_rett.mod
    • Anat2015
      SERVO TURRET AMP ALARM - переводить не пробовали, учитывая, что AMP - это и есть сервопривод? Что тут не ясного, язык не знаком?
    • Viktor2004
      я бы попробовал позиционировать шпиндель сигналом G70.6 Ответный сигнал, шпиндель сориентирован F45.7
    • gudstartup
      Здравствуйте уважаемые участники форума попрошу вас поделится мнениями. Имеется фрезерный шпиндель с позиционным кодером (зубчатое колесо и считывающая головка). Ошибку 9047 выдает после того отрабатывает M19 c последующим включением оборотов Sxxx M3 причем обороты могут быть любыми от 100-макс. Если задавать обороты без предварительной ориентации то 9047 не возникает и шпиндель нормально вращается на любой скорости при этом обороты не гуляют точность поддержания +-1 об\мин. Кабель заменен заменить сам датчик нет возможности. Мой вопрос в следующем почему шпиндель после M19 и включении управления скоростью Sxxx M3 остается в режиме позиционирования а не сразу переходит врежим управления скоростью и даже явная отмена ориентации кодом M20 ничего не меняет ведь режим позиционного управления должен быть отменен или он остается и отменяется не сразу после подачи команды M3 и кода скорости? Станок достаточно старый 2009 года плату управления шпинделем меняли на соседний ничего не изменилось. Это может быть логика PMC или скорее всего это датчик?   PMC-SB.001
    • Killerchik
      А мы-то думали, что пропорциональный гидравлический.....
    • Grig75
    • EnKlL
      Я написал полностью ошибку, серво привод
    • gudstartup
      https://www.aggsoft.com/cnc-dnc/maho-cnc-432-1000.htm
×
×
  • Создать...