الخميس، 12 أكتوبر 2017

تحليل القواعد بإستعمال Linear spring على ال ETABS

      في هذا المقال سوف نتحدث عن الكيفية التي يمكن من خلالها نمذجة و تحليل القواعد الخرسانية بإستعمال برنامج ETABS .

القواعد هي ذلك الجزء من المبنى التي تنقل كافة الأحمال إلى التربة , و لها عدة أنواع في هذا المقال سوف نتعامل مع  القواعد منفردة single footings والقواعد مشتركة combined footing.
    تكمن الإشكالية في نمذجة القواعد بشكل رئيسي في نمذجة التفاعل بين القاعدة و التربة soil-base interaction , أحد النماذج التي يمكن إتباعها هو إستعمال النوابض الخطية Linear springs لتمثل التشوه Deformation الذي يحصل في التربة نتيجة الأحمال.
   تخضع النوابض الخطية إلى ما يسمى بـ Hooks law حيث ينص هذا القانون على أن القوة في النابض تتناسب خطيا مع التشوه الحاصل في النابض , و ثات التناسب في هذا القانون يسمى K : Stiffness coefficient و هو مقدار القوة اللازمة لإحداث إزاحة بمقدار 1m  , و تمثل هذه القيمة ما يسمى بالـ subgrade modulus في برامج مثل ال CSI-SAFE , و يمكن كتابة القانون على الصورة :
F = Ku 
حيث F : القوة بوحدة نيوتن N 
u : الإزاحة بالمتر m 
 K : Stiffness coefficient N/m 

   يمثل الرسم البياني الاتي التناسب الخطي بين القوة و الإزاحة . 

تتلخص عملية النمذجة بشكل رئيسي في اعتبار القاعدة بأنها Shell و تقسيمها إلى أجزاء صغيرة , ثم تعيين إسناد من نوع spring عند كل Node , و كل spring منها يمتلك قيمة K مختلفة حسب المساحة التي تحملها ال Node كما سوف نبين لاحقا. 
و يبن الشكل تلخيص لعملية النمذجة : 
لتوضيح ذلك , سوف نقوم بنمذجة قاعدة أبعادها 2x2 m و سماكة 0.4 m و تحمل في مركزها عامود بأبعاد 0.4x0.4 الذي يحمل بدوره حملا مركزا مقداره 400 kN , كما هو موضح في الشكل : 
1. بعد تشغيل برنامج ETABS 2015 , أبدأ مشروع جديد , ثم أختر Use Built in settings with >> Metric SI 
2. تظهر نافذة Quick Templates قم بإدخال الآتي : 
لتعيين حدود القاعدة :
Number of grid lines in x direction = 2  
Number of grid lines in Y direction = 2
لتعيين أبعاد القاعدة : 
Spacing of grids in X direction  = 2 m 
Spacing of grids in Y direction = 2 m
و لأننا سوف نحلل في بعدين فقط فإننا لن نحتاج لإضافة طوابق 
Number of stories = 1 

3. لتظهر مساحة العمل بالشكل الإتي : 

4. تعيين المقاطع و المواد : 
- المواد 
في هذا المثال سوف نستعمل المواد الإفتراضية الموجودة في البرنامج , من قائمة define >> Materials properties , سوف نستعمل مادة 4000 psi , يمكنك الدخول إليها لمعرفة خصائصها كما في الصورة .
- المقاطع 
يمكن تمثيل القاعدة كعنصر Shell و يمتلك thickness مقدارها 0.4 m , 
- من قائمة Define >> Section properties أختر Slab section 
- في نافذة slab properties اختر add new property 

- في نافذة Slab property data , قم بتغيير modeling type إلى shell-thick , و إدخال قيمة ال thickness بمقدار  400 mm , كما في الصورة . 

5. رسم القاعدة 
قم برسم القاعدة بإستعمال أداة draw rectangular floor , و في نافذة Properties of object قم بإختيار المقطع الذي قمنا بتعريفه سابقا


