Ansys для проектировщика КЖ-КМ

[Fedor 1 597]

"Ни одна не открывается!" - попросите Валерия Московского он поможет открыть

[ValeryMoscow 3]

"Ни одна не открывается!" - попросите Валерия Московского он поможет открыть

статьи можно посмотреть на нашем сайте

cae-services

в разделе ГАЛЕРЕЯ => Наши расчеты

[Chardash 306]

статьи можно посмотреть на нашем сайте

cae-services

в разделе ГАЛЕРЕЯ => Наши расчеты

спасибо, посмотрел. еще бы увидеть сравнение результатов расчета ветра в ансис (или тп) со СНиП или Еврокод. Изменено 16 мая 2011 пользователем Chardash

[Chardash 306]

Здравствуйте, считаю безбалочное перекрытие в ansys . Подскажите пожалуйста аналог АЖТ лира-скад-робот? Необходимо для моделирования капители в оболочечной модели.

[sergeyd 3]

насчет аналога точного не знаю, но капитель моделируется с эксцентриситетом

берете 181 элемент. задаете ему сечение (все равно, начиная с 13.0, real не применяется)

капители одно, перекрытию другое.

изучаете команду secoffset

проверяете визуально /eshap,1

[Chardash 306]

насчет аналога точного не знаю, но капитель моделируется с эксцентриситетом

берете 181 элемент. задаете ему сечение (все равно, начиная с 13.0, real не применяется)

капители одно, перекрытию другое.

изучаете команду secoffset

проверяете визуально /eshap,1

спасибо!

[Fedor 1 597]

Наверное лучше 281. Хотя в чем его отличие от 93 для однородного не понял. В задачах ведут себя одинаково. Расхождение по напряжениям 0.33 %

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

[Chardash 306]

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

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

В Лира с эти проще, выбираем КЭ или вообще, делаем экспорт из Лира Арм. Как тут обстоит с этим вопрос, пока не ясно.

Или, первое, что приходит в голову, отдельно моделировать эти зоны элементами.

Изменено 20 июня 2011 пользователем Chardash

[Fedor 1 597]

Я определяю нормальные силы и моменты и подбираю арматуру в Mathematica. Все равно больше двух, максимум 3 видов армирования на чертеже не делается. Обычно делю на области с напряжениями больше определенного уровня и меньше, Потом вспавляю изображение в автокад масштабирую и поверх картинки ставлю арматуру, чтобы охватывала необходимые области.

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

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

[Chardash 306]

Я определяю нормальные силы и моменты и подбираю арматуру в Mathematica. Все равно больше двух, максимум 3 видов армирования на чертеже не делается. Обычно делю на области с напряжениями больше определенного уровня и меньше, Потом вспавляю изображение в автокад масштабирую и поверх картинки ставлю арматуру, чтобы охватывала необходимые области.

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

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

спасибо, Федор! попробую найти эти сообщения. Жалко, что в WB это все не реализованно. до этого просто понижал модуль для прогибов.в лира, и в солидах/оболочках ансиса. но, пока что солидами и напряжениями мой расчет в ансис заканчивается. все же хочется увидеть картину, как можно более похожую на действительность. Особенно после расчетов балок с большими сечениями по СП, когда армирование получается больше СНиПовского в 2 раза.

"А как это оценивать с учетом ползучести тоже обсуждалось где-то здесь"

в лира это удобно, просто пару цифр в установке нелинейного расчета

Изменено 21 июня 2011 пользователем Chardash

[Fedor 1 597]

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

В Ansys много разных нелинейностей. Пишешь код если понимаешь что хочешь, дальше контрал це и контрал вэ

[Chardash 306]

Пишешь код если понимаешь что хочешь

постараюсь разобраться. Изменено 21 июня 2011 пользователем Chardash

[Fedor 1 597]

Это в Ansys самое главное - понимание механики. Остальное мелочи.

[Chardash 306]

согласен. кнопки жать и потом смотреть картинку - не самое сложное. но и это в ansys как полоса препятствий. в abaqus интерфейс подружелюбней

[Fedor 1 597]

Что может быть проще нормального языка программирования? :)

[KostyaPC 0]

сильно не бить. что-то когда-то делал и бросил. в каком состоянии сейчас - не помню

FINISH



/CLEAR



/CONFIG,NRES,200000



/FILNAME,rc



! Redirect text output

! 0 - to the screen 

! 1 - to a file

CreateLog = 0



*IF,CreateLog,EQ,0,THEN

	/OUTPUT

*ELSE

	/OUTPUT,log,txt

*ENDIF



CreateReport = 0



E_c = 3e10

nu_c = 0.2

R_c_c_u = 1.45e7

R_c_t_u = 1.05e6



! Strain-stress diagram for concrete

! 1 - SP 52-101 bilinear

! 2 - SP 52-101 trilinear

! 3 - EN 1992-1-1 bilinear

! 4 - EN 1992-1-1 parabolic-rectangular



