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

Изгиб пластины

Kultura

Автор: Kultura,
в ANSYS Mechanical APDL

 Рассказать Подписчики 1

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

Kultura

Kultura 1

Опубликовано:
Опубликовано: (изменено)

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

post-31346-0-42429100-1409574895_thumb.png

post-31346-0-40381600-1409574915_thumb.png

post-31346-0-36093400-1409574956_thumb.png

post-31346-0-85636100-1409574975_thumb.png

post-31346-0-88536900-1409574987_thumb.png

test.txt

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


soklakov

soklakov 1 475

Опубликовано:
Опубликовано:
корректные результаты только с приложенными уравнениями связи на краю помимо момента

Вы про CP?

CE пробовали?

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

Kultura 1

Опубликовано:
  • Автор
Опубликовано:

 

корректные результаты только с приложенными уравнениями связи на краю помимо момента

Вы про CP?

CE пробовали?

 

 

Да, про CP. Пробовал СЕRIG, а вот CE не пробовал

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

soklakov 1 475

Опубликовано:
Опубликовано:

А на что равняетесь-то? Дайте, пожалуйста, пластину с размерами, закреплениями и нагрузками + правильный результат.

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

Kultura 1

Опубликовано:
  • Автор
Опубликовано:

А на что равняетесь-то? Дайте, пожалуйста, пластину с размерами, закреплениями и нагрузками + правильный результат.

 

Для сложного нагружения так сходу не отвечу, домой хочется. Завтра с утра пересчитаю и вышлю

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

Fedor 1 597

Опубликовано:
Опубликовано:

"Видно, что без задания связей пластину ведёт" - Это при симметричных условиях без всяких связей ?

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

Kultura 1

Опубликовано:
  • Автор
Опубликовано: (изменено)

"Видно, что без задания связей пластину ведёт" - Это при симметричных условиях без всяких связей ?

 

да

А на что равняетесь-то? Дайте, пожалуйста, пластину с размерами, закреплениями и нагрузками + правильный результат.

 

Пластина 500х500 мм толщиной 30мм, Модуль упругости 2е11. Пластина свободная, для четвертины. Моменты Мх=1000 Н*мм, Му=2000 Н*мм. Определяем перемещения крайних точек пластины по формуле для чистого изгиба. Если считать, что одна из 4-х точек находится посередине пластины, ее перемещение - 0, то перемещения смежных точек будут u1=0,278е-3 и u2=0,556е-3 мм, а перемещение противоположного узла, соответственно, равны сумме перемещений смежных точек для случая сложного изгиба u3=u1+u2

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

_Fedor 0

Опубликовано:
Опубликовано:

" Это при симметричных условиях без всяких связей ?

 

да

"  так не бывает. Где-то напахали. Покажите разбиение и где условия отметьте... :)

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

soklakov 1 475

Опубликовано:
Опубликовано:

так не бывает.

 

Из того что удается наблюдать...

чем короче сторона, к которой приложен момент(при постоянной длине второй стороны), тем чище получающийся изгиб.

Для примера гнём линейку, как привыкли:

post-33642-0-99921300-1409662834_thumb.png

Все ок, изгиб чистый.

Попробуем погнуть "вдоль":

post-33642-0-97795200-1409662971_thumb.png

Её не то, что ведёт, вообще довольно забавно выворачивает.

Ну а в случае равенства сторон получаем седлообразную форму, как было показано выше.

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

Fedor 1 597

Опубликовано:
Опубликовано:

Приложите пару сил по разным концам там где момент и все должно быть нормально. Симметрия должна соблюдаться, коль есть исходная в условиях и нагрузках, как ни крути :)

Чтобы задрать вверх туда и сила должна быть, а при моменте ее нет, а мелочи принцип Сен-Венана сглаживает.  Или момент надо приложить с двух сторон пластины, чтобы симметрия соблюдалась :)

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

Kultura 1

Опубликовано:
  • Автор
Опубликовано: (изменено)

" Это при симметричных условиях без всяких связей ?

 

да

"  так не бывает. Где-то напахали. Покажите разбиение и где условия отметьте... :)

 

Да вроде там все просто, модель твердотельная, цп-шка закрыта:

FINISH
/CLEAR,START
/FILNAME,texp6

mx=1000
mz=2000

a=45 	!mm
b=17
l=250 	!mm
h=30.23 !mm

/prep7
MP,EX,2,200000	!MPa
MP,PRXY,2,0.3
MP,DENS,2,7800e-9 !kg/mm3

