Алгоритм метода перемещений с параметрическим изменением внешних силового и кинематического воздействий

Бесплатный доступ

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

Строительная механика, метод перемещений, односторонние связи, линейная задача дополнительности, метод лемке

Короткий адрес: https://sciup.org/142142423

IDR: 142142423

Текст научной статьи Алгоритм метода перемещений с параметрическим изменением внешних силового и кинематического воздействий

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

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

Задача расчета систем с идеальными (без учета трения) односторонними связями после дискретизации сводится к условной задаче квадратичного программирования [1]:

П (z) = — zTRz - zTR F ^ min,

2                                                             (1)

z > 0.

Здесь z – взаимное перемещение контактных пар точек по направлению односторонних связей; R – контактная матрица жесткости; R F – вектор сил взаимодействия контактных пар точек по направлению односторонних связей от заданной нагрузки (грузовой вектор). Перемещения принято считать положительными Z i >  0 , если они не ограничиваются односторонней связью.

Обозначим ~ = z + п , где П — вектор известных обобщенных зазоров по направлению z, и запишем условия Куна-Таккера для задачи (1):

~

1) x = R ~ + R F ; 2) x 0; 3) ~ 0; 4) ~ Tx = 0,

~ где x - вектор сил взаимодействия контактных пар точек; RF = RF - Rn.

Система уравнений 1, условий неотрицательности 2, 3 и условий дополняющей нежесткости 4 из (2) представляют собой линейную задачу дополнительности (ЛЗД) [3, 4]. Алгоритм решения данной задачи можно интерпретировать как метод перемещений строительной механики с применением терминов: основная система метода перемещений, единичные состояния (дислокации); взаимные перемещения контактных пар точек; контактная матрица жесткости и т. п.

Для решения задачи (2) предлагается алгоритм с параметрическим изменением внешнего воздействия, являющийся модификацией алгоритма Лемке [3, 4]. При параметрическом изменении внешнего воздействия R F = pR V +R C, где pR V , R C – изменяющаяся и постоянная составляющие внешнего воздействия; p – параметр изменения внешнего воздействия. Тогда ЛЗД (2), расширенная введением параметра p , перепишется в виде:

1) x - R~ - pR V = R C ; 2) x 0; 3) ~ 0; 4) ~ T x = 0; 5) p 0. (3)

Расширенная задача (3) будет эквивалентна исходной задаче (2) при p = 1.

Ясно, что только одна переменная из каждой пары xi , zi в силу 2) – 4) из (3) будет отлична от нуля. Условие 4) из (3) устанавливает правило дополняющей нежесткости для определения индекса переменной, вводимой в базис. А именно в базис вводится переменная, взаимодополнительная к переменной, вышедшей из базиса. Если из базиса вышла (стала равной нулю) переменная xi (что соответствует выключению связи в контактной паре i – «момент отрыва»), то в базис должна войти переменная zi , определяющая взаимное удаление точек контактной пары i . Если из базиса вышла переменная zi (включение связи, «момент касания»), то в базис, естественно, войдет xi – сила взаимодействия вступивших в контакт точек.

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

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

Компоненты контактной матрицы жесткости R представляют собой реакции во введенных связях от «единичных» дислокаций опор 1, 2, 3 (рис. 1 а , в , г , д ). С целью определения значений этих компонент построены «единичные» эпюры M 1 , M 2 , M 3 (рис. 1 е , ж , з ).

Балка загружена двумя сосредоточенными силами F 1 и F 2 . Силы прикладываются к конструкции поочередно: нагрузка F возрастает до значения 2,5 кН , затем при действующей силе F 1 = 2,5 кН начинает возрастать нагрузка F 2 до значения 7,5 кН .

Грузовой вектор метода перемещений формируется из реакций во введенных связях основной системы от заданной нагрузки. С целью формирования грузовых векторов R F1 и R F2 построены эпюры M 01 и M 02 соответственно от сил F 1 и F 2 (рис. 1 и , к ).

В результате имеем:

" 4548.48 -3979.92 -1137.12 "

'- 0 . 375

' 1 . 125 '

' 0 1

R =

