[rafitop]

Кто нибудь решал задачи связанные с механикой грунтов в Ansys?
В постаноме Друкера -Прагера?
Например Задача о стене в грунте, т.е с одной стороны котлован, с другой грунт, который соответсвенно давит на стену под действием силы тяжести и она должна отклониться.


Пробовал решить в различных постановках и перепробовал все типы контактов и сопряжений..не сходится...
!!!

Материал песок E=4e7 Па с=5000 Па Phi=32 (угол трения)
Плотность 1700 коэфф Пуассона 0,35
Элемент Plane 82

[Хворобьев]

Пробовал СергейД
Хотя это вам лучше на форум www.ansys.ru
И код конечно же...
И не сходится если то с чем (ответ и как получен).

[rafitop]

решеие не сходится..само с собой...т.е просто не решается
решается нелинейная задача соответсвенно

[Хворобьев]

Сообщение от rafitop решеие не сходится..само с собой...т.е просто не решается решается нелинейная задача соответсвенно

Ну, возможно вы где то что то не так ввели. По идее он выдает вам сообщение (соотв. нам неведомые). Обычно удобно использовать верификационные примеры из справки - для освоения. Но по этой модели я не нашел такого.
В любом случае - вам лучше обратится на форум по ансис, где с вас затребуют текстовый лог задачи.
Ибо в противном случае узнать в чем вы накосячили не представляется возможным....

[rafitop]

Код:

[Выделить все]

!*  
/PREP7  
!*  
ET,1,PLANE82
!*  
KEYOPT,1,3,2
KEYOPT,1,5,0
KEYOPT,1,6,0
!*  
!*  
MPTEMP,,,,,,,,  
MPTEMP,1,0  
MPDATA,EX,1,,40e6   
MPDATA,PRXY,1,,0.35 
TB,DP,1,,,  
TBMODIF,1,1,5000
TBMODIF,1,2,32  
TBMODIF,1,3,
MPTEMP,,,,,,,,  
MPTEMP,1,0  
MPDATA,DENS,1,,1700 
MPTEMP,,,,,,,,  
MPTEMP,1,0  
MPDATA,DENS,2,,8000 
MPTEMP,,,,,,,,  
MPTEMP,1,0  
MPDATA,EX,2,,2.1e11 
MPDATA,PRXY,2,,0.3  
K, ,,,, 
K, ,,20,,   
K, ,9.5,20,,
K, ,9.5,10,,
K, ,10,10,, 
K, ,10,40,, 
K, ,40,40,, 
K, ,40,0,,  
FLST,2,8,3  
FITEM,2,1   
FITEM,2,2   
FITEM,2,3   
FITEM,2,4   
FITEM,2,5   
FITEM,2,6   
FITEM,2,7   
FITEM,2,8   
A,P51X  
RECTNG,9.5,10,10,40,
CM,_Y,AREA  
ASEL, , , ,       1 
CM,_Y1,AREA 
CMSEL,S,_Y  
!*  
CMSEL,S,_Y1 
AATT,       1, ,   1,       0,  
CMSEL,S,_Y  
CMDELE,_Y   
CMDELE,_Y1  
!*  
CM,_Y,AREA  
ASEL, , , ,       2 
CM,_Y1,AREA 
CMSEL,S,_Y  
!*  
CMSEL,S,_Y1 
AATT,       2, ,   1,       0,  
CMSEL,S,_Y  
CMDELE,_Y   
CMDELE,_Y1  
!*  
FLST,5,12,4,ORDE,2  
FITEM,5,1   
FITEM,5,-12 
CM,_Y,LINE  
LSEL, , , ,P51X 
CM,_Y1,LINE 
CMSEL,,_Y   
!*  
LESIZE,_Y1,0.5, , , , , , ,1
!*  
MSHAPE,0,2D 
MSHKEY,0
!*  
CM,_Y,AREA  
ASEL, , , ,       1 
CM,_Y1,AREA 
CHKMSH,'AREA'   
CMSEL,S,_Y  
!*  
AMESH,_Y1   
!*  
CMDELE,_Y   
CMDELE,_Y1  
CMDELE,_Y2  
!*  
CM,_Y,AREA  
ASEL, , , ,       2 
CM,_Y1,AREA 
CHKMSH,'AREA'   
CMSEL,S,_Y  
!*  
AMESH,_Y1   
!*  
CMDELE,_Y   
CMDELE,_Y1  
CMDELE,_Y2  
!*  
FINISH  
/SOL
FLST,2,2,4,ORDE,2   
FITEM,2,1   
FITEM,2,7   
!*  
/GO 
DL,P51X, ,UX,   
FLST,2,1,4,ORDE,1   
FITEM,2,8   
!*  
/GO 
DL,P51X, ,UY,   
FINISH  
/PREP7  
!*  
!*  
CM,_NODECM,NODE 
CM,_ELEMCM,ELEM 
CM,_KPCM,KP 
CM,_LINECM,LINE 
CM,_AREACM,AREA 
CM,_VOLUCM,VOLU 
/GSAV,cwz,gsav,,temp
MP,MU,1,
MAT,1   
R,3 
REAL,3  
ET,2,169
ET,3,172
KEYOPT,3,9,0
KEYOPT,3,10,2   
R,3,
RMORE,  
RMORE,,0
RMORE,0 
! Generate the target surface   
LSEL,S,,,3  
LSEL,A,,,4  
LSEL,A,,,5  
CM,_TARGET,LINE 
TYPE,2  
NSLL,S,1
ESLN,S,0
ESURF   
CMSEL,S,_ELEMCM 
! Generate the contact surface  
LSEL,S,,,9  
LSEL,A,,,10 
LSEL,A,,,12 
CM,_CONTACT,LINE
TYPE,3  
NSLL,S,1
ESLN,S,0
ESURF   
!* Create Companion Pair - Start
R,4 
REAL,4  
ET,4,169
ET,5,172
KEYOPT,5,9,0
KEYOPT,5,10,2   
R,4,
RMORE,  
RMORE,,0
RMORE,0 
TYPE,4  
ESEL,S,TYPE,,3  
NSLE,S  
ESLN,S,0
ESURF   
TYPE,5  
ESEL,S,TYPE,,2  
NSLE,S  
ESLN,S,0
ESURF   
!* Create Companion Pair - End  
ALLSEL  
ESEL,ALL
ESEL,S,TYPE,,2  
ESEL,A,TYPE,,3  
ESEL,R,REAL,,3  
/PSYMB,ESYS,1   
/PNUM,TYPE,1
/NUM,1  
EPLOT   
ESEL,ALL
ESEL,S,TYPE,,2  
ESEL,A,TYPE,,3  
ESEL,R,REAL,,3  
ESEL,A,TYPE,,4  
ESEL,A,TYPE,,5  
CMSEL,A,_NODECM 
CMDEL,_NODECM   
CMSEL,A,_ELEMCM 
CMDEL,_ELEMCM   
CMSEL,S,_KPCM   
CMDEL,_KPCM 
CMSEL,S,_LINECM 
CMDEL,_LINECM   
CMSEL,S,_AREACM 
CMDEL,_AREACM   
CMSEL,S,_VOLUCM 
CMDEL,_VOLUCM   
/GRES,cwz,gsav  
CMDEL,_TARGET   
CMDEL,_CONTACT  
/MREP,EPLOT 
FINISH  
/SOL
ACEL,0,9.8,0,   
EPLOT   
ANTYPE,0
NLGEOM,1
DELTIM,0.01,0.001,0.5   
NEQIT,100   
TIME,1  
/REPLOT,RESIZE

