[derik]

Здраствуйте уважаемые форумчане! Просьба тех кто считает в Ансисе подсказать - где у меня ошибка в схеме! Хелп перечитал, он мне уже скоро сниться будет :-( помогите пожалуйста!
вот код:

Код:

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

/prep7
a=4
h=2.5
b=.4
t=.2
q=100
/view,1,1,1,1
/ang,1
/rep,fast
blc5,a/2,a/2,a,a,
k,,b/2,b/2,0
k,,a-b/2,b/2,0
k,,b/2,a-b/2,0
k,,a-b/2,a-b/2,0
k,,b/2,b/2,-h
k,,a-b/2,b/2,-h
k,,b/2,a-b/2,-h
k,,a-b/2,a-b/2,-h
lstr,5,9
lstr,6,10
lstr,7,11
lstr,8,12
nummrg,all
boptn,,keep
boptn,keep,1
boptn,nwar,0
boptn,vers,rv52
btol,1e-005,
lina,5,1
lina,6,1
lina,8,1
lina,7,1
mptemp,,,,,,,,
mptemp,1,0
mpdata,ex,1,,3e8
mpdata,prxy,1,,.3
mpdata,dens,1,,2500
et,1,shell63
r,1,t,,,,,,
et,2,beam4
keyopt,2,9,9
r,2,b*b,(b*b**3)/12,(b*b**3)/12,b,b,
lsel,s,,,5
lsel,a,,,6
lsel,a,,,7
lsel,a,,,8
latt,1,2,2,
asel,s,,,1
aatt,1,1,1
allsel,all
lesize,5,,,b,1
lesize,6,,,b,1
lesize,7,,,b,1
lesize,8,,,b,1
lmesh,5,6,
lmesh,7,8,
aesize,1,b/2,
amesh,1
LGEN,3,all, , , , ,h, ,0   
AGEN,3,all, , , , ,h, ,0
nummrg,all
allsel,all
/eshape,1
/replot
finish
/sol
nsel,s,loc,x,b/2
nsel,r,loc,y,b/2
nsel,r,loc,z,-h
D,all, , , , , ,ALL, , , , , 
nsel,s,loc,x,b/2
nsel,r,loc,y,a-b/2
nsel,r,loc,z,-h
D,all, , , , , ,ALL, , , , , 
nsel,s,loc,x,a-b/2
nsel,r,loc,y,b/2
nsel,r,loc,z,-h
D,all, , , , , ,ALL, , , , , 
nsel,s,loc,x,a-b/2
nsel,r,loc,y,a-b/2
nsel,r,loc,z,-h
D,all, , , , , ,ALL, , , , ,
allsel,all
acel,0,0,9.8,
solve
finish
/post1
set,first

[СергейД]

перед решением
solve
поставь allsel

[derik]

Спасибо!
Все равно пишет, что 18 узел превысел максимальное перемещение заложенное в программе.
Какой командой нужно "склеивать" стержни и пластины? Или, возможно, я что-то не так делаю?

[СергейД]

нет рядом ансиса...
значит, не соединены части модели
поставь
nummrg,all
и посмотри на диагностику.
или лучше для начала нарисуй узлы с номерами
где будет двоиться, там и разбирайся.

[derik]

Спасибо Сергей! Эта команда помогла.
Чувствовал что нужно "склеить" модель, но не знад как и этой команды нигде не встречал! спасибо еще раз, штурмуем Ансис дальше :-D

[СергейД]

склеить нужно бы линии до генерации узлов,
ИНАЧЕ могут быть несклеены узлы в середине линий.
совет- все пересечь, lptn,all
а уж потом присвыаивать атрибуты и генеировать узля.

[derik]

Попробую как Вы говорите - все пересечь, но и так посчиталось. Правда я код в первом посте вставил не тот но это уже мелочи.

Сегрей, а в Ансис, для моделирования ригелей с помощью "жестких вставок" как в лире, используются элементы типа conta171-178 и targe169-170 или есть иные способы?

[derik]

Код:

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

/prep7
a=4
b=6
c=.4
t=.2
blc5,a/2,b/2,a,b
k,,c/2,0,-(t/2+c/2)
k,,a-c/2,0,-(t/2+c/2)
k,,a-c/2,b,-(t/2+c/2)
k,,c/2,b,-(t/2+c/2)
lstr,5,8
lstr,6,7
mptemp,,,,,,,,
mptemp,1,0
mpdata,ex,1,,3e8
mpdata,prxy,1,,.3
mpdata,dens,1,,2500
et,1,beam4
r,1,c*c,(c**3)/6,(c**3)/6,c,c,0
et,2,shell63
r,2,t,,,,,,
lsel,s,,,5
lsel,a,,,6
latt,1,1,1,
allsel,all
/replot
asel,s,,,1
aatt,1,2,2,
allsel,all
lesize,5,,,c/2,1
lesize,6,,,c/2,1
lmesh,5,6,
aesize,1,c/2,
amesh,1,
allsel,all
/eshape,1
/replot
finish
/sol
nsel,s,loc,x,c/2
nsel,r,loc,y,0
nsel,r,loc,z,-(c/2+t/2)
d,all,,,,,,all,,,,,
nsel,s,loc,x,a-c/2
nsel,r,loc,y,0
nsel,r,loc,z,-(c/2+t/2)
d,all,,,,,,all,,,,,
nsel,s,loc,x,c/2
nsel,r,loc,y,b
nsel,r,loc,z,-(c/2+t/2)
d,all,,,,,,all,,,,,
nsel,s,loc,x,a-c/2
nsel,r,loc,y,b
nsel,r,loc,z,-(c/2+t/2)
d,all,,,,,,all,,,,,
allsel,all
/replot
acel,0,0,9.8,
finish

как мне соединить плиту и балки, какими элементами или командами, как мне учесть эксцентриситет балки? элементы cont и targe это нелин., как я понял из хелпа, до него еще не дорос поделитесь пжалуйста опытом :-)

[aka raduntsev-nias]

6pa4o, с рсу как то решил вопрос?

[derik]

Мне до РСУ как ежу до Пекина я пока учусь моделировать и задавть свойства материалов и нагрузки на простых примерах - рамы, каркасы, фермы. В своем предыдущем вопросе разобрался, понял почему не решается, спасибо СергеюД! Вот дошел до ребристой плиты и стопарнулся, не пойму как задать эксцентриситет ребра плиты:
1 - вроде нужно через conta и targe, но не пойму как;
2 - через задание secdata, но там, наверно, не учитывается контакт ребра и плиты.
Считать солидами тоже пока рано, надо со стержнями и пластинами разобраться

[СергейД]

secoffs
secdata

связь произойдет из общности узлов.
почитай теперь хелп хотя бы пару дней. многое прояснится.
все эти вставки= абстракция, сложная и ненужная по сути вещь
(но лучше всего MPC184)

а солиды чем раньше, тем лучше.
набери балку балками, потом шеллами, потом солидами.
сравни.
и работай солидами.

[aka raduntsev-nias]

СергейД, но ведь в солиды еще арматуру придется вставлять(солидами). но как узнать, какая нужна до расчета?

без шуток, для меня пока это секрет.

[СергейД]

линейно решайте 190 или 186 или 45 типами. ничего туда вставлять не обязательно.

а подбор арматуры не нужно усложнять....
поставил по опыту в запас. решил нелинейно и убрал, где лишняя (если хочется)
успехов.

[СергейД]

линейно решайте 190 или 186 или 45 типами. ничего туда вставлять не обязательно.

а подбор арматуры не нужно усложнять....
поставил по опыту в запас. решил нелинейно и убрал, где лишняя (если хочется)
успехов.

[СергейД]

линейно решайте 190 или 186 или 45 типами. ничего туда вставлять не обязательно.

а подбор арматуры не нужно усложнять....
поставил по опыту в запас. решил нелинейно и убрал, где лишняя (если хочется)
успехов.

[СергейД]

линейно решайте 190 или 186 или 45 типами. ничего туда вставлять не обязательно.

а подбор арматуры не нужно усложнять....
поставил по опыту в запас. решил нелинейно и убрал, где лишняя (если хочется)
успехов.

[СергейД]

==============

[derik]

Спасибо СергейД! я вот склонялся к 184, но решил его не трогать и переключился на 170 и 175, показалось что будет проще, а оказывается проще - моделировать солидами буду изучать!

[aka raduntsev-nias]

согласен, что проще. есть для этого адванс, ревит и тд. сложнее в cad системе нарисовать плоскими объемную схему и сделать импорт в ансис. хотя плоскими привычнее. но сравнить результат интересно. не думаю, что будет сильно отличаться.

[МишаИнженер]

Ребята честно Вам скажу: я с вас тащюсь (еле хожу). Как Вы можете считать конструкции на основе каких-то записей в блокноте? Неужели нельзя предложить файл который может нарисовать что-то похожее на конструкцию на экране. Мне просто за Вас обидно. Честное слово.

[aka raduntsev-nias]

МишаИнженер, главное научиться. и сделать расчет. а там решим и посмотрим. и сравним.

[derik]

Offtop: А ты попробуй! меня, лично, очень сильно наколяло в лире и, особенно, в скаде стотыщьпятьсотиконок, а тут свобода действий! Нужно, к примеру, изменить все данные конструкции: высота этажа, размеры в плане, сечение колонн и тд. в Ансисе просто меняешь значения параметров конструкций и все! а в лире или старке все нужно заново отрисовывать!

[aka raduntsev-nias]

Offtop: ансис - для людей! свободных и достойных.

[derik]

Offtop: +1

[aka raduntsev-nias]

Offtop: и не для тупых карандашников.

[Кулик Алексей aka kpblc]

C флудом завязывайте.

[derik]

Код:

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

/prep7
a=8
b=3.5
c=.3
h=.2
k=c+h
q=5
et,1,solid186
mptemp,,,,,,,,
mptemp,1,0
mpdata,ex,1,,3e8
mpdata,prxy,1,,.3
mpdata,dens,1,,2500
blc5,c/2,c/2,c,c
blc5,c+h/2,b/2,h,b
agen,2,1,,,,b-c,,,0
aadd,1,2
aadd,4,3
vext,all,,,0,0,a,,,,
vsweep,1
/view,1,1,1,1
/ang,1
/rep,fast
finish
/sol
asel,s,,,1
da,all,all
asel,s,,,2
da,all,all
acel,9.8,0,0
sfa,8,1,pres,q
allsel,all
/replot
solve
finish
/post1
set,first

все это набрал за 30 мин, разбирался с солидом. все перехожу на солиды - быстрее и меньше телодвижений :-)

[aka raduntsev-nias]

дак такую штуковину еще быстрее сделать в 3д потом закинуть в ансис.

[An2]

какую книжку для начала почитать, посоветуйте пожалуйста.

[СергейД]

браво, bra40! ты на верном пути.
raduncev! а если изменится пара размера, снова станешь перерисовывать и импортировать?
чем сложнее объект, тем его тяжелее импортировать. уж послушайте опытного человека. и на импортированном хорошую сетку создать сложнее. (импорт = 10% расчета разве что)
проще бывает перечертить в ансис. да я уж не раз писал это на форуме. повторять неловко.

2 an2
уж не раз эту тему поднимали. поищи по форуму советы+ в даунлоаде есть огромная подборка.
я проводил курсы в МГСУ. и осенью (не раньше) буду.

[derik]

Спасибо, Сергей, стараюсь а моделировать в среде Ансис проще и быстрее чем в CAD программах...жалко кнопки "назад" нет как будет время выложу Басова "Справочник пользователя" 2005г. почти хелп на русском.

[derik]

Код:

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

finish
~cfclear,nostart
/prep7
a=5
b=.4
c=.3
h=3.2
q=1500
et,1,solid186
mptemp,,,,,,,,
mptemp,1,0
mpdata,ex,1,,.3e9
mpdata,prxy,1,,.2
mpdata,dens,1,,2200
blc5,b/2,b/2,b,b,-h
vgen,3,1,,,a,,,,0
vgen,5,all,,,,a,,,0
/view,1,1,1,1,
/ang,1
/rep,fast
blc5,b/2,(4*a)/2,b,4*a,c
vgen,3,16,,,a,,,,0
blc5,(2*a+b)/2,b/2,2*a+b,b,c
vgen,5,19,,,,a,,,,0
blc5,(2*a+b)/2,(4*a+b)/2,2*a+b,4*a+b,c
vgen,,24,,,,,c,,,1
blc5,(2*a+b)*3/4,(4*a+b)*7/8,a/2,a/2.5,c
vgen,,25,,,,,c,,,1
vsbv,24,25
vovlap,all
smrtsize,on
mshape,1,3d
mshkey,0
desize,b,b,b,,,,,,
vsweep,all
vsel,s,,,58
mshkey,0
aesize,all,b,
vmesh,58
allsel,all
/replot
finish
/sol
asel,s,loc,z,-h
da,all,all
acel,0,0,9.8
allsel,all
solve
finish
/post1
set,first
plnsol,u,sum,0,1

Вот. Жалко на Ансис времени мало остается

[МишаИнженер]

Что это такое?
А можно в ANSYSe начертить балку с поясами разной толщины и стенкой с переменной толщиной, назначить сварные швы и посчитать их?

[derik]

Можно. а вот как это сделать я, по крайней мере, не знаю
код над твоим постом это рассчитывает такой каркас (выдуманный мной )

[aka raduntsev-nias]

Сообщение от 6pa4o Можно. а вот как это сделать я, по крайней мере, не знаю код над твоим постом это рассчитывает такой каркас (выдуманный мной )

еще бы рсу туда. но рсн возможно. по крайней мере в сивилфеме.

[derik]

ручками вводить комбинации, а кому сейчас легко? тем более со временем это будет делаться пару секунд! Хотя СергейД, скажет более точно.

[СергейД]

Миша! (насчет Инженер, не знаю...) твоя задачка элементарно решается с помощью "блокнота", над которым ты (не понимая, о чем речь) "тащился".
про лисп слышал? это типа того... только проще и удобнее. именно этот файл и генерирует модель и рисунок. параметрически.
разумеется, есть и обычное рисование...

бра40! UNdo есть. его просто нужно заранее заказать
/undo,on
но есть куда лучше.
session editor найди в меню задач.
открывается блокнот с командами, выполненными со времени последнего SAve.
и можно исправить любые команды из них.
исправляешь радиус (например) и все операции повторяются заново (пока они исполнимы)
если считаешь, что в ансис чего-то нет, ты сначала спроси...
или лучше почитай на форуме мои прежние ответы.
возможно, что-то пригодится. мне, согласись, повторять неловко.

РСН с трудом подвластно логике. я отставил это все.
НИИМосстрой составил макрос РСН и сертифицировал его в Госстрое....

[derik]

session editor,ом пользуюсь а вот про undo, не знал, спасибо! Сергей, а где этот макрос можно получить? пока он не нужен - хочу, сначала, уверенно моделировать конструкции и "разбивать" их. но макрос пусть будет, кушать не просит ведь!

[aka raduntsev-nias]

точно. и мне макрос нужен. могу даже в ответ что-нить подарить.

[derik]

а что у тебя есть, может быть и мне пригодиться?

[aka raduntsev-nias]

смею тебя заверить много чего, что будет интересно

[derik]

огласить весь список пжалуйста
З.Ы. я тоже кой что накопал

[master_luc]

Давно для определения РСУ в комплексе Femap/Nastran применяю интересную программу - комбинатор.
Суть - РСУ строяться в виде графа, результаты берутся из текстового файла, потом происходит расчет комбинаций усилий. Результаты расчета импортирую обратно также через текстовый файл.
Посмотрите сайт автора программы http://sites.google.com/site/40ushek/.
Думаю, что ее свободно можно адаптировать под Ansys.
А если по богатому - то CivilFem, к сожалению реализации РСУ там нет, только РСН, которые пользователь определеяет сам (в последнем не уверен).

[СергейД]

макрос РСН в ниимосстрое... думаю, так просто не дадут.

[derik]

aka raduntsev-nias, а ты упругое основание не пробывал моделировать в ANSYS?

[aka raduntsev-nias]

Сообщение от 6pa4o aka raduntsev-nias, а ты упругое основание не пробывал моделировать в ANSYS?

нет, 6pa4o, к сожалению нет. это пока в мечтах. ansys для меня пока очень сложен и труден.

[derik]

пнятно, а что там насчет списка литературы, которую ты насобирал?

[aka raduntsev-nias]

есть и потихоньку читается. халтуру взял для пополнения финансов. иначе мне быстро прокручу, потом продолжу. поэтому на время исчезаю.

[IMG]http://s003.***********/i202/1005/c4/dd2228da0786t.jpg[/IMG]

[derik]

а можно всю книжку посмотреть А по поводу РСН: либо макрос написать, либо, для начала, рассчитать все РСН вручную и каждую комбинацию прикладывать к модели, записывая рассчитанные нагружения через LSWRITE, апотом запустить на пошаговый расчет. СергейД, такой вариант имеет право на жизнь?

[aka raduntsev-nias]

это Перельмутер. "Расч модели сооруж и возм их анализа_2002"
в дауне дб

[СергейД]

именно посмотрев справку по скад, я и понял, что у меня есть более важные и более
легкие дела.
была ситуация. посмотрели коэфф в скад.
задали их в ансис lcoper. все.

[derik]

Вот, вроде маленькая модель, а считает долго.

Код:

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

finish
~cfclear,nostart
/prep7
a=5
b=.4
c=.3
h=3.2
q=1500
et,1,solid186
mptemp,,,,,,,,
mptemp,1,0
mpdata,ex,1,,.3e9
mpdata,prxy,1,,.2
mpdata,dens,1,,2200
blc5,b/2,b/2,b,b,-h
vgen,3,1,,,a,,,,0
vgen,5,all,,,,a,,,0
/view,1,1,1,1,
/ang,1
/rep,fast
blc5,b/2,(4*a)/2,b,4*a,c
vgen,3,16,,,a,,,,0
blc5,(2*a+b)/2,b/2,2*a+b,b,c
vgen,5,19,,,,a,,,,0
blc5,(2*a+b)/2,(4*a+b)/2,2*a+b,4*a+b,c
vgen,,24,,,,,c,,,1
blc5,(2*a+b)*3/4,(4*a+b)*7/8,a/2,a/2.5,c
vgen,,25,,,,,c,,,1
vsbv,24,25
vgen,3,all,,,,,2*c+h,,,0
vovlap,all
smrtsize,on
mshape,1,3d
mshkey,0
desize,b,b,b,,,,,,
vsweep,all
vsel,s,,,194
vsel,a,,,195
vsel,a,,,196
aesize,all,b,
vmesh,all
allsel,all
/replot
finish
/sol
asel,s,loc,z,-h
da,all,all
acel,0,0,9.8
allsel,all
solve
finish
/post1
set,first
plnsol,u,sum,0,1

абидно

[Дмитрий Дмитриев]

Бра4о! Поменяйте решатель с прямого SPARSE(по умолчанию) на итерационный, к примеру PCG(см. команду EQSLV)

[derik]

Спасибо, помогло!

[СергейД]

а теперь
msave,on
eqslv,pcg,1e-5
еще быстрее. вероятно, раза в полтора.

кстати, этот решатель вполне годится и для балочно-оболочечных моделей. но не любит жестких вставок и тп. при нем матрица упаковывается и занимает на порядок меньше места (и помещается в оперативке обычно)
решили недавно 42 миллиона уравнений (на кластере, правда)

[derik]

кхе...такое ощущение, что я плиту считаю, а не каркас спасибо Сергей! 42 миллиона а покажите в картинке, что это за конструкция такая? если не секрет

[aka raduntsev-nias]

веселенький каркасик!

[derik]

башня федераций или буржуй тауэр

[aka raduntsev-nias]

в сивиле есть глава о комбинациях

[IMG]http://s43.***********/i100/1005/86/d4d722c03075t.jpg[/IMG]

"ANSYS allows the results combinations fulfillment operating with different results sets
(Load Case) in a explicit way, that is, the user must define which results sets should
be operated and in which way must be combined in order to obtain the expected
result. The way followed for obtaining the result must be the same for the whole
model. Therefore, the user could do a combination of the following type:
“LoadCase1 = 1.20 * LoadStep1 + 0.95 * LoadStep3” "

[СергейД]

42 миллиона=
мы тестировали возможности кластера
и проверяли параллельный решатель ансис (вдруг там баги-глюки...)
и макросом последовательно сгущали сетку для здания (балки-оболочки) на грунте (солиды).
от 100 тыс уравнений до 7 миллионов уравнений.
все совпало. потом начали копировать модель и решать.
6 одинаковых зданий. тоже проверили.
с виду ничего особенного.
НО Теперь можно считать большие реальные задачи уверенно. (не с них же начинать)

в прочности решаю нелин= там задачи не слишком большие по объему, но долгие по врмени.

я сейчас занимаюсь, в основном расчетами на ветер CFX
там другие критерии. максимум 83 миллиона ячеек
в задаче 24 одинаковых скопированных модели застроек, но при разных направлениях ветра.
удобно графики и огибающие строить. все варианты- в одной задаче (поочередно их запускать неэффективно).
мои картинки про ветер = в теме про мост волгограда.

[Vlamos]

Сообщение от СергейД Теперь можно считать большие реальные задачи уверенно.

правильно ли я понимаю,что большие реальные задачи (немаленькие здания) Вы будете моделировать

Сообщение от СергейД балки-оболочки

а не

Сообщение от СергейД солиды

из-за размерности задачи?

[СергейД]

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

[Vlamos]

хотелось бы узнать - есть ли практика моделирования реальных немаленьких зданий солидами или нет?

[СергейД]

разумно верхние и типовые этажи смоделировать для экономии шеллами
а если и так считается= зачем усложнять.

Сити_палас в Сити (45 этажей+ грунт+ сваи) подойдет? (3 млн уравнений)
гостиница у стадиона Цска строящегося (40 эт+округа+грунт= 5 милл ур.). 10-15 этажные здания все набираются солидами.
проблема в привычке и в определении моментов для армирования.
сейчас объекты у нас в основном не высотные.

нелином, солидами считал 12-этажку (1 млн уравнений= 14 часов)
огнестойкость пробую.

от задачи зависит и от здания. и от тех, кто подбирает арматуру.
модель из ревита проще решить солидами, чем превращать ее в стерженьки и пластинки.

повторю, я сейчас занимаюсь немного иным.
а в чем вопрос? помочь посчитать солидами?

[derik]

Я знаю, я читать не умею но.у меня не получается разбить плиту из 186го правильными кирпичами, получаются только треугольники. пытался vsweep'ом: разбивал плиту на объемы, чтобы были ограничемы 4 плоскостями. соединял поверхности, затем линии lccat и accat, desize задавал. плиту без отверстия разбивал без проблем, а вот с отверстием - стопарнулся. направьте если где то не туда свернул.
З.Ы.:мой пост № 52, там есть эта плита, она вся в треугольниках

Все разобрался, оказывается нужно было, просто идти спать, а проснувшись и сев за ансис все получилось даже не пришлось лезть в хелп

Код:

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

finish
~cfclear,nostart
/prep7
/units,si
a=5
b=.4
c=.4
t=.2
h=3.2
q=1500
msave,on
eqslv,pcg,1e-5
et,1,solid186
mptemp,,,,,,,,
mptemp,1,0
mpdata,ex,1,,3e8
mpdata,prxy,1,,.2
mpdata,dens,1,,2200
blc5,b/2,b/2,b,b,-h
vgen,3,1,,,a,,,,0
vgen,5,all,,,,a,,,0
/view,1,1,1,1,
/ang,1
/rep,fast
blc5,b/2,(4*a+b)/2,b,4*a+b,c+t
vgen,3,16,,,a,,,,0
blc5,(2*a+b)/2,b/2,2*a+b,b,c+t
vgen,5,19,,,,a,,,,0
blc5,(a+b)/2,(a+b)/2,a-b,a-b,t
vgen,,24,,,,,b,,,1
vgen,2,24,,,a,,,0
vgen,4,24,25,,,a,,,0
blc5,(2*a+b)*3/4,(4*a+b)*7/8,a/2,a/2.5,t
vgen,,32,,,,,c,,,1
vsbv,31,32
vgen,6,all,,,,,h+c,,,0
vovlap,all
mshape,1
mshkey,0
esize,b,
vsweep,all
/replot
finish
/sol
asel,s,loc,z,-h
da,all,all
acel,0,0,9.8
allsel,all
solve
finish
/post1
set,first
plnsol,u,sum,0,1

Сергей, а можно попросить, при наличии у Вас свободного времени, покритиковать эту модель :-) на что мне следует обратить внимание. Заранее спасибо!

[Dron_629]

Сообщение от 6pa4o Сергей, а можно попросить, при наличии у Вас свободного времени, покритиковать эту модель :-) на что мне следует обратить внимание. Заранее спасибо!

Я хоть и не СергейД, но покритикую. При расчёте solid элементами по толщине нельзя ставить 1 элемент. Лучше от 3. А так вроде бы нормально.
Смотрю у тебя CivilFEM установлен (~cfclear - его команда). Тогда арматуру можешь посчитать через эквивалентные оболочки. Параграф 8.3 Help CivilFEM 12.0.

[derik]

Спасибо! Всмысле перекрытие и ригеля моделировать не одним а 3мя солидами?
З.Ы. а не могли бы Вы выложить хелп по цивилфему. в моем нет 8 и 9 главы. почему не знаю

[Dron_629]

По сечению должно быть не менее 2 элементов (а лучше от 3). Для ригеля в каждом направлении (по высоте и ширине сечения).
А какая версия CivilFEM?

[derik]

Попробуем :-) а у Вас хелп есть по цивилфему, Вам не сложно его выложить здесь? :-)
версия 11.0

[Dron_629]

Если у Вас версия CivilFEM старая, то такой функции там нет. А если 12.0 то она уже и в хелпе есть. Здесь выкладывать наверное некорректно, могу на e-mail.
Ищите новую версию .

[derik]

Я имел ввиду в даунлоад скинуть, пригодиться ;-) а так - [email protected], заранее спасибо! а скажите как правильно перекрытие задать: единой плитой( к примеру если здание 25х25 то и плиту перекрытия задавать цельным солидом(в 3 слоя) или плитами 5х5, к примеру?
З.Ы. новая есть как и ансис, но...он требователен к железу, чем 11. к сожалению :-(

[Dron_629]

Сообщение от 6pa4o перекрытия задавать цельным солидом(в 3 слоя) или плитами 5х5, к примеру?

Просто не могу понять что имели ввиду? Геометрическую модель или кэ после разбивки?
Запускал 12 версию на P4 одноядерник - 512 оперативки, и ничего работает. Возьмите хелп оттуда.

[derik]

геом. модель :-) просто при задании цельного перекрытия, 25х25 к примеру, у меня были проблемы с построением правильной сетки - vsweep'у не нравилось, что я пересек объемы, vovlap, и получилдось множество поверхностей. картинку прилагаю.
и правильная сетка не строилась, а вот когда между ригелями вставил плиту все получилось :-)

[Dron_629]

Есть множество путей постройки данной геом. модели. Главное чтобы после разбивки сетка была нормальная.

[derik]

а какой проще? и праильней, на Ваш взгляд?

[Cfytrr]

Расчетная программа другая, но вопрос из той же оперы когда нельзя моделить пластину одним элементом по высоте?
1. Тетра-10 1 по высоте
2 Гекса-8 1 по высоте
3. Гекса-8 5 по высоте

[Dron_629]

Моё мнение что 8 узловым никогда, а 20 узловым в принципе можно. Но это повторюсь только моё мнение без притензии на абсолютную истину.
бра4о, проще тот которым умеешь пользоваться. Я бы бил бы одну ячейку Вашего здания а затем копировал вместе с сеткой и объеденил бы совпадающие узлы (а дальше бы смотрел - если не нравится другой). Пробуйте и найдёте удобный Вам способ.

[Cfytrr]

Dron_629,Значит приведенные эскизы не верны? Можете скинуть скрин примера где применение одного элемента по высоте пластины выдает ложный (отклонение больше чем 5-10%) результат.

[Dron_629]

А зачем скидывать? На Ваших всё видно. Обратите внимание на распределение напряжений по толщине Ваших пластин.

[Cfytrr]

Сообщение от Dron_629 Обратите внимание на распределение напряжений по толщине

Полностью с вами согласен, распределение по толщине корректней в многослойной сетке, только вот в повседневной инженерной жизни как часто вам нужно распределение напряжений по толщине в пластине?

[Dron_629]

Вопрос был в правильности или в нужности? Как я помню в правильности. Ответ Вами получен, из Ваших же примеров.

[Cfytrr]

Сообщение от Dron_629 правильности или в нужности

Ну тут тонкая гань... можно договорится до того что результаты расчета никогда не бывают правильными....

[Dron_629]

Тогда зачем Вам объёмники? если не интерено распределение по толщине?

[Cfytrr]

Сообщение от Dron_629 Тогда зачем Вам объёмники? если не интерено распределение по толщине?

1. Иногда нужно создать модель в которой есть "пластинчатая" часть и есть "массивная" очень неудобно одну часть делать пластинами а другую солидами.
2. Иногда нужно решать контакт с торцом пластины (при соотв. сгущении сетки в области контакта)
3. Геометрию узлов где сходятся много листов гораздо проще и быстрей экспортировать и решать в солидах, чем предварительно тщательно и подробно прорисовывать "срединные поверхности" для решения в пластинах.
во всех решениях интересуют напряжения в крайних фибрах

[СергейД]

dron_629
вы тестировали сами 1 элемент по толщине?
решите перекрытие 186 типом (опции Enh strain). 1,2,3 эл для перемещений и частот идентичны практически. по напряжениям разница незначительна. а начинающим- хорошее упражнение.
+есть 190 тип который всегда 1 по толшине...
и т.д.
у меня на эту тему много ПРАКТИЧЕСКОГО опыта.

насчет разбивки- советую отрезать (vsbw) топологически простые объемы, пригодные для VSWEEP, а остальное бить тетраэдрами.
в 12 версии появились специальные 4узловые тетры тип 285= рекомендую.

к сожалению, на анализ модели нет пока врмени...и другой софт перед глазами.

[derik]

Сколько советов. а на проверку нет времени, только вечером. жалко. а руки то чешуться сесть и попробывать Спасибо Всем ответившим!

[Dron_629]

