Добрый вечер, техники и механики.
Возникла необходимость считать диссипацию энергии колеблющейся системы в контактной паре. Пока, решил научиться просто решать контакт. Взял простую задачу брусочек по брусочку. В статике все хорошо - задача сходится. Но в дальнейшем необходимо будет смотреть отклик колеблющейся системы, то есть возбуждать импульсным воздействием, а затем проводить Фурье анализ для определения собственных частот и амплитуд различных форм колебания. Вот и пытаюсь посчитать теперь эту задачу в транзиенте, чтобы убедиться что все работает и потом уже переходить на реальную конструкцию. Но когда я считаю промежуток в 1 секунду скольжения с шагом 0.01 - на 24 шаге, после смещения верхнего бруска на 4,8 мм устойчивость схемы теряется и расчет расходится. Ошибка "ERROR IN ELEMENT FORMULATION". У меня вопрос, неужели это контактный элемент эррорится? Как исключить возможность данной ошибки и добиться устойчивого сходящегося решения?
ниже привожу ЛОГ файл задачи с брикетами. заранее спасибо, тем кто откликнется)
Код:
/COM,ANSYS RELEASE 11.0SP1 UP20070830 13:56:04 06/24/2011
/input,menust,tmp,'',,,,,,,,,,,,,,,,1
/GRA,POWER
/GST,ON
/PLO,INFO,3
/GRO,CURL,ON
/CPLANE,1
/REPLOT,RESIZE
WPSTYLE,,,,,,,,0
/COM, Structural
FINISH
/PREP7
!*
ET,1,SOLID186
!*
ET,2,TARGE170
!*
ET,3,CONTA174
!*
!*
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,EX,1,,2000
MPDATA,PRXY,1,,0.3
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,EX,2,,2e12
MPDATA,PRXY,2,,0.3
BLOCK,0,4,0,1,0,0.5,
BLOCK,0,1,0,1,0.5,1,
BLOCK,3.5,4,0,1,0.5,1,
!!!! МЕШИНГ ДВИЖУЩЕЙСЯ ЧАСТИ И ДЛИННОЙ СТАТИЧЕСКОЙ ПОВЕРХНОСТИ
CM,_Y,VOLU
VSEL, , , , 1
CM,_Y1,VOLU
CMSEL,S,_Y
!*
CMSEL,S,_Y1
VATT, 1, , 1, 0
CMSEL,S,_Y
CMDELE,_Y
CMDELE,_Y1
!*
CM,_Y,VOLU
VSEL, , , , 2
CM,_Y1,VOLU
CMSEL,S,_Y
!*
CMSEL,S,_Y1
VATT, 2, , 1, 0
CMSEL,S,_Y
CMDELE,_Y
CMDELE,_Y1
!*
!!!!
ESIZE,0.1,0, ! elemetn size
FLST,5,2,6,ORDE,2
FITEM,5,1 ! volume 1 medium
FITEM,5,2 ! and volume 2 long
CM,_Y,VOLU
VSEL, , , ,P51X
CM,_Y1,VOLU
CHKMSH,'VOLU'
CMSEL,S,_Y
!*
MSHAPE,0,3d
MSHKEY,1
VMESH,_Y1
MSHKEY,0
!*
CMDELE,_Y
CMDELE,_Y1
CMDELE,_Y2
!*
!!!! МЕШИНГ МАЛЕНЬКОГО "УПОРА" - тормозного башмака=)
ESIZE,0.1,0,
MSHKEY,1
!*
CM,_Y,VOLU
VSEL, , , , 3
CM,_Y1,VOLU
CHKMSH,'VOLU'
CMSEL,S,_Y
!*
VMESH,_Y1
!*
CMDELE,_Y
CMDELE,_Y1
CMDELE,_Y2
!*
!!!!! Displacement LOADS
FLST,2,1,5,ORDE,1 ! основание
FITEM,2,1
!*
/GO
DA,P51X,ALL,0 !!! dof = 0
FLST,2,6,5,ORDE,4 ! все сотальное
FITEM,2,3
FITEM,2,-6
FITEM,2,15
FITEM,2,-16
!*
/GO
DA,P51X,ALL,0 !!! dof = 0
!*
KEYOPT,2,1,1
KEYOPT,2,2,0
KEYOPT,2,3,0
!*
KEYOPT,2,4,0
KEYOPT,2,5,0
!*
KEYOPT,3,1,0
KEYOPT,3,2,1
KEYOPT,3,4,0
KEYOPT,3,5,0
KEYOPT,3,6,0
KEYOPT,3,7,0
KEYOPT,3,8,2
KEYOPT,3,9,0
KEYOPT,3,10,0
KEYOPT,3,11,0
KEYOPT,3,12,0
!*
/MREP,EPLOT
!!!! СОЗДАНИЕ КОНТАКТНОЙ ПАРЫ
/COM, CONTACT PAIR CREATION - START
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,0.12
MAT,1
MP,EMIS,1,7.88860905221e-031
R,3
REAL,3
ET,4,170
ET,5,174
R,3,,,1.0,0.1,0,
RMORE,,,0.07,0.0,1.0,
RMORE,0,0,1.0,,1.0,0.5
RMORE,0,1.0,1.2,0.0,,1.0
KEYOPT,5,4,0
KEYOPT,5,5,0
KEYOPT,5,7,0
KEYOPT,5,8,0
KEYOPT,5,9,0
KEYOPT,5,10,2
KEYOPT,5,11,0
KEYOPT,5,12,0
KEYOPT,5,2,0
KEYOPT,4,5,0
! Generate the target surface
ASEL,S,,,2
CM,_TARGET,AREA
TYPE,4
NSLA,S,1
ESLN,S,0
ESLL,U
ESEL,U,ENAME,,188,189
NSLE,A,CT2 ! CZMESH patch (fsk qt-40109 8/2008)
ESURF
CMSEL,S,_ELEMCM
! Generate the contact surface
ASEL,S,,,7
CM,_CONTACT,AREA
TYPE,5
NSLA,S,1
ESLN,S,0
NSLE,A,CT2 ! CZMESH patch (fsk qt-40109 8/2008)
ESURF
ALLSEL
ESEL,ALL
ESEL,S,TYPE,,4
ESEL,A,TYPE,,5
ESEL,R,REAL,,3
/PSYMB,ESYS,1
/PNUM,TYPE,1
/NUM,1
EPLOT
ESEL,ALL
ESEL,S,TYPE,,4
ESEL,A,TYPE,,5
ESEL,R,REAL,,3
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
/COM, CONTACT PAIR CREATION - END
/MREP,EPLOT
/COM, CONTACT PAIR PROPERTIES - START
MP,MU,1,0.12
RMODIF,3,9,0.07
RMODIF,3,21,1.4
RMODIF,3,22,100
KEYOPT,5,2,1
/COM, CONTACT PAIR PROPERTIES - END
NROPT,UNSYM
/MREP,EPLOT
/REPLOT
FLST,2,1,5,ORDE,1
FITEM,2,8
/GO
!*
SFA,P51X,1,PRES,2 ! вертикальная "прижимающая" составляющая нагрузки на движущийся брусок
FLST,2,1,5,ORDE,1
FITEM,2,11
/GO
FLST,2,1,5,ORDE,1
FITEM,2,11
/GO
!*
SFA,P51X,1,PRES,0.6 ! Горизонтальная составляющая нагрузки на движущийся брусок
!*
ANTYPE,4
!*
TRNOPT,FULL
LUMPM,1
!*
FINISH
!!!! РЕШАЛКА
/SOL
ANTYPE,4
NLGEOM,1
DELTIM,0.01,0,0
OUTRES,ERASE
OUTRES,ALL,ALL
AUTOTS,0
NCNV,2,0,0,0,0
TIME,1
