Главная » Просмотр файлов » Hutton - Fundamentals of Finite Element Analysis

Hutton - Fundamentals of Finite Element Analysis (523155), страница 45

Файл №523155 Hutton - Fundamentals of Finite Element Analysis (Hutton - Fundamentals of Finite Element Analysis) 45 страницаHutton - Fundamentals of Finite Element Analysis (523155) страница 452013-09-15СтудИзба
Просмтор этого файла доступен только зарегистрированным пользователям. Но у нас супер быстрая регистрация: достаточно только электронной почты!

Текст из файла (страница 45)

Ifthe lateral surfaces do not exhibit convection (i.e., the surfaces are insulated), theconvection terms are removed by setting h = 0 . Note that, in many finite elementsoftware packages, the convection portion of the conductance matrix is not automatically included in element matrix formulation. Instead, lateral surface (aswell as edge) convection effects are specified by applying convection “loads” tothe surfaces as appropriate. The software then modifies the element matrices asrequired.The element forcing functions are described in column matrix (vector)form as (e) TQ[N ] t d A =Q{N } t d AfQ =Aff(e) h(e) gA= 2h Ta=−[N ] d A = 2h Ta{N } d ATAqs n s [N ] T t dS = −S(7.38)Aqs n s {N } t dSSwhere [N ] = {N } is the M × 1 column matrix of interpolation functions.Equations 7.36–7.38 represent the general formulation of a finite element fortwo-dimensional heat conduction with convection from the surfaces.

Note inparticular that these equations are valid for an arbitrary element having M nodesand, therefore, any order of interpolation functions (linear, quadratic, cubic, etc.).In following examples, use of specific element geometries are illustrated.T7.4.2 Boundary ConditionsThe boundary conditions for two-dimensional conduction with convection maybe of three types, as illustrated by Figure 7.8 for a general two-dimensionaldomain. On portion S1 of the boundary, the temperature is prescribed as a knownHutton: Fundamentals ofFinite Element Analysis7. Applications in HeatTransferText© The McGraw−HillCompanies, 20047.4 Heat Transfer in Two Dimensionsh(T Ta)S3T T *S1q*S2Figure 7.8 Types of boundaryconditions for two-dimensionalconduction with convection.constant value TS1 = T ∗ .

In a finite element model of such a domain, every element node located on S1 has known temperature and the corresponding nodalequilibrium equations become “reaction” equations. The reaction “forces” arethe heat fluxes at the nodes on S1 . In using finite element software packages, suchconditions are input data; the user of the software (“FE programmer”) enterssuch data as appropriate at the applicable nodes of the finite element model (inthis case, specified temperatures).The heat flux on portion S2 of the boundary is prescribed as q S2 = q ∗ . Thisis analogous to specified nodal forces in a structural problem.

Hence, for all elements having nodes on S2 , the third of Equation 7.38 gives the correspondingnodal forcing functions as (e) f g = − q ∗ n S2 {N }t dS(7.39)S2Finally, a portion S3 of the boundary illustrates an edge convection condition. In this situation, the heat flux at the boundary must be equilibrated by theconvection loss from S3 .

For all elements having edges on S3 , the convectioncondition is expressed as (e) f g = − q S3 n S3 {N }t dS = − h(T (e) − Ta ){N }t dS(7.40)S3S3Noting that the right-hand side of Equation 7.40 involves the nodal temperatures, we rewrite the equation as (e) Tf g = − h[N ] [N ]{T }t dS3 + h Ta {N }t dS3(7.41)S3S3and observe that, when inserted into Equation 7.36, the first integral term on theright of Equation 7.41 adds stiffness to specific terms of the conductance matrix241Hutton: Fundamentals ofFinite Element Analysis2427. Applications in HeatTransferCHAPTER 7Text© The McGraw−HillCompanies, 2004Applications in Heat Transferassociated with nodes on S3 .