Cергей Иванович, не рассмотрел что там в коде 20 узловой элемент был, так что вопрос с конкретной моделью снят по части элементов на толщину, но я всё равно уверен что 8 узловой элемент один на толщину брать не следует.

[СергейД]

специальный 190 есть, а 45 и тем более 185 действительно не стоит.
хотя для верхних этажей почему бы и нет...
допустим получается тяжелая модель из 186.
в верхних 2/3 высотки можно смело менять 186 на 45 через Emodif.
и EMID,remove

[derik]

Эт Вы про мою модель или нет?

[Dron_629]

190 элемент не пробовал, буду иметь ввиду.
Только для придания весовых и жёсткостных характеристик один слой объёмников будет ни чем не хуже оболочек. Но всё же для интересных мест использовать стандартные 8 узловые элементы я бы не стал.
В тех конструкциях которые я обычно считаю толщина значительная по крайней мере больше 0,5 м. и мне спокойнее когда 3 элемента.

[СергейД]

не спорю. речь о перекрытиях. они по 0.5м бывают крайне редко.
в линейных задачах жб я бы о 3% не думал
и разбивал бы экономно, просто чтобы получить общую картину
. Ошибка ведь в другом.
(в трещинах, в монтаже)
нужно примерно определить моменты. подобрать арматуру
вырезать типовой фрагмент и решить нелинейно

в нелинейных задачах 65 обычно 3-4 элемента по толщине (по крайним размазывается арматура)

но если говорить о точности, то нужно начинать с грунта.
скажем, какой смысл ставить несколько эл по толщине, но применять коэфф постели равномерный по основанию и брать его из кросса.
нужно тогда и грунт представлять друкер-прагером и тп.
все должно быть "соразмерно"

[Dron_629]

Полностью согласен.

[derik]

Вот пример из Лиры №9 решил в ансисе. :-) где наш МишаИнженер, пусть сравнит "блокнотный" расчет с лирой :-D

Код:

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

finish
~cfclear,nostart
/prep7
/units,si
a=5.2
b=7.6
c=.4
d=.6
f=.8
h=3
t=.5
t2=.15
p=1.96e4
p2=.98e6
q=.98e4
q2=.49e4
msave,on
eqslv,pcg,1e-5
et,1,186
mptemp,,,,,,,,
mptemp,1,0
mpdata,ex,1,,3e9
mpdata,prxy,1,,.2
mpdata,dens,1,,2750
blc5,b/2,a/2,b,a,t
blc5,f,d,f,c,h
vgen,2,2,,,,a-2*d,,,,0
vgen,2,2,3,,b-2*f,,,,0
vgen,,2,,,,,t,,,1
vgen,,3,,,,,t,,,1
vgen,,4,,,,,t,,,1
vgen,,5,,,,,t,,,1
blc5,b/2,a/2,b-c,a-c,t2
vgen,,6,,,,,h+t,,,1
mshape,1
mshkey,0
esize,c/2
vsweep,all
nummrg,all
finish
/sol
asel,s,,,1
da,all,all
acel,0,0,9.8
allsel,all
lswrite,1,
sfa,2,1,pres,q
sfa,32,2,pres,q2
nsel,s,loc,x,2
nsel,r,loc,y,1.6
nsel,r,loc,z,t
f,all,fz,-p2
nsel,s,loc,x,2
nsel,r,loc,y,3.4
nsel,r,loc,z,t
f,all,fz,-p2
nsel,s,loc,x,5.6
nsel,r,loc,y,1.6
nsel,r,loc,z,t
f,all,fz,-p2
nsel,s,loc,x,5.6
nsel,r,loc,y,3.4
nsel,r,loc,z,t
f,all,fz,-p2
allsel,all
lswrite,2,
f,4973,fx,p
allsel,all
lswrite,3,
lssolve,1,3,1,
finish
/post1
set,first

[aka raduntsev-nias]

6pa4o, грунтового основания не видно.

[derik]

вот теперь займусь грунтом :-) не все сразу ;-)

[derik]

Код:

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

finish
~cfclear,nostart
/prep7
a=3
b=9
c=.3
h=1.2
h1=2.4
msave,on
eqslv,pcg,1e-5
et,1,45
keyotp,1,2,1
et,2,186
MPTEMP,,,,,,,,  
MPTEMP,1,0  
MPDATA,DENS,1,,2000 
MPTEMP,,,,,,,,  
MPTEMP,1,0  
MPDATA,EX,1,,2.5e8 
MPDATA,PRXY,1,,.3   
TB,DP,1,,,  
TBMODIF,1,1,65000   
TBMODIF,1,2,20  
TBMODIF,1,3,0   
MPTEMP,,,,,,,,  
MPTEMP,1,0  
MPDATA,DENS,2,,1900 
MPTEMP,,,,,,,,  
MPTEMP,1,0  
MPDATA,EX,2,,3e8  
MPDATA,PRXY,2,,.3   
TB,DP,2,,,  
TBMODIF,1,1,1000
TBMODIF,1,2,37  
TBMODIF,1,3,7   
MPTEMP,,,,,,,,  
MPTEMP,1,0  
MPDATA,DENS,3,,2500 
MPTEMP,,,,,,,,  
MPTEMP,1,0  
MPDATA,EX,3,,3e8
MPDATA,PRXY,3,,.2
blc5,b/2,a/2,b,a,c
vatt,3,,2,
blc5,b/2,-a/2,b,a,-h
vgen,2,2,,,b,,,,,0
vgen,2,2,,,-b,,,,,0
vsel,s,volu,,2,4,
vgen,3,all,,,,a,,,0
vatt,1,,1,
blc5,b/2,-a/2,b,a,-h1
vsel,s,,,11
vgen,2,11,,,b,,,,,0
vgen,2,11,,,-b,,,,,0
vsel,s,volu,,11,13,
vgen,3,all,,,,a,,,0
vgen,,all,,,,,-h,,,1
vatt,2,,1,
allsel,all

Добрый день! Вот пытаюсь смоделировать грунт через Д-П. чтот у меня плита висит в воздухе после расчета. что здесь не так?

[СергейД]

сначала нужно научиться решать контактную задачу
грунт-фундамент. это важнее. и сложнее.

а уж потом пробовать DP.
проверять файл не могу.
но замечу, что 45 и 186 не стыкуются
(8 и 20 узлов) без ухищрений.
для 186 dp нельзя.

начинать нелин с Dp не стоит. нужно научиться контакту, а потом пластике металлов.
потом Dp
потом бетон

у меня нелин сопровождется 10-15 строками параметров. у тебя ничего не вижу.
как решать собираешься?

[shtoss]

Ну наконец-то подняли тему расчета жбк в ANSYS я уж думал это вообще не произойдет. Можно пожалуйста хоть какаю-нибудь книжку для студентов чтобы можно было как лиру освоить, именно к расчету жбк+грунт, тут вроде мелькали страницы из какого-то руководства по CIVILFEM можно его полностью почитать.

[aka raduntsev-nias]

shtoss, в самом ansys есть справка(и в сивиле). учебную версию возможно найти тоже можно. литературу можно найти в http://dwg.ru/dnl/.

[derik]

спасибо Сергей Иванович! буду дальше разбираться :-)
то shtoss, такой книжки нет :-) вот только если Сергей Иванович не собереться что нибудь написать (нужно его дружно попросить) а те страницы - скачай цивилфем и будет тебе счастье :-)

[СергейД]

попробую к осени что-нибудь составить. Буду(надеюсь) курсы проводить.
а bra40 молодец. прогресс очень наглядный.
теперь советую освоить пластику металла -дело несложное. какой-нибудь узел из снипа посчитать
в больших перемещениях.
потом порешать с учетом контакта классическое- типа задачи Герца.
потом порешать контакт фундамент-грунт. освоить хорошо опции сходимости.
потом друкер-прагера понять.
сначала обычного. потом EDP.
а тогда с ЖБ можно уверенно работать.
все просто.
делайте. разбирайтесь.
в макросе команды подсказать правильные время всегда найду.

[derik]

Спасибо ! Вот все таки подумываю на Ваши курсы все таки подъехать если ни чего не случиться :-) теперь полез в дебри контактных задач Кстати, Сергей Иванович как там поживает книга http://forum.dwg.ru/showthread.php?p=352600&#post352600?

[derik]

Добрый вечер! Таки вот:

Код:

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

finish
~cfclear
/prep7
msave,on
eqslv,pcg,1e-5
a=.296
b=.141
t=.0125
t1=.018
l=10
p=200
et,1,186
mptemp,,,,,,,
mptemp,1,0
mpdata,ex,1,,2.06e8
mpdata,prxy,1,,.3
mpdata,dens,1,,7850
blc5,b/2,t1/2,b,t1,
blc5,t/2,a/2,t,a
agen,2,1,,,,a-t1,,,0
agen,,2,,,b,,,,,1
agen,2,1,,,b,,,,0
agen,2,3,,,b,,,,0
aadd,all
vext,all,,,0,0,l,,,,
/view,,1,1,1
/ang,1
/rep,fasta
mshape,1
mshkey,0
lesize,all,t*10,
vsweep,all
nummrg,all
finish
/sol
asel,s,,,1
asel,a,,,6
da,all,all
allsel,all
asel,s,area,,10,12,1
sfa,all,,pres,p,
acel,0,0,9.8
allsel,all
solve
finish
/post1
set,first
/efacet,1
plnsol,u,sum,0,1

Код:

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

finish
~cfclear
/prep7
msave,on
eqslv,pcg,1e-5
/view,,1,1,1
/ang,1
/rep,fasta
p=200
et,1,188
mptemp,,,,,,,
mptemp,1,0
mpdata,ex,1,,2.06e8
mpdata,prxy,1,,.3
mpdata,dens,1,,7850
k,,0,0,0
k,,5,0,0
lstr,1,2
sectype,1,beam,I,,0
secoffset,cent
secdata,.3,.3,.296,.018,.018,.0125,0,0,0,0
lesize,all,.2
lmesh,all
/eshape,1
finish
/sol
acel,0,0,9.8
ksel,s,,,1
ksel,a,,,2
dk,all,all
sfbeam,all,1,pres,p,p,,,,,
solve
finish
/post1
set,first
/efacet,1
plnsol,u,sum,0,1

это балкой
Ваши комментарии

[viking1963]

Доброго времени суток, уважаемые участники темы.
Хотелось бы получить небольшую консультацию.
Делаем проект в сейсмоопасном районе.
Я работаю над Административным зданием (монолит, 5 этажей). Расчетную часть я уже закончил (SCAD v.11.1), и если у меня будет время то продублирую расчеты в Лире 9.6 или в Старке, приготовил документацию в виде расчетно-пояснительной записки к отправке в Экспертизу. Рабочие чертежи конструкций по моему зданию готовы.
Работа у нас поставлена так, что нет главного конструктора (просто что бы исключить ненужные словопрения и давления). Всего два ведущих, и у каждого свой объект и объем работы. Так вот второй ведущий получил задачу от руководства спроектировать забор-ограждение. Второй ведущий конструктор - это дедок лет 59, всю жизнь крутившийся в машиностроении. Не мудрствуя лукаво, он загнал геометрию забора в ANSYS и стал получать там напряжения и перемещения забора от сейсмовоздействия...........
Чертежи забора он взял из типовика советских времен........
Там армирование забора - арматура диаметром 8 мм. Дедок прикинул - ускорения забора от сейсмовоздействия возрастают ровно в 8 раз, и решил армировать забор арматурой 36-го-диаметра....
Я сижу работаю над Административным зданием и на подходе у меня следующий блок - Спорткомплекс, так что мне его забор по-барабану, что зайцу балалайка..........
То есть у меня своих дел хватает, и мне не до него....... Дедку 1 год до пенсии......
Так вот, сегодня поймал себя на мысли - че он делает то ?
И есть ли смысл выполнять расчет армирования забора в программе ANSYS ? (а самое забавное - как эт так у него это армирование то получается?) .

P.S. : до этого дедок пытался получить напряжения в лестничном монолитном марше (длина марша 3 метра, ширина 1 метр) в программе ANSYS, и соответственно заармировать лестничный марш - к сожалению здесь только фото забора - и то немножко размытое.........
лесенку я у него сфоткать не успел..........

[СергейД]

умеючи-то в чем угодно можно.
получил моменты в ансис и армируешь по снипу.
в чем сомнения-то?

но строго говоря.
в программе ансис (в основном модуле) нет в принципе подбора врмирования.
можно задать какое-нибудь и проверить нелинейно.
потом при желании исправлять и пересчитывать.
модели доложны быть из объемных элементов.
в идеале 3-4 по толщине.
нелинейный расчет подразумевает ввод акселерограммы...
не уверен, подойдет ли вам это все.
вы хотите сами научиться или кого-либо попросить?
я лично еще месяц занят. да и задача трудоемкая. надо всерьез.

возьмите что-нибудь полегче- типа лиры...

[viking1963]

СергейД,
Доброго времени суток,

я лично еще месяц занят. да и задача трудоемкая

?
Не понял - Вам трудоемко посчитать армирование забора? - да оно мне на ф*г не надо, пусть дедок сам свой забор в программе ANSYS считает, а по поводу Лиры и Скада - зря Вы это.
Мои объекты, которые реально были построены были расчитаны именно с применением Лиры и Скада, просто какой смысл использовать машиностроительную программу ANSYS в строительстве ? На скринах - мои объекты

в программе ансис (в основном модуле) нет в принципе подбора армирования. можно задать какое-нибудь и проверить нелинейно. потом при желании исправлять и пересчитывать.

[СергейД]

2viking
я могу вывесить картинки посчитанного в ансис.
и моими коллегами, тем более. ну и что?
посчитали моменты, подобрали арматуру экселем..

кто чем умеет.. и кому что нравится.
просто считает раз в десять быстрее... и нелинейные возможности несопоставимы.
он просто не только для заборов...
не ансис виноват, что товарищ выбрал такую "армутуру".

уточните, в чем идея вашего поста в теме, где собрались им интересующиеся? та же самая, что и "МИшиИнженера", который "тащился", нас читая?

[igr]

Зачем достаточно простые конструкции( забор, 5 этажное здание) считать на такой крутой программе как Ansys, когда можноэто выполнить в болеепростых прогах?

[viking1963]

СергейД,
Доброго времени суток,

я могу вывесить картинки посчитанного в ансис.

уточните, в чем идея вашего поста в теме

Просто я пытаюсь разобраться - начать изучение ANSYS и PLAXIS, или нет.
Наверное, как только появится хоть немного свободного времени, то помаленьку начну.
Если можно, приведите скрины зданий, посчитанных в ANSYS, и допустим картинку усилий, в колоннах, балках (стержневых элементах), плитах, диафрагмах (в пластинчатых элементах) : (N, Q, Mкр, Мx, Му) - то есть исходников по которым можно выполнить подбор арматуры.

С уважением, viking1963
P.S.: не обижайтесь...........

[igr]

Наверное времени невпроворот и зарплату платят просто так

[aka raduntsev-nias]

Сообщение от igr Зачем достаточно простые конструкции( забор, 5 этажное здание) считать на такой крутой программе как Ansys, когда можноэто выполнить в болеепростых прогах?

когда я ездил на тюнингованной наглухо тонированной ВАЗ 12шке я все равно хотел машину лучше. Сейчас у меня машина более лучше. Но все равно.

Offtop: нет ничего лучше BMW x7!

[Vlamos]

2ceргейД - извините,не понял

получил моменты в ансис и армируешь по снипу. ... модели доложны быть из объемных элементов.

реальные здания по акселерограмме все-же как вы считаете в Ansys ?

[СергейД]

2 viking igr
если вам хватает того, что уже знаете, значит ансис вам ни к чему.

считаем в основном объемниками.снимаем моменты= подбирать арматуры
достаточно экселя (я этим не занимаюсь)
лично я делаю поверочные расчеты нелинейные подобранного.
прогрессирующее, огнестойкость

2 vlamos
.
резервуара сейсмику считал объемниками.
из оболочек здание считал на акс-мму и не одно.
прогрессирующее для здания считал объемниками. (фрагмент)
огнестойкость фрагмента из нескольких этажей считал.
считал по акселерограмме покрытие ДС Мегаспорт (стержни-оболочки)
сейсмику зданий из объемников не считал. закажут- посчитаю.
проблем не вижу.
можно и LS-dyna привлечь (есть легальная и владею почти как ансисом)

2 viking
картинки вывешивать повода не вижу.
очень занят аэродинамикой.
приезжайте еще= побеседуем (я ж Вам предлагал уже)
plaxis чем больше изучать, тем больше к нему вопросов...
он для геотехников, которые и без него знают, сколько должно быть..
трата времени.
уж поверьте.
грунты и в ансис есть. можно все в связке считать.

[derik]

Вечер добрый! Я конечно дико извиняюсь, но как мне повернуть арию/волюмэ (уже созданную) вокруг какой нибудь из осей? Спасибо заранее!

[СергейД]

local,11,1,0,0,0,0,0,90
vgen,2,all,,,0,45,0,,,1
csys,0
будет поворот вокруг х на 45

csys,1
vgen,2,all,,,0,120,0,,,1
csys,0
поворот вокруг z на 120

это если они не связаны ни с чем.

иначе
нужно создать копию

vgen,2,7,,,0,120,0
vclear,7
vdel,7,,,1

для обхема 7

agen для поверхностей

[derik]

Спасибо!

[derik]

Вечер добрый! Вот на очередную порцию критики выставляю:

Код:

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

fini
~cfclear
/undo,on
/prep7
msave,on
eqslv,pcg,1e-5
a=1
b=.5
r=.3
t=.02
q=500
et,1,186
mptemp,,,,,,,,
mptemp,1,0
mpdata,ex,1,,2.06e8
mpdata,prxy,1,,.3
mpdata,dens,1,,7850
blc5,,,a,a,t
cyl4,,,r,,r-t,,b
blc5,0,0,t,a,t
vext,14,,,0,0,b-t,,.7,,
csys,1
vgen,4,3,4,,,90,,,0
csys,0
/view,1,1,1,1
/ang,1
/rep,fast
vptn,all
vdele,23,24,,1
vdele,47,48,,1
vdele,36,,,1
allsel,all
asel,s,area,,131,134,3
asel,a,area,,150,153,3
asel,a,area,,158,161,3
asel,a,area,,226,227,1
vext,all,,,0,0,10*b,,,0
allsel,all
vplot
mshape,1
mshkey,0
esize,all,b/10
vsweep,all
nummrg,all
finish
/sol
asel,s,area,,57,82,5
asel,a,area,,90,95,5
asel,a,area,,125,128,3
asel,a,area,,135,139,4
asel,a,area,,107,111,4
asel,a,area,,206,210,4
asel,a,area,,222,224,2
asel,a,area,,170,176,6
asel,a,area,,191,194,3
asel,a,area,,198,203,5
asel,a,,,146
da,all,all
asel,s,area,,1,8,7
asel,a,area,,17,35,6
asel,a,area,,39,44,5
sfa,all,,pres,q
allsel,all
acel,0,0,9.8
solve
finish
/post1
set,first
plnsol,u,sum,0,1

[derik]

Добрый день! вот, никак н могу разобраться...вроде все как по мануалу...но что то закралась ошибка или жара так действует...
почему не генерятся контактная и ответная поверхности? заранее спасибо!

Код:

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

fini
/clear,nostart
/undo,on
/prep7
msave,on
eqslv,pcg,1e-5
/view,1,1,1,1
/ang,1
/rep,fast
a=4
b=8
t=.5
t1=1
mptemp,,,,,,,
mptemp,1,0
mpdata,ex,1,,3e8
mpdata,prxy,1,,.2
mpdata,dens,1,,2500
et,1,185
et,2,170
et,3,174
blc5,b/2,b/2,a,a,-t1
blc5,b/2,b/2,a,a,t
esize,t/2
vsweep,all
et,2,170
et,3,174
r,3
real,3
r,3,,,.1,0.1,0,
rmore,,,1.0E20,0.0,1.0,
rmore,0.0,0,1.0,,1.0,0.5
rmore,0,1.0,1.0,0.0,,1.0
keyopt,3,4,0
keyopt,3,5,0
keyopt,3,7,0
keyopt,3,8,0
keyopt,3,9,0
keyopt,3,10,2
keyopt,3,11,0
keyopt,3,12,0
keyopt,3,2,0
keyopt,2,5,0
asel,s,,,7
nsla,s,0
esln,s,0
esll,u
cm,contact,area
esurf
type,3
asel,s,,,2
nsla,s,0
esln,s,0
esll,u
cm,target,area
type,2
esurf

все открывает :-)

[aka raduntsev-nias]

Offtop: даже не открывается полностью.
пишет
The ESURF command can produce only line or area elements.

наверное из за этого

[СергейД]

друзья! ансис под рукой нет.
но type,3
надо поставить перед esurf.
иначе она контакты не сгенерирует (type,1 по умолч)

контакты дело серьезное.
советую несколько раз перечитать хелп...
вообще применяю
keyopt,,2,1- штрафной алгоритм-
в моих задачах он лучше сходится

[derik]

спасибо! Сергей Иванович, а что это за команда в мануале по контакту -

Использовать команду GCGEN (маршрут меню Main Menu>Preprocessor> Create>Elements>At ContactSrf) для идентификации исходной и ответной поверхностей с помощью имен и для генерации контактных элементов между этими именованными поверхностями

это цитата из 9 главы хелпа. чтот такой команды я не найду в описании. она вообще существует?

[aka raduntsev-nias]

Мне конечно Очень Далеко до уважаемого Сергея Иановича. Но я хочу поделиться тем, что нашел по этому вопросу.

ps все таки не зря убил семь месяцев на сбор информации!!)

Точно, Брачо. Пишет
GCGEN is not a recognized RUNSTAT commandj abbreviation or macro This command will be ignored.

Спасибо, Сергей Иванович.

[derik]

вот вот, я это тоже ьнашел. а теперь попробуй ввести эту команду в ансисе или найти ее в хелпе или перечне команд. будешь приятно удивлен

[СергейД]

это 2000 год. все отменено...
читайте хотя бы К.Басова.

контакт узел-узел делается типом 177
eintf и тп.

[derik]

Код:

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

fini
/clear,nostart
/undo,on
/prep7
msave,on
eqslv,pcg,1e-5
/view,1,1,1,1
/ang,1
/rep,fast
a=4
b=8
t=.5
t1=1
mptemp,,,,,,,
mptemp,1,0
mpdata,ex,1,,3e8
mpdata,prxy,1,,.2
mpdata,dens,1,,2500
et,1,185
et,2,170
et,3,174
blc5,b/2,b/2,a,a,-t1
blc5,b/2,b/2,a,a,t
esize,t/2
vsweep,all
et,2,170
et,3,174
r,3
real,3
r,3,,,.1,0.1,0,
rmore,,,1.0E20,0.0,1.0,
rmore,0.0,0,1.0,,1.0,0.5
rmore,0,1.0,1.0,0.0,,1.0
keyopt,3,4,0
keyopt,3,5,0
keyopt,3,7,0
keyopt,3,8,0
keyopt,3,9,0
keyopt,3,10,2
keyopt,3,11,0
keyopt,3,12,0
keyopt,3,2,0
keyopt,2,5,0
asel,s,,,7
nsla,s,1
cm,contact,area
type,3
esurf
asel,s,,,2
nsla,s,1
cm,target,area
type,2
esurf
fini
/sol
acel,0,0,9.8
da,1,all,
allsel,all
solve
finish
/post1
plnsol,u,sum,0,1

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

[СергейД]

нужно отправить вывод в файл
/outp,,outp
и внимательно изучить потом протокол.
стараясь разобраться в каждой строчке.
там сплошные советы будут...
+
все же
keyopt,3,4,1 не помешает.
185 ЛУЧШЕ ЗАМЕНИТЬ НА 186
ИЛИ 174 НА 173- они должны соответствовать по числу узлов.

[derik]

день добрый!
вот жду замечаний:

Код:

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

fini
/clear,nostart
/undo,on
/outp,123,txt,C:\ansys\output,
/cwd,'C:\ansys\out' 
/prep7
msave,on
eqslv,pcg,1e-5
a=12
b=4
h=3
t=.3
q=200
blc5,a/2,-h/2,a,-h,
blc5,a/2,t/2,b,t,
et,1,82
keyopt,1,3,2
mptemp,,,,,,,
mptemp,1,0
mpdata,ex,1,,3e7
mpdata,prxy,1,,.2
!mpdata,dens,1,,2500
mptemp,1,0
mpdata,ex,2,,2.1e6
mpdata,prxy,2,,.2
!mpdata,dens,2,,1350
asel,s,,,1
aatt,2,,1,
asel,s,,,2
aatt,1,,1,
allsel,all
mshkey,0
mshape,0,2d
esize,t/5
amesh,all
et,2,targe169
et,3,conta172
real,10
type,2
lsel,s,,,3
nsll,s,1
esurf
type,3
lsel,s,,,5
nsll,s,1
esurf
allsel
keyopt,3,5,3
fini
/sol
lsel,s,,,1
lsel,a,line,,2,4,2
dl,all,,all,
allsel,all
!acel,0,9.8,0
sfl,7,pres,q,
antype,static
autots,on
lnsrch,auto
solve
fini
/post1
set,first
plnsol,u,sum,0,1

[aka raduntsev-nias]

Offtop: уже и до основания добрался

[aka raduntsev-nias]

Тоже сижу решаю какой путь выбрать. Или построение модели средствами ansys. Но больше склоняюсь к построению с помощью поддерживаемых им Cad системами с последующим экспортом из них в ansys workbench, расчете в нем и последующей передачей в классический для более быстрого с ним ознакомления, ведь это все реализовано. И плюсы налицо, хотите параметров, ansys импортирует из Cad system параметры. Даже сварку можно передать.

[derik]

лучше классику сначала освоить, а потом уже ВБ.
в плане создания геометрии ансис самодостаточен, может построение и заморочено, но это только по началу, поверь мне :-), а потом все легче и легче. в конструкции, для анализа, нужно изменять геометрию, как писал Сергей Иванович, будешь каждый раз перерисовывать? проще изменить параметр и пустить на расчет заново - 5 сек!
а сварка...можно и в классике. там все можно. мне, лично, ансис нравится все больше и больше!!
все сказанное сугубо ИМХО.

[aka raduntsev-nias]

Сообщение от 6pa4o в конструкции, для анализа, нужно изменять геометрию, как писал Сергей Иванович, будешь каждый раз перерисовывать? проще изменить параметр и пустить на расчет заново - 5 сек!

хорошо, я понял тебя. согласен. но изменять геометрию и в workbench несложно с кнопки update(быстро изменив ее в модели) и нагрузки на месте вроде остаются. конечно, если у тебя плита. а если какой нибудь узел структурой плиты МАРХИ или подобное с подвыподвертом. тут наверное сложнее описать параметрически. и большой плюс, что имеется возможность все это в классику перекинуть.

[derik]

дело не в кнопках, дело в прогармировании. ты сам выбираешь как и что будет генериться. про ВБ молчу - открыл и закрыл, мне не опнравилось
а параметрическое описание..я уже модель для акада, после ночи, сел писать в блокноте, а потом долго втыкал - почему он не хочет ее строить помучаешься недельку, а потом все войдет в колею

З.Ы. попробуй в любом из вывешенных мною макросов изменить параметр, геометрии или материала, - получишь новую модель и новые результаты. попробуй :-)

[aka raduntsev-nias]

это конечно. пробовал, спасибо за коды. Но конструкции на скринах, думается мне, быстрей получить в кадах. Я их уже знаю тем более. Fortran забыл почти( А конструкции срочно считать надо.

[IMG]http://s02.***********/i175/1006/6e/8e1b86f23af1t.jpg[/IMG]

[IMG]http://i065.***********/1006/1a/bf4a0f122cc5t.jpg[/IMG]

[derik]

а это не текла случаем??? да нет построить ее в ансисе можно, считай чертишь в каде каком нить. просто через параметры удобней ИМХО.
на вкус и цвет фломастеры разные.

[aka raduntsev-nias]

да, она. это точно!

[derik]

Доброе утро! подскажите пожалуйста, почему моменты у меня выводятся в такой форме? что я не так указал?
вот код:

Код:

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

fini
/clear,nostsrt
/undo,on
!/outp,123,txt,C:\ansys\output,
!/cwd,'C:\ansys\out' 
/prep7
msave,on
eqslv,pcg,1e-5
a=.4
b=.5
c=.8
k=.1
h=.3
t=.03
e=3
q=500
et,1,43
keyopt,1,5,1
r,1,t,0,0,0,0,
mptemt,,,,,,,
mptemp,1,0
mpdata,ex,1,,2.06e8
mpdata,prxy,1,,.2
mpdata,dens,1,,7850
k,,0,0,0
k,,0,c,0
k,,a,c,0
k,,b+k,c-h,0
k,,b,0,0
a,1,2,3,4,5,
blc5,a/4,c/2,a/2,e,
blc5,-a/4,c/2,-a/2,e,
blc5,0,c/2,h,e,
csys,5
agen,,1,,,,90,,,,1
agen,,4,,,,90,,,,1
csys,0
agen,,4,,,a/2,,,,,1
agen,2,4,,,-a,,,,,0
aovlap,all
aesize,all,a/16
amesh,all
nummrg,all
/eshape,1
allsel,all
fini
/sol
acel,0,9.8,0
lsel,s,line,,6,32,26
lsel,a,line,,22,28,2
dl,all,,all
lsel,s,line,,8,33,25
lsel,a,line,,23,29,2
sfl,all,pres,q,
lsel,s,,,3
sfl,all,pres,q/2,
allsel,all
solve
fini
/post1
set,first
etable,,m,x
avprin,0,,
etable,,m,y
avprin,0,,
etable,,m,z
etable,,f,x
avprin,0,,
etable,,f,y
avprin,0,,
etable,,f,z
pletab,mx,avg

разобрался :-D

[aka raduntsev-nias]

У меня тоже неясности.

[IMG]http://s06.***********/i179/1007/c4/3404b225b68bt.jpg[/IMG]

[derik]

воркбенч?
а ты контакт болтов с уголком задал?

[aka raduntsev-nias]

Сообщение от 6pa4o воркбенч? а ты контакт болтов с уголком задал?

да. задал . и уголка с балкой трение. болты без трения

[derik]

тогда хз...считай в класике
я вот с моментом разобрался как задать и как посмотреть. все оказалось просто - апдл рулит
а покажи скрин сбоку - вид на болты во фронте

[aka raduntsev-nias]

Просмотрел файл урнала. много преупреждений. буду разбираться

[IMG]http://s44.***********/i103/1007/da/0d29bc1ad752t.jpg[/IMG]

[IMG]http://s47.***********/i115/1007/76/1e5a9fceb9e3t.jpg[/IMG]

[IMG]http://s51.***********/i131/1007/9a/23b0d42ce76dt.jpg[/IMG]

[IMG]http://s54.***********/i143/1007/cf/68de8bac2e25t.jpg[/IMG]

узелок маленький, считается долго. и решатель менял и все. что такое

[derik]

еще больше не понятно что и как исправлять. попробуй болты с трением поставить, те с проскальзываньем
и уголок не передает усилие на двутавр..
а выложи логи расчета

[aka raduntsev-nias]

Сообщение от 6pa4o еще больше не понятно что и как исправлять. попробуй болты с трением поставить, те с проскальзываньем и уголок не передает усилие на двутавр.. а выложи логи расчета

ошибки расчета. если так интересны скину в PM

болты с bolt load и lock/

ладно. поразбираюсь еще. потом отпишусь, когда решу.

[aka raduntsev-nias]

обидно. при обновлении геометрии нагрузки улетели. может что не нажал)