!Построение модели
kd=20
es=l/kd
block,-a/2,a*(b-1)/2,0,l,-h/2,0
block,-a/2,a*(b-1)/2,0,l,0,h/2
vglue,all

ET,1,SOLID186
ESIZE,es

allsel
type,1
mat,2
vsweep,all

/SOLU
NSEL,S,LOC,y,0  
dsym,symm,y
NSEL,S,LOC,X,a*(b-1)/2 
dsym,symm,x 
NSEL,S,,,n0
D,n0,UZ,0 

! Моменты:
!NSEL,S,LOC,y,0
!D,all,UZ,0 
NSEL,S,LOC,y,l
NSEL,U,LOC,x,a*(b-1)/2
NSEL,U,LOC,x,-a/2
!CP,1,UZ,ALL
LSEL,S,,,19
nsll
*GET,nodecount,NODE,,COUNT
F,all,FY,-Mx/h/nodecount
LSEL,S,,,2
nsll
*GET,nodecount,NODE,,COUNT
F,all,FY,Mx/h/nodecount

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

Fedor 1 597

Опубликовано:
Опубликовано:

Не считает. Только строит. 

/SOLU NSEL,S,LOC,y,0 dsym,symm,y NSEL,S,LOC,X,a*(b-1)/2 dsym,symm,x NSEL,S,,,n0 D,n0,UZ,0
   что-то тут напахано, лень разбираться
Ссылка на сообщение
Поделиться на других сайтах
Kultura

Kultura 1

Опубликовано:
  • Автор
Опубликовано:

 

Не считает. Только строит. 

 

 

/SOLU NSEL,S,LOC,y,0 dsym,symm,y NSEL,S,LOC,X,a*(b-1)/2 dsym,symm,x NSEL,S,,,n0 D,n0,UZ,0
   что-то тут напахано, лень разбираться

 

 

Извиняюсь, лишнего удалил, вот рабочий вариант:

FINISH
/CLEAR,START
/FILNAME,test

mx=1000
mz=2000

a=500 	!mm
l=500 	!mm
h=30	!mm

/prep7
MP,EX,1,200000		!MPa
MP,DENS,1,7800e-9 	!kg/mm3

!Построение модели
block,0,a/2,0,l/2,-h/2,0
block,0,a/2,0,l/2,0,h/2
vglue,all

ET,1,SOLID186
ESIZE,l/20

allsel
type,1
mat,1
vsweep,all

/SOLU
NSEL,S,LOC,y,0  
dsym,symm,y
NSEL,S,LOC,X,a/2 
dsym,symm,x 

n0=NODE(a/2,0,0)
NSEL,S,,,n0
D,n0,UZ,0 

! Моменты:
!NSEL,S,LOC,y,0
!D,all,UZ,0 
!NSEL,S,LOC,y,l/2
!CP,1,UZ,ALL
LSEL,S,loc,y,l/2
LSEL,R,loc,z,h/2
nsll
*GET,nodecount,NODE,,COUNT
F,all,FY,-Mx/h/nodecount
LSEL,S,loc,y,l/2
LSEL,R,loc,z,-h/2
nsll
*GET,nodecount,NODE,,COUNT
F,all,FY,Mx/h/nodecount

ALLSEL
SOLVE
FINISH

/POST1
PLNSOL, U,Z, 0,1.0
Ссылка на сообщение
Поделиться на других сайтах
Fedor

Fedor 1 597

Опубликовано:
Опубликовано:
NSEL,S,LOC,y,0 dsym,symm,y NSEL,S,LOC,X,a/2 dsym,symm,x
  и какой механический смысл в таких условиях требовать симметрии при несимметричном приложении условий и нагрузок  ? 
Ссылка на сообщение
Поделиться на других сайтах
Fedor

Fedor 1 597

Опубликовано:
Опубликовано:

FINISH

 

/CLEAR,START

 

/FILNAME,test

 

 

mx=1000

 

mz=2000

 

 

a=500  !mm

 

l=500  !mm

 

h=30 !mm

 

 

/prep7

 

MP,EX,1,200000  !MPa

 

MP,DENS,1,7800e-9  !kg/mm3

 

 

!Построение модели

 

block,0,a/2,0,l/2,-h/2,0

 

block,0,a/2,0,l/2,0,h/2

 

vglue,all

 

 

ET,1,SOLID186

 

ESIZE,l/20

 

 

allsel

 

type,1

 

mat,1

 

vsweep,all

 

 

/SOLU

 

NSEL,S,LOC,y,0 

 

