فهرست مطالب

مکانیک هوافضا - سال سیزدهم شماره 1 (پیاپی 47، بهار 1396)

فصلنامه مکانیک هوافضا
سال سیزدهم شماره 1 (پیاپی 47، بهار 1396)

  • رفتار مکانیک مواد و سازه ها
  • تاریخ انتشار: 1396/02/12
  • تعداد عناوین: 9
|
  • مرتضی دزیانی، حسین دلایلی، شهرام یوسفی*، حمید نخ فروش صفحه 1
    مراحل مختلف در فرایند طراحی سازه های هوایی، اثرات متقابل موضوعات مختلف طراحی و لزوم در اختیار داشتن یک پیش بینی اولیه در مراحل مختلف طراحی، نیاز به وجود یک روش طراحی پارامتریک کلی برای اجزاء مختلف سازه ای را ضروری می سازد. روش طراحی بر اساس مفهوم شناسه سازه یک روش طراحی و تحلیل همزمان بوده که با دریافت پارامتر کلیدی شناسه سازه و خصوصیات مواد می تواند طراحی کاملی از مقطع پنل تقویت شده فشاری را در اختیار طراح قرار دهد. در این مقاله ضمن بررسی و تحلیل روش طراحی پنل تقویت شده با مفهوم شناسه سازه با پیکربندی کلی، نتایج آن از دو راهبرد مختلف استخراج شده و با انجام اصلاحاتی، نتایج روش شناسه سازه به نتایج تحلیل های چرخه ای مرسوم نزدیکتر شده است. در نهایت نتایج حاصل با نتایج روش های جاری طراحی مقایسه شده اند. روش شناسه سازه علاوه بر اینکه خود می تواند یک روش طراحی محسوب شود، به عنوان ابزاری برای تعیین حدود ناحیه بهینه و انتخاب نقاط اولیه طراحی در روش های بهینه سازی عددی بر پایه روابط تحلیلی و یا روش های اجزاء محدود قابل استفاده می باشد.
    کلیدواژگان: پنل تقویت شده فشاری، شناسه سازه، طراحی بهینه
  • محمد اسماعیل گلمکانی*، مصطفی امامی صفحه 13
    هدف از این تحقیق تحلیل غیر خطی خمش و کمانش قطاع صفحات مدرج تابعی در راستای شعاعی تحت بار مکانیکی می باشد. خواص مکانیکی صفحه تابعی در جهت شعاع با استفاده از مدل موری- تاناکا تعیین شده است. معادلات غیر خطی بصورت نموی و با استفاده از تئوری تغییر شکل های برشی مرتبه اول و فرضیه تغییر شکل های بزرگ فون کارمن ارائه گردیده اند. برای حل معادلات غیرخطی نموی از روش رهایی پویا و روش تفاوت های محدود استفاده شده است. با توجه به عدم وجود تحقیقی مشابه برای کمانش قطاع صفحات گرد تابعی با تغییرات خواص در راستای شعاعی تحت بارهای مکانیکی فشاری، برخی از نتایج با پاسخ های بدست آمده از نرم افزار المان محدود آباکوس برای صفحات تابعی مقایسه شده است و مطابقت خوب بدست آمده حاکی از صحت و دقت روش عددی بکار رفته می باشد. در انتها اثرات پارامترهای مختلفی همچون شرایط مرزی گوناگون، تاثیر توان ماده تابعی، نسبت ضخامت به شعاع خارجی، نسبت شعاع داخلی به شعاع خارجی، زاویه قطاع برای خمش و کمانش صفحات مورد بررسی قرار گرفته است.
    کلیدواژگان: خمش، کمانش، تئوری تغییر شکل برشی مرتبه اول، مواد تابعی
  • جواد شهبازی کرمی*، کرامت ملک زاده فرد، محمد احمدی صفحه 29
    در این مقاله توزیع تنش میان لایه ای در اتصالات چسبی تک لبه و جفت لبه مورد بررسی قرار گرفته است. یک معادله ی جابجایی درجه دوم برای چسب در راستای ضخامت در نظر گرفته شده است. برای تحلیل تنش این اتصال، نمودار جسم آزاد المانی از الیاف، در نواحی چسبی و غیر چسبی اتصال رسم شده و با نوشتن معادلات تعادل برای این المان ها، معادلات دیفرانسیل حاکم بر اتصال استخراج شده است. سپس، با حل معادلات حاکم و با استفاده از معادلات شرایط مرزی، مقادیر جابجایی در الیاف و در نتیجه تنش برشی ایجاد شده در چسب و بار محوری ایجاد شده در الیاف بدست آمده است. علاوه بر این موارد، اثرات انبساط گرمایی و نم گرمایی نیز در نظر گرفته شده است. نتایج نشان می دهند که توزیع تنش علاوه بر مشخصات هندسی هم چون ضخامت و طول ناحیه ی چسبی، به مشخصات فیزیکی اتصال از جمله مدول یانگ، ضریب پواسون و ضریب انبساط حرارتی لایه نیز بستگی دارد. برخی از انواع اتصالات توزیع تنش یکنواخت تری نسبت به بقیه دارند. مثلا تنش برشی بیشینه در اتصال جفت لبه تقریبا نصف اتصال تک لبه است. نتایج حل های تحلیلی با نتایج به دست آمده از حل اجزا محدود به وسیله ی نرم افزار ANSYS صحت سنجی شده است و تطابق خوبی بین این نتایج مشاهده شده است.
    کلیدواژگان: توزیع تنش، اتصال چسبی، مواد مرکب، ماتریس فلزی
  • فرهاد مرادی، سعید جعفری مهرآبادی* صفحه 39
    دراین مقاله، کمانش حرارتی پانل استوانه ای نازک از جنس ماده تابعی مدرج دوبعدی که در ساختار آن یک لایه ماده چندفازی مگنتو-الکترو-الاستیک باخواص الاستیکی خاص به منظوراثرپذیری سازه از میدان های مغناطیسی و الکتریکی استفاده شده، بررسی می شود. خواص ماده تابعی مدرج درجهات طول و محیط با کسرحجمی هایی از توابع نمایی درنظرگرفته شده اند. معادلات تعادل سازه مذکور براساس تئوری مرتبه اول تغییرشکل برشی و با درنظرگرفتن ترم های غیرخطی کرنش در روابط کرنش- تغییر مکان بدست آمده و سپس با ایجاد نمو در مولفه های تغییرمکان و اعمال معادلات اویلر بر تابع تغییرات دوم انرژی پتانسیل، معادلات پایداری بدست آمده اند. درادامه با نوشتن معادلات پایداری برحسب تغییر مکان ها و اعمال شرایط مرزی، با استفاده از روش مربعات دیفرانسیلی تعمیم یافته دو بعدی، دمای بحرانی کمانش سازه بدست آمده و با مقایسه نتایج حاصل از این روش با نتایج بدست آمده از روش المان محدود، اطمینان کافی از دقت روش حل حاصل شده است. درپایان اثرات خصوصیات هندسی سازه، ولتاژ مغناطیسی اعمالی به سطح خارجی پانل و تغییر خواص مواد بر اختلاف دمای بحرانی کمانش حرارتی پوسته مورد بررسی قرار گرفته است.
    کلیدواژگان: پیزو مغناطیس، ماده تابعی مدرج دوبعدی، کمانش حرارتی، روش مربعات دیفرانسیلی تعمیم یافته، پانل استوانه ای نازک و کوتاه
  • فاروق بابایی بیردایی*، بشیر بهجت صفحه 53
    نانو تیر یکی از مهمترین نانو ساختارها با کاربرد در زمینه های مختلف از جمله ساخت سیستم های نانو الکترومکانیکال با وزن کم، مصرف انرژی پایین و حساسیت بالا برای استفاده در پزشکی، رایانه، بیوحسگرها و... می باشد. از آنجا که رفتار مواد در مقیاس نانو تفاوت آشکاری با رفتار در مقیاس معمول دارد، روش های متنوع و جدیدی برای مطالعه رفتار مکانیکی مواد در مقیاس نانو ابداع شده اند؛ یکی از این روش ها استفاده از مدل غیرموضعی می باشد. در این مقاله رفتار تیر بر روی بستر الاستیک غیرخطی با در نظر گرفتن اثرات غیر موضعی در ابعاد نانو بررسی شده است. از تئوری تیر اویلر-برنولی برای مدلسازی تیر استفاده شده است. معادلات حاکم بر رفتار تیر با استفاده از روش جابجایی مجازی به دست آمده است و از روش تحلیلی برای حل معادلات حاکم بر مسئله استفاده شده است. مقادیر خیز بیشینه، بار کمانش بحرانی و فرکانس طبیعی برای شرط مرزی تکیه گاه ساده به ازای مقادیر مختلف ثابت فنری و پارامتر غیر موضعی به دست آمده است. نتایج به دست آمده نشان می دهد که لحاظ کردن اثرات غیر موضعی موجب افزایش خیز و کاهش نیروی کمانش بحرانی و فرکانس طبیعی می شود. همچنین بستر الاستیک موجب کاهش خیز و افزایش نیروی کمانش بحرانی و فرکانس طبیعی می شود.
    کلیدواژگان: تئوری غیرموضعی، بستر الاستیک پسترناک، تحلیل استاتیکی، کمانشی و ارتعاشی
  • بهزاد مشیری اول، محمد جعفری* صفحه 61
    در این مقاله، برای بررسی تمرکز تنش صفحات چندلایه ی کامپوزیتی متقارن دارای گشودگی مثلثی از روش متغیر مختلط لخنیتسکی استفاده شده است. روش ارائه شده توسط لخنیتسکی به حل صفحات ناهمسانگرد حاوی گشودگی دایره ای و بیضی شکل محدود می شود. برای استفاده از این روش در تحلیل صفحات دارای گشودگی مثلثی، با استفاده از تابع نگاشت همنوا، ناحیه خارج گشودگی مثلثی به محدوه خارج دایره ای به شعاع واحد نگاشته می شود. تاثیر زاویه چرخش گشودگی، چیدمان لایه ها، کشیدگی گشودگی، شعاع انحنای گوشه گشودگی و زاویه بار بر توزیع تنش اطراف گشودگی مثلثی بررسی می شود. در مواردی، به منظور بررسی درستی روش حل تحلیلی، نتایج حاصل از روش تحلیلی با نتایج اجزای محدود مقایسه می شود. نتایج ارائه شده در این تحقیق، نشان می دهد که روش حل حاضر برای تعیین تمرکز تنش در چندلایه های متقارن دارای گشودگی هایی با اشکال مختلف قابل ارائه است.
    کلیدواژگان: روش متغیر مختلط، حل تحلیلی، چندلایه متقارن، گشودگی مثلثی
  • جمشید فضیلتی*، علی قمریان، نادر آریایی فر صفحه 73
    در مقاله حاضر فرایند بهینه سازی یک هندسه جدارنازک فروریزشی جاذب انرژی بر مبنای روش طراحی آزمایش ها ارائه شده است. بهینه سازی برای بیشترین قابلیت جذب انرژی و کمترین شوک برخورد انجام شده است. بر این اساس تعداد چهار پارامتر کنترل هندسه انتخاب شده است که برای هر یک دو سطح مقداری پیش فرض مورد استفاده قرار داده شده است. در مجموع 16 مدل عددی با هندسه متفاوت تهیه و رفتار فروریزشی این هندسه ها به صورت مدل پوسته سه بعدی در نرم افزار LS-Dyna شبیه سازی شده است. مدل ماده مورد استفاده برای جاذب آلیاژ آلومینیمی با رفتار الاستیک-پلاستیک است که مشخصات آن از آزمون کشش استخراج شده است. مدل عددی جاذب در مجاورت سطوح تکیه گاهی و نگه دارنده و در برابر برخورد دینامیکی جرم متحرک با انرژی مشخص در راستای محوری تحلیل شده است. مقایسه نتایج شبیه سازی و بررسی بر اساس روش طراحی آزمایش ها به تعیین شکل و میزان اثرگذاری پارامترهای تغییرداده شده و ترکیب اثر آن ها بر رفتار نمونه جاذب منتهی شده است. نشان داده شد که اهمیت ضخامت سازه جاذب در رفتار آن قابل توجه است. همچنین مشاهده شد که علاوه بر این که شعاع انحنا در انتهای بسته هندسه عاملی مهم در تعیین مشخصات رفتار سازه جاذب است، بلکه بیشترین تاثیر را دارا است. علاوه بر این، برای تایید مدل عددی استفاده شده یک مساله شبه استاتیکی مشابه مورد تحلیل قرار داده شده است و با نتایج حاصل از آزمون آزمایشگاهی مقایسه شده است که همخوانی مناسب نتایج کیفیت مناسب شبیه سازی را تایید می کند.
    کلیدواژگان: بهینه سازی، بارگذاری محوری دینامیکی، سازه جدار نازک، انرژی جذب شده، طراحی آزمایش ها
  • هاشم بابایی*، توحید میرزابابای مستوفی، ابوالفضل درویزه، مجید علی طاولی صفحه 87
    استفاده روز افزون از امواج ضربه ای در صنعت، سبب رواج آزمایش های تجربی، مدل سازی های تحلیلی و عددی مختلف برای بررسی تاثیر امواج ضربه ای بر تغییر شکل فلزات شده است. در این مقاله به بررسی تجربی و عددی تغییر شکل پلاستیک ورق های مستطیلی گیردار فولادی و آلومینیومی تحت بارگذاری هیدرودینامیکی توسط سامانه چکش ضربه ای برای بیان خواص فلزات تحت اثر امواج ضربه ای پرداخته شده است. شبیه سازی المان محدود توسط نرم افزار تحلیلی ABAQUS/EXPLICIT انجام گرفته که به منظور اعتبار سنجی نتایج بدست آمده از آن، آزمایش هایی بر روی ورق های مستطیلی با ضخامت و سطوح انرژی متفاوت انجام شده است. این مقایسه نشان می دهد که نتایج عددی با نتایج تجربی سازگاری قابل قبولی دارد، بنابراین استفاده از مدل عددی برای پیش بینی خیز مرکز ورق مستطیلی تحت بارگذاری هیدرودینامیکی مطلوب است.
    کلیدواژگان: تغییر شکل پلاستیک، بارگذاری هیدرودینامیکی، چکش ضربه ای، شبیه سازی عددی، ورق مستطیلی کاملا گیردار
  • بیژن محمدی*، مرتضی حلمی، داود سلیمی مجد صفحه 95
    هدف از این تحقیق، تخمین عمر خستگی کامپوزیت های چندلایه با استفاده از یک مدل جدید مبتنی بر مکانیک خرابی محیط پیوسته است. بدین منظور یک معیار به فرم بسته بر پایه روش انرژی و در چارچوب مکانیک خرابی ارائه شده است، به گونه ای که تفاوت آسیب در سه جهت مختلف را لحاظ می نماید و در عین حال محدود به لایه چینی خاصی نمی باشد. مدل پیشنهادی با استفاده از ماده نوشته شده توسط کاربر، در نرم افزار اجزاء محدود ANSYS (یوزرمت) پیاده سازی شده است. در این مقاله نحوه مشخصه سازی ثوابت مدل بیان شده و برای کامپوزیت کربن- اپکسی AS4/3501-6، ثوابت محاسبه شده اند و در نهایت اعتبارسنجی مدل برای کامپوزیت های لایه ای تک جهته و چندجهته مورد ارزیابی قرار گرفته است. در چندلایه های چندجهته، نتایج عمر خستگی برای چند لایه های متعامد دارای سوراخ [04/904]s و [904/04]s با نتایج آزمایشگاهی مقایسه شده است. نتایج نشان می دهد که مدل ارائه شده می تواند با تقریب بسیار قابل قبولی عمر چندلایه های کامپوزیتی را تنها با استفاده از ثوابت به دست آمده از آزمایش چندلایه های تک جهته در دو نسبت تنش محاسبه نماید.
    کلیدواژگان: عمر خستگی، کامپوزیت چند لایه، مکانیک خرابی محیط پیوسته، متغیر آسیب، چگالی انرژی کرنشی
