2014 dxdy logo

Научный форум dxdy

Математика, Физика, Computer Science, Machine Learning, LaTeX, Механика и Техника, Химия,
Биология и Медицина, Экономика и Финансовая Математика, Гуманитарные науки




 
 сходимость разностной схемы в расчете скольжения. ansys.
Сообщение28.06.2011, 15:34 
Добрый вечер, техники и механики.
Возникла необходимость считать диссипацию энергии колеблющейся системы в контактной паре. Пока, решил научиться просто решать контакт. Взял простую задачу брусочек по брусочку. В статике все хорошо - задача сходится. Но в дальнейшем необходимо будет смотреть отклик колеблющейся системы, то есть возбуждать импульсным воздействием, а затем проводить Фурье анализ для определения собственных частот и амплитуд различных форм колебания. Вот и пытаюсь посчитать теперь эту задачу в транзиенте, чтобы убедиться что все работает и потом уже переходить на реальную конструкцию. Но когда я считаю промежуток в 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

 
 
 [ 1 сообщение ] 


Powered by phpBB © 2000, 2002, 2005, 2007 phpBB Group