dsym,symm,y

 

NSEL,S,LOC,X,a/2

 

!dsym,symm,x

 

 

allsel

 

DK,6, ,0, ,0,ALL, , , , , , 

 

 

DK,5, ,0, ,0,ALL, , , , , , 

 

 

 

 

n0=NODE(a/2,0,0)

 

NSEL,S,,,n0

 

!!!!!D,n0,UZ,0

 

 

! Моменты:

 

 

LSEL,S,loc,y,l/2

 

LSEL,R,loc,z,h/2

 

nsll

 

*GET,nodecount,NODE,,COUNT

 

F,all,FY,-Mx/h/nodecount

 

LSEL,S,loc,y,l/2

 

LSEL,R,loc,z,-h/2

 

nsll

 

*GET,nodecount,NODE,,COUNT

 

F,all,FY,Mx/h/nodecount

 

 

ALLSEL

 

SOLVE

 

FINISH

 

 

/POST1

 

PLNSOL, U,Z, 0,1.0

 

 

Вот так все симметрично получается...

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

Kultura 1

Опубликовано:
  • Автор
Опубликовано: (изменено)

 

NSEL,S,LOC,y,0 dsym,symm,y NSEL,S,LOC,X,a/2 dsym,symm,x
  и какой механический смысл в таких условиях требовать симметрии при несимметричном приложении условий и нагрузок  ? 

 

 

Ну а если теперь увеличить ширину пластины, например, раза в 3, прогиб останется постоянным по ширине пластины или изменится? Да и при квадратной пластине  прогиб не совпадает с теоретическим значением. Пластину ведёт

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

Fedor 1 597

Опубликовано:
Опубликовано:

n0=NODE(a/2,0,0)

 

NSEL,S,,,n0

 

D,n0,UZ,0

 

Выкиньте это и замените симметричным и ничего никуда вести не будет.

Будьте проще не мудрите и все будет хорошо :)

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

Kultura 1

Опубликовано:
  • Автор
Опубликовано:

n0=NODE(a/2,0,0)

 

NSEL,S,,,n0

 

D,n0,UZ,0

 

Выкиньте это и замените симметричным и ничего никуда вести не будет.

Будьте проще не мудрите и все будет хорошо :)

 

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

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

Fedor 1 597

Опубликовано:
Опубликовано:

Я же поправил код, чтобы все было нормально. Просто запустите и убедитесь :)

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

Kultura 1

Опубликовано:
  • Автор
Опубликовано:

Я же поправил код, чтобы все было нормально. Просто запустите и убедитесь :)

 

Да, вы правы, стало нормально не зависимо от размеров пластины. Хотя по краям значения все равно отличаются.

 

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

FINISH
/CLEAR,START
/FILNAME,test

mx=10000
mz=1000
a=5000 	!mm
l=500 	!mm
h=30	!mm

/prep7
MP,EX,1,200000		!MPa
MP,DENS,1,7800e-9 	!kg/mm3

!Построение модели
block,0,a/2,0,l/2,-h/2,0
block,0,a/2,0,l/2,0,h/2
vglue,all

ET,1,SOLID186
ESIZE,l/20
allsel
type,1
mat,1
vsweep,all

/SOLU
NSEL,S,LOC,y,0 
D,all,UY,0 
NSEL,S,LOC,X,0
D,all,UX,0 

nsel,s,loc,y,0
nsel,r,loc,z,0
nsel,r,loc,x,0
d,all,all

! Моменты:
LSEL,S,loc,y,l/2
LSEL,R,loc,z,h/2
nsll
*GET,nodecount,NODE,,COUNT
F,all,FY,-Mx/h/nodecount
LSEL,S,loc,y,l/2
LSEL,R,loc,z,-h/2
nsll
*GET,nodecount,NODE,,COUNT
F,all,FY,Mx/h/nodecount

LSEL,S,loc,x,a/2
LSEL,R,loc,z,-h/2
nsll
*GET,nodecount,NODE,,COUNT
F,all,FX,Mz/h/nodecount
LSEL,S,loc,x,a/2
LSEL,R,loc,z,h/2
nsll
*GET,nodecount,NODE,,COUNT
F,all,FX,-Mz/h/nodecount

ALLSEL
SOLVE
FINISH

/POST1
PLNSOL, U,Z, 0,1.0
Ссылка на сообщение
Поделиться на других сайтах

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

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

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

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

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

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

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

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

×
 Рассказать
Подписчики 1
Перейти к списку тем

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

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


×
×
  • Создать...