ConcreteDiagram = 4



E_s = 2e11

nu_s = 0.3

R_s_t_u = 4.35e8



! Strain-stress diagram for prestressing steel

! 1 - SP 52-102 trilinear

! 2 - EN 1992-1-1 trilinear (not implemented)



PrestressingSteelDiagram = 1



! Strain-stress diagram for steel

! 1 - SP 52-101 bilinear (not implemented)

! 2 - EN 1992-1-1 bilinear (not implemented)



SteelDiagram = 1



L = 4

b = 0.3

h = 0.5



d_s = 2e-2



q = 5.05e4



epsilon_c1 = 1.5e-3

epsilon_c_u1 = 3.5e-3



epsilon_c2 = 2e-3

epsilon_c_u2 = 3.5e-3



epsilon_c3 = 1.75e-3

epsilon_c_u3 = 3.5e-3



epsilon_c4 = 2e-3

epsilon_c_u4 = 3.5e-3



epsilon_s_u1 = 2.5e-2



Concrete = 1

Steel = 2

PrestressingSteel = 3



/UNITS,SI



/PREP7



	PSCONTROL,ALL,ON



	ET,1,SOLID65

	ET,2,LINK180



	!KEYOPT,1,3,2

	!KEYOPT,1,5,1

	!KEYOPT,1,6,3

	KEYOPT,1,7,1



	R,2,ACOS(-1)*d_s**2/4



	MP,NUXY,Concrete,nu_c

	

	*IF,ConcreteDiagram,EQ,1,THEN

		MP,EX,Concrete,R_c_c_u/epsilon_c1

		TB,MISO,Concrete,1,2

		TBPT,,epsilon_c1,R_c_c_u

		TBPT,,epsilon_c_u1,R_c_c_u

	*ENDIF



	*IF,ConcreteDiagram,EQ,2,THEN

		MP,EX,Concrete,E_c

		TB,MISO,Concrete,1,3

		TBPT,,0.6*R_c_c_u/E_c,0.6*R_c_c_u

		TBPT,,epsilon_c2,R_c_c_u

		TBPT,,epsilon_c_u2,R_c_c_u

	*ENDIF



	*IF,ConcreteDiagram,EQ,3,THEN

		MP,EX,Concrete,R_c_c_u/epsilon_c3

		TB,MISO,Concrete,1,2

		TBPT,,epsilon_c3,R_c_c_u

		TBPT,,epsilon_c_u3,R_c_c_u

	*ENDIF



	! EN 1992-1-1 parabolic-rectangular

	*IF,ConcreteDiagram,EQ,4,THEN

		MP,EX,Concrete,R_c_c_u*(1-(1-0.1)**2)/(0.1*epsilon_c4)

		TB,MISO,Concrete,1,11

		TBPT,,0.1*epsilon_c4,R_c_c_u*(1-(1-0.1)**2)

		TBPT,,0.2*epsilon_c4,R_c_c_u*(1-(1-0.2)**2)

		TBPT,,0.3*epsilon_c4,R_c_c_u*(1-(1-0.3)**2)

		TBPT,,0.4*epsilon_c4,R_c_c_u*(1-(1-0.4)**2)

		TBPT,,0.5*epsilon_c4,R_c_c_u*(1-(1-0.5)**2)

		TBPT,,0.6*epsilon_c4,R_c_c_u*(1-(1-0.6)**2)

		TBPT,,0.7*epsilon_c4,R_c_c_u*(1-(1-0.7)**2)

		TBPT,,0.8*epsilon_c4,R_c_c_u*(1-(1-0.8)**2)

		TBPT,,0.9*epsilon_c4,R_c_c_u*(1-(1-0.9)**2)

		TBPT,,epsilon_c4,R_c_c_u

		TBPT,,epsilon_c_u4,R_c_c_u

	*ENDIF



	!TBLIST

	TBPLOT

	/AXLAB,X,STRAIN

	/AXLAB,Y,STRESS

	/REPLOT

	*IF,CreateLog,EQ,1,THEN

		/SHOW,PNG,,0

		/GFILE,1600

		/REPLOT

		/SHOW,CLOSE

	*ENDIF



	TB,CONCR,1

	TBDATA,1,0.4,0.9,R_c_t_u,R_c_c_u



	MP,NUXY,Steel,nu_s



	*IF,SteelDiagram,EQ,1,THEN

		MP,EX,Steel,2e11

		TB,MISO,Steel,1,2

		TBPT,,R_s_t_u/E_s,R_s_t_u

		TBPT,,epsilon_s_u1,R_s_t_u

	*ENDIF



	*IF,SteelDiagram,EQ,2,THEN

		MP,EX,Steel,2e11

		TB,MISO,Steel,1,2

		TBPT,,R_s_t_u/E_s,R_s_t_u

		TBPT,,epsilon_s_u2,R_s_t_u

	*ENDIF



	!TBLIST

	TBPLOT

	/AXLAB,X,STRAIN

	/AXLAB,Y,STRESS

	/REPLOT

	*IF,CreateLog,EQ,1,THEN

		/SHOW,PNG,,0

		/GFILE,1600

		/REPLOT

		/SHOW,CLOSE

	*ENDIF



	/VIEW,,1,1,1

	!/REPLOT



	K,1,b/2-0.05,h-0.05

	K,2,b/2-0.05,h-0.05,L

	L,1,2



	BLC5,b/4,h/2,b/2,h,L



	MSHKEY,1

	MSHAPE,0,3D

	ESIZE,0.05



	LSEL,S,LENGTH,,L

	LESIZE,ALL,0.2

	LSEL,ALL



	TYPE,1

	MAT,Concrete



	VSWEEP,ALL



	TYPE,2

	MAT,2

	REAL,2



	LMESH,1



	NSEL,ALL

	NUMMRG,NODE



	/ESHAPE,1



	ESEL,,ENAME,,SOLID65

	/TRLCY,ELEM,0.75

	ALLSEL



	/REPLOT



