4.4: Ламіновані композитні плити
Вступ
Цей документ призначений для окреслення механіки армованих волокном ламінованих плит, що призводить до обчислювальної схеми, яка пов'язує внутрішньоплощинну деформацію і кривизну ламінату з накладеними на нього тягами і згинальними моментами. Хоча це невелика частина загальної області композитів, армованих волокном, або навіть теорії ламінату, це важлива техніка, яку повинні розуміти всі інженери композитів.
У наступних розділах ми розглянемо конституційні відносини для ізотропних матеріалів у матричній формі, а потім покажемо, що розширення до поперечно-ізотропних композитних пластин дуже просте. Оскільки кожен шар в ламінаті може бути орієнтований довільно, то ми потім покажемо, як еластичні властивості окремих ламінат можуть трансформуватися в загальний напрямок. Нарешті, ми збалансуємо напруження окремих шарів проти застосованих тяг і моментів, щоб розробити матричні відносини, що регулюють для ламінату в цілому.
Розрахунки для механіки ламінату найкраще проводити за допомогою комп'ютера, а алгоритми окреслені для еластичних ламінатів, ламінатів, що демонструють ефекти теплового розширення, та ламінатів, що демонструють в'язкопружну реакцію.
Ізотропні лінійні еластичні матеріали
Як показано в елементарних текстах з механіки матеріалів (пор. Roylance 1996 (див. Посилання, наведені в кінці цього документа)), декартові деформації, що виникають внаслідок стану плоского напруження (σz=τxz=τyz=0), є
γxy=1Gτxy
У плоскому напруженні також є деформація вz напрямку через ефект Пуассона:ϵz=−ν(σx+σy); цей компонент деформації буде ігноруватися в секціях, які слід дотримуватися. У вищезазначених співвідношеннях є три константи пружності: модуль ЮнгаEν, коефіцієнт Пуассона та модуль зсувуG. Однак для ізотропних матеріалів існує лише дві незалежні пружні константи, і, наприклад,G можуть бути отримані зE іν як
G=E2(1+ν)
Використовуючи матричні позначення, ці відносини можуть бути записані як
{ϵxϵxγxy}=[1/E−ν/E0−ν/E1/E0001/G]{σxσyτxy}
Величина в дужках називається матрицею відповідності матеріалу, позначаєтьсяS абоSij. Важливо зрозуміти фізичну значимість різних його термінів. Безпосередньо з правил множення матриці елемент вith рядку іjth стовпціSij - це внесокjth напруги вith деформацію. Наприклад, компонент у положенні 1,2 є внеском напруги -напрямку до деформаціїy -напрямку:σy множення на1/E дає деформаціюx -direction, породженуσy, а потім множення цього на−ν дає деформацію Пуассонаy індукований вx напрямку. Нульові елементи показують відсутність зчеплення між нормальними та зсувними компонентами.
Якщо ми хочемо записати напруження через деформації, рівняння 4.4.1 можна інвертувати, щоб дати:
{σxσyτxy}=E1−ν2[1ν0ν1000(1−ν)/2]{ϵxϵyγxy}
де тутG був заміненийE/2(1+ν). Це відношення може бути скорочено далі як:
σ=Dϵ
деD=S−1 - матриця жорсткості. Відзначимо, що модуль Юнга можна відновити, приймаючи зворотний елемент 1,1 матриці відповідностіS, але що положення 1,1 матриці жорсткостіD містить ефекти Пуассона і не дорівнюєE.
анізотропні матеріали
Якщо матеріал має текстуру, таку як дерево або однонаправлено армовані волокнисті композити, як показано на малюнку 1, модуль пружностіE1 у напрямку волокна, як правило, буде більшим, ніж у поперечних напрямках (E2іE3). КолиE1≠E2≠E3, матеріал, як кажуть, ортотропний. Однак загальним є те, що властивості в площині, поперечній до напрямку волокна, щоб бути ізотропними до хорошого наближення (E2=E3); такий матеріал називається поперечно-ізотропним. Пружні конститутивні закони повинні бути змінені з урахуванням цієї анізотропії, а наступна форма є продовженням звичайних рівнянь ізотропної пружності до поперечно-ізотропних матеріалів:
{ϵ1ϵ2γ12}=[1/E1−ν21/E20−ν12/E11/E20001/G12]{σ1σ2τ12}

