Гармоническая сила

3 марта 2011

Всем привет!

Вопрос в следующем, как задать на некоторые узлы гармоническую силу (закон вида Fcos(wt), либо Fsin(wt)).

Есть некоторая кольцевая пластинка. На внутренний радиус действует данная гармоническая сила. Внешний радиус закреплен.

В итоге требуется получить перемещения.

Сколько не пыталась, что-то не получается ((

3 марта 2011

задаете амплитуду силы и решаете harmonic response на необходимой частоте.

3 марта 2011

задаете амплитуду силы и решаете harmonic response на необходимой частоте.

Приложила силу и задала амплитуду силы.

Затем задала Analys Type - Harmonic

А где задается частота? ((

И ещё ругается на некие Master Degree of freedom are needed for harmonic analysis

Изменено 3 марта 2011 пользователем Mosjka

4 марта 2011

Вы модель то не забыли закрепить?

Частота в LoadStepOpt - time/freq, там же и демпфирование!

4 марта 2011

пример решения классической задачи на гармонический отклик.

finish
/clear,start

! задание мех. характеристик и геометрии
/PREP7  
ET,1,BEAM189
SECTYPE,   1, BEAM, RECT, , 0   
SECOFFSET, CENT 
SECDATA,10,10,10,10,0,0,0,0,0,0 
MPTEMP,,,,,,,,  
MPTEMP,1,0  
MPDATA,EX,1,,20000  
MPDATA,PRXY,1,,0.3  
MPDATA,DENS,1,,8e-10
MPDATA,DAMP,1,,0.05 

K, ,,,, 
K, ,50,,,   

LSTR,1,2  

ESIZE,2,0,  
LMESH,1  

/eshape,1,on
EPLOT   

! задание ГУ с нагрузкой
D,1,,,,,,ALL
F,2,FY,10

! выполнение модального анализа (нужен для суперпозиции мод в гармоническом анализе)
/SOL
! 20 форм и частот, включая симметричные
ANTYPE,2
MODOPT,LANB,20  
SOLVE   
FINISH  

! диапазон частот и количество шагов
f1=0	
f2=8000
stp=20

! гармонический анализ
/SOL
ANTYPE,3
HROPT,MSUP,20, ,0   
HROUT,ON,ON,0   
LVSCALE,1,  ! обязательная команда
HARFRQ,f1,f2 ! диапазон частот
NSUBST,stp,  ! кол. расчетных точек
KBC,1   
DMPRAT,0.02,	! модальное демпфирование
SOLVE  
Finish

/SOLU   ! теперь это все нужно разложить
EXPASS,1
NUMEXP,stp,f1,f2,1 ! в диапазоне с шагом со сгущением на собственных частотах
HREXP,ALL,  
SOLVE   
FINISH  

/POST26 ! классический график
FILE,'file','rfrq','.'  
NSOL,2,2,U,Y,   
XVAR,1  
PLVAR,2,

4 марта 2011

Вы модель то не забыли закрепить?

Частота в LoadStepOpt - time/freq, там же и демпфирование!

Модель закреплена.

Да, спасибо огромное, нашла!

пример решения классической задачи на гармонический отклик.

Мне вроде как не нужен пункт

"выполнение модального анализа (нужен для суперпозиции мод в гармоническом анализе)".

А вот, что такое Master DOFs? Потому что то появляется этот пункт, то нет.

И какая разница между harmonic анализами: full, reduced?

Изменено 4 марта 2011 пользователем Mosjka

4 марта 2011

Вроде считает, но почему-то никаких перемещений нет.

Посмотрите, пожалуйста!

!* Задание элементов
/PREP7  
!*  
ET,1,LINK1  
!*  
ET,2,PLANE42
!*  

!* Задание реальных констант для LINK1 - всё равно какие
R,1,20,0,   
!*  
!*Задание свойств материалов
!*LINK1 Стержни из очень жесткого материала  
MPTEMP,,,,,,,,  
MPTEMP,1,0  
MPDATA,EX,1,,2e20   
MPDATA,PRXY,1,,0.3  
MPTEMP,,,,,,,,  
MPTEMP,1,0  
!*PLANE42 Плоские элементы с к-либо свойствами
MPDATA,EX,2,,2e11   
MPDATA,PRXY,2,,0.3  
MPTEMP,,,,,,,,  
MPTEMP,1,0  
MPDATA,DENS,2,,8000 

!* Задание параметров 
*SET,r1,0.2 
*SET,r2,0.3 
*SET,F0,2 

!* Построение геометричесакой модели
!* Два круга
FLST,2,2,8  
FITEM,2,0,0,0   
FITEM,2,r1,0,0 
CIRCLE,P51X 
FLST,2,2,8  
FITEM,2,0,0,0   
FITEM,2,r2,0,0 
CIRCLE,P51X 
!* Cоздание линий для дальнейшей разбивки
LSTR,       1,       5  
LSTR,       2,       6  
LSTR,       3,       7  
LSTR,       4,       8  
!* Создание площадей
FLST,2,4,4  
FITEM,2,6   
FITEM,2,10  
FITEM,2,2   
FITEM,2,11  
AL,P51X 
FLST,2,4,4  
FITEM,2,1   
FITEM,2,10  
FITEM,2,5   
FITEM,2,9   
AL,P51X 
FLST,2,4,4  
FITEM,2,9   
FITEM,2,8   
FITEM,2,12  
FITEM,2,4   
AL,P51X 
FLST,2,4,4  
FITEM,2,7   
FITEM,2,3   
FITEM,2,11  
FITEM,2,12  
AL,P51X 

!* Задание аттрибутов для разбивки
TYPE,   2   
MAT,       2
REAL,       1   
ESYS,       0   
SECNUM, 
!* Задание разметки для разбивки и разбивка
FLST,5,4,4,ORDE,2   
FITEM,5,9   
FITEM,5,-12 
CM,_Y,LINE  
LSEL, , , ,P51X 
CM,_Y1,LINE 
CMSEL,,_Y   
!*  
LESIZE,_Y1, , ,6, , , , ,1  
!*  
FLST,5,8,4,ORDE,2   
FITEM,5,1   
FITEM,5,-8  
CM,_Y,LINE  
LSEL, , , ,P51X 
CM,_Y1,LINE 
CMSEL,,_Y   
!*  
LESIZE,_Y1, , ,15, , , , ,1 
!*  
FLST,5,4,5,ORDE,2   
FITEM,5,1   
FITEM,5,-4  
CM,_Y,AREA  
ASEL, , , ,P51X 
CM,_Y1,AREA 
CHKMSH,'AREA'   
CMSEL,S,_Y  
!*  
AMESH,_Y1   
!*  
CMDELE,_Y   
CMDELE,_Y1  
CMDELE,_Y2  
!*  
!* Создание узла в центре координат
N, ,,,,,,,  
!* Создание жестких балок - для сохранения жесткости внутренней окружности
!*Аттрибуты
TYPE,   1   
MAT,       1
REAL,       1   
ESYS,       0   
SECNUM, 
TSHAP,LINE  
!*  Создание элементов
!*  
FLST,2,2,1  
FITEM,2,421 
FITEM,2,23  
E,P51X  
FLST,2,2,1  
FITEM,2,421 
FITEM,2,127 
E,P51X  
FLST,2,2,1  
FITEM,2,421 
FITEM,2,126 
E,P51X  
FLST,2,2,1  
FITEM,2,421 
FITEM,2,125 
E,P51X  
FLST,2,2,1  
FITEM,2,421 
FITEM,2,124 
E,P51X  
FLST,2,2,1  
FITEM,2,421 
FITEM,2,123 
E,P51X  
FLST,2,2,1  
FITEM,2,421 
FITEM,2,122 
E,P51X  
FLST,2,2,1  
FITEM,2,421 
FITEM,2,121 
E,P51X  
FLST,2,2,1  
FITEM,2,421 
FITEM,2,120 
E,P51X  
FLST,2,2,1  
FITEM,2,421 
FITEM,2,119 
E,P51X  
FLST,2,2,1  
FITEM,2,421 
FITEM,2,118 
E,P51X  
FLST,2,2,1  
FITEM,2,421 
FITEM,2,117 
E,P51X  
FLST,2,2,1  
FITEM,2,421 
FITEM,2,116 
E,P51X  
FLST,2,2,1  
FITEM,2,421 
FITEM,2,115 
E,P51X  
FLST,2,2,1  
FITEM,2,421 
FITEM,2,114 
E,P51X  
FLST,2,2,1  
FITEM,2,421 
FITEM,2,113 
E,P51X  
FLST,2,2,1  
FITEM,2,421 
FITEM,2,24  
E,P51X  
FLST,2,2,1  
FITEM,2,421 
FITEM,2,25  
E,P51X  
FLST,2,2,1  
FITEM,2,421 
FITEM,2,26  
E,P51X  
FLST,2,2,1  
FITEM,2,421 
FITEM,2,27  
E,P51X  
FLST,2,2,1  
FITEM,2,421 
FITEM,2,28  
E,P51X  
FLST,2,2,1  
FITEM,2,421 
FITEM,2,29  
E,P51X  
FLST,2,2,1  
FITEM,2,421 
FITEM,2,30  
E,P51X  
FLST,2,2,1  
FITEM,2,421 
FITEM,2,31  
E,P51X  
FLST,2,2,1  
FITEM,2,421 
FITEM,2,32  
E,P51X  
FLST,2,2,1  
FITEM,2,421 
FITEM,2,33  
E,P51X  
FLST,2,2,1  
FITEM,2,421 
FITEM,2,34  
E,P51X  
FLST,2,2,1  
FITEM,2,421 
FITEM,2,35  
E,P51X  
FLST,2,2,1  
FITEM,2,421 
FITEM,2,36  
E,P51X  
FLST,2,2,1  
FITEM,2,421 
FITEM,2,37  
E,P51X  
FLST,2,2,1  
FITEM,2,421 
FITEM,2,17  
E,P51X  
FLST,2,2,1  
FITEM,2,421 
FITEM,2,337 
E,P51X  
FLST,2,2,1  
FITEM,2,421 
FITEM,2,338 
E,P51X  
FLST,2,2,1  
FITEM,2,421 
FITEM,2,233 
E,P51X  
FLST,2,2,1  
FITEM,2,421 
FITEM,2,350 
E,P51X  
FLST,2,2,1  
FITEM,2,421 
FITEM,2,349 
E,P51X  
FLST,2,2,1  
FITEM,2,421 
FITEM,2,348 
E,P51X  
FLST,2,2,1  
FITEM,2,421 
FITEM,2,347 
E,P51X  
FLST,2,2,1  
FITEM,2,421 
FITEM,2,346 
E,P51X  
FLST,2,2,1  
FITEM,2,421 
FITEM,2,345 
E,P51X  
FLST,2,2,1  
FITEM,2,421 
FITEM,2,344 
E,P51X  
FLST,2,2,1  
FITEM,2,421 
FITEM,2,343 
E,P51X  
FLST,2,2,1  
FITEM,2,421 
FITEM,2,342 
E,P51X  
FLST,2,2,1  
FITEM,2,421 
FITEM,2,341 
E,P51X  
FLST,2,2,1  
FITEM,2,421 
FITEM,2,340 
E,P51X  
FLST,2,2,1  
FITEM,2,421 
FITEM,2,339 
E,P51X  
FLST,2,2,1  
FITEM,2,421 
FITEM,2,239 
E,P51X  
FLST,2,2,1  
FITEM,2,421 
FITEM,2,240 
E,P51X  
FLST,2,2,1  
FITEM,2,421 
FITEM,2,242 
E,P51X  
FLST,2,2,1  
FITEM,2,421 
FITEM,2,241 
E,P51X  
FLST,2,2,1  
FITEM,2,421 
FITEM,2,243 
E,P51X  
FLST,2,2,1  
FITEM,2,421 
FITEM,2,244 
E,P51X  
FLST,2,2,1  
FITEM,2,421 
FITEM,2,245 
E,P51X  
FLST,2,2,1  
FITEM,2,421 
FITEM,2,246 
E,P51X  
FLST,2,2,1  
FITEM,2,421 
FITEM,2,247 
E,P51X  
FLST,2,2,1  
FITEM,2,421 
FITEM,2,248 
E,P51X  
FLST,2,2,1  
FITEM,2,421 
FITEM,2,249 
E,P51X  
FLST,2,2,1  
FITEM,2,421 
FITEM,2,250 
E,P51X  
FLST,2,2,1  
FITEM,2,421 
FITEM,2,251 
E,P51X  
FLST,2,2,1  
FITEM,2,421 
FITEM,2,252 
E,P51X  

!* Приложение силы к центральному узлу
/SOL
FLST,2,1,1,ORDE,1   
FITEM,2,421 
!*  
/GO 
F,P51X,FY,F0

!* Закрепление внешнего радиуса по "нормальному" направлению
csys,1  
FINISH  
/PREP7  
FLST,2,60,1,ORDE,8  
FITEM,2,1   
FITEM,2,-16 
FITEM,2,128 
FITEM,2,-142
FITEM,2,218 
FITEM,2,-232
FITEM,2,323 
FITEM,2,-336
NROTAT,P51X 
FINISH  
/SOL
FLST,2,60,1,ORDE,8  
FITEM,2,1   
FITEM,2,-16 
FITEM,2,128 
FITEM,2,-142
FITEM,2,218 
FITEM,2,-232
FITEM,2,323 
FITEM,2,-336
!*  
/GO 
D,P51X, , , , , ,UX, , , , ,

!*Закрепление центарльного узла по перемещениям X
FLST,2,1,1,ORDE,1   
FITEM,2,421 
!*  
/GO 
D,P51X, , , , , ,UX, , , , ,


!* Переход в декартовы координаты и решение
csys,0  
ANTYPE,3
!*  
!*  
HROPT,FULL  
HROUT,ON
LUMPM,0 
!*  
EQSLV, ,1e-008, 
PSTRES,0
!*  
harfrq,0,3  
NSUBST,10
SOLVE   
FINISH

Изменено 4 марта 2011 пользователем Mosjka

4 марта 2011

Сложно понять что вы пытаетесь смоделировать. Не могли бы вы объяснить?

4 марта 2011

Сложно понять что вы пытаетесь смоделировать. Не могли бы вы объяснить?

Да, конечно.

Есть некоторая трубка с радиусом R1.

Данная трубка вибрируюет в цилиндрическом канале с R2, заполненном металлом/жидкостью.

Соотвественно в ansys моделирую зазор между трубкой и каналом.

Центр трубки соединен с границей нерастяжимыми стержневыми элементами (т.е. трубка колеблется как жесткое тело в этом канале). На внешней границе цилиндрического канала запрещены нормальные перемещения. К центру трубки прикладывалась сосредоточенная сила, направленная в положительном направлении оси y.

В результате необходимо получить амплитуду колебаний.

Изменено 4 марта 2011 пользователем Mosjka

4 марта 2011

Понял. Завтра на работе посмотрю.

4 марта 2011

Понял. Завтра на работе посмотрю.

спасибо огромное!)

И ещё небольшой вопрос. Может делать не harmonic analys, а таблично задавать силу? или это одно и тоже?

Изменено 4 марта 2011 пользователем Mosjka

4 марта 2011

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

4 марта 2011

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

Хм. Сила гармоническая с законом F0cos(wt), колебание системы соотвестенно будут по закону Acos(wt). Вот так вроде...

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

Ещё раз спасибо за такую консультацию)

Изменено 4 марта 2011 пользователем Mosjka

5 марта 2011

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

finish
/clear,start

/PREP7  
/nerr,,,,,,10000
! задание центрального узла с нулевой массой для приложения в нем силы
ET,1,MASS21 
R,1,0,  
n,,,,,  
e,1  

! модель зазора
ET,2,PLANE42
KEYOPT,2,1,0
KEYOPT,2,2,0
KEYOPT,2,3,0
KEYOPT,2,5,0
KEYOPT,2,6,0

type,2

MPTEMP,,,,,,,,  
MPTEMP,1,0  
MPDATA,EX,1,,2e11   
MPDATA,PRXY,1,,0.3  
MPDATA,DENS,1,,8000 


r1=0.2 
r2=0.3 
F0=2 

CYL4,0,0,r1, ,r2

ESIZE,(r2-r1)/10,0,   
MSHAPE,0,2D 
MSHKEY,0
AMESH,All

! задание через CERIG жесткости внутренего цилиндра
lsel,,,,5,8 
nsll,,1 
nsel,a,,,1  

CERIG,1,All,ALL, , , ,  

! закрепление внешнего цилиндра
lsel,,,,1,4
dl,All,,All

! амплитуда силы
F,1,FY,F0

allsel,all



! вырезка из моего предыдущего мароса

! выполнение модального анализа (нужен для суперпозиции мод в гармоническом анализе)
/SOL
! 20 форм и частот, включая симметричные
ANTYPE,2
MODOPT,LANB,20  
SOLVE   
FINISH  

! диапазон частот и количество шагов
f1=0	
f2=30000
stp=50

! гармонический анализ
/SOL
ANTYPE,3
HROPT,MSUP,20, ,0   
HROUT,ON,ON,0   
LVSCALE,1,  ! обязательная команда
HARFRQ,f1,f2 ! диапазон частот
NSUBST,stp,  ! кол. расчетных точек
KBC,1   
DMPRAT,0.02,	! модальное демпфирование
SOLVE  
Finish

/SOLU   ! теперь это все нужно разложить
EXPASS,1
NUMEXP,stp,f1,f2,1 ! в диапазоне с шагом со сгущением на собственных частотах
HREXP,ALL,  
SOLVE   
FINISH  

/POST26 ! классический график
FILE,'file','rfrq','.'  
NSOL,2,1,U,Y,   
XVAR,1  
PLVAR,2,

5 марта 2011

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

Стараюсь быть аккуратной, но ошибки же у всех случаются?)

Попыталась запустить ваш макрос...

Ansys выдал, что около 2000 тысяч warning. И затем вылетел :(( Попробую потом дома на 12 ансисе запустить. может 11 глючит.

Алгоритм дал мне мой преподаватель. И сомневаться в правильности не приходится вроде как)

Вопросы по вашему макросу:

1) CERIG - данная функция работает аналогично тому, как я раньше закрепляла балками? Но где там задается жесткость, чему она равна?

2) Внешний цилиндр Вы закрепили по всем степеням свободы. А мне требуется только по нормали. Но это я изменю.