[derik]

:-) считать учись в классике

[aka raduntsev-nias]

так и знал, что так скажешь)))) да, придется поднапрячься.

[derik]

скинь файл и журнал. помочуюсь и я :-)

[aka raduntsev-nias]

попробуй параметрами в apdl. я например собрал узелок и задал нагрузки в течении часа (Спасибо, Willy!) (не спеша и аккуратно)+ минут 5 на редактирование после update. + ко всему этому имею еще и файл, для классика. кроме этого можно получить импорт параметров из Cad модели при желании. интересно, сколько времени это займет у тебя? и представим, что еще и рама с узлами опирания на фундамент.
И, например, как поступать в том случае, если у нас изменится сечение колонны например на трубу. Мне понадобится не больше 5 минут на исправление.
Конечно, при больших объектах с несколькими этажами удобно иметь параметры для изменения. Но если сравнивать параметрическое моделирование в Tekla и задание параметров с помощью apdl, и возможность постоянного его применения, то скорее всего, на практике не так часто будет возможно его применять (пока не накопится много макросов и будет спроектирован не один объект с помощью классика). Конечно, исследуя определенную группу конструкций, без них никуда. Но при гонке со сроками с индивидуальной архитектурой, возможно не всегда найдется время под нее подстраивать свои конструкции.
Я конечно же за apdl. И тоже изучаю его.

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

[derik]

я все размеры через параметры, так что мне тоже только исправить значение

[aka raduntsev-nias]

Еще есть abaqus, краем глаза удалось увидеть, что в нем делают. Прямо таки управление природой. Тоже классный продукт.

[derik]

угу...наслышан...но сначала АНСИС! там тоже люди чудеса творят...а мы пока так, страдаем ерундой

[aka raduntsev-nias]

Сообщение от 6pa4o угу...наслышан...но сначала АНСИС! там тоже люди чудеса творят...а мы пока так, страдаем ерундой

Offtop: это точно!

[aka raduntsev-nias]

можно ли в ansys увидеть разрушение материала визуально? например в abaqus видел скрин, как конструкция рассыпается. тоже бы хотелось увидеть груды искореженного метала, последствия большого тоннажа

[СергейД]

расчет по неявной схеме подразумевает целостность конструкции- отсутствие нулевых частот.
для рассыпания нужна ls-dyna или autodyn.
работающие по явной схеме интегрирования
или abaqus/explicit (а не обычный /standard)=
та картинка получена в последнем.

просто нужно иметь в виду, что правильно в этом подходе считается
обычно только доли секунды. далее идет накопление погрешностей.
и просто красивая мультипликация.
чуть подробнее об этом в моей статье в САПР и ГРафика 2004 год№3

[aka raduntsev-nias]

Спасибо, Сергей Иванович. Будем разбираться.

[derik]

Добрый вечер! Сергей Иванович, не могли бы Вы сказать свое мнение? :-)

Код:

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

fini
/clear,nostart
/undo,on
/outp,123,txt,C:\ansys\output,
/cwd,'C:\ansys\out' 
/prep7
msave,on
eqslv,pcg,1e-5
/view,1,1,1,1
/ang,1
/rep,fast
r=2
p=500
cyl4,,-r,r,90,
cyl4,,r,r,-90,,
lccat,1,2
lccat,4,6
sfl,5,press,p
et,1,82
mptemp,,,,,,,
mptemp,1,0
mpdata,ex,1,,1.85e8
mpdata,prxy,1,,.3
mptemp,1,0
mpdata,ex,2,,2e8
mpdata,prxy,2,,.2
asel,s,,,1
aatt,2,,1,
asel,s,,,2
aatt,1,,1,
allsel,all
mshkey,0
mshape,0,2d
esize,r/20
amesh,all
eref,736,,,1,2,,
eref,338,,,1,2,,
et,2,targe169
et,3,conta172
real,10
r,10,0,0,1.23,0,0,1
keyopt,2,5,0
keyopt,3,2,0
keyopt,3,4,1
keyopt,3,5,0
keyopt,3,7,0
keyopt,3,8,0
keyopt,3,9,0
keyopt,3,10,2
keyopt,3,11,0
keyopt,3,12,6
type,2
lsel,s,,,1
nsll,s,1
esurf
type,3
lsel,s,,,4
nsll,s,1
esurf
allsel,all
fini
/sol
/gst,1
antype,0
nlgeom,1
nsubst,20,0,0
outres,erase
outres,all,all
time,1
dl,3,,all,
lsel,s,line,,2,6,4
dl,all,,symm
allsel,all
!cncheck
solve
fini
/post1
set,first
etable,,f,x
avprin,0,,
etable,,f,y
avprin,0,,
plnsol,s,1,0,1
!/expand,9,axis,,,10

[СергейД]

извините. отпуск

[aka raduntsev-nias]

И это хорошо! Желаем набраться сил и быстрей к нам вернуться!

to 6pa4o/ Нашел видеоролик построения и расчета узла. Если интересно, скину в личку.

[derik]

Хорошего отдыха и приятных эмоций! угу ^_^

[aka raduntsev-nias]

эти ролики есть в сети. кому интерсно, пишите, скину ссылку

[aka raduntsev-nias]

Сообщение от aka raduntsev-nias 6pa4o, с рсу как то решил вопрос?

читай в хэлпе Creating and Combining Load Case - РСН

РСУ

[derik]

спасибо! я где о так и думал, но особенно не вдавался - с контактом разбираюсь :-)

Сергей Иванович, если у Вас уже закончился отпуск, не могли бы Вы посмотреть код в #156. Жду критики :-)

[СергейД]

посмотрите =исправил основное.
разбирайтесь.

[derik]

________

[soppr]

А чем открывать файл mac?

[derik]

любым текстовым редактором. scite пользуюсь

[aka raduntsev-nias]

Все признаки несходимости контактной задачи. ответ можно найти на соседнем форуме, передирать не буду, не я отвечал, да и тренируюсь пока в WB
может Сергей Иванович повторится.

[derik]

я в курсе что несходиться... :-) а что за соседний форум?

[aka raduntsev-nias]

Форум САПР

[СергейД]

форум fsapr2000.ru
но там что-то стало малоинтересно.
советы такие.
хотя все же что-то и забыл....
1. контакт узел в узел идеален и может быть реализован через тип 178
2. если сгущать, то только контакт. частая причина проблем=
когда таргет кое-где (и тем более, если много где) мельче , чем контакт.
Нужно понимать, что контакт проверяется в точках интегрирования (для прямоугольнкиа
их четыре и они по совету Гаусса в 29% от центра, если 100% размер линейный)
то есть если узел таргета попадает между точками интегрирования (ТИ) контактной поверхнгости-
он будет внедряться. и картина исказится.
3. главный критерий =величина внедрения и ее нужно всегда рисовать
etab,cont,pene. она должна быть раз в 100 минимум меньше перемещений соседних узлов.
4. контактная жесткость по умолчанию= слишком "жесткая".
вполне достаточно 0.1 и даже 0.01 для улучшения сходимости
если внимательно проверять после этого результаты по п.3
5. опция 10, 5 сглаживает сходимость
6. всегда стараюсь для отладки сначала успешно решить задачу при 12,5 bonded
6. Всегда применяю штрафной алгоритм penalty 2,1
7. нельзя решать без трения вообще.
но при реальном трении 0.2-0.3 необходим несимметричный решатель
nropt,unsym,,on
8. нужно решить с долей нагрузки (максимум 1%)
проанализировать результаты и убедиться в установлении контакта и физичности результатов
9. а уж потом приложить полную минимум за 10 шагов. или сразу ставьте Nsubst,1000,1000000,100
10. невязки часто допустимы много больше принимаемых по умолчанию CNVTOL,f,
11. всегда рисуйте мультики и графикик и записывайте промежуточные шаги.
12. обязательно весь вывод направляйте в текстовый файл. Там зачастую достаточно
подсказок в диагностике.
13. ПЕРЕЧИТАЙТЕ хелп прежде чем решать сколько нибудь новую задачу.
14. элемент можно считать плохим, если в нем неудачно расположились ТИ
скажем острые треугольники нежелательны. ТИ в них далеко от острого угла.
форма солидов здесь вторична, но чем они лучше, тем лучше для сходимотси
15. если задача не сошлась, нужно нарисовать картину распределения невязок.
(но их запись нужно заказывать nldiag,nrre,on)
иногда положение мест несходимости бывает очень неожиданным. (где-либо на границе контактной зоны)

[derik]

это все + к мануалу по контакту. Сергей Иванович, а как мне сгенерерить таргет и контакт возле отверстия в пластине (болтовое соединение, контакт головки болта с пластиной) - выбираю линии круга lsel,s,,,l1,l2,n
nsll,s,1
esln
esurf

ансис пишет, что нет поверхности для контакта...т.о. выделяются не все элементы пластины, которые контактируют с головкой болта. как вариант сделать очень грубую сетку пластины, но это исказит результат. выделять всю поверхность (я сейчас так и делаю, т.к. другова способа не нашел) расточительно - мощностя компьютера не безграничны...макрос прикладываю, просьба сильно не пинать и камни не кидать - он еще сыроват :-) первоначально попробывал решить четвертинку - тот же результат - не геерерятся элементы, только если выделить поверхность.

[СергейД]

друзья!
очень занят.

[aka raduntsev-nias]

попробую решить в WB. если разберусь с исходником в твоем mac файле. если мозгов хватит)

единственное, не пойму, как задать контакт узел в узел в WB для пластин!

[derik]

^_^ там все просто...разбей пластинки на одинаковое кол-во элементов. мб нужно пластинки, которые будут контактировать, вообще моделировать отдельно, а потом сшивать узлы. тут Сергей Иванович подсобил бы...

[aka raduntsev-nias]

Если сравнивать с лирой, то у тебя оси пластин не выровнены

[derik]

ээээ..нормали, вроде, направлены как по мануалу...а вот про рластины забыл начисто....чорт...ночью разберемся пасиб за "тропинку"

[aka raduntsev-nias]

)))ну и ж этот apdl
6 часов убил и мало что понял

и еще. если сравнивать с wb, то болтовое задаю в 2 загружениях.
1 load
2 lock

здесь пока не понял что.
пока вижу только одну силу. 1e3

[derik]

не, апдл это круть :-) Сергей Иванович прав. с ним проще...я уже мышкой практически не пользуюсь...ток чтоб покрутить модель диван...музыка спокойная, кофе с печенькой и бук на коленках...
тебе просто сначала нужно было апдл учить, потом команды сами запоминаются, я вообще был в шоке, как можно запонмить столько команд и столько параметров каждой команды...ниче запомнил и в хелп дажо лезть не нужно

нагрузка здесь одна sfa,1e3...на головку болта.самое простое

а здесь все через

type,4 и type,3
меняются только лишь реальные константы - r,10....14
выбор каждой контактной пары заканчивается esurf

[aka raduntsev-nias]

хорошая ссылочка

[derik]

тот же хелп ток забитый в инет...имхо, встроенная справка удобней - тут непревычно! :-) вот если бы она на русском была...

[aka raduntsev-nias]

вопрос. как избавиться от этих дряней.

наверное из за этого
SHPP,OFF

????

[derik]

это контактный(целевой элемнт) т ыот нео никак, по мойму не избавишься, могу ошибатся.

[aka raduntsev-nias]

не от контактного. нужно чтобы поверхность на месте оставалась.

[derik]

поверхности контакта и целевые поверхности должны деформироваться вместе с моделью...у тебя, судя по картинке, контактная поверхность остается на месте...

[aka raduntsev-nias]

я про то же. это при самом большом увеличении визуализации деформаций.
при среднем и малом все ok

[derik]

Нужна помощь Сергея Ивановича...а он занят...

[aka raduntsev-nias]

может ошибаюсь, но при трении 0,2 несимметричный контакт

[derik]

отдельно болтик решил :-) потом пластики, а потом все вместе :-) Макс можошь посмотреть в ВБ, результаты совпадут? :-)

[aka raduntsev-nias]

хорошо. если разберусь в этом коде . не успеваю сильно вникнуть, урывками apdl изучаю. на wb. с геометрией вроде ясно. и с нагузкой. с граничными условиями не совсем
и не пойму еще реультат в классике. вроде тоже выбираю, что и в вб, но цифры какие то другие. в общем пока.
не обещаю что скоро. еще свою первую не закончил чтобы совсем, но почти. с контактами только более менее разобрался. пропали тянучки наконец.
кстати болты сделал с трением. в разной литературе по разному с болтами.

скинешь данные на бумажке быстрей решу) признаюсь. я жесткий чайник.

хотя ладно, решу. потом сравним

тоже файлик. классик не решает. wb пропустил. не знаю пока, как эти ошибки в wb исправить

[derik]

если трудности с командами - спрашивай, что знаю расскажу а чем открыть присланный файл? ни ВБ ни ансис его не читают...а вот и сам болтик, что считаю. примитивный, но ведь с примитивов и нужно начинать

по поводу журнала ошибок:
pcg решатель не подходит, убери его, а в /sol встать следующие строчки:

Код:

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

/gst,1
nldiag,nrre,on
nropt,unsym,,on
antype,0
nlgeom,1
nsubst,100,5000,10
outres,erase
outres,all,all
time,10
solve
fini

попробуй сменить 186 на 187 тип
и что там там у тебя не то в задании таблиц материала, хотя стоит попробовать сменить тип элемента...может быть и ошибаюсь, пускай меня поправят старшие товарищи

[aka raduntsev-nias]

решатель видимо таким из WB передался. выбирал интерактивный, тк считал нессиметричный контакт. файл можно переименовать в .txt. или открывать как коммандный файл в ansys

сменить элементы тоже не понятно как в WB
вот так. тоже не все просто

Сообщение от 6pa4o ведь с примитивов и нужно начинать )

согласен. но сроки жмут. мне сдавать работу скоро руководителю.

[derik]

тада хз....считай в классике ИМХО! брось ВБ переходи на АНСИС!!!! даешь стране АПДЛ!!!!

[aka raduntsev-nias]

Сообщение от 6pa4o тада хз....считай в классике ИМХО! брось ВБ переходи на АНСИС!!!! даешь стране АПДЛ!!!!

обязательно! потихоньку перейду. к тебе буду обращаться!

[derik]

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


попробуй пустить на расчет

[aka raduntsev-nias]

Сообщение от 6pa4o попробуй пустить на расчет

заработало? спасиб! после работы проверю

[derik]

не я не знаю...счас ансиса нет под рукой, я там поправил тип элемента и pcg убрал...авось поможет :-D

Код:

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

fini
/clear,nostart
/outp,outfile,txt,C:\ansys\output
/cwd,'C:\ansys\out'

a=100
b=150
c=150
h=40
t=20
r=20
r2=35
p=5e4

/prep7

et,1,187
et,2,185

mptemp,,,,,,,
mptemp,1,0
mpdata,ex,1,,2.06e8
mpdata,prxy,1,,.3

cylind,r,,0,a,,90,
cylind,r,r2,0,h,,90,
cylind,r2,,0,h,,90,
vgen,,1,,,,,-a,,,1
vgen,,2,,,,,-2*h,,,1
vadd,1,3
blc5,b/2,c/2,b,c,-t
cylind,r,,0,-t,,90
vsbv,1,3
vgen,2,5,,,,,-t,,,0
sfa,13,1,pres,-p

mshkey,0
mshape,0,3d
esize,r/4
vsel,s,,,1,5,4
vsweep,all
vsel,inve
vmesh,all

et,3,170
et,4,174
keyopt,4,2,1
keyopt,4,5,1
keyopt,4,9,3
keyopt,4,10,5
keyopt,4,12,5

mp,mu,3,.3
mat,3

r,3,,,.1,.1,,40
r,4,,,.1,.1,,40
r,5,,,.1,.1,,40
r,6,,,.1,.1,,40
r,7,,,.1,.1,,40

real,3
type,3
asel,s,,,3
nsla,s,1
esln
esurf
type,4
asel,s,,,9
nsla,s,1
esln
esurf

real,4
type,3
asel,s,,,2,24,22
nsla,s,1
esln
esurf
type,4
asel,s,,,3
nsla,s,1
esln
esurf

real,5
type,3
asel,s,,,15
nsla,s,1
esln
esurf
type,4
asel,s,,,29
nsla,s,1
esln
esurf

real,6
type,3
asel,s,,,30
nsla,s,1
esln
esurf
type,4
asel,s,,,17
nsla,s,1
esln
esurf

real,7
type,3
asel,s,,,22
nsla,s,1
esln
esurf
type,4
asel,s,,,7
nsla,s,1
esln
esurf
cncheck,auto

asel,s,,,4,5
asel,a,,,10,11
asel,a,,,18,19
asel,a,,,12,20,8
asel,a,,,27,28
da,all,symm
asel,s,,,16,21,5
asel,a,,,23,25,2
da,all,all
!da,13,uz,1.8
allsel,all
fini

/sol

/gst,1
nldiag,nrre,on
nropt,unsym,,on
antype,0
nlgeom,1
nsubst,100,10000,50
outres,erase
outres,all,all
time,5
solve
fini

/post1
set,first

и все таки я ее решил :-D
2 гектара занимает решение...

[aka raduntsev-nias]

Сообщение от 6pa4o;620780 /post1 set,first[/code и все таки я ее решил :-D 2 гектара занимает решение...

видимо это не предел.
у меня сейчас получился 38,1 Гб

[derik]

ты шо там такое решал?

[aka raduntsev-nias]

Сообщение от 6pa4o ты шо там такое решал?

контактная. "узелок" посчитал) помельче сетку побил)

[derik]

Ы...а скинь файлик охота глянуть

[aka raduntsev-nias]

позже выложу.
после реальных испытаний.
лучше твой сегодня завтра добью и скину

[Cfytrr]

Сообщение от aka raduntsev-nias после реальных испытаний

Вы случаем не фланец будете испытывать?

[aka raduntsev-nias]

Cfytrr, да, выше выкладывал скрин) единственное, не совсем такой. начинаю работать с узлами стержневых конструкций. Пока обрабатываю старые наработки института(оформленные изобретения), где числюсь аспирантом, если есть какие -то предложения, замечания, поделитесь.

[Cfytrr]

aka raduntsev-nias, Тут в одной теме долго и с увлечением рассуждали о работе фланцевых соединений но без доказательной базы в виде реальных экспериментов разговор как то увял.
Хотелось бы получить исходные данные по проверяемому узлу и результаты испытаний, тогда параллельно с ANSYS можно проверить и в Femap (Nstran NX)

[aka raduntsev-nias]

Cfytrr, хорошо - так даже качественнее результат получится. Я свяжусь с Вами (скину геометрию в step_iges), но позже немного
Кстати я немного(совсем пока немного) освоил Siemens NX UG(чем то похож на Tekla, отличная программа), использую иногда для связки с ansys.
Проверил в ansys некоторые фланцевые соединения, что можно найти в примерах учебника Кузина Н.Я. "Проектирование и расчет стальных ферм", результат очень даже совпадает.

в WB
в классике у тебя скорей всего эту геометрию править придется
это пример из Кузина, если кто непонял

[derik]

Макс, а ты это в ВБ решаешь или классике?

скинь геометрию в акаде или номер страницы кузина...шас посморим

[aka raduntsev-nias]

Сообщение от 6pa4o Макс, а ты это в ВБ решаешь или классике? скинь геометрию в акаде или номер страницы кузина...шас посморим

стр 150 файл выше

[derik]

см. # 217

[aka raduntsev-nias]

6pa4o,
классно!
повешу код на стену)
только ребра можно было оборвать на торцевой пластине и отверстие в ней не нужно, видимо я геометрию с ошибкой кинул
и сетка не везде связная у ребер с круглой пластиной.
те
убрать отверстие в центре круглой, так, чтобы трубы доходили до нее
ребра одовести до круглой, дальше не надо, у тебя пока ее пересекают
ну и сетка, что то не то
а в общем классно. можно потом провети оптимизацию.
весь мой узелок уместился на листочек)

[derik]

Исправил. Макс забирай с поста выше. контактные поферхости и сделаю..но потом, сейчас некада...тебе это вчера нужно было?

[aka raduntsev-nias]

6pa4o, спасибо! да - это.

[Cfytrr]

Сообщение от aka raduntsev-nias ...это пример из Кузина

Так Вы именно этот фланец будете испытывать?
Если да, то какая нагрузка в расчетной схеме, есть ли предварительное натяжение, а то из фрагмента примера как то не очень понятно...

[derik]

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

Макс, нет, нужно пересекать объемы...но там контактных пар будет много...эх...прикольно...за ночь можно все сделать! нужно будет эту ночь выкроить.

[aka raduntsev-nias]

Cfytrr, Да нет, не этот. Но он ничем не отличается от всех подобных фланцевых. Ведь принцип решения тот-же.
В этом узле есть предварительное натяжение. В болтах. Болты М24
To Брачо
наверное скоро будешь

[Cfytrr]

Сообщение от aka raduntsev-nias Да нет, не этот. Но он ничем не отличается от всех подобных...

Тогда может быть "моделить" именно то, что будут испытывать? Принципы могут быть и одинаковыми, но количественно решения будут отличаться...

[aka raduntsev-nias]

я уже сделал тот, который надо. (я умею быстро моделить, поэтому сделать геометрию почти любого узла для меня не составляет трудностей)
но пока не могу его выложить, тк он Разработка. не моя даже полностью. сами понимаете в общем.
работать над ним дальше нужно. как только смогу, выложу.
Результаы расчета в МКЭ для этого, если хотите, можно сравнить.
Данные в учебнике.

[derik]

Cfytrr а картинка из какой прграммы? результаты можно расчета?

Код:

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

fini
/clear,nostart
/outp,outfile,txt,C:\ansys\output
/cwd,'C:\ansys\out'

*set,a,.8
*set,b,.2
*set,t1,.05
*set,t2,.01
*set,r1,.5
*set,r2,.19
*set,r3,.04
*set,r4,.06
*set,p,5e4

/prep7

et,1,95

mptemp,,,,,,,
mptemp,1,0
mpdata,ex,1,,2.06e8
mpdata,prxy,1,,.3

cyl4,,,r2,0,r2+t2,360,a+t1
k,,-r1,
k,,r1,
k,,r1,,t1+t1/2
k,,r2+t1*1.5,,a+t1
k,,-r2-t1*1.5,,a+t1
k,,-r1,,t1+t1/2
a,17,18,19,20,21,22,17
vext,7,,,,t2,,,
vgen,,2,,,,-t2/2,,,,1
csys,1
vgen,2,2,,,,90,,,0
csys,0
vptn,all
vsel,s,,,9,10
vsel,a,,,13,14
vsel,a,,,4
vdele,all,,,1
asel,s,,,30,48,6
asel,a,,,64,67,3
asel,a,,,102,104,2
vext,all,,,,,a/3,,,
larc,29,18,73,r1
al,185,169,161,105,
vext,109,,,,,t1,,,
cyl4,r1/2,-r1/2,r3,0,0,360,t1
vsbv,21,22
csys,1
vgen,4,23,,,,90,,,0
csys,0
cyl4,,,r2,0,,360,t1
allsel,all
nummrg,kp
vsymm,z,all,,,,0,0
cyl4,r1/2,r1/2,r3,0,0,360,t1*2
vgen,,51,,,,,-t1,,,1
cyl4,r1/2,r1/2,r4,0,0,360,t1*1.5
vgen,,52,,,,,t1,,,1
cyl4,r1/2,r1/2,r4,0,0,360,t1
vgen,,53,,,,,-t1*2,,,1
vsel,s,,,51,53,1
csys,1
vgen,4,all,,,,90,,,0
csys,0
allsel,all

mshkey,1
mshape,0,3d
esize,r3/2
vsweep,all

воть, теперь все как нада
осталось наложить контактные пары и ограничить узел! считайте, мой комп расчет не потянет, слабоват для такого

[aka raduntsev-nias]

несовсем. опять фасонками стенку трубы пробил

[derik]

нет, подругому получится, что элементы несвязаны будут с собой...а так получилось, что это единый узел( верхний и нижхий) скрепленый только болтовым соединением.

а вот и полный макрос
пусть если что неправильно поправте, плиз...я только учусь

Код:

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

fini
/clear,nostart
/outp,outfile,txt,C:\ansys\output
/cwd,'C:\ansys\out'

*set,a,.8
*set,b,.2
*set,t1,.05
*set,t2,.01
*set,r1,.5
*set,r2,.19
*set,r3,.04
*set,r4,.06
*set,p,5e4

/prep7

et,1,95

mptemp,,,,,,,
mptemp,1,0
mpdata,ex,1,,2.06e8
mpdata,prxy,1,,.3

cyl4,,,r2,0,r2+t2,360,a+t1
k,,-r1,
k,,r1,
k,,r1,,t1+t1/2
k,,r2+t1*1.5,,a+t1
k,,-r2-t1*1.5,,a+t1
k,,-r1,,t1+t1/2
a,17,18,19,20,21,22,17
vext,7,,,,t2,,,
vgen,,2,,,,-t2/2,,,,1
csys,1
vgen,2,2,,,,90,,,0
csys,0
vptn,all
vsel,s,,,9,10
vsel,a,,,13,14
vsel,a,,,4
vdele,all,,,1
asel,s,,,30,48,6
asel,a,,,64,67,3
asel,a,,,102,104,2
vext,all,,,,,a/3,,,
larc,29,18,73,r1
al,185,169,161,105,
vext,109,,,,,t1,,,
cyl4,r1/2,-r1/2,r3,0,0,360,t1
vsbv,21,22
csys,1
vgen,4,23,,,,90,,,0
csys,0
cyl4,,,r2,0,,360,t1
allsel,all
nummrg,kp
vsymm,z,all,,,,0,0
cyl4,r1/2,r1/2,r3,0,0,360,t1*2
vgen,,51,,,,,-t1,,,1
cyl4,r1/2,r1/2,r4,0,0,360,t1*1.5
vgen,,52,,,,,t1,,,1
cyl4,r1/2,r1/2,r4,0,0,360,t1
vgen,,53,,,,,-t1*2,,,1
vsel,s,,,51,53,1
csys,1
vgen,4,all,,,,90,,,0
csys,0
allsel,all

mshkey,1
mshape,0,3d
esize,r3/2
vsweep,all

et,2,170
et,3,174
keyopt,3,2,1
keyopt,3,5,1
keyopt,3,9,3
keyopt,3,10,5
keyopt,3,12,5
r,3,,,.1,.1,,10
r,4,,,.1,.1,,10
r,5,,,.1,.1,,10
r,6,,,.1,.1,,10
r,7,,,.1,.1,,10
mp,mu,3,.2

mat,3

!фланцы
real,3
type,2
vsel,s,loc,z,0,a
vsel,u,,,51,60,3
aslv
asel,r,loc,z,0
nsla,s,1
esln
esurf
type,3
vsel,s,loc,z,0,-a
vsel,u,,,51,60,3
aslv
asel,r,loc,z,0
nsla,s,1
esln
esurf

!болтовое соединение
!лев. нижний
real,4
type,2
asel,s,,,133,134
asel,a,,,277,278
nsla,s,1
esln
esurf
type,3
vsel,s,,,57
aslv
asel,r,loc,z,0
nsla,s,1
esln
esurf

!прав. нижний
real,5
type,2
asel,s,,,117,118
asel,a,,,269,270
nsla,s,1
esln
esurf
type,3
vsel,s,,,60
aslv
asel,r,loc,z,0
nsla,s,1
esln
esurf

!прав. верхний
real,6
type,2
asel,s,,,109,110
asel,a,,,253,254
nsla,s,1
esln
esurf
type,3
vsel,s,,,51
aslv
asel,r,loc,z,0
nsla,s,1
esln
esurf

!лев. верхний
real,7
type,2
asel,s,,,125,126
asel,a,,,261,262
nsla,s,1
esln
esurf
type,3
vsel,s,,,54
aslv
asel,r,loc,z,0
nsla,s,1
esln
esurf
cncheck,auto

asel,s,loc,z,t1+t1*1.5
da,all,all
asel,s,loc,z,a+t1+a/3
sfa,all,1,pres,-p
asel,s,loc,z,-(a+t1+a/3)
sfa,all,2,pres,-p
allsel,all
fini

/sol

/gst,1
nldiag,nrre,on
nropt,unsym,,on
antype,0
nlgeom,1
nsubst,100,10000,50
deltim,0.1,0.1,0.1
outres,erase
outres,all,all
time,5
solve
fini

/post1
set,first

[aka raduntsev-nias]

для нессиметр контакта нужен комп-р неслабый. мой с 8 Ггб оперативки 3Ггц процессор 2 ядра и то тормозит
по моему ты болты ненапряг.

[derik]

я не умею этго делать, пока это следующий шаг...

[aka raduntsev-nias]

В ANSYS сечение предварительного натяжения задано командой SECDATA (используя опцию типа PRETENSION) во входном файле ANSYS. Это задает часть модели, через которую программа применяет предварительное нагружение. Сечение предварительного натяжения состоит:

