[lapkin]

Здравствуйте!
Я рассчитываю в ансис небольшой объект на сейсмическое воздействие.
Вроде бы посчитала, но есть небольшие сомнения в построенном алгоритме, так как впервые решаю подобную задачу. Прочитала литературу по переходному анализу, протестировала ряд простых примеров. Остановилась пока на упругом материале.
Большая просьба - высказать свои мысли и замечания по этому поводу. Заранее спасибо!
Еще есть вопрос: Можно ли так смоделировать вибрирующее основание?
Постановка задачи: металлоконструкция на консоле нагружена силой, связана с основанием, которое смоделировано как сосредоточенные массы большой величины, к которым приложены внешние ускорения. Во вложении алгоритм с исходными данными.
Алгоритм

Код:

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

/PREP7 
!Тип элемента
!Балочная конструкция
ET,1,BEAM188
KEYOPT,1,1,1!(1)Признак степени свободы депланации сечения при стесненном кручении, 7 степеней свободы
KEYOPT,1,2,0
KEYOPT,1,3,2
KEYOPT,1,4,2
KEYOPT,1,6,0
KEYOPT,1,7,2
KEYOPT,1,8,3
KEYOPT,1,9,3
KEYOPT,1,10,0   
KEYOPT,1,11,0   
KEYOPT,1,12,0
!*
!Основание - земля
ET,2,MASS21 
KEYOPT,2,1,0
KEYOPT,2,2,0
KEYOPT,2,3,2
!*
!Свойства материала
!Балочная конструкция
MPTEMP,1,0  
MPDATA,EX,1,,2e+11  
MPDATA,PRXY,1,,0.3  
MPDATA,DENS,1,,7850 
!*
!Основание - земля
MPTEMP,1,0  
MPDATA,EX,2,,48e+6  
MPDATA,PRXY,2,,0.35  
MPDATA,DENS,2,,2500
!*
!Постоянные
!Пояс
SECTYPE,   1, BEAM, I, h-beam, 0
SECOFFSET, CENT 
SECDATA,152e-3,152e-3,180e-3,21e-3,21e-3,13e-3,0,0,0,0  
!*
!Земля
R,1,20000, 
!*
!ПОСТРОЕНИЕ РАСЧЁТНОЙ МОДЕЛИ
base=1.5
h=2
Length=3
!
K,1,base/2,0,base/2,
K,2,-base/2,0,base/2,   
K,3,-base/2,0,-base/2,  
K,4,base/2,0,-base/2,   
K,5,base/2,h,base/2,
K,6,-base/2,h,base/2,   
K,7,-base/2,h,-base/2,  
K,8,base/2,h,-base/2,   
K,9,0,2*h,0,
K,10,0,2*h,Length,  
K,11,0,2*h,1.2*Length,  
!
FLST,3,4,3,ORDE,2   
FITEM,3,1   
FITEM,3,-4  
FLST,3,4,3,ORDE,2   
FITEM,3,1   
FITEM,3,-4  
KGEN,2,P51X, , ,0,0,0, ,0   
FLST,3,4,3,ORDE,2   
FITEM,3,5   
FITEM,3,-8  
KGEN,4,P51X, , ,0,0,0, ,0   
FLST,3,1,3,ORDE,1   
FITEM,3,9   
KGEN,5,P51X, , ,0,0,0, ,0   
FLST,3,1,3,ORDE,1   
FITEM,3,10  
KGEN,2,P51X, , ,0,0,0, ,0  
!
LSTR,      13,       6  
LSTR,      12,       5  
LSTR,      15,       8  
LSTR,      14,       7  
LSTR,      17,      16  
LSTR,      20,      19  
LSTR,      23,      18  
LSTR,      22,      21  
LSTR,      25,       9  
LSTR,      24,      28  
LSTR,      27,      29  
LSTR,      26,      30  
LSTR,      31,      10  
LSTR,      32,      11  
!
FLST,5,14,4,ORDE,2  
FITEM,5,1   
FITEM,5,-14 
CM,_Y,LINE  
LSEL, , , ,P51X 
CM,_Y1,LINE 
CMSEL,S,_Y  
CMSEL,S,_Y1 
LATT,1,1,1, , , ,1  
CMSEL,S,_Y  
CMDELE,_Y   
CMDELE,_Y1  
!*  
FLST,5,4,3,ORDE,2   
FITEM,5,1   
FITEM,5,-4  
CM,_Y,KP
KSEL, , , ,P51X 
CM,_Y1,KP   
CMSEL,S,_Y  
CMSEL,S,_Y1 
KATT,       2,       1,   2,       0
CMSEL,S,_Y  
CMDELE,_Y   
CMDELE,_Y1  
!*   
LESIZE,ALL, , ,20, ,1, , ,1,
ESIZE,0,1,  
FLST,2,4,3,ORDE,2   
FITEM,2,1   
FITEM,2,-4  
KMESH,P51X  
!*
FLST,2,14,4,ORDE,2  
FITEM,2,1   
FITEM,2,-14 
LMESH,P51X 
!
FLST,4,4,1,ORDE,4   
FITEM,4,6   
FITEM,4,89  
FITEM,4,153 
FITEM,4,173 
CP, ,ALL,P51X
!
FLST,4,4,1,ORDE,4   
FITEM,4,69  
FITEM,4,132 
FITEM,4,152 
FITEM,4,236 
CP, ,ALL,P51X  
!
FLST,4,4,1,ORDE,4   
FITEM,4,48  
FITEM,4,111 
FITEM,4,131 
FITEM,4,215 
CP, ,ALL,P51X 
!
FLST,4,4,1,ORDE,4   
FITEM,4,27  
FITEM,4,90  
FITEM,4,110 
FITEM,4,194 
CP, ,ALL,P51X
!
FLST,4,5,1,ORDE,5   
FITEM,4,174 
FITEM,4,195 
FITEM,4,216 
FITEM,4,237 
FITEM,4,257 
CP, ,ALL,P51X
!
FLST,4,2,1,ORDE,2   
FITEM,4,258 
FITEM,4,278 
CP, ,ALL,P51X
!
FLST,4,2,1,ORDE,2   
FITEM,4,2   
FITEM,4,5   
CP, ,ALL,P51X  
!
FLST,4,2,1,ORDE,2   
FITEM,4,1   
FITEM,4,26   
CP, ,ALL,P51X  
!
FLST,4,2,1,ORDE,2   
FITEM,4,4   
FITEM,4,47   
CP, ,ALL,P51X
!
FLST,4,2,1,ORDE,2   
FITEM,4,3   
FITEM,4,68   
CP, ,ALL,P51X  
!Создание таблиц
stp=3!Количество шагов
Tstp=0.01!Шаг по времени
*DIM,A_X,TABLE,stp,1,1,Time
*DIM,A_Y,TABLE,stp,1,1,Time 
*DIM,A_Z,TABLE,stp,1,1,Time
!*
!Присвоение нулевому столбцу значения времени
*CREATE,ansuitmp
*VREAD,A_X(1,0),'Time','txt',' ',IJK,stp,1,1, ,
(e12.6) 
*END
/INPUT,ansuitmp

