CC BY
37
Dual Matrix and Biquaternion Methods of Solving Direct and Inverse Kinematics Problems of Manipulators for Example Stanford Robot Arm. II
E. I. Lomovtseva, Yu. N. Chelnokov
Saratov State University, 83, Astrahanskaya str., 410012, Saratov, Russia, [email protected], [email protected]
The methology of solving the inverce kinematics problem of manipulators by using biquaternion theory of kinematics control is shown on the example of Stanford robot arm. Solving of the inverce kinematics problem of Stanford robot arm is performed using the simplest control law. The analysis of numerical solution results is made. The efficacy of applying the theory of kinematics control for solving the inverce kinematics problem of manipulators is proved. Dual matrix and biquaternion methods of solving direct kinematics problem of manipulators were considered in [1].
Keywords: robot-manipulator, direct kinematics problem, dual direction cosine matrix, biquaternion, quaternion, kinematic equations.
References
1. Lomovtseva E. I., Chelnokov Yu. N. Dual Matrix and Biquaternion Methods of Solving Direct and Inverse Kinematics Problems of Manipulators, for Example Stanford Robot Arm. I. Izv. Saratov. Univ. (N. S.), Ser. Math. Mech. Inform., 2013, vol. 13, no. 4. pp. 82-89 (in Russian).
2. Chelnokov Yu. N. Biquaternion Solution of the Kinematic Control Problem for the Motion of a Rigid Body and Its Application to the Solution of Inverse Problems of Robot-Manipulator Kinematics. Mechanics of
УДК 539.3
Solids [Izv. RAN. Mehanika tverdogo tela], 2013, vol. 48, no. 1. pp. 31-46.
3. Fu K. S., Gonzalez R. C. ,Lee C. S. G. Robotics : Control, Sensing, Vision, and Intelligence. McGraw-Hill, Inc, 1987, 580 p.
4. Chelnokov Yu. N. Kvaternionnye i bikvaternionnye modeli i metody mehaniki tverdogo tela i ih prilozhenija. Geometrija dvizhenija [Quaternion and Biquaternion Models and Methods of Mechanics of a Rigid Body and their Applications. Geometry of Motion.] Saratov, Saratov Univ. Press, 2006, 236 p. (in Russian).
УСТОЙЧИВОСТЬ КОНСТРУКТИВНО-ОРТОТРОПНОЙ НЕОДНОРОДНОЙ ЦИЛИНДРИЧЕСКОЙ ОБОЛОЧКИ ОТ НЕРАВНОМЕРНОЙ РАДИАЛЬНОЙ НАГРУЗКИ
А. А. Мочалин
Кандидат физико-математических наук, доцент кафедры прикладной математики и системного анализа, Cаратовский государственный технический университет имени Гагарина Ю. А., [email protected]
На базе полубезмоментной теории В. З. Власова рассматривается задача об устойчивости цилиндрической конструктивно-ортотропной оболочки переменной вдоль образующей толщины при действии осесимметричного изменяющегося вдоль оси оболочки радиального давления. При одном соотношении изменения толщины и давления получено точное решение для нахождения одной из величин в законе изменения давления, при которой происходит потеря устойчивости оболочки.
Ключевые слова: цилиндрическая оболочка, теория оболочек, устойчивость оболочек, радиальная нагрузка, критическое давление, толщина оболочка, теория упругости.
1. В линейной теории цилиндрических оболочек широкое применение нашла полубезмоментная теория В. З. Власова, учитывающая особенности напряженного состояния оболочек, длина которых находится в пределах Б < I < 8Б, где Б — диаметр оболочки, I — её длина. В основе этой теории лежат две гипотезы — статическая и геометрическая, позволяющие существенно упростить уравнения, описывающие состояние устойчивости оболочки. Считается, что удлинение в окружном направлении и углы сдвига в срединной поверхности равны нулю, полагаются равными нулю перерезывающая сила и изгибающий момент в осевом направлении, а также крутящий момент.
На основе гипотез полубезмоментной теории цилиндрических оболочек В. З. Власова [1] рассмотрим конструктивно-ортотропную неоднородную оболочку под действием неравномерной радиальной
© Мочалин Л. Л., 2014
95
Изв. Сарат. ун-та. Нов. сер. Сер. Математика. Механика. Информатика. 2014. Т. 14, вып. 1
нагрузки. Уравнение совместности деформаций и уравнение равновесия после несложных преобразований перепишется в виде
Rd4ei + 1 д2ei _ ду4 + R ду2 _
d4 Mo 1 д2 M2 , R——;—г — ^ 0—+
д2
К2
дх2 д 2 Ti
ду4
R ду2
дх2
= P,
(1) (2)
р = _ дз! + дз! _
дх ду ду2
где х — координата, откладываемая по образующей, у — координата, откладываемая по дуге поперечного круга цилиндра.
Соотношения упругости с учетом допущений полубезмоментной теории В. З. Власова запишутся в форме
Т ^ ЕНтхе1, М2 = 5П2П2к2.
Упругие усилия и моменты связаны с деформациями срединной поверхности соотношениями (2).
Введем безразмерные переменные С = х/1, в = у/Я (1 — длина оболочки) и функцию Ф, связанную с компонентами поля перемещений и (и, V, т) равенствами
u = —
1 д3Ф / дС дв2'
v =
i д2Ф
R дв2'
w =
i д4 ф R дв4"
Тогда
£1 =
Ti _
i д4 ф
I2 дС2дв2' Eh д 4Ф
Y дС2дв2'
= 1 /д4Ф д6 Ф
К2 = - R3 V дв4 + дв6
D ( д6 Ф д4 Ф M2 =----1--
M2 R3 V дв6 + дв4
(3)
Подстановка (3) тождественно удовлетворяет уравнению (1), а уравнение (2) перепишется в сле дующем виде:
д^ дС2
д2Ф
МС) w
D/4 /4
+^ E
в в
P * dtfdtf = 0.
(4)
о о
д2
1 +
дифференциальный оператор
о 1 1 / д^2 д^„\ д2
Здесь P* = -ТаС + Пае - «V• а п = №
В. З. Власова.
Уравнением (4) будем пользоваться при исследовании устойчивости цилиндрических оболочек средней и большой длины под действием радиальной нагрузки и интегрировать при граничных условиях, зависящих от способа закрепления торцов оболочки.
2. Пусть оболочка находится под действием неравномерной радиальной нагрузки q = q0(1 + аС)-6. Компонента нагрузки qn в этом случае определится выражением qn = —q0RK2(1 + аС)-6, где параметр q0 подлежит определению.
Выбирая функцию Ф в форме Ф = X(С) sin кв, где k — число волн вдоль окружности оболочки при потере устойчивости, и задавая закон изменения толщины в виде Л,(С) = h0(1 + аС)-2 , приходим к уравнению
¿С2
(1 + аС)
-2
d2 X
¿С2
где
Л4 = q0k4/4 (k2 _ 1) _
Л = ER3h0 (k 1)
— Л4(1 + аС)-6 X = 0, /4h0 (k2 — 1)k4
(5)
12R6 (1 — v2)
Уравнение (5) после введения новой переменной г = 1п(1 + а£) и новой функции Ш(г) соотношением X = Шв-2/6 может быть приведено к дифференциальному уравнению четвертого порядка с
2
d
постоянными коэффициентами:
^ (4) -10 ^ (3) + 25 ^ (2) + 89 ^ -(Л4 - ^ ^ = о
3 6 54 V«4 1296 У '
характеристическое уравнение для которого запишется в виде
(52 + о5 + с)(52 + Ь5 + й) = О,
где величины о, Ь, с, й связаны равенствами
0 + Ь = -1°, оЬ + с + й =25, Ьс + йо = |, сй =±1 - *
3 6 54 1296 а4
и равны
7 5 25а2 + /6ООа4 + 1296Л4 , 25а2 ± ^6о0а4ТТ296Л4
о = Ь = —, с =---, а =---.
3 36а2 36а2
Выражение для функции X(С) примет вид
X(С) = (1 + аС)2/3(Сх йп(5 ■ 1п(1 + аС)) + С2 еов(5 ■ 1п(1 + аС))+ +С3 бЬ(5 ■ 1п(1 + аС)) + С4 еЬ(5 ■ 1п(1 + аС)),
Р / 6°°а4 + 1296Л4
где 5 =---.
6а
Рассмотрим два вида граничных условии: жесткое закрепление и свободное опирание краёв оболочки.
В случае жесткого закрепления краёв граничные условия имеют вид
и = О, у = О при С = О, С = 1
или
х (С) = о, аХСС) = ° при С = °, С = 1.
Если края оболочки свободно оперты, то следует удовлетворять условиям вида
и = О, Тх = О при С = О, С = 1
или
X (С) = О, = О при С = О, С = 1.
а2Х(С)
йС2
Подчиняя общее решение уравнения граничным условиям и требуя совместности системы однородных уравнении, получаем соотношения для определения собственных значении для жесткого закрепления в форме
51п(1 + а) =4.73, а в случае свободного опирания краев — в форме
Бт(51п(1 + а)) = О.
Для очень длинных оболочек (Я/1 ^ О) величина до нечувствительна к виду закрепления и равна
= (к2 - 1) до = Л312(1 - V2) •
После минимизации по к2 (к2 ^ 1) для оболочек средней длины получаем формулу для вычисления критического давления:
*
Iт (6)
Механика
97
Изв. Сарат. ун-та. Нов. сер. Сер. Математика. Механика. Информатика. 2014. Т. 14, вып. 1
где величина во в случае жесткого закрепления равна
во =
а в случае свободного опирания краев —
во =
4.73а ln(1 + а)
75 162
а ,
па
ln(1 + а) J 162
75
а,
д* = д/п.
Формулу (6) можно переписать для жесткого закрепления краев в форме
q* = 4.73ЕК (а)(Но/R)5/2 R ^
4 п2
l 9(1 - д*2)3/4 \/ mi' a для свободного опирания торцов соответственно в виде
,R V6
= пЕК (а)(Но/R)5/2-
42
I 9(1 - д*2)3/4 у шГ Значения коэффициента К(а) для некоторых значений а приведены в таблице.
(7)
(8)
4
4
3
a -0.5 -0.4 -0.3 -0.2 -0.1 0 0.1 0.2 0.3 0.4 0.5
K(a) 0.750 0.798 0.847 0.888 0.950 1 1.049 1.093 1.136 1.173 1.20
На основании полученных результатов можно сделать следующий вывод.
Значения критического давления неоднородной цилиндрической оболочки под действием неравномерного давления, изменяющегося по закону изменения жесткости, пропорциональны критическому давлению оболочки постоянной толщины с коэффициентом пропорциональности К (а), и при увеличении отношения толщин Н\/Н0 на концах оболочки критическое давление уменьшается, а при его уменьшении критическое давление увеличивается.
Библиографический список
1. Власов В. З. Общая теория оболочек. М. : ГИТТЛ, 1949. 784 с.
2. Бурмистров Е. Ф. Симметричная деформация кон-структивно-ортотропных оболочек вращения. Саратов : Изд-во Сарат. ун-та, 1962. 108 с.
3. Мочалин А. А. Устойчивость полубезмоментной ци-
линдрической оболочки переменной толщины // Изв. вузов. Машиностроение. 1975. № 11. С. 27-31. 4. Штейнберг М. В. Расчет круговых цилиндрических оболочек с толщиной переменной в направлении образующей // Прикладная механика. 1965. Т. 1, вып. 7. С. 67-72.
The Stability of the Constructive-orthotropic Heterogeneous Cylindrical Shell
under Uneven Radial Load
A. A. Mochalin
Saratov State Technical University, 77, Polytechnicheskaya str., 410054, Saratov, Russia, [email protected]
On the base haft-momentum Vlasov theory the problem of stability of cylindrical homogeneas shell with variation of thicknees atv radial symmetrical ractial pressure variated onalong axe distance. At one reletion between thickness and pressure values the accurate solution was produced for one values in pressure variation law when stability of shell is sailed.
Key words: cylindrical shell, theory of shells, stability of shells, radial loading, critical pressure, shell thickness, the theory of elasticity.
Е. В. Мурашкин. Расчет параметров нагружения полого шара
References
1. Vlasov V. Z. Obshchaia teoriia obolochek [General theory of shells]. Moscow, Gos. izd. teh. teor. lit., 1949, 784 p. (in Russian).
2. Burmistrov E. F. Simmetrichnaia deformatsiia kon-struktivno-ortotropnykh obolochek vrashcheniia [Symmetric deformation of structurally orthotropic shells of revolution]. Saratov, Saratov. Univ. Press, 1962. 108 p. (in Russian).
УДК 539.274
3. Mochalin A. A. Sustainability semi-momentless cylindrical shell of variable thickness. Izv. vuzov. Mashino-stroenie, 1975, no. 11, pp. 27-31 (in Russian).
4. Shteinberg M. V. Calculation of circular cylindrical shells with variable thickness in the direction of forming. Prikladnaia mekhanika, 1965, vol. 1, iss. 7, pp. 67-72 (in Russian).
РАСЧЕТ ПАРАМЕТРОВ НАГРУЖЕНИЯ ПОЛОГО ШАРА В УСЛОВИЯХ БОЛЬШИХ УПРУГОПОЛЗУЧИХ ДЕФОРМАЦИЙ
Е. В. Мурашкин
Кандидат физико-математических наук, научный сотрудник, Институт проблем механики им. А. Ю. Ишлинского РАН, Москва, [email protected], [email protected]
Предложена модель больших упругоползучих деформаций. Разделение тензора полных деформаций альманси определяется уравнениями изменения обратимой и необратимой его составляющих. Рассмотрено сферически симметричное деформирование полого шара в процессе установившейся ползучести. Получена разрешающая система уравнений рассматриваемой краевой задачи. Предложен способ определения нагружающего усилия вызывающего заданное деформированное состояние. По заданным законам изменения поля перемещений построены функции внешнего нагружающего усилия.
Ключевые слова: большие деформации, ползучесть, упругость, релаксация. ВВЕДЕНИЕ
Необходимость повышения точности математического описания процессов, происходящих при технологической обработке и эксплуатации металлоизделий, вынуждает учитывать упругие свойства материалов на всех стадиях жизненного цикла изделия. Рассмотрение задач в классических моделях малых деформаций невозможно, когда относительное изменение формы рассматриваемого тела велико. Одной из таких характерных задач, где нельзя обойтись без применения модели больших деформаций, является задача о моделировании процессов в окрестности микропоры в металле, происходящих под действием интенсивного давления. Актуальность данной задачи обусловлена обнаруженным на опыте эффектом существенного повышения эксплуатационных характеристик металла при интенсивном всестороннем сжатии образцов [1] «залечивания» микродефектов сплошности. Попытки смоделировать процесс залечивания микропоры в металле делались неоднократно, в том числе и на основе модели больших упругопластических деформаций [2], обладающей эффектом приспосабливаемости к периодическим нагружениям по циклу «нагрузка - разгрузка» [3].
В настоящей работе решена задача о сферически симметричном сжатии шара с микропорой в центре. Условие несжимаемости среды определяет кинематику среды с точностью до неизвестной функции времени, что позволяет по известному закону деформирования определить процесс нагру-жения, вызывающий заданное деформированное состояние.
1. ОСНОВНЫЕ МОДЕЛЬНЫЕ ЗАВИСИМОСТИ
За основу возьмем модель больших упругопластических деформаций [3], основные кинематические соотношения которой в прямоугольной декартовой системе пространственных (эйлеровых)
© Мурашкин Е. В., 2014
99