Элементы PRETS179, которые задают соединения или соединительные детали непосредственно.

Узел "предварительного натяжения", который программа использует для управления и контроля полных растягивающих нагрузок.

В ANSYS нагрузки применяются к сечению предварительного натяжения командой SLOAD.

[derik]

либо через psmesh...проще и быстрее только нужно посмотреть как она будет себя вести на контактной поверхности...чую, psmesh и esurf не дружат меж собой

[soppr]

Два брата-акробата (Raduntsev и bra4o), вы куда пропали?))
Темка интересная.

[aka raduntsev-nias]

добрый день.
немного изменил код из хелпа ansys

Код:

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

/prep7   
/title,Sample application of PSMESH 
 
et,1,92
mp,ex,1,3e7  
mp,alpx,1,8.4e-6 
mp,prxy,1,0.30 
tref,70  
 
/foc,,-.09,.34,.42  
/dist,,.99 
/ang,,-55.8  
/view,,.39,-.87,.31 
/pnum,volu,1  
/num,1 
cylind,0.5,,  -0.25,0,  0,180 
cylind,0.5,,  1,1.25,   0,180 
cylind,0.25,, 0,1,      0,180 
vglue,all  
numc,all  
vplot 
mat,1  
smrt,off  
vmesh,1,3
/pnum,mat,1  
eplot 
psmesh,,example,,volu,1,0,z,0.5,,,,elems    
CM,lines,LINE  
/dist,,1.1  
cmplot 
/solu 
eqslve,pcg,1e-8 
asel,s,loc,y
da,all,symm  
asel,all 
dа,1,ux  
da,3,ux  
da,1,uz 
da,3,uz
da,1,uy  
da,3,uy 
sload,1,9,,force,100,1,2 
/title,Sample application of PSMESH - preload only   
solve

при решении ansys ругается на команду da
что в ней не так?
сравнил решение с abaqus. возможно такая разница из-за ограничения студенческой версии. результат меняется при изменении сетки

to soppr: дел больше стало . Наверное и у Брачо тоже

[derik]

а попробуй:

asel,s,,,1,3,2
da,all,all

вместо

dа,1,ux
da,3,ux
da,1,uz
da,3,uz
da,1,uy
da,3,uy

З,Ы,:а еще лучше не привязываться к номерам поверхностей, а выбирать через asel,s,loc,z(x, y), параметр выбора(высота, ширина, длина)

[aka raduntsev-nias]

пошло дело. результат только сильно поменялся.

[derik]

ансис от абакуса?

[aka raduntsev-nias]

Сообщение от 6pa4o ансис от абакуса?

ну да. в abaqus то же задать и все нормально будет. но дальше 1000 узлов непрыгнешь, поэтому все равно разница будет.
кстати, появилось другое предупреждение. когда открываешь результаты.
это мой первый код. видимо еще где-то напортачил.

[aka raduntsev-nias]

mp,ex,1,1
mp,prxy,1,0.01
поставим нули и получим то же практически.

abaqus (у меня) вообще нереагирует на изменения материалов(в самой простой задаче из хелпа ansys). при разных характеристиках ответ один. неправильный.


У кого-нибудь есть какие-то представления о связи Fortran Compiler и Abaqus. Например установлен Abaqus на моем компьютере, а также Intel Visual Fortran. Как запустить UMAT пользователем подпрограммы? Как правильно связать compier и Abaqus?

как подключить material library?

[derik]

а попробуй:
mptemp,,,,,,,
mptemp,1,0
mpdata,ex,1,,3e6
mpdata,prxy,1,,.2

[shtoss]

Добрый день всем.
Вопрос от чайника
Где в ANSYS Mechanical команда отменить/вернуть т.е команда undo всю программу облазил не могу найти пробовал в командной строке вводить UNDO /UNDO,ON ни какой реакции ansys. А то как странно в воркбенче есть, а в обычном нет в TOOLBAR. А то очень не удобно приходится каждый раз файл сохранять.

[СергейД]

есть способ гибче.
в session editor (в меню)
изменяете Любые команды
оставшиеся после save

/undo,on нужно заказывать заранее.
если работаете макросами такая проблема почти неактуальна.

[iliorik]

Добрый день!
Я вижу вы тут самые рулевые по Ансис.
Подкажите пожалуста как здесь все колонны мне симметрично скопировать, и возможно ли верх колонн с контуром жестко завязать (что б он никуда не перемещался)?

Код:

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

/FILNAME,PIONERSK_70
/PREP7
ET,1,SHELL63
 
          !Первая плита
          k,1,-0.0500,13.5000, 0.0000
          k,2, 3.4500,13.5000, 0.0000
          k,3, 3.4500,15.0000, 0.0000
          k,4,13.9000,15.0000, 0.0000
          k,5,13.9000, 7.4500, 0.0000
          k,6,21.0000, 7.4500, 0.0000
          k,7,21.0000,15.0000, 0.0000
          k,8,31.4500,15.0000, 0.0000
          k,9,31.4500,13.5000, 0.0000
          k,10,34.9500,13.5000, 0.0000
          k,11,34.9500, 1.5000, 0.0000
          k,12,27.9000, 1.5000, 0.0000
          k,13,27.9000, 0.0000, 0.0000
          k,14, 7.0000, 0.0000, 0.0000
          k,15, 7.0000, 1.5000, 0.0000
          k,16,-0.0500, 1.5000, 0.0000
          a,1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16
          !Плита корридора
          k,17,13.9000,16.2000, 0.0000
          k,18,18.0500,16.2000, 0.0000
          k,19,18.0500,14.6500, 0.0000
          k,20,16.0000,14.6500, 0.0000
          k,21,16.0000, 9.6250, 0.0000
          k,22,18.0500, 9.6250, 0.0000
          k,23,18.0500,12.0900, 0.0000
          k,24,21.0000,12.0900, 0.0000
          k,25,21.0000,10.6200, 0.0000
          k,26,20.1310,10.6200, 0.0000
          k,27,20.1310, 9.2500, 0.0000
          k,28,21.0000, 9.2500, 0.0000
          k,29,21.0000, 7.4500, 0.0000
          k,30,13.9000, 7.4500, 0.0000
          a,17,18,19,20,21,22,23,24,25,26,27,28,29,30
!Колонны
!1
COL_X=-0.05                        !X координата верха колонны
COL_Y=1.75                        !Y координата верха колонны
COL_Z=0                            !Z координата верха колонны
COL_L=3                            !Длина колонны
COL_B=0.3                        !Ширина сечения колонны
COL_H=0.8                          !Высота сечения колонны
BLC5,COL_X,COL_Y,COL_B,COL_H    !Рисуем проекцию колонны (верх)
K,,COL_X,COL_Y,COL_Z            !Точка верха колонны
K,,COL_X,COL_Y,COL_Z-COL_L        !Точка низа колонны
*GET,MAXKP,kp,,num,max        !Находим последнее значение точки
LNODE=MAXKP-1
L,LNODE,MAXKP                    !Рисуем колонну
*GET,MAXAREA,area,,num,max        !Находим последнее значение плиты
ASBA,ALL,MAXAREA,,,KEEP            !Вырезаем сечение колонны
NUMCMP,AREA                        !Перенумерация поверхностей
 
!2
COL_X=-0.05                        !X координата верха колонны
COL_Y=1.75+5.75                        !Y координата верха колонны
COL_Z=0                            !Z координата верха колонны
COL_L=3                            !Длина колонны
COL_B=0.3                        !Ширина сечения колонны
COL_H=0.8                          !Высота сечения колонны
BLC5,COL_X,COL_Y,COL_B,COL_H    !Рисуем проекцию колонны (верх)
K,,COL_X,COL_Y,COL_Z            !Точка верха колонны
K,,COL_X,COL_Y,COL_Z-COL_L        !Точка низа колонны
*GET,MAXKP,kp,,num,max        !Находим последнее значение точки
LNODE=MAXKP-1
L,LNODE,MAXKP                    !Рисуем колонну
*GET,MAXAREA,area,,num,max        !Находим последнее значение плиты
ASBA,ALL,MAXAREA,,,KEEP            !Вырезаем сечение колонны
NUMCMP,AREA                        !Перенумерация поверхностей
 
!3
COL_X=-0.05                        !X координата верха колонны
COL_Y=1.75+5.75*2                !Y координата верха колонны
COL_Z=0                            !Z координата верха колонны
COL_L=3                            !Длина колонны
COL_B=0.3                        !Ширина сечения колонны
COL_H=0.8                          !Высота сечения колонны
BLC5,COL_X,COL_Y,COL_B,COL_H    !Рисуем проекцию колонны (верх)
K,,COL_X,COL_Y,COL_Z            !Точка верха колонны
K,,COL_X,COL_Y,COL_Z-COL_L        !Точка низа колонны
*GET,MAXKP,kp,,num,max        !Находим последнее значение точки
LNODE=MAXKP-1
L,LNODE,MAXKP                    !Рисуем колонну
*GET,MAXAREA,area,,num,max        !Находим последнее значение плиты
ASBA,ALL,MAXAREA,,,KEEP            !Вырезаем сечение колонны
NUMCMP,AREA                        !Перенумерация поверхностей
 
!4
COL_X=3.5                       !X координата верха колонны
COL_Y=1.5                         !Y координата верха колонны
COL_Z=0                            !Z координата верха колонны
COL_L=3                            !Длина колонны
COL_B=0.8                        !Ширина сечения колонны
COL_H=0.3                          !Высота сечения колонны
BLC5,COL_X,COL_Y,COL_B,COL_H    !Рисуем проекцию колонны (верх)
K,,COL_X,COL_Y,COL_Z            !Точка верха колонны
K,,COL_X,COL_Y,COL_Z-COL_L        !Точка низа колонны
*GET,MAXKP,kp,,num,max        !Находим последнее значение точки
LNODE=MAXKP-1
L,LNODE,MAXKP                    !Рисуем колонну
*GET,MAXAREA,area,,num,max        !Находим последнее значение плиты
ASBA,ALL,MAXAREA,,,KEEP            !Вырезаем сечение колонны
NUMCMP,AREA                        !Перенумерация поверхностей
 
!5
COL_X=10.2                       !X координата верха колонны
COL_Y=0                         !Y координата верха колонны
COL_Z=0                            !Z координата верха колонны
COL_L=3                            !Длина колонны
COL_B=0.8                        !Ширина сечения колонны
COL_H=0.3                          !Высота сечения колонны
BLC5,COL_X,COL_Y,COL_B,COL_H    !Рисуем проекцию колонны (верх)
K,,COL_X,COL_Y,COL_Z            !Точка верха колонны
K,,COL_X,COL_Y,COL_Z-COL_L        !Точка низа колонны
*GET,MAXKP,kp,,num,max        !Находим последнее значение точки
LNODE=MAXKP-1
L,LNODE,MAXKP                    !Рисуем колонну
*GET,MAXAREA,area,,num,max        !Находим последнее значение плиты
ASBA,ALL,MAXAREA,,,KEEP            !Вырезаем сечение колонны
NUMCMP,AREA                        !Перенумерация поверхностей
 
!6
COL_X=13.7                       !X координата верха колонны
COL_Y=0                         !Y координата верха колонны
COL_Z=0                            !Z координата верха колонны
COL_L=3                            !Длина колонны
COL_B=0.8                        !Ширина сечения колонны
COL_H=0.3                          !Высота сечения колонны
BLC5,COL_X,COL_Y,COL_B,COL_H    !Рисуем проекцию колонны (верх)
K,,COL_X,COL_Y,COL_Z            !Точка верха колонны
K,,COL_X,COL_Y,COL_Z-COL_L        !Точка низа колонны
*GET,MAXKP,kp,,num,max        !Находим последнее значение точки
LNODE=MAXKP-1
L,LNODE,MAXKP                    !Рисуем колонну
*GET,MAXAREA,area,,num,max        !Находим последнее значение плиты
ASBA,ALL,MAXAREA,,,KEEP            !Вырезаем сечение колонны
NUMCMP,AREA                        !Перенумерация поверхностей
 
!7
COL_X=10.45                       !X координата верха колонны
COL_Y=1.5                         !Y координата верха колонны
COL_Z=0                            !Z координата верха колонны
COL_L=3                            !Длина колонны
COL_B=0.3                        !Ширина сечения колонны
COL_H=0.3                          !Высота сечения колонны
BLC5,COL_X,COL_Y,COL_B,COL_H    !Рисуем проекцию колонны (верх)
K,,COL_X,COL_Y,COL_Z            !Точка верха колонны
K,,COL_X,COL_Y,COL_Z-COL_L        !Точка низа колонны
*GET,MAXKP,kp,,num,max        !Находим последнее значение точки
LNODE=MAXKP-1
L,LNODE,MAXKP                    !Рисуем колонну
*GET,MAXAREA,area,,num,max        !Находим последнее значение плиты
ASBA,ALL,MAXAREA,,,KEEP            !Вырезаем сечение колонны
NUMCMP,AREA                        !Перенумерация поверхностей
 
!8
COL_X=13.95                       !X координата верха колонны
COL_Y=1.5                         !Y координата верха колонны
COL_Z=0                            !Z координата верха колонны
COL_L=3                            !Длина колонны
COL_B=0.3                        !Ширина сечения колонны
COL_H=0.3                          !Высота сечения колонны
BLC5,COL_X,COL_Y,COL_B,COL_H    !Рисуем проекцию колонны (верх)
K,,COL_X,COL_Y,COL_Z            !Точка верха колонны
K,,COL_X,COL_Y,COL_Z-COL_L        !Точка низа колонны
*GET,MAXKP,kp,,num,max        !Находим последнее значение точки
LNODE=MAXKP-1
L,LNODE,MAXKP                    !Рисуем колонну
*GET,MAXAREA,area,,num,max        !Находим последнее значение плиты
ASBA,ALL,MAXAREA,,,KEEP            !Вырезаем сечение колонны
NUMCMP,AREA                        !Перенумерация поверхностей
 
!9
COL_X=3.7                       !X координата верха колонны
COL_Y=7.5                         !Y координата верха колонны
COL_Z=0                            !Z координата верха колонны
COL_L=3                            !Длина колонны
COL_B=0.8                        !Ширина сечения колонны
COL_H=0.3                          !Высота сечения колонны
BLC5,COL_X,COL_Y,COL_B,COL_H    !Рисуем проекцию колонны (верх)
K,,COL_X,COL_Y,COL_Z            !Точка верха колонны
K,,COL_X,COL_Y,COL_Z-COL_L        !Точка низа колонны
*GET,MAXKP,kp,,num,max        !Находим последнее значение точки
LNODE=MAXKP-1
L,LNODE,MAXKP                    !Рисуем колонну
*GET,MAXAREA,area,,num,max        !Находим последнее значение плиты
ASBA,ALL,MAXAREA,,,KEEP            !Вырезаем сечение колонны
NUMCMP,AREA                        !Перенумерация поверхностей
 
!10
COL_X=7.2                        !X координата верха колонны
COL_Y=7.5                         !Y координата верха колонны
COL_Z=0                            !Z координата верха колонны
COL_L=3                            !Длина колонны
COL_B=0.8                        !Ширина сечения колонны
COL_H=0.3                          !Высота сечения колонны
BLC5,COL_X,COL_Y,COL_B,COL_H    !Рисуем проекцию колонны (верх)
K,,COL_X,COL_Y,COL_Z            !Точка верха колонны
K,,COL_X,COL_Y,COL_Z-COL_L        !Точка низа колонны
*GET,MAXKP,kp,,num,max        !Находим последнее значение точки
LNODE=MAXKP-1
L,LNODE,MAXKP                    !Рисуем колонну
*GET,MAXAREA,area,,num,max        !Находим последнее значение плиты
ASBA,ALL,MAXAREA,,,KEEP            !Вырезаем сечение колонны
NUMCMP,AREA                        !Перенумерация поверхностей
 
!11
COL_X=13.95                        !X координата верха колонны
COL_Y=7.5                         !Y координата верха колонны
COL_Z=0                            !Z координата верха колонны
COL_L=3                            !Длина колонны
COL_B=0.3                        !Ширина сечения колонны
COL_H=0.3                          !Высота сечения колонны
BLC5,COL_X,COL_Y,COL_B,COL_H    !Рисуем проекцию колонны (верх)
K,,COL_X,COL_Y,COL_Z            !Точка верха колонны
K,,COL_X,COL_Y,COL_Z-COL_L        !Точка низа колонны
*GET,MAXKP,kp,,num,max        !Находим последнее значение точки
LNODE=MAXKP-1
L,LNODE,MAXKP                    !Рисуем колонну
*GET,MAXAREA,area,,num,max        !Находим последнее значение плиты
ASBA,ALL,MAXAREA,,,KEEP            !Вырезаем сечение колонны
NUMCMP,AREA                        !Перенумерация поверхностей
 
!12
COL_X=10.95                       !X координата верха колонны
COL_Y=7.5                         !Y координата верха колонны
COL_Z=0                            !Z координата верха колонны
COL_L=3                            !Длина колонны
COL_B=0.3                        !Ширина сечения колонны
COL_H=0.3                          !Высота сечения колонны
BLC5,COL_X,COL_Y,COL_B,COL_H    !Рисуем проекцию колонны (верх)
K,,COL_X,COL_Y,COL_Z            !Точка верха колонны
K,,COL_X,COL_Y,COL_Z-COL_L        !Точка низа колонны
*GET,MAXKP,kp,,num,max        !Находим последнее значение точки
LNODE=MAXKP-1
L,LNODE,MAXKP                    !Рисуем колонну
*GET,MAXAREA,area,,num,max        !Находим последнее значение плиты
ASBA,ALL,MAXAREA,,,KEEP            !Вырезаем сечение колонны
NUMCMP,AREA                        !Перенумерация поверхностей
 
!13
COL_X=10.45                        !X координата верха колонны
COL_Y=9.3                        !Y координата верха колонны
COL_Z=0                            !Z координата верха колонны
COL_L=3                            !Длина колонны
COL_B=0.3                        !Ширина сечения колонны
COL_H=0.3                          !Высота сечения колонны
BLC5,COL_X,COL_Y,COL_B,COL_H    !Рисуем проекцию колонны (верх)
K,,COL_X,COL_Y,COL_Z            !Точка верха колонны
K,,COL_X,COL_Y,COL_Z-COL_L        !Точка низа колонны
*GET,MAXKP,kp,,num,max        !Находим последнее значение точки
LNODE=MAXKP-1
L,LNODE,MAXKP                    !Рисуем колонну
*GET,MAXAREA,area,,num,max        !Находим последнее значение плиты
ASBA,ALL,MAXAREA,,,KEEP            !Вырезаем сечение колонны
NUMCMP,AREA                        !Перенумерация поверхностей
 
!14
COL_X=6.95                       !X координата верха колонны
COL_Y=9.3                         !Y координата верха колонны
COL_Z=0                            !Z координата верха колонны
COL_L=3                            !Длина колонны
COL_B=0.3                        !Ширина сечения колонны
COL_H=0.3                          !Высота сечения колонны
BLC5,COL_X,COL_Y,COL_B,COL_H    !Рисуем проекцию колонны (верх)
K,,COL_X,COL_Y,COL_Z            !Точка верха колонны
K,,COL_X,COL_Y,COL_Z-COL_L        !Точка низа колонны
*GET,MAXKP,kp,,num,max        !Находим последнее значение точки
LNODE=MAXKP-1
L,LNODE,MAXKP                    !Рисуем колонну
*GET,MAXAREA,area,,num,max        !Находим последнее значение плиты
ASBA,ALL,MAXAREA,,,KEEP            !Вырезаем сечение колонны
NUMCMP,AREA                        !Перенумерация поверхностей
 
!15
COL_X=6.95                       !X координата верха колонны
COL_Y=14.75                         !Y координата верха колонны
COL_Z=0                            !Z координата верха колонны
COL_L=3                            !Длина колонны
COL_B=0.3                        !Ширина сечения колонны
COL_H=0.8                          !Высота сечения колонны
BLC5,COL_X,COL_Y,COL_B,COL_H    !Рисуем проекцию колонны (верх)
K,,COL_X,COL_Y,COL_Z            !Точка верха колонны
K,,COL_X,COL_Y,COL_Z-COL_L        !Точка низа колонны
*GET,MAXKP,kp,,num,max        !Находим последнее значение точки
LNODE=MAXKP-1
L,LNODE,MAXKP                    !Рисуем колонну
*GET,MAXAREA,area,,num,max        !Находим последнее значение плиты
ASBA,ALL,MAXAREA,,,KEEP            !Вырезаем сечение колонны
NUMCMP,AREA                        !Перенумерация поверхностей
 
!16
COL_X=10.45                        !X координата верха колонны
COL_Y=15.0                        !Y координата верха колонны
COL_Z=0                            !Z координата верха колонны
COL_L=3                            !Длина колонны
COL_B=0.3                        !Ширина сечения колонны
COL_H=0.3                          !Высота сечения колонны
BLC5,COL_X,COL_Y,COL_B,COL_H    !Рисуем проекцию колонны (верх)
K,,COL_X,COL_Y,COL_Z            !Точка верха колонны
K,,COL_X,COL_Y,COL_Z-COL_L        !Точка низа колонны
*GET,MAXKP,kp,,num,max        !Находим последнее значение точки
LNODE=MAXKP-1
L,LNODE,MAXKP                    !Рисуем колонну
*GET,MAXAREA,area,,num,max        !Находим последнее значение плиты
ASBA,ALL,MAXAREA,,,KEEP            !Вырезаем сечение колонны
NUMCMP,AREA                        !Перенумерация поверхностей
 
!17
COL_X=10.45                       !X координата верха колонны
COL_Y=16.2                         !Y координата верха колонны
COL_Z=0                            !Z координата верха колонны
COL_L=3                            !Длина колонны
COL_B=0.3                        !Ширина сечения колонны
COL_H=0.3                          !Высота сечения колонны
BLC5,COL_X,COL_Y,COL_B,COL_H    !Рисуем проекцию колонны (верх)
K,,COL_X,COL_Y,COL_Z            !Точка верха колонны
K,,COL_X,COL_Y,COL_Z-COL_L        !Точка низа колонны
*GET,MAXKP,kp,,num,max        !Находим последнее значение точки
LNODE=MAXKP-1
L,LNODE,MAXKP                    !Рисуем колонну
*GET,MAXAREA,area,,num,max        !Находим последнее значение плиты
ASBA,ALL,MAXAREA,,,KEEP            !Вырезаем сечение колонны
NUMCMP,AREA                        !Перенумерация поверхностей

[СергейД]

некогда все читать...

vgen agen lgen - копирование
arsym vsymm lsymm - симметричное отражение
хелп открывал?

про жестко связать- непонятно. способов много

[iliorik]

Спасибоза ответы. Уже разберусь. Хелп открывал, но там как то туго искать все это.

Надо на подобии как в лире жесткое тело сделать. Только в лире это по точкам делается. А здесь есть ли вариант грани элемента жестко связать с концом колонны, или хотя бы точки на этих гранях?

Как вижу в данной теме человек перешел из 2д к 3д сеткам. В связи с этим встает вопрос - возможно ли как то будет в лиру перекинуть все это как 2д сетку? Если да то тогда половина запарок сразу ненужна.
Что думаете по этому поводу?

[derik]

Нужно колонны моделировать конечным элементом элементом beam, а сечение через sectype, secdata. жесткое закрепление - все пересечь - lptn,all...а вообще стоит почитать эту тему с перывых сообщений. Сергей Иванович все это уже рассказывал!
не пойму причем здесь лира? вы хотите результаты полученные в ансис проверить в лире?!

[iliorik]

Да. Хотелось бы в лиру перегонять потом все это и там решать. Хотя бы потому что лира как произведственная, а ансис лишь студенческая есть. Но ансис сетки хорошие строит, лире пока до этого еще расти и расти. И вот хотелось бы строить сетку ансис и передавать ее в лиру.

[aka raduntsev-nias]

Для этого конвертер нужен.
Вообще лучше сначала в лира, потом еще раз посчитать в ansys

[derik]

а можно сразу считать в ансис и не париться :-D зачем такие сложности? конвертер...перепроверки...а снип и сп накой? по ним, с помощью экселя, можно и сделать поверочный расчет

[aka raduntsev-nias]

да. главное не ошибиться в ansys/ а потом в ексель.
программа правильно посчитает, что в ней только могут шаловливые ручки натворить.

[таи]

Позвольте полюбопытствовать, а как решается вопрос с сочетанием нагрузок в Ансисе?
Я Ансис не знаю - меня интересует сам алгоритм действий.
С одной нагрузкой более менее все ясно.

Но предположим имеется три вида нагружений постоянное - коэфф. 1, временное длительное - коэфф.0.95, кратковременная - коэфф. 0.9.
Итого 1*(1) 2*(0,9) 3 (0,95).
Просветите, если не трудно.

[derik]

поэтому сейчас штудирую мат. часть
с помощью команды loadcase, точнее lcdef...эт так навскидку, до комбинаций пока не добрался...а проще сделать как советовали знающие люди: нагрузил, решил; нагрузил, решил .
т.к. в ансисе присутствуют макросы, это даже упрощает дело

[aka raduntsev-nias]

_____________

[таи]

Сообщение от 6pa4o проще сделать как советовали знающие люди: нагрузил, решил; нагрузил, решил

Можно и так.

Сообщение от 6pa4o в ансисе присутствуют макросы, это даже упрощает дело

Понятно, надо тогда макрос найти будет или самому, что труднее, написать со временем, Offtop: (речь идет с моей стороны о code-aster).

Сообщение от aka raduntsev-nias да. главное не ошибиться в ansys/ а потом в ексель. программа правильно посчитает, что в ней только могут шаловливые ручки натворить.

А в экселе как считаете? Крутящий момент учитываете? Диаграммы деформирования железобетона? Вопрос сугубо практический, интересен подобный опыт - ту же Лиру проверять по ее же усилиям...

[СергейД]

если (как на юге, западе и востоке= словом везде) задают конкретные коэффициенты, то просто
можете суммировать либо загружения (LS), либо результаты от них (lcase)
операций достаточно.
+в civilfem есть процедуры.

с нашими РСУ кроме нас никто не работает,
я повозился немного с макросом и забросил= логике моей оказалось непосильно.
В НИИмосстрое такое что-то сделали и даже сертифицировали макрос в Госстрое.
когда нам очень нужно было= в скаде запомнили множители и
просто задали их в ансис.

про проверку лиры:
каждой программе свое. куда интереснее просчитать нелинейно то,
что наармировали в лире.
Жаль (слегка), что на прочность у меня сейчас времени нет совершенно
(на всякий случай извиняюсь)= на аэродинамику заказов полно.

[derik]

проще действительно так:
нагрузка*коэффциент - 1 загружение
нагрузка*коэффциент - 2-е загружение
и т.д....
потом смотрим по истрории какое самое "опасное" нагружение и по нему ведем расчет дальше...
данные экспортируем в эксель, считаем и фсе

Сообщение от 6pa4o анные экспортируем в эксель, считаем и фсе

В таком случае и к Вам вопрос из поста № 246...

[derik]

формулы из снипа и сп+ результаты из ансиса. в чем вопрос?

Сообщение от 6pa4o формулы из снипа и сп+ результаты из ансиса. в чем вопрос?

Ну хотя бы в подходе - Вы плоский элемент делите на слои или пользуетесь обобщенными уравнениями предельного равновесия..? (Сп-52-103-2007, п.8.4). И каков подход к выбору точек для проверки..?
Ну и как к работающим в Ansys - возможно ли изополя усилий преобразовывать в нужные мне изополя...? Т.е. я домножаю на какие либо нужные мне цифры значения изополей и получаю графически эту картинку - так можно?

[derik]

а как Вы проверяете скажем армирование полученные в программе, незнаю в какой Вы считаете?
про изополя ничего не могу сказать.

Проверяем по формулам СП-52-103-2007 и иногда по формулам как для изгибаемых элементов - проверяем результаты Лиры, собираемся познакомиться со Статикой и Микрофе для этих целях (в прошлом я успешно считал в них металл). Просто раньше в экселе считали исключительно стержни (тем более что лировские результаты для сжатых элементов у нас сильно расходились с ручным счетом, балки более менее сходились), а вот плиты проверяли как изгибаемые элементы... Но позже вышел сп 2007-го года где четко всплыла необходимость учета крутящих моментов... Написать то напишем в экселе, но теперь уже не так все просто - раньше мы могли рассчитать предельный момент для различных типов армирования и настроить шкалу усилий в Лире и вперед, а теперь результат зависит от нескольких шкал (моменты в двух направлениях, кручение..) Он и раньше зависел =), свойства железобетона не поменялись с тех пор, но раньше рассматривать плиты как полосовой элемент допускалось, а сейчас вот думаем...
И тогда уж еще один вопрос - а Вы пользуетесь Ansys ED? А где можно достать официально его для диссертации? На сайте ихнем что-то глухо... И не знаете ли - если я скачаю образовательную версию "где-то" - я смогу ей пользоваться легально?

[aka raduntsev-nias]

Сообщение от Ал-й Ну и как к работающим в Ansys - возможно ли изополя усилий преобразовывать в нужные мне изополя...? Т.е. я домножаю на какие либо нужные мне цифры значения изополей и получаю графически эту картинку - так можно?

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

"А где можно достать официально его для диссертации?"

Бесплатной версии для ansys теперь нет. Можно найти организацию имеющую на него лицензию и поработать так.

[таи]

Сообщение от Ал-й И тогда уж еще один вопрос - а Вы пользуетесь Ansys ED? Offtop: Извините, что вмешиваюсь в разговор. А где можно достать официально его для диссертации?

Мне кажется, Вы не слышали про code-aster - бесплатный открытый пакет - аналог ансиса по возможностям или достаточно близко по возможностям.

Поплутав по европейским сайтам, сделал вывод, что многие фрилансеры и просто специалисты пересаживаются с ансиса на code-aster и не особенно жалеют, - по своим возможностям он практически может равнять с ним.

Разумеется, требуется довольно значительное время для освоения, не меньше, видимо, чем и в вопросе с освоениея ансиса.

Сообщение от таи Мне кажется, Вы не слышали про code-aster - бесплатный открытый пакет - аналог ансиса по возможностям или достаточно близко по возможностям.