*CREATE,ansuitmp
*VREAD,A_Y(1,0),'Time','txt',' ',IJK,stp,1,1, ,
(e12.6) 
*END
/INPUT,ansuitmp 
  
*CREATE,ansuitmp
*VREAD,A_Z(1,0),'Time','txt',' ',IJK,stp,1,1, ,
(e12.6) 
*END
/INPUT,ansuitmp
!Присвоение первому столбцу значения ускорений
*CREATE,ansuitmp
*VREAD,A_X(1,1),'a_x','txt',' ',IJK,stp,1,1, ,  
(e12.6) 
*END
/INPUT,ansuitmp

*CREATE,ansuitmp
*VREAD,A_Y(1,1),'a_y','txt',' ',IJK,stp,1,1, ,  
(e12.6) 
*END
/INPUT,ansuitmp
  
*CREATE,ansuitmp
*VREAD,A_Z(1,1),'a_z','txt',' ',IJK,stp,1,1, ,  
(e12.6) 
*END
/INPUT,ansuitmp

!*vplot,A_Z(1,0),A_Z(1,1) ! вывод графика  нагружения
!Нагрузка и решение
!Начальные условия
FLST,2,4,3,ORDE,2   
FITEM,2,1   
FITEM,2,-4  
DK,P51X, , , ,0,UX,UY,UZ, , , , 
!Нагрузка на конце балки
FLST,2,1,3,ORDE,1   
FITEM,2,11  
FK,P51X,FY,-500
ACEL,0,9.81,0,  
!
FINISH
/SOLU
*do,i,1,stp
ANTYPE,TRANS
TRNOPT,full
OUTRES,ERASE
OUTRES,ALL,ALL 
AUTOTS,1
NSUBST,1
KBC,0
NLGEOM,ON
PSTRES,0
TIME,i*Tstp
D,1,ACCX,A_X(i,1)
D,2,ACCX,A_X(i,1)
D,3,ACCX,A_X(i,1)
D,4,ACCX,A_X(i,1)
!
D,1,ACCY,A_Y(i,1)
D,2,ACCY,A_Y(i,1)
D,3,ACCY,A_Y(i,1)
D,4,ACCY,A_Y(i,1)
!
D,1,ACCZ,A_Z(i,1)
D,2,ACCZ,A_Z(i,1)
D,3,ACCZ,A_Z(i,1)
D,4,ACCZ,A_Z(i,1)
!
NEQIT,100  
Solve
*enddo
FINISH