-3979.92 4548.48  1705.68

, R F 1 = '

1 . 813

- , R F 2 ='

- 5 . 44

- , з = ■

0 . 001

-1137.12 1705.68   758.08

. - 1

. 3

. 0

а)

I M 20

F 1 = 2.5 кН

' Г

1 мм

2    F 2

Основная сиситема

= 7.5 кН

д )

F. = 1705.68 r = 1137.12

2274.24

3411.36

2274.24

0.125

r F1 = 0.375   F = 1.813

1516.16

r F 1 = 1

Рис. 1

На первом этапе алгоритма снимаем преднапряжение основной системы метода перемещений, вызванное замкнутыми зазорами, т. е. полагаем RV = 0, что означает отсутствие силовой нагрузки. Ис- пользуя алгоритм с искусственно введенным обжатием [2], формируем начальную таблицу ЛЗД и про-< 3,98 )

изводим расчет конструкции на воздействие R C = — R n =

- 4,549

(— 1,706 )

Грузовые векторы R F1 и R F2 располагаются в разных столбцах и преобразуются вместе со всей таблицей в результате выполнения шагов алгоритма на этапе I.

Таблица 1

Начальная таблица этапа I

Базис

x1

x2

x3

z1

z2

z3

p1

p2

p3

R3

x1

1

0

0

-4548,5

3979,9

1137,1

-1

-0,375

1,125

3,98

x2

0

1

0

3979,9

-4548,5

-1705,7

-1

1,813

-5,438

-4,549

x3

0

0

1

1137,1

-1705,7

-758,08

-1

-1,0001

3,0004

-1,706

Базисное решение x i = 3 . 98 , x 2 = - 4 . 549 , Х 3 = - 1 . 706 недопустимо из-за нарушения переменными x 2 и x 3 условия неотрицательности [4]. Для решения ЛЗД на шаге 1 в базис следует ввести параметр p1 такой, чтобы все базисные переменные оказались большими нуля. Физически это означает введение в опорных сечениях такой прижимающей силы p1 (одинаковой для всех опор), которая для одной из опор обеспечит «момент отрыва».

В результате выполнения шага 1 имеем следующее.

Таблица 2

Шаг 1

Базис

x1

x2

x3

z1

z2

z3

p1

p2

p3

- Кз

min

x1

1

-1

0

-8528,4

8528,4

2842,8

0

-2,1876

6,563

8,5284

0,001

p1

0

-1

0

-3979,9

4548,48

1705,68

1

-1,8126

5,438

4,5485

0,001

x3

0

-1

1

-2842,8

2842,8

947,6

0

-2,8128

8,4383

2,8428

0,001

Базисные переменные: x i = 8 . 528 , pl = 4 . 549 , Х 3 = 2 . 843 . Величина обжатия pl такова, что деформированная балка отрывается от опоры 2 (рис. 2 а ). Поскольку из базиса вышла переменная x 2 , то на следующем шаге согласно правилу дополняющей нежесткости x 2 z 2 = 0 в базис должна войти переменная z 2 . Для выбора ведущей строки воспользуемся правилом минимального отношения. Поскольку величина отношения одинакова во всех трех случаях, выбираем ведущую строку так, чтобы переменная p1 вышла из базиса, и проводим шаг однократного замещения Жордана-Гаусса. В результате имеем следующее.

Таблица 3

Шаг 2

Базис

x1

x2

x3

z1

z2

z3

p1

p2

p3

- Кз

x1

1

0,875

0

-1066

0

-355,35

-1,875

1,211

-3,633

0

z2

0

-0,0002

0

-0,875

1

0,375

0,0002

-0,0004

0,0012

0,001

x3

0

-0,375

1

-355,4

0

-118,45

-0,625

-1,6799

5,0397

0

Базисное решение: x i = 0 , z 2 = 0 . 001 м , x 3 = 0 . Полученный результат показывает, что смещение односторонней опоры 2 вертикально вниз и, следовательно, выключение ее из работы при отсутствии внешней нагрузки не влияет на напряженное состояние балки.

Параметр обжатия pl вышел из базиса. Следовательно, pl = 0 и этап I завершен (рис 2 б).