Слышал, и про Calculix тоже, просто пока равняюсь на преимущества Ansys в количестве ресурсов, специалистов, курсов, а раз речь об образовательной версии, то и цена не важна... А вот если речь зайдет о практическом применении, то возможно придется знакомиться с другими пакетами - я думаю знание Ansys (даже начальное) только поможет в их освоении - ведь логика программ наверняка похожа...
з.ы. Самые "крутые" на мой взгляд спецы Европы, те, кто реализовывал удивительные конструкции Калатравы, считали судя по всему в CivilFem, раз официальный сайт CivilFem'а весь в фотках его объектов =)
з.ы. з.ы. aka raduntsev-nias, спасибо за вариант, кажется такая организация на примете есть...

[aka raduntsev-nias]

Сообщение от Ал-й з.ы. з.ы. aka raduntsev-nias, спасибо за вариант, кажется такая организация на примете есть...

Не за что) Это вариант предложенный Сергей Ивановичем.
Я честно повторил.
Но, впрочем, так и поступил.

Сообщение от aka raduntsev-nias Не за что) Это вариант предложенный Сергей Ивановичем. Я честно повторил. Но, впрочем, так и поступил.

Раз уж Вы разбираетесь, то позволю себе еще один вопросик - вроде от Вас слышал, что есть образовательная версия Abaqus... Это так? И если да, то где можно ее получить..? И как много проблем при освоении - неужели совсем нет литературы и пользователей-советчиков? Моя предполагаемая модель - балка ж/б, усиленная углепластиковой арматурой, причем как нагруженная, так и разгруженная (т.е. для первого случая еще и учет последовательности загружения нужен) - необходимо будет сравнивать численный и натурный эксперимент... Abaqus поможет, или все же изыскивать друзей с Ansys?

[aka raduntsev-nias]

Сообщение от Ал-й Раз уж Вы разбираетесь

немного пока

Сообщение от Ал-й есть образовательная версия Abaqus.

да, заполните анкету и Вам ее вышлют

Сообщение от Ал-й Моя предполагаемая модель - балка ж/б, усиленная углепластиковой арматурой, причем как нагруженная, так и разгруженная (т.е. для первого случая еще и учет последовательности загружения нужен) - необходимо будет сравнивать численный и натурный эксперимент

Давно не занимаюсь жб, может вопросы и ответы в этой теме помогут. Если еще не видели. От себя добавлю лишь пару скринов расчета жб балки сделанного в ansys и один из примеров abaqus
Выбирать Вам, что я могу сказать. Сам до сих пор нахожусь в процессе выбора. Конечно, имея многое переведенное по ansys и примеры выполненных задач по специальности в нем же, проще то же воспроизвести в abaqus/ Который имеет неплохую справку и набор примеров. Единственный минус все на анг. Методичек и литературы по abaqus на русском очень мало.
Из плюсов ansys - есть wb, который проще классик, те можно начать с него, но потом все равно вернетесь к классике, для более серъезной обработки задачи. wb без знания APDL - слишком куцо пока. и еще плюс - есть сивилфем с набором готовых макросов для проектирования гражданских зданий.
те проще все таки в ansys начать работать

[Kostya_PC]

http://www.simulia.com/academics/student.html

aka raduntsev-nias, Kostya_PC - спасибо, подумаю над выбором.

[Cfytrr]

aka raduntsev-nias,
Насколько я понял, в этой модели бетон работает упруго даже в растянутой зоне

Потом, автор предлагает натянуть арматуру на соответствующее усилие и вполне закономерно считает, что после этого бетон действительно "трещать" не должен.
А нет ли у Вас подобных отчетов для простого железобетона, где бетон вынужден в сжатой и растянутой зоне работать по-разному?

[aka raduntsev-nias]

Cfytrr,
Вас интересует процесс развития трещин в растянутых зонах бетона? Те сравнение аналитического метода расчета жб балки с расчетом в мкэ программах и сравнимость с результатами эксперимента? те увидеть виртуальное образование трещин, разрушение балки и тд.
Такого отчета у меня нет. Поискал в сети, нашел наверное похожий
Тоже бы с удовольствием посмотрел на технику моделирования таких процессов.
Нашел. Сильно не вникал. Может не совсем то, что Вам нужно. )
Искал целенаправленно abaqus/ То же можно и в ansys найти думаю.

[derik]

Код:

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

fini
/clear,nostart
/outp,output,txt,C:\ansys\output
/cwd,'C:\ansys\out'
/units,mks

/prep7

*set,a,4
*set,b,.4
*set,c,.2
*set,q,5e4
*set,p,2e5

et,1,185
mptemp,,,,,,,
mptemp,1,0
mpdata,ex,1,,2.04e5
mpdata,prxy,1,,.3
tb,biso,1,1,2,
tbtemp,0
tbdata,,.84e5,3044

blc5,,,a,b,c
csys,5
vgen,2,all,,,,90,,,0
csys,0
vgen,,2,,,-b/2,,c/2,,,1
vovlap,all

mshkey,0
mshape,2,3d
esize,b/4
vsweep,all
nummrg,all

asel,s,,,5,6
asel,a,,,11,12
da,all,all,
sfa,34,1,press,q
alls
fini

/sol

/gst,1
nldiag,nrre,on
antype,0
nlgeom,1
nsubst,250,5500,75
outres,erase
outres,all,all
time,1
solve
fini

/post1

set,first
plnsol,uy,sum,0,1

Добрый вечер!
на критику первую нелинейную модель выставляю

[aka raduntsev-nias]

Сообщение от 6pa4o первую нелинейную модель выставляю

круто! у меня в вб так пока и не получилось с нелином

[st2008]

Сообщение от 6pa4o на критику первую нелинейную модель выставляю

Выдает такое сообщение...

One or more elements have become highly distorted. Excessive distortion of elements is usually a symptom indicating the need for corrective action elsewhere. Try incrementing the load more slowly (increase the number of substeps or decrease the time step size). You may need to improve your mesh to obtain elements with better aspect ratios. Also consider the behavior of materials, contact pairs, and/or constraint equations. If this message appears in the first iteration of first substep, be sure to perform element shape checking.

[derik]

вообщем нужно увеличить число подшагов или понизить временной шаг. с этой командой у меня пока туго не могу найти как подбирать эти самые шаги, их количество и мин. и макс. число шагов

[aka raduntsev-nias]

Сообщение от 6pa4o найти как подбирать эти самые шаги, их количество и мин. и макс. число шагов

или может с материалами еще что то не то

[derik]

материал брал из верификационного примера...то кол-во шагов которое там же указано, эту ошибку выдавало раз 5.

[soppr]

Маленькая контактная задачка для понимания контакта в ansys
Два кубика, один скользит по другому. Верхний кубик имеет модуль упругости 2000 и коэф. Пуассона 0.3. Нижний модуль 2e12 и Пуассона 0.3.
Верхний кубик нагружен на верхней площади распр. нагрузкой 2 и на боковой грани 0.5. Нижний кубик заделан со всех сторон (кроме контактной). Задано трение 0.3, отношение статического к динамическому трению - 1.2. Трение покоя (cohesion) - 0.4.

Код:

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

/COM,  Structural
!*  
/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,1,0,1,0,0.5,
BLOCK,0,1,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,
FLST,5,2,6,ORDE,2
FITEM,5,1
FITEM,5,-2  
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  
!*  

/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.3
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,,,1.0E20,0.0,1.0, 
RMORE,0.4,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

/SOL
FLST,2,5,5,ORDE,2
FITEM,2,8
FITEM,2,-12 
!*  
/GO 
DA,P51X,ALL,
FLST,2,1,5,ORDE,1
FITEM,2,1
/GO 
!*  
SFA,P51X,1,PRES,2
FLST,2,1,5,ORDE,1
FITEM,2,5
/GO 
!*  
SFA,P51X,1,PRES,0.5

Solve
/post1

Лог файл конечно кривоват)). Но работает)
Комментарии господа)

[aka raduntsev-nias]

Сообщение от soppr Комментарии господа)

решилась. но трение 0,3 - нессиметричный решатель нужно. или с меньшим решить

[soppr]

Сообщение от aka raduntsev-nias решилась. но трение 0,3 - нессиметричный решатель нужно. или с меньшим решить

От величины трения зависит решатель? Не мог бы подсказать, в каком разделе найти эту информацию - в решателях или контактной?

Возникли еще ряд вопрсов:
1. При выборе “жесткой” и “поддатливой” поверхностей руководствовался принципами, написанными в ANSYS help, те контактной поверхностью выбрал наименее жесткую (верхний куб). При этом расчет не производился. Как только поменял наименование поверхностей местами (нижний куб с большим модулем упругости стал контактным), расчет прошел. С чем связано?
2. При вводе коэффициента трения он присваивается какому то из материалов. Влияет ли на результаты задание трения разным материалам? (как я полагаю - нет)
3. Контактная жесткость. В helpANSYSсказано, что этот параметр вычисляется автоматически (как я полагаю в зависимости от модулей упругости, геометрии и других). Тем не менее есть рекомендация по начальному заданию контактной жесткости равной 1 для объемных задач. То есть его нужно подбирать вручную? Что это за параметр и как он влияет на результаты?

[aka raduntsev-nias]

сколько вопросов. да я то откуда все это знаю.
могу только посоветовать.

Сообщение от soppr 1. При выборе “жесткой” и “поддатливой” поверхностей

решение контактных задач, CADFEM

Сообщение от soppr При вводе коэффициента трения он присваивается какому то из материалов. Влияет ли на результаты задание трения разным материалам

на САПРе это обсуждали. так понял, что задавать ДВА разных коэффициента. Но не уверен.

Сообщение от soppr Контактная жесткость. В helpANSYSсказано, что этот параметр вычисляется автоматически

да, если установить опцию. имеет и другие значения. для трения можно 0,1 поставить. тут есть про это в теме.

[soppr]

1. CADFEM читал и сделал как там написано и именно так она не решалась) Поэтому попоробовал изменить CONTACT и TARGET местами.
2. В контактном мастере задается friction только для одного ID материала. Хотя если контактный мастер не использовать то можно задать каждому материалу. Поэтому и вопрос возник.

[aka raduntsev-nias]

NROP,UNSYM нужно добавить
и субстепов поболее. в контакте сетку побейте сильней.(на контакте)
еще у тебя временного интервала нет

[soppr]

raduntsev, большое спасибо за советы.

[aka raduntsev-nias]

soppr, если помогло, хорошо. не за что.
Сравнение

[derik]

погрешности бывают и у машин это абакус, Макс?

[aka raduntsev-nias]

Сообщение от 6pa4o это абакус, Макс?

да, Ярослав. если есть интерес, вот простой пример из хелпа abaqus
твоя задачка очень помогла, выручил. всю страну научил apdl у)

Код:

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

"""
beamExample.py

Reproduce the cantilever beam example from the
Appendix of the Getting Started with
Abaqus: Interactive Edition Manual.
"""

from abaqus import *
import testUtils
testUtils.setBackwardCompatibility()
from abaqusConstants import *

# Create a model.

myModel = mdb.Model(name='Beam')

# Create a new viewport in which to display the model
# and the results of the analysis.

myViewport = session.Viewport(name='Cantilever Beam Example',
    origin=(20, 20), width=150, height=120)
    
#-----------------------------------------------------

import part

# Create a sketch for the base feature.

mySketch = myModel.ConstrainedSketch(name='beamProfile',
    sheetSize=250.)

# Create the rectangle.

mySketch.rectangle(point1=(-100,10), point2=(100,-10))

# Create a three-dimensional, deformable part.

myBeam = myModel.Part(name='Beam', dimensionality=THREE_D,
    type=DEFORMABLE_BODY)

# Create the part's base feature by extruding the sketch 
# through a distance of 25.0.

myBeam.BaseSolidExtrude(sketch=mySketch, depth=25.0)

#-----------------------------------------------------

import material

# Create a material.

mySteel = myModel.Material(name='Steel')

# Create the elastic properties: youngsModulus is 209.E3
# and poissonsRatio is 0.3

elasticProperties = (209.E3, 0.3)
mySteel.Elastic(table=(elasticProperties, ) )

#-------------------------------------------------------

import section

# Create the solid section.

mySection = myModel.HomogeneousSolidSection(name='beamSection',
    material='Steel', thickness=1.0)

# Assign the section to the region. The region refers 
# to the single cell in this model.

region = (myBeam.cells,)
myBeam.SectionAssignment(region=region,
    sectionName='beamSection')

#-------------------------------------------------------

import assembly

# Create a part instance.

myAssembly = myModel.rootAssembly
myInstance = myAssembly.Instance(name='beamInstance',
    part=myBeam, dependent=OFF)

#-------------------------------------------------------

import step

# Create a step. The time period of the static step is 1.0, 
# and the initial incrementation is 0.1; the step is created
# after the initial step. 

myModel.StaticStep(name='beamLoad', previous='Initial',
    timePeriod=1.0, initialInc=0.1,
    description='Load the top of the beam.')

#-------------------------------------------------------

import load

# Find the end face using coordinates.

endFaceCenter = (-100,0,12.5)
endFace = myInstance.faces.findAt((endFaceCenter,) )

# Create a boundary condition that encastres one end
# of the beam.

endRegion = (endFace,)
myModel.EncastreBC(name='Fixed',createStepName='beamLoad',
    region=endRegion)

# Find the top face using coordinates.

topFaceCenter = (0,10,12.5)
topFace = myInstance.faces.findAt((topFaceCenter,) )

# Create a pressure load on the top face of the beam.

topSurface = ((topFace, SIDE1), )
myModel.Pressure(name='Pressure', createStepName='beamLoad',
    region=topSurface, magnitude=0.5)

#-------------------------------------------------------

import mesh

# Assign an element type to the part instance.

region = (myInstance.cells,)
elemType = mesh.ElemType(elemCode=C3D8I, elemLibrary=STANDARD)
myAssembly.setElementType(regions=region, elemTypes=(elemType,))

# Seed the part instance.

myAssembly.seedPartInstance(regions=(myInstance,), size=10.0)

# Mesh the part instance.

myAssembly.generateMesh(regions=(myInstance,))

# Display the meshed beam.

myViewport.assemblyDisplay.setValues(mesh=ON)
myViewport.assemblyDisplay.meshOptions.setValues(meshTechnique=ON)
myViewport.setValues(displayedObject=myAssembly)

#-------------------------------------------------------

import job

# Create an analysis job for the model and submit it.

jobName = 'beam_tutorial'
myJob = mdb.Job(name=jobName, model='Beam',
    description='Cantilever beam tutorial')

# Wait for the job to complete.

myJob.submit()
myJob.waitForCompletion()

#-------------------------------------------------------

import visualization

# Open the output database and display a
# default contour plot.

myOdb = visualization.openOdb(path=jobName + '.odb')
myViewport.setValues(displayedObject=myOdb)
myViewport.odbDisplay.display.setValues(plotState=CONTOURS_ON_DEF)

myViewport.odbDisplay.commonOptions.setValues(renderStyle=FILLED)

[derik]

Код:

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

fini
/clear,nostart
/outp,output,txt,C:\ansys\output
/cwd,'C:\ansys\out'
/units,mks

/prep7

*set,w1,.18
*set,w2,.18
*set,w3,.443
*set,t1,.011
*set,t2,.011
*set,t3,.0078
*set,a,3
*set,p,5e3

et,1,188
keyopt,1,8,1
keyopt,1,9,1
mptemp,,,,,,,
mptemp,1,0
mpdata,ex,1,,2.15e8
mpdata,prxy,1,,.3
tb,biso,1,1,2,
tbtemp,0
tbdata,,1.95e8,2540

sectype,1,beam,i,,0
secoffset,cent
secdata,w1,w2,w3,t1,t2,t3,
k,,
k,,a,
lstr,1,2

lesize,all,a/40,
lmesh,all

sfbeam,all,1,pres,p,p,

ksel,s,loc,x,0
dk,all,,,,,all,
alls
fini
/eshape,1

/sol

/gst,1
nldiag,nrre,on
antype,0
nlgeom,1
nsubst,20,100,10
outres,erase
outres,all,all
time,.5
solve
fini

/post1

set,first
etable,fx,smisc,1
etable,fx,smisc,14
etable,my,smisc,3
etable,mz,smisc,24
etable,tq,smisc,4
etable,sfz,smisc,5
plls,my,my,1,0

вот еще один пример - балочка обыкновенная!

[aka raduntsev-nias]

Вот еще интересная
ее в ansys. я сейчас это делаю. на задании законов нелин пока тормознулось
Вот, тут человек поделился файлом. На форуме нашел где-то
сетку и геометрию добавил.
координаты точек на кривых деформаций пока не нашел. это и есть моя проблема
вообще нелин(несчитая контактных), да еще в железобетоне ) делаю в первый раз.

[derik]

подобное и я моделировал, но без учета нелина даже поперечку накладывал, ну это так - от нечего делать а потом, по совету Сергея Ивановича, пересел на металлы, для начала...сижу читаю теорию пластичности и хелп ансиса с верификационными примерами...до бетона еще огого как далеко

[aka raduntsev-nias]

_удалил_чтоб Брачо недулся.
Позже свой пример может подготовлю.

[derik]

А где фундаментную плиту потерял?
грунт нужно через моедль друкера-прагера задавать, а бетон через et,65 и concrete...

[aka raduntsev-nias]

Сообщение от 6pa4o А где фундаментную плиту потерял? грунт нужно через моедль друкера-прагера задавать, а бетон через et,65 и concrete...

да знаю я. это же просто, пример . ответы разные в вб и классике.

С плитой на 3д грунте вроде ясно, не потому что самый умный, а только смысл. но как быть в случае столбчатых фундаментов? Как с ними каркас проверять? Кто нибудь пробовал проверить каркас на столбчатых ф-ах в ansys ? Наверное так же, нелин объемниками

Ув. специалисты!
Объясните - есть ли принципиальные отличия между различными по сроку давности версиями Ansys (9-13 или более ранние) с точки зрения обучения им. Меняются ли команды? Насколько эффективно заниматься с литературой по одной версии, а программой - по другой..? Пока примерно так =)

[derik]

разница есть - добавились новые КЭ, скорость расчетов уменьшилась...

[aka raduntsev-nias]

Сообщение от Ал-й Меняются ли команды

да, в новой немного изменились. добавились элементы даже.

Сообщение от 6pa4o разница есть - добавились новые КЭ, скорость расчетов уменьшилась...

А именно команды? Вот, скажем, буду я учить команды APDL по литературе для версии 9, а практически воспользуюсь 13-й версией... Я столкнусь с тем, что команды (способ их записи) поменяются..?

[СергейД]

все обиходные команды ( даже начиная с версии 5.0) еще работоспособны.
добавляются только новые.
с 9.0 проблем точно не будет.

СергейД, aka raduntsev-nias, 6pa4o, спасибо!

[Dakar]

Сообщение от 6pa4o скорость расчетов уменьшилась...

Извините за подколку, но может наоборот?

Хочуууу ANSYS изучаааааать!!!

[derik]

Offtop: точно

[derik]

Код:

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

fini
/clear,nostart
/outp,output,txt,C:\ansys\output
/cwd,'C:\ansys\out'
/units,mks

/prep7

*set,a,5
*set,w1,.18
*set,w2,.18
*set,w3,.443
*set,t1,.011
*set,t2,.011
*set,t3,.0078
*set,q,5e4
*set,E,206e6 !модуль упругости, кПа
*set,s02,345000 !условный предел текучести, кПа
*set,fs,490000 !предел прочности, кПа
*set,psi,63 !относительное сужение при разрыве, %
*set,nuxy,.3 !коэффициент Пуассона

! показатель упрочнения
m=0.75*log(fs*(1+1.4*(psi/100))/s02)/log(log(100/(100-psi))/((s02/E)+0.002))
! "истинный" предел текучести
sy=(s02/(E*0.002+s02)**m)**(1/(1-m))
! напряжение, соответствующее предельной деформации
sm=sy*(log(100/(100-psi))/(sy/E))**m

mp,ex,1,E
mp,nuxy,1,nuxy
tb,kinh,1,1,20 ! закон упрочнения материала
tbpt,,(sy/E),sy ! начальная точка графика sigma=f(e)
*do,i,1,18
tbpt,,(sy/E)*(1+((sm/sy)-1)*(i/19))**(1/m),sy+(sm-sy)*(i/19)
*enddo
tbpt,,log(100/(100-psi)),sm ! конечная точка графика sigma=f(e)
/go

et,1,185

k,,
k,,,,t1
k,,,w1/2-t3/2,t1
k,,,w1/2-t3/2,w3-t2
k,,,,w3-t2
k,,,,w3
k,,,w2,w3
k,,,w2,w3-t2
k,,,w1/2+t3/2,w3-t2
k,,,w1/2+t3/2,t1
k,,,w1,t1
k,,,w1,
a,1,2,3,4,5,6,7,8,9,10,11,12,1
vext,all,,,a,
asel,s,loc,z,w3
sfa,all,1,press,q
alls
asel,s,,,1,2
da,all,all,
alls

mshkey,0
mshape,2,3d
esize,t1*5
vsweep,all
alls
fini

/sol

/gst,1
nldiag,nrre,on
antype,0
nlgeom,1
nsubst,85,8500,15
outres,erase
outres,all,all
time,.1
solve
fini

/post1

set,first

вот...еще одна балка

[derik]

Код:

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

fini
/clear,nostart
/outp,output,txt,C:\ansys\output
/cwd,'C:\ansys\output'
/units,mks

/prep7

*set,a,4!м
*set,b,.2!м
*set,h,.4!м
*set,q,.3!МПа
*set,Eb,2.75e4!МПа
*set,Rb,11.5!МПа
*set,Rbt,.9!МПа
*set,ny,.2
*set,alfa,1.5e-6

et,1,65
keyopt,1,3,1
keyopt,1,6,3
keyopt,1,7,1
mptemp,,,,,,,
mptemp,1,0
mpdata,ex,1,,Eb
mpdata,nuxy,1,,ny
mpdata,alpx,1,,alfa
tb,conc,1,1,9,1
tbtemp,0
tbdata,1,.2,.6,Rbt,Rb,
tb,kinh,1,1,3,0 
tbtemp,0
tbpt,,0.6*Rb/eb,0.6*Rb
tbpt,,0.002,Rb 
tbpt,,0.0048,Rb

blc5,b/2,h/2,b,h,a
asel,s,loc,y,h
sfa,all,1,pres,q

mshkey,0
mshape,2,3d
esize,b/4
vsweep,all
asel,s,loc,z,a
asel,a,loc,z,0
da,all,all
alls
fini

/sol

cnvtol,f,,.8
/gst,1
nldiag,nrre,on
antype,0
nlgeom,1
nsubst,50,700,15
outres,erase
outres,all,all
time,.1
solve
fini

/post1

set,first
!set,last
!/device,vector,1
!plcrack,1,0

вот бетонная балка, прошу покритиковать

[derik]

___

[soppr]

Это понятно, вот поэтому я и выбрал балку 600*300*6000 для примера.

=))))))))))) Без палева, как говорится)

[derik]

Всмысле?

[СергейД]

обучение провожу... занят.
разве что как учебный пример рассмотреть попробую.

[derik]

пожалуйста только, если можно с исправлениями и пояснениями

СергейД, здравствуйте! У меня появились экспериментальные данные по испытанию кирпичных простенков на срез с различными вариантами усиления и коллеги, желающие узнать - можно ли все это смоделировать в Ansys (конечно можно, мы же видели...).
Могу я принести исходные данные на курсы - как пример порассматривать подобную задачу...?
Единственное, что в этот понедельник меня точно не отпустят... - буду только в пятницу...

[СергейД]

приносите. составьте в ансис модель. обсудим.посчитаем

[Sxela]

В посте http://forum.dwg.ru/showpost.php?p=658641&postcount=259 приведены картинки из примера расчета железобетонной балки в модуле Explicit Dynamics (ANSYS Worbench). Хотелось бы посмотреть сам пример. В tutorial его не нашел, есть только ряд аналогичных картинок без пояснений.
Интересует вопрос по подготовке геометрии элементов (тел) и установлении связи между ними, ну и разбиении на сетку.

В help написано (раздел Reinforcement Type), что при наличии линейных объектов внутри объема тела связь между ними устанавливается автоматически, но у меня ничего не получается (хочу реализовать данный пример с целью обучения Worcbench). Создал объем (solid), создал линейные элементы (line body), присвоил им требуемый тип сечения (cross section) и результат ноль. Пробовал объеденить элементы геометрии в один блок (part), но выдает ошибку о несовместимости типов (solid и line), что соответствует написанному в справке.

Помогите пожалуйста разобраться с данным примером.

P.S. В классическом ANSYS подобную балку с некоторыми упрощениями я реализовал, но нужен конкретно Worcbench

[Chardash]

Сообщение от Sxela Хотелось бы посмотреть сам пример

пожалуйста

насчет
P.S. "В классическом ANSYS подобную балку с некоторыми упрощениями я реализовал, но нужен конкретно Worcbench", тсссс, сейчас СергейД придет....

[Sxela]

и еще вопрос созрел: как назначить в WB тип КЭ Solid65 и возможно ли это вообще?
все что на шел по данной теме - это подключение mehanical APDL, но связь устанавливается односторонняя, только на выгрузку данных в классический ANSYS, изменения обратно в WB не переходят.

[Chardash]

Сообщение от Sxela возможно ли это вообще?

возможно. читайте матчасть.

подсказка: с помощью команд APDL
ps прав был Сергей Иванович.

[СергейД]

ну кто вас заставляет применять WB??? не для этих задач он разрабатывался.

[Sxela]

Сообщение от Chardash возможно. читайте матчасть. подсказка: с помощью команд APDL ps прав был Сергей Иванович.

не знаю кто такой Сергей Иванович, тем более не могу оценить степень его правоты по интересующему меня вопросу. Могу только догадываться по постам СергейД .

С APDL знаком достаточно хорошо и успешно использую с 2000 г., но в классическом ANSYS.
По изучению "матчасти" для WB есть существенные ограничения из-за незнания английского, по Ansys хоть рускоязычной литературы достаточно да еще с примерами кода APDL.

За сужение области поиска был бы благодарен, а то из "матчасти" в распоряжении только куцый help для v13.0, и тот на 80% под классику. Tutorial по WB отсутсвует.

За пример с балкой большое спасибо.

p.s. WB изучаю всего неделю, для реализации большинства текущих задач вполне работоспособная система, интуитивно понятна, особенно для пользователей, владеющих CAD системами. По большому счету остались вопросы только по програмированию и реализации двустаронней связи Ansys-WB, но наверное не для данного форума .

[Chardash]

Сообщение от Sxela Могу только догадываться по постам СергейД

Сообщение от Sxela С APDL знаком достаточно хорошо и успешно использую с 2000 г., но в классическом ANSYS.

тогда тем более, команды в вб и в классик практически не отличаются. Для их ввода есть клавиша, см скрин
матчасть, при желании, можно отыскать в инете . если что, пишите в личку, подмогу.

[derik]

___

[Chardash]

Брачо, в последней книге, что выложил в download, есть пример

[derik]

всю пересмотрел, честно, ничего не нашел! ткните носом пжалуйста

пункт 5 читал, вдоль и поперек

[Chardash]

6pa4o, я думал там (гл5) есть

[derik]

ты меня понял
как все таки создать это сечение и увидить картину моментов в эквивалентном сечении! неужели здесь нет людей которые могут подсказать?

[Chardash]

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

недавно был на презентации ansys для строительства.
вышел в твердой уверенности, что продолжу и дальше с ним работать.

ansys масштабная система. и научиться считать строительные задачи в нем, было бы очень полезно.

ИМХо использовать лучше всего как пару для обычной строительной расчетной(Лира, Скад, Робот)

[derik]

угу....лиру или старк, как сертифицированную систему в России и ансис

[derik]

___

[Chardash]

еще пример

Код:

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

!EX7.9
FINISH
/CLEAR
/PREP7
R1=250
R2=500
C=20
B=62.16
ET,1,LINK8
MP,EX,1,3030
MP,GXY,1,1096
R,1,317
! ET,1,BEAM4
! R,1,317,23770,2950,1,1
! RMORE,,9180
CSYS,1
K,1,R2,30 
KGEN,6,1,1,,,60
K,7,R1,,B 
KGEN,6,7,7,,,60
K,13,,,C + B 
CSYS,0 
*DO,I,1,5
L,I,I + 6
L,I,I + 7
*ENDDO
L,6,7
L,6,12
*DO,I,7,11 
L,I,I + 1
*ENDDO
L,12,7 
*DO,I,7,12
L,I,13
*ENDDO
LESIZE,ALL,,,1
!LESIZE,ALL,,,5
LMESH,ALL
KSEL,S,LOC,Z,0
DK,ALL,ALL
KSEL,ALL
FK,13,FZ, -1
/SOLU
ANTYPE,0
PSTRES,ON
SOLVE
FINISH
/SOLU
ANTYPE,1
BUCOPT,LANB,1
SOLVE
/POST1
SET,LAST
PLNSOL,U,Z,0,1.0

и свая.

Код:

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

!EX8.31
FINISH$/CLEAR$/PREP7
apile=1.0$hpil1=2.0$hpil2=8.0$asoil=7.0$hsoil=24
et,1,solid45$et,2,170$et,3,173
keyopt,3,9,0$keyopt,3,12,2
mp,ex,1,2.6e10$mp,prxy,1,0.167$mp,dens,1,2500
mp,ex,2,2.6e8$mp,prxy,2,0.42$mp,dens,2,1900
tb,dp,2$tbdata,,19000,31,29
mp,mu,3,0.2$r,1$r,2$r,3
blc4,,,asoil,asoil,-hsoil$blc4,,,1,1,-hpil2$vsbv,1,2
wpoff,apile$wprota,,,90$vsbw,all
wpoff,apile$wprota,,,90$vsbw,all$wpcsys,-1
wpoff,,,-hpil2$vsbw,all$numcmp,all
vatt,2,2,1$esize,0.5$mshkey,1$vmesh,all
vsel,none$blc4,,,apile,apile,hpil1 + hpil2
vatt,1,1,1$vmesh,all
asel,s,,,37,39,2$asel,a,,,34$nsla,s,1
real,3$type,2$mat,3$esurf,all$allsel,all
asel,s,,,6,10,2$nsla,s,1$real,3$type,3$esurf,all
asel,s,loc,x,0$asel,a,loc,y,0$da,all,symm
asel,s,loc,x,asoil$asel,a,loc,y,asoil$da,all,all
asel,s,loc,z,-hsoil$da,all,uz
asel,s,loc,z,hpil1$sfa,all,1,pres,50e4$allse,all
/solu$antype,static$nlgeom,on$autots,on
acel,,,9.8$time,1$nsubst,100$outres,all,all
solve$finish
/post1$set,last
esel,s,real,,1$plnsol,u,z$plnsol,s,z
esel,s,real,,2$plnsol,u,z$plnsol,s,eqv
esel,s,type,,3$etable,st1,nmisc,41$pletab,st1
etable,sz1,nmisc,45$pletab,cz1