To generalize, we rewrite Equation 7.41 asf(e) g (e) (e) = − k h S {T } + f h S(7.42)wherek(e) hS=h[N ] T [N ]t dS(7.43)Sis the contribution to the element conductance matrix owing to convection onportion S of the element boundary and (e) f hS =h Ta {N }t dS(7.44)Sis the forcing function associated with convection on S.Incorporating Equation 7.42 into Equation 7.36, we have (e) (e) (e) (e) k (e) {T } = f Q + f h + f g + f h Swhere the element conductance matrix is now given byT T (e) ∂N∂N∂N∂Nk=+ kyt dAkx∂x∂x∂y∂yAT+ 2h[N ] [N ] d A + h [N ] T [N ] t dSA(7.45)(7.46)Swhich now explicitly includes edge convection on portion(s) S of the elementboundary subjected to convection.EXAMPLE 7.443sr1Determine the conductance matrix (excluding edge convection) for a four-node, rectangular element having 0.5 in.

thickness and equal sides of 1 in. The material has thermalproperties k x = k y = 20 Btu/(hr-ft-◦ F ) and h = 50 Btu/(hr-ft2-◦ F ).■ Solution2Figure 7.9 Elementnode numberingfor Example 7.4;the length of eachedge is 1 in.The element with node numbers is as shown in Figure 7.9 and the interpolation functions,Equation 6.56, areN 1 (r, s) =1(1 − r )(1 − s)4N 2 (r, s) =1(1 + r )(1 − s)4Hutton: Fundamentals ofFinite Element Analysis7. Applications in HeatTransferText© The McGraw−HillCompanies, 20047.4 Heat Transfer in Two DimensionsN 3 (r, s) =1(1 + r )(1 + s)4N 4 (r, s) =1(1 − r )(1 + s)4in terms of the normalized coordinates r and s.

For the 1-in. square element, we have2a = 2b = 1 and d A = dx dy = ab dr ds . The partial derivatives in terms of the nor-malized coordinates, via the chain rule, are∂ Ni∂ N i ∂r1 ∂ Ni==∂x∂r ∂ xa ∂ri = 1, 4∂ Ni∂ Ni ∂ s1 ∂ Ni==∂y∂s ∂yb ∂si = 1, 4Therefore, Equation 7.37 becomesk(e)T T 1 1 ∂N∂N 1∂N∂N 1kxt ab dr ds=+ ky∂r∂r a 2∂s∂ s b2−1 −11+ 2h[N ] T [N ]ab dr ds−1or, on a term by term basis,1 1 ki j =−1 −1∂ Ni ∂ N j 1∂ Ni ∂ N j 1kxt ab dr ds+ ky∂r ∂r a 2∂ s ∂ s b21 1+ 2hN i N j ab dr dsi, j = 1, 4−1 −1or1 1 ki j =kx−1 −1∂ Ni ∂ N j b∂ Ni ∂ N j at dr ds+ ky∂r ∂r a∂s ∂s b1 1+ 2hN i N j ab dr dsi, j = 1, 4−1 −1Assuming that k x and k y are constants, we haveki j = k x tba1 1−1 −11∂ Ni ∂ N jadr ds + k y t∂r ∂rb1 1−1 −11+ 2 habN i N j dr ds−1 −1i, j = 1, 4∂ Ni ∂ N jdr ds∂s ∂s243Hutton: Fundamentals ofFinite Element Analysis2447. Applications in HeatTransferCHAPTER 7Text© The McGraw−HillCompanies, 2004Applications in Heat TransferThe required partial derivatives are∂ N11= (s − 1)∂r4∂ N11= (r − 1)∂s4∂ N21= (1 − s)∂r4∂ N21= − (1 + r )∂s4∂ N31= (1 + s)∂r4∂ N31= (1 + r )∂s41∂ N4= − (1 + s)∂r4∂ N41= (1 − r )∂s4Substituting numerical values (noting that a = b ), we obtain, for example,1 1 k 11 = 20−1 −111(s − 1) 2 +(r − 1) 216161 11(1 − r ) 2 (1 − s) 216+ 2(50)−1 −1Integrating first on r,20(0.5)=16(12)k 11100−161−10.5120.512dr ds2dr ds11(r − 1) 3 (s − 1) 2 r −1 + ds3−10.5122 11(1 − r ) 3 (1 − s) ds3−12−1ork 1120(0.5)=16(12)1 8(s − 1) (2) +32−1100ds +160.5122 1(1 − s) 2−18ds3Then, integrating on s, we obtaink 1120(0.5)=16(12)ork 11 =20(0.5)16(12) 2 112(s − 1) 3(1 − s) 38100 0.58 + s −3316 1233−1−11616+33+100160.512 2 88= 0.6327 Btu/(hr-◦ F )33The analytical integration procedure just used to determine k 11 is not the method used byfinite element software packages; instead, numerical methods are used, primarily theGauss quadrature procedure discussed in Chapter 6.

