1503_Конструкция и расчет самолета на прочность крыло
.pdfДля расчета сечения крыла на изгиб вычерчивается профиль расчетного сечения крыла, на котором размещаются пронумерованные стрингеры и лонжероны (рис.12). В носике и хвостике профиля следует располагать стрингеры с большим шагом, чем между лонжеронами. Расчет сечения крыла на изгиб проводится методом редукционных коэффициентов и последовательных приближений.
y
|
|
y 1 |
|
xi |
Fi |
|
|
|
|
3 |
i |
1 |
2 |
x1 |
|
yi O |
|
O |
|
x n
Рис. 12
6.1. Порядок расчета первого приближения
Определяются в первом приближении приведенные площади поперечного сечения продольных ребер (стрингеров, поясов лонжеронов) с присоединенной обшивкой
Fiпр. = Fi Σ × j1i , Fi Σ = Fi + Fiобшприс. . , (9)
где F - действительная площадь сечения i-го ребра; F обш. |
- |
присоединенная |
||
i |
i прис. |
|
|
|
площадь обшивки ( F обш |
= 0,8 × d × ℓ- для растянутой панели, |
F обш |
= 30 × d2 - |
|
i прис. |
|
|
i прис. |
|
для сжатой панели); j1i - редукционный коэффициент первого приближения. Если материал полок лонжеронов и стрингеров разный, то следует сде-
лать приведение к одному материалу через редукционный коэффициент по модулю упругости
ji E = E i ,
Eo
где Ei - модуль материала i-го элемента; Eo - модуль материала, к которому приводится конструкция (как правило, это материал пояса самого нагруженного лонжерона). Тогда
j1i Σ = j1i × ji E
21
Вслучае разных материалов поясов лонжеронов и стрингеров в формулу
(9)вместо j1i подставляется j1i Σ .
Определяем координаты xi и yi центров тяжести сечений продольных элементов профиля относительно произвольно выбранных осей x и y (рис. 12)
и вычисляем статические моменты элементов Fiпр. × xi и Fiпр. × yi .
Определяем координаты центра тяжести сечения первого приближения по формулам
|
|
n |
|
|
|
|
∑ Fiпр. × xi |
|
|
x1 |
= |
i =1 |
, |
|
n |
||||
c |
|
|
||
|
|
∑ Fiпр. |
|
|
|
|
i =1 |
|
|
|
n |
|
|
|
|
∑ Fiпр. × yi |
|
|
y1 |
= |
i =1 |
. |
|
n |
||||
c |
|
|
||
|
|
∑ Fiпр. |
|
|
|
|
i =1 |
|
Через найденный центр тяжести проводим оси x1и y1 (ось x1 удобно выбрать параллельной хорде сечения) и определяем координаты центров тяжести всех элементов сечения относительно новых осей.
Вычисляем моменты инерции (осевые и центробежный) приведенного сечения относительно осей x1 и y1 :
n |
|
n |
I пр1 . |
n |
× yi1 . |
I пр1 . = ∑ Fiпр. × yi12 , |
I пр1 . = ∑ Fiпр. × xi12 , |
1 = ∑ Fiпр. × xi1 |
|||
x |
y |
i =1 |
x y |
i =1 |
|
i =1 |
|
|
|
Определяем угол поворота главных центральных осей сечения:
|
2 × I |
пр1 . |
1 |
|
tg 2a = |
|
x y |
|
. |
I пр1 . - I пр1 . |
||||
|
y |
|
x |
Если угол α будет больше 5о, то оси x1и y1 следует повернуть на этот угол (положительное значение угла соответствует вращению осей по часовой стрелке) и далее вести расчет относительно главных центральных осей. В целях упрощения расчета угол α рекомендуется вычислять только при расчетах последнего приближения. Обычно, если ось x1выбрана параллельно хорде сечения , угол α оказывается незначительным и им можно пренебречь.
Определяем напряжения в элементах сечения в первом приближении
si |
= |
M |
изг. |
× yi |
× ji × jE . |
1 |
|
1 |
1 |
||
|
|
I |
пр1 . |
|
|
|
|
|
x |
|
|
Полученные напряжения s1i сравниваем с σкр.л−н и σкр.стр. для сжатой панели и с σв.л−н и σв.стр. - для растянутой панели.
22
6.2. Определение критических напряжений стрингеров
Критическое напряжение стрингера σкр.стр. вычисляется из условия об-
щей и местной форм потери устойчивости. Для вычисления σкр.стр. общей фор-
мы потери устойчивости используем выражение |
|
||||||||||||
|
|
s |
|
= s |
|
× |
|
1 + ν |
, |
(10) |
|||
|
кр.стр. |
|
|
|
|
||||||||
|
ув |
|
|
|
в |
1 + n + n2 |
|
|
|||||
гден= |
. Здесь σкр. э - критическое напряжение, вычисленное по формуле |
||||||||||||
|
|||||||||||||
|
укр. э |
|
|
|
|
|
|
|
|
|
|
||
Эйлера: |
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
sкр. э |
= |
p2 |
× E |
, l = |
n × a |
, |
(11) |
|||
|
|
|
l2 |
|
|
ix |
|||||||
|
|
|
|
|
|
|
|
|
|
где μ - коэффициент, зависящий от условий опирания концов стрингера; a - шаг нервюр; λ - гибкость стрингера с присоединенной обшивкой; ix - радиус инер-
ции относительно центральной оси сечения.
В формуле (11) под ix следует понимать imin , но в целях упрощения по-
ложение главной инерциальной оси считаем совпадающим с осью x. В свою очередь
i x |
= |
I |
x |
, |
|
F пр. |
|||||
|
|
|
|||
|
|
|
стр. |
|
где I x - момент инерции стрингера с присоединенной обшивкой относительно оси x (рис.13); Fстрпр.. - площадь сечения стрингера с присоединенной обшивкой. Ширина присоединенной обшивки берется рав-
ной 30 δ (рис.13).
F пр. |
= F |
+ F прис. |
стр. |
стр. |
обш. |
При этом
I |
x |
= I |
x |
+ F |
прис. × a2 |
+ I |
x |
+ F |
× a |
2 |
, |
|
|
обш. 1 |
|
стр. |
|
2 |
|
||||
|
|
|
|
1 |
|
|
|
2 |
|
|
|
23
где I x1 - момент инерции присоединенной обшивки относительно собственной центральной оси x1 (обычно значения I x1 -малы); I x2 - момент инерции стрин-
гера относительно собственной центральной оси x2 .
Для вычисленияσкр.стр. местной формы потери устойчивости рассмотрим
потерю устойчивости свободной полки стрингера как пластины, шарнирно опертой по трем сторонам (рис.14). На рис. 14 обозначено: а – шаг нервюр; b1 – высота свободной полки стрингера (рис.13). Для рассматриваемой пластинки σкр. вычисляется по асимптотической формуле (10), в которой
sкр. э |
|
δ |
|
2 |
|
= 0,9 × kσ × E |
|
с |
|
, |
|
|
|
b1 |
|
|
где kσ – коэффициент, зависящий от условий нагружения и опирания пластины, dс – толщина свободной полки стрингера.
Для рассматриваемого случая
b |
|
2 |
||
kσ = 0,425 + |
1 |
. |
||
a |
||||
|
|
|
Для сравнения с действительными напряжениями, полученными в результате редуцирования, выбирается меньшее напряжение, найденное из расчетов общей и местной потери устойчивости.
В процессе редуцирования необходимо обратить внимание на следующее: если напряжения в сжатой полке лонжерона окажутся больше или равными разрушающим в любом из приближений, то конструкция крыла не способна выдержать расчетную нагрузку и ее надо усилить. Дальнейшие приближения в этом случае делать не следует. Если в каком-либо сжатом стрингере с номером
"k" (с присоединенной обшивкой) напряжение sik окажется меньше σкр.стр. k , то редукционный коэффициент для него и в последующем приближении следует оставить прежним jik+1 = jik ; если в каком-либо сжатом стрингере (с присоединенной обшивкой) с номером "m" напряжение s im окажется больше σ кр.стр. m то в последующем
приближении редукционный коэффициент следует вычислять по формуле
ji +1 = sкр.стр. m × yлi − на m × ymi ;
если ни в одном стрингере напряжение sij не пре-
высит σкр.стр. j , то конструкция явно перетяже-
лена и требует облегчения.
24
В растянутой зоне уточнение редукционных коэффициентов в процессе последовательных приближений ведется так же, но сравнение расчетных напряжений ведется не с σкр.стр. , а с σв .
В результате мы получаем новые уточненные редукционные коэффициенты последующего приближения jij+1 . Далее рассчитываем следующее приб-
лижение в том же порядке и снова уточняем редукционные коэффициенты. Расчет продолжается до тех пор, пока редукционные коэффициенты двух последующих приближений практически совпадут (в пределах 5%).
7. Расчет сечения крыла на сдвиг.
Расчет сечения крыла на сдвиг ведется без учета влияния кручения (поперечная сила QΣ считается приложенной в центре жесткости сечения, полагая, что на сдвиг работают стенки лонжеронов и обшивка).
7.1. Порядок расчета
Для расчета многоконтурного сечения на сдвиг делаются продольные разрезы в панелях таким образом, чтобы контур стал открытым. Для сечения крыла разрезы удобно делать в плоскости хорд в носке крыла и в стенках лонжеронов (рис. 15). В местах разрезов прикладываются неизвестные замыкающие погонные касательные усилия qI , qII , qIII .
Погонные касательные усилия qi в обшивке панелей сечения крыла определяются как сумма погонных касательных усилий qoi в незамкнутом контуре и замыкающих усилий qI , qII , qIII . Усилия qoi определяются формулой
|
|
Q |
× S |
отс. |
|
|
qoi |
= - |
Σ |
|
1, i −1 |
, |
(12) |
I |
|
|
||||
|
|
xред. |
|
гдеQΣ -расчетная перерезываю-
|
k =i−1 |
|
|
|
щая сила; S1отс. = |
∑ y |
k |
× F пр. |
- |
,i−1 |
k =1 |
k ред. |
|
|
|
|
|
|
статический момент площади части сечения, ограниченного 1-м и (i-1) – м ребрами (принятый порядок нумерации ребер очевиден из
рис. 14); I xред. - главный момент инерции всего сечения, причем положение центра тяжести берется
25
из последнего приближения расчета на изгиб.
В формуле (12) направление поперечной силы QΣ считается положительным при его совпадении с положительным направлением оси y, т.е. вверх. Положительные направления потоков касательных усилий совпадают с направлением обхода начала координат по часовой стрелке.
Для определения замыкающих потоков погонных касательных усилийqI , qII , qIII составляем канонические уравнения
[A]×{q}+ {Ao } = {0}
Коэффициенты канонических уравнений (элементы матрицы [A]и вектора {Ao }) определяются выражениями:
a11 |
= |
1 |
× ∑ |
li |
, a22 |
= |
1 |
× ∑ |
li |
, a33 |
= |
1 |
× ∑ |
li |
, |
|
Go |
di ред. |
Go |
di ред. |
Go |
di ред. |
|||||||||||
|
|
ABCDEFA |
|
|
BCDEFB |
|
|
CDEC |
|
(здесь суммирование ведется по панелям, где qI , qII , qIII не равны нулю соответственно),
a12 |
= a21 |
= |
1 |
× ∑ |
|
li |
|
|
, |
a23 |
= a32 |
= |
1 |
× ∑ |
li |
, |
|||
Go |
|
di ред. |
Go |
di ред. |
|||||||||||||||
|
|
|
BCDEF |
|
|
|
|
|
|
|
CDE |
|
|||||||
|
|
|
|
a13 |
= a31 |
= |
1 |
× ∑ |
|
li |
, |
|
|
|
|
||||
|
|
|
|
|
G |
|
d |
i ред. |
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
CDE |
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
o |
|
|
|
|
|
|
|
|
(здесь суммирование ведется по панелям, где qI , qII ; qII , qIII ; qI , qIII - соответственно не равны нулю),
a10 |
= |
1 |
× ∑ |
li × qoi |
, a20 |
= |
1 |
× ∑ |
li × qoi |
, a30 |
= |
1 |
× ∑ |
li × qoi |
, |
|
|
|
|||||||||||||
|
|
|
di ред. |
|
|
||||||||||
|
|
Go ABCDEFA |
di ред. |
|
Go BCDEF |
|
|
Go CDE di ред. |
(здесь суммирование ведется по панелям, где qI ,q0i ; qII , q0i ; qIII , q0i - соответственно не равны нулю). Здесь li -длина i-той панели; Go - приведенный модуль
сдвига (для обшивки из дюраля G = 105 |
кг × см− 2 ) ; d |
i ред. |
- |
редуцированная тол- |
o |
|
|
|
|
|
щина |
обшивки di ред. = di × Y ; |
||
|
Ψ - редукционный коэффи- |
|||
|
циент обшивки. |
|||
|
|
Модуль сдвига обшивки |
||
|
панели крыла не равен модулю |
|||
|
сдвига материала обшивки, а |
|||
|
зависит еще от ее кривизны, |
|||
|
толщины, шага нервюр и |
|||
|
стрингеров (размеров подкреп- |
|||
|
ляющей |
|
клетки), подкрепля- |
ющих профилей, характера нагружения пластины. Значения модуля сдвига бо-
26
лее или менее точно определяются опытным путем для данной конструкции. В расчете приходится большей частью пользоваться средними величинами G, полученными из испытаний аналогичных конструкций. Так как
Y = G ,
Go
то при вычислении мы будем пользоваться значениями редукционных коэффициентов, приведенными на рис. 15. Значения коэффициента Ψ для обшивки из другого материала следует умножить на ϕE .
Определяем погонные сдвиги в панелях как соответствующую сумму погонных касательных усилий в открытом контуре и замыкающих интенсивностей:
qΣсдв. = qoi + qI + qII + qIII
По результатам расчета строим схему потоков погонных касательных усилий по контуру сечения.
8.Расчет сечения крыла на кручение
8.1.Определение положения центра жесткости сечения крыла
Положение центра жесткости определяется по формуле
|
2 |
n |
|
|
xцуточн.ж. . = |
× (qI × FI + qII × FII + qIII × FIII + ∑qoi |
× wi ), |
||
QΣ |
||||
|
i =1 |
|
где FI , FII , FIII - площади контуров ABCDEFA, BCDEFB, CDEC соответственно
ℓi |
Уточненное |
QΣ положение ц.ж. |
i-1 i
ri
|
|
|
|
|
ц. ж. |
A |
|
x |
уточ. |
|
x |
|
|
||||
|
ц.ж. |
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
Рис. 16
(рис. 15), которые подсчитываются по чертежу сечения крыла; qI , qII , qIII - потоки погонных касательных усилий, полученные в результате расчета сечения на сдвиг от силы QΣ ; qoi - потоки погонных касательных усилий в открытом контуре сечения крыла от сдвига; ωi - секториальная площадь, соответствующая i – той панели (рис. 16). Приближенно значения ωi
можно вычислить как площадь треугольника
27
wi » 1 × ℓi × ri ,
2
где ri - длина перпендикуляра, опущенного из произвольно выбранного полюса А (в качестве такого полюса можно взять координату предварительно определенного центра жесткости) на хорду дуги элемента контура сечения ℓi .
|
|
8.2. Определение потока касательных усилий от кручения |
|
|
При |
расчете на кручение замыкающие потоки |
касательных усилий |
qкр. |
, qкр. , q |
кр. определяются из системы уравнений |
|
I |
II |
III |
|
|
|
[ A]×{q} = 2 ×x{F} , |
|
|
|
2 ×{F}т ×{q}= M кр.уточн. |
(13) |
|
|
|
|
Здесь компоненты вектора {F} - FI , FII , FIII - площади |
контуров ABCDEFA, |
||
BCDEFB, CDEC (см. рис. 15); ξ - относительный угол закручивания сечения; |
|||
M кр.уточн. - крутящий момент относительно уточненного положения центра же- |
|||
сткости сечения, определяемый из выражения |
|
||
|
|
M кр.уточн. = M кр. ± QΣ × d |
|
Здесь d – |
расстояние между приближенным и уточненным положениями цен- |
||
тров жесткости сечения крыла. |
|
||
|
Значения коэффициентов aij канонических уравнений (13) те же, что при |
расчете на сдвиг. После определения потоков замыкающих касательных усилий qIкр. , qIIкр. , qIIIкр. при кручении, суммарные погонные сдвиговые усилия находим подобно расчету на сдвиг, положив qoi = 0 , т.е.
qiкрΣ . = qIкр. + qIIкр. + qIIIкр.
По результатам расчета строим суммарную эпюру потоков погонных ка-
сательных усилий от сдвига и кручения qiсдвΣ + кр по контуру расчетного сечения крыла. При построении суммарной эпюры положительные значения потоков откладываем внутрь контура сечения.
28
9.Проверка обшивки и стенок лонжеронов на прочность
иустойчивость
Врезультате проверочного расчета должно быть дано заключение о прочности подобранного сечения крыла. Для этого обшивка и стенки лонжеронов проверяются на прочность и устойчивость.
Максимальные нормальные напряжения, действующие на соответствую-
щую панель обшивки (или стенки лонжерона) с учетом ϕE
smax = s |
× |
E |
|
обш |
, |
||
|
|||
обш |
i max |
E |
|
|
|
||
|
|
стр |
а значения редукционного коэффициента обшивки находятся по выражению
ϕобш = bприс. ,
b
где bприс . = 30 × d - для сжатой зоны; bприс. = 0,8 × ℓ- для растянутой зоны; ℓ - ширина рассматриваемой панели обшивки (шаг стрингеров). Тогда средние нормальные напряжения в панелях обшивки
sср. обш. = jобш × smax . обш
Касательные напряжения, действующие в обшивке (или стенке лонжерона) от сдвига и кручения, вычисляются как
tiΣ |
= |
qiсдвΣ |
+ кр |
|
|
|
. |
||
|
|
|||
|
|
di |
Критические касательные напряжения tкр вычисляются аналогично sкр по
формуле |
|
|
|
|
|
|
|
|
|
ℓ |
σср.обш. |
|
1 + n* |
|
|
|
|
|
|
tкр |
|
|
|
|||
|
|
|
= tв ×1 + n* + n* 2 , |
|
||||
τiΣ |
|
|
где tв |
= 0,7 × sв , n* |
= |
τв |
, |
|
|
a |
|
|
|
|
tкр.э |
|
|
|
|
|
tкр.э |
d |
|
2 |
|
|
|
|
|
= 0,9 × E × kτ × |
|
i . |
|
||
|
|
|
|
|
b |
|
||
|
|
|
Значения коэффициента kτ берется из ра- |
|||||
|
|
|
боты [8] для пластины шарнирно опертой по |
|||||
|
Рис. 17 |
|
контуру. |
|
|
|
|
|
|
|
Для оценки устойчивости элементов кры- |
||||||
|
|
|
||||||
|
|
|
29 |
|
|
|
|
|
ла вычисляется коэффициент
|
|
sср.обш. |
|
t |
iΣ |
2 |
||
k |
|
= |
|
|
+ |
|
. |
|
|
|
|
|
|||||
|
кр |
|
sкр |
|
|
|
|
|
|
|
|
|
tкр |
Значения коэффициента kкр позволяют судить о работе обшивки (стенки
лонжерона) на устойчивость при одновременном действии сжатия и сдвига (рис. 17). Растянутая панель находится в этом случае в облегченных условиях и в данном случае не рассматривается. Если kкр < 1, то потери устойчивости не
произойдет, при kкр > 1 пластинка потеряет устойчивость. Потеря устойчивости
не является критерием потери несущей способности конструкции.
При проверке обшивки на прочность вычисляются значения коэффициента kпр с использованием четвертой теории прочности:
kпр |
= |
s2max + 3 × ti2Σ |
|
обш |
|||
sв |
|||
|
|
£ 1 , |
(14) |
где s2max + 3 × ti2Σ = sэквIV . .
обш
Для стенок лонжеронов σmax = 0 (чистый сдвиг) коэффициент kпр вычис-
ляется по формуле (14) до потери устойчивости и, если стенка потеряла устойчивость, то по выражению
kпр = |
2 × tΣ |
£ 1. |
|
sв |
|||
|
|
Значения коэффициента kпр < 1 позволяют сделать вывод о том, что условие прочности соблюдается.
Ли т е р а т у р а
1.Бадягин А.А., Егер С.М. и др. Проектирование самолетов. – М.: Машиностроение, 1972. 516с.
2.Зайцев В.Н. Рудаков В.Л. Конструкция и прочность самолетов. – Киев: Вища школа, 1976. 400с.
3.Шульженко М.Н. Конструкция самолетов. – М.: Машиностроение, 1971.
4.Кравец А.С. Характеристики авиационных профилей. – М.: Оборонгиз, 1939.
5.Макаревский А.И., Корчемкин Н.Н., Француз Т.А., Чижов В.М. Прочность самолета. – М.: Машиностроение, 1975. 280с.
30