|
  • M. Dezyani, H. Dalayeli, Sh. Yousefi*, H. Nakhforoosh Page 1
    Different steps of design process in aerospace structures, interaction of design parameters in the complete design process and the necessity of the preliminary prediction of structural geometrydictate the presence of a general parametric design method for preliminary sizing of different structural components. The design method of structural index concept is a simultaneous design and analysis process that can provide complete details of panel geometry by having the structural index and material properties of a stiffened panel. In this paper, the design and analysis process of compressed stiffened panel is described by structural index concept with general configuration and the results are extracted with two different approaches. Then the results are modified to confirm the common iterative methods. Finally, the results are compared with the results of current methods. In addition to its application as a design method, structural index concept also can be used to determine optimal region of design space and choose initial design points which required for numerical optimization methods based on analytical or F.E.M methods.
    Keywords: Structural Index Method, Structural Optimization, Stiffened Panel, Dimensional Sizing
  • M.E. Golmakani*, M. Emami Page 13
    The main purpose of this study is to investigate the nonlinear bending and buckling analysis of radially functionally graded sector plates subjected to uniform in-plane compressive loads. The mechanical properties of plates assumed to vary continuously along the radial direction by the Mori–Tanaka distribution. The incremental form of nonlinear formulations are derived based on first order shear deformation theory (FSDT) and large deflection von Karman equations. The dynamic relaxation (DR) method combined with the finite difference discretization technique is employed to solve the equilibrium equations. Also, due to the lack of similar research for the buckling of functionally graded sector plates with material variation in the radial direction, some results are compared with the ones reported by the ABAQUS finite element software. The achieved good agreements between the results indicate the accuracy of the present numerical method. Finally, numerical results for the maximum displacement and critical buckling load for various boundary conditions, effects of grading index, thickness-to-radius ratio, sector angle and inner radius-to-outer radius ratio are presented.
    Keywords: Bending, Buckling, First Order Shear Deformation Theory, Functionally Graded Material
  • J. Shahbazi Karami*, K. Malekzade Fard, M. Ahmadi Page 29
    In this article, the interfacial stresses in single and double lap joints are studied. A quadratic displacement field is assumed in the adhesive layer. For Analyzing stress in joint the free body diagram element of fiber in adhesive and unadhesive part of join has been drawn. By writing the static equation for this elements, the differential equation for joint has been derived. Then by solving this differential equation and using the boundary condition equation, displacement in fibers and so the shear stress in adhesive and axial force in fibers has been derived. The effect of thermal and hygrothermal stresses are also included. It is observed that the peak shear stress developed within the adhesive layer is function of mechanical and physical parameters as Young’s module, Poisson’s ratio, thickness, thermal and hygrothermal coefficients of expansion. The results show that the double lap joint appears to reduce the peak shear stress within the joint significantly. For example maximum shear stress in double lap joints are nearly the half of single lap joints. The results are validated by FEM solution. An excellent agreement is observed between analytical and FE methods.
    Keywords: Stress Distribution, Composite, Adhesive Joint, Metal Matrix
  • F. Moradi, S. Jafari* Page 39
    In this paper thermal buckling of a thin cylindrical panel made of piezo-magnetic two dimensional functionally graded materials (2D-PFGM) subjected to magnetic field have been investigated. The material properties of structure assumed varying by exponential functions of volume fraction in longitudinal and circumferential directions. In order to solving the problem at the first, the equilibrium equations have been derived by the first order shear deformation theory by concerning the nonlinear terms from the strain-displacement relations. In the next step, by giving the incremental changes on the displacement components and calculating the resultant forces and moments, using Lagrang,s equations on the second functional of potential energy the stability equations have been derived. After solution the mentioned equations by generalized differential quadrature method based on the simply supported boundary conditions the critical buckling temperature has been determined. By solving the numerical example in special cases and comparing the obtained results we confident from the method of analysis. At the end, effect of geometrical features, applied voltage on the external surface of panel, and changing the material properties on the critical buckling load have been investigated.
    Keywords: Piezo, Magnetic, 2D, FGM, Thermal Buckling, GDQ Method, Thin, Short Cylindrical Panel
  • F. Babaei Birdaei*, B. Behjat Page 53
    Nano-beam is one of the most important nano-structures with applications in different fields such as NEMS (Nano electromechanical systems) considering their low weight, low energy consumption and high sensitivity for use in medical, computer, bio-sensors and etc. Since the behavior of material in nano scale is different from usual scales, new methods are innovated to study the mechanical behavior of material in nano scale which one of them is nonlocal model. In this paper, nonlocal elasticity theory is used to describe the behavior of beam on the nonlinear elastic foundation. The Euler-Bernoulli beam theory is used to modeling the beam. The governing equations are obtained using principle of virtual displacements and analytical solution is used to solve governing equations. Maximum deflection, critical buckling load and natural frequencies for simply supported boundary condition and different spring constant, and nonlocal parameter is presented. The obtained results show that deflection is increased and critical buckling load and natural frequencies are decreased by considering nonlocal effects. Also elastic foundation decreases deflection and increases critical buckling load and natural frequencies.
    Keywords: Nonlocal Theory, Pasternak Elastic Foundation, Vibration, Buckling, Static Analysis
  • B. Moshiri Aval, M. Jafari* Page 61
    In this paper, the complex variable method of Lekhnitskii has been applied for solving stress concentration problems in symmetric laminates with triangular hole. Lekhnitskii solution is limited to circular and elliptical holes in anisotropic plates. In order to use this approach to triangular hole, by means of conformal mapping, the area external to the hole can be represented by the area outside the unit circle. The effect of rotation angle of hole, stacking sequence, hole curvature and load angle on stress distribution around triangular hole are considered. The results based on analytical solution are compared with the results obtained using finite element methods. Good agreement is observed and provides confidence in the accuracy of the present results. The results presented in this paper, Indicated that the presented method can be used to determine accurately the stress concentration in symmetric laminates with special shape holes.
    Keywords: Complex Variable Method, Analytical Solution, Symmetric Laminate, Triangular Hole
  • J. Fazilati*, A. Ghamarian, N. Ariaeifar Page 73
    In the present paper the optimization process of a thin-walled collapsible shock absorber geometry is performed using the well-known design of experiments method. The optimizing target parameters are the maximum energy absorption and the low impact shock. A number of four geometrical control parameters have been chosen in two different levels making a totally 16 different geometries. The collapse behavior of the absorbent geometries then has been simulated using a full 3D shell element LS-Dyna numerical model with an elastic-plastic material behavior and some main absorbing characteristics extracted. The collapse behavior of the geometries is simulated under axial impact of a traveling mass with a presumed kinetic energy. A comparison between the different geometries can be performed to obtain the effectiveness of the control parameters and decide about the most effective ones. It has been shown that while the thickness of the crushable structure is an important factor, an increase in the radius of curvature at the end of the presumed geometry is the most effective parameter in absorber efficiency. A study also is performed to validate the numerical simulation process with conical absorber geometry under axial quasi-static loading. The comparison exhibited a very good agreement between the numerical finite element results and the data acquired from the experimental test that showed the validity of the simulation procedure.
    Keywords: Optimization, Dynamic Axial Loading, Thin, Walled Structure, Absorbed Energy, Design of Experiment
  • H. Babaei*, T. Mirzababaie Mostofi, A. Darvizeh, M. Alitavoli Page 87
    Nowadays due to increasing industrial utilization of shock waves, it became popular to use analytical, numerical modeling and experimental tests to evaluate effect of shock waves on metal deformation. To explain this effect on metal properties, current paper considers both experimental and numerical plastic deformation of clamped rectangular steel and aluminum plates, under hydrodynamic loading of drop hammer system. To validate the explicit finite element simulation by solver ABAQUS, experimental test on rectangular plates with various thickness and energy levels have been performed. Obtained results by the numerical simulation are in a very good agreement with those in the experimental tests. Hence, it is desirable to utilize the numerical modeling to predict mid-point deflection of rectangular plate under hydrodynamic loading.
    Keywords: Plastic Deformation, Hydrodynamic Loading, Drop Hammer, Numerical Simulation, Fully Clamped Rectangular Plate
  • B. Mohammadi*, M. Helmi, D. Salimi, Majd Page 95
    The purpose of this study is to estimate the fatigue life of laminated composites using a novel model based on continuum damage mechanics. For this purpose, a closed form criterion based on the energy method has been presented in the context of damage mechanics such that considers the difference of damage in three directions and moreover is not limited to a special layup. The proposed model has been implemented in the ANSYS finite element software using the material coded by the user (Usermat). The method of characterization of the model constants has been presented in this paper and the constants have been determined for AS4/3501-6 Carbon-Epoxy composite and finally, validation of the model for unidirectional and multidirectional layered composites has been evaluated. For multidirectional laminates, results of fatigue life for cross-ply laminates with hole with the layups [04/904]s and [904/04]s compared with experimental results. The results show that the developed model can predict the fatigue life of laminated composites only by using the material constants which is obtained from experiments of unidirectional laminates in two stress levels.
    Keywords: Fatigue Life, Laminated Composite, Continuum Damage Mechanics, Damage Variable, Strain Energy Density