Новые операционные соотношения и их применение к решению нестационарных задач для стержней на основе теории С.П. Тимошенко

Полный текст

Введение 1 Одной из задач динамики сооружений является разработка принципов и методов расчета соору- жений при действии на конструкцию сил, приложенных лишь на короткий промежуток времени или быстро изменяющихся. Это явление надо рассматривать с точки зрения распространения волн напряжения [1]. Проблема распространения волн приводит к отысканию решения нестационарной задачи теории волн. International License указанных воздействиях классическая теория Бер- https://creativecommons.org/licenses/by/4.0/ В случае поперечных колебаний стержней при нулли - Эйлера во многих отношениях становится непригодной. Классическое уравнение относится к параболическому типу, что физически означает бесконечную скорость распространения возмущения. Теория С.П. Тимошенко [2], учитывающая в отличие от классической теории также инерцию вращения поперечных сечений стержня и деформацию сдвига, дает результаты, очень близкие к точным результатам, вытекающим из теории упругости. Уравнение Тимошенко является вполне гиперболическим и описывает распространение двух слабых разрывов (изгибного и сдвигового) с конечными скоростями. Для исследования переходных волновых процессов деформации в стержнях на основе теории Тимошенко необходимо иметь точные аналитические решения нестационарных задач в общем виде. Каждое точное аналитическое решение имеет большую значимость, так как в рамках данной аналитической модели является точным описанием реального процесса, служит критерием при оценке точности приближенных решений, позволяет отлаживать и верифицировать программы, в которых реализованы численные методы, является фундаментальным теоретическим фактом. Наиболее ясное представление о физической сущности нестационарных задач обычно дают решения в бегущих волнах. В общем виде такие решения отсутствуют. При использовании интегрального преобразования Лапласа по времени для построения решения в виде бегущих волн наибольшие трудности представляет именно этап перехода от изображения к начальной функции (оригиналу). Анализ работ по теории Тимошенко, в том числе упомянутых в обзоре [3] и современных [4-8], показал, что имеющиеся решения в бегущих волнах некоторых частных задач либо очень громоздки и имеют структуру, которая неясно выражает наиболее существенные черты описываемого процесса, либо их эффективность при вычислениях достигается только в некоторых довольно ограниченных областях значений координаты и времени. Формы представления этих решений не удовлетворяют, например, общим условиям, упомянутым в работе [9], для форм, в которых желательно представить решение. Исключениями являются исследования [10-12], в которых рассмотрены лишь частные задачи. Таким образом, построение в общем виде точных аналитических решений нестационарных задач теории С.П. Тимошенко для стержней в физически наглядной и удобной для практических расчетов форме является актуальной проблемой. В данной работе в общем виде сформулированы три типа нестационарных задач для полубесконечного стержня и описан метод получения указанных решений. Рассматриваются однородные призматические стержни. Решения строятся в виде бегущих волн методами операционного исчисления на основе интегрального преобразования Лапласа - Карсона по времени. Выводятся новые операционные соотношения, на основе которых показан способ нахождения оригиналов в наглядном и компактном виде, без использования общей формулы обращения. 1. Постановка задач и изображение их решений Для постановки исследуемых задач принимаем следующие условные обозначения: Ox, Oy, Oz - оси прямоугольной системы координат (x - координата вдоль оси стержня); t - время; Ab - площадь поперечного сечения стержня; Jb - момент инерции поперечного сечения относительно оси, проходящей через центр тяжести и параллельной оси Oy; k′ - коэффициент сдвига;  - объемная плотность; E - модуль Юнга; G - модуль сдвига; q - интенсивность внешней поперечной нагрузки; w - поперечное перемещение центра тяжести сечения стержня от положения его статического равновесия;  - угол поворота поперечного сечения, обусловленный изгибом, относительно нейтральной оси; M - изгибающий момент; Q - поперечная сила. Уравнения Тимошенко имеют следующий вид: 2ψ w  2ψ  EJb x2  k A G b x  ψ ρJb t2  0  (1) 2w ψ  2w  k A G b x2 x   ρAb t2  q В качестве основы для определения и исследования изображений решений нестационарных задач теории Тимошенко могут служить соотношения, полученные при решении в пространстве изображений нестационарной задачи для бесконечного стержня (  x  ) при нулевых начальных условиях, находящегося под действием поперечной сосредоточенной импульсной нагрузки, приложенной в момент времени t = 0 (q x t( , )Rδ1xδ1t, где R = const, 1 - единичная импульсивная функция). Точку приложения силы примем за начало координат. Вследствие симметрии рассмотрим только часть стержня, расположенную правей точки приложения нагрузки. Таким образом, математическая формулировка этой задачи следующая: найти решение системы уравнений (1) (q(x, t) = 0) при 0  x  , t  , удовлетворяющее граничным условиям ψ 0 ,t  0   R  (t  0) (2) Q0,t   δ1 t  2  и начальным условиям w x ,0  0, ψx,0  0,   w ψ 0  x . (3) x,0  0, x,0  0 t t  Ищем решение, обращающееся вместе со своими производными по x в нуль при x  . Аналогично рассматриваются три типа нестационарных задач для полубесконечного стержня (0  x  ) при нулевых начальных условиях и следующих условиях закрепления конца x = 0 и внешних нестационарных воздействиях, приложен- ных к концу x = 0: - задача I0 - конец x = 0 свободно оперт, опора смещается по закону w0(t), на конце x = 0 приложен изгибающий момент M0(t); - задача II0 - конец x = 0 защемлен, опора смещается по закону w0(t) и поворачивается по закону 0(t); - задача III0 - на свободном конце x = 0 приложены изгибающий момент M0(t) и поперечная сила Q0(t). Здесь w0(t), 0(t), M0(t), Q0(t) - заданные функции времени t, причем исходя из физического смысла будем считать функции w0(t) и 0(t) непрерывными при t  0 и равными нулю при t = 0. Математическая формулировка этих задач аналогична формулировке задачи, приведенной выше.  Пусть F p   p e pt f t dt   f t  представляет 0 изображение начальной функции (оригинала) f(t), где p  s iω - комплексная переменная (параметр интегрального преобразования Лапласа - Карсона). К сожалению, объем статьи не позволяет привести здесь все изображения решений сформулированных задач. Но, чтобы получить представление о структуре изображений в этих задачах, рассмотрим, например, изображение функции  (x,t) в задаче II0 при 0(t) = 0 . Оно записано так, чтобы выделить указанные ниже представители изображений двух типов: 1 x p,   C W0  p p p  p  p  p2 r2 erx1  p2 r2 er x2   c22 1 2  2 erx1   p r22 r12  r1  er x2    r2  . (4)      p2 rr r22 r12  c22 1  r1  c22 2  r2   Здесь r1,2(p) - непрерывные однозначные ветви многозначной функции r p( )  A1 , (5) выделенные на плоскости p с разрезами  A3 Re p  A3 ,  λ Im p λ и на плоскости p с разрезом  A3 Re p  A3 соответственно (рис. 1). Рис. 1. Области, в которых выделены однозначные ветви r1(p), r2(p) [Figure 1. Define of single-valued functions r1(p), r2(p) on the complex plane] зуемых здесь и далее: Приведем выражения для постоянных, исполь- A1 , A2  c122  c222 , A3 , c1  c2  c1 , c22 c12 c1 , v  c2 A4 A v3 1, λ  c c1 2 C. (6) Ветви r1 и r2 удовлетворяют на действительной оси следующим асимптотическим формулам: p r1 p ~ A1 1 A p2  , c1 p r2  p ~ A1 1 A p2  Re p ,Im p  0 . c2 В этом случае Rer p1   0 и Rer2 p 0 при Re p  A3. Изображения искомых функций содержат два типа функций. В настоящей статье будут рассмотрены основные представители этих типов: - изображения 1-го типа: n 2p2 2  cp22 r1p2  er1r x1  φ(1)n x,t, -  r2  r1  n 2p2 2  cp22 rp22  er2r x2  φ(2)n x,t, - r2  r1  где n = 0, 1; 2) изображение 2-го типа: p  κ t . p  rr1 2 2 c2 Изображения 1-го типа не содержат произведения r1r2, то есть не содержат радикал и характеризуются наличием радикала p  A . Изображения 2-го типа содержат произведения r1r2, то есть радикал p2  λ2 , и не содержат радикал p2  A32 . Отметим, что из рассмотренных задач изображения 2-го типа имеют место только в задачах II0 и III0. 2. Проблема определения начальных функций. Вывод новых операционных соотношений В практической реализации метода интегральных преобразований наиболее трудным этапом обычно является построение начальной функции (оригинала) по полученному изображению. Для этой цели можно было бы воспользоваться сразу формулой обращения Римана - Меллина [13]: f  t  1 ept F p( ) dp , 2πi p L где F(p) - изображение начальной функции f(t), L - контур Римана - Меллина в плоскости p, лежащий справа от всех особых точек функции F(p)/p. Однако очевидно, что такого способа определения оригинала следует по возможности избегать и вместо этого пытаться применить теоремы операционного исчисления [13], с помощью которых полученные изображения сводятся к имеющимся в справочниках по операционному исчислению (например, [14]). Перечислим следующие основные способы определения оригиналов, рассматриваемые в литературе [9; 15; 16]: 1) использование операционных теорем и табличных соотношений; 2) непосредственное использование формулы обращения; 3) представление изображения в виде ряда, допускающего достаточно простой почленный переход от изображения к оригиналу; 4) асимптотические разложения; 5) численное обращение. Каждый из перечисленных методов обладает достоинствами и недостатками. При этом, как было отмечено во введении, важно, чтобы формы представления оригиналов удовлетворяли некоторым условиям [9]. Это - обозримость результата, позволяющая судить об основных чертах исследуемого процесса, и возможность вычисления значений оригинала с заданной точностью. Итак, возвратимся к рассматриваемым задачам. В данной статье применяется первый способ. Следует отметить, что полученные изображения имеют весьма сложную структуру из-за сложности выражения (5). Видимо, поэтому в известных нам работах, в которых уравнения Тимошенко решались с помощью преобразования Лапласа по t (см. введение), не был использован первый способ. Исключениями являются работы [10-12], в которых рассмотрены лишь частные задачи. Однако нам удалось в общем случае пойти по этому пути - пути построения оригиналов с помощью операционных теорем и соотношений. Но для этого выводятся новые, не имеющиеся в литературе по операционному исчислению, операционные соотношения. Одно из этих соотношений (назовем его основным) дает выражение начальной функции для изображения следующего вида: F a p 1  a p  α2  , (7) 1 где F p( )  f t  - некоторое известное (табличное) соотношение; a1 a2 0 (а1 и a2 не равны нулю одновременно); , m - комплексные числа, Re m > -1. Таким образом, таблица операционных соотношений [14] будет дополнена новыми соотношениями. Оригиналы для всех изображений, полученных при решении рассматриваемых задач, можно найти с помощью этой расширенной таблицы и некоторых теорем операционного исчисления. Как будет показано далее, при таком способе нахождения оригиналов решения поставленных задач, полученные в явной аналитической форме, представляются в виде интегралов, но имеют ясно выраженный волновой характер и компактный вид. Очень важно также, что и для стержней конечной длины с помощью перечисленных средств можно получить в явной аналитической форме оригиналы для изображений прямых и отраженных волн. Отметим, что при построении оригиналов рассматриваемым способом ограничиться чисто операционными методами не удается. Связано это с тем, что, хотя операционное решение более коротко, более просто, в нем явно не выражается наличие условий, достаточных для законности той или иной операции. Поэтому в сложных случаях (на- пример, при наличии многозначных функций, как в нашем случае) для того, чтобы аналитический смысл каждой операции был совершенно ясен, подробно исследовано поведение изображений в комплексной плоскости. Основное операционное соотношение найдем, применяя преобразование А.М. Эфроса (обобщенную теорему Бореля) [17]. Укажем здесь два варианта вывода этого соотношения. 1. Непосредственно используем преобразование Эфроса при a1  a2  0 ,  a2  0. Приводим окончательные результаты. При a1  a2  0,  a2  0 получим m2    a 0 a2  0 , Rem  1. (8) Здесь Im - модифицированная функция Бесселя первого рода порядка m [18]. При a2 = 1 мы получаем имеющийся в литературе по операционному исчислению (например, в [14]) частный случай преобразования Эфроса. При a1  a2  0 ,  a2  0 , обозначив a  a1  a2 , (9) a1  a2 Rem  1, (10) где a1 и a2 удовлетворяют системе a1  a2  0 ,  a2  0. 2. Другой вариант вывода искомого операционного соотношения основан на использовании известных из литературы частных случаев преобразования Эфроса. При этом принимаем, что a1  a2  a1. Тем самым будет рассмотрен также случай a1  a2  0 (a1 и a2 не равны нулю одновременно), не охваченный в первом варианте вывода. Данный вариант вывода искомого операционного соотношения рассмотрим подробно. Допустим сначала, что a1  a2  a1. Изображение в левой части соотношения (10) можно представить как результат следующей цепочки преобразований: a1  a2 α2 1 F p  F  a1  a2 p  2  p  a p 2α αm1p 1 Для получения соответствующей каждому этапу преобразования начальной функции используем известные из литературы следующие частные случаи преобразования Эфроса [14]: если F p( ) f t( ), то F1 p  F a1  a2 p f  2 t f t1( ) , (11)  2   a1  a2  где J0 - бесселева функция первого рода нулевого порядка, 0 Rem  2 . (13) Следовательно, для получения искомого ори- 2α гинала надо найти функцию f3( ).t a При -a1 = a2 < a1 получим соотношение a1  0, Rem  1, (14) которое следует из соотношения (13) при  t  F2( )p  F a p( 1 )  f   , где вместо m взято m - 1.  a1  Соотношение (13) является, таким образом, частным случаем соотношения (14). Обратимся снова к случаю -a1 < a2 < a1. Подставим в (13) выражение для функции f2(τ) из (12) и переменим порядок интегрирования, используя формулу Дирихле для двойных интегралов [19]. Следовательно, m1   f3( )t  t 2 0 f1 1(τ ) τd 11(t 2τ)Im1α t t( 2τ) J02α a1τ (τ τ )1  1 dτ,  Rem 2. Преобразуем внутренний интеграл, введя но- Перейдем в (15) от модифицированной бесселевой функции первого рода Im+1 к бесселевой функции первого рода Jm+1. При этом удобней рассматривать не функцию вида Im(z) (z - комплексная переменная), а функцию Im( )z . Испольzm зуя [18], приходим к следующим соотношениям: Im( )z  Jm(iz) , (16) zm (iz)m Jm ( )z  Im(iz). (17) zm (iz)m Отметим, что функции Im( )z и Jm( )z являzm zm ются четными. Тогда выражение (15) согласно (16) преобразуется к виду, содержащему разрывный интеграл Сонина [18]:  1 0, t  2a τ ,1   αt 2 1  Jm α 2τ1 tt 2a1 τ1, t  21a τ .1 Для обеспечения сходимости этого интегра- 1 ла под величинами α 2t и α τ1 подразумева- a ются положительные величины, то есть α 0 , а Re (m+1) > 0 [18]. Ограничение α 0 вызвано только способом вывода, а не существом дела: как сказано выше, основное операционное соотношение, которое будет получено далее, справедливо и для комплексных . Далее получим m at f3( )t  1 2 t  21 τ1  2   1   2   a Im α t  2τ1t  2 τ1  f  τ1dτ .1  α  t  2τ1    a   a1 a2  0   Переходя к новой переменной интегрирова- 1 ния τ  2 τ1, получим соотношение (10), но при a поставленном условии a1  a2  a1. Отметим, что второй вариант вывода основного операционного соотношения послужил нам методической основой при первоначальных попытках нахождения оригиналов для частных видов изображения (7). Объединяя условия для a1 и a2, поставленные при обоих вариантах вывода, получаем искомое операционное соотношение (основное) t  τ где, в соответствии с обобщенной теоремой Бореля, F(p) удовлетворяет условиям обратной теоремы относительно формул Римана - Меллина [17] (в дальнейшем об этих условиях для F(p) упоминать не будем; они будут всегда выполнены), a  a1  a2 , a1  a2  0, Rem  1, a1  a2 а также соотношения (8) и (14) при a1  a2  0 , Re m > -1. Смысл верхнего индекса в обозначении Sm(2) станет понятен далее. Из соотношения (18) при a1 = 0, a2 = 1, а также при a1 > 0, a2 = 0 следуют имеющиеся в литературе частные случаи преобразования Эфроса [14]. Далее нам потребуются некоторые дополнительные соотношения. 1) Заменяем в соотношении (18) a2 на -a2 и переходим к новой переменной интегрирования τ τ1  (ниже обозначим ее опять через  ), получим a Sm(1)(a a1, 2,α,F p( ), p)  Sm(2)(a1,a2,α,F p( ), p)  a a1 2αm  dτ, (19) t  τ  a1 a2  a  a1  a2 , a1  a2  0, Rem  1. a1  a2 2) Получим оригинал для случая, когда m = -1 ( S1(2)( )p является изображением): S1(2)( )p  p (2)( )p  1  (2)(2)   α 2pS0 ( )p  p  S0 ( ) .p  Здесь не перечислены для сокращения записи все параметры, от которых зависит функция Sm(2) . Оригинал найдем из (18): для первого изображения в квадратных скобках - по теореме дифференцирования оригинала [13] (оригинал в правой части соотношения (18) при t = 0 равен нулю); для второго изображения - из соотношения  p  p2  α2 S0(2)( )p  αS1(2)( ).p В результате имеем a  a1  a2 , a1  a 0. (20) 2 a1  a2 Заменяя здесь a2 на -a2 и переходя к новой переменной интегрирования, получим a  1  a2 , a1  a2  0. (21) a1  a2 3) Сравнивая (18) и (19) при m = 0, а также (18) при m = 1 и (21), отметим, что соответствующие интегралы отличаются только верхним пределом интегрирования. Тогда напишем в обобщенном виде очень полезное соотношение: Sn(a a1, 2,α,F p( ), p)  Sn(1)(a a1, 2,α,F p( ), p)      f  dτ  n f  . (22)  a1  a2  aα  a1  a2  Здесь должно быть a1 - a2 > 0 и a1 + a2 > 0, то есть a1  a2 ; n = 0, 1. Используя (19) при m = 1 и (20), получим S a a%1( 1, 2,α,F p( ), p)  S1(1)(a a1, 2,α,F p( ), p)  t 2  t   f  , a1 a2 . (23) α  a1a2 Легко заметить, что S a a%1 1 2( , ,α,F p p( ), )  1  S a1 1( ,a2,α,F p( ), p) . a 4) Заменяя в соотношении (18)  на i и используя (17), получим t m  t  τ 2 α t  aτt  τ f  τ    t  aτ  Jm    a1 a2 dτ, 0 a1 a2  0, Rem 1. (24) При a1 = 0, a2 = 1 из (24) следует имеющийся, например, в [14], частный случай преобразования Эфроса. Заменяя  на i в соотношениях (8) и (14), получим новые операционные соотношения, которые при a2 = 1 и a1 = 1 соответственно сводятся к имеющимся в [14]. Таким образом, найден ряд новых операционных соотношений, с помощью которых будут определены начальные функции в рассматриваемых задачах теории Тимошенко. 3. Определение начальных функций с помощью новых операционных соотношений Возвратимся к нестационарным задачам теории Тимошенко для бесконечного и полубесконечного стержней. Изображения решений этих задач не- обходимо представить в таком виде, чтобы применить новые операционные соотношения и некоторые известные теоремы операционного исчисления к отдельным составляющим изображений (см., например, формулу (4)). Если эти составляющие относятся к разным типам (1-му или 2-му), то применяются разные новые операционные соотношения. Покажем, как определить начальные функции с помощью новых операционных соотношений. 1) Начальные функции для изображений 1-го типа. Легко доказать, что r1 и r2 получаются из выражения при замене p на  1  1  p 1  v   и 1 1 p 1  v   соответственно. Отсюда следует, что при использовании соотношений (19), (21) и (18), (20) соответственно надо принять 1 1 a1 1 , a2  1 , α  A3 . v v Тогда, согласно (9), 1 a  . v При этом за исходное соотношение F p( ) f t( ) будем принимать следующее табличное соотношение [14]: F p( )  2c2   x 0,   f t( )   (25) 2c J2 0 A3 1.   Допустим, требуется найти оригинал для p2 er x1 изображения . Используя соотношение r22 r12 r1 (19), получим где f   τ принимается согласно (25).  2 В результате получим p2 erx1 (1) r22 r12 r1 φ0 x,t   x 0, t  c1 ,   (26) .    c1 x x При t  имеем  τ t , таким образом, выражеc1 c1 ние t vτ может быть больше, меньше или равно нулю. В зависимости от знака этого выражения имеем или функцию I0 A3 (t  τ)(t vτ) , или функцию J0 A3 ( τv t t)(  τ) (от действительного аргумента). Следовательно, рассмотрим два случая, x t определяемые соотношением и : c1 v t x x x а)   t , или  t ; тогда на всем инv c1 c1 c2 t тервале интегрирования τ  ; v x t x x б)   t , или t   ; тогда интервал инc1 v c2 c1 тегрирования следует разбить на две части: x t t  τ  и  τ  t . c1 v v p2 er x2 Оригинал для изображения мож- r22 r12 r2 но получить, применяя соотношение (18), так как Тогда p2 er x2 (2) r22 r12 r2 φ0 (x,t)   x 0, t  c2 ,    (27) .    c2 Теперь становится понятен смысл верхнего индекса в обозначениях Sm(1) и Sm(2) : изображения рассматриваемого типа, содержащие множитель erx1 (1) , приводятся к изображению Sm , а изобраr1 er x2 жения, содержащие множитель , приводятся r2 к изображению Sm(2) , где по-прежнему F(p) принимается согласно (25). Для сравнения приведем теперь в наших обозначениях оригинал, например, для изображения p2 er x1 , найденный в работе [20] с использоr22 r12 r1 ванием общей формулы обращения для преобразования Лапласа:  x     1 1   πA1 C 0     sh A 3 1ρ t  A A x  r x Сравнение при t > x/c1 (26) с (28) показывает, что оригинал в соотношении (26), в отличие от оригинала в соотношении (28), имеет компактный, наглядный вид, что способствует более простому анализу оригинала. Мы не имеем еще полного решения конкретной задачи, однако вид оригиналов в соотношениях (26) и (27) подтверждает известный факт: теория Тимошенко описывает распространение двух волн со скоростями фронтов c1 и c2; структура указанных оригиналов (нижний и верхний пределы интегралов, вид аргумента функ-  2  ции типа J0 A4 ) имеет ясно выражен-  ный волновой характер, что нельзя сказать об оригинале в соотношении (28). Характерным обстоятельством является то, что в изображениях, полученных при решении задач для бесконечного и полубесконечного стержней, одна половина слагаемых (при er x2 ) получается из другой (при er x1 ) заменой r1 на r2, и наоборот. В формуле (4) это изображения в квадратных скобках. Например, следующие два изображения всегда находятся вместе в такой комбинации: p2 er x1 p2 er x2  . (29) r22 r12 r1 r22 r12 r2 x Если при t  в первом слагаемом оригинала c2 φ0(1)(x,t) перейти к новой переменной интегрирования τ1  vτ , то видно, что это слагаемое взаимно x (2) уничтожается при t  с оригиналом φ0 (x,t) . c2 В связи с этим лучше сразу использовать совместные соотношения (22) и (23) при нахождении оригиналов для изображений типа (29). Например, из (22) в обобщенном виде получим n n 2 2 rx1 2 2 r x2 p  p r1  e p  p r2  e         φ (x,t)  2 2 2 2 2 2 n r2  r1  c2 p  r1 r2  r1  c2 p  r2   x x x , t   1  , t  ,  c2 n  0,1 . 2) Начальные функции для изображений 2-го типа. Допустим, требуется найти оригинал для p изображения . Так как p 2  r r1 2 c2 p  c12 p 1 c c1 2 , p 2  r r1 2 c2 то оригинал для первого слагаемого найдем, применяя соотношение (24) для c c1 22  f t( )   c c1 22 F p( )   c1  c2 c1  c2 при a1 v a, 2 1, α  λ, m  0; оригинал для второго слагаемого находим по [14]; тогда получим: p  p 2  rr12 c2  c c1 22 t J0λ t  c1c2 τt  τdτ  c c1 2t J0λτdτ  c1 c2 0   c1 c2   0  κ( )t . (31) Итак, оригиналы для изображений рассмотренных типов представляются в виде интегралов от бесселевых функций и определяются практически с помощью двух операционных соотношений (18) и (25). Аналогично определяются все оригиналы в рассматриваемых задачах. Точные аналитические выражения для перемещений и усилий в бесконечном и полубесконечном стержнях сведены в таблицы. Ограничение объема статьи позволяет привести здесь только выражения для угла поворота и изгибающего момента в рассмотренной далее тестовой задаче. 4. Решение тестовой задачи для полубесконечного стержня Рассмотрим на основе теории Тимошенко частную нестационарную задачу о действии на защемленный полубесконечный стержень ( x0 ) равномерно распределенной поперечной импульсной нагрузки, приложенной в момент времени t = 0; начальные условия - нулевые. Решение этой задачи, найденное с использованием общей формулы обращения для преобразования Лапласа, имеется в работе [21] (решение дано без вывода). Поставим перед собой целью сравнить числовые результаты, полученные по формулам данной работы и в работе [21]. При решении этой задачи была отработана методика вычисления значений искомых функций на ЭВМ. Математическая формулировка этой задачи сле- дующая: найти решение системы уравнений (1) при q x t( , )  Rδ ( )δ ( )0 x 1 t ( R  const; δ ( ),δ ( )0 x 1 t - единичная функция Хевисайда и единичная импульсивная функция), 0   x , t 0 , удовлетворяющее граничным условиям w(0,t)  0, ψ(0,t)  0 (t  0) и начальным условиям (3). Ищем решение, ограниченное при x . Будем искать функцию w (x, t) в виде R w x t( , )  w x t1( , )  t. (32) ρAb В результате приходим к задаче II0 относительно функций w1 (x, t) и  (x, t) при R w0( )t   t , ψ ( )0 t  0. ρAb В работе [21] дано выражение для безразмерного изгибающего момента M(x t, ) . Чтобы полу- R чить выражение для этой функции, используя новые операционные соотношения, удобней найти сначала оригинал для изображения Ψ (x, p) (4), а затем - ψ изгибающий момент из соотношения M  EJb . x Используя (30), (31) и теорему свертыва- ния [13], получим  EJb (x, p)  EJb ψ(x t, )  R R t t φ ( ,τ)κ(1 x t  τ) τd φ ( ,τ) τ.0 x d (33) 0 0 x x В зависимости от соотношения t, , интерваc1 c2 лы интегрирования в (33) можно разбить на интервалы, в которых подынтегральная функция имеет конкретное аналитическое выражение. Доказав, что эти интегралы имеют частную производную по x в соответствующих областях, получим CM x t( , )   EJb C ψ ( x,t)  R R     [1]0 , t  x ,  c1      t    % (x t, ,τ) τd  κ(t  x ) , cx1  t  cx2 , (34)   C  f c1   x  c1    x   c2 t   C  f% (x t, ,τ) τd  x f% (x t, ,τ) τ κ(d  t  x ) , t  cx2 , c1 x  f% (x t, ,τ)  φ ( ,τ)κ(1 x t  τ) φ0 x,τ . На защемленной опоре: M (0,t) C  Cκ( ).t R В сравниваемой работе [21] безразмерные из- CM гибающие моменты вычислялись для стерж- R ня с прямоугольным поперечным сечением высотой h. Введем, как и в [21], безразмерную координату точки оси стержня и безразмерное время: x c t1 . η  , T  h h Анализ работы [21] показал, что значение v было принято равным 1,9. Построенный по результатам вычислений гра- CM фик функции (рис. 2 ) при η = 2 совпадает R с графиком, приведенным в работе [21]. T Рис. 2. Действие на защемленный полубесконечный стержень ( x  0 ) равномерно распределенной поперечной импульсной нагрузки q x t( , )  Rδ ( )δ ( )0 x 1 t . Изгибающий момент в сечении η = 2 при v = 1,9 [Figure 2. A clamped semi-infinite beam ( x  0 ) which is subjected to a uniform lateral impulse q xt( , )Rδ ( )δ ( )0 x 1 t . Bending moment at η = 2, v = 1,9] Выводы Итак, в ходе исследования показан метод получения точных аналитических решений задач о распространении упругих волн в однородных приз- матических стержнях при поперечных кратковременных и быстроизменяющихся воздействиях на основе теории Тимошенко. В общем виде сформулированы три типа нестационарных задач для полубесконечного стержня. Рассмотрены бегущие волны. Использовано операционное исчисление. Отмечено, что полученные изображения имеют очень сложную структуру. Изучена проблема определения начальных функций (оригиналов) при решении нестационарных задач теории Тимошенко. Выводятся новые, не имеющиеся в литературе по операционному исчислению, операционные соотношения, на основе которых показан способ нахождения оригиналов без использования общей формулы обращения. Решения задач записываются в виде интегралов от бесселевых функций и, в отличие от решений, имеющихся в литературе, ясно показывают волновой характер изучаемых процессов, имеют компактный вид и удобны для практических вычислений. Ограничение объема статьи позволило привести здесь только некоторые точные аналитические выражения для различных функций. В общем виде такие выражения для различных типов стержней получены и сведены в таблицы, которые могут быть представлены автором статьи. Полученные новые операционные соотношения позволяют дополнить таблицы формул операционного исчисления и могут быть использованы в дальнейшем при решении различных проблем.

Об авторах

Александр Леонидович Зоненберг

АО «ЦНИИЭП жилища - Институт комплексного проектирования жилых и общественных зданий»

Автор, ответственный за переписку.
Email: zonenberg@list.ru

главный специалист-инженер, отдел конструкций жилых и общественных зданий

Российская Федерация, 127434, Москва, Дмитровское шоссе, д. 9, стр. 3

Список литературы