Удобнее писать после $

[derik]

Здрасте!
помогите через патч получить момент по длине пути.
Заранее спасибо!

Код:

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

finish
/clear
b0=.2
h0=.300
l0=3.000
ec=3e6
p0=100
/prep7
et,1,solid95
mp,ex,1,ec
mp,prxy,1,.2
mp,dens,2.5
blc4,,,b0,h0,l0
/view,1,1,1,1
/ang,1
vplot
asel,s,loc,z,0
asel,a,loc,z,l0
da,all,all
asel,s,loc,y,h0
sfa,all,1,pres,p0
alls
esize,.1
vsweep,all
finish
/solu
solve
finish
/post1
padel,all
/pbc,path,1
path,mx,2,60,40
ppath,1,,b0/2,0,0
ppath,2,,b0/2,0,l0
pdef,sigmx,s,x,avg
pcalc,intg,mx,sigmx,s

[Chardash]

под другим углом нужно сечение.
Точнее, на САРПе есть про это, с отрисовкой графиков.

[derik]

немного топик аппа

смоделировал балку из диссертации Зиганшина Х.А.. Балка под номером Б-II-2а, все характеристики материалов, закрепления взяты из данной работы. В аттаче картинка трещинообразования по опытным данным и смоделированная в ансис.

[MrWhite]

Картинки отличные, а критическая сила насколько отличается от опыта?

[tqthanh]

Здравствуйте! 6ра4о
Интересно, а какие характеристики материалов, закрепления из диссертации Зиганшина, и какой тип элементов Вы использовали при моделировании в ансис, если не секрет? Меня интересует проблема моделированния железжбетонных конструкций в анснс, однако у меня нет доступа к российской базе диссертаций.
С интересом следил за этой темой. Кстати, я иностранец, но это неважно.
Заранее спасибо.

[derik]

материалы из диссера кроме модуля упругости бетона и арматуры, элемент в ансисе 65 для бетона и 264 для арматуры;concr модель материала бетона с кинематическим упрочнением, арматура - просто кинематическое упрочнение . данный диссер есть на сайте жбк.рф в разделе библиотека, есть не сможете скачать напишите, "намылю".

[tqthanh]

6ра4о!
Большое спасибо за наводку. Стараюсь все материалы данные Вами изучать.

[derik]

Пожалуйста, но Вы не на меня ориентируйтесь, я только учусь! MrWhite и СергейД вот чьи сообщения бесценны!

[tqthanh]

М... да!
Ещё раз спасибо. Не волнуйтесь! Я тайно слежу за их дисскурсей.

[ruslannur]

Подскажите, что я делаю не так:

/UNITS, SI !единицы измерения СИ
s1=0.010 !толщина позиции 1
s2=0.010 !толщина позиции 2
Dotv=0.009 !диаметр отверстий
H=0.098 !Длина
B=0.030 !ширина
D1=0.170 !высота поз1
D2=0.125 !высота поз2

a=0.070 !расстояние между центрами отверстий
c1=0.015 !отступ до отверстия по вертикали
c2=0.019 !отступ до отверстия по горизонтали

F1=D1-D2+0.010 !размеры выреза
F2=0.038
R1=0.005

CH=0.010 !РАЗМЕР ФАСКИ

PI=ACOS(-1)
G=400*0.5*1.1 !Вес блока, Н
MOM=G*0.12 !Пара сил (момент)
PRESM=(mom/2*a)/(PI*s1*(Dotv/2)) !давление на стенку отверстия от пары сил
PRES=(G/3)/(PI*s1*(Dotv/2)) !давление на стенку отверстия от веса блока
hr=0.030 !положение верха пластины отноcительно
sp=0.012 !толщина пластины

W=7000 !сила ветра на одну опору, Н
PRESW=(W/3)/(PI*s1*(Dotv/2)) !давление на стенку отверстия от ветра


/prep7

eqslv,pcg,1e-5

MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,EX,1,,2e11
MPDATA,PRXY,1,,0.28
MPDATA,DENS,1,,7850


ET,1,solid186

BLC4,0,-S1/2,H-S2,S1,-D1
BLC4,H-S2,-B/2,S2,B,-D2

k,101,H-S2-F2, -s1/2, -D1
k,102,H-S2-F2, -s1/2, -D1+F1
k,103,H-s2, -s1/2, -D1+F1
k,104,H-s2, -S1/2, -D1

L,101,102
L,102,103
L,103,104
L,104,101

LFILLT,25,26,R1

AL,25,26,27,28,29

VEXT,13,,,0,S1,0,,,,

VSBV,1,3

k,201,c2,-S1/2,-c1 ! центры отверстий
k,202,c2,-S1/2,-(c1+a) ! центры отверстий
k,203,c2,-S1/2,-(c1+2*a) ! центры отверстий

k,204,c2,S1/2,-c1 ! центры отверстий
k,205,c2,S1/2,-(c1+a) ! центры отверстий
k,206,c2,S1/2,-(c1+2*a) ! центры отверстий

circle,201,Dotv/2,204
circle,202,Dotv/2,205
circle,203,Dotv/2,206

AL,2,3,4,10
AL,11,27,28,33
AL,34,39,44,45
VEXT,1,,,0,S1,0,,,,
VEXT,3,4,,0,S1,0,,,,

VSBV,4,1
VSBV,6,3
VSBV,1,5

k,,H-S2, -B/2+CH, -D2
k,,H-s2, -B/2, -(D2-CH)
k,,H-s2, -B/2, -D2
A,45,46,47
VEXT,1,,,s2,0,0,,,,
VSYMM,Y,1,,,
VSBV,2,1
VSBV,5,4

k,251,H, B/2, -hr !пятно контакта с пластиной
k,252,H, -B/2, -hr
k,253,H, -B/2, -(hr+sp)
k,254,H, B/2, -(hr+sp)

A,251,252,253,254

VEXT,10,,,-S2,0,0,,,,

VOVLAP,1,2

ESIZE,S1/3
SMRT,1
MSHAPE,1,3D
MSHKEY,0

CM,_Y,VOLU
VSEL,S,VOLU,,1,6
CM,_Y1,VOLU
CHKMSH,'VOLU'
CMSEL,S,_Y
VMESH,_Y1
CMDELE,_Y
CMDELE,_Y1
CMDELE,_Y2

ASEL,S,AREA,,13,14
ASEL,A,AREA,,18,19
ASEL,A,AREA,,25,28
ASEL,A,AREA,,30,33
AREFINE,ALL,,,1

/SOL

DA,10,UX,0
DA,10,UY,0

DA,4,UX,0
DA,4,UZ,0

asel,s,area,,13 !нагрузка в отверстиях от веса блока
asel,a,area,,14
asel,a,area,,25
asel,a,area,,26
asel,a,area,,30
asel,a,area,,31
sfa,all,,pres,PRES

allsel,all

asel,s,area,,14 !нагрузка от момента в верхнем отверстии
asel,a,area,,18
sfa,all,,pres,PRESM

allsel,all

asel,s,area,,30 !нагрузка от момента в нижнем отверстии
asel,a,area,,33
sfa,all,,pres,PRESM

allsel,all

asel,s,area,,13 !ветровая нагрузка
asel,a,area,,19
asel,a,area,,25
asel,a,area,,28
asel,a,area,,30
asel,a,area,,33
sfa,all,,pres,PRESW

allsel,all

solve

[MrWhite]

Ну 186 без опций оставить уже не правильно.. А в чем проблема можете описать, запустить сейчас не могу?

[derik]

в месте соединения ребра и пластины у Вас несогласованная сетка - Вам необходимо пересечь эти два тела. Сообщение говорит что получаются перемещения узла такого то больше чем допустимо. ребро удлините до наружной плоскости пластины и введите vptn,all.

[ruslannur]

Спасибо, с vptn все получилось.
Посмотрите еще, правильно ли задан контакт?

/UNITS, SI !переход к системе единиц измерения СИ

s1=0.010 !толщина средний стенки
s2=0.010 !толщина стенки основания
Dotv=0.009 !диаметр отверстий
H=0.098 !Длина
B=0.030 !ширина
D1=0.170 !высота поз1
D2=0.100 !высота поз2

a=0.070 !расстояние между центрами отверстий
c1=0.015 !отступ до отверстия по вертикали
c2=0.019 !отступ до отверстия по горизонтали

F1=D1-D2+0.010 !размеры выреза
F2=0.038
R1=0.005

CH=0.010 !РАЗМЕР ФАСКИ

PI=ACOS(-1)
G=4000*0.5*1.1 !Вес блока, Н
MOM=G*0.12 !Пара сил
PRESM=(mom/2*a)/(PI*s1*(Dotv/2))!давление на стенку отверстия от пары сил
PRES=(G/3)/(PI*s1*(Dotv/2)) !давление на стенку отверстия от веса блока
hr=0.030 !положение верха пластины отноcительно
sp=0.012 !толщина пластины

W=7000 !сила ветра на одну опору
PRESW=(W/3)/(PI*s1*(Dotv/2)) !давление на стенку отверстия от ветра

PB=0.036 !ширина пластины
PH=0.110 !длина пластины
PC=0.060 !расстояние от края пластины до центра детали
PS=0.012 !толщина пластины

/prep7

eqslv,pcg,1e-5

MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,EX,1,,2e11
MPDATA,PRXY,1,,0.28
MPDATA,DENS,1,,7850

ET,1,solid186

BLC4,0,-S1/2,H-S2,S1,-D1
BLC4,H-S2,-B/2,S2,B,-D2

k,101,H-S2-F2, -s1/2, -D1
k,102,H-S2-F2, -s1/2, -D1+F1
k,103,H-s2, -s1/2, -D1+F1
k,104,H-s2, -S1/2, -D1

L,101,102
L,102,103
L,103,104
L,104,101

LFILLT,25,26,R1

AL,25,26,27,28,29

VEXT,13,,,0,S1,0,,,,

VSBV,1,3

k,201,c2,-S1/2,-c1 ! центры отверстий
k,202,c2,-S1/2,-(c1+a) ! центры отверстий
k,203,c2,-S1/2,-(c1+2*a) ! центры отверстий

k,204,c2,S1/2,-c1 ! центры отверстий
k,205,c2,S1/2,-(c1+a) ! центры отверстий
k,206,c2,S1/2,-(c1+2*a) ! центры отверстий

circle,201,Dotv/2,204
circle,202,Dotv/2,205
circle,203,Dotv/2,206

AL,2,3,4,10
AL,11,27,28,33
AL,34,39,44,45
VEXT,1,,,0,S1,0,,,,
VEXT,3,4,,0,S1,0,,,,

VSBV,4,1
VSBV,6,3
VSBV,1,5

k,,H-S2, -B/2+CH, -D2
k,,H-s2, -B/2, -(D2-CH)
k,,H-s2, -B/2, -D2
A,45,46,47
VEXT,1,,,s2,0,0,,,,
VSYMM,Y,1,,,
VSBV,2,1
VSBV,5,4

k,251,H, B/2, -hr !пятно контакта с пластиной
k,252,H, -B/2, -hr
k,253,H, -B/2, -(hr+sp)
k,254,H, B/2, -(hr+sp)

A,251,252,253,254

VEXT,10,,,-S2,0,0,,,,

VOVLAP,1,2

nummrg,ALL

!mshkey,1
!mshape,0,3d
!esize,S1/9
!vsweep,all

VPTN,ALL

BLOCK,H-S2/2-PB/2,H-S2/2+PB/2,PC,PC-PH,-D2,-D2-PS, !пластина

ESIZE,S1/3
SMRT,1
MSHAPE,1,3D
MSHKEY,0

CM,_Y,VOLU
VSEL,S,VOLU,,1,8
CM,_Y1,VOLU
CHKMSH,'VOLU'
CMSEL,S,_Y
VMESH,_Y1
CMDELE,_Y
CMDELE,_Y1
CMDELE,_Y2

ASEL,S,AREA,,13,14
ASEL,A,AREA,,18,19
ASEL,A,AREA,,25,28
ASEL,A,AREA,,30,33
AREFINE,ALL,,,1

ASEL,S,AREA,,9
ASEL,A,AREA,,45
ASEL,A,AREA,,49
AREFINE,ALL,,,1

ASEL,S,AREA,,1
ASEL,A,AREA,,7
AREFINE,ALL,,,1

ALLSEL, ALL

!генерация контакта через contact wizard

CM,_TARGET,AREA
!*
CM,_CONTACT,AREA
!*
/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.2
MAT,1
MP,EMIS,1,7.88860905221e-031
R,3
REAL,3
ET,2,170
ET,3,174
R,3,,,1.0,0.1,0,
RMORE,,,1.0E20,0.0,1.0,
RMORE,0.0,0,1.0,,1.0,0.5
RMORE,0,1.0,1.0,0.0,,1.0
KEYOPT,3,4,0
KEYOPT,3,5,1
KEYOPT,3,7,0
KEYOPT,3,8,0
KEYOPT,3,9,0
KEYOPT,3,10,2
KEYOPT,3,11,0
KEYOPT,3,12,0
KEYOPT,3,2,0
KEYOPT,2,5,0
! Generate the target surface
ASEL,S,,,3
ASEL,A,,,8
ASEL,A,,,11
ASEL,A,,,12
ASEL,A,,,21
ASEL,A,,,22
CM,_TARGET,AREA
TYPE,2
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,,,4
ASEL,A,,,6
ASEL,A,,,35
ASEL,A,,,40
ASEL,A,,,41
ASEL,A,,,43
ASEL,A,,,44
ASEL,A,,,49
ASEL,A,,,50
CM,_CONTACT,AREA
TYPE,3
NSLA,S,1
ESLN,S,0
NSLE,A,CT2 ! CZMESH patch (fsk qt-40109 8/2008)
ESURF
ALLSEL
ESEL,ALL
ESEL,S,TYPE,,2
ESEL,A,TYPE,,3
ESEL,R,REAL,,3
/PSYMB,ESYS,1
/PNUM,TYPE,1
/NUM,1
EPLOT
ESEL,ALL
ESEL,S,TYPE,,2
ESEL,A,TYPE,,3
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

ALLSEL, ALL

/SOL

!DA,10,ALL,0
DA,1,UX,0
DA,7,UX,0
!DA,34,ALL,0
!DA,43,ALL,0
!DA,4,ALL,0
!DA,4,UX,0
!DA,4,UZ,0

asel,s,area,,13 !нагрузка в отверстиях от веса блока
asel,a,area,,14
asel,a,area,,25
asel,a,area,,26
asel,a,area,,30
asel,a,area,,31
sfa,all,,pres,PRES

allsel,all

asel,s,area,,14 !нагрузка от момента в верхнем отверстии
asel,a,area,,18
sfa,all,,pres,PRESM

allsel,all

asel,s,area,,30 !нагрузка от момента в нижнем отверстии
asel,a,area,,33
sfa,all,,pres,PRESM

allsel,all

asel,s,area,,13 !ветровая нагрузка
asel,a,area,,19
asel,a,area,,25
asel,a,area,,28
asel,a,area,,30
asel,a,area,,33
sfa,all,,pres,PRESW

allsel,all


DA,12,ALL !Закрепление пластины
DA,11,SYMM


solve


/post1
PLESOL, S,INT, 0,1.0

[Stanum]

Доброго времени суток, только начинаю разбираться с solid65 и бетоном в ansys, и у меня такой вопрос возник, при расчете, какие необходимо вводить данные: кубиковая прочность или призменная? Насколько сильно влияет ввод таких данных как прочность на двухосное сжатие и прочность при гидростатическом давлении, и как следствие где можно найти эти данные о материале?

[СергейД]

советов универсальных дать трудно.
если собираетесь сравнивать со снип (чтобы показывать практикам) нужно брать из снип два предела, на растяжение и на сжатие.
остальное по умолчанию. (там есть значения- разберитесь в логике их назначения)
у меня с арбатом отличия получаются не более 5%. далее можно считать то, что арбату не по силам.

если сравниваете с экспериментом- надо подбирать характеристики и параметры, чтобы совпало с замерами. начиная с тех, что назначены по итогам эксперимента.

все дополнительные параметры лично я толком не варьировал.
Думаю, их влияние более важно для управления сходимостью в существенно трехосных задачах (вроде плотин).

надеюсь, через месяц-два я смогу добавить информации.

[frostyfrost]

Сергей Иванович, вопрос к Вам есть небольшой. Давили ли Вы кубики, цилиндры и призмы в ANSYS, чтобы на выходе иметь при одинаковых свойствах материалов идентичную кубиковую, цилиндрическую и призменную прочность? И если давили, то какие результаты?

[Stanum]

СергейД, спасибо за ответ, и если можно еще такой вопрос:
В ansys необходимо вводить истинные напряжения - деформации в качестве исходных данных, для бетонного образца во всех виденных мной источниках присутствует ниспадающая ветвь, как я понял это связано с тем, что происходит разрушения сечения образца (образование пирамид), а следовательно диаграммы в нормах представлены в инженерных напряжениях деформациях, переход от которых к истинным через логарифм не сделать (как сделать этот переход если честно, в литературе пока не нашел), отсюда вопрос:

Следует ли задавать материалу ниспадающую ветвь, или лучше заканчивать прямым участком?
Проверил оба варианта на задаче по раздавливанию кубика (с ниспадающей ветвью и без). С ниспадающей ветвью задача решилась лучше ( эквивалентные относительные пластических деформации + упругие относительные деформации к моменту обрыва решения(оборвал сам, так как пошло много кратное деление шага безрезультатное) соответствовали исходной диаграмме.

Без ниспадающей ветви задача разошлась по достижению предельных напряжений.

[СергейД]

задаю ниспадающую ветвь kinh

кубики-призмы не давил. делал подобное только для грунтов.
напишите в личку данные= я решу, как привык. едва только найду время.

[frostyfrost]

Сергей Иванович, предлагаю для вас и коллег следующее описание верификационного теста. Интересно насколько модель бетона для 65-х сможет его выполнить. Давятся образцы со следующими параметрами:
- цилиндр 150х300 мм,
- куб - 150х150х150 мм,
- призма - 150х150х400 мм.
Параметры материалов одинаковые для всех вариантов. Варианта испытанию как я понимаю возможно два.
1-й вариант - крепление нижнего торца по XYZ, верхнего по XY, но с приложением перемещений эквивалентных предельной деформации сжатия. По суммарной опорной реакции по Z определение графика деформирования с определением прочности на сжатие (кубической, призменной и цилиндрической, которые связаны известными зависимостями).
2-й вариант - контактная задача с теми же параметрами.

[derik]

Код:

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

!/OUTPUT,NULL$/NOPR
/PREP7
/GRAPHICS,POWER
!ПАРАМЕТРЫ
*SET,A,.15 			!М
*SET,B,.15 			!М
*SET,H,.15 			!М
*SET,Q,21E6			!ПА
*SET,N,8			!ЭЛЕМЕНТОВ ПО ШИРИНЕ
*SET,ST,200			!МИН КОЛ-ВО ШАГОВ
*SET,STM,1E6		!МАКС КОЛ-ВО ШАГОВ
!ХАРАКТЕРИСТИКИ БЕТОНА В20
*SET,EB,3.75E10 !ПА
*SET,D_CONCRET,2500 !КГ/М3
*SET,RB,20.5E6		!ПА
*SET,RBT,2.2E6 		!ПА
*SET,N_CONCR,.2
*SET,ALFA_BR,1E-5
!НАЗНАЧЕНИЕ ЭЛЕМЕНТА
ET,1,65,,,2,,1,3$KEYOPT,1,7,1$KEYOPT,1,8,1
!НАЗНАЧЕНИЕ МАТЕРИАЛА
UIMP,1,EX,DENS,NUXY,EB,D_CONCRET,N_CONCR
!НЕЛИНЕЙНЫЙ БЕТОН
TB,CONC,1,1,9,1
TBDATA,,.3,.6,RBT,RB,,,
TBDATA,,,,.6,,,
TB,KINH,1,1,3
TBPT,,0.000328,RB*0.6
TBPT,,0.0034,RB
TBPT,,0.0048,RB
!ГЕОМЕТРИЯ
BLOCK,0,A/2,0,B/2,0,H
!СЕТКА
MSHKEY,2
MSHAPE,0,3D
ESIZE,A/2/N
VSWEEP,ALL
!ЗАКРЕПЛЕНИЯ И НАГРУЗКИ
ASEL,S,LOC,Z,0
DA,ALL,ALL
ASEL,S,LOC,X,0
ASEL,A,LOC,Y,0
DA,ALL,SYMM
ASEL,S,LOC,Z,H
SFA,ALL,1,PRES,Q
NC=NODE(A/2,B/2,H/2)
ALLS$ACEL,,,9.8$FINI
!НАСТРОЙКИ РЕШАТЕЛЯ
/SOL
	/GST,1,1
	SOLCONTROL,ON
	NLDIAG,NRRE,ON
	EQSLV,PSG,1E-6
	NROPT,UNSYM,,ON
	OUTRES,ALL,ALL
	LNSRCH,ON
	NLGEOM,ON
	PRED,ON
	AUTOTS,ON
	CNVTOL,F,,.05
	CNVTOL,U,,.09
	CUTCON,NOSHAPE,1
	CUTCON,PLSLIMIT,1
	NSUBST,ST,STM,ST,1
	TIME,30!Q/1000
	NEQIT,50
	SOLVE,,,,,NOCHECK
!РЕЗУЛЬТАТЫ
FINI$/POST1

[frostyfrost]

6pa4o, и как в итоге? Получилось получить призменную, кубическую и цилиндрическую прочность?

[derik]

я тока макрос сделал, нет времени на анализ. я бегло смотрел, графики не строил - но выглядит правдоподобно

[frostyfrost]

Там важно получить конкретную прочность. Поскольку для куба и призмы при одинаковых площадях опирания должна получиться разная прочность.

[СергейД]

немного уточнений=
EQSLV,PSG,1E-6! надо pcg
NROPT,UNSYM,,ON! не активируется, поскольку только со sparse
CNVTOL,F,,.05
CNVTOL,U,,.09
NEQIT,50!

[derik]

ставте свои данные и будет вам и призма и куб сейчас куб в макросе забит, точнее 1/4 куба.
Сергей Иванович, спасибо!!! изменения внесу

[СергейД]

немного уточнений=
EQSLV,PSG,1E-6! надо pcg
NROPT,UNSYM,,ON! не активируется, поскольку только со sparse
CNVTOL,F,,.05! начинаю с .50 потом уменьшаю. обычно уменьшение на результаты влияния не оказывает= просто решение дольше и может не сойтись.
CNVTOL,U,,.09! отключаю
NEQIT,50! 7 ставлю= много итераций нет резона=лучше шаг меньше делать
! решать нужно в контактной постановке.

[derik]

Сообщение от СергейД немного уточнений= EQSLV,PSG,1E-6! надо pcg NROPT,UNSYM,,ON! не активируется, поскольку только со sparse CNVTOL,F,,.05! начинаю с .50 потом уменьшаю. обычно уменьшение на результаты влияния не оказывает= просто решение дольше и может не сойтись. CNVTOL,U,,.09! отключаю NEQIT,50! 7 ставлю= много итераций нет резона=лучше шаг меньше делать

стояло 10, и шагов побольше, но бывает он шаг делит до бесконечности, и несколькими итерациями довел что более менее решение сходится/расходится. сходимость тоже варьирую, это уже "нормально" (вроде как), если ставить чуть больше и + + + + то решение проскакивает

[Stanum]

6pa4o, а здесь лучше грузить деформацией, или всё же нагрузкой?
equivalent plastic strain+ elastic strain/equivalent plastic stress = входной диаграмме напряжения деформации для точки центра куба - говорит ли это о приемлимости решения?

[frostyfrost]

Если лабораторные испытания моделируете, то перемещениями вы сможете в ряде случаев ниспадающую ветвь поймать, тогда как для нагрузки будет "прощелкивание". Хотя все зависит от задачи.

[Stanum]

Продолжаю разбираться с solid65, читая хелп у меня возникли вопросы по таким параметрам как:
real constant CSTIF(жесткость раздавленных элементов) имеет ли смысл увеличивать ее исходное значение 1е-6?
Параметр в solcontrol INCP - роль этого параметра непонятна, по нему практически нет описания кроме как того что он добавляет жесткость от нагрузки давлением, с чем и зачем его есть непонятно?
Также вычитал в интернете, что есть по solid65 не документированные возможности, где о них можно узнать поподробней?

[retal]

Обычно при сжатии призмы разрушение начинается с образованием продольных трещин, в конце вообще вырисовывается "огрызок" в виде 2-х пирамидок. Запускал код 6pa4o такой картинки не увидел.

[frostyfrost]

Сообщение от retal Обычно при сжатии призмы разрушение начинается с образованием продольных трещин, в конце вообще вырисовывается "огрызок" в виде 2-х пирамидок

Зависит от закрепления по концам. Возможны 3 варианта. Продольные трещины, одна наклонная и огрызок. Достаточно специфическая задача моделирования бетона и в частности моделирования достоверного процесса разрушения.

[СергейД]

там еще трение о плиты может повлиять.

[Stanum]

Решал ту же задачу (модель куб между двух плит(контактная задача с трением)) столкнулся с проблемой не сходится реакция для кубика. Решил проверить в бетоне ли проблема. В результате прогнал 2 задачи:
1) элементы куба solid 185(пробовал с разными опциями, полное интегрирование, улучшенные относительные деформации)
2) элементы куба solid 65
материал для них одинаковый: только tb,kinh без tb, conc/ В результате две абсолютно разные реакции к тому же задача с 65 расходится на самом последнем шаге.

вот сам текст задачи

Код:

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

 finish
/CLEAR,NOSTART  
/units, si
/filnam, kinh
/title, kinh
/OUTPUT, kinh_1, txt, , D:\ANSYS files
!бетон не обжатый
!параметры
fc=40e6!прочность не обжатого бетона максимальная (М)Па
fct=2.1e6! прочность бетона на растяжение (М)Па
ec1=0.00186! напряжения соответствующие макс напряжениям
ecu=0.003! предельные напряжения в необжатом бетоне
Eb=32e9! модуль упругости бетона (Г)Па
nub=0.2
an=15
*dim, eps_b, array, 20
*dim, sig_b, array, 20
eps_b(1)=0.3*fc/Eb
eps_b(an)=ec1
eps_b(20)=ecu

*do,i,2,an
eps_b(i)=eps_b(i-1)+(eps_b(an)-eps_b(1))/(an-1)
*enddo
*do,i,an+1,20
eps_b(i)=eps_b(i-1)+(eps_b(20)-eps_b(an))/(20-an)
*enddo

sig_b(1)=0.3*fc
*do,i,2,20
x=eps_b(i)/ec1
Ebc1=5e6*(fc**(1/2))
Esec1=fc/ec1
r1=Ebc1/(Ebc1-Esec1)
sig_b(i)=fc*x*r1/(r1-1+x**r1)
*enddo

!бетон без ниспадающей ветви
*dim, eps_bn, array, an+1
*dim, sig_bn, array, an+1
*do,i,1,an
eps_bn(i)=eps_b(i)
sig_bn(i)=sig_b(i)
*enddo
eps_bn(an+1)=eps_b(20)
sig_bn(an+1)=sig_b(an)*1.01

!сталь
nus=0.3! коэффициент пуассона стали
Es=210e9!модуль упругости стали, (Г)Па


/prep7
!тип элемента
et,1,185
et,2,65
et,3,185
keyopt,1,2,3
!keyopt,1,1,1
!keyopt,2,3,2
!keyopt,2,6,3
!keyopt,2,7,1
!keyopt,2,10,1
keyopt,3,2,0

!материала
mp,ex,1,Eb
mp,nuxy,1,nub
mp,mu,1,0.3
TB, conc, 1, 1, 9, 0, 
tbtemp,0
tbdata,1,0.2,0.6,fct,fc,
tb,kinh,1,1,20,0 
tbtemp,0
*do, i, 1, 20
tbpt, ,Eps_B(i), Sig_B(i)
*enddo

mp,ex,5,Es
mp,nuxy,5,nus
mp,mu,5,0.3						! коэффициент трения

mp,ex,6,Eb
mp,nuxy,6,nub
mp,mu,6,0.3
!TB, conc, 6, 1, 9, 0, 
!tbtemp,0
!tbdata,1,0.2,0.6,fct,fc,
tb,kinh,6,1,an+1,0 
tbtemp,0
*do, i, 1, an+1
tbpt, ,Eps_Bn(i), Sig_Bn(i)
*enddo


!построение модели
a=0.15
a1=0.151
hpl=0.05
apl=0.2
apl1=0.201
block,0,-a1/2,0,a,0,-a/2
vsel,u,,,all
block,-a1/2,0,-hpl,0,-a/2,0
block,-apl1/2,-a1/2,-hpl,0,-a/2,0
block,-a1/2,0,-hpl,0,-apl/2,-a/2
block,-apl1/2,-a1/2,-hpl,0,-apl/2,-a/2
vglue,all
vsel,u,,,all
block,-a1/2,0,a,a+hpl,-a/2,0,
block,-apl1/2,-a1/2,a,a+hpl,-a/2
block,-a1/2,0,a,a+hpl,-apl/2,-a/2
block,-apl1/2,-a1/2,a,a+hpl,-apl/2,-a/2
vglue,all
alls
!vglue,all

!граничные условия
asel,s,loc,y,-hpl,-hpl
Da,all,Uy,0
asel,s,loc,x,0,0
asel,a,loc,z,0,0
DA,all,Symm

alls