Параметрν12 є основним коефіцієнтом Пуассона; це відношення деформації, індукованої у напрямку 2 деформацією, застосованою в 1-напрямку. Цей параметр не обмежується значеннями менше 0,5, як в ізотропних матеріалах. І навпаки,ν21 дає деформацію, індуковану в 1-напрямку деформацією, застосованою у напрямку 2. Оскільки 2-напрямок (поперечний до волокон), як правило, має набагато меншу жорсткість, ніж 1-напрямок, дана деформація в напрямку 1 зазвичай розвиває набагато більшу деформацію в напрямку 2, ніж така ж деформація в напрямку 2 індукує деформацію в напрямку 1. Отже, ми, як правило, матимемоν12>ν21. У вищенаведеному рівнянні є п'ять констант (E1,E2,ν12,ν21іG12. Однак тільки чотири з них незалежні; оскількиS матриця симетрична, ми маємоν21/E2=ν12/E1.
Проста форма Рівняння 4.4.4, з нулями в термінами, що представляють зв'язок між нормальною та зсувною складовими, виходить лише тоді, коли осі вирівняні по основним напрямкам матеріалу; тобто уздовж і поперечно до осей волокон. Якщо осі орієнтовані по якомусь іншому напрямку, всі терміни матриці відповідності будуть заповнені, і симетрія матеріалу не буде кидатися в очі. Якщо, наприклад, напрямок волокна знаходиться поза віссю від напрямку завантаження, матеріал буде розвивати деформацію зсуву, коли волокна намагаються зорієнтуватися вздовж напрямку завантаження. Тому буде зв'язок між нормальним напруженням і деформацією зсуву, яка не відбувається в ізотропному матеріалі.
Трансформація осей
Важливо вміти перетворювати осі в «лабораторний»x−y каркас і з нього на каркас з натурального матеріалу, в якому осі можуть бути марковані 1 - 2, що відповідають волокну та поперечним напрямкам, як показано на малюнку 2.

Як показано в елементарних підручниках, закон трансформації для декартового наголосу Коші можна записати:
σ1=σxcos2θ+σysin2θ+2τxysinθcosθσ2=σxsin2θ+σycos2θ−2τxysinθcosθτ12=(σy−σx)sinθcosθ+τxy(cos2θ−sin2θ)
Деθ - кут відx осі до осі 1 (волокна). Ці відносини можуть бути записані в матричному вигляді як
{σ1σ2τ12}=E1−ν2[c2s22scs2c2−2sc−scscc2−s2]{σxσyτxy}
деc=cosθ іs=sinθ. Це може бути скорочено як
σ′=Aσ
деA матриця перетворення в дужках вище. Цей вираз може бути застосований як до тривимірних, так і двовимірних напружених станів, хоча конкретна форма,A наведена в рівнянні 4.4.6, дійсна лише у двох вимірах (плоска напруга) та для декартових координат.
Використовуючи математичні або геометричні аргументи, можна показати, що компоненти нескінченно малої деформації трансформуються майже однаковими співвідношеннями:
{ϵ1ϵ212γ12}=A{ϵxϵy12γxy}
Коефіцієнт 1/2 на складових зсуву виникає з класичного визначення деформації зсуву, яка вдвічі перевищує тензоріальну деформацію зсуву. Це вносить деяку незручність в трансформаційні відносини, які можна зменшити, вводячи матрицю Рейтера, визначену як
[R]=[100010002] or [R]−1=[1000100012]
Тепер ми можемо написати:
{ϵ1ϵ2γ12}=R{ϵ1ϵ212γ12}=RA{ϵxϵy12γxy}=RAR−1{ϵxϵyγxy}
або
ϵ′=RAR−1ϵ
Закон трансформації для комплаєнсу тепер може бути розроблений з законів трансформації для деформацій і стресів. Послідовними перетвореннями деформація вx−y довільному напрямку пов'язана з деформацією в 1-2 (основний матеріал) напрямках, потім до напружень в 1-2 напрямках і, нарешті, до напружень вx−y напрямках. Остаточною групуванням матриць перетворення, що стосуютьсяx−y деформацій доx−y напружень, є перетворена матриця відповідності вx−y напрямку:
{ϵxϵyγxy}=R{ϵxϵy12γxy}=RA−1{ϵ1ϵ212γ12}=RA−1R−1{ϵ1ϵ2γ12}
=RA−1R−1S{σ1σ2τ12}=RA−1R−1SA{σxσyτxy}≡ˉS{σxσyτxy}
деˉS - перетворена матриця відповідності щодоx−y осей. Зворотний відˉS isˉD, матриця жорсткості щодоx−y осей:
ˉS=RA−1R−1SA, ˉD=ˉS−1
Приклад4.4.1
Розглянемо шар кевлар-епоксидного композиту з жорсткістюE1=82,E2=4,G12=2.8 (всі ГПа) іν12=0.25. орієнтований на30∘ відx осі. Жорсткість уx напрямку можна знайти як зворотну елементу 1,1 перетвореної матриці відповідностіS, як дано рівнянням 4.4.11. Нижче показано, як це можна зробити за допомогою символічного математичного програмного забезпечення Maple (відредаговано для стислості):
Read linear algebra package > with(linalg): Define compliance matrix > S:=matrix(3,3,[[1/E[1],-nu[21]/E[2],0],[-nu[12]/E[1],1/E[2],0],[0,0,1/G[12]]]); Numerical parameters for Kevlar-epoxy > Digits:=4;unprotect(E);E[1]:=82e9;E[2]:=4e9;G[12]:=2.8e9;nu[12]:=.25; nu[21]:=nu[12]*E[2]/E[1]; Compliance matrix evaluated > S2:=map(eval,S);
S2:=[.122010−10−.305010−110−.304910−11.250010−9000.357110−9]
Transformation matrix > A:=matrix(3,3,[[c^2,s^2,2*s*c],[s^2,c^2,-2*s*c],[-s*c,s*c,c^2-s^2]]); Trigonometric relations and angle > s:=sin(theta);c:=cos(theta);theta:=30*Pi/180; Transformation matrix evaluated > A2:=evalf(map(eval,A));
A2:=[.7500.2500.8660.2500.7500−.8660−.4330.4330.5000]
Reuter's matrix > R:=matrix(3,3,[[1,0,0],[0,1,0],[0,0,2]]); Transformed compliance matrix > Sbar:=evalf(evalm( R &* inverse(A2) &* inverse(R) &* S2 &* A2 ));
Sbar:=[.882810−10−.196810−10−.122210−9−.196910−10.207110−9−.837010−10−.122210−9−.837710−10.290510−9]
Stiffness in x-direction > 'E[x]'=1/Sbar[1,1];
Ex=.11331011
Зверніть увагу, що перетворена матриця відповідності симетрична (до числової помилки округлення), але існують ненульові значення зв'язку. Користувач, який не знає про внутрішній склад матеріалу, вважав би його повністю анізотропним.
Ламіновані композитні плити
Однією з найбільш поширених форм армованих волокном композитних матеріалів є перехресний ламінат, в якому виробник «закладає» послідовність односпрямованих армованих «шарів», як зазначено на малюнку 3. Кожен шар, як правило, являє собою тонкий (приблизно 0,2 мм) лист колімованих волокон, просочений неотвержденной епоксидною смолою або іншим термореактивним полімерним матричним матеріалом. Орієнтація кожного шару довільна, а послідовність укладання підібрана для досягнення бажаних властивостей ламінату. У цьому розділі ми описуємо, як розробляються та аналізуються такі ламінати.

