[sookin]

Уважаемые участники данного форума, прошу опять ткнуть меня носом где я не прав!
Озадачился я вопросом нахождения критической силы при потере устойчивости по изгибно-крутильной форме (LTB). Определился с габаритами двутавра, задался мощным моментов в силовой плоскости и мааааленьким моментиком из плоскости. Закрепления шарнирное; концевые сечения не могут крутиться, но могут депланировать
Решил задачу тремя путями:
1) beam188 с использованием SECOFFSET; получил критическую силу что-то вроде 570 кН
2) beam188 с использованием моментов и центрального сжатия; получил критическую силу что-то около 960кН
3) shell181 с использованием моментов и центрального сжатия; получил критическую силу примерно 960 кН


Формы потери устойчивости во 2-м и 3-м случае бьют между собой и если верить статистике и внутреннему чутью, то этот вариант верный.
Но! тогда я не понимаю смысла в SECOFFSET, потому что решая такую же по смыслу задачу, но использованием швеллера, верная критическая сила нашлась при использовании SECOFFSET, а не варианта с моментами (но у швелллера центр тяжести не совпадает с центром изгиба, а у двутавра совпадает)
Отсюда вопрос - какая хитрость имеется при решении задач об LTB с заданием эксцентриситетов с помощью SECOFFSET?
Коды первых двух вариантов прикладываю для облегчения взаимопонимания

FINISH
/CLEAR
/FILNAME, LTB_secoffset,1
/TITLE, A1
tf = 0.014
tw = 0.01
bf = 0.32
H = 0.5
L = 10.08
E = 2.06e11
m = 0.3
or_k = H/2 ! Точка ориентации сечения
P = 1 ! Значение силы (н)
ey = 0.25 ! Эксцентриситет (м) в плоскости действия момента
ex = 0.017179 ! Эксцентриситет (м) из плоскости действия момента
/PREP7
ET,1,BEAM188
KEYOPT,1,1,1
KEYOPT,1,2,1
KEYOPT,1,4,2
sectype,1,beam,I,SCT1,2
secdata,bf,bf,H,tf,tf,tw
SECOFFSET,USER,-ex0,(ey+H/2), , , ,
MP,EX,1,E
MP,PRXY,1,m
K,1,0,0,0
K,2,0,0,L
k,3,0,or_k,0
L,1,2 !
LSEL,S, LINE, ,1
LATT,1, , 1, ,3, ,1
LESIZE,ALL,0.05,,,,1,,,1
LSEL,ALL
LMESH,ALL
/SOLU
ANTYPE, STATIC
PSTRES, ON
DK,1, , , ,0,UX,UY,UZ,ROTZ , , ,
DK,2, , , ,0,UX,UY,ROTZ , , , ,
FK,2,FZ, -P
SOLVE
FINISH
/SOLU
ANTYPE, BUCKLE
PSTRES,ON
BUCOPT,LANB,1
SOLVE
FINISH
/SOLU
EXPASS, ON
MXPAND, 1,,, 0
SOLVE
FINISH
/POST1
SET, FIRST
PLDISP, 1

FINISH
/CLEAR
/FILNAME, LTB_moments,1
/TITLE, A1
tf = 0.014 !Толщина полок (м)
tw = 0.01 !Толщина стенки (м)
bf = 0.32 ! Ширина полок (м)
H = 0.5 ! Высота сечения (м)
L = 10.08 ! Длина элемента (м)
E = 2.06e11 !Модуль упругости (Па)
m = 0.3 ! Коэффициент Пуассона
or_k = H / 2 ! Точка ориентации сечения
ey = 0.25
ex = 0.017179
P = 1 ! Значение силы (н)
My = P*ex
Mx = P*ey
/PREP7
ET,1,BEAM188
KEYOPT,1,1,1
KEYOPT,1,2,1
KEYOPT,1,4,2
sectype,1,beam,I,SCT1,2
secdata,bf,bf,H,tf,tf,tw
SECOFFSET,CENT,,, , , ,
MP,EX,1,E
MP,PRXY,1,m
K,1,0,0,0
K,2,0,0,L
k,3,0,or_k,0
L,1,2
LSEL,S, LINE, ,1
LATT,1, , 1, ,3, ,1
LESIZE,ALL,0.05,,,,1,,,1
LSEL,ALL
LMESH,ALL
/SOLU
ANTYPE, STATIC
PSTRES, ON
DK,1, , , ,0,UX,UY,UZ,ROTZ , , ,
DK,2, , , ,0,UX,UY,ROTZ , , , ,
FK,1,MX,Mx
FK,1,MY,My
FK,2,FZ, -P
FK,2,MX,-Mx
FK,2,MY,-My
SOLVE
FINISH
/SOLU
ANTYPE, BUCKLE
PSTRES,ON
BUCOPT,LANB,1
SOLVE
FINISH
/SOLU
EXPASS, ON
MXPAND, 1,,, 0
SOLVE
FINISH
/POST1
SET, FIRST
PLDISP, 1

[MrWhite]

Ничего не понял... Зачем вы задаете SECOFFSET и одновременно эксцентриситеты и еще и пытаетесь сравнить с шелами?

[sookin]

Ничего не понял... Зачем вы задаете SECOFFSET и одновременно эксцентриситеты и еще и пытаетесь сравнить с шелами?

три разные задачи:
1-я: сила приложена через SECOFFSET (просто в коде переменные для secoffset названы как эксцентриситеты в двух плоскостях)
2-я: приложена сила плюс моменты (соответственно M=P*e в двух плоскостях)
3-я: исходные данные второй задачи, но двутавр сделан из оболочек (в этой задаче как раз многое прозрачно и она не вызывает много вопросов)

[MrWhite]

Прошу прощения, не разобрался.