If we examine the terms in the integrands of the equation defining k i j , we find that the integrands are quadratic functionsHutton: Fundamentals ofFinite Element Analysis7. Applications in HeatTransferText© The McGraw−HillCompanies, 20047.4 Heat Transfer in Two Dimensionsof r and s.

Therefore, the integrals can be evaluated exactly by using two Gauss pointsin r and s. Per Table 6.1, the required Gauss points and weighting factors are r i , s j =± 0.57735 and W i , W j = 1.0 , i, j = 1, 2 . Using the numerical procedure for k 11 , we writek 11b= kx ta1 1−1 −11+ 2hab−1= kx tb1(s − 1) 2 dr ds + k y t16a1 1−1 −11(r − 1) 2 dr ds161(r − 1) 2 (s − 1) 2 dr ds16222 2 b1a1W i W j (s j − 1) 2 + k y tW i W j (r i − 1) 2a i =1 j =1 16b i =1 j =1 16+ 2hab22 1W i W j (1 − r i ) 2 (1 − s j ) 216i =1 j =1and, using the specified integration points and weighting factors, this evaluates tok 11b= kx ta 14a 1+ kyt+ 2hab3b 39It is extremely important to note that the result expressed in the preceding equation is thecorrect value of k 11 for any rectangular element used for the two-dimensional heat conduction analysis discussed in this section.

The integrations need not be repeated for eachelement; only the geometric quantities and the conductance values need be substituted toobtain the value. Indeed, if we substitute the values for this example, we obtaink 11 = 0.6327 Btu/(hr-◦ F )as per the analytical integration procedure.Proceeding with the Gaussian integration procedure (calculation of some of theseterms are to be evaluated as end-of-chapter problems), we findk 11 = k 22 = k 33 = k 44 = 0.6327 Btu/(hr-◦ F )Why are these values equal?The off-diagonal terms (again using the numerical integration procedure) are calculated ask 12 = −0.1003k 13 = −0.2585k 14 = −0.1003k 23 = −0.1003k 24 = −0.2585k 34 = −0.1003245Hutton: Fundamentals ofFinite Element Analysis2467.

Applications in HeatTransferCHAPTER 7Text© The McGraw−HillCompanies, 2004Applications in Heat TransferBtu/(hr-◦ F ), and the complete element conductance matrix is0.6327 −0.1003 −0.2585 −0.1003 (e)  −0.1003 0.6327 −0.1003 −0.2585 ◦k= −0.2585 −0.1003 0.6327 −0.1003  Btu/(hr- F )−0.1003 −0.2585 −0.1003 0.6327EXAMPLE 7.5Figure 7.10a depicts a two-dimensional heating fin. The fin is attached to a pipe on itsleft edge, and the pipe conveys water at a constant temperature of 180 ◦ F .

Характеристики

Тип файла
PDF-файл
Размер
5,36 Mb
Тип материала
Учебное заведение
Неизвестно

Список файлов книги