[lapkin]

Неужели, никому не интересна эта тема?

[viking1963]

А ANSYS у Вас лицензионный ?

[lapkin]

рассматриваю задачу не в коммерческих целях, а сугубо для самообразования, освоения переходного анализа.

[Фролов Дмитрий]

для того, чтобы в конструкции не возникали паразитные колебания от нагрузки собственным весом, обычно в начале добавляют шаг с этим нагружением с командой TIMINT,OFF...затем TIMINT,ON и начинают трясти

[lapkin]

Уважаемый Дмитрий, спасибо за подсказку. Если я правильно поняла, сначала решаем задачу с собственным весом с эффектом преднапряжения, потом удаляем ускорения g, решаем задачу с ускорением основания?

[Фролов Дмитрий]

Ничего удалять не надо! Сначала загружаем собственным весом, без учета инерции, чтобы не возникало колебаний от этого нагружения. Потом включаем инерцию и начинаем трясти загруженную модель.

[lapkin]

Уважаемый Дмитрий, благодарствую за помощь! У меня к Вам маленькая просьба, посмотрите, пожалуйста, блог solution с учетом Ваших рекомендаций
/SOLU
ANTYPE,TRANS
TRNOPT,full
TIMINT,OFF
NLGEOM,ON
ACEL,0,9.81,0,
TIME, 0.001
LSWRITE,1, ! Записать, как первый шаг нагружения
*do,i,1,stp
TIMINT,ON
NLGEOM,ON
TIME,i*Tstp
ALPHAD,0
BETAD,0
D,1,ACCX,A_X(i,1)
D,2,ACCX,A_X(i,1)
D,3,ACCX,A_X(i,1)
D,4,ACCX,A_X(i,1)
!
D,1,ACCY,A_Y(i,1)
D,2,ACCY,A_Y(i,1)
D,3,ACCY,A_Y(i,1)
D,4,ACCY,A_Y(i,1)
!
D,1,ACCZ,A_Z(i,1)
D,2,ACCZ,A_Z(i,1)
D,3,ACCZ,A_Z(i,1)
D,4,ACCZ,A_Z(i,1)
!
LSWRITE,(i+1), ! Записать, как последующие шаги нагружения
*enddo
LSSOLVE,1,(stp+1),1, ! Решать, используя пошаговое нагружение с 1 по (stp+1) с шагом 1
FINISH
Буду очень благодарна

[Фролов Дмитрий]

Проверьте главное, чтобы переменные stp и Tstp давали верные промежутки времени...а так в принципе все верно.

И ещё..проверьте, чтобы TIMINT,OFF не сбивалось при LSWRITE, т.к. она зависит от команды SOLCONTROL.

[lapkin]

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