6. تقسيم القاعدة meshing و تعيين الإسناد restrains 
    في حالة الإستناد على التربة فإن الإستناد يكون مستمر على كامل مساحة القاعدة , في ال ETABS لا يمكن تعيين مسند لمساحة, (فقط يمكن تعينه عند نقاط Nodes) لذلك سوف نقوم بتقسيم القاعدة إلى مساحات صغيرة بحيث نحصل على عدد ملائم من نقاط الإسناد. 
     أيضا , في حالة التربة فإن أستعمال إسناد من نوع Pinned لن يمثل السلوك الحقيقي للتربة حيث أن هذا الإسناد يمنع الحركة بالأتجاه العامودي و في المقابل فإن التربة مادة قابلة للإنضغاط يمكن ان تتحرك نتيجة الاحمال نحو الأسفل.
   بدلا من إستعمال مسند تقليدي "إن جاز التعبير" سوف نستعمل مسند من نوع spring , حيث أن قابلية النابض للإنضغاط يمكن أن تتشابه مع سلوك التربة . 
  - تقسيم القاعدة و علاقته بالحمل المركز القادم من العامود : 
في الحقيقة الحمل المركز القادم من العامود هو حمل غير مركز , و لكن حمل موزع على مساحة العامود ( 0.4X0.4) لذلك عند تقسيم القاعدة ينبغي أن نحاول تقسيمها بحيث نحصل على مساحات كاملة تحت مكان العامود , و يمثل الشكل التالي مثالا على تقسيم غير منتظم للقاعدة حتى تتناسب مع المساحة المحملة:
- لتقسيم القاعدة في مثالنا , قم بتحديد القاعدة , ثم من قائمة Edit >> Edit shell اختر Divide shell , لتظهر نافذة divide selected shell . 

- سوف نقوم بالتقسيم إلى 5 أجزاء بكل اتجاه , كل جزء منها بطول 0.4 ( و يمكن تقسيمها كل 0.2 للحصول على 10 اجزاء بكل اتجاه ). 
قم بتحديد الخيار Divide Quadrilaterals و إدخل في كلا المربعين الرقم 5 . 

لتحصل على قاعدة مقسمة كما في الصورة 
تعيين الإسناد 