3) Я правильно понимаю, что мне необязательно проводить модальный анализ? И затем вводить модальное демпфирование? И получать отклик? Мне же просто нужно определить амплитуду... Или я чего-то не понимаю :((

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

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

Изменено 5 марта 2011 пользователем Mosjka

5 марта 2011

В макрос вставте /nerr,,,,,,10000 и он не будет вылетать.

1) жесткость = бесконечности.

2) ок

3) это один из вариантов - разложение по собственным формам.

и литература в помощь.

5 марта 2011

Отсутствующие степени свободы соединяю поэтому столько предупреждений.

В вашей задаче задана гармоническая нагрузка. Какой тип задачи тогда нужно решить?

5 марта 2011

В макрос вставте /nerr,,,,,,10000 и он не будет вылетать.

1) жесткость = бесконечности.

2) ок

3) это один из вариантов - разложение по собственным формам.

и литература в помощь.

Что-то всё равно вышибает (( Но попробую дома. У меня тут ансис часто тупит.

Спасибо за литературу - буду изучать.

В вашей задаче задана гармоническая нагрузка. Какой тип задачи тогда нужно решить?

извините, просто в последнее время слышала много советов. И такой совет тоже был. А поскольку я с ансисом ещё на Вы, поэтому возникают такие странные вопросы((

В общем буду пробовать запускать, когда приду домой.

Изменю закрепление и поменяю часть кода в решении:

вместо

NSOL,2,1,U,Y,

XVAR,1

PLVAR,2,

_linenums:0'>! выполнение модального анализа (нужен для суперпозиции мод в гармоническом анализе)/SOL! 20 форм и частот, включая симметричныеANTYPE,2MODOPT,LANB,20SOLVEFINISH! диапазон частот и количество шаговf1=0f2=30000stp=50! гармонический анализ/SOLANTYPE,3HROPT,MSUP,20, ,0HROUT,ON,ON,0LVSCALE,1, ! обязательная командаHARFRQ,f1,f2 ! диапазон частотNSUBST,stp, ! кол. расчетных точекKBC,1DMPRAT,0.02, ! модальное демпфированиеSOLVEFinish/SOLU ! теперь это все нужно разложитьEXPASS,1NUMEXP,stp,f1,f2,1 ! в диапазоне с шагом со сгущением на собственных частотахHREXP,ALL,SOLVEFINISH/POST26 ! классический графикFILE,'file','rfrq','.'NSOL,2,1,U,Y,XVAR,1PLVAR,2,

напишу так:

f1=0f2=30000stp=50! гармонический анализ/SOLANTYPE,3HROPT,FULL  HROUT,ONLUMPM,0 EQSLV, ,1e-008, PSTRES,0LVSCALE,1, ! обязательная командаHARFRQ,f1,f2 ! диапазон частотNSUBST,stp, ! кол. расчетных точекKBC,1SOLVEFinish

И покажу перемещения. Вроде всё правильно? Главное чтобы перемещения опять не по нулям получились))

И ещё. Сколько не читала хелп, моего английского не хватает понять, зачем вот эти строчки :helpsmilie: ?

LVSCALE,1, ! обязательная команда

KBC,1

P.s. Ещё раз спасибо! Я вам так благодарна!!!! :blush:

Изменено 5 марта 2011 пользователем Mosjka

5 марта 2011

Ваш код работает.

LVSCALE - это масштабный коэффициент нагрузки. При суперпозиции мод нужно его задавать.

КВС,1 - ступенчатое приложение нагрузки. Если посмотрите на график, то при 0Гц увидете ненулевую амплитуду. При КВС,0 она будет нулевой.

5 марта 2011

Ваш код работает.

LVSCALE - это масштабный коэффициент нагрузки. При суперпозиции мод нужно его задавать.

КВС,1 - ступенчатое приложение нагрузки. Если посмотрите на график, то при 0Гц увидете ненулевую амплитуду. При КВС,0 она будет нулевой.

получилось запустить. но результат следующий:

Т.е. нет результата... Перемещения по Y - нулевые. И остальные нулевые... ((( :angel_sadangel:

/clear,start/PREP7/nerr,,,,,,10000! задание центрального узла с нулевой массой для приложения в нем силыET,1,MASS21R,1,0,n,,,,,e,1! модель зазораET,2,PLANE42KEYOPT,2,1,0KEYOPT,2,2,0KEYOPT,2,3,0KEYOPT,2,5,0KEYOPT,2,6,0type,2MPTEMP,,,,,,,,MPTEMP,1,0MPDATA,EX,1,,2e11MPDATA,PRXY,1,,0.3MPDATA,DENS,1,,8000r1=0.2r2=0.3F0=2CYL4,0,0,r1, ,r2ESIZE,(r2-r1)/10,0,MSHAPE,0,2DMSHKEY,0AMESH,All! задание через CERIG жесткости внутренего цилиндраlsel,,,,5,8nsll,,1nsel,a,,,1CERIG,1,All,ALL, , , ,allsel,all! амплитуда силыF,1,FY,F0!* Закрепление внешнего цилиндра по нормалиcsys,1 lsel,,,,1,4DL,all, ,UX,   csys,0f1=0f2=30stp=10! гармонический анализ/SOLANTYPE,3HROPT,FULLHROUT,ONLUMPM,0EQSLV, ,1e-008,PSTRES,0HARFRQ,f1,f2 ! диапазон частотNSUBST,stp, ! кол. расчетных точекKBC,1SOLVEFinish

Изменено 5 марта 2011 пользователем Mosjka

Войти

Гармоническая сила

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

Mosjka 0

Ссылка на сообщение

Поделиться на других сайтах

Влад. 6

Ссылка на сообщение

Поделиться на других сайтах

Mosjka 0

Ссылка на сообщение

Поделиться на других сайтах

es-master 0

Ссылка на сообщение

Поделиться на других сайтах

Влад. 6

Ссылка на сообщение

Поделиться на других сайтах

Mosjka 0

Ссылка на сообщение

Поделиться на других сайтах

Mosjka 0

Ссылка на сообщение

Поделиться на других сайтах

Влад. 6

Ссылка на сообщение

Поделиться на других сайтах

Mosjka 0

Ссылка на сообщение

Поделиться на других сайтах

Влад. 6

Ссылка на сообщение

Поделиться на других сайтах

Mosjka 0

Ссылка на сообщение

Поделиться на других сайтах

Влад. 6

Ссылка на сообщение

Поделиться на других сайтах

Mosjka 0

Ссылка на сообщение

Поделиться на других сайтах

Влад. 6

Ссылка на сообщение

Поделиться на других сайтах

Mosjka 0

Ссылка на сообщение

Поделиться на других сайтах

Влад. 6

Ссылка на сообщение

Поделиться на других сайтах

Влад. 6

Ссылка на сообщение

Поделиться на других сайтах

Mosjka 0

Ссылка на сообщение

Поделиться на других сайтах

Влад. 6

Ссылка на сообщение

Поделиться на других сайтах

Mosjka 0

Ссылка на сообщение

Поделиться на других сайтах

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

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

Сообщения