FINISH



/SOLUTION



	ASEL,,LOC,Z,0

	DA,ALL,ALL

	ALLSEL



	ASEL,,LOC,Y,h

	SFA,ALL,1,PRES,q

	ALLSEL



	NSEL,,LOC,X,0

	DSYM,,X

	ALLSEL



	/REPLOT



	ANTYPE,STATIC



	!NLGEOM,ON

	NROPT,FULL,,ON

	!STAOPT,VT

	NSUBST,50000,100000,1000

	NEQIT,500

	!NLDIAG,NRRE,ON

	!NLDIAG,EFLG,ON

	AUTOTS,ON



	!RESCONTROL,,ALL,1



	OUTRES,ERASE

	OUTRES,ALL,ALL



	!CNvTOL,F,,0.05

	!CNvTOL,U,,0.05



	EQSLV,PCG,1e-5

	PCGOPT,,,YES

	!BCSOPTION,,,,,,PERFORMANCE



	*GET,StartTime,ACTIVE,0,TIME,WALL



	SOLVE



	*GET,StopTime,ACTIVE,0,TIME,WALL

	SolutionTime=(StopTime-StartTime)*3600

	!*MSG,INFO,SolutionTime

	!Solution time: %G seconds



FINISH



/POST1



	SET,LAST

	/DEVICE,VECTOR,1

	PLCRACK,0,1



	!/EXPAND,1,RECT,HALF,1e-8

	!/AUTO



	!ETABLE,IFORCE,LEPPL,1

	!ETABLE,JFORCE,LEPPL,2



	!ESEL,S,ENAME,,LINK180



	!PLLS,IFORCE,JFORCE



	! Turn contour values on

	!/PNUM,SVAL,ON



	PATH,SIGMA,2

	NSEL,,LOC,X,0

	NSEL,R,LOC,Y,0

	NSEL,R,LOC,Z,0

	*GET,NodeNumber,NODE,0,NUM,MIN

	ALLSEL

	PPATH,1,NodeNumber

	NSEL,,LOC,X,0

	NSEL,R,LOC,Y,h

	NSEL,R,LOC,Z,0

	*GET,NodeNumber,NODE,0,NUM,MIN

	ALLSEL

	PPATH,2,NodeNumber

	PDEF,SIGMA,S,Z

	PLPAGM,SIGMA

	/AUTO



!FINISH



!/POST26



!	NSEL,,LOC,X,0

!	NSEL,R,LOC,Y,0

!	NSEL,R,LOC,Z,L

!	*GET,NodeNumber,NODE,0,NUM,MIN

!	NSOL,2,NodeNumber,U,Y,UY

!	STORE,MERGE

!	!XVAR,1

!	PLVAR,2

!	*IF,CreateLog,EQ,1,THEN

!		/SHOW,PNG,,0

!		/GFILE,1600

!		/REPLOT

!		/SHOW,CLOSE

!	*ENDIF

[pavlik 1]

сильно не бить. что-то когда-то делал и бросил. в каком состоянии сейчас - не помню

Спасибо за интересный материал!

Можно только один вопрос? Всегда было интересно зачем в код вставляют закоментированные команды - вот как у вас:

!TBLIST

Ведь они все равно не выполняются?

Изменено 21 июня 2011 пользователем pavlik

[Fedor 1 597]

Экспериментировал человек с разными вариантами. Чтобы с ифами не связываться поотключал :)

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

[Chardash 306]

"сильно не бить. что-то когда-то делал и бросил. в каком состоянии сейчас - не помню"

Это ничего. Спасибо, конечно. В первом приближении не пошло, а дальше посмотрим.

[KostyaPC 0]

Всегда было интересно зачем в код вставляют закоментированные команды - вот как у вас:

!TBLIST

Ведь они все равно не выполняются?

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

В первом приближении не пошло, а дальше посмотрим.

что именно не пошло?