Вот как-то так..

[Хворобьев]

У меня образовательная версия, лимит на вашу задачу превышен.
Попробуйте обратится на...

http://www.cadfem.ru/

http://www.emt.ru/forum/viewforum.ph...bf17b6e4a16ccc

Кстати, директор cadfem вроде бы самолично занимался задачами устойчивости откосов много лет назад.

[ASЪ]

А что нить подобное есть по скаду?

[Constantin Shashkin]

Я не знаток Ansys. Вернемся к смыслу задачи: может стенка и должна падать? И в этом случае отсутствие сходимости итераций закономерный результат? Какой глубины котлован? Какие распорные системы? А руками считать пробовали?

[rafitop]

на emt шном сайте меня убеждают, что она сходится..
надо сетку поменять..и почитать хелп....
нет стенка не падает это точно гнется но не падает)

[Constantin Shashkin]

Все равно руками посчитайте! Какой глубины котлован-то? Или это секрет?

[zai]

Сообщение от rafitop Кто нибудь решал задачи связанные с механикой грунтов в Ansys? Элемент Plane 82

Ansys осуществляет поиск решения с критериями сходимости по умолчанию, которые могут быть очень консервативными. Попробуйте загрубить их. Элемент 82 с промежуточными узлами. У него квадратичные функции формы и в нелинейных задачах вначале неплохо решить задачу с элементом без промежуточных узлов (42). Контакт между стенкой и нелинейной средой также лучше задавать после решения задачи с полной агдезией.

[rafitop]

задача с 4х узловыми элементами решается, однако при угле внутреннего трения порядка 30 градусов, при изменении его допустим до 10 она снова не сходится(

[СергейД]

коэфф трения нулевой...
под стенкой в зоне пластичности нет сгущения
нелин задачи решаются поэтапно
геостатика и лучше несколько слоев грунта последовательно удалять

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

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

реши сначала линейно
ЧИТАЙ ХЕЛП!!!!!!

[СергейД]

в физлинейной постановке у тебя отклонение 23 см!!!
так что советую сначала искать ошибки в линейном расчете
а уж потом браться за нелин

СОВЕТОВАТЬ ПО нелину не буду. разберись сначала с вышесказанным.

я уж не говорю о том, что И по ручному счету (ЭТО ПРАВИЛЬНО ТЕБЕ РЕКОМЕНДОВАЛИ...) ОНА ВПОЛНЕ МОЖЕТ ЗАВАЛИТЬСЯ

[rafitop]

Я просто задал вопрос: как лучше решать такие задачи!.
Использовать 4 или 8 узловые элементы
Каков должен быть размер элементов, так как в плаксисе сетка крупная

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

НО все же я проси ВАС несколько раз показать ВАшу сетку..
Подсказать параметры контакта, ВАС, СЕРГЕЙ, посоветовали как специалиста, а Вы же просто раскритиковали..за что вам конечно же спасибо))

Я знаю что задачи решаются в несколько Этапов и неплохо бы учитывать и поэтапную выемку грунта (можно через Кил Элемент), и сначала решить задачу о том что грунт осел, а потом учитывать стену подпорную..да...не спорю...но думаю ВЫ как и я знаете что информации и примеров решения подобных задач в Ансисе в Интернете нет...(только Сивил ФЭм )

Вот и разбираюсь!...Если вы обладаете какими-либо примерами или у вас есть ссылки на статьи или журналы.. где рассматриваются задачи о гркнтах сваях...и их взаимодействии ..я был бы ВАМ признателен.если ВЫ поделились!