!сетка
elem=a/20
lsel,s,length,,hpl
LESIZE,ALL, , ,6, ,1, , ,1,
lsel,inve
LESIZE,ALL,elem, , , ,1, , ,1,  
alls
vsel,s,loc,y,0,a
VATT, 6, 1, 2, !(3)
vsel,s,loc,y,-hpl,0
VATT, 5, 1, 1, 
vsel,s,loc,y,a,a+hpl
VATT, 5, 1, 1, 
alls
MSHAPE,0,3d 
MSHKEY,1
VMESH,all
alls

!контакт
MAT,6							!задание номера материала 
R,101							!константы для третего типа кэ
REAL,101							!номер типа констант
R,102							!константы для третего типа кэ
REAL,102							!номер типа констант
ET,100,170						!создаёт тип кэ цель
ET,101,173						!создаёт тип кэ контактор

	FKN=10
	FTOLN=0.01
	!ICONT=0.2
	PINB	=0.1
	!PMAX=0
	!PMIN=0
	TAUMAX=1.00E+20
	!CNOF=0.001
	FKT=1
	!COHE=0
	!FHTG=0
	!FACT=1
	!DC=1

!Real Constants

R,102,, , FKN, FTOLN, , PINB,
RMORE,, , TAUMAX, , , FKT,

R,101,, , FKN, FTOLN, , PINB,
RMORE,, , TAUMAX, , , FKT,

!Real Constants
!R,3,R1, R2, FKN, FTOLN, ICONT, PINB,
!RMORE,PMAX, PMIN, TAUMAX, CNOF, FKOP, FKT,
!RMORE,COHE, TCC, FHTG, SBCT, RDVF, FWGT,
!RMORE,ECC, FHEG, FACT, DC, SLTO, TNOP,
!RMORE,TOLS, MCC, PPCN, FPAT, COR, STRM
!RMORE,FDMN, FDMT, FDMD, FDMS, TBND

KEYOPT,101,2,1						!Контактный алгоритм
KEYOPT,101,5,0						!автоматическое закрытие открытых контактов
!KEYOPT,101,6,2						!?
KEYOPT,101,7,0						!доп контроль шага нагрузки, уменьшение шага при высоком проникновении
KEYOPT,101,9,1						!начальный пенетрации борьба
KEYOPT,101,10,2   						!когда обновлять контактную жесткость
KEYOPT,101,12,0
KEYOPT,100,5,0						!Определяет тип контактируемых поверхностей объём/объём, оболочка/оболочка и т.д

! Generate the target surface1   
vsel,s,loc,y,a,hpl+a

eslv, s
aslv,s
asel,r,loc,y,a,a
NSLA,S,1							!Выбор узлов принадлежащих area 
ESLN,r,0							!Выбор элементов которым принадлежат выбранные узлы
TYPE,100							!
REAL,101
ESURF							!Создаём элементы контакта на поверхности

! Generate the contact surface 1 
vsel,s,loc,y,0,a
eslv, s
aslv,s
asel,r,loc,y,a,a
NSLA,S,1							!Выбор узлов принадлежащих area
ESLN,r,0							!Выбор элементов которым принадлежат выбранные узлы
TYPE,101							
REAL, 101
ESURF							!Создаём элементы контакта на поверхности
ALLSEL							

! Generate the target surface2   
vsel,s,loc,y,-hpl,0
eslv, s
aslv,s
asel,r,loc,y,0,0
NSLA,S,1							!Выбор узлов принадлежащих area 
ESLN,r,0							!Выбор элементов которым принадлежат выбранные узлы
TYPE,100							!
REAL,102
ESURF							!Создаём элементы контакта на поверхности

! Generate the contact surface 2 
vsel,s,loc,y,0,a
vsel,r,loc,x,-a1/2,0
vsel,r,loc,z,-a/2,0
eslv, s
aslv,s
asel,r,loc,y,0,0
NSLA,S,1							!Выбор узлов принадлежащих area
ESLN,r,0							!Выбор элементов которым принадлежат выбранные узлы
TYPE,101							
REAL, 102
ESURF							!Создаём элементы контакта на поверхности
ALLSEL							

CNCHECK,AUTO
fini

/sol
/gst,1
/NERR, 10000,
nlgeom,on
nldiag,nrre,on
antype,0
autots,on
OUTRES,ALL,ALL  
!EQSLV,PCG,1E-5 
nropt,unsym,,on
lnsrch,on
!pred,on
cutcon,noshape,1
cutcon,plslimit,0.3
neqit,10
nsubst,400,10000,100
cnvtol,f,,0.005
cnvtol,u,,0.002

asel,s,loc,y,a+hpl,a+hpl
Da,all,Uy,-a*0.005
alls
time, 100
solve, , , , , nocheck

!*do,i,1,100
!asel,s,loc,y,a+hpl,a+hpl
!Da,all,Uy,-a*0.005*i/100
!cnvtol,f,,0.05+i/100*0.45
!time, i
!solve, , , , , nocheck
!*GET, dal, ACTIVE, 0, solu, cnvg,
!*IF, dal,eq, 0, exit,:label1
!*enddo
!:label1

fini
/post1
alls
set,last
*GET, sets, ACTIVE, 0, set, nset
sets=sets-1
*dim, react_sum1, array, sets,,,Time
*do,i,1,sets
SET,,, ,,, ,i 
nsel,s,loc,y,-hpl
*get,node_count,node,0,count
*do,j,1,node_count
*get, ni, node, 0, num, min
*get, react, node, ni, rf, fy
react_sum1(i)=react_sum1(i)+react
nsel,u,node,,ni
*enddo
*enddo
alls

fini

на схемах представлены значения реакции. величина не сошедшехся сил и их расположение

[Stanum]

Прошу прощения за двойной пост, все приложения не влезали
Прикладываю output для 2 задач и ещё одну картинку.
здесь на мой взгляд возможный источник проблемы напряжения по элементам имеют странный вид: они меняются очень резко от элемента к элементу

[MrWhite]

материал для них одинаковый: только tb,kinh без tb, conc/

Запустил макрос, просматриваю кривые strain-stress существенно отличаются для материалов 1 и 6.
По элементному решению которое вы считаете корнем проблемы не беспокойтесь, это жертва скорости решения.

[Stanum]

MrWhite, материал оба раза 6. 1 с ниспадающей ветвью пока не использовал.

Сообщение от MrWhite Цитата: По элементному решению которое вы считаете корнем проблемы не беспокойтесь, это жертва скорости решения.

Вы имеете ввиду, что у меня сетка крупная?

Просто, при использовании 185, таких скачков нет.

[frostyfrost]

Такой вопрос, а какую диаграмму KINH деформирования задавать для цилиндров, призм и т.д.? Ведь в РФ - диаграмма задается в призменной прочности, в ЕС - цилиндрическая. По идее должна быть одна диаграмма общая, поскольку в зависимости от того, что считаем (призму, куб, цилиндр), на выходе должны иметь соответствующую прочность, как этот момент можно увязать. Плюс вопрос в том, что по KINH и т.п. условие пластичности не совсем применимо к бетону.

[MrWhite]

Я простите не понимаю суть Вашей проблемы.
Вы сравниваете одинаковый материал используя разные конечные элементы? Тогда зачем Вам 65-ый? Он не фига не аналог 185-ого. И зачем Вам менять опции 185-ого, когда в хелпе ясно описано какая для чего нужна. Какие вы там реакции на графике смотрите я тоже не понимаю. Опишите как Вы ее находите.

[Stanum]

MrWhite,спасибо за участие в теме, считал, нелинейно с 65 бетон, после разрушения, решил проверить, при какой нагрузке элемент разрушился, для каждого set, суммировал реакции в закреплении нижней пластины, в итоге имею зависимость суммарной реакции от прилагаемого смещения верхней пластины(времени). Умножав пиковую реакцию на 4 и поделив на площадь куба хотел бы получить ту прочность, что задал в начале. Не получил. Решил проверить где ошибся задал кубу 185 элемент, убрал tb, conc, проверил нагрузку снова - получил. Проверил тоже самое но с 65, не получил + вдобавок значение реакции менялось в зависимости от количества шагов.

Опции менял, для того, чтоб попытаться получить результат как с 65, не смог этого добиться. Я просто прочитав одну статью почему-то решил, что это из-за shear locking в элементах, хотя изгиба у меня нет.

[MrWhite]

А... вот теперь более понятно (вы меня простите, Ваша задача очень долго считается на моем калькуляторе). Тогда такой вопрос вы когда 65-ым делали смотрели картину трещинообразования?
И да, как и правильно отметил frostyfrost, закон KINH работает по критерию мизеса, закон CONC (для 65-ого) это сочетания 4-х отличающихся законов прочности. Т.е. как Вы их так просто собираетесь сравнить - непонятно. Если пластика в бетоне наступит раньше чем в KINH (что очень вероятно) то и реакция будет конечно другой сразу (вы же нагрузку перемещениями прикладываете?).

[Stanum]

MrWhite, да трещины смотрел (когда использовал материал бетон), общий вид получал 2 пирамиды НО, что мне не понравилось, трещины на первом этапе проявлялись :
А) очень резко = одна трещинка на грани куба следующий шаг почти весь куб в трещинах
Б) проскакивали узлы, то есть появлялись в тех же местах, что и скачки на картинке, напряжений, что я выложил, то есть сначала к примеру во всех 65 трещит к примеру 8 узел, а, потом трещит и все остальное,
Поэтому и решил проверить, не ошибся ли где-то. 65 без законов , по идее должен быть = обычному 8 узловому объемнику, а не получился.

На ваш вопрос по совмещенную законов, пока не могу ответить, я считал, что поверхность Вилама - Варнаке, определит поверхность разрушения, а kinh, путь к этой поверхности, наверно, я здесь не прав, но это не должно влиять на задачу тестовую, так как поверхность Варнаке, я в ней не задал. И 65 должен работать только по закону kinh.

Да: нагрузка перемещением верхней пластины.

[MrWhite]

Stanum, я попробую воспроизвести то что Вы делали:
- сделали 65 + tb,conc: получили трещины, но Вам трещинообразование не понравилось;
- сделали 65 c tb,kinh: не получили реакции какие ожидали и дальше пошли проверять с 185-ыми.

Ну тут вообщем 3 замечания:
1) трещино-образование нужно обязательно смотреть совместно с относительными деформациями (в случае с 65 линейными, т.е. трещинообразование он как пластику не понимает). Оценивать развитие трещин нужно именно по второму... Трещина в 65 означает наступление предельного состояние в элементе, но вовсе не означает, что Вы визуально эту трещину сможете лицезреть. Т.е. еще раз трещина в 65 это отображение наступления момента критерия прочности + площадка по которой оно происходит, не более того.
2) 65 без законов должен вести себя как родоначальный 45. По идее... Но я в хелпе Ansys не нашел прямой информации об этом. Поэтому это утверждение нужно верифицировать вообще-то. И в общем случае эта логика может быть не верной.
3) Однозначно Вам могу сказать что поверхность прочности бетона нельзя выразить через kinh. Если бы все так просто было . Вам достаточно посмотреть вид этих двух поверхностей в 3D чтобы понять что они совершенно не идентичны. Бетон еще можно попытаться выразить через DP или EDP - там хотя бы есть некое подобие у критериев прочности.

Ну и так по мелочи: советую в настройках 65-ого элемента порыться. Некоторые кейопты лучше изменить с дефолтных, какие - думаю сами разберетесь.

[frostyfrost]

Попробуйте без контакта решить. Низ закрепляете по xyz, верх - xy+перемещение по z. Потом просуммируйте реакции по нижней грани цилиндра. Для трещин хорошо бы добавить шагов. Если верх крепить не будете, то вместо пирамид должна быть наклонная трещина

[Stanum]

MrWhite,

Сообщение от MrWhite 65 без законов должен вести себя как родоначальный 45. По идее... Но я в хелпе Ansys не нашел прямой информации об этом. Поэтому это утверждение нужно верифицировать вообще-то. И в общем случае эта логика может быть не верной.

из справки второй абзац: solid 65= structural solid но с возможностями трещать и раздавливаться.
Справка после таблицы где описываются константы для таблицы tb, conc : если таблицы нет: трещать и раздавливаться не будет. То есть суммируя должны получать 45.

Сообщение от MrWhite 1) трещино-образование нужно обязательно смотреть совместно с относительными деформациями (в случае с 65 линейными, т.е. трещинообразование он как пластику не понимает). Оценивать развитие трещин нужно именно по второму... Трещина в 65 означает наступление предельного состояние в элементе, но вовсе не означает, что Вы визуально эту трещину сможете лицезреть. Т.е. еще раз трещина в 65 это отображение наступления момента критерия прочности + площадка по которой оно происходит, не более того.

Спасибо, за этот момент не знал, смотрел значит не туда.

По поводу поверхностей: edp и варнаке очень похожи по виду (нечто похожее на конус, ось - гидростатическое давление).
Edp и miso можно комбинировать если верить справке, + я среди макросов, даже в этой ветке видел комбинирование tb, conc и tb, kinh

И все же, я не понимаю как задав прочность бетона x, я не могу получить реакцию ей соответствующую при использовании 65, хоть с tb, conc хоть tb, kinh, хоть с их совмещением. У меня все таки не сложное НДС , конструкции.

По поводу кейоптов 65: спасибо, покручу, может, что и поможет.
frostyfrost, так тоже пробовал, сделаю ещё раз выложу результаты.
Мне не понравились, так-как вылазили пики напряжений на рёбрах куба.

[frostyfrost]

Реакция может не совпадать по соедующим моментам:
- плохо контакт считается (поэтому можно тестово через dof считать)
- модель не способна решить данный пример корректно
-в качестве fc задаете призменную или кубиковую прочность, хотя надо циллиндрическую (для большинства моделей американского происхождения базовая прочность цилиндрическая)
Несколько моделей, которые смотрел не в Ansys давали одинаковую кубиковую и призменную прочность при моделировании.

[Stanum]

frostyfrost,

Сообщение от frostyfrost плохо контакт считается (поэтому можно тестово через dof считать)

то есть контакт плохо считается только для solid65? Тот же контакт на 185 и всё идёт, но я все равно проверю, спасибо

Сообщение от frostyfrost модель не способна решить данный пример корректно

Не понял. немного, разъясните, пожалуйста, что Вы подразумевали?

Сообщение от frostyfrost в качестве fc задаете призменную или кубиковую прочность, хотя надо цилиндрическую

Опять же solid185 и 65 для которых задан только tb, kinh не дали одинакового результата. То есть если б 65 и 185 мне дали в tb, kinh, равенство, хоть приблизительное +/- 5-7 % я б даже не обратил внимание, но разница то большая 20-30%. И тут поверхность разрушения для бетона не причем, так как она даже не задана.

[MrWhite]

Сообщение от Stanum из справки второй абзац: solid 65= structural solid но с возможностями трещать и раздавливаться. Справка после таблицы где описываются константы для таблицы tb, conc : если таблицы нет: трещать и раздавливаться не будет. То есть суммируя должны получать 45.

Я же говорю - прямой отсылки к тому что это эквивалент нету. Кроме того для 65 родной является только модель conc - он и создан под нее, поэтому изменение 45-ого могли зайти достаточно далеко. Например в линейных задачах он может вести себя нормально (это я проверял), а вот при использовании сторонней модели нелинейности уже далеко не факт. Вы в справке читали о совместимости нелинейных моделей tb,conc и tb,kinh?

Реакция может не совпадать по соедующим моментам: - плохо контакт считается (поэтому можно тестово через dof считать) - модель не способна решить данный пример корректно -в качестве fc задаете призменную или кубиковую прочность, хотя надо циллиндрическую (для большинства моделей американского происхождения базовая прочность цилиндрическая) Несколько моделей, которые смотрел не в Ansys давали одинаковую кубиковую и призменную прочность при моделировании.

- не в контакте точно, там по другому проблемы выглядят;
- вероятно (например так описано выше);
- ну точно не в этом дело. В тестовом примере простейший кубик, при чем тут прочность цилиндра?

[frostyfrost]

Попробуйте решить модель без контакта плюс еще одну из одного кэ из 65го и 185го, чтобы понять, что это не влияние контакта.

[Stanum]

frostyfrost, попробовал, вот код. результат тот же.

Код:

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

 finish
/CLEAR,NOSTART  
/units, si
/filnam, kinh
/title, kinh
/OUTPUT, kinh_1, txt, , D:\ANSYS files
!бетон не обжатый
!параметры
fc=40e6!прочность не обжатого бетона максимальная (М)Па
fct=2.1e6! прочность бетона на растяжение (М)Па
ec1=0.00186! напряжения соответствующие макс напряжениям
ecu=0.003! предельные напряжения в необжатом бетоне
Eb=32e9! модуль упругости бетона (Г)Па
nub=0.2
an=15
*dim, eps_b, array, 20
*dim, sig_b, array, 20
eps_b(1)=0.3*fc/Eb
eps_b(an)=ec1
eps_b(20)=ecu

*do,i,2,an
eps_b(i)=eps_b(i-1)+(eps_b(an)-eps_b(1))/(an-1)
*enddo
*do,i,an+1,20
eps_b(i)=eps_b(i-1)+(eps_b(20)-eps_b(an))/(20-an)
*enddo

sig_b(1)=0.3*fc
*do,i,2,20
x=eps_b(i)/ec1
Ebc1=5e6*(fc**(1/2))
Esec1=fc/ec1
r1=Ebc1/(Ebc1-Esec1)
sig_b(i)=fc*x*r1/(r1-1+x**r1)
*enddo

!бетон без ниспадающей ветви
*dim, eps_bn, array, an+1
*dim, sig_bn, array, an+1
*do,i,1,an
eps_bn(i)=eps_b(i)
sig_bn(i)=sig_b(i)
*enddo
eps_bn(an+1)=eps_b(20)
sig_bn(an+1)=sig_b(an)*1.01

!сталь
nus=0.3! коэффициент пуассона стали
Es=210e9!модуль упругости стали, (Г)Па


/prep7
!тип элемента
et,1,185
et,2,65
et,3,185
keyopt,1,2,3
keyopt,2,3,2
keyopt,2,5,2
keyopt,2,6,3
keyopt,2,7,1

keyopt,3,2,0

!материала
mp,ex,5,Es
mp,nuxy,5,nus

mp,ex,6,Eb
mp,nuxy,6,nub
tb,kinh,6,1,an+1,0 
tbtemp,0
*do, i, 1, an+1
tbpt, ,Eps_Bn(i), Sig_Bn(i)
*enddo

!построение модели
a=0.15
a1=0.151
hpl=0.05
apl=0.2
apl1=0.201
block,0,-a1/2,0,a,0,-a/2

!граничные условия
asel,s,loc,y,0,0
Da,all,uy,0
asel,s,loc,x,0,0
asel,a,loc,z,0,0
DA,all,Symm
alls

!сетка
elem=a/20
LESIZE,ALL,elem, , , ,1, , ,1,  
alls
vsel,s,loc,y,0,a
VATT, 6, 1, 2, !(3)

alls
MSHAPE,0,3d 
MSHKEY,1
VMESH,all
alls
/sol
/gst,1
/NERR, 10000,
nlgeom,on
nldiag,nrre,on
antype,0
autots,on
OUTRES,ALL,ALL  
EQSLV,PCG,1E-5 
!nropt,unsym,,on
lnsrch,on
pred,on
cutcon,noshape,1
cutcon,plslimit,0.3
neqit,10
nsubst,400,10000,100
cnvtol,f,,0.05
cnvtol,u,,0.02
asel,s,loc,y,a,a
Da,all,Uy,-a*0.005
alls
time, 100
solve, , , , , nocheck

fini
/post1
alls
set,last
*GET, sets, ACTIVE, 0, set, nset
sets=sets-1
*dim, react_sum1, array, sets,,,Time
*do,i,1,sets
SET,,, ,,, ,i 
nsel,s,loc,y,0
*get,node_count,node,0,count
*do,j,1,node_count
*get, ni, node, 0, num, min
*get, react, node, ni, rf, fy
react_sum1(i)=react_sum1(i)+react
nsel,u,node,,ni
*enddo
*enddo
alls

fini

[frostyfrost]

Stanum, посчитал свои схемки без контакта и без KINH и аналогичных зависимостей, аналогично получилась несущая способность ниже F_c (правда и 65-ми я особо никогда не пользовался, поскольку считаю в основном сталь в ANSYS). Для куба, цилиндра и призмы получились разные несущие способности, что хорошо, но по значениям ниже F_c. Несущую способность определял по суммированию опорной реакции по Z с последующим делением на площадь опирания. Буду смотреть дальше в чем мой косяк или не косяк.

[Stanum]

frostyfrost,спасибо, что продолжаете участвовать в теме. Я пока проверил свою схема но не с kinh, а с bkin результат совпадает с желаемым.

[frostyfrost]

Сообщение от Stanum а с bkin результат совпадает с желаемым.

Для куба, призмы и цилиндра результаты совпадают с соответствующими прочностями при одной и той же кривой напряжения-деформации?

[Stanum]

frostyfrost, я считал только куб, говорю о совпадении 185/65
проверил kinh, с меньшим количеством точек, результат тоже совпадает.
Проблема в диаграмме выходит.

[frostyfrost]

А))) Понял))))

[Stanum]

Сравнил результат для задачи из solid65 с kinh при разном количестве позиций в таблице (от 6 до 11).
Получил следующий результат : при количестве точек в таблице меньше либо равном 6 решение получаем верное.
При увеличении количества шагов результат меняется, постепенно уменьшаясь.
В чём причина такого поведения элементов не знаю, может какие-то не обозначенные ограничения в solid65?

При использовании изотропного упрочнения (MISO) результат остаётся правильным независимо от количества точек на кривой.

Вот листинг тестовой задачи:

Код:

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

 finish
/CLEAR,NOSTART  
/units, si
/filnam, kinh
/title, kinh
/OUTPUT, kinh_1, txt, , D:\ANSYS files
!бетон не обжатый
!параметры
fc=40e6!прочность не обжатого бетона максимальная (М)Па
fct=2.1e6! прочность бетона на растяжение (М)Па
ec1=0.00186! напряжения соответствующие макс напряжениям
ecu=0.003! предельные напряжения в необжатом бетоне
Eb=32e9! модуль упругости бетона (Г)Па
nub=0.2
kolmax=20
kolmin=6
ele=kolmax-kolmin+1

*dim, eps_b, array, kolmax,ele
*dim, sig_b, array, kolmax,ele

*do,kol,kolmin,kolmax
j=kol-kolmin+1
an=nint(kol/2)+1

eps_b(1,j)=0.3*fc/Eb
eps_b(an,j)=ec1
eps_b(kol,j)=ecu

*do,i,2,an
eps_b(i,j)=eps_b(i-1,j)+(eps_b(an,j)-eps_b(1,j))/(an-1)
*enddo

*do,i,an+1,kol
eps_b(i,j)=eps_b(i-1,j)+(eps_b(kol,j)-eps_b(an,j))/(kol-an)
*enddo

sig_b(1,j)=0.3*fc

*do,i,2,kol
x=eps_b(i, j)/ec1
Ebc1=5e6*(fc**(1/2))
Esec1=fc/ec1
r1=Ebc1/(Ebc1-Esec1)
sig_b(i, j)=fc*x*r1/(r1-1+x**r1)
*enddo

*enddo

!бетон без ниспадающей ветви
anmax=nint(kolmax/2+1)
*dim, eps_bn, array, anmax+1,ele
*dim, sig_bn, array, anmax+1,ele

*do,kol,kolmin,kolmax
j=kol-kolmin+1
an=nint(kol/2)+1

*do,i,1,an
eps_bn(i, j)=eps_b(i, j)
sig_bn(i, j)=sig_b(i, j)
*enddo

eps_bn(an+1, j)=eps_b(kol, j)
sig_bn(an+1, j)=sig_b(an, j)*1.001

*enddo

!сталь
nus=0.3! коэффициент пуассона стали
Es=210e9!модуль упругости стали, (Г)Па

/prep7
!тип элемента
et,1,185
et,2,65
et,3,185
keyopt,1,2,3
keyopt,2,1,1
keyopt,2,3,2
keyopt,2,5,1
keyopt,2,6,3
keyopt,2,7,1
keyopt,3,2,0

!материала
mp,ex,6,Eb
mp,nuxy,6,nub

!построение модели
a=0.15
a1=0.151
hpl=0.05
apl=0.2
apl1=0.201
block,0,-a1/2,0,a,0,-a/2
alls

!граничные условия
asel,s,loc,y,0,0
Da,all,uy,0
asel,s,loc,x,0,0
asel,a,loc,z,0,0
DA,all,Symm
alls

!сетка
elem=a/20
LESIZE,ALL,elem, , , ,1, , ,1,  
alls
VATT, 6, 1, 2, !
alls
MSHAPE,0,3d 
MSHKEY,1
VMESH,all
alls

*dim, react_sum1, array, 120,ele,

*do, kol, kolmin, kolmax
j=kol-kolmin+1
an=nint(kol)/2+1

/prep7
TBDE,KINH,6 
tb,kinh,6,1,an+1,0 
tbtemp,0
*do, i, 1, an+1
tbpt, ,Eps_Bn(i, j), Sig_Bn(i, j)
*enddo

/sol
SOLCONTROL,,,INCP
/gst,1
/NERR, 10000,
nlgeom,on
nldiag,nrre,on
antype,0
autots,on
OUTRES,ALL,ALL  
!EQSLV,PCG,1E-5 
nropt,unsym,,on
lnsrch,on
pred,on
cutcon,noshape,1
cutcon,plslimit,0.3
neqit,10
nsubst,400,10000,100
cnvtol,f,,0.05
cnvtol,u,,0.02

asel,s,loc,y,a,a
Da,all,Uy,-a*0.005
alls
time, 100
solve, , , , , nocheck

fini
/post1
alls
set,last
*GET, sets, ACTIVE, 0, set, nset
*do,ik,1,sets
SET,,, ,,, ,ik 
nsel,s,loc,y,0
*get,node_count,node,0,count
*do,jk,1,node_count
*get, ni, node, 0, num, min
*get, react, node, ni, rf, fy
react_sum1(ik, j)=react_sum1(ik, j)+react
nsel,u,node,,ni
*enddo
*enddo
alls
fini

*enddo

[miko2009]

Пробовал кто нибудь моделировать в ANSYS worckbench кусок перекрытия с колоннами + вся арматура не дискретная, а физическая + моделирование контакта арматуры и бетона ?
Для бетона принял:

ET,matid,SOLID65
MPDATA,EX,matid,,2e4
MPDATA,PRXY,matid,,0.25
TB,CONC,matid,1,9,
TBTEMP,0
TBDATA,,0.2,0.8,3,30,,
TBDATA,,,,,,,
R,matid, , , , , , ,
RMORE, , , , , , ,
RMORE, ,

арматура = сталь по умолчанию + нелинейные внутренние параметры из библиотеки.

как замоделить контакт ?

[derik]

а какой тип элемента: link/beam или reinf?

[miko2009]

solid185
с link180 проблем нету.
Когда бетон SOLID65 ругается:
Unconverged решение (определены в качестве подшаге 999999) выводится для анализа отладочных целей. Результаты в это время не должны использоваться для любых других целей.
КОгда бетон тот же 185 или 186 все отлично, но без развития трещин не итересно

[derik]

попробуйте узу-узел, сделав сетку соответствующую. а зачем Вам контакт бетона и арматуры, мышей ловите? только сейчас увидел про ВБ :/ попробуйте так, а там не знаю, ВБ не открываю

[miko2009]

Сообщение от 6pa4o а зачем Вам контакт бетона и арматуры, мышей ловите?

лень двигатель прогресса, и именно прогресс двигает меня учитывать контакт ,благо при помощи Revit-a на модель уходят считанные минуты и потом экспорт в SAT формат со всей арматурой. Вот и задался вопросом как сделать контакт. это же не просто коэффициент трения, и анзис меня материт что контакт не возможен при таких настройках материалов , и у меня закрадывается ощущение что SOLID65 не работает с другими объемными элементами, надеюсь что я не прав, я же не единственный кто задался таким вопросом ? кто-то хотя бы для себя уже решил такую задачу, вот хочется услышать мнение этих людей.

[derik]

я в 65-й свободно вставляю 264, макросом. контактом никогда в данной постановке не пользовался. и сколько читал диссеров не делают этого, мб не так перевел (я слаб в англицком)

[Stanum]

miko2009,
Вот статейка, и вот гугль

[miko2009]

так я же не про контакт между solid (объемным НДС) и link (линейным НДС) а про solid\solid в частности solid65\solid185
(какой КЭ по умолчанию в ВБ ? ) 186

[miko2009]

6pa4o Теперь лучше ?

[derik]

да а кстати - было бы проще по времени расчета и по наложению сетки - моделировать арматуру линком/бимом/реинфом. а то получаеться что все объемниками, тем более арматура
я что то сразу не сообразил что арматура объемником....

[miko2009]

Так это понятно что проще у меня задачка сегодня весила 2 000 000 КЭ но ее дополнил еще маленьким куском здания, ожидаю в районе 7 000 000-10 000 000, но в течении часа думаю уложится посчитать , на создание модели у меня уходят минуты ...... так что пока вы в классике строите модель я ее уже обрабатываю по моим прикидкам когда я закончу экспериментировать с различными настройками сеточного препроцессора и настройками то на такую задачу будет уходить не более 3 часов с полным анализом. Пока что отдельные узлы и конструкции здорово получаются а вот с кусками здания все под вопросом , много чего нужно учесть заранее чего нету ни в одной книге......Кстати даже контакт замоделить в такой задаче как само собой разумеющиеся, просто вычесть из модели бетона всю арматуру чтобы получить грани по которым происходит CONTA и TARGE и это все достигается пару движениями в самом ВБ. Я все же надеюсь что у такого подхода намного больше арсенал выходных данных.
А если еще кто нибудь поделится макросом для ВБ по crack вообще расцелую

[Stanum]

miko2009, а разве нормы допускают проскальзывание арматуры?
И еще вопрос что у вас за зверь, что вы за три часа такую громаду считать планируете? У меня значительно меньшая модель до 300к элементов считалась при учете только геом нелинейности около 10 минут.

[derik]

miko2009, а что трещинообразование не показывается?

[miko2009]

Сообщение от Stanum miko2009, а разве нормы допускают проскальзывание арматуры?

прочность сцепления

Сообщение от 6pa4o miko2009, а что трещинообразование не показывается?

