Архитектору и проектировщику|Тепло- и звукоизоляция URSA.RU
Показать сообщение отдельно
Непрочитано 06.05.2018, 13:12 #1
Почему не работает KINH или что я делаю не правильно?
smirnovnn
 
Регистрация: 29.11.2017
Сообщений: 40

Для изучения нелинейных моделей материалов MISO, BISO, BKIN, KINH взял самую простую задачу: стальной стержень (С245) диаметром 10мм длиной 100мм жестко закрепил с одной стороны, а с другой нагрузил перемещениями величиной 15 мм, которые растягивают стержень. При анализе результатов модели MISO, BISO, BKIN дают адекватные результаты. Но модель KINH дает очень большие напряжения, на последнем шаге аж 2811 МПа, по сути в 10 раз больше чем должно быть. Ошибку найти не могу. При замене в коде KINH на MISO дает нормальный результат. В чем может быть ошибка?
Код задачи следующий:

Код:
[Выделить все]
FINI
/CLEAR,STARt
/outp,,outp
/PREP7

!Ввод геометрических данных
D=10 !мм, диаметр образца
L=100 !мм

!Задание круга
CYL4,0,0,D/2

!Разбивка круга на элементы
ET,2,200,6 
SMRT,4
ESIZ,10 ! 
AMESH,ALL

!Присваивание типов элементов

ET,1,SOLID185,,2

!Выдавливание элементов в цилиндр
TYPE,1
EXTOPT,ESIZ,80
VEXT,ALL,,,0,0,L

!Закрепление узлов опоры - жесткая заделка
NSEL,S,LOC,z,0
D,ALL,ALL

! Задание перемещений в 15мм на другом конце
NSEL,S,LOC,Z,L
CP,1,UZ,ALL
nn=ndnext(0)
D,nn,UZ,15


!Присваивание свойств для стали

EMET=206E3 !Н/мм2 модуль деформации
SIGT=245 !Н/мм2 предел текучести для стальной трубы расчет в маткаде (пониженное значение)

!Задаем упругие характеристики металлического стержня
MP,EX,1,EMET
MP,DENS,1,7850E-9
MP,NUXY,1,0.30

TB,PLASTIC,1,1,3,KINH     
TBTEMP,20.0                
TBPT,DEFI,0.0108,196    
TBPT,DEFI,0.02295,245
TBPT,DEFI,0.2,245


ALLSEL
EPLOT

/SOLU !Вход  в решатель 
ANTYPE,STATIC !Тип задачи -статическая 
NROPT,UNSYM,,on !Метод Ньютона-Рафсона с несиметричной матрицей элементов
	NLGEOM,on !Включен учет больших деформаций и перемещений
	LNSRCH,on

CNVTOL,F,,0.01,,0.01 ! Сходимость по силе
CNVTOL,U,,0.01,,0.01 ! Сходимость по перемещению
NEQIT,300
OUTRES,ALL,ALL
NLDIAG,NRRE,ON	
NSUBST,50,100000,50
	CUTCONTROL,PLSLIMIT,0.15
        CUTCONTROL,NOSHAP,1
	PIVCHECK,OFF
	TIME,1

SOLVE,,,,,nocheck

FINI
Просмотров: 4116
 
Размещение рекламы