- من قائمة Define >> Spring properties ثم أختر Point spring , في الشاشة التي تظهر أختر add new property لتظهر النافذة Point spring properties data .
هذه النافذة تحتوي على Simple Spring stiffness in global directions , بتعيين كافة القيم إلى الصفر بإستثناء قيمة Transition in Z , سوف نقوم بحسابها من خلال العلاقة F=Ku . 
هنا سوف نحتاج إلى قوة تحمل التربة Bearing capacity , و يمكن الحصول عليها من تقرير مختبر التربة لطبقة التأسيس المطلوبة .  و بإستعمال أحد معادلات حساب ال settlement S لحساب قيمة u ( يمكن إعتبار قيمة S في معظم الاحيان كـ 0.1 cm (
يمكن استعمال المعادلة الآتية لحساب S : ( و تستعمل هذه المعادلة فقط لحساب ال settlement تحت مركز القاعدة فقط ) 
حيث :
q : bearing capacity 
B : footing width
v : Poisson ratio
Es : Modulus of elasticity
Iw : Influance factor , depending on the shape and rigidity of footing (0.95)

و يمكن بدلا من حساب قيمة K , إستعمال بعض القيم الإفتراضية لمختلف أنواع التربة كالأتي : 
losse sand : 4800 - 16000 kN/m
Medium dense sand : 4600 - 80000 kN/m
Dense sand : 6400 - 128000 kN/m
Clayey medium dense sand : 32000-80000 kN/m
Silty medium dense sand : 24000 - 48000
Clay soil 
  - q <= 200 kPa : 12000-24000 kN/m
  - 200 < q <=800 : 24000 - 48000 kN/m
  - q > = 800 : K > 48000 kN/m

في مثالنا سوف نعتبر تربة التأسيس من نوع silty clay , من تقرير مختبر التربة حصلنا على : 
Es = 600 Kg/cm.sqr
q = 2.5 kg/cm.sqr ( 25000 kg/m.sqr = 245 kN/m.sqr) 
v = 0.15 
بالتعويض في المعادلة السابقة نحصل على قيمة S = 0.774 cm أو ( 0.00774 m ) .


- تختلف المساحة التي يدعمها كل نابض , لذلك سوف نقوم بتعريف ثلاث قيم كما في الشكل  :
من قائمة define >> spring properties أختر point spring , قم بتعريف ثلاث نوابض بإستعمال قيم K في الصورة السابقة و ذلك بإدخالها في مربع Transition in Z كما في الصورة : 




بذلك قمنا بتعريف spring property

- إضافة نوابض إلى القاعدة . 
* قم بتحديد كل ال nodes التي تمتلك نفس ال K ( مثلا قم يتحديد نقاط الزوايا فقط التي لها K1 و من قائمة assign >> joint أختر spring ثم أختر spring corner K1 كرر العلمية لباقي النوابض 


* فيقوم البرنامج بتعييبن إسناد من نوع spring لجميع النقاط كما في الصورة 


7. تعيين حمل العامود 

لتعيين حمل العامود على شكل Area load ( مساحة العامود = 0.4*0.4 = 0.16 متر مربع , فالحمل الموزع على هذه المساحة يكون 400/0.16 = 2500 kN/m , قم بتحديد المساحة في منتصف القاعدة كما في الشكل 
من قائمة Assign>> shell load اختر Uniform , ثم أدخل القيمة 2500 

استعراض نتائج التحليل : 

يمكن الحصول على نتائج أدق و أفضل من خلال تقليل مساحة ال mesh 





الخميس، 11 أغسطس 2016

مشكلة إنتهاء ترخيص برامج CSI.

من المشاكل الشائعة في مجموعة برامج CSI هي تعذر تحميل الرخصة  , حيث تظهر نافذة عند تشغيل البرنامج و عند إغلاقها فإن البرنامج يتوقف.
و المتسبب بهذه المشكلة هو إنتهاء ترخيص البرنامج , أو أن الكراك المستعمل قديم.
تعتمد برامج CSI في ترخيصها على ملف يسمى Lservrc موجود في مسار البرنامج , يحتوي هذا الملف على تاريخ أنتهاء الترخيص , لذلك عندما يكون الكراك قديما فإنه سوف يولد تاريخا يأتي بعد التاريخ الفعلي لترخيص البرنامج.

و يظهر في الصورة ملف Lservrc , و يحتوي على تاريخ الأنتهاء . 
و بعد معرفة تاريخ إنتهاء البرنامج فإنه يمكنك تغيير تاريخ الجهاز ليكون قبل تاريخ الإنتهاء , أو بإمكانك إستعمال أحد الأدوات التي تقوم بتشغيل البرامج الأخرى بالتاريخ الذي ترغب دون الحاجة لتغيير وقت الجاهز مثل Run As Date ( أو لنسخة 64 هنا ).

يمكن فتح الملف بإستعمال ال Notepad .
توجه إلى مسار تثبيت البرنامج المطلوب مثلا Program Files >> Computer and structure >> CSiBridge و ابحث عن ملف Lservrc .

أضغط على الملف بزر الماوس الأيمن و أختر Open , لتظهر النافذة الآتية :

أختر منها Notepad .
لتظهر محتويات الملف على الشاشة .
حرك شريط التصفح حتى تصل إلى الأسطر التي تحتوي كلمة Expires , ثم أقرأ التاريخ الموجود .


الجمعة، 29 يوليو 2016

Shell and membrane in ETABS

يوفر برنامج ETABS و برامج ال Finite elements بشكل عام , ما يسمى بالعناصر (القشرية Shell) و العناصر ( الغشائية Membrane) , و يختلف سلوك كل من هما عن الآخر. و يستعمل برنامج ETABS هذه العناصر لتمثيل الأسقف Slabs و الجدران Walls , و يتم إختيارهما عند تعيين مقطع خرساني للسقف أو الجدار.

يأتي الأختلاف الأساسي بين ال Shell و ال Membrane من طريقة مقاومة و نقل كل منهما للأحمال المسلطة عليهما.

خصائص العناصر القشرية Shell 
- تمتلك العناصر القشرية  In plane , Out of plane stiffness . و نعني ب In و Out, أتجاه تحميل العنصر , ففي الأسقف يتم تحميلها بحمل موزع بشكل عامودي على مستوى السقف , و في الجدران الأستنادية يتم تحميلها بحمل يمثل ضغط التربة عامودي على مستوى الجدار , أما في حالة الجدران الحاملة Bearing walls فيتم تحميلها بحمل موازي لمستواها .
لذلك العناصر القشرية ملائمة لتمثيل هذه الأسقف و الجدران , لأنها تستطيع مقاومة كل من Bending و Axial .
- تحتاج لعمل Mesh للعنصر, و ذلك حتى يمتلك العنصر bending stiffness ملائمة , مما يمكن هذه العناصر من مقاومة جزء من القوى المنقولة إلى العناصر الحاملة مثل الجسور Beams . لذلك يكون الحمل المنقول من عنصر قشري إلى جسر أقل من الحمل المنقول من عنصر غشائي إلى نفس الجسر.
- إمكانية عمل Mesh في إتجاه واحد لتمثيل One way slab.
- الأقرب إلى الواقع (من الغشائية) في تمثيل الأسقف و الجدران .

لتعيين مقطع خرساني كعنصر قشري نتبع الخطوات التالية . ( بإستخدام نسخة 9.7)
1. من قائمة Define نختار Wall/slab/Deck section
2. من النافذة نختر Add slab أو Add wall.
3. لتظهر نافذة Wall/slab section و منها نختار المقطع خرساني CONC و السماكة Thickness , و في التبويب Type نختار Shell.
سلوك العنصر القشري . 
لتوضيح سلوك العنصر القشري قم برسم مبنى يتكون من طابق واحد و يرتكز على أربعة أعمدة و أربعة جسور beams . 




لإضافة حمل إلى السقف , قم بتحديد السقف , ثم من قائمة Assign نختار shell/area Load ثم Uniform لتظهر نافذة Uniform Surface Load
قم بعمل Run

و لمشاهدة تأثيرات الأحمال , من قائمة Display نختار show deformed shape و نحدد حالة التحميل Live

لاحظ أن التشوه عكس المتوقع , فالمفترض بالجسور أن تنحني بتأثير ال bending moment المنتقل من السقف , كما يفترص بالأعمدة أن تتشوه بتأثير ال buckling , و لكن ما حدث أن كامل السقف مع الجسور قد أنزاح إلى الأسفل و سبب تشوه الأعمدة بالأتجاه الرأسي فقط.

و من قائمة Display أختر Show member Forces/stress diagram ثم أختر Support/spring reaction , ثم حدد حالة التحميل Live

ليقوم البرنامج بعرض قيم ردود الفعل عند المساند . و هي أيضا عكس المتوقع , حيث أن ال Reaction الوحيد هو بإتجاه z في حين باقي القيم x,y عبارة عن أصفار. 

ملاحظة
_________________________________________________________________________________
هذه القيمة - في هذه الحالة فقط - , عبارة عن حاصل ضرب (الحمل الموزع * مساحة السقف) مقسوما على عدد الاعمدة , و يمكن معرفة مساحة العناصر و الأحمال الموجودة عليه بالضغط بزر الماوس الأيمن على العنصر , لتظهر نافذة Area information و من تبويب Location يمكن الحصول على المساحة Area و من تبويب Loads يمكن الحصول على الأحمال . 
10*36 \ 4 = 90 Kn

_________________________________________________________________________________

و حتى تظهر هذه النافذة يجب مراعاة أن البرنامج لا يقوم بعرض أي من Deformed shape أو Forces , ففي حالة تفعيل عرض ال Forces فإن الضغط بزر الماوس الأيمن سوف يظهر نافذة تحتوي على تفاصيل ما هو معروض , مثلا عند عرض ال Moment فإنه سوف يعرض تفاصيل ال Moment , Shear , axial .
و لعرض قيم ال moment  من قائمة Display أختر Show member Forces/stress diagram ثم أختر Frame/Pier/spandrel forces , ثم أختر Moment-33  ثم حدد حالة التحميل Live .
في هذه الحالة لن يظهر أن شيء على الشاشة , قم بالضغط بزر الماوس الأيمن على أي من الجسور  لتظهر لك شاشة Diagram for , و نفس الأمر للأعمدة. 



لماذا لم تظهر على الجسور أي قوى ؟
ذكرنا أنه من خصائص ال Shell انها بحاجة لعمل Mesh حتى تستطيع نقل الأحمال بشكل صحيح
يتعامل البرنامج مع العنصر الخطي بين كل Node و الأخرى على أنها عنصر لا يمكنه التشوه و الأنزياح , فالذي ينزاح هي ال Nodes و ليس العناصر الرابطة بينها , في حالتنا فإن جوانب السقف تربط بين ال Nodes المرتكزة على الأعمدة , لذلك فالاحمال سوف تنتقل مباشرة إلى الأعمدة , في حين انه في منطقة ألتقاءه بالجسر فإن السقف لا يمتلك أي Nodes يمكنها نقل الاحمال إلى الجسر , لذلك فإن عمل Mesh سوف يوفر المزيد من ال Nodes ترتكز على الجسر . 
و لعمل Mesh للسقف : قم بتحديد السقف , ثم من قائمة Edit أختر Mesh Areas  لتظهر شاشة Mesh Selected Areas 
ثم قم بأختيار عدد ال Meshes المطلوبة , و كلما زاد عددها كلما زادت دقة النتائج في المقابل يزدا الوقت المستغرق في التحليل


قم بعمل Run .

ثم لاحظ شكل التشوه الناتج . حيث أن السقف أصبح يمتلك Nodes يمكنه التشوه من خلالها (إن جاز التعبير ) و كل Node منها عليها حمل مختلف عن الأخرى . 

و يمكن الحصول على قيمة إنزياح Displacement لكل Node بالضغط بر الماوس الأيمن على ال Node 

قم بعرض قيم ال moment  من قائمة Display أختر Show member Forces/stress diagram ثم أختر Frame/Pier/spandrel forces , ثم أختر Moment-33  ثم حدد حالة التحميل Live .

بالضغط بزر الماوس الأيمن على أحد الجسور , تظهر نافذة Diagram for beam .

في تبويب Dist Load لاحظ وجود أحمال مركزة Concentrated  على الجسر , مصدر هذه الأحمال هي ال Nodes الخاصة بال Shell.
كما يمكن عرض ال Stresses في الShell على شكل خريطة كنتورية ملونة , من قائمة Display ختر Show member Forces/stress diagram ثم Shell stresses/forces . و من النافذة التي تظهر يمكن الأختيار بين Forces او Stresses 


و لجعل لعنصر القشري يعمل كـ One way element نقوم تقسيم ال Mesh بإتجاه واحد فقط . و للقيام بذلك قم بحذف كل ال Meshes التي نتجت من التقسيم في الخطوات السايقة , و ذلك بتحديدها من خلال القائمة Select ثم إختيار by Area Object type لتظهر نافذة نختار منها Floor . ثم أضغط delete. ثم قم برسم سقف جديد . و إضافة حمل موزع عليه.


قم بتحديد السقف الجديد , ثم من قائمة Edit أختر Mesh Areas . و في النافذة التي تظهر قم بتغيير أحد الأرقام إلى 1 . (حيث يمثل المربع الأول عدد ال meshes بإتجاه X و الآخر بإتجاه Y ).


قم بعمل Run . و لاحظ النتائج .
شكل التشوه :