все показывается , с трещинами я еще не пробовал но вот пользователь Дмитрий Фролов показал что это возможно в ВБ
http://www.cae-club.ru/content/napryazheniya-v-betone
http://forum.dwg.ru/showthread.php?t=76274
offtop: старый уже зверь, 16 ГБ, GTX780, i7 -4 ГГц.

[miko2009]

Пришел к выводу , моделировать ЖБ элементы Solid65 (бетон) и solid187(арматура) невозможно с учетом развития трещин (точнее можно но с ОЧЕНЬ грубой сеткой как на 1 скрине) так как при разбиении тетраэдрами всего 5 узлов , для моделирования SOLID65 необходима 8 узловая фигура. Такая модель работает но без развития трещин когда бетон замоделирован так же 187 солидом. Так же и по времени и по размерности задача все же не реальная , проверял обычную пустотку (1 000 0003 700 000 элементов) хотя результаты правдивые но без трещин . Печально.........

[miko2009]

Сообщение от 6pa4o miko2009, а что трещинообразование не показывается?

Нажмите для просмотра

http://www.youtube.com/watch?v=OAfZQsiXw5gэто пока что предварительно....

[derik]

Сообщение от miko2009 https://www.youtube.com/watch?v=OAfZ...ature=youtu.be это пока что предварительно....

пишет что видео с ограниченным доступом

[miko2009]

исправил

[miko2009]

Итак прошло уже несколько месяцев как я взялся за Жб в WB. Пытался классику передрать в надстройку, столкнулся с интересной задачей что возможно считать любые конструкции при помощи ANSYS WB, но вопрос был в том как задавать арматуру в конструкциях различной формы , как считать узлы ЖБ, решил это просто, все считать солидами , так как справка рекомендует 185(186) в качестве бетона и арматуры я решил что почему бы и нет.Долго 65 и 185 конфликтовали друг с другом в одной модели (по крайней мере у меня).Но все же получились вот такие результаты. Итак исходные данные:
Есть балка пролетом 2770, армирована только 2 стержнями 12д А500С, сечение балки 200х300 В25, рас.до ц.т.арматуры 50 мм. нагрузка 4Т на расстоянии 1.1м от опоры , балка шарнирная.
1. Расчет в арбате показал 3.5 Т если игнорировать ширину раскрытия трещин, прогиб 7.62мм
2. Модель классическая бетон SOLID65, арматура LINK180 результаты такие:

Нажмите для просмотра

http://www.youtube.com/watch?v=OAfZQsiXw5g&list=UUb76hQTe1cgAMa8pVyjfwsgпрогиб:

Нажмите для просмотра

http://www.youtube.com/watch?v=3ZlC7oHptfY&feature=youtu.be3. Собственно обсуждаемая модель(хотя насколько я понял и первая тоже будет критиковаться ) бетон SOLID65, арматура SOLID185 :

Нажмите для просмотра

http://www.youtube.com/watch?v=pWsiMcWpRsEпрогиб:

Нажмите для просмотра

http://www.youtube.com/watch?v=l63_vgxrhDc&feature=youtu.beнапряжение в арматуре:

Нажмите для просмотра

http://www.youtube.com/watch?v=JQCPzZ_Kh2M&feature=youtu.beво всех моделях каждая секунда = 0.5 Тс , 8 секунд (счетчик с верху в углу) = 4 Тс.
Во всех моделях нету TB,KINH, арматура BISO во всех случаях. Результат 3 модели ОЧЕНЬ близок к результату АРБАТа.

Так же в качестве интереса нашел испытания плиты ПК60.15 , по описанию выдерживает 10 Тс , при натуральных испытаниях нагрузкой 9.8 Тс дала прогиб 70мм, замоделировал эту же плиту получил 68 мм, напряжения в арматуре 458 МПа. (забыл на работе)
Кто что думает по этому поводу ?
P.S. если честно возможно есть агрехи в схемах или в работе конструкций так как у меня их уже сотня обьемом 150 ГБ, пробовал различные сеточные модели , различные настройки , различные граничные условия и тд и тп что сам начал запутываться.

[miko2009]

что совсем тема тухлая ? неужели некому поделится своим опытом ? лично у меня такие ощущения что я постиг темную сторону силы
я тут кусочек здания считал (одна колонна и перекрытие) замороченный узел стыка, в калькуляторах ни как не давался. 1 200 000 элементов , сходимость по силе 5%, считался 6 суток на офисном ПК в не лине и не досчитал , свет вырубали днем:
зато упругий 10 минут (бетон на изображениях скрыл).
Теперь я понял что для такой работы инженеру нужен нормальный ПК.........SSD накопители , огромные вычислительные мощности (i7ext.,xeon на базе материнок двухпроцессорных/трех/четырех........),64+оперативы,quadro K5000\K6000 или эквивалент tesla. На всем что ниже только только такие не хитрые вещи считать , аля балка , плита и тд. Может есть способы ускорить работу ? когда WB считает нелин он похоже напрягает не ЦП а веник ?! SSD винчестер я так понимаю в корне ускорит процесс ?
6pa4o ну что там с пустоткой ? не пробовал считать ? я кстати ошибся, прогиб 48 мм у меня и в испытаниях нагружали плиту с "дефектами" , так что 70 явный перебор.

[genya]

Вопрос не по теме, как учесть депланацию при расчете тонкостенных стержнeй?

[Stanum]

genya, если вы о beam элементах посмотрите в кейоптах, там есть пункт отвечающий за стесненное кручение
Keyopt(1)=1

[genya]

Извиняюсь, что не уточнил, Ansys WB? С классикой пока не очень!

[miko2009]

тоже самое и в WB , добавить commands для элемента:

ET,1,BEAM189
KEYOPT,1,1,1
ну и дальше по желанию......

[genya]

Спасибо, вставил эту команду, но что-то ни чего не изменяется, для просмотра напряжений в дерево проекта вставил Linearized Equivalent Stress в окне Details of выбрал Membrane+Bending. Пытаюсь сравнить расчет по методичке и МКЭ. Скажите, а что еще нужно учитывать при расчете тонкостенных стержне? Потому как, что-то не срастается!

Сообщение от genya Спасибо, вставил эту команду, но что-то ни чего не изменяется, для просмотра напряжений в дерево проекта вставил Linearized Equivalent Stress в окне Details of выбрал Membrane+Bending.

Сообщение от genya Потому как, что-то не срастается!

А что Вы хотите увидеть? Что там в методичке такого, чего не видно в WB?

[genya]

Хотелось бы увидеть сходимость решения по нормальным напряжениям, хотя бы с небольшим отклонением!

[genya]

Вот эпюра напряжений нижнего пояса, так-же прикрепляю методичку, на странице 27 построена эпюра нормальных напряжений:
в)

Сообщение от genya на странице 27 построена эпюра нормальных напряжений:

а Вы смотрите эквивалентные. где логика?

[genya]

Спасибо palexxvlad! Был не внимателен, учусь!
Теперь разница составила 4,5% все еще не удовлетворительно, грубая ошибка найдена, осталось уточнить результат.
В чем еще может быть проблема?

А вот по примеру №2 сходимость по нормальным напряжениям на 100%!

Начал проверять касательные напряжения, указав путь в месте наибольших поперечных сил, но что-то не сходится, возможно, нужно применять другой подтип и Membrane+Bending здесь не подходит, потому, как Peak дает более близкий результат.

Кто-то может объяснить применение именно ПОДТИПА относительно нормальных и касательных напряжений?

Можно ли в Ansys WB посмотреть угол закручивания балки?

Сообщение от genya В чем еще может быть проблема?

Может быть неточность в определении геометрических и секториальных характеристиках путем перемножения "ниточных" эпюр координат по методичке. В местах стыковки элементов сечения материал там учитывается дважды.

[genya]

Я думаю с неточностями я сам разберусь, больше интересует остальные вопросы. Еще один момент, а можно ли в Ansys WB посмотреть геометрические характеристики сечения, я имею ввиду не только момент инерции или момент сопротивления сечения, хотя и их я найти не могу?

Сообщение от genya Еще один момент, а можно ли в Ansys WB посмотреть геометрические характеристики сечения, я имею ввиду не только момент инерции или момент сопротивления сечения, хотя и их я найти не могу?

Скорее всего нет.
Вы бы выложили свою схемку в формате WB.

[genya]

Пожайлуста!

Мда-а-а... Так Вы что же, на солидах считаете "тонкостенку"? Причем же здесь beam189 со всеми "кейоптами"? на объемниках достаточно трудно будет получить сходимость результатов по касательным напряжениям при такой сетке КЭ.

[genya]

Можно чуть поподробнее с типом конечных элементов, так как с этим я еще не совсем разобрался. Уменьшить сетку не получается, оперативной памяти при расчете не хватает!

[miko2009]

сколько элементов ?

[GVI]

Добрый день. Увидел Вашу задачу, решил проверить в SOFiSTiK попутно с вопросом о программном обеспечении, позволяющем получить бимомент.
Пожалуйста, смотрите результаты.

Сообщение от genya Можно чуть поподробнее с типом конечных элементов, так как с этим я еще не совсем разобрался.

Нельзя. Чуть подробнее никак не получится. Это не табличка начальных параметров в методичке.

Сообщение от genya Уменьшить сетку не получается, оперативной памяти при расчете не хватает!

Еще бы. genya, дорогой, не майтесь дурью. Прошу прошения за мой французский.

----- добавлено через ~1 мин. -----
GVI, а софистик позволяет прикладывать нагрузки в произвольных местах линейного КЭ?

[genya]

Сообщение от palexxvlad Мда-а-а... Так Вы что же, на солидах считаете "тонкостенку"?

Прокомментируйте свой ответ.
Я же прописал в командной вставке:
ET,1,BEAM189
KEYOPT,1,1,1

Я, так понял, что Я использовал тип конечных элементов SOLID"..." а нужно использовать BEAM189!
Вы пишите, что я "считаю на солидах" Я понятия не имею какие элементы использовались, если я бы знал как их посмотреть и поменять то проблем не было.
Как-же исправить, где можно посмотреть используемый тип конечного элемента в Ansys WB?

genya, ну нельзя же так. Раз Вы понятия не имеете о том, с чем работаете, что хорошего из Вашей работы может получится? Ansys не дорогая иномарка, которую можно купить(или украсть) и ездить абсолютно не ведая как она устроена и какие принципы заложены в основу ее работы. Ищите курсы по Ansys или читайте книжки по МКЭ - станет понятно.

[genya]

Ваша позиция, стала мне понятна! Могли бы и сразу сказать я не хочу, я не знаю, я не буду. Я тогда зря бы время не тратил и отложил бы до лучших времен, но спасибо и на этом.

Сообщение от genya Ваша позиция, стала мне понятна!

Да, я не обучаю ANSYSу через форум, я лишь пытаюсь подсказывать там, где это уместно. Ваш случай не подходит для подсказок. genya, Вам нужно учится, учится и еще раз учится. Удачи.

[Stanum]

genya, начните с классики , там очень и очень подробный хелп. Откройте раздел basic analysis, и прочтите. После прочтения сделаете вашу задачу.

[GVI]

Сообщение от palexxvlad Нельзя. Чуть подробнее никак не получится. Это не табличка начальных параметров в методичке. Еще бы. genya, дорогой, не майтесь дурью. Прошу прошения за мой французский. ----- добавлено через ~1 мин. ----- GVI, а софистик позволяет прикладывать нагрузки в произвольных местах линейного КЭ?

Да, позволяет. Так и сделано.

[genya]

Спасибо конечно за совет Stanum, но если я сейчас начну решать эту задачу в классике то время пойдет не на часы, а на недели, а у меня столько времени нет. К тому же я уже решил эту задачу в WB, ошибки были пустяковые, конечно от незнания самого WB, так как изучать Ansys я начал совсем недавно. Так-же разобрался и с углом закручивания, у WB существует специальный инструмент для этого.
Тема закрыта.

[fasol]

Господа, а кто ходил в экспертизу с расчетами железобетона в нелинейной постановке?
Вы учитывали нормативные коэффициенты условия работ для бетона и арматуры? Разделяли первую и вторую группу предельных состояний с соответствующими расчетными сопротивлениями?

[miko2009]

ANSYS не совсем для экспертизы , там строят численные эксперименты которые желательно проверить натурными.

[MrWhite]

Сообщение от miko2009 ANSYS не совсем для экспертизы , там строят численные эксперименты которые желательно проверить натурными.

Вы серьезно? И что значит "не совсем для экспертизы"?

[miko2009]

Совсем серьезно, экспертиза смотрит 87 постановление как в настольную книгу и им все равно что там задано в ANSYS, им нужна бумажка. Увы в РФ ANSYS живет совместно с СТУ.

[MrWhite]

Вы же о составе проектной документации и их содержание? Что там ограничивает использование Ansys?

[derik]

Offtop:

Сообщение от MrWhite Вы же о составе проектной документации и их содержание? Что там ограничивает использование Ansys?

к сожалению эксперДы. им подавай лиц программы сертифицированные и верифицированные в РФ. хотя некоторых из них можно склонить на свою сторону доказав что объект нетривиальный и обычные РК здесь малоприменимы. но таких случаев я знаю всего 1-2

[MrWhite]

Сообщение от derik Offtop: к сожалению эксперДы. им подавай лиц программы сертифицированные и верифицированные в РФ. хотя некоторых из них можно склонить на свою сторону доказав что объект нетривиальный и обычные РК здесь малоприменимы. но таких случаев я знаю всего 1-2

Из всего этого Ansys только что не сертифицирован. Сертификации соответствия никак не относится к численным методам, только к методике расчетов, материалам и пр. Если мы говорим про нелинейное моделирование ж/б то нет ни одной программы сертифицированной под это дело, по причине отсутствия методики в нормах. Т.е. не на что сертифицировать. Суть то какая, нет никаких ограничений в применении Ansys и других подобных программ, просто очень кривая у нас нормативная база.
По своему опыту прохождения экспертизы (небольшому) по расчетам в Ansys на порядок меньше замечаний, чем по расчетам в Лире или ручным расчетам. Правда у меня не ПГС.

[miko2009]

Ну у вас же с гидросооружениями работа как я понял, а Лира и остальные калькуляторы геотехнику считают очень ....... точнее можно сказать пальцем в небо на коэф. постели и с объёмом там все очень туго, даже в Лира 10.4 в которой есть DP.
С ANSYS вопрос сертификации очень сложный , там же нету так таковых методик СНиП, только сотни и тысячи теорий. Вот вам для примера работа которая предстоит у меня на одну крупную строительную компанию, есть новая технология, методики расчета нету в СНиП, но у компании есть свой ЦНИИ, есть патенты на нее и есть материальная база для испытаний конструкций и не только. В итоге будем экспериментировать с материалами в ANSYS , проводить численные эксперименты , потом натурные испытания принципиальных ячеек здания. По результатам вносить изменения в настройки материалов (что бы получить такой же характер разрушения только численно), выводить из этого всего дела законы упрочнения экспериментальных материалов. Будем подключать НИИЖБ для того что бы они узаконили через СТУ и описали расчетную методику и только потом будем заходить в экспертизу.

[MrWhite]

miko2009, ну гидротехника это не геотехника, но вы правы в том, что реализовать необходимые расчеты с помощью Лиры и пр - невозможно. Как впрочем невозможно все сделать в том же Плаксисе, хотя он то сертифицирован на расчеты бетонных и железобетонных конструкций . Это к вопросу о статусе этих самых сертификатов - я лично их и в грош не ценю.

Изначально вопрос от fasol звучал так: кто ходил в экспертизу с расчетами железобетона в нелинейной постановке?
Я его воспринял (может ошибочно), как моделирование железобетона. Вот если говорить только про моделирование, то нет не ходил. Я выполнял расчеты в Ansys в соответствии с СП в том числе и армирование, по методике которую сам реализовал. Юридически это шло как ручной расчет в дополнение к определению НДС в Ansys.
Моделирование если и идет, то как дополнение уточняющее или проверяющее точные места. Но не как единственное обоснование.
В настоящее время, с текущими нормативами не вижу возможности отказаться от расчета армирования и перейти к его моделированию. Но это не отменяет того факта, что моделирование делать не нужно, во-первых это требование ФЗ №384, а во вторых дополнительный уровень самопроверки и надо сказать очень эффективный.

Ну и касательно экспертизы: для ГТС и объектов атомной энергетики - расчеты в Ansys это уже скорее обычное явление, чем исключение.

[miko2009]

MrWhite просто в калькуляторах нелинейный расчет сертифицирован а в ANSYS нет( я понимаю что это все ерунда но факт остается фактом). А вообще да можно взять пособие к СП 2003 года и посчитать армирование по допускаемым напряжениям но по закону мы и этого делать не можем так как вы можете ошибиться в алгоритме, а сертификация этот момент и проверяет. Но насколько я знаю скоро у нас СП и СНиП уйдут на второй план, все же чаша весов в последнее время склоняется к тому что ни Мельников ни Кучеренко ни Герсанов ни Промздания ни Гвоздев и др. не хотят разрабатывать нормативы , они хотят отдавать эту работу на плечи заинтересованных лиц и это далеко не проектировщики.

[MrWhite]

Сообщение от miko2009 MrWhite просто в калькуляторах нелинейный расчет сертифицирован а в ANSYS нет( я понимаю что это все ерунда но факт остается фактом). А вообще да можно взять пособие к СП 2003 года и посчитать армирование по допускаемым напряжениям но по закону мы и этого делать не можем так как вы можете ошибиться в алгоритме, а сертификация этот момент и проверяет. Но насколько я знаю скоро у нас СП и СНиП уйдут на второй план, все же чаша весов в последнее время склоняется к тому что ни Мельников ни Кучеренко ни Герсанов ни Промздания ни Гвоздев и др. не хотят разрабатывать нормативы , они хотят отдавать эту работу на плечи заинтересованных лиц и это далеко не проектировщики.

Ну если подходить дотошно, расчет по диаграммам (гипотеза плоского поворота) и расчет по теории предельного равновесия в одинаковой степени являются нелинейными, а точнее псевдонелинейными. Федеральный закон требует от нас гораздо большего.
Ручной расчет нормами не запрещен. Т.е. если весь алгоритм описан, показаны усилия, к-ты, формулы - проверить результаты не представляется сложным. В случае с гипотезой плоского поворота - это не только формулы, но и диаграммы напряженного состояния по сечению. В конечно счете ошибки возможны, но они так же возможны и при сертификации. Просто ответственность на разные плечи ложиться . Хотя опять же не знаю как юридически решаются вопросы связанные с ошибками при сертификации. Кто в конечном счете несет наказание? Авторы программы, органы сертификации или расчетчик?
Кстати сертифицированы не все действующие СП и СНиПы по расчетам железобетона. Например СНиП 2.06.08 и его дальнейшие вариации в виде СП. Не знаю ни одной программы которая прошла сертификацию по ним.

Ну а уж если говорить вообще про философию СНиП, то она конечно на мой взгляд устарела. Вполне достаточно ограничиваться общими положениями. Методы расчета лучше оставлять в ответственности исполнителя целиком и полностью.

[fasol]

Я под своим вопросом подразумевал следующее:
На территории РФ действуют посвященные железобетону СП 63 (ПГС) и СП 41 (ГТС), в которых однозначно определены все параметры расчета по теории предельных состояний.
У меня, допустим, есть конструкция, которую существующими методиками не осилить. И именно поэтому я расчехляю solid65 и concr (ручками пока что приятнее и понятнее).
Но куда сувать коэффициенты?
На коэффициент надежности по ответственности и на коэффициент сочетания нагрузок я могу домножить полезную нагрузку, окей. Тут вроде особых вопросов нет.
А коэффициент условия работ? Если он больше единицы, разумнее его проигнорировать, а если меньше? В сталежелезобетонных водоводах, например, я уменьшал расчетное сопротивление арматуры и бетона на 10%, потому что коэффициент по нормам 0,9. Насколько это адекватно?
А расчет по раскрытию трещин? Нагрузка должна быть нормативная, а расчетные сопротивления? Все по второй группе?
Непонятно.

----- добавлено через ~7 мин. -----
А сниповская диаграмма сжатия/растяжения бетона с выпуклостью ВНИЗ, которую ансис физически не может схавать? Зачем так?

[fasol]

Не могу понять, почему 10-я версия на i7+16gb считает solid65+conc в два раза быстрее, чем 14.5 и 15 на 2хXeon+128gb (которые частенько вылетают при этом)?
Элементы 181, 65, 8.

Настройки решателя
NROPT,UNSYM,,1
NLGEOM,1
NSUBST,20,1E5,20
NEQIT,7
OUTRES,ALL,ALL
CNVTOL,F,,0.5
CUTCON,NOSHAPE,1
CUTCON,PLSLIMIT,1

[derik]

а зачем настройки когда нужно лог смотреть.

[fasol]

А что конкретно там смотреть?

----- добавлено через ~2 мин. -----
И где этот лог найти?

[derik]

эм.....кхм....даже не знаю....в папке с рабочими файлами должен быть log файл. если что то завершается с ошибками он пишет почему эррор вышел, ну или после чего произошел фатал. тут и нужно смотреть и копать. настройки решателя тут влияют меньше, только если скорость счета
мб поменять нужно решатель. вы только несимметричным пользовались? другие не подставляли?

[fasol]

Ах, этот лог. Я думал, у решателя какой-то свой есть. Расчет абортится чаще всего из-за отсутствия сходимости. Бисекш, бисекшн, бисекшн и привет. Иногда dof limit вылетает. Что вообще очень странно, учитывая, что исполнительный макрос один и тот же.
Другой решатель попробую.

[MrWhite]

fasol это как раз очень возможно. Особенно если расчет прерывается. Это значит что в лог спамит какая то ошибка или варнинг. Спам отображения в окне этого события и задерживает расчет.

[fasol]

Конструкция: сталежелезобетонный напорный водовод. То есть металлическая облицовка, одетая в обойму из железобетона с двумя рядами кольцевых каркасов. Нагрузка - давление воды на облицовку, поэтому напряженное состояние соответствует внецентренному растяжению.
Неожиданно обнаружил, что изменение начального модуля упругости бетона в диапазоне от 0,3Е до 1,0Е (путем варьирования относительных деформаций, соответствующих 0,6Rb) вызывают изменение напряжений:
в облицовке и внутренней арматуре +/-60%
во внешней арматуре - в несколько раз.
Хотя по идее работа CONCR на растяжение билинейна - сначала монотонно до Rbt, потом полочка. Почему характеристики бетона на сжатие так сильно влияют?
Solid65, трехлинейная диаграмма по СП63.

[miko2009]

Сообщение от fasol Solid65, трехлинейная диаграмма по СП63.

так CONCR или трехлинейная диаграмма по СП63 ?
Offtop: Solid65 прошлый век

[MrWhite]

Сообщение от fasol Хотя по идее работа CONCR на растяжение билинейна - сначала монотонно до Rbt, потом полочка.

Никакой полочки там нет. Монотонно до Rbt, затем вертикально до 0.
И присоединюсь к вопросу miko2009: так 65-ый или трехлинейная диаграмма по СП63?

[Бахил]

Сообщение от MrWhite Монотонно до Rbt, затем вертикально до 0.

Да ладно. Ансиска так не умеет.

[fasol]

TB,CONC,2
TBTEMP,0
TBDATA,,0.3,0.8,RBTN,-1
TB,MISO,2,1,4
TBTEMP,0
TBPT,,ESB1,0.6*RBN
TBPT,,3.4E-3,RBN
TBPT,,4.8E-3,RBN


----- добавлено через ~10 мин. -----


----- добавлено через ~13 мин. -----
Если модуль на растяжение равен модулю на сжатие, тогда вопросов нет, потому что варьирование положения _b1 автоматически изменяет наклон в четвертом квадранте графика. Но так ли это? Мне казалось, бетон там живет совершенно обособленной от сжатия жизнью.

----- добавлено через ~27 мин. -----
То есть в моем понимании CONC (с параметрами трения в трещинах и Rbt) отвечает за зону растяжения, а MISO (с трехлинейкой по СП63) - за зону сжатия. Вроде всегда так делал, с натурой даже сходилось. Криминал штоле?

[MrWhite]

Сообщение от fasol То есть в моем понимании CONC (с параметрами трения в трещинах и Rbt) отвечает за зону растяжения, а MISO (с трехлинейкой по СП63) - за зону сжатия. Вроде всегда так делал, с натурой даже сходилось. Криминал штоле?

Ну вообще Вы условно конус с цилиндром пересекаете. Да это будет работать так как Вы хотите, но только до окружности их пересечения. В растяжении действительно CONC будет перебивать MISO, а вот в сжатии MISO будет перебивать CONC. И чем больше будет расти обжатие, тем сильнее будет врать MISO. Вплоть до полного бреда.

Вполне допускаю, что конкретно в Вашем расчете именно это и происходит.

[fasol]

Вот оно что...
А как правильно-то? Чтоб и трещины в 65м сохранить, и нелинейную диаграмму сжатия оставить?

[miko2009]

Сообщение от fasol Вроде всегда так делал, с натурой даже сходилось. Криминал штоле?

у вас ядро ANSIS при расчете никаких предупреждений не выдает ? По идее он должен вам сказать фразу на подобие "модель CONC отключена, в дальнейшем расчете применяется MISO" и вы считаете чистое MISO, трещинообразованием там и не пахнет.

Сообщение от fasol А как правильно-то? Чтоб и трещины в 65м сохранить, и нелинейную диаграмму сжатия оставить?

забыть про SOLID65 и использовать Menetrey-Willam с 18 версии ANSYS
https://repository.lib.ncsu.edu/bits...=1&isAllowed=y
http://cae-club.ru/publications/novy...mechanical-r18

[fasol]

Нет, 65й работает как надо, трещины раскрываются, закрываются и вообще красота (на первый взгляд).
Видать, действительно, конус с цилиндром пересекается и кто первее, того и тапки.

Сообщение от miko2009 18 версии ANSYS

Дальше 17й не обновлялись и вряд ли будем. А лично меня вообще 11я устраивает.

[miko2009]

Сообщение от fasol трещины раскрываются, закрываются

интересно как вы видите это ? трещинообразование происходит за счет выключение из работы просто элементов при достижении максимального сжатия/растяжения ,контролируется соответствующим KEYOPT (не помню какой именно), поэтому закрываться нечему, происходит эрозия конечных элементов.

[fasol]

Вижу как и все, либо через conrete plot 1-2-3-все трещины, либо через EPTO1, где картинка намного нагляднее.
При приложении нагрузки, трещины раскрываются, при снятии закрываются, кроме тех мест, где металл потек.

[MrWhite]

Сообщение от fasol Дальше 17й не обновлялись и вряд ли будем. А лично меня вообще 11я устраивает.

Ну с трещинами только Solid65 работает, остальные их только в виде нелинейных деформаций показывают. А вообще в 18.* реально хорошие материалы по бетону добавили. Мне кажется отказ от обновления до 18-ой версии не рационален, если занимаетесь бетоном и грунтами. Уж на что я консервативен, но и то сразу поставил
В 17-ой можно microplane попробовать, но лично мне его реализация в Ansys совсем не нравится.

[fasol]

Окей, а по-старинке 65м-то можно в принципе адекватно посчитать? Без MISO.

[MrWhite]

Solid65 нормально подходит для моделирования сооружения с проектной арматурой на проектные нагрузки. Особенно гидротехнических. Потому что там напряжения сжатия в бетоне редко когда превышают 0,65 от предела прочности. А если превышают то находят в состоянии трехосного сжатия, что сильно повышает предел прочности. Solid65 хорошо подходит под динамические задачи.
Solid65 плохо подходит под задачи подбора арматуры (если на основании НДС пытаться определить нужную арматуру). Задач разрушения - типа прогрессирующего. Ну и вообще любых других задач на которых нужно замоделировать уже серьезно текущую арматуру. Так как там факторы нелинейного деформирования бетона становятся ключевыми.

Кстати диаграмма которую вы пытаетесь применить описывает одноосное НДС поэтому худо бедно подходит под конструктивные элементы типа колон, балок и пр. хуже для элементов типа плит, и совсем плохо для сложных по форме объемных тел. И по своей сути это расчетный инструмент, а не инструмент для моделирования.

[fasol]

Как бы то ни было, CONC у меня лучше всего срабатывает в связке с MISO. В одиночку или с KINH получается нефизичная ерунда с раскрытием трещин там, где их совершенно не должно быть. И эту комбинацию CONC+MISO я не сам придумал, а у Волански подсмотрел. Да, у него там была однопролетная балка.

[miko2009]

Сообщение от fasol В одиночку или с KINH получается нефизичная ерунда с раскрытием трещин там, где их совершенно не должно быть.

KINH + CONC не работают, это же совершенно разные мат. теории , почитайте файл "solution information", вам расчетное ядро заменит solid65 на solid185\186 (в зависимости от формы КЭ) и соответственно никакого CONC не будет, будет просто обычный изотропный материал

[MrWhite]

На однопролетной балке даже материал чугуна с забитыми характеристиками бетона себя вести будет очень похоже на бетон. Балка не показатель совсем..

[fasol]

Сообщение от miko2009 KINH + CONC не работают, это же совершенно разные мат. теории , почитайте файл "solution information", вам расчетное ядро заменит solid65 на solid185\186 (в зависимости от формы КЭ) и соответственно никакого CONC не будет, будет просто обычный изотропный материал

У нас с вами, видимо, какие-то совсем разные ансисы, потому у меня тип элемента не меняется, трещины рисуются, все как надо.
К тому же в ломалке бетоно у MrWhite KINH есть:

[miko2009]

Сообщение от fasol У нас с вами, видимо, какие-то совсем разные ансисы

это не ANSYS разный, это разные математические модели, и я 100% уверен в этом, да и не я один

Сообщение от fasol К тому же в ломалке бетоно у MrWhite KINH есть:

KINH это закон деформирования и если я не ошибаюсь то в ломалке MrWhite все же заложил их порознь, или все же у MrWhite "вера" пошатнулась ?

[MrWhite]

Сообщение от miko2009 KINH это закон деформирования и если я не ошибаюсь то в ломалке MrWhite все же заложил их порознь, или все же у MrWhite "вера" пошатнулась ?

Практика критерий истины жеж . Я проверил (вдруг я чего то недопонимал) и сделал правильные выводы.