\ Класична теорія ламінату» є продовженням теорії згинання однорідних пластин, але з урахуванням внутрішньоплоских тяг на додаток до згинальних моментів, і для різної жорсткості кожного шару при аналізі. У загальних випадках визначення тяг і моментів в заданому місці зажадає вирішення загальних рівнянь для рівноваги і зміщення сумісності пластин. Ця теорія розглядається в ряді стандартних текстів (пор. С.Тимошенко і С.Войновський-Крігер, Теорія плит і раковин, Макгроу-Хілл, Нью-Йорк, 1959 р.), і тут обговорюватися не буде.
Почнемо з того, що припускаємо знання тяг N і моментів M, застосованих до пластини в положенніx,y, як показано на малюнку 4:
N={NxNyNxy} M={MxMyMxy}

Зручно буде нормалізувати ці тяги і моменти по ширині плити, так вони мають одиниціN/m іN−m/m, або простоN відповідно. Координатиx іy є напрямками в площині плити, іz прийнято приймати як позитивні вниз. Прогин вz напрямку називаєтьсяw, також приймається як позитивний вниз.

Аналогічно з припущенням Ейлера для балок, припущення Кіршкофа для згинання пластини приймає спочатку прямі вертикальні лінії, щоб залишатися прямими, але обертатися навколо середньої площини (z=0). Як показано на малюнку 5, горизонтальні зсувиu іv вy напрямкахx і за рахунок обертання можуть бути прийняті до розумного наближення від кута повороту і відстані від середньої площини, і це обертальне переміщення додається до зміщення середньої площини ( u0,v0):
u=u0−zw0,x
v=v0−zw0,y
ϵ={ϵxϵyγxy}={u,xv,yu,y+v,x}={u0,x−zw0,xxv0,y−zw0,yy(u0,y+v0,x)−2zw0,xy}=ϵ0+zκ
Штами є лише градієнтами переміщень; використовуючи матричні позначення їх можна записати, деϵ0 знаходиться деформація середньої площини іκ є вектором других похідних зміщення, що називається кривизною:
κ={κxκyκxy}={−w0,xx−w0,yy−2w0,xy}
Компонентκxy являє собою скручувальну кривизну,x вказуючи, як змінюється нахил середньої площиниy -directiony (або еквівалентно, як змінюється нахил напрямкуx).
Напруження щодоx−y осей тепер визначаються з деформацій, і при цьому потрібно враховувати, що кожен шар взагалі матиме різну жорсткість, залежно від власних властивостей, а також своєї орієнтації щодоx−y осей. Це враховується шляхом обчислення перетвореної матриці жорсткостіˉD, як описано в попередньому розділі (Equation 4.4.11). Нагадаємо, що жорсткість шару, задана рівнянням 4.4.4, - це ті, що вздовж волокна та поперечні напрямки цього конкретного шару. Властивості кожного шару необхідно трансформувати в загальніx−y осі, підібрані довільно для всього ламінату. Напруги при будь-якому вертикальному положенні тоді:
σ=ˉDϵ=ˉDϵ0+zˉDκ
де тутˉD трансформована жорсткість шару в положенні, при якому обчислюються напруги.
Кожен з цих шарів напружень повинен додати, щоб збалансувати тягу на одиницю шириниN:
N=∫+h/2−h/2σdz=N∑k=1∫zk+1zkσkdz
деσk напруга вk шарі іzk відстань від середини ламінату доk нижньої частини шару. Використання рівняння 4.4.16 для запису напружень через деформації та кривизни середньої площини:
N=N∑k=1(∫zk+1zkˉDϵ0dz+∫zk+1zkˉDκzdz)
Кривизнаκ і деформація середньої площиниϵ0 є постійними протягом усьогоz, і перетворена жорсткістьˉD не змінюється в межах заданого шару. Видалення цих величин зсередини інтегралів:
N=N∑k=1(ˉDϵ0∫zk+1zkdz+ˉDκ∫zk+1zkzdz)
Після оцінки інтегралів цей вираз можна записати в компактному вигляді:
N=Aϵ0+Bκ
деA - «матриця розтягнутої жорсткості», яка визначається як:
A=N∑k=1ˉD(zk+1−zk)
іB являє собою «матрицю жорсткості зчеплення», яка визначається як:
B=12N∑k=1ˉD(z2k+1−z2k)
Обґрунтування назв «extensional» і «couption» запропоновано рівнянням 4.4.20. AМатриця дає вплив розтягнутої середньої площини деформаціїϵ0 на внутрішню тягуN, аB матриця дає внесок кривизниκ в тягу. Може бути не очевидно, чому згинання плити зажадає внутрішньоплощинной тяги, або навпаки, чому витягування пластини в її площині призведе до її згинання. Але візуалізуйте пластину, що містить шари все тієї ж жорсткості, за винятком деяких дуже низькомодульних шарів десь над її середньою площиною. Коли пластина витягується, більш сумісні шари над середньою площиною, як правило, розтягуються більше, ніж жорсткіші шари нижче середньої площини. Верхня половина ламінату тягнеться більше, ніж нижня половина, тому вона приймає увігнуто-вниз кривизну.
Аналогічно, результуючі моменти на одиницю ширини повинні врівноважуватися моментами, що сприяють внутрішнім напруженням:
M=∫+h/2−h/2σzdz=Bϵ0+Dκ
деD «матриця жорсткості на вигин» визначається як:
D=13N∑k=1ˉD(z3k+1−z3k)
Повний набір зв'язків між прикладними силами та моментами, а також отримані деформації та кривизни середньої площини можна узагальнити як єдине матричне рівняння:
{NM}=[ABBD]{ϵ0κ}
A/B/B/DМатриця в дужках - це матриця жорсткості ламінату, а її зворотним буде матриця відповідності ламінату.
Наявність ненульових елементів в матриці зчепленняB вказує на те, що застосування в площині тяги призведе до викривлення або викривлення пластини, або що прикладений згинальний момент також буде генерувати подовжувальну деформацію. Ці ефекти, як правило, небажані. Однак їх можна уникнути, зробивши ламінат симетричним щодо середньої площини, як може виявити експертиза Рівняння 4.4.22. (У деяких випадках ця муфта подовжувач-кривизна може бути використана як цікава конструктивна особливість. Наприклад, можна сконструювати композитну лопатку гвинта, кут атаки якого автоматично змінюється зі швидкістю обертання: збільшена швидкість збільшує доцентрове навантаження в площині, що викликає скручування в лопаті.)
Наведені вище співвідношення забезпечують прямий (хоча і виснажливий, якщо не використовується комп'ютер) засіб визначення напружень і переміщень в ламінованих композиціях, що піддаються внутрішньоплощинним тяговим або згинальним навантаженням:
- Для кожного типу матеріалу в послідовності укладання отримують шляхом вимірювання або мікромеханічної оцінки чотири незалежних анізотропних параметра, що відображаються в Рівнянні 4.4.4: (E1,E2,ν12, іG12).
- Використовуючи Рівняння 4.4.11, перетворіть матрицю відповідності для кожного шару з основних напрямків матеріалу шару на деякі зручні опорні осі, які будуть використовуватися для ламінату в цілому.
- Інвертувати перетворену матрицю відповідності для отримання перетвореної (щодоx−y осей) матриці жорсткості. ˉD.
- Додайте внесок кожного шару доD матрицьA,B і, як це передбачено Eqns. 4.4.21, 4.4.22 та 4.4.24.
- Введіть задані тягиN і згинальні моментиM, і сформуйте системні рівняння, задані рівнянням 4.4.25.
- Розв'яжіть отриману систему для невідомих значень внутрішньоплоского деформаціїϵ0 та кривизниκ.
- Використовуйте Рівняння 4.4.16 для визначення напружень шару для кожного шару в ламінаті в перерахунку наϵ0,κ іz. Це будуть напруги щодоx−y осей.
- Використовуйте рівняння 4.4.6 для перетворенняx−y напружень назад до основних осей матеріалу (паралельних і поперечних волокнам).
- При бажанні окремі напруження шару можуть бути використані у відповідному критерії відмови для оцінки ймовірності виходу цього шару з ладу. Критерій Цай-Хілла широко використовується для цієї мети:
(σ1ˆσ1)2−σ1σ2ˆσ21+(σ2ˆσ2)2+(τ12ˆτ12)2=1
ˆσ2 Осьˆσ1 і міцність на розрив шарів паралельно і вздовж напрямку волокна,ˆτ12 а також міцність внутрішньомісячного шару. Цей критерій прогнозує невдачу, коли ліва сторона вищевказаного рівняння дорівнює або перевищує одиницю.
Приклад4.4.2
Аналіз ламінату, описаний вище, був реалізований в коді з назвою пластини, і цей приклад демонструє використання цього коду при визначенні жорсткості двошарового композиту 0/90 графіту/епоксидного композиту. Тут кожному з двох шарів дається товщина 0,5, тому загальна висота ламінату складе одиницю. Теорія ламінату передбачає одиницю ширини, тому загальна жорсткість і відповідність матриць будуть засновані на одиничному перерізі.
> plate assign properties for lamina type 1... enter modulus in fiber direction... (enter -1 to stop): 230e9 enter modulus in transverse direction: 6.6e9 enter principal Poisson ratio: .25 enter shear modulus: 4.8e9 enter ply thickness: .5 assign properties for lamina type 2... enter modulus in fiber direction... (enter -1 to stop): -1 define layup sequence, starting at bottom... (use negative material set number to stop) enter material set number for ply number 1: 1 enter ply angle: 0 enter material set number for ply number 2: 1 enter ply angle: 90 enter material set number for ply number 3: -1 laminate stiffness matrix: 0.1185E+12 0.1653E+10 0.2942E+04 -0.2798D+11 0.0000D+00 0.7354D+03 0.1653E+10 0.1185E+12 0.1389E+06 0.0000D+00 0.2798D+11 0.3473D+05 0.2942E+04 0.1389E+06 0.4800E+10 0.7354D+03 0.3473D+05 0.0000D+00 -0.2798E+11 0.0000E+00 0.7354E+03 0.9876D+10 0.1377D+09 0.2451D+03 0.0000E+00 0.2798E+11 0.3473E+05 0.1377D+09 0.9876D+10 0.1158D+05 0.7354E+03 0.3473E+05 0.0000E+00 0.2451D+03 0.1158D+05 0.4000D+09 laminate compliance matrix: 0.2548E-10 -0.3554E-12 -0.1639E-16 0.7218D-10 0.7125D-19 -0.6022D-16 -0.3554E-12 0.2548E-10 -0.2150E-15 0.3253D-18 -0.7218D-10 -0.1228D-15 -0.1639E-16 -0.2150E-15 0.2083E-09 -0.6022D-16 -0.1228D-15 0.2228D-19 -0.6022E-16 -0.1228E-15 0.2228E-19 -0.1967D-15 -0.2580D-14 0.2500D-08
Зверніть увагу, що цей несиметричний ламінат генерує ненульові значення в матриці зв'язкуB, як і очікувалося. Жорсткість однакова вx іy напрямках, в чому можна переконатися, розглянувши 1,1 і 2,2 елементи матриці відповідності ламінату. Ефективним модулем єEx=Ey=1/0.2548×10−10=39.2 ГПа. Однак ламінат не є ізотропним, як це можна виявити шляхом перекатки плити з укладанням 0/90, орієнтованим під іншим кутом відx−y осей.
Температурні ефекти
Існує ряд поліпшень, які можна розглянути для описаного вище коду пластини: він може бути розширений, щоб включити міжламінарні напруги зсуву між шарами, він може включати базу даних комерційно доступних препрегів і основних матеріалів, або користувальницький інтерфейс може бути зроблений\ більш дружнім» і графічно орієнтований. Багато таких функцій доступні в комерційних кодах, або можуть бути додані користувачем, і далі тут не будуть обговорюватися. Однак ефекти теплового розширення настільки важливі в застосуванні, що код ламінату майже повинен мати цю особливість, щоб бути корисним, і загальний підхід буде викладено тут.
Загалом, підвищення температуриΔT викликає теплове розширення, яке дається відомим співвідношеннямϵT=αΔT, деϵT знаходиться термічно індукована деформація іα коефіцієнт лінійного теплового розширення. Ця теплова деформація отримується без необхідності застосування напруги, так що коли закон Гука використовується для обчислення напруги від деформації, спочатку віднімається теплова складова:σ=E(ϵ−αΔT). Теплове розширення викликає лише нормальну деформацію, тому зсувні компоненти деформації не впливають. Таким чином, рівняння 4.4.3 може бути розширено як
де вектор теплової деформації в координатній рамці 1 - 2 дорівнює
ϵT={α1α20}ΔT
Осьα1 іα2 наведені анізотропні коефіцієнти теплового розширення в волоконному і поперечному напрямках. Трансформуючись в загальніx−y осі, це відношення стає:
{σxσyτxy}=[ˆD11ˆD12ˆD13ˆD21ˆD22ˆD23ˆD31ˆD32ˆD33]{[ϵxϵyϵxy]−[αxαyαxy]ΔT}
Індекси наˉD елементах стосуються позицій рядків та стовпців у матриці жорсткості, а не координатних напрямків; надпланка служить нагадуванням про те, що ці елементи стосуютьсяx−y осей. Вектор теплового розширення на правій стороні (α=αx,αy,αxy) по суті є вектором деформації, і тому його можна отримати з (α1,α2,0), як у Рівнянні 4.4.10:
α={σxσyσxy}=RA−1R−1{α1α20}
Зверніть увагу, що в загальномуx−y напрямку теплове розширення викликає як нормальні, так і зсувні деформації.
Попереднє незалежне від температури розвиток тепер можна повторити, модифікувавши лише шляхом проведення термінів теплового розширення. Як і раніше, вектор деформації для будь-якого положенняz від середньої площини задається через деформаціюϵ0 та кривизнуκ середньої площини
ϵ=ϵ0+zκ
Відповідний стрес тоді
Врівноваження напружень від застосованих тяг і моментів, як і раніше:
Цей результат ідентичний результату рівняння 4.4.20 та рівняння 4.4.23, крім додавання інтегралів, що представляють «теплові навантаження». Це дозволяє вирішувати проблеми, що залежать від температури, за допомогою «еквівалентної механічної формули»; загальні керуючі рівняння можуть бути записані як
{ˉNˉM}=[ABBD]{ϵ0κ}, or {ϵ0κ}=[ABBD]−1{ˉNˉM}
де\ еквівалентні теплові навантаження» задаються як
ˉN=N+∫ˉDαΔTdz
ˉM=M+∫ˉDαΔTzdz
Розширення коду пластини для розміщення теплових ефектів, таким чином, складається з модифікації вектора навантаження 6× 1 шляхом додавання двох векторних інтегралів 3× 1 у вищезгаданому виразі.
В'язкопружні ефекти
Оскільки матриця багатьох композитних ламінатів полімерна, конструктору може знадобитися розглянути можливість в'язкопружного розслаблення напружень або повзучості під час навантаження. Будь-який такий ефект, ймовірно, не буде великим, так як волокна, які несуть більшу частину навантаження, зазвичай не є в'язкоеластичними. Далі матеріал матриці зазвичай використовується значно нижче його температури склопереходу, і буде діяти в склоподібному еластичному режимі.
Однак деякі програми можуть бути не такими простими. Якщо ламінат використовується при підвищеній температурі, і якщо напруги діють в напрямках, не підтримуваних армуючими волокнами, можуть спостерігатися релаксаційні ефекти. На малюнку 6 показана повзучість, виміряна в однонаправленому графіто-епоксидному ламінаті T300/5208 (M.E. Tuttle and H.F. Brinson, «Прогнозування довгострокової відповідності повзучості загальних композитних ламінатів,» Експериментальна механіка, стор. 89, березень 1986.), завантажений поперечно до волокон на 149∘C. Навіть у цьому майже найгіршому випадку деформації повзучості відносно невеликі (менше 10% еластичної деформації), але малюнок 6 показує, що ефекти релаксації можуть бути важливими в деяких ситуаціях.

Наведена вище стаття Таттла-Брінсона описує покрокову обчислювальну схему, яка може бути використана для моделювання цих в'язкопружних ефектів ламінату, і тут буде викладена спрощена форма їх методу. В'язкопружна деформація повзучості, що виникає в даному шарі протягом певного приросту часу,dt може бути розрахована на основі напруги в шарі в той час, припускаючи, що шар вільний від сусідніх шарів; це дає незалежний шар деформації повзучості. Цей штам буде діяти, щоб розслабити шаровий стрес.
Звичайно, шари не можуть вільно напружуватися довільно, і належну сумісність деформацій можна відновити шляхом розрахунку зовнішніх навантажень, які б створювали пружні деформації, рівні незалежним шарам повзучості. Ці навантаження підсумовуються по всіх шарах в ламінаті, щоб дати еквівалентну навантаження на повзучість ламінату. Це навантаження застосовується до ламінату для обчислення набору сумісних деформацій і кривизн, які називаються еквівалентним деформацією повзучості ламінату. Ця деформація додається до початкової пружної деформації при обчисленні напруги на заданому шарі, тоді як незалежний шар деформації повзучості віднімається.
Наступний список більш детально розробляє ці кроки:
1. Пружні деформації та кривизни середньої площини вирішуються для заданих згинальних моментів і тяг, використовуючи склоподібні модулі різних шарів. З рівняння 4.4.25:
{ϵ0κ}=[ABBD]−1{NM}
2. Потім пружна деформація в кожному шарі отримується з Рівняння 4.4.15. Для шаруkth, з центром у координатіz, це:
ϵp_e_xy=ϵ0+zκ
деp_e_xy індекси позначають шар, пружний, напруга вx−y напрямку. Штами пружного шару відносно 1-2 (волоконно-поперечних) напрямків задаються перетворенням Рівняння 4.4.10:
ϵp_e_12=RAR−1ϵp_e_xy
Ці перші два кроки виконуються кодом пружної пластини, а адаптація до в'язкопружної реакції складається з додавання наступних кроків.
3. Поточна напруга шаруσk_12 в 1-2 напрямках становить:
Величинаϵp_lc_12−ϵp_c_12 - це різниця між еквівалентною деформацією повзучості ламінату та незалежним шаром деформації повзучості. Величиниϵp_lc_12 іϵp_c_12 встановлюються на нуль спочатку, але оновлюються в кроках 4 і 8 нижче, щоб врахувати в'язкопружну релаксацію.
4. Поточне напруження шару в потім використовується у відповідній в'язкопружній моделі для обчислення повзучості, яка виникла б, якби шар вільно деформувався незалежно від сусідніх шарів; це називається незалежною деформацією повзучості шару. Для простої моделі Voigt поточне значення деформації повзучості може бути оновлено з його значення на попередньому часовому кроці як:
де надскрипти на деформації вказують значення на поточному та попередньому кроках часу. CvОсь в'язкопружна відповідність повзучості іτ час релаксації. Штам повзучості, рівний,Cv буде розвиватися на додаток до початкової еластичної деформації в ламінаті, і частка цієї2/e деформації повзучості буде розвиватися через часτ. Різні значенняCv будуть використані для волокон, поперечних і зсувних деформаційних компонентів через анізотропії шару.
5. Напруження в напрямах 1-2 та\ (x-y), які були б необхідні для розвитку незалежних деформацій повзучості, якщо пласт був пружним, складають
σk_12=Dϵp_c_12
σk_xy=A−1σk_12
6. Ці еквівалентні напруження пружного шару підсумовуються над усіма шарами в ламінаті, щоб створити еквівалентну навантаження на повзучість ламінату. Внесокkth п'єса становить:
Nc=Nc+tkσk_xy
Mc=Mc+tkzσk_xy
деkth - товщина шару іkth координата йогоz осьової лінії.
7. Еквівалентна деформація повзучості ламінату потім обчислюється з еластичної матриці відповідності та еквівалентних навантажень на повзучість ламінату як
{ϵ0lcκlc}=[ABBD]−1{NcMc}
8. Шарова деформація повзучості ламінату вx−y і 1-2 напрямках
ϵp_pc_xy=ϵ0lc+zκlc
ϵp_lc_12=RAR−1ϵp_lc_xy
9. Нарешті, час збільшується,(t←t+dt) і інший часовий цикл обчислюється, починаючи з кроку 3.
Приклад4.4.3
В якості ілюстрації наведеного алгоритму розглянемо просту модель ламінату з одним ізотропним шаром. Пружні константи єE=100 (довільні одиниці) іν=0.25, а одиничне напруження прикладається вx -напрямку. Отже,x початкова деформація напрямкуϵx,0=σx/E=0.01. У цьому ізотропному тесті код обчислює модуль зсуву якG=E/2(1+ν). Деформація повзучості регулюється параметромvfrac, який встановлює відповідністьGv повзучості Voigtvfrac/E2 у поперечному напрямкуvfrac/G12 для компонентів зсуву та нуль у напрямку волокна (припускаючи лише еластичну реакцію вздовж волокон).
На малюнку 7 показана історія деформації повзучості цього ламінату для часу релаксаціїτ=1000 s Код кроків лінійно в журналі часу, в даному випадку з чотирма часовими кроками на десятиліття. Деформація повзучості - це деформація над початковою пружною деформацією, яка переходить від нуля до того,Cvϵx,0=5×10−4 як час прогресує через час релаксації.

Посилання
- Ештон, JE, J.C. Halpin і П.Х. Петі, грунтовка з композитних матеріалів: аналіз, Technomic Press, Вестпорт, Коннектикут, 1969.
- Джонс, Р.М., Механіка композиційних матеріалів, Макгроу-Хілл, Нью-Йорк, 1975.
- Пауелл, P.C., Інжиніринг з полімерами, Чепмен і Холл, Лондон, 1983.
- Ройланс, Д., Механіка матеріалів, Wiley & Sons, Нью-Йорк, 1996.
Вправа4.4.1
Випишітьx−y двовимірну матрицю відповідностіˉS та матрицю жорсткостіˉD (Рівняння 4.4.11) для одного шару графіту/епоксидного композиту з вирівняними поx−y осях волокнами.
Вправа4.4.2
Випишітьx−y двовимірну матрицю відповідностіˉS та матрицю жорсткостіˉD (Рівняння 4.4.11) для одного шару графіту/епоксидного композиту з вирівняними волокнами 30∘ відx осі.
Вправа4.4.3
Побудуйте ефективний модуль Юнга, виміряний вздовжx− осі, одного односпрямованого шару графіто-епоксидного композиту як функцію кута між напрямком волокна шару таx− віссю.
Вправа4.4.4
Використовуючи мову програмування за вашим вибором, напишіть код ламінату, подібний до коду пластини, згаданого в тексті, і перевірте його, обчисливши жорсткість ламінату і відповідність матриць, наведених у прикладі4.4.2.
Вправа4.4.5
(60∘/0∘/−60∘)SРозкладка (Sверхній індекс вказує, що шари повторюються для отримання симетричного ламінату) є прикладом того, що називаються «квазі-ізотропними» ламінатами, що мають рівні жорсткості вy напрямкахx і, незалежно від орієнтації ламінату. Переконайтеся, що це так для двох орієнтацій ламінату, один має 0∘ шарів, орієнтованих вздовжx осі, а інший з 0∘ шарів, орієнтованих на 30∘ відx осі.
Механічні властивості композиційних матеріалів
У наступній таблиці наведено значення фізичних та механічних властивостей для репрезентативних матеріалів шарів та сердечників, широко використовуваних у композитних ламінатах, армованих волокнами. Властивості шару взяті з F.P.Gerstle, «Композити», Енциклопедії полімерних наук та інженерії, Wiley, Нью-Йорк, 1991, з якими слід звернутися до даних з більш широкого спектру матеріалів. Див. також Г.Любін, Довідник з композитів, Ван Ностранд, Нью-Йорк, 1982.
S-скло/епоксидна смола | Кевлар/епоксидна смола | HM графіт/епоксидна смола | Сосна | Жорсткий пінопласт Rohacell 51 | |
Еластичні властивості: E1, ГПа |
|
|
|
|
|
Міцність на розрив: σ1, МПа |
1800 |
2000 |
1100 |
78 |
1,9 |
Міцність на стиск: σ1, МПа |
690 |
280 |
620 |
33 |
0,9 |
Фізичні властивості: α1,10−6/∘C |
2,1 |
-4,0 |
-0,7 |
|
33 |