|
||
| Правила | Регистрация | Пользователи | Поиск | Сообщения за день | Все разделы прочитаны | Справка по форуму | Файлообменник | |
|
Поиск в этой теме |
|
||||
немного топик аппа
смоделировал балку из диссертации Зиганшина Х.А.. Балка под номером Б-II-2а, все характеристики материалов, закрепления взяты из данной работы. В аттаче картинка трещинообразования по опытным данным и смоделированная в ансис. |
||||
|
||||
Регистрация: 18.10.2012
Сообщений: 5
|
Здравствуйте! 6ра4о
Интересно, а какие характеристики материалов, закрепления из диссертации Зиганшина, и какой тип элементов Вы использовали при моделировании в ансис, если не секрет? Меня интересует проблема моделированния железжбетонных конструкций в анснс, однако у меня нет доступа к российской базе диссертаций. С интересом следил за этой темой. Кстати, я иностранец, но это неважно. Заранее спасибо. |
|||
|
||||
материалы из диссера кроме модуля упругости бетона и арматуры, элемент в ансисе 65 для бетона и 264 для арматуры;concr модель материала бетона с кинематическим упрочнением, арматура - просто кинематическое упрочнение . данный диссер есть на сайте жбк.рф в разделе библиотека, есть не сможете скачать напишите, "намылю".
|
||||
|
||||
Регистрация: 23.04.2013
Сообщений: 2
|
Подскажите, что я делаю не так:
/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 |
|||
|
||||
в месте соединения ребра и пластины у Вас несогласованная сетка - Вам необходимо пересечь эти два тела. Сообщение говорит что получаются перемещения узла такого то больше чем допустимо. ребро удлините до наружной плоскости пластины и введите vptn,all.
|
||||
|
||||
Регистрация: 23.04.2013
Сообщений: 2
|
Спасибо, с 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 |
|||
|
||||
Dr Регистрация: 16.10.2011
Днепропетровск
Сообщений: 1,020
|
Доброго времени суток, только начинаю разбираться с solid65 и бетоном в ansys, и у меня такой вопрос возник, при расчете, какие необходимо вводить данные: кубиковая прочность или призменная? Насколько сильно влияет ввод таких данных как прочность на двухосное сжатие и прочность при гидростатическом давлении, и как следствие где можно найти эти данные о материале?
|
|||
|
||||
расчеты МКЭ и CFD. ктн Регистрация: 11.05.2005
Подмосковье
Сообщений: 2,170
|
советов универсальных дать трудно.
если собираетесь сравнивать со снип (чтобы показывать практикам) нужно брать из снип два предела, на растяжение и на сжатие. остальное по умолчанию. (там есть значения- разберитесь в логике их назначения) у меня с арбатом отличия получаются не более 5%. далее можно считать то, что арбату не по силам. если сравниваете с экспериментом- надо подбирать характеристики и параметры, чтобы совпало с замерами. начиная с тех, что назначены по итогам эксперимента. все дополнительные параметры лично я толком не варьировал. Думаю, их влияние более важно для управления сходимостью в существенно трехосных задачах (вроде плотин). надеюсь, через месяц-два я смогу добавить информации. |
|||
|
||||
Юмористические расчеты, комедийные диссертаций, цирковые статьи Регистрация: 12.12.2012
Москва
Сообщений: 892
|
Сергей Иванович, вопрос к Вам есть небольшой. Давили ли Вы кубики, цилиндры и призмы в ANSYS, чтобы на выходе иметь при одинаковых свойствах материалов идентичную кубиковую, цилиндрическую и призменную прочность? И если давили, то какие результаты?
|
|||
|
||||
Dr Регистрация: 16.10.2011
Днепропетровск
Сообщений: 1,020
|
СергейД, спасибо за ответ, и если можно еще такой вопрос:
В ansys необходимо вводить истинные напряжения - деформации в качестве исходных данных, для бетонного образца во всех виденных мной источниках присутствует ниспадающая ветвь, как я понял это связано с тем, что происходит разрушения сечения образца (образование пирамид), а следовательно диаграммы в нормах представлены в инженерных напряжениях деформациях, переход от которых к истинным через логарифм не сделать (как сделать этот переход если честно, в литературе пока не нашел), отсюда вопрос: Следует ли задавать материалу ниспадающую ветвь, или лучше заканчивать прямым участком? Проверил оба варианта на задаче по раздавливанию кубика (с ниспадающей ветвью и без). С ниспадающей ветвью задача решилась лучше ( эквивалентные относительные пластических деформации + упругие относительные деформации к моменту обрыва решения(оборвал сам, так как пошло много кратное деление шага безрезультатное) соответствовали исходной диаграмме. Без ниспадающей ветви задача разошлась по достижению предельных напряжений. Последний раз редактировалось Stanum, 18.11.2013 в 20:28. |
|||
|
||||
Юмористические расчеты, комедийные диссертаций, цирковые статьи Регистрация: 12.12.2012
Москва
Сообщений: 892
|
Сергей Иванович, предлагаю для вас и коллег следующее описание верификационного теста. Интересно насколько модель бетона для 65-х сможет его выполнить. Давятся образцы со следующими параметрами:
- цилиндр 150х300 мм, - куб - 150х150х150 мм, - призма - 150х150х400 мм. Параметры материалов одинаковые для всех вариантов. Варианта испытанию как я понимаю возможно два. 1-й вариант - крепление нижнего торца по XYZ, верхнего по XY, но с приложением перемещений эквивалентных предельной деформации сжатия. По суммарной опорной реакции по Z определение графика деформирования с определением прочности на сжатие (кубической, призменной и цилиндрической, которые связаны известными зависимостями). 2-й вариант - контактная задача с теми же параметрами. |
|||
|
Опции темы | Поиск в этой теме |
|
|
Похожие темы | ||||
Тема | Автор | Раздел | Ответов | Последнее сообщение |
Нужна помощь с расчетом ж/б "клюшки" по серии 1.822.1-2/82 вып.4 | Lexzus | Железобетонные конструкции | 4 | 11.12.2011 11:10 |
Помогите с расчетом ж/б стены прямоугольного в плане бассейна | Den_13 | Железобетонные конструкции | 2 | 18.01.2010 10:41 |
Скатная кровля для ж/б каркаса ИИ-04 | Antony206 | Железобетонные конструкции | 6 | 14.12.2009 23:18 |
Наращивание существующего ж/б каркаса | Белозерова Татьяна | Железобетонные конструкции | 17 | 16.02.2009 14:23 |
узлы сопряжения ж/б ригелей и стальных стоек каркаса | ТОЛЬКО ПРИШЕДШИЙ | Поиск литературы, чертежей, моделей и прочих материалов | 0 | 10.12.2007 12:02 |