CC BY
33
Теория расчета строительных конструкций
УДК 624.048
СПОСОБ ФОРМИРОВАНИЯ НЕЛИНЕЙНЫХ ФИЗИЧЕСКИХ СООТНОШЕНИЙ В ПРЯМЫХ И ОБРАТНЫХ ЗАДАЧАХ РАСЧЕТА МНОГОФАЗНЫХ СТЕРЖНЕЙ
А.В. Мищенко
Предложен алгоритмичный способ формирования физических соотношений для составных многофазных стержней, основанный на аппроксимации диаграмм деформирования фазовых материалов целыми рациональными полиномами произвольной степени. Применены жесткостные характеристики высших порядков.
Ключевые слова: многофазный стержень, армированный стержень, слоистая структура, физическая нелинейность, жесткостные характеристики высших порядков, предельные состояния.
Рассматриваются стержни, составленные из квазиоднородных частей (фаз) произвольной геометрической формы. При условии рационального проектирования подобные структуры позволяют создавать весьма эффективные несущие конструкции, адаптированные к возникающему в них напряженно-деформированному состоянию, обладающими повышенной несущей способностью, жесткостью, трещиностойкостью при одновременном снижении расхода материалов [1]. Номенклатура используемых в них материалов может включать металлы, бетоны, конструкционные пластмассы, древесину, искусственные дисперсные композиты и т.п.
1. В системе координат xyz многофазный стержень (МФС) обладает свойствами симметрии относительно плоскости ху, в которой стержень испытывает прямой продольно-поперечный изгиб, продольная ось совмещена с осью х. Будем считать, что в стержне: обеспечены условия идеального межфазного контакта при отсутствии разрывов в деформациях на границах фаз; материал каждой фазы представляет собой сплошную, однородную, изотропную нелинейно упругую среду со значениями физических характеристик одинакового порядка; длина и характерный поперечный размер стержня удовлетворяют условию I > (5 ^ 6) h .
В зависимости от геометрических параметров фаз может быть выделен ряд частных типов структур МФС (рис. 1): поперечно- и продольно-слоистая, полярно- и радиально-слоистая, поперечно-продольно-слоистые структуры. Отличаются они формами межфазных поверхностей и направлением реализации неоднородности. Так, для схем «б», «в», «г», «д» такими направлениями соответственно являются: ось у, луч ст, радиус г, ось х.
Примем вариант кинематических соотношений и условий равновесия в форме
еx +1 и 1, sin8 »8, cos8» 1, ех(xy) = 8о-к-y , Sy(xy) = 0
Уyx(X У) = -Уо(x) = v'-8 , е0(x) = u', к(х) = 8',
(1) (2)
(3)
N' + (бе)' = -?х, Q'- (N8)' = ду,
М' = Q + т2.
Здесь () ' = д() / дх, уух, ех, еу - сдвиг, продольная и поперечная деформация; и, V - продольное и поперечное перемещение точек отсчет-ной оси х; 8 - угол поворота поперечного сечения;
[N,Q,М](х) = £Ц^,т£>,-ахк)у]М - (4)
к=! Ак
продольная и поперечная силы, изгибающий момент.
Основой расчетной модели МФС является способ представления физических соотношений (определяющих соотношений). Из трех компонент тензора напряжения в данной задаче при учете принятых допущений основной компонентой будем считать стх, а тух , сту - второстепенными [2].
Последней из них в данной задаче будем пренебрегать. Для основной компоненты будут выполнены все группы соотношений, а для второстепенных - все, кроме физических, что является результатом введения в расчетную модель гипотез (1).
Компоненту ст^) тензора напряжений для материала к-й фазы представим в виде суммы
(5)
стхк )(е) = 1 р/ (е), 1=0
где /е1 (1 = 0,1,..., г) - набор гладких базисных
функций, определенных либо на всей области допустимых деформаций, либо на непересекающихся ее подобластях, что дает возможность описывать особенности законов деформирования в виде
Мищенко А.В.
Способ формирования нелинейных физических соотношений в прямых и обратных задачах расчета многофазных стержней
д)
б)
1 С
й
Ь1(х,у)
Л
У1
в)
- у
Гк
к\П/-
1д\\ф
И/ ¡1 дт\«1 1
(к кс
-¿г) z
- у
Гк
>д\ ф
11/1 /лк\п
.у !! z
X
Рис. 1. Структуры многофазных стержней: общего вида (а), поперечно-слоистая (б), полярно-слоистая (в), радиально-слоистая (г), продольно-слоистая (д), продольно-поперечно-слоистая (е)
разрывов функции напряжения и ее производных; Рк = [Рк,+, ра,7 ] - числовые параметры 1-й аппроксимации для областей растяжения и сжатия в к-м материале.
В данной статье рассмотрим способ построения физических соотношений, основанный на использовании в качестве базисных функций степеней с натуральными показателями [3]
/.дв) = в1, (1 = 0,1,...,г). (6)
В отличие от множества других вариантов представления законов деформирования способ (5), (6) отличается: наличием физического смысла и совпадением с классическими методами в случае применения первых двух членов ряда (6), дающих модели жестко-пластического и линейно-упругого деформирования [4]; возможностью последовательного уточнения при увеличении членов разложения [5].
Подставив (5), (6), (1) в первые два выражения (4), получим
£ Вт (Во ,к) = N , £ ВМ1 (Во ,к) = М . (7)
1=0 1=0
Здесь Вм, ВМ1 - полные алгебраические полиномы 1-й степени обобщенных деформаций вп, к :
В№ = £ CijDijВ0 7=0
(8)
ВМ/ = £ СуГ17+1В0 7 к^ , 7=0
записанные с использованием жесткостных характеристик Г)^ многофазного сечения из
нелинейно-упругих материалов 1 -го физического и 7 -го геометрического порядков
Г (х)=£ рк// л у^А=£ Г
(к) 7
(9)
к=1
к=1
Су (7 = 0,..., 1) - коэффициенты биномов
Ньютона 1 -й степени (в0 7 к)1. Вычисляя
жесткости (9), следует использовать параметры
материалов р^ = [ рк,+, рк,7 ], соответствующие
знакам деформации.
Для степеней 1 = 0, 1, 2, 3 левые части уравнений (7) принимают вид
|вт = Г10В0 7Г11к ,
IВМ1 = 7Г11В0 + Г12 к ,
IBN2 = Г20В2 7 2Г21В0к + Г22к , ВМ2 =7Г21В2 + 2Г22В0к 7 Г23к ,
I BN0 = Г00,
IВМ 0 = 7 Г01/
I^3 = А0В3 7 3Г31В0к + 3Г32В0к2 7 Гззк3,
1ВМ3 = 7Г31В0 + 3Г32В2к 7 3Г33В0к2 + Г34к3.
Для конкретных неоднородных структур в зависимости от свойств фазовых материалов могут применяться различные (необязательно последовательные) наборы блоков (8). Физическое соотношение для сдвигов может быть представлено в виде у0DQ = Q с использованием сдвиговой жесткости многофазного сечения DQ.
2. Рассмотрим постановки обратных задач на основе соотношений (8).
Для замыкания краевой задачи (2), (3), (7) при соответствующих граничных условиях используем ограничение деформации
|в х (х, у) <в®т =вХк)/ К, (к = 1,..., * ), (10)
где вХ^ - допустимый уровень деформации при растяжении и сжатии при оценке предельных состояний материала к-й фазы; кВ> 1 - коэффици-
z
к
Теория расчета строительных конструкций
ент запаса по отношению к уровню предельной деформации е® .
Способ записи расчетного критерия через деформации обладает определенными преимуществами по отношению к способу записи через напряжения. Среди них: линейность критерия относительно величин обобщенных деформаций; отсутствие необходимости использования физических соотношений при формулировке левой части критерия. Применение же в расчетном критерии разрывных нелинейных функций напряжений существенно усложняет анализ возникающих предельных состояний.
В правой части критерия (10) могут использоваться предельные деформации (как таковые), превышение которых недопустимо по тем или иным соображениям. Кроме того, предельные деформации могут определяться на основе специальных зависимостей, описывающих некоторые физические предельные состояния. Например, предельный уровень деформации, вводимой по условию прочности, потребует формулировки обобщенного эквивалентного напряжения ст®, = ^(стХк), тУ?) <ст(1к) согласно принятой теории прочности и, как следствие, - возникновение зависимости правой части (10) от различных компонент тензора напряжения [4]. Таким образом, несмотря на принятую для удобства квазиодноосную форму записи, критерий (10) отражает и сложные напряженные состояния.
В обратных задачах рационального проектирования, связанных с поиском эффективной структуры стержня, потребуем реализации условия (10) в форме строгого равенства в максимально возможном числе точек пространства X, Y, занятого стержнем. В 5 -фазном стержне условие предельного состояния сформулируем в виде непрерывно-дискретного многоточечного расчетного критерия (НДРК) на основе следующего определения.
Предельное состояние на участке многофазного стержня, занимающего пространство X, Y, при нелинейно-упругом деформировании реализуется в том случае, если условие (10) выполняется в форме строгого равенства
|е0(х) - к( х) у\ = 8®+ (11)
непрерывно по продольной координате х е X* е X в дискретном числе уровней пу по
поперечной координате у = = Y* е Y. Для
всех иных точек х е X*, у г Y*. условие (10) реализуется в форме неравенства.
Кратность критерия пу, в зависимости от условий задачи, принимает значения от 0 до 5 +1. В общем случае удовлетворение указанных равенств кроме варьирования геометрических функций требует и изменения физической структуры стержня. Для стержней с фиксированной структурой максимальная кратность .. равна двум.
В поперечно-слоистых структурах выявление серии точек |1,...,удобно выполнять с использованием полигона предельных деформаций „¡¿Ш (У). На рис. 2 показано предельное состояние для двухточечного НДРК при осуществлении внутреннего касания эпюрой деформаций е*(у) полигона допустимых деформаций в двух точках. Это позволяет сразу определить расчетную деформацию и кривизну
+Р(к1)+£ + Р(к2)+е -ьаат Ъ2 + ьаат Ъ1
80* = "
к*(х)=
^2 -^1
+ №)+ + „(¿2)+ —ьааш + ьаат
^2 -^1
(12)
для дальнейшего использования при решении краевой задачи (2), (3), (7). Здесь при растяжении берутся верхние, а при сжатии - нижние знаки.
5 с
у
8*( У)
В^) Г
8^- [
( 5) + адш
] 8(х)+
Рис. 2. Геометрическая интерпретация расчетного критерия (11) для поперечно-слоистой структуры при пу = 2
В сечении многофазной структуры может существовать несколько допустимых вариантов расчетных полей деформаций, отличающихся координатами точек касания допустимого полигона линией е*(у). Для выбора наиболее рационального варианта может использоваться оценка
1 5 1
К8 = - Т , (к)+ („0Ак ШaX, (13)
5 к=1 Ак8ааш
основанная на суммарной интегральной величине 0 < К8 < 1 относительного уровня деформации. Применение (13) дает возможность: а) для заданных неоднородных структур осуществить отбор наилучших деформированных состояний 8* ( у) ; б) осуществить отбор физических структур, наиболее эффективных для деформированных состояний определенных типов. Для структур с прямоугольными сечениями слоев выражение (13) принимает вид
К8 = ^ Т -?1
' к=1 8,
(к )+ а^
2 2 Л . У2+1 - Ук
' 2hk
^ шах.
1
г
8
0
Мищенко А.В.
Способ формирования нелинейных физических соотношений в прямых и обратных задачах расчета многофазных стержней
Жесткостные характеристики записываются для конкретной структуры при помощи выражений (9). Например, для структуры, показанной на рис. 2, получим выражения
Dj (x) = ¿ pkb (yj+i - yj+1) / (j +1), (14)
k=i
которые позволяют организовать численную или аналитическую процедуру нахождения ширины bki (x), bj2 (x) или толщины hji (x), hj2 (x) расчетных (варьируемых) слоев.
Рассмотрим стержень, армированный продольными волокнами (рис. 3). Выделим в нем две фазы: арматуру (k = 1), деформирующуюся линейно ctx1^ = Ea6, и связующее (k = 2) при квадратичном законе деформирования
42)+ = P2,+ + Р2> + Р2>2 , ^ = Plfi + P2> + РУ2 .
Систему уравнений (7) сформируем из блоков (8) 0, 1 и 2 степени
+D00 + AoSo - D11К + D206(2 - 2D2160K + D22K = N,
-/3,160 + D2К - d21sg + 2d22S0K -d23K2 = m. (15) 1 1 " 2 ' Жесткостные характеристики (9) примут вид
D00( x) = p2 ,0 A2 + p2 ,0A2 ,
D01(x) = p2,0S2 + p2 ,0S2 ,
D10 (x) = EA + p£A2+ + Р2,- A2- ,
Dn(x) = EAya + ^2+ + Р2,-^2 ,
D12 (x) = EAy2a + P2,+12+ + Р2,-/2 ,
D20( x) = p2 ,2A2 + p2 ,2A2 ,
D21(x) = p2 ,2 S2 + p2 ,2 S2 ,
D22 (x) = P2+2I+ + P2,2I2 ,
D23( x) = p2 ,2 J2 + p2 ,2J2 ,
где A1 = Aaw - площадь сечения w армирующих
элементов; A2, S2, /2 , J - геометрические характеристики областей растяжения и сжатия вто-
рой фазы - обобщенные моменты инерции соответственно 0, 1, 2 и 3-го порядков, в соответствии с (9) при к > 0 заданные выражениями
У0( х)
[А+,^,/2+,J2+ ](х) = | ¿2(х)[1,у,У2,У3Му ,
Л( х)
У2( х)
[А2-, ^2 ,12 , J2 ](х) = | Ь (х)[1, У, у2, уЪ]йу .
У0( х)
Для сформированной системы (15) расчетные деформация и кривизна устанавливаются по выражениям типа (12) при анализе полигона предельных деформаций двухфазного сечения. В результате их нахождения из двух уравнений (15) могут быть найдены геометрические параметры, например, А1 и Ь2 или А1 и ^ = у2 - у1.
Из рассмотренной постановки (15) как частный случай, вытекает методика проектирования железобетонного сечения при следующих расчетных предпосылках: замена диаграммы деформирования бетона «нулевой» аппроксимацией (6)
стХ2)~ = />20 = -ст- ; пренебрежение сопротивлением растянутого бетона = 0. В таком случае вышерассмотренные уравнения принимают вид
-2О (У2 - У0 )b2 + EaA1 (б0* - УаК*) = N,
0,5ст0 (У22 -
У0 = 6* / К*,
0 520 (У2 - Уо )b2 + EaA1 Уа (б0* - УаК*) = M,
-0*
6аУ2 -6ъУа У2 - ya
К*( x) = -
У2 - Уа
где еа, е- - допустимые деформации арматуры и бетона. Из них легко находится площадь арматуры А1 и ширина прямоугольного сечения Ь2 .
Изложенный способ формирования физических зависимостей нелинейно деформируемых многофазных стержней на основе аппроксимации нормального напряжения степенным рядом позволяет алгоритмизировать запись физических соот-
a
Рис. 3. Расчетные и предельные деформации в сечении армированного стержня
Теория расчета строительных конструкций
ношений для интегральных силовых факторов с применением жесткостных характеристик сечений высших порядков; обоснованно, в соответствии с реальными свойствами используемых конструкционных материалов, формировать блочную структуру систем уравнений при построении алгоритмов решения обратных задач.
Литература
1. Немировский, Ю.В. Рациональное и оптимальное проектирование слоистых стержневых систем / Ю.В. Немировский, А.В. Мищенко, И.Т. Вохмянин. - Новосибирск: НГАСУ, 2004. - 488 с.
2. Васильев, В.В. Механика конструкций из композитных материалов / В.В. Васильев. - М.: Машиностроение, 1988. - 272 с.
3. Мищенко, А.В. Прямая и обратная задачи деформирования слоистых стержней с учетом физической нелинейности / А.В. Мищенко // Численные методы решения задач теории упругости и пластичности: труды XIX Всероссийской конференции, Бийск, 28-31 августа 2005 г. - Новосибирск: Параллель, 2005. - С. 184-188.
4. Мищенко, А.В. Напряжения в слоистых стержнях переменного сечения / А.В. Мищенко // Механика композиционных материалов и конструкций. - 2007. - Т. 13. - № 4. - С. 537-547.
5. Мищенко, А.В. Расчет и проектирование деревянных стержневых систем с учетом физической нелинейности / А.В Мищенко, Ю.В. Немировский // Строительная механика и расчет сооружений. - 2007. - № 6. - С. 46-52.
Мищенко Андрей Викторович, кандидат технических наук, доцент кафедры строительной механики, Новосибирский государственный архитектурно-строительный университет (Сибстрин), mavr@hnet.ru
Поступила в редакцию 14 января 2014 г.
Bulletin of the South Ural State University Series "Construction Engineering and Architecture" _2014, vol. 14, no. 2, pp. 12-16
METHOD OF FORMING NON-LINEAR PHYSICAL EQUATIONS IN THE DIRECT AND INVERSE PROBLEMS OF MULTIPHASE RODS ANALYSIS
A.V. Mishchenko, Novosibirsk State University of Architecture and Civil Engineering (Sibstrin), Novosibirsk, Russian Federation, mavr@hnet.ru.
The algorithmic method of forming non-linear physical equations for multiphase rods analysis is suggested. The method is based on the approximation of the diagrams of phase materials deformation by integer rational polynomial in arbitrary degree. The higher order rigidity characteristics are used.
Keywords: multiphase rod, armature rod, layered structure, physical non-linearity, higher order rigidity characteristics, limiting states.
References
1. Nemirovsky Yu.V., Mishchenko A.V., Vokhmianin I.T. Ratsional'noe i optimal'noe proektirovanie slois-tykh sterzhnevykh sistem [Rational and Optimal Design of Layered Rod Structures]. Novosibirsk, NSACEU Publ., 2004. 488 p.
2. Vasiliev V.V. Mekhanika konstruktsiy iz kompozitnykh materialov [Mechanics of the Composite Material Structures]. Moscow, Machinery construction Publ., 1988. 272 p.
3. Mishchenko A.V. [Direct and Inverse Problem of Deformation of Layered Rods with the Regard for Physical Nonlinear]. Chislennie metodi reshenia zadach teorii uprugosti i plastichnosti: Trudy 19 Vserossiiskoy konfe-rencii [Numerical Methods Solution of Elasticity and Plasticity Theory Problems: Transactions of XIX All-Russian Conference]. Novosibirsk, Parallel Publ., 2005, pp. 184-188 (in Russ.).
4. Mishchenko A.V., Nemirovsky Yu.V. [Analysis and Rational Designing of Wood Rod Systems with the Regard for Physical Nonlinear]. Stroitelnaia mekhanika i raschet sooruseniy [Structural Mechanics and Analysis of Building], 2007, no. 6, pp. 46-52 (in Russ.).
5. Mishchenko A.V. [Stresses on the layered rods of variable cross sections]. Mekhanika kompozitsionnykh materialov i konstruktsii [Composite Mechanics and Design], 2007, vol. 13, no. 4, pp. 537-547 (in Russ.).
Received 14 January 2014