Свежие статьи
Популярно сейчас
А знаете ли Вы, что из года в год задания практически не меняются? Математика, преподаваемая в учебных заведениях, никак не менялась минимум 30 лет. Найдите нужный учебный материал на СтудИзбе!
Ответы на популярные вопросы
Да! Наши авторы собирают и выкладывают те работы, которые сдаются в Вашем учебном заведении ежегодно и уже проверены преподавателями.
Да! У нас любой человек может выложить любую учебную работу и зарабатывать на её продажах! Но каждый учебный материал публикуется только после тщательной проверки администрацией.
Вернём деньги! А если быть более точными, то автору даётся немного времени на исправление, а если не исправит или выйдет время, то вернём деньги в полном объёме!
Да! На равне с готовыми студенческими работами у нас продаются услуги. Цены на услуги видны сразу, то есть Вам нужно только указать параметры и сразу можно оплачивать.
Отзывы студентов
Ставлю 10/10
Все нравится, очень удобный сайт, помогает в учебе. Кроме этого, можно заработать самому, выставляя готовые учебные материалы на продажу здесь. Рейтинги и отзывы на преподавателей очень помогают сориентироваться в начале нового семестра. Спасибо за такую функцию. Ставлю максимальную оценку.
Лучшая платформа для успешной сдачи сессии
Познакомился со СтудИзбой благодаря своему другу, очень нравится интерфейс, количество доступных файлов, цена, в общем, все прекрасно. Даже сам продаю какие-то свои работы.
Студизба ван лав ❤
Очень офигенный сайт для студентов. Много полезных учебных материалов. Пользуюсь студизбой с октября 2021 года. Серьёзных нареканий нет. Хотелось бы, что бы ввели подписочную модель и сделали материалы дешевле 300 рублей в рамках подписки бесплатными.
Отличный сайт
Лично меня всё устраивает - и покупка, и продажа; и цены, и возможность предпросмотра куска файла, и обилие бесплатных файлов (в подборках по авторам, читай, ВУЗам и факультетам). Есть определённые баги, но всё решаемо, да и администраторы реагируют в течение суток.
Маленький отзыв о большом помощнике!
Студизба спасает в те моменты, когда сроки горят, а работ накопилось достаточно. Довольно удобный сайт с простой навигацией и огромным количеством материалов.
Студ. Изба как крупнейший сборник работ для студентов
Тут дофига бывает всего полезного. Печально, что бывают предметы по которым даже одного бесплатного решения нет, но это скорее вопрос к студентам. В остальном всё здорово.
Спасательный островок
Если уже не успеваешь разобраться или застрял на каком-то задание поможет тебе быстро и недорого решить твою проблему.
Всё и так отлично
Всё очень удобно. Особенно круто, что есть система бонусов и можно выводить остатки денег. Очень много качественных бесплатных файлов.
Отзыв о системе "Студизба"
Отличная платформа для распространения работ, востребованных студентами. Хорошо налаженная и качественная работа сайта, огромная база заданий и аудитория.
Отличный помощник
Отличный сайт с кучей полезных файлов, позволяющий найти много методичек / учебников / отзывов о вузах и преподователях.
Отлично помогает студентам в любой момент для решения трудных и незамедлительных задач
Хотелось бы больше конкретной информации о преподавателях. А так в принципе хороший сайт, всегда им пользуюсь и ни разу не было желания прекратить. Хороший сайт для помощи студентам, удобный и приятный интерфейс. Из недостатков можно выделить только отсутствия небольшого количества файлов.
Спасибо за шикарный сайт
Великолепный сайт на котором студент за не большие деньги может найти помощь с дз, проектами курсовыми, лабораторными, а также узнать отзывы на преподавателей и бесплатно скачать пособия.
Популярные преподаватели
Добавляйте материалы
и зарабатывайте!
Продажи идут автоматически
6455
Авторов
на СтудИзбе
305
Средний доход
с одного платного файла
Обучение Подробнее