На втором этапе решаем ЛЗД с R C = - R n( I 1 и R V = R F 11 1 1 . Начальная таблица этапа II получается из таблицы, соответствующей концу первого этапа. Векторы R F /11 и R F2W переносятся влево за знак равенства с противоположным знаком. Параметр p2 характеризуют изменение внешней силовой нагрузки F 1 . Параметр p3 на данном этапе считается равным 0 в связи с отсутствием силового воздействия F 2 .

Таблица 4

Начальная таблица этапа II

Базис

x1

x2

x3

z1

z2

z3

p1

p2

p3

- Кз $

min

x1

1

0,875

0

-1066

0

-355,4

-1,875

-1,211

3,6332

0

z2

0

-0,0002

0

-0,875

1

0,375

0,0002

0,0004

-0,001

0,001

2,5093

x3

0

-0,375

1

-355,4

0

-118,5

-0,625

1,6799

-5,04

0

0,0000

При p2 = 0 базисное решение: x i = 0 , z 2 = 0 . 001 м , Х 3 = 0 , что соответствует отсутствию силовой нагрузки. Для начала итерационного процесса вводим в базис параметр нагрузки p2 (таким образом, ведущий столбец соответствует p2 ). По правилу минимального отношения определяем номер ведущей строки и проводим шаг однократного замещения Жордана-Гаусса. В результате получим следующее.

Таблица 5

Шаг 3

Базис

x1

x2

x3

z1

z2

z3

p1

p2

p3

- R3 $

min

x1

1

0,6047

0,721

-1322

0

-440,7

-2,326

0

0

0

z2

0

-0,0001

-0,0002

-0,791

1

0,4031

0,0004

0

0

0,001

0,002

p2

0

-0,223

0,595

-211,5

0

-70,51

-0,372

1

-3

0

В таблицах 6-8 (этап II: шаги 4 – 6) показано выполнение шагов алгоритма с параметрическим изменением нагрузки F 1 . На рисунке 2 в , г , д , е представлены рабочие схемы, соответствующие каждому шагу алгоритма.

Таблица 6

Шаг 4

Базис

x1

x2

x3

z1

z2

z3

p1

p2

p3

- R3 $

min

x1

1

0,4615

0,4615

-2186,8

1093,4

0

-1,9231

0

0

1,0934

2,37

z3

0

-0,0003

-0,0006

-1,9616

2,4808

1

0,0009

0

0

0,0025

p2

0

-0,2461

0,5538

-349,84

174,92

0

-0,3077

1

-3

0,1749

Таблица 7

Шаг 5

Базис

x1

x2

x3

z1

z2

z3

p1

p2

p3

- R3

min

x2

2,1667

1

1

-4738

2369

0

-4,1667

0

0

2,369

z3

0,0007

0

-0,0003

-3,5

3,25

1

-0,0004

0

0

0,0033

p2

0,5333

0

0,79991

-1516

757,99

0

-1,3332

1

-3

0,758

Таблица 8

Шаг 6

Базис

x1

x2

x3

z1

z2

z3

p1

p2

p3

- R3

x2

0,5

1

-1,5

0

0

0

0

-3,1254

9,376

0

z3

-0,0005

0

-0,0021

0

1,5

1

0,00264

-0,0023

0,0069

0,0015

z1

-0,0004

0

-0,0005

1

-0,5

0

0,00088

-0,0007

0,002

-0,0005

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

Базисные переменные x 2 , z 3 , z 1 получаются вычитанием столбца, соответствующего p2 , из столбца, соответствующего - R3 ^ :

B = •

f 1 2

z 3

>

f 0

= 0 . 0015

М

' - 3 . 125 '

- 0 . 0023

> = <

3 . 125

0 . 0038

. Z 1

_- 0 . 0005

_- 0 . 0007

0 . 0002

Начальная таблица этапа III формируется из результатов, полученных по окончании этапа II, с учетом того, что R C = B и R V = R F 2 ( II > .

В таблице 9-14 (этап III: шаги 7 – 11) представлены результаты выполнения шагов алгоритма с параметрическим изменением нагрузки F 2 с учетом приложенной силы F = 2 . 5 кН . На рисунке 2 ж , з , и , к , л приведены рабочие схемы, соответствующие каждому шагу алгоритма.

Таблица 9

Начальная таблица этапа III

Базис

x1

x2

x3

z1

z2

z3

p1

p2

p3

B

min

x2

0,5

1

-1,5

0

0

0

0

-3,1254

9,3761

3,1254

0,333

z3

-0,0005

0

-0,0021

0

1,5

1

0,0026

-0,0023

0,0069

0,0038

0,55

z1

-0,0004

0

-0,0005

1

-0,5

0

0,0009

-0,0007

0,002

0,0002

0,1

Таблица 10

Шаг 7

Базис

x1

x2

x3

z1

z2

z3

p1

p2

p3

B

min

x2

2,1667

1

1

-4738

2369

0

-4,1667

0

0

2,369

1,093

z3

0,0007

0

-0,0003

-3,5001

3,25

1

-0,0004

0

0

0,0033

4,62

p3

-0,1778

0

-0,2666

505,328

-252,66

0

0,4444

-0,3333

1

0,0807

Таблица 11

Шаг 8

Базис

x1

x2

x3

z1

z2

z3

p1

p2

p3

B

min

x1

1

0,4615

0,4615

-2186,8

1093,4

0

-1,9231

0

0

1,0934

0.001

z3

0

-0,0003

-0,0006

-1,9616

2,4808

1

0,0009

0

0

0,0025

0.001

p3

0

0,082

-0,1846

116,614

-58,307

0

0,1026

-0,3333

1

0,275

Таблица 12

Шаг 9

Базис

x1

x2

x3

z1

z2

z3

p1

p2

p3

B

min

x1

1

0,6047

0,7209

-1322,2

0

-440,74

-2,3256

0

0

0

0

z2

0

-0,0001

-0,0002

-0,7907

1

0,4031

0,0004

0

0

0,001

p3

0

0,0744

-0,1984

70,511

0

23,5036

0,124

-0,3333

1

0,3333

Таблица 13

Шаг 10

Базис

x1

x2

x3

z1

z2

z3

p1

p2

p3

B

min

x3

1,3871

0,8387

1

-1834,1

0

-611,36

-3,2258

0

0

0

z2

0,0003

0,00007

0

-1,2258

1

0,2581

-0,0004

0

0

0,001

p3

0,2752

0,2408

0

-293,42

0

-97,806

-0,5161

-0,3333

1

0,3333

Шаг 11

Базис

x1

x2

x3

z1

z2

z3

p1

p2

p3

B

x3

-0,3333

-0,6667

1

0

0

0

0

2,0836

-6,2507

-2,0836

z2

-0,0008

-0,0009

0

0

1

0,6667

0,0018

0,0014

-0,0042

-0,0004

z1

-0,0009

-0,0008

0

1

0

0,3333

0,0018

0,0011

-0,0034

-0,0011

в ) ШагЗ

F 1 = 0

и) Шаг 9

p3 = 0.333

p2 = 0

z 2 = 1 мм $

Момент

отрыва

x 1 = 0

F1=2.5кН z2=1ммф     '

Момент касания

2 . 5 кН

x 1 = 0

к) Шаг 10

p3 = 0.333

F 1 = 2 . 5 кН

x 3 = 0

= 1 мм

Момент

2.5кН отрыва

Рис. 2

Этап III также завершается получением лучевого решения (шаг 13). Значения базисных переменных x 3 , z 2 , z 1 получаются вычитанием столбца, соответствующего p3 , из столбца, соответствующего B :

' x 3 '

' - 2 . 084 '

' - 6 . 251 '

' 4 . 167 '

z z 2

> = <

- 0 . 0004

- 0 . 0042

> = <

0 . 0038

, z

_- 0 . 0011

_- 0 . 0034

0 . 0023

Таким образом, в результате произведенных расчетов установлена действительная рабочая схема конструкции, определены «включенные» и «выключенные» из работы односторонние связи, получены численные значения реакций связей от действия приложенных сил.

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

Статья научная