ĐẠI HỌC QUỐC GIA TP HCM TRƯỜNG ĐẠI HỌC BÁCH KHOA NGUYỄN TẤN PHÁT PHÂN TÍCH PHI TUYẾN CỘT ỐNG THÉP NHỒI BÊ TÔNG DƯỚI TÁC ĐỘNG CỦA TẢI TRỌNG VÀ NHIỆT ĐỘ Chu yên n gành : Xây d ựng cơng trình dân dụng cơng nghiệp Mã số : 60.58.20 LUẬN VĂN THẠC SĨ TP HỒ CHÍ MINH, tháng năm 2013 CÔNG TRÌNH ĐƯC HOÀN THÀNH TẠI TRƯỜNG ĐẠI HỌC BÁCH KHOA-ĐHQG TP HCM Người hướng dẫn khoa học: Người chấm nhận xét 1: Người chấm nhận xét 2: Luận văn thạc sĩ bảo vệ Trường Đại học Bách Khoa, ĐH Quốc gia Tp HCM ngày tháng năm Thành phần Hội đồng đánh giá luận văn thạc sĩ gồm: Xác nhận Chủ tịch Hội đồng đánh giá luận văn Trưởng Khoa quản lý chuyên ngành sau luận văn sửa chữa (nếu có) CHỦ TỊCH HỘI ĐỒNG TRƯỞNG KHOA ĐẠI HỌC QUỐC GIA TP.HCM TRƯỜNG ĐẠI HỌC BÁCH KHOA CỘNG HÒA XÃ HỘI CHỦ NGHĨA VIỆT NAM Độc lập - Tự - Hạnh phúc NHIỆM VỤ LUẬN VĂN THẠC SĨ Họ tên học viên : NGUYỄN TẤN PHÁT MSHV : 11210242 Ngày, tháng, năm sinh : 26/09/1986 Nơi sinh : Bà Rịa – Vũng Tàu Chuyên ngành: Xây dựng Cơng trình DD & CN Mã số : 605820 I TÊN ĐỀ TÀI: Phân tích phi tuyến cột ống thép nhồi bê tông tác động tải trọng nhiệt độ II NHIỆM VỤ VÀ NỘI DUNG: - Xây dựng phần tử hữu hạn dầm-cột mơ ứng xử phi tuyến hình học vật liệu cột ống thép nhồi bê tông tác dụng tải trọng nhiệt độ - Thiết lập thuật tốn giải phi tuyến để phân tích kết cấu có kể đến giai đoạn gia tăng nhiệt ngồi giai đoạn gia tăng tải trọng phát triển chương trình máy tính nhằm tự động hóa q trình phân tích ngơn ngữ lập trình C++ - Kiểm tra độ tin cậy chương trình cách so sánh kết đạt với kết nghiên cứu sẵn có tác giả khác giới - Rút nhận xét kết luận khối lượng công việc thực được; Đề xuất hướng phát triển đề tài III NGÀY GIAO NHIỆM VỤ : 02/07/2012 IV NGÀY HOÀN THÀNH NHIỆM VỤ : 30/11/2012 V : TS TRƯƠNG HỒI CHÍNH CÁN BỘ HƯỚNG DẪN TS NGÔ HỮU CƯỜNG Tp HCM, ngày tháng năm 20 CÁN BỘ HƯỚNG DẪN (Họ tên chữ ký) BAN QUẢN LÝ CHUYÊN NGÀNH (Họ tên chữ ký) TS TRƯƠNG HỒI CHÍNH TS NGƠ HỮU CƯỜNG TRƯỞNG KHOA KỸ THUẬT XÂY DỰNG (Họ tên chữ ký) Lời Cảm Ơn “ Công cha nghĩa mẹ ơn thầy” Lời em xin tỏ lòng biết ơn sâu sắc đến thầy hướng dẫn luận án, thầy TS Trương Hồi Chính thầy TS Ngơ Hữu Cường, người thầy mẫu mực uyên bác, người cố vấn đầy kinh nghiệm, tận tình hướng dẫn, định hướng khoa học động viên tinh thần cho em vượt qua khó khăn suốt q trình nghiên cứu Đạo đức tri thức thầy gương sáng cho em noi theo Tôi xin cảm ơn Ban giám hiệu, Phòng Quản lý khoa học Sau đại học, Khoa Kỹ thuật xây dựng, thầy cô giảng dạy cao học Trường Đại học Bách khoa TP Hồ Chí Minh tất thầy cô dạy từ trước đến kiến thức quý báu truyền đạt Một lòng biết ơn vô hạn xin gửi tới cha mẹ Cha mẹ dạy cho điều hay lẽ phải, an ủi động viên nâng đỡ con, giúp nghị lực để vượt qua khó khăn sống Cha mẹ gương sáng để học tập noi theo Tôi xin cảm ơn Giám đốc đồng nghiệp Công ty cổ phần tư vấn kỹ thuật xây dựng Kiến An, xin cảm ơn tất bạn bè và người thân động viên tạo điều kiện giúp hoàn thành luận án Xin chân thành cảm ơn! Tp.HCM, ngày 26 tháng 11 năm 2012 Nguyễn Tấn Phát TĨM TẮT LUẬN VĂN THẠC SĨ PHÂN TÍCH PHI TUYẾN CỘT ỐNG THÉP NHỒI BÊ TÔNG DƯỚI TÁC ĐỘNG CỦA TẢI TRỌNG VÀ NHIỆT ĐỘ Tóm tắt Luận án nghiên cứu ứng xử phi tuyến cấu kiện cột ống thép nhồi bê tông tác dụng tải trọng nhiệt độ có kể đến ứng xử phi tuyến hình học vật liệu Ma trận độ cứng phần tử hữu hạn phi tuyến thiết lập việc áp dụng nguyên lý toàn phần dừng có kể đến tác động tải trọng nhiệt Tiết diện gồm ống thép lõi bê tông phần tử hữu hạn chia thành nhiều thớ quan hệ ứng suất biến dạng thớ thép bê tơng cập nhật suốt q trình phân tích để mơ tác động phi đàn hồi qua mặt cắt ngang dọc theo chiều dài cấu kiện Một chương trình phân tích ngơn ngữ lập trình C++ dùng thuật tốn giải phi tuyến Euler đơn giản hiệu phát triển kết so sánh với kết nghiên cứu sẵn có khác để minh họa độ tin cậy chương trình phân tích đề xuất NONLINEAR ANALYSIS OF CONCRETE FILLED STEEL TUBE UNDER MECHANICAL AND FIRE LOADINGS Abstract This thesis presents a numerical procedure for analysis of concrete-filled steel tube column members under mechanical and fire loading considering geometric and material nonlinear effects The stiffness matrix of the nonlinear finite element is developed by the application of the principle of stationary total energy The crosssection comprised of steel tube and concrete core at the element mid-length is divided into many fibers and the strain-stress relation of each fiber is updated during the analysis process to simulate the inelastic effects thoughout the cross-section and along the member length An analysis program written by C++ programming language using the simple but effective Euler nonlinear algorithm is developed and its results are compared with those of existing studies to verify the reliability of the proposed analysis program Mục lục DANH MỤC HÌNH VẼ DANH MỤC BẢNG BIỂU DANH MỤC CÁC KÍ HIỆU VÀ CHỮ VIẾT TẮT CHƯƠNG TỔNG QUAN 11 1.1 GIỚI THIỆU 11 1.2 TÌNH HÌNH NGHIÊN CỨU ĐỀ TÀI 11 1.2.1 Tình hình nghiên cứu giới 11 1.2.2 Tình hình nghiên cứu nước 13 1.3 MỤC TIÊU CỦA LUẬN ÁN 13 1.4 Ý NGHĨA THỰC TIỄN CỦA ĐỀ TÀI 13 CHƯƠNG CƠ SỞ LÝ THUYẾT 14 2.1 MƠ HÌNH VẬT LIỆU 14 2.1.1 Mơ hình vật liệu bê tông 14 2.1.2 Mơ hình vật liệu thép 15 2.2 ẢNH HƯỞNG CỦA NHIỆT ĐỘ LÊN VẬT LIỆU 16 2.2.1 Ảnh hưởng nhiệt độ lên vật liệu thép 16 2.2.2 Ảnh hưởng nhiệt độ lên vật liệu bê tông 19 2.3 SỰ PHÂN PHỐI NHIỆT TRONG TIẾT DIỆN 21 2.3.1 Phân phối nhiệt vỏ thép 21 2.3.2 Phân phối nhiệt lõi bê tông 22 2.4 TIÊU CHUẨN ISO-834 25 2.5 MƠ HÌNH PHẦN TỬ HỮU HẠN 26 2.5.1 Các giả thiết 26 2.5.2 Phương pháp RayLeigh-Ritz 26 2.5.3 Thế lực tác dụng 35 2.5.4 Nguyên lý toàn phần dừng 35 2.6 CÁC PHƯƠNG TRÌNH CÂN BẰNG GIA TĂNG 60 CHƯƠNG CHƯƠNG TRÌNH ỨNG DỤNG 62 3.1 Giới thiệu 62 3.2 Chi tiết q trình phân tích 62 3.2.1 Mô hình phần tử 62 3.2.2 Xác định trạng thái phần tử thớ 64 3.2.3 Chuyển từ hệ trục tọa độ địa phương sang toàn cục 67 3.3 LƯU ĐỒ GIẢI THUẬT 67 3.3.1 Vòng lặp gia tăng tải đứng 68 3.3.2 Vòng lặp gia tăng nhiệt độ 69 3.3.3 Vòng lặp gia tải ngang 70 3.4 THUẬT GIẢI PHI TUYẾN 71 3.4.1 Giới thiệu 71 3.4.2 Thuật toán Euler 71 3.4.3 Kỹ thuật Newton-Raphson 72 3.4.4 Kỹ thuật điều chỉnh chuyển vị 74 3.4.5 Kết luận 75 3.5 CHƯƠNG TRÌNH ỨNG DỤNG 75 CHƯƠNG VÍ DỤ ÁP DỤNG CHƯƠNG TRÌNH CFT-FIRE 78 4.1 GIỚI THIỆU 78 4.2 BÀI TỐN GIA TẢI Ở NHIỆT ĐỘ PHỊNG 78 4.2.1 Bài toán Sangdo Hong 78 4.2.2 Bài toán dầm cột Shakir 80 4.3 BÀI TOÁN GIA TẢI VÀ NHIỆT KẾT HỢP CỦA SANGDO HONG 82 4.4 KHẢO SÁT SỰ ẢNH HƯỞNG CỦA NHIỆT ĐỘ 93 4.4.1 Mẫu CFT 254x254x8-16.4 93 4.4.2 Mẫu CFT 305x305x8-16.4 94 4.5 KHẢO SÁT SỰ HỘI TỤ CỦA CHƯƠNG TRÌNH 95 CHƯƠNG KẾT LUẬN 97 5.1 TÓM TẮT LUẬN VĂN 97 5.2 KẾT LUẬN VÀ HƯỚNG PHÁT TRIỂN CỦA ĐỀ TÀI 98 5.2.1 Kết luận 98 5.2.2 Hướng phát triển đề tài 98 TÀI LIỆU THAM KHẢO 100 DANH MỤC HÌNH VẼ Hình 2.1 Quan hệ ứng suất-biến dạng bê tông theo Euro code 2: “Design of concrete struture” 14 Hình 2.2 Quan hệ ứng suất-biến dạng thép có xét tái bền 16 Hình 2.3 Biểu đồ quan hệ hệ số giảm cường đô (K) thép theo nhiệt độ 17 Hình 2.4 Hệ số giảm cường độ (K) thép theo Herzt, Xu, Eurocode 18 Hình 2.5 Mơ hình độ giảm cường độ bê tông theo nhiệt độ 20 Hình 2.6 Quan hệ λ-T vật liệu thép 21 Hình 2.7 Quan hệ λ-T vật liệu bê tông 22 Hình 2.8 Phân phối nhiệt theo khoảng cách A.Y NASSIF 22 Hình 2.9 Phân phối nhiệt theo chiều cao-R.Kichuki 23 Hình 2.10 Phân phối nhiệt theo khoảng cách –G.M Freskakis 23 Hình 2.11 Phân phối nhiệt theo khoảng cách –Josephine V.C 24 Hình 2.12 Phân phối nhiệt theo khoảng cách –Euro code 24 Hình 2.13 Phân phối nhiệt theo khoảng cách –Ellingwood 25 Hình 2.14 Đường cong quan hệ nhiệt độ thời gian theo ISO-834 26 Hình 3.2 Phần tử hữu hạn dầm cột CFT 26 Hình 3.3 Biểu đồ quan hệ ứng suất-biến dạng bê tơng thép 28 Hình 3.4 Phần tử điển hình chịu tải trọng 35 Hình 4.1 Phần tử hữu hạn dùng mô chảy dẻo theo chiều dài cấu kiện 62 Hình 4.2 Chia phần tử hứu hạn thành phần tử 63 Hình 4.3 Chia thớ tiết diện cột CFT 63 Hình 4.4 Lưu đồ gia tải đứng 68 Hình 4.5 Lưu đồ gia nhiệt 69 Hình 4.6 Lưu đồ gia tải ngang 70 Hình 3.5 Ứng xử tải trọng-chuyển vị khung cổng chịu tải trọng phân bố tải trọng ngang 71 Hình 3.6 Sơ đồ minh họa thuật toán Euler đơn giản 72 Hình 3.7 Sơ đồ minh họa kỹ thuật Newton-Raphson 73 Hình 3.8 Sơ đồ minh họa kỹ thuật Newton-Raphson hiệu chỉnh 74 Hình 4.7 Các giai đoạn q trình phân tích 78 Hình 4.8 Quan hệ Lực-Chuyển vị ngang (CFT-254x254x8-16.4-20) 79 Hình 4.9 Quan hệ Moment -độ cong (CFT-254x254x8-16.4-20) 80 Hình 4.10 Sơ đồ hình học phần trăm chảy dẻo (254x254x8-16-20) 80 Hình 4.11 Sơ đồ phần tử quan hệ lực-chuyển vị mẫu (Shakir) 81 Hình 4.12 Sơ đồ phần tử quan hệ lực-chuyển vị mẫu (Shakir) 82 Hình 4.13 Các giai đoạn q trình phân tích 83 Hình 4.14 Quan hệ tải ngang chuyền vị ngang (CFT-254x254x8-16.4-300) 84 Hình 4.15 Quan hệ Moment -độ cong (CFT-254x254x8-16.4-300) 84 Hình 4.16 Sơ đồ hình học phần trăm chảy dẻo (CFT-254x254x8-16.4-300) 85 Hình 4.17 Quan hệ tải ngang chuyền vị ngang (CFT-254x254x8-16-500) 85 Hình 4.18 Quan hệ Moment -độ cong (CFT-254x254x8-16-500) 86 Hình 4.19 Sơ đồ hình học phần trăm chảy dẻo (CFT-254x254x8-16-500) 86 Hình 4.20 Quan hệ Moment -độ cong (CFT-254x254x8-17-500) 87 Hình 4.21 Quan hệ Moment -độ cong (CFT-254x254x8-17-500) 87 Hình 4.22 Quan hệ Moment -độ cong (CFT-254x254x8-36-300) 88 Hình 4.23 Sơ đồ hình học phần trăm chảy dẻo (CFT-254x254x8-36-300) 88 Hình 4.24 Quan hệ Moment -độ cong(CFT-254x254x8-36-500) 89 Hình 4.25 Sơ đồ hình học phần trăm chảy dẻo (CFT-254x254x8-36-500) 89 CFT-Fire Code ktc216= -(6*Mc_n[i*N_SELE + k]*anpC*del_T*tc2)/sel; ktc221= 0; ktc222= (18*Nc_n[i*N_SELE + k]*(anpC*anpC)*(del_T*del_T)*tc2)/(5*sel) + (12*Nc_n[i*N_SELE + k]*anpC*del_T*tc1)/(5*sel) - (72*Ic_ze[i*N_SELE + k]*anpC*del_T*ftb_c*tc2)/(sel*sel*sel); ktc223= (3*Nc_n[i*N_SELE + k]*(anpC*anpC)*(del_T*del_T)*tc2)/10 + (Nc_n[i*N_SELE + k]*anpC*del_T*tc1)/5 - (36*Ic_ze[i*N_SELE + k]*anpC*del_T*ftb_c*tc2)/(sel*sel); ktc224= 0; ktc225= (72*Ic_ze[i*N_SELE + k]*anpC*del_T*ftb_c*tc2)/(sel*sel*sel) - (12*Nc_n[i*N_SELE + k]*anpC*del_T*tc1)/(5*sel) - (18*Nc_n[i*N_SELE + k]*(anpC*anpC)*(del_T*del_T)*tc2)/(5*sel); ktc226= (3*Nc_n[i*N_SELE + k]*(anpC*anpC)*(del_T*del_T)*tc2)/10 + (Nc_n[i*N_SELE + k]*anpC*del_T*tc1)/5 - (36*Ic_ze[i*N_SELE + k]*anpC*del_T*ftb_c*tc2)/(sel*sel); ktc231= (6*Mc_n[i*N_SELE + k]*anpC*del_T*tc2)/sel; ktc232= (3*Nc_n[i*N_SELE + k]*(anpC*anpC)*(del_T*del_T)*tc2)/10 + (Nc_n[i*N_SELE + k]*anpC*del_T*tc1)/5 - (36*Ic_ze[i*N_SELE + k]*anpC*del_T*ftb_c*tc2)/(sel*sel); ktc233= (4*sel*Nc_n[i*N_SELE + k]*anpC*del_T*tc1)/15 + (2*sel*Nc_n[i*N_SELE + k]*(anpC*anpC)*(del_T*del_T)*tc2)/5 - (24*Ic_ze[i*N_SELE + k]*anpC*del_T*ftb_c*tc2)/sel; ktc234= -(6*Mc_n[i*N_SELE + k]*anpC*del_T*tc2)/sel; ktc235= (36*Ic_ze[i*N_SELE + k]*anpC*del_T*ftb_c*tc2)/(sel*sel) - (Nc_n[i*N_SELE + k]*anpC*del_T*tc1)/5 - (3*Nc_n[i*N_SELE + k]*(anpC*anpC)*(del_T*del_T)*tc2)/10; ktc236= - (sel*Nc_n[i*N_SELE + k]*anpC*del_T*tc1)/15 - (sel*Nc_n[i*N_SELE + k]*(anpC*anpC)*(del_T*del_T)*tc2)/10 - (12*Ic_ze[i*N_SELE + k]*anpC*del_T*ftb_c*tc2)/sel; ktc241= (6*Nc_n[i*N_SELE + k]*anpC*del_T*tc2)/sel; ktc242= 0; ktc243= -(6*Mc_n[i*N_SELE + k]*anpC*del_T*tc2)/sel; ktc244= -(6*Nc_n[i*N_SELE + k]*anpC*del_T*tc2)/sel; ktc245= 0; ktc246= (6*Mc_n[i*N_SELE + k]*anpC*del_T*tc2)/sel; ktc251= 0; ktc252= (72*Ic_ze[i*N_SELE + k]*anpC*del_T*ftb_c*tc2)/(sel*sel*sel) - (12*Nc_n[i*N_SELE + k]*anpC*del_T*tc1)/(5*sel) - (18*Nc_n[i*N_SELE + k]*(anpC*anpC)*(del_T*del_T)*tc2)/(5*sel); ktc253= (36*Ic_ze[i*N_SELE + k]*anpC*del_T*ftb_c*tc2)/(sel*sel) - (Nc_n[i*N_SELE + k]*anpC*del_T*tc1)/5 - (3*Nc_n[i*N_SELE + k]*(anpC*anpC)*(del_T*del_T)*tc2)/10; ktc254= 0; ktc255= (18*Nc_n[i*N_SELE + k]*(anpC*anpC)*(del_T*del_T)*tc2)/(5*sel) + (12*Nc_n[i*N_SELE + k]*anpC*del_T*tc1)/(5*sel) - (72*Ic_ze[i*N_SELE + k]*anpC*del_T*ftb_c*tc2)/(sel*sel*sel); ktc256= (36*Ic_ze[i*N_SELE + k]*anpC*del_T*ftb_c*tc2)/(sel*sel) - (Nc_n[i*N_SELE + k]*anpC*del_T*tc1)/5 - (3*Nc_n[i*N_SELE + k]*(anpC*anpC)*(del_T*del_T)*tc2)/10; ktc261= -(6*Mc_n[i*N_SELE + k]*anpC*del_T*tc2)/sel; ktc262= (3*Nc_n[i*N_SELE + k]*(anpC*anpC)*(del_T*del_T)*tc2)/10 + (Nc_n[i*N_SELE + k]*anpC*del_T*tc1)/5 - (36*Ic_ze[i*N_SELE + k]*anpC*del_T*ftb_c*tc2)/(sel*sel); ktc263= - (sel*Nc_n[i*N_SELE + k]*anpC*del_T*tc1)/15 - (sel*Nc_n[i*N_SELE + k]*(anpC*anpC)*(del_T*del_T)*tc2)/10 - (12*Ic_ze[i*N_SELE + k]*anpC*del_T*ftb_c*tc2)/sel; ktc264= (6*Mc_n[i*N_SELE + k]*anpC*del_T*tc2)/sel; ktc265= (36*Ic_ze[i*N_SELE + k]*anpC*del_T*ftb_c*tc2)/(sel*sel) - (Nc_n[i*N_SELE + k]*anpC*del_T*tc1)/5 - (3*Nc_n[i*N_SELE + k]*(anpC*anpC)*(del_T*del_T)*tc2)/10; ktc266= (4*sel*Nc_n[i*N_SELE + k]*anpC*del_T*tc1)/15 + (2*sel*Nc_n[i*N_SELE + k]*(anpC*anpC)*(del_T*del_T)*tc2)/5 - (24*Ic_ze[i*N_SELE + k]*anpC*del_T*ftb_c*tc2)/sel; ktc311= 0; ktc312= (36*Nc_n[i*N_SELE + k]*anpC*d2*del_T*tc2)/(5*(sel*sel)) + (3*Nc_n[i*N_SELE + k]*anpC*d3*del_T*tc2)/(5*sel) - (36*Nc_n[i*N_SELE + k]*anpC*d5*del_T*tc2)/(5*(sel*sel)) + (3*Nc_n[i*N_SELE + k]*anpC*d6*del_T*tc2)/(5*sel); ktc313= (4*Nc_n[i*N_SELE + k]*anpC*d3*del_T*tc2)/5 - (Nc_n[i*N_SELE + k]*anpC*d6*del_T*tc2)/5 + (3*Nc_n[i*N_SELE + k]*anpC*d2*del_T*tc2)/(5*sel) - (3*Nc_n[i*N_SELE + k]*anpC*d5*del_T*tc2)/(5*sel); ktc314= 0; ktc315= (36*Nc_n[i*N_SELE + k]*anpC*d5*del_T*tc2)/(5*(sel*sel)) - (3*Nc_n[i*N_SELE + k]*anpC*d3*del_T*tc2)/(5*sel) - (36*Nc_n[i*N_SELE + k]*anpC*d2*del_T*tc2)/(5*(sel*sel)) (3*Nc_n[i*N_SELE + k]*anpC*d6*del_T*tc2)/(5*sel); ktc316= (4*Nc_n[i*N_SELE + k]*anpC*d6*del_T*tc2)/5 - (Nc_n[i*N_SELE + k]*anpC*d3*del_T*tc2)/5 + (3*Nc_n[i*N_SELE + k]*anpC*d2*del_T*tc2)/(5*sel) - (3*Nc_n[i*N_SELE + k]*anpC*d5*del_T*tc2)/(5*sel); ktc321= (36*Nc_n[i*N_SELE + k]*anpC*d2*del_T*tc2)/(5*(sel*sel)) + (3*Nc_n[i*N_SELE + k]*anpC*d3*del_T*tc2)/(5*sel) - (36*Nc_n[i*N_SELE + k]*anpC*d5*del_T*tc2)/(5*(sel*sel)) + (3*Nc_n[i*N_SELE + k]*anpC*d6*del_T*tc2)/(5*sel); ktc322= (36*Nc_n[i*N_SELE + k]*anpC*d1*del_T*tc2)/(5*(sel*sel)) - (36*Nc_n[i*N_SELE + k]*anpC*d4*del_T*tc2)/(5*(sel*sel)); ktc323= (3*Nc_n[i*N_SELE + k]*anpC*d1*del_T*tc2)/(5*sel) - (3*Nc_n[i*N_SELE + k]*anpC*d4*del_T*tc2)/(5*sel); ktc324= (36*Nc_n[i*N_SELE + k]*anpC*d5*del_T*tc2)/(5*(sel*sel)) - (3*Nc_n[i*N_SELE + k]*anpC*d3*del_T*tc2)/(5*sel) - (36*Nc_n[i*N_SELE + k]*anpC*d2*del_T*tc2)/(5*(sel*sel)) (3*Nc_n[i*N_SELE + k]*anpC*d6*del_T*tc2)/(5*sel); ktc325= (36*Nc_n[i*N_SELE + k]*anpC*d4*del_T*tc2)/(5*(sel*sel)) - (36*Nc_n[i*N_SELE + k]*anpC*d1*del_T*tc2)/(5*(sel*sel)); ktc326= (3*Nc_n[i*N_SELE + k]*anpC*d1*del_T*tc2)/(5*sel) - (3*Nc_n[i*N_SELE + k]*anpC*d4*del_T*tc2)/(5*sel); ktc331= (4*Nc_n[i*N_SELE + k]*anpC*d3*del_T*tc2)/5 - (Nc_n[i*N_SELE + k]*anpC*d6*del_T*tc2)/5 Page: 53 CFT-Fire Code + (3*Nc_n[i*N_SELE + k]*anpC*d2*del_T*tc2)/(5*sel) - (3*Nc_n[i*N_SELE + k]*anpC*d5*del_T*tc2)/(5*sel); ktc332= (3*Nc_n[i*N_SELE + k]*anpC*d1*del_T*tc2)/(5*sel) - (3*Nc_n[i*N_SELE + k]*anpC*d4*del_T*tc2)/(5*sel); ktc333= (4*Nc_n[i*N_SELE + k]*anpC*d1*del_T*tc2)/5 - 6*Mc_n[i*N_SELE + k]*anpC*d3*del_T*tc2 (4*Nc_n[i*N_SELE + k]*anpC*d4*del_T*tc2)/5; ktc334= (Nc_n[i*N_SELE + k]*anpC*d6*del_T*tc2)/5 - (4*Nc_n[i*N_SELE + k]*anpC*d3*del_T*tc2)/5 - (3*Nc_n[i*N_SELE + k]*anpC*d2*del_T*tc2)/(5*sel) + (3*Nc_n[i*N_SELE + k]*anpC*d5*del_T*tc2)/(5*sel); ktc335= (3*Nc_n[i*N_SELE + k]*anpC*d4*del_T*tc2)/(5*sel) - (3*Nc_n[i*N_SELE + k]*anpC*d1*del_T*tc2)/(5*sel); ktc336= (Nc_n[i*N_SELE + k]*anpC*d4*del_T*tc2)/5 - (Nc_n[i*N_SELE + k]*anpC*d1*del_T*tc2)/5; ktc341= 0; ktc342= (36*Nc_n[i*N_SELE + k]*anpC*d5*del_T*tc2)/(5*(sel*sel)) - (3*Nc_n[i*N_SELE + k]*anpC*d3*del_T*tc2)/(5*sel) - (36*Nc_n[i*N_SELE + k]*anpC*d2*del_T*tc2)/(5*(sel*sel)) (3*Nc_n[i*N_SELE + k]*anpC*d6*del_T*tc2)/(5*sel); ktc343= (Nc_n[i*N_SELE + k]*anpC*d6*del_T*tc2)/5 - (4*Nc_n[i*N_SELE + k]*anpC*d3*del_T*tc2)/5 - (3*Nc_n[i*N_SELE + k]*anpC*d2*del_T*tc2)/(5*sel) + (3*Nc_n[i*N_SELE + k]*anpC*d5*del_T*tc2)/(5*sel); ktc344= 0; ktc345= (36*Nc_n[i*N_SELE + k]*anpC*d2*del_T*tc2)/(5*(sel*sel)) + (3*Nc_n[i*N_SELE + k]*anpC*d3*del_T*tc2)/(5*sel) - (36*Nc_n[i*N_SELE + k]*anpC*d5*del_T*tc2)/(5*(sel*sel)) + (3*Nc_n[i*N_SELE + k]*anpC*d6*del_T*tc2)/(5*sel); ktc346= (Nc_n[i*N_SELE + k]*anpC*d3*del_T*tc2)/5 - (4*Nc_n[i*N_SELE + k]*anpC*d6*del_T*tc2)/5 - (3*Nc_n[i*N_SELE + k]*anpC*d2*del_T*tc2)/(5*sel) + (3*Nc_n[i*N_SELE + k]*anpC*d5*del_T*tc2)/(5*sel); ktc351= (36*Nc_n[i*N_SELE + k]*anpC*d5*del_T*tc2)/(5*(sel*sel)) - (3*Nc_n[i*N_SELE + k]*anpC*d3*del_T*tc2)/(5*sel) - (36*Nc_n[i*N_SELE + k]*anpC*d2*del_T*tc2)/(5*(sel*sel)) (3*Nc_n[i*N_SELE + k]*anpC*d6*del_T*tc2)/(5*sel); ktc352= (36*Nc_n[i*N_SELE + k]*anpC*d4*del_T*tc2)/(5*(sel*sel)) - (36*Nc_n[i*N_SELE + k]*anpC*d1*del_T*tc2)/(5*(sel*sel)); ktc353= (3*Nc_n[i*N_SELE + k]*anpC*d4*del_T*tc2)/(5*sel) - (3*Nc_n[i*N_SELE + k]*anpC*d1*del_T*tc2)/(5*sel); ktc354= (36*Nc_n[i*N_SELE + k]*anpC*d2*del_T*tc2)/(5*(sel*sel)) + (3*Nc_n[i*N_SELE + k]*anpC*d3*del_T*tc2)/(5*sel) - (36*Nc_n[i*N_SELE + k]*anpC*d5*del_T*tc2)/(5*(sel*sel)) + (3*Nc_n[i*N_SELE + k]*anpC*d6*del_T*tc2)/(5*sel); ktc355= (36*Nc_n[i*N_SELE + k]*anpC*d1*del_T*tc2)/(5*(sel*sel)) - (36*Nc_n[i*N_SELE + k]*anpC*d4*del_T*tc2)/(5*(sel*sel)); ktc356= (3*Nc_n[i*N_SELE + k]*anpC*d4*del_T*tc2)/(5*sel) - (3*Nc_n[i*N_SELE + k]*anpC*d1*del_T*tc2)/(5*sel); ktc361= (4*Nc_n[i*N_SELE + k]*anpC*d6*del_T*tc2)/5 - (Nc_n[i*N_SELE + k]*anpC*d3*del_T*tc2)/5 + (3*Nc_n[i*N_SELE + k]*anpC*d2*del_T*tc2)/(5*sel) - (3*Nc_n[i*N_SELE + k]*anpC*d5*del_T*tc2)/(5*sel); ktc362= (3*Nc_n[i*N_SELE + k]*anpC*d1*del_T*tc2)/(5*sel) - (3*Nc_n[i*N_SELE + k]*anpC*d4*del_T*tc2)/(5*sel); ktc363= (Nc_n[i*N_SELE + k]*anpC*d4*del_T*tc2)/5 - (Nc_n[i*N_SELE + k]*anpC*d1*del_T*tc2)/5; ktc364= (Nc_n[i*N_SELE + k]*anpC*d3*del_T*tc2)/5 - (4*Nc_n[i*N_SELE + k]*anpC*d6*del_T*tc2)/5 - (3*Nc_n[i*N_SELE + k]*anpC*d2*del_T*tc2)/(5*sel) + (3*Nc_n[i*N_SELE + k]*anpC*d5*del_T*tc2)/(5*sel); ktc365= (3*Nc_n[i*N_SELE + k]*anpC*d4*del_T*tc2)/(5*sel) - (3*Nc_n[i*N_SELE + k]*anpC*d1*del_T*tc2)/(5*sel); ktc366= 6*Mc_n[i*N_SELE + k]*anpC*d6*del_T*tc2 + (4*Nc_n[i*N_SELE + k]*anpC*d1*del_T*tc2)/5 (4*Nc_n[i*N_SELE + k]*anpC*d4*del_T*tc2)/5; ktc411= 0; ktc412= 0; ktc413= 0; ktc414= 0; ktc415= 0; ktc416= 0; ktc421= 0; ktc422= (1296*Nc_n[i*N_SELE + k]*anpC*d2*d5*del_T*tc2)/(35*(sel*sel*sel)) - (27*Nc_n[i*N_SELE + k]*anpC*(d3*d3)*del_T*tc2)/(35*sel) - (648*Nc_n[i*N_SELE + k]*anpC*(d5*d5)*del_T*tc2)/(35*(sel*sel*sel)) - (27*Nc_n[i*N_SELE + k]*anpC*(d6*d6)*del_T*tc2)/(35*sel) - (162*Nc_n[i*N_SELE + k]*anpC*d2*d3*del_T*tc2)/(35*(sel*sel)) - (648*Nc_n[i*N_SELE + k]*anpC*(d2*d2)*del_T*tc2)/(35*(sel*sel*sel)) - (162*Nc_n[i*N_SELE + k]*anpC*d2*d6*del_T*tc2)/(35*(sel*sel)) + (162*Nc_n[i*N_SELE + k]*anpC*d3*d5*del_T*tc2)/(35*(sel*sel)) + (162*Nc_n[i*N_SELE + k]*anpC*d5*d6*del_T*tc2)/(35*(sel*sel)); ktc423= (9*Nc_n[i*N_SELE + k]*anpC*(d3*d3)*del_T*tc2)/140 - (9*Nc_n[i*N_SELE + k]*anpC*(d6*d6)*del_T*tc2)/140 - (81*Nc_n[i*N_SELE + k]*anpC*(d2*d2)*del_T*tc2)/(35*(sel*sel)) - (81*Nc_n[i*N_SELE + k]*anpC*(d5*d5)*del_T*tc2)/(35*(sel*sel)) - (9*Nc_n[i*N_SELE + k]*anpC*d3*d6*del_T*tc2)/70 - (54*Nc_n[i*N_SELE + k]*anpC*d2*d3*del_T*tc2)/(35*sel) + (162*Nc_n[i*N_SELE + k]*anpC*d2*d5*del_T*tc2)/(35*(sel*sel)) + (54*Nc_n[i*N_SELE + k]*anpC*d3*d5*del_T*tc2)/(35*sel); ktc424= 0; Page: 54 CFT-Fire Code ktc425= (648*Nc_n[i*N_SELE + k]*anpC*(d2*d2)*del_T*tc2)/(35*(sel*sel*sel)) + (27*Nc_n[i*N_SELE + k]*anpC*(d3*d3)*del_T*tc2)/(35*sel) + (648*Nc_n[i*N_SELE + k]*anpC*(d5*d5)*del_T*tc2)/(35*(sel*sel*sel)) + (27*Nc_n[i*N_SELE + k]*anpC*(d6*d6)*del_T*tc2)/(35*sel) + (162*Nc_n[i*N_SELE + k]*anpC*d2*d3*del_T*tc2)/(35*(sel*sel)) - (1296*Nc_n[i*N_SELE + k]*anpC*d2*d5*del_T*tc2)/(35*(sel*sel*sel)) + (162*Nc_n[i*N_SELE + k]*anpC*d2*d6*del_T*tc2)/(35*(sel*sel)) - (162*Nc_n[i*N_SELE + k]*anpC*d3*d5*del_T*tc2)/(35*(sel*sel)) - (162*Nc_n[i*N_SELE + k]*anpC*d5*d6*del_T*tc2)/(35*(sel*sel)); ktc426= (9*Nc_n[i*N_SELE + k]*anpC*(d6*d6)*del_T*tc2)/140 - (9*Nc_n[i*N_SELE + k]*anpC*(d3*d3)*del_T*tc2)/140 - (81*Nc_n[i*N_SELE + k]*anpC*(d2*d2)*del_T*tc2)/(35*(sel*sel)) - (81*Nc_n[i*N_SELE + k]*anpC*(d5*d5)*del_T*tc2)/(35*(sel*sel)) - (9*Nc_n[i*N_SELE + k]*anpC*d3*d6*del_T*tc2)/70 + (162*Nc_n[i*N_SELE + k]*anpC*d2*d5*del_T*tc2)/(35*(sel*sel)) (54*Nc_n[i*N_SELE + k]*anpC*d2*d6*del_T*tc2)/(35*sel) + (54*Nc_n[i*N_SELE + k]*anpC*d5*d6*del_T*tc2)/(35*sel); ktc431= 0; ktc432= (9*Nc_n[i*N_SELE + k]*anpC*(d3*d3)*del_T*tc2)/140 - (9*Nc_n[i*N_SELE + k]*anpC*(d6*d6)*del_T*tc2)/140 - (81*Nc_n[i*N_SELE + k]*anpC*(d2*d2)*del_T*tc2)/(35*(sel*sel)) - (81*Nc_n[i*N_SELE + k]*anpC*(d5*d5)*del_T*tc2)/(35*(sel*sel)) - (9*Nc_n[i*N_SELE + k]*anpC*d3*d6*del_T*tc2)/70 - (54*Nc_n[i*N_SELE + k]*anpC*d2*d3*del_T*tc2)/(35*sel) + (162*Nc_n[i*N_SELE + k]*anpC*d2*d5*del_T*tc2)/(35*(sel*sel)) + (54*Nc_n[i*N_SELE + k]*anpC*d3*d5*del_T*tc2)/(35*sel); ktc433= (9*Nc_n[i*N_SELE + k]*anpC*d2*d3*del_T*tc2)/70 - (27*Nc_n[i*N_SELE + k]*anpC*(d5*d5)*del_T*tc2)/(35*sel) - (27*Nc_n[i*N_SELE + k]*anpC*(d2*d2)*del_T*tc2)/(35*sel) - (9*Nc_n[i*N_SELE + k]*anpC*d2*d6*del_T*tc2)/70 - (9*Nc_n[i*N_SELE + k]*anpC*d3*d5*del_T*tc2)/70 + (9*Nc_n[i*N_SELE + k]*anpC*d5*d6*del_T*tc2)/70 (18*sel*Nc_n[i*N_SELE + k]*anpC*(d3*d3)*del_T*tc2)/35 - (3*sel*Nc_n[i*N_SELE + k]*anpC*(d6*d6)*del_T*tc2)/70 + (54*Nc_n[i*N_SELE + k]*anpC*d2*d5*del_T*tc2)/(35*sel) + (9*sel*Nc_n[i*N_SELE + k]*anpC*d3*d6*del_T*tc2)/70; ktc434= 0; ktc435= (9*Nc_n[i*N_SELE + k]*anpC*(d6*d6)*del_T*tc2)/140 - (9*Nc_n[i*N_SELE + k]*anpC*(d3*d3)*del_T*tc2)/140 + (81*Nc_n[i*N_SELE + k]*anpC*(d2*d2)*del_T*tc2)/(35*(sel*sel)) + (81*Nc_n[i*N_SELE + k]*anpC*(d5*d5)*del_T*tc2)/(35*(sel*sel)) + (9*Nc_n[i*N_SELE + k]*anpC*d3*d6*del_T*tc2)/70 + (54*Nc_n[i*N_SELE + k]*anpC*d2*d3*del_T*tc2)/(35*sel) (162*Nc_n[i*N_SELE + k]*anpC*d2*d5*del_T*tc2)/(35*(sel*sel)) - (54*Nc_n[i*N_SELE + k]*anpC*d3*d5*del_T*tc2)/(35*sel); ktc436= (9*Nc_n[i*N_SELE + k]*anpC*d3*d5*del_T*tc2)/70 - (9*Nc_n[i*N_SELE + k]*anpC*d2*d6*del_T*tc2)/70 - (9*Nc_n[i*N_SELE + k]*anpC*d2*d3*del_T*tc2)/70 + (9*Nc_n[i*N_SELE + k]*anpC*d5*d6*del_T*tc2)/70 + (9*sel*Nc_n[i*N_SELE + k]*anpC*(d3*d3)*del_T*tc2)/140 + (9*sel*Nc_n[i*N_SELE + k]*anpC*(d6*d6)*del_T*tc2)/140 (3*sel*Nc_n[i*N_SELE + k]*anpC*d3*d6*del_T*tc2)/35; ktc441= 0; ktc442= 0; ktc443= 0; ktc444= 0; ktc445= 0; ktc446= 0; ktc451= 0; ktc452= (648*Nc_n[i*N_SELE + k]*anpC*(d2*d2)*del_T*tc2)/(35*(sel*sel*sel)) + (27*Nc_n[i*N_SELE + k]*anpC*(d3*d3)*del_T*tc2)/(35*sel) + (648*Nc_n[i*N_SELE + k]*anpC*(d5*d5)*del_T*tc2)/(35*(sel*sel*sel)) + (27*Nc_n[i*N_SELE + k]*anpC*(d6*d6)*del_T*tc2)/(35*sel) + (162*Nc_n[i*N_SELE + k]*anpC*d2*d3*del_T*tc2)/(35*(sel*sel)) - (1296*Nc_n[i*N_SELE + k]*anpC*d2*d5*del_T*tc2)/(35*(sel*sel*sel)) + (162*Nc_n[i*N_SELE + k]*anpC*d2*d6*del_T*tc2)/(35*(sel*sel)) - (162*Nc_n[i*N_SELE + k]*anpC*d3*d5*del_T*tc2)/(35*(sel*sel)) - (162*Nc_n[i*N_SELE + k]*anpC*d5*d6*del_T*tc2)/(35*(sel*sel)); ktc453= (9*Nc_n[i*N_SELE + k]*anpC*(d6*d6)*del_T*tc2)/140 - (9*Nc_n[i*N_SELE + k]*anpC*(d3*d3)*del_T*tc2)/140 + (81*Nc_n[i*N_SELE + k]*anpC*(d2*d2)*del_T*tc2)/(35*(sel*sel)) + (81*Nc_n[i*N_SELE + k]*anpC*(d5*d5)*del_T*tc2)/(35*(sel*sel)) + (9*Nc_n[i*N_SELE + k]*anpC*d3*d6*del_T*tc2)/70 + (54*Nc_n[i*N_SELE + k]*anpC*d2*d3*del_T*tc2)/(35*sel) (162*Nc_n[i*N_SELE + k]*anpC*d2*d5*del_T*tc2)/(35*(sel*sel)) - (54*Nc_n[i*N_SELE + k]*anpC*d3*d5*del_T*tc2)/(35*sel); ktc454= 0; ktc455= (1296*Nc_n[i*N_SELE + k]*anpC*d2*d5*del_T*tc2)/(35*(sel*sel*sel)) - (27*Nc_n[i*N_SELE + k]*anpC*(d3*d3)*del_T*tc2)/(35*sel) - (648*Nc_n[i*N_SELE + k]*anpC*(d5*d5)*del_T*tc2)/(35*(sel*sel*sel)) - (27*Nc_n[i*N_SELE + k]*anpC*(d6*d6)*del_T*tc2)/(35*sel) - (162*Nc_n[i*N_SELE + k]*anpC*d2*d3*del_T*tc2)/(35*(sel*sel)) - (648*Nc_n[i*N_SELE + k]*anpC*(d2*d2)*del_T*tc2)/(35*(sel*sel*sel)) - (162*Nc_n[i*N_SELE + k]*anpC*d2*d6*del_T*tc2)/(35*(sel*sel)) + (162*Nc_n[i*N_SELE + k]*anpC*d3*d5*del_T*tc2)/(35*(sel*sel)) + (162*Nc_n[i*N_SELE + k]*anpC*d5*d6*del_T*tc2)/(35*(sel*sel)); ktc456= (9*Nc_n[i*N_SELE + k]*anpC*(d3*d3)*del_T*tc2)/140 - (9*Nc_n[i*N_SELE + k]*anpC*(d6*d6)*del_T*tc2)/140 + (81*Nc_n[i*N_SELE + k]*anpC*(d2*d2)*del_T*tc2)/(35*(sel*sel)) + (81*Nc_n[i*N_SELE + k]*anpC*(d5*d5)*del_T*tc2)/(35*(sel*sel)) + (9*Nc_n[i*N_SELE + Page: 55 CFT-Fire Code k]*anpC*d3*d6*del_T*tc2)/70 - (162*Nc_n[i*N_SELE + k]*anpC*d2*d5*del_T*tc2)/(35*(sel*sel)) + (54*Nc_n[i*N_SELE + k]*anpC*d2*d6*del_T*tc2)/(35*sel) - (54*Nc_n[i*N_SELE + k]*anpC*d5*d6*del_T*tc2)/(35*sel); ktc461= 0; ktc462= (9*Nc_n[i*N_SELE + k]*anpC*(d6*d6)*del_T*tc2)/140 - (9*Nc_n[i*N_SELE + k]*anpC*(d3*d3)*del_T*tc2)/140 - (81*Nc_n[i*N_SELE + k]*anpC*(d2*d2)*del_T*tc2)/(35*(sel*sel)) - (81*Nc_n[i*N_SELE + k]*anpC*(d5*d5)*del_T*tc2)/(35*(sel*sel)) - (9*Nc_n[i*N_SELE + k]*anpC*d3*d6*del_T*tc2)/70 + (162*Nc_n[i*N_SELE + k]*anpC*d2*d5*del_T*tc2)/(35*(sel*sel)) (54*Nc_n[i*N_SELE + k]*anpC*d2*d6*del_T*tc2)/(35*sel) + (54*Nc_n[i*N_SELE + k]*anpC*d5*d6*del_T*tc2)/(35*sel); ktc463= (9*Nc_n[i*N_SELE + k]*anpC*d3*d5*del_T*tc2)/70 - (9*Nc_n[i*N_SELE + k]*anpC*d2*d6*del_T*tc2)/70 - (9*Nc_n[i*N_SELE + k]*anpC*d2*d3*del_T*tc2)/70 + (9*Nc_n[i*N_SELE + k]*anpC*d5*d6*del_T*tc2)/70 + (9*sel*Nc_n[i*N_SELE + k]*anpC*(d3*d3)*del_T*tc2)/140 + (9*sel*Nc_n[i*N_SELE + k]*anpC*(d6*d6)*del_T*tc2)/140 (3*sel*Nc_n[i*N_SELE + k]*anpC*d3*d6*del_T*tc2)/35; ktc464= 0; ktc465= (9*Nc_n[i*N_SELE + k]*anpC*(d3*d3)*del_T*tc2)/140 - (9*Nc_n[i*N_SELE + k]*anpC*(d6*d6)*del_T*tc2)/140 + (81*Nc_n[i*N_SELE + k]*anpC*(d2*d2)*del_T*tc2)/(35*(sel*sel)) + (81*Nc_n[i*N_SELE + k]*anpC*(d5*d5)*del_T*tc2)/(35*(sel*sel)) + (9*Nc_n[i*N_SELE + k]*anpC*d3*d6*del_T*tc2)/70 - (162*Nc_n[i*N_SELE + k]*anpC*d2*d5*del_T*tc2)/(35*(sel*sel)) + (54*Nc_n[i*N_SELE + k]*anpC*d2*d6*del_T*tc2)/(35*sel) - (54*Nc_n[i*N_SELE + k]*anpC*d5*d6*del_T*tc2)/(35*sel); ktc466= (9*Nc_n[i*N_SELE + k]*anpC*d2*d6*del_T*tc2)/70 - (27*Nc_n[i*N_SELE + k]*anpC*(d5*d5)*del_T*tc2)/(35*sel) - (9*Nc_n[i*N_SELE + k]*anpC*d2*d3*del_T*tc2)/70 (27*Nc_n[i*N_SELE + k]*anpC*(d2*d2)*del_T*tc2)/(35*sel) + (9*Nc_n[i*N_SELE + k]*anpC*d3*d5*del_T*tc2)/70 - (9*Nc_n[i*N_SELE + k]*anpC*d5*d6*del_T*tc2)/70 (3*sel*Nc_n[i*N_SELE + k]*anpC*(d3*d3)*del_T*tc2)/70 - (18*sel*Nc_n[i*N_SELE + k]*anpC*(d6*d6)*del_T*tc2)/35 + (54*Nc_n[i*N_SELE + k]*anpC*d2*d5*del_T*tc2)/(35*sel) + (9*sel*Nc_n[i*N_SELE + k]*anpC*d3*d6*del_T*tc2)/70; // thiet lap ma tran cung cua phan tu sub_sm[0][0] = +Etb*A_e[i*N_SELE + k]/sel +(2.0*Nc_n[i*N_SELE + k]*tc1)/sel -(6.0*tc2*(Mc_n[i*N_SELE + k]*d3 - Mc_n[i*N_SELE + k]*d6 Nc_n[i*N_SELE + k]*d1 + Nc_n[i*N_SELE + k]*d4))/(sel*sel) -(Nc_n[i*N_SELE + k]*tc2*(2.0*d3*d3 - d3*d6 + 2.0*d6*d6))/(5.0*sel) (18.0*Nc_n[i*N_SELE + k]*tc2*(d2 - d5)*(d2 - d5))/(5.0*(sel*sel*sel)) (3.0*Nc_n[i*N_SELE + k]*tc2*(d3 + d6)*(d2 - d5))/(5.0*(sel*sel)) + heso1*(kts11) + heso2*(ktc211+ktc311+ktc411)+ heso3*(kps211+kps311+kps411)+heso4*(ktp11); sub_sm[0][1] = +Etb*A_e[i*N_SELE + k]*(6/(5*sel*sel)*(d5 - d2)- 1/(10*sel)*(d3 + d6)) +(36.0*Ic_ze[i*N_SELE + k]*ftb_c*tc2*(d3 + d6))/(sel*sel*sel) (Nc_n[i*N_SELE + k]*tc1*(d3 + d6))/(5.0*sel) - (12.0*Nc_n[i*N_SELE + k]*tc1*(d2 - d5))/(5.0*(sel*sel)) + (72.0*Ic_ze[i*N_SELE + k]*ftb_c*tc2*(d2 - d5))/(sel*sel*sel*sel) -(3.0*Nc_n[i*N_SELE + k]*tc2*(d1 - d4)*(12.0*d2 - 12.0*d5 + sel*d3 + sel*d6))/(5.0*(sel*sel*sel)) +tc2*((216.0*Nc_n[i*N_SELE + k]*(d2 - d5)*(d2 - d5)*(d2 d5))/(35.0*(sel*sel*sel*sel)) + (27.0*Nc_n[i*N_SELE + k]*(2*d2 2*d5)*(d3*d3 + d6*d6))/(70.0*(sel*sel)) + (81.0*Nc_n[i*N_SELE + k]*(d3 + d6)*(d2 - d5)*(d2 - d5))/(35.0*(sel*sel*sel)) - (3*Nc_n[i*N_SELE + k]*(d3 + d6)*(d3*d3 - 4*d3*d6 + d6*d6))/(140.0*sel)) + heso1*(kts12) + heso2*(ktc212+ktc312+ktc412)+ heso3*(kps212+kps312+kps412)+heso4*(ktp12); sub_sm[0][2] = -Etb*S_e[i*N_SELE + k]/sel +Etb*A_e[i*N_SELE + k]*( 1/(10*sel)*(d5 - d2) + 1.0/30*(d6 - 4*d3)) -(2.0*Mc_n[i*N_SELE + k]*tc1)/sel +(Nc_n[i*N_SELE + k]*d6*tc1)/15.0 - (4.0*Nc_n[i*N_SELE + k]*d3*tc1)/15.0 + (6.0*tc2*(Mc_n[i*N_SELE + k]*d4 - Mc_n[i*N_SELE + k]*d1 + 4.0*Ic_ze[i*N_SELE + k]*d3* ftb_c + 2.0*Ic_ze[i*N_SELE + k]*d6* ftb_c))/(sel*sel) - (Nc_n[i*N_SELE + k]*tc1*(d2 - d5))/(5.0*sel) + (36.0*Ic_ze[i*N_SELE + k]* ftb_c*tc2*(d2 - d5))/(sel*sel*sel) +(tc2*(15.0*Mc_n[i*N_SELE + k]*d3*d3 - 4.0*Nc_n[i*N_SELE + k]*d1*d3 + Nc_n[i*N_SELE + k]*d1*d6 + 4.0*Nc_n[i*N_SELE + k]*d3*d4 - Nc_n[i*N_SELE + k]*d4*d6))/(5.0*sel) - (3.0*Nc_n[i*N_SELE + k]*tc2*(d1 - d4)*(d2 d5))/(5.0*(sel*sel)) +(6.0*Nc_n[i*N_SELE + k]*d3*d3*d3*tc2)/35.0 - (3.0*Nc_n[i*N_SELE + k]*d6*d6*d6*tc2)/140.0 + (3.0*Nc_n[i*N_SELE + k]*d3*d6*d6*tc2)/70.0 (9.0*Nc_n[i*N_SELE + k]*d3*d3*d6*tc2)/140.0 + (27.0*Nc_n[i*N_SELE + k]*tc2*(d2 - d5)*(d2 - d5)*(d2 - d5))/(35.0*(sel*sel*sel)) + (27.0*Nc_n[i*N_SELE + k]*d3*tc2*(d2 - d5)*(d2 - d5))/(35.0*(sel*sel)) + (9.0*Nc_n[i*N_SELE + k]*tc2*(d2 - d5)*(2*d3*d6 - d3*d3 + Page: 56 CFT-Fire Code d6*d6))/(140.0*sel) + heso1*(kts13) + heso2*(ktc213+ktc313+ktc413)+ heso3*(kps113+kps213+kps313+kps413)+heso4*(ktp13); sub_sm[0][3] = -1*A_e[i*N_SELE + k]*Etb/sel -(2.0*Nc_n[i*N_SELE + k]*tc1)/sel +(6.0*tc2*(Mc_n[i*N_SELE + k]*d3 - Mc_n[i*N_SELE + k]*d6 Nc_n[i*N_SELE + k]*d1 + Nc_n[i*N_SELE + k]*d4))/(sel*sel) +(Nc_n[i*N_SELE + k]*tc2*(2.0*d3*d3 - d3*d6 + 2.0*d6*d6))/(5.0*sel) + (18.0*Nc_n[i*N_SELE + k]*tc2*(d2 - d5)*(d2 - d5))/(5.0*(sel*sel*sel)) + (3.0*Nc_n[i*N_SELE + k]*tc2*(d3 + d6)*(d2 - d5))/(5.0*(sel*sel)) + heso1*(kts14) + heso2*(ktc214+ktc314+ktc414)+ heso3*(kps214+kps314+kps414)+heso4*(ktp14); sub_sm[0][4] = +Etb*A_e[i*N_SELE + k]*(-6/(5*sel*sel)*(d5 - d2)+ 1/(10.0*sel)*(d3 + d6)) +(12.0* Nc_n[i*N_SELE + k]*tc1*(d2 - d5))/(5.0*(sel*sel)) + ( Nc_n[i*N_SELE + k]*tc1*(d3 + d6))/(5.0*sel) - (36.0*Ic_ze[i*N_SELE + k]*ftb_c*tc2*(d3 + d6))/(sel*sel*sel) - (72.0*Ic_ze[i*N_SELE + k]*ftb_c*tc2*(d2 - d5))/(sel*sel*sel*sel) +(3.0*Nc_n[i*N_SELE + k]*tc2*(d1 - d4)*(12.0*d2 - 12.0*d5 + sel*d3 + sel*d6))/(5.0*(sel*sel*sel)) -tc2*((216.0* Nc_n[i*N_SELE + k]*(d2 - d5)*(d2 - d5)*(d2 d5))/(35.0*(sel*sel*sel*sel)) + (27.0* Nc_n[i*N_SELE + k]*(2*d2 2*d5)*(d3*d3 + d6*d6))/(70.0*(sel*sel)) + (81.0* Nc_n[i*N_SELE + k]*(d3 + d6)*(d2 - d5)*(d2 - d5))/(35.0*(sel*sel*sel)) - (3.0* Nc_n[i*N_SELE + k]*(d3 + d6)*(d3*d3 - 4*d3*d6 + d6*d6))/(140.0*sel)) + heso1*(kts15) + heso2*(ktc215+ktc315+ktc415)+ heso3*(kps215+kps315+kps415)+heso4*(ktp15); sub_sm[0][5] = +Etb*S_e[i*N_SELE + k]/sel +Etb*A_e[i*N_SELE + k]*(1/(10*sel)*(d5 - d2) + 1.0/30*(d3 - 4*d6)) +(2.0*Mc_n[i*N_SELE + k]*tc1)/sel +( Nc_n[i*N_SELE + k]*d3*tc1)/15.0 - (4.0*Nc_n[i*N_SELE + k]*d6*tc1)/15.0 + (6.0*tc2*(Mc_n[i*N_SELE + k]*d1 - Mc_n[i*N_SELE + k]*d4 + 2.0*Ic_ze[i*N_SELE + k]*d3*ftb_c + 4.0*Ic_ze[i*N_SELE + k]*d6*ftb_c))/(sel*sel) - ( Nc_n[i*N_SELE + k]*tc1*(d2 - d5))/(5.0*sel) + (36.0*Ic_ze[i*N_SELE + k]*ftb_c*tc2*(d2 - d5))/(sel*sel*sel) -(tc2*(15.0*Mc_n[i*N_SELE + k]*d6*d6 - Nc_n[i*N_SELE + k]*d1*d3 + 4.0*Nc_n[i*N_SELE + k]*d1*d6 + Nc_n[i*N_SELE + k]*d3*d4 4.0*Nc_n[i*N_SELE + k]*d4*d6))/(5.0*sel) - (3.0* Nc_n[i*N_SELE + k]*tc2*(d1 - d4)*(d2 - d5))/(5.0*(sel*sel)) +(6.0*Nc_n[i*N_SELE + k]*d6*d6*d6*tc2)/35.0 - (3.0*Nc_n[i*N_SELE + k]*d3*d3*d3*tc2)/140.0 - (9.0*Nc_n[i*N_SELE + k]*d3*d6*d6*tc2)/140.0 + (3.0*Nc_n[i*N_SELE + k]*d3*d3*d6*tc2)/70.0 + (27.0*Nc_n[i*N_SELE + k]*tc2*(d2 - d5)*(d2 - d5)*(d2 - d5))/(35.0*(sel*sel*sel)) + (27.0* Nc_n[i*N_SELE + k]*d6*tc2*(d2 - d5)*(d2 - d5))/(35.0*(sel*sel)) + (9.0*Nc_n[i*N_SELE + k]*tc2*(d2 - d5)*(d3*d3 + 2*d3*d6 d6*d6))/(140.0*sel) + heso1*(kts16) + heso2*(ktc216+ktc316+ktc416)+ heso3*(kps216+kps316+kps416)+heso4*(ktp16); sub_sm[1][1] = 12*Etb*I_ze[i*N_SELE + k]/(sel*sel*sel) + Etb*A_e[i*N_SELE + k]*(6/(5*sel*sel)*(d4 - d1) + 1.0/140*(18/sel*(d3*d3 + d6*d6) + 432/(sel*sel*sel)*(d5 - d2)*(d5 - d2) - 108/(sel*sel)*(d5 - d2)*(d3 + d6))) + 6/(5*sel)*N_p[i*N_SELE + k] +(24.0*Ic_ze[i*N_SELE + k]*ftb_c*tc1)/(sel*sel*sel) +(6.0*Nc_p[i*N_SELE + k])/(5.0*sel) -(12.0*Nc_n[i*N_SELE + k]*tc1*(d1 - d4))/(5.0*(sel*sel)) (72.0*ftb_c*tc2*(Cc_n[i*N_SELE + k]*d3 - Cc_n[i*N_SELE + k]*d6 Ic_ze[i*N_SELE + k]*d1 + Ic_ze[i*N_SELE + k]*d4))/(sel*sel*sel*sel) +(9.0*Nc_n[i*N_SELE + k]*tc1*(d3*d3 + d6*d6))/(35.0*sel) (126.0*Nc_n[i*N_SELE + k]*tc2*d1*d1 - 252.0*Nc_n[i*N_SELE + k]*tc2*d1*d4 216.0*Nc_n[i*N_SELE + k]*tc1*d2*d2 + 432*Nc_n[i*N_SELE + k]*tc1*d2*d5 + 216.0*Ic_ze[i*N_SELE + k]*ftb_c*tc2*d3*d3 + 180.0*Ic_ze[i*N_SELE + k]*ftb_c*tc2*d3*d6 + 126.0*Nc_n[i*N_SELE + k]*tc2*d4*d4 - 216.0*Nc_n[i*N_SELE + k]*tc1*d5*d5 + 216.0*Ic_ze[i*N_SELE + k]*ftb_c*tc2*d6*d6)/(35.0*(sel*sel*sel)) - (3888.0*Ic_ze[i*N_SELE + k]*ftb_c*tc2*(d2 - d5)*(d2 - d5))/(35.0*(sel*sel*sel*sel*sel)) + (54.0*Nc_n[i*N_SELE + k]*tc1*(d3 + d6)*(d2 - d5))/(35.0*(sel*sel)) (324.0*Ic_ze[i*N_SELE + k]*ftb_c*tc2*(d3 + d6)*(d2 d5))/(7.0*(sel*sel*sel*sel)) +(27.0*Nc_n[i*N_SELE + k]*tc2*(d1 - d4)*((sel*sel)*d3*d3 + (sel*sel)*d6*d6 + 6.0*sel*d2*d3 + 6.0*sel*d2*d6 - 6.0*sel*d3*d5 - 6.0*sel*d5*d6 + 24.0*d2*d2 48.0*d2*d5 + 24.0*d5*d5))/(35.0*(sel*sel*sel*sel)) -tc2*((14580.0*Nc_n[i*N_SELE + k]*(d2 - d5)*(d2 - d5)*(d2 - d5)*(d2 Page: 57 CFT-Fire Code d5))/(1001.0*(sel*sel*sel*sel*sel)) + (3.0*Nc_n[i*N_SELE + k]*(29.0*d3*d3*d3*d3 - 14.0*d3*d3*d3*d6 + 57.0*d3*d3*d6*d6 - 14.0*d3*d6*d6*d6 + 29.0*d6*d6*d6*d6))/(2002.0*sel) + (324.0*Nc_n[i*N_SELE + k]*(d2 - d5)*(d2 d5)*(9.0*d3*d3 + 5.0*d3*d6 + 9.0*d6*d6))/(1001.0*(sel*sel*sel)) + (8100.0*Nc_n[i*N_SELE + k]*(d3 + d6)*(d2 - d5)*(d2 - d5)*(d2 d5))/(1001.0*(sel*sel*sel*sel)) + (27.0*Nc_n[i*N_SELE + k]*(d3 + d6)*(2.0*d2 2.0*d5)*(10.0*d3*d3 + 7.0*d3*d6 + 10.0*d6*d6))/(2002.0*(sel*sel))) + heso1*(kts22) + heso2*(ktc222+ktc322+ktc422)+ heso3*(kps222+kps322+kps422)+heso4*(ktp22); sub_sm[1][2] = 6*Etb*I_ze[i*N_SELE + k]/(sel*sel) + Etb*A_e[i*N_SELE + k]*( 1/(10*sel)*(d4 - d1) + 1.0/280*(3*(d6*d6 - d3*d3) + 6*d3*d6 - 72/sel*d3*(d5 - d2) + 108/(sel*sel)*(d5 - d2)*(d5 - d2))) + 1.0/10*N_p[i*N_SELE + k] +(12.0*Ic_ze[i*N_SELE + k]*ftb_c*tc1)/(sel*sel) +Nc_p[i*N_SELE + k]/10.0 -(Nc_n[i*N_SELE + k]*tc1*(d1 - d4))/(5.0*sel) - (36.0*ftb_c*tc2*(2.0*Cc_n[i*N_SELE + k]*d3 Ic_ze[i*N_SELE + k]*d1 + Ic_ze[i*N_SELE + k]*d4))/(sel*sel*sel) - (72.0*Cc_n[i*N_SELE + k]*ftb_c*tc2*(d2 - d5))/(sel*sel*sel*sel) +(3.0*Nc_n[i*N_SELE + k]*d6*d6*tc1)/140.0 - (3.0*Nc_n[i*N_SELE + k]*d3*d3*tc1)/140.0 (42.0*Nc_n[i*N_SELE + k]*tc2*d1*d1 - 84.0*Nc_n[i*N_SELE + k]*tc2*d1*d4 - 108.0*Nc_n[i*N_SELE + k]*tc1*d2*d2 + 216.0*Nc_n[i*N_SELE + k]*tc1*d2*d5 + 1440.0*Ic_ze[i*N_SELE + k]*ftb_c*tc2*d3*d3 + 192.0*Ic_ze[i*N_SELE + k]*ftb_c*tc2*d3*d6 + 42.0*Nc_n[i*N_SELE + k]*tc2*d4*d4 108.0*Nc_n[i*N_SELE + k]*tc1*d5*d5 + 96.0*Ic_ze[i*N_SELE + k]*ftb_c*tc2*d6*d6)/(140.0*(sel*sel)) + (3.0*Nc_n[i*N_SELE + k]*d3*d6*tc1)/70.0 (162.0*Ic_ze[i*N_SELE + k]*ftb_c*tc2*(d2 - d5)*(d2 - d5))/(7.0*(sel*sel*sel*sel)) + (18.0*Nc_n[i*N_SELE + k]*d3*tc1*(d2 - d5))/(35.0*sel) - (36.0*Ic_ze[i*N_SELE + k]*ftb_c*tc2*(d2 - d5)*(12.0*d3 + 5.0*d6))/(35.0*(sel*sel*sel)) +(9.0*Nc_n[i*N_SELE + k]*tc2*(d1 - d4)*(2.0*(sel*sel)*d3*d6 - (sel*sel)*d3*d3 + (sel*sel)*d6*d6 + 24.0*sel*d2*d3 - 24.0*sel*d3*d5 + 36.0*d2*d2 - 72.0*d2*d5 + 36.0*d5*d5))/(140.0*(sel*sel*sel)) -tc2*((57.0*Nc_n[i*N_SELE + k]*d6*d6*d6*d6)/8008.0 - (255.0*Nc_n[i*N_SELE + k]*d3*d3*d3*d3)/8008.0 - (9.0*Nc_n[i*N_SELE + k]*d3*d3*d6*d6)/8008.0 + (2025.0*Nc_n[i*N_SELE + k]*(d2 - d5)*(d2 - d5)*(d2 - d5)*(d2 - d5))/(1001.0*(sel*sel*sel*sel)) - (3.0*Nc_n[i*N_SELE + k]*d3*d6*d6*d6)/4004.0 + (57.0*Nc_n[i*N_SELE + k]*d3*d3*d3*d6)/2002.0 + (27.0*Nc_n[i*N_SELE + k]*(d2 - d5)*(d2 - d5)*(10.0*d3*d3 + 7.0*d3*d6 + 10.0*d6*d6))/(2002.0*(sel*sel)) + (108.0*Nc_n[i*N_SELE + k]*(d2 - d5)*(d2 - d5)*(d2 - d5)*(18.0*d3 + 5.0*d6))/(1001.0*(sel*sel*sel)) + (3.0*Nc_n[i*N_SELE + k]*(2.0*d2 - 2.0*d5)*(116.0*d3*d3*d3 42.0*d3*d3*d6 + 114.0*d3*d6*d6 - 14.0*d6*d6*d6))/(4004.0*sel) + (27.0*Nc_n[i*N_SELE + k]*(d3 + d6)*(d2 - d5)*(d2 - d5)*(20.0*d3 + 7.0*d6))/(2002.0*(sel*sel))) - heso1*(kts23) + heso2*(ktc223+ktc323+ktc423)+ heso3*(kps223+kps323+kps423)+heso4*(ktp23); sub_sm[1][3] =+Etb*A_e[i*N_SELE + k]*( -6/(5*sel*sel)*(d5 - d2) + 1/(10*sel)*(d3 + d6)) +(12.0*Nc_n[i*N_SELE + k]*tc1*(d2 - d5))/(5.0*(sel*sel)) + (Nc_n[i*N_SELE + k]*tc1*(d3 + d6))/(5.0*sel) - (36.0*Ic_ze[i*N_SELE + k]*ftb_c*tc2*(d3 + d6))/(sel*sel*sel) - (72.0*Ic_ze[i*N_SELE + k]*ftb_c*tc2*(d2 - d5))/(sel*sel*sel*sel) +(3.0*Nc_n[i*N_SELE + k]*tc2*(d1 - d4)*(12.0*d2 - 12.0*d5 + sel*d3 + sel*d6))/(5.0*(sel*sel*sel)) -tc2*((216.0*Nc_n[i*N_SELE + k]*(d2 - d5)*(d2 - d5)*(d2 d5))/(35.0*(sel*sel*sel*sel)) + (27.0*Nc_n[i*N_SELE + k]*(2.0*d2 2.0*d5)*(d3*d3 + d6*d6))/(70.0*(sel*sel)) + (81.0*Nc_n[i*N_SELE + k]*(d3 + d6)*(d2 - d5)*(d2 - d5))/(35.0*(sel*sel*sel)) - (3.0*Nc_n[i*N_SELE + k]*(d3 + d6)*(d3*d3 - 4.0*d3*d6 + d6*d6))/(140.0*sel)) + heso1*(kts24) + heso2*(ktc224+ktc324+ktc424)+ heso3*(kps224+kps324+kps424)+heso4*(ktp24); sub_sm[1][4] = -12*Etb*I_ze[i*N_SELE + k]/(sel*sel*sel) + Etb*A_e[i*N_SELE + k]*(-6/(5*sel*sel)*(d4 - d1) + 1.0/140*(-18/sel*(d3*d3 + d6*d6) - 432/(sel*sel*sel)*(d5 - d2)*(d5 - d2) + 108/(sel*sel)*(d5 - d2)*(d3 + d6))) - 6/(5*sel)*N_p[i*N_SELE + k] -(24.0*Ic_ze[i*N_SELE + k]*ftb_c*tc1)/(sel*sel*sel) -(6.0*Nc_p[i*N_SELE + k])/(5.0*sel) +(12.0*Nc_n[i*N_SELE + k]*tc1*(d1 - d4))/(5.0*(sel*sel)) + (72.0*ftb_c*tc2*(Cc_n[i*N_SELE + k]*d3 - Cc_n[i*N_SELE + k]*d6 Ic_ze[i*N_SELE + k]*d1 + Ic_ze[i*N_SELE + k]*d4))/(sel*sel*sel*sel) +(126.0*Nc_n[i*N_SELE + k]*tc2*d1*d1 - 252.0*Nc_n[i*N_SELE + k]*tc2*d1*d4 216.0*Nc_n[i*N_SELE + k]*tc1*d2*d2 + 432.0*Nc_n[i*N_SELE + k]*tc1*d2*d5 + 216.0*Ic_ze[i*N_SELE + k]*ftb_c*tc2*d3*d3 + 180.0*Ic_ze[i*N_SELE + k]*ftb_c*tc2*d3*d6 + 126.0*Nc_n[i*N_SELE + k]*tc2*d4*d4 - 216.0*Nc_n[i*N_SELE + k]*tc1*d5*d5 + 216.0*Ic_ze[i*N_SELE + k]*ftb_c*tc2*d6*d6)/(35.0*(sel*sel*sel)) - (9.0*Nc_n[i*N_SELE + k]*tc1*(d3*d3 + d6*d6))/(35.0*sel) + (3888.0*Ic_ze[i*N_SELE + k]*ftb_c*tc2*(d2 - d5)*(d2 Page: 58 CFT-Fire Code d5))/(35.0*(sel*sel*sel*sel*sel)) - (54.0*Nc_n[i*N_SELE + k]*tc1*(d3 + d6)*(d2 - d5))/(35.0*(sel*sel)) + (324.0*Ic_ze[i*N_SELE + k]*ftb_c*tc2*(d3 + d6)*(d2 d5))/(7.0*(sel*sel*sel*sel)) -(27.0*Nc_n[i*N_SELE + k]*tc2*(d1 - d4)*((sel*sel)*d3*d3 + (sel*sel)*d6*d6 + 6.0*sel*d2*d3 + 6.0*sel*d2*d6 - 6.0*sel*d3*d5 - 6.0*sel*d5*d6 + 24.0*d2*d2 48.0*d2*d5 + 24.0*d5*d5))/(35.0*(sel*sel*sel*sel)) +tc2*((14580.0*Nc_n[i*N_SELE + k]*(d2 - d5)*(d2 - d5)*(d2 - d5)*(d2 d5))/(1001.0*(sel*sel*sel*sel*sel)) + (3.0*Nc_n[i*N_SELE + k]*(29.0*d3*d3*d3*d3 - 14.0*d3*d3*d3*d6 + 57.0*d3*d3*d6*d6 - 14.0*d3*d6*d6*d6 + 29.0*d6*d6*d6*d6))/(2002.0*sel) + (324.0*Nc_n[i*N_SELE + k]*(d2 - d5)*(d2 d5)*(9.0*d3*d3 + 5.0*d3*d6 + 9.0*d6*d6))/(1001.0*(sel*sel*sel)) + (8100.0*Nc_n[i*N_SELE + k]*(d3 + d6)*(d2 - d5)*(d2 - d5)*(d2 d5))/(1001.0*(sel*sel*sel*sel)) + (27.0*Nc_n[i*N_SELE + k]*(d3 + d6)*(2.0*d2 2.0*d5)*(10.0*d3*d3 + 7.0*d3*d6 + 10.0*d6*d6))/(2002.0*(sel*sel))) + heso1*(kts25) + heso2*(ktc225+ktc325+ktc425)+ heso3*(kps225+kps325+kps425)+heso4*(ktp25); sub_sm[1][5] = 6*Etb*I_ze[i*N_SELE + k]/(sel*sel) + Etb*A_e[i*N_SELE + k]*(1/(10*sel)*(d4 - d1) + 1.0/280*(3*(d3*d3 - d6*d6) + 6*d3*d6 - 72/sel*d6*(d5 - d2) + 108/(sel*sel)*(d5 - d2)*(d5 - d2))) + 1.0/10*N_p[i*N_SELE + k] +(12.0* Ic_ze[i*N_SELE + k]*ftb_c*tc1)/(sel*sel) +Nc_p[i*N_SELE + k]/10.0 +(36.0*ftb_c*tc2*(2.0*Cc_n[i*N_SELE + k]*d6 + Ic_ze[i*N_SELE + k]*d1 Ic_ze[i*N_SELE + k]*d4))/(sel*sel*sel) - (Nc_n[i*N_SELE + k]*tc1*(d1 d4))/(5.0*sel) + (72.0*Cc_n[i*N_SELE + k]*ftb_c*tc2*(d2 d5))/(sel*sel*sel*sel) +(3.0*Nc_n[i*N_SELE + k]*d3*d3*tc1)/140.0 - (42.0*Nc_n[i*N_SELE + k]*tc2*d1*d1 - 84.0*Nc_n[i*N_SELE + k]*tc2*d1*d4 - 108.0*Nc_n[i*N_SELE + k]*tc1*d2*d2 + 216.0*Nc_n[i*N_SELE + k]*tc1*d2*d5 + 96.0*Ic_ze[i*N_SELE + k]*ftb_c*tc2*d3*d3 + 192.0* Ic_ze[i*N_SELE + k]*ftb_c*tc2*d3*d6 + 42.0*Nc_n[i*N_SELE + k]*tc2*d4*d4 - 108.0*Nc_n[i*N_SELE + k]*tc1*d5*d5 + 1440.0* Ic_ze[i*N_SELE + k]*ftb_c*tc2*d6*d6)/(140.0*(sel*sel)) - (3.0*Nc_n[i*N_SELE + k]*d6*d6*tc1)/140.0 + (3.0*Nc_n[i*N_SELE + k]*d3*d6*tc1)/70.0 - (162.0* Ic_ze[i*N_SELE + k]*ftb_c*tc2*(d2 - d5)*(d2 - d5))/(7.0*(sel*sel*sel*sel)) + (18.0*Nc_n[i*N_SELE + k]*d6*tc1*(d2 - d5))/(35.0*sel) - (36.0* Ic_ze[i*N_SELE + k]*ftb_c*tc2*(d2 - d5)*(5.0*d3 + 12.0*d6))/(35.0*(sel*sel*sel)) +(9.0*Nc_n[i*N_SELE + k]*tc2*(d1 - d4)*((sel*sel)*d3*d3 + 2.0*(sel*sel)*d3*d6 - (sel*sel)*d6*d6 + 24.0*sel*d2*d6 - 24.0*sel*d5*d6 + 36.0*d2*d2 - 72.0*d2*d5 + 36.0*d5*d5))/(140.0*(sel*sel*sel)) -tc2*((57.0*Nc_n[i*N_SELE + k]*d3*d3*d3*d3)/8008.0 - (255.0*Nc_n[i*N_SELE + k]*d6*d6*d6*d6)/8008.0 - (9.0*Nc_n[i*N_SELE + k]*d3*d3*d6*d6)/8008.0 + (2025.0*Nc_n[i*N_SELE + k]*(d2 - d5)*(d2 - d5)*(d2 - d5)*(d2 d5))/(1001.0*(sel*sel*sel*sel)) + (57.0*Nc_n[i*N_SELE + k]*d3*d6*d6*d6)/2002.0 - (3.0*Nc_n[i*N_SELE + k]*d3*d3*d3*d6)/4004.0 + (27.0*Nc_n[i*N_SELE + k]*(d2 d5)*(d2 - d5)*(10.0*d3*d3 + 7.0*d3*d6 + 10.0*d6*d6))/(2002.0*(sel*sel)) + (108.0*Nc_n[i*N_SELE + k]*(d2 - d5)*(d2 - d5)*(d2 - d5)*(5.0*d3 + 18.0*d6))/(1001.0*(sel*sel*sel)) - (3.0*Nc_n[i*N_SELE + k]*(2.0*d2 2.0*d5)*(14.0*d3*d3*d3 - 114.0*d3*d3*d6 + 42.0*d3*d6*d6 116.0*d6*d6*d6))/(4004.0*sel) + (27.0*Nc_n[i*N_SELE + k]*(d3 + d6)*(d2 d5)*(d2 - d5)*(7.0*d3 + 20.0*d6))/(2002.0*(sel*sel))) + heso1*(kts26) + heso2*(ktc226+ktc326+ktc426)+ heso3*(kps226+kps326+kps426)+heso4*(ktp26); sub_sm[2][2] = 4*Etb*I_ze[i*N_SELE + k]/sel +Etb*S_e[i*N_SELE + k]*d3 + Etb*A_e[i*N_SELE + k]*(2.0/15*(d4 - d1) + 1.0/140*(12*sel*d3*d3 + sel*d6*d6 - 3*sel*d3*d6 + 18/sel*(d5 - d2)*(d5 - d2) + 3*(d5 - d2)*(d3 - d6))) + 2*sel/15*N_p[i*N_SELE + k] +(8.0*Ic_ze[i*N_SELE + k]*ftb_c*tc1)/sel +(2.0*sel*Nc_p[i*N_SELE + k])/15.0 +2.0*Mc_n[i*N_SELE + k]*d3*tc1 - (4.0*Nc_n[i*N_SELE + k]*d1*tc1)/15.0 + (4.0*Nc_n[i*N_SELE + k]*d4*tc1)/15.0 - (12.0*ftb_c*tc2*(5.0*Cc_n[i*N_SELE + k]*d3 + Cc_n[i*N_SELE + k]*d6 - 2.0*Ic_ze[i*N_SELE + k]*d1 + 2.0*Ic_ze[i*N_SELE + k]*d4))/(sel*sel) - (72.0*Cc_n[i*N_SELE + k]*ftb_c*tc2*(d2 - d5))/(sel*sel*sel) +sel*((6.0*Nc_n[i*N_SELE + k]*tc1*d3*d3)/35.0 - (3.0*Nc_n[i*N_SELE + k]*tc1*d3*d6)/70.0 + (Nc_n[i*N_SELE + k]*tc1*d6*d6)/70.0) (28.0*Nc_n[i*N_SELE + k]*tc2*d1*d1 - 420.0*Mc_n[i*N_SELE + k]*tc2*d1*d3 56.0*Nc_n[i*N_SELE + k]*tc2*d1*d4 - 18.0*Nc_n[i*N_SELE + k]*tc1*d2*d2 + 36.0*Nc_n[i*N_SELE + k]*tc1*d2*d5 + 1536.0*Ic_ze[i*N_SELE + k]*ftb_c*tc2*d3*d3 + 420.0*Mc_n[i*N_SELE + k]*tc2*d3*d4 + 384.0*Ic_ze[i*N_SELE + k]*ftb_c*tc2*d3*d6 + 28.0*Nc_n[i*N_SELE + k]*tc2*d4*d4 - 18.0*Nc_n[i*N_SELE + k]*tc1*d5*d5 + 24.0*Ic_ze[i*N_SELE + k]*ftb_c*tc2*d6*d6)/(70.0*sel) Page: 59 CFT-Fire Code (3.0*Nc_n[i*N_SELE + k]*d2*d3*tc1)/70.0 + (3.0*Nc_n[i*N_SELE + k]*d2*d6*tc1)/70.0 + (3.0*Nc_n[i*N_SELE + k]*d3*d5*tc1)/70.0 (3.0*Nc_n[i*N_SELE + k]*d5*d6*tc1)/70.0 - (216.0*Ic_ze[i*N_SELE + k]*ftb_c*tc2*(d2 - d5)*(d2 - d5))/(35.0*(sel*sel*sel)) - (48.0*Ic_ze[i*N_SELE + k]*ftb_c*tc2*(d2 - d5)*(15.0*d3 + d6))/(35.0*(sel*sel)) +(18.0*Nc_n[i*N_SELE + k]*d1*d3*d3*tc2)/35.0 - 3.0*Mc_n[i*N_SELE + k]*d3*d3*d3*tc2 + (3.0*Nc_n[i*N_SELE + k]*d1*d6*d6*tc2)/70.0 (18.0*Nc_n[i*N_SELE + k]*d3*d3*d4*tc2)/35.0 - (3.0*Nc_n[i*N_SELE + k]*d4*d6*d6*tc2)/70.0 + (27.0*Nc_n[i*N_SELE + k]*tc2*(d1 - d4)*(d2 - d5)*(d2 d5))/(35.0*(sel*sel)) - (9.0*Nc_n[i*N_SELE + k]*d1*d3*d6*tc2)/70.0 + (9.0*Nc_n[i*N_SELE + k]*d3*d4*d6*tc2)/70.0 - (9.0*Nc_n[i*N_SELE + k]*tc2*(d1 d4)*(d2 - d5)*(d3 - d6))/(70.0*sel) -tc2*(sel*((285.0*Nc_n[i*N_SELE + k]*d3*d3*d3*d3)/2002.0 - (95.0*Nc_n[i*N_SELE + k]*d3*d3*d3*d6)/2002.0 + (3.0*Nc_n[i*N_SELE + k]*d3*d3*d6*d6)/182.0 (3.0*Nc_n[i*N_SELE + k]*d3*d6*d6*d6)/572.0 + (Nc_n[i*N_SELE + k]*d6*d6*d6*d6)/364.0) + (486.0*Nc_n[i*N_SELE + k]*(d2 - d5)*(d2 - d5)*(d2 d5)*(d2 - d5))/(1001.0*(sel*sel*sel)) - (255.0*Nc_n[i*N_SELE + k]*d2*d3*d3*d3)/2002.0 - (3.0*Nc_n[i*N_SELE + k]*d2*d6*d6*d6)/4004.0 + (255.0*Nc_n[i*N_SELE + k]*d3*d3*d3*d5)/2002.0 + (3.0*Nc_n[i*N_SELE + k]*d5*d6*d6*d6)/4004.0 + (3.0*Nc_n[i*N_SELE + k]*(d2 - d5)*(d2 d5)*(348.0*d3*d3 - 84.0*d3*d6 + 114.0*d6*d6))/(4004.0*sel) + (9.0*Nc_n[i*N_SELE + k]*(d2 - d5)*(d2 - d5)*(d2 - d5)*(20.0*d3 + 7.0*d6))/(1001.0*(sel*sel)) - (9.0*Nc_n[i*N_SELE + k]*d2*d3*d6*d6)/4004.0 + (171.0*Nc_n[i*N_SELE + k]*d2*d3*d3*d6)/2002.0 + (9.0*Nc_n[i*N_SELE + k]*d3*d5*d6*d6)/4004.0 - (171.0*Nc_n[i*N_SELE + k]*d3*d3*d5*d6)/2002.0 + (90.0*Nc_n[i*N_SELE + k]*(d3 + d6)*(d2 - d5)*(d2 - d5)*(d2 d5))/(1001.0*(sel*sel))) + heso1*(kts33) + heso2*(ktc233+ktc333+ktc433)+ heso3*(kps233+kps333+kps433)+heso4*(ktp33); sub_sm[2][3] = Etb*S_e[i*N_SELE + k]/sel + Etb*A_e[i*N_SELE + k]*(-1/(10*sel)*(d5 - d2) - 1.0/30*(d6 - 4*d3)) +(2.0* Mc_n[i*N_SELE + k]*tc1)/sel +(4.0*Nc_n[i*N_SELE + k]*d3*tc1)/15.0 - (Nc_n[i*N_SELE + k]*d6*tc1)/15.0 (6.0*tc2*( Mc_n[i*N_SELE + k]*d4 - Mc_n[i*N_SELE + k]*d1 + 4.0*Ic_ze[i*N_SELE + k]*d3*ftb_c + 2.0*Ic_ze[i*N_SELE + k]*d6*ftb_c))/(sel*sel) + (Nc_n[i*N_SELE + k]*tc1*(d2 - d5))/(5.0*sel) - (36.0*Ic_ze[i*N_SELE + k]*ftb_c*tc2*(d2 d5))/(sel*sel*sel) +(3.0*Nc_n[i*N_SELE + k]*tc2*(d1 - d4)*(d2 - d5))/(5.0*(sel*sel)) (tc2*(15.0* Mc_n[i*N_SELE + k]*d3*d3 - 4.0*Nc_n[i*N_SELE + k]*d1*d3 + Nc_n[i*N_SELE + k]*d1*d6 + 4.0*Nc_n[i*N_SELE + k]*d3*d4 - Nc_n[i*N_SELE + k]*d4*d6))/(5.0*sel) +(3.0*Nc_n[i*N_SELE + k]*d6*d6*d6*tc2)/140.0 - (6.0*Nc_n[i*N_SELE + k]*d3*d3*d3*tc2)/35.0 - (3.0*Nc_n[i*N_SELE + k]*d3*d6*d6*tc2)/70.0 + (9.0*Nc_n[i*N_SELE + k]*d3*d3*d6*tc2)/140.0 - (27.0*Nc_n[i*N_SELE + k]*tc2*(d2 - d5)*(d2 - d5)*(d2 - d5))/(35.0*(sel*sel*sel)) - (27.0*Nc_n[i*N_SELE + k]*d3*tc2*(d2 - d5)*(d2 - d5))/(35.0*(sel*sel)) - (9.0*Nc_n[i*N_SELE + k]*tc2*(d2 - d5)*(2.0*d3*d6 - d3*d3 + d6*d6))/(140.0*sel) + heso1*(kts34) + heso2*(ktc234+ktc334+ktc434)+ heso3*(kps234+kps334+kps434)+heso4*(ktp34); sub_sm[2][4] = -6*Etb*I_ze[i*N_SELE + k]/(sel*sel) + Etb*A_e[i*N_SELE + k]*(-1/(10*sel)*(d4 - d1) + 1.0/280*(-3*(d6*d6 - d3*d3) - 6*d3*d6 + 72/sel*d3*(d5 - d2) - 108/(sel*sel)*(d5 - d2)*(d5 - d2))) - 1.0/10*N_p[i*N_SELE + k] -(12.0*Ic_ze[i*N_SELE + k]* ftb_c*tc1)/(sel*sel) -Nc_p[i*N_SELE + k]/10.0 +(Nc_n[i*N_SELE + k]*tc1*(d1 - d4))/(5.0*sel) + (36.0* ftb_c*tc2*(2*Cc_n[i*N_SELE + k]*d3 - Ic_ze[i*N_SELE + k]*d1 + Ic_ze[i*N_SELE + k]*d4))/(sel*sel*sel) + (72.0*Cc_n[i*N_SELE + k]* ftb_c*tc2*(d2 d5))/(sel*sel*sel*sel) +(42.0*Nc_n[i*N_SELE + k]*tc2*d1*d1 - 84.0*Nc_n[i*N_SELE + k]*tc2*d1*d4 108.0*Nc_n[i*N_SELE + k]*tc1*d2*d2 + 216.0*Nc_n[i*N_SELE + k]*tc1*d2*d5 + 1440.0*Ic_ze[i*N_SELE + k]* ftb_c*tc2*d3*d3 + 192.0*Ic_ze[i*N_SELE + k]* ftb_c*tc2*d3*d6 + 42.0*Nc_n[i*N_SELE + k]*tc2*d4*d4 - 108.0*Nc_n[i*N_SELE + k]*tc1*d5*d5 + 96.0*Ic_ze[i*N_SELE + k]* ftb_c*tc2*d6*d6)/(140.0*(sel*sel)) + (3.0*Nc_n[i*N_SELE + k]*d3*d3*tc1)/140.0 - (3.0*Nc_n[i*N_SELE + k]*d6*d6*tc1)/140.0 - (3.0*Nc_n[i*N_SELE + k]*d3*d6*tc1)/70.0 + (162.0*Ic_ze[i*N_SELE + k]* ftb_c*tc2*(d2 - d5)*(d2 d5))/(7.0*(sel*sel*sel*sel)) - (18.0*Nc_n[i*N_SELE + k]*d3*tc1*(d2 d5))/(35.0*sel) + (36.0*Ic_ze[i*N_SELE + k]* ftb_c*tc2*(d2 - d5)*(12.0*d3 + 5.0*d6))/(35.0*(sel*sel*sel)) +(9.0*Nc_n[i*N_SELE + k]*tc2*((sel*sel)*d3*d3 - 2.0*(sel*sel)*d3*d6 (sel*sel)*d6*d6 - 24.0*sel*d2*d3 + 24.0*sel*d3*d5 - 36.0*d2*d2 + 72.0*d2*d5 36.0*d5*d5)*(d1 - d4))/(140.0*(sel*sel*sel)) Page: 60 CFT-Fire Code +tc2*((57.0*Nc_n[i*N_SELE + k]*d6*d6*d6*d6)/8008.0 - (255.0*Nc_n[i*N_SELE + k]*d3*d3*d3*d3)/8008.0 - (9.0*Nc_n[i*N_SELE + k]*d3*d3*d6*d6)/8008.0 + (2025.0*Nc_n[i*N_SELE + k]*(d2 - d5)*(d2 - d5)*(d2 - d5)*(d2 d5))/(1001.0*(sel*sel*sel*sel)) - (3.0*Nc_n[i*N_SELE + k]*d3*d6*d6*d6)/4004.0 + (57.0*Nc_n[i*N_SELE + k]*d3*d3*d3*d6)/2002.0 + (27.0*Nc_n[i*N_SELE + k]*(d2 - d5)*(d2 - d5)*(10.0*d3*d3 + 7.0*d3*d6 + 10.0*d6*d6))/(2002.0*(sel*sel)) + (108.0*Nc_n[i*N_SELE + k]*(d2 - d5)*(d2 - d5)*(d2 - d5)*(18.0*d3 + 5.0*d6))/(1001.0*(sel*sel*sel)) + (3.0*Nc_n[i*N_SELE + k]*(2.0*d2 2.0*d5)*(116.0*d3*d3*d3 - 42.0*d3*d3*d6 + 114.0*d3*d6*d6 14.0*d6*d6*d6))/(4004.0*sel) + (27.0*Nc_n[i*N_SELE + k]*(d3 + d6)*(d2 d5)*(d2 - d5)*(20.0*d3 + 7.0*d6))/(2002.0*(sel*sel))) + heso1*(kts35) + heso2*(ktc235+ktc335+ktc435)+ heso3*(kps235+kps335+kps435)+heso4*(ktp35); sub_sm[2][5] = 2*Etb*I_ze[i*N_SELE + k]/sel + Etb*A_e[i*N_SELE + k]*(-1.0/30*(d4 - d1) + 1.0/280*(4*sel*d3*d6 - 3*sel*(d3*d3 + d6*d6) - 6*(d5 - d2)*(d3 + d6))) - sel/30*N_p[i*N_SELE + k] +(4.0*Ic_ze[i*N_SELE + k]*ftb_c*tc1)/sel -(sel*Nc_p[i*N_SELE + k])/30.0 +(Nc_n[i*N_SELE + k]*d1*tc1)/15.0 - (Nc_n[i*N_SELE + k]*d4*tc1)/15.0 (12.0*ftb_c*tc2*(Cc_n[i*N_SELE + k]*d3 - Cc_n[i*N_SELE + k]*d6 Ic_ze[i*N_SELE + k]*d1 + Ic_ze[i*N_SELE + k]*d4))/(sel*sel) +(3.0*Nc_n[i*N_SELE + k]*d2*d3*tc1)/70.0 - (tc2*(14.0*Nc_n[i*N_SELE + k]*d1*d4 - 7.0*Nc_n[i*N_SELE + k]*d1*d1 + 192.0*Ic_ze[i*N_SELE + k]*ftb_c*d3*d3 + 48.0*Ic_ze[i*N_SELE + k]*ftb_c*d3*d6 - 7.0*Nc_n[i*N_SELE + k]*d4*d4 + 192.0*Ic_ze[i*N_SELE + k]*ftb_c*d6*d6))/(70.0*sel) - sel*((3.0*Nc_n[i*N_SELE + k]*tc1*d3*d3)/140.0 - (Nc_n[i*N_SELE + k]*tc1*d3*d6)/35.0 + (3.0*Nc_n[i*N_SELE + k]*tc1*d6*d6)/140.0) + (3.0*Nc_n[i*N_SELE + k]*d2*d6*tc1)/70.0 (3.0*Nc_n[i*N_SELE + k]*d3*d5*tc1)/70.0 - (3.0*Nc_n[i*N_SELE + k]*d5*d6*tc1)/70.0 - (18.0*Ic_ze[i*N_SELE + k]*ftb_c*tc2*(d2 - d5)*(d2 d5))/(7.0*(sel*sel*sel)) - (48.0*Ic_ze[i*N_SELE + k]*ftb_c*tc2*(d3 + d6)*(d2 d5))/(35.0*(sel*sel)) -(3.0*Nc_n[i*N_SELE + k]*tc2*(d1 - d4)*(6.0*d3*d5 - 6.0*d2*d6 - 6.0*d2*d3 + 6.0*d5*d6 + 3.0*sel*d3*d3 + 3.0*sel*d6*d6 - 4.0*sel*d3*d6))/(140.0*sel) +(95.0*sel*Nc_n[i*N_SELE + k]*d3*d3*d3*d3*tc2)/8008.0 + (95.0*sel*Nc_n[i*N_SELE + k]*d6*d6*d6*d6*tc2)/8008.0 - (57.0*Nc_n[i*N_SELE + k]*d2*d3*d3*d3*tc2)/2002.0 - (57.0*Nc_n[i*N_SELE + k]*d2*d6*d6*d6*tc2)/2002.0 + (57.0*Nc_n[i*N_SELE + k]*d3*d3*d3*d5*tc2)/2002.0 + (57.0*Nc_n[i*N_SELE + k]*d5*d6*d6*d6*tc2)/2002.0 - (135.0*Nc_n[i*N_SELE + k]*tc2*(d2 - d5)*(d2 d5)*(d2 - d5)*(d2 - d5))/(1001.0*(sel*sel*sel)) + (9.0*Nc_n[i*N_SELE + k]*tc2*(d2 - d5)*(d2 - d5)*(7.0*d3*d3 - 38.0*d3*d6 + 7.0*d6*d6))/(2002.0*sel) - (sel*Nc_n[i*N_SELE + k]*d3*d6*d6*d6*tc2)/91.0 - (sel*Nc_n[i*N_SELE + k]*d3*d3*d3*d6*tc2)/91.0 + (9.0*Nc_n[i*N_SELE + k]*d2*d3*d6*d6*tc2)/4004.0 + (9.0*Nc_n[i*N_SELE + k]*d2*d3*d3*d6*tc2)/4004.0 - (9.0*Nc_n[i*N_SELE + k]*d3*d5*d6*d6*tc2)/4004.0 - (9.0*Nc_n[i*N_SELE + k]*d3*d3*d5*d6*tc2)/4004.0 (153.0*Nc_n[i*N_SELE + k]*tc2*(d3 + d6)*(d2 - d5)*(d2 - d5)*(d2 d5))/(1001.0*(sel*sel)) + (9.0*sel*Nc_n[i*N_SELE + k]*d3*d3*d6*d6*tc2)/1144.0 + heso1*(kts36) + heso2*(ktc236+ktc336+ktc436)+ heso3*(kps236+kps336+kps436)+heso4*(ktp36); sub_sm[3][3] =+A_e[i*N_SELE + k]*Etb/sel +(2.0*Nc_n[i*N_SELE + k]*tc1)/sel -(6.0*tc2*(Mc_n[i*N_SELE + k]*d3 - Mc_n[i*N_SELE + k]*d6 - Nc_n[i*N_SELE + k]*d1 + Nc_n[i*N_SELE + k]*d4))/(sel*sel) -(Nc_n[i*N_SELE + k]*tc2*(2.0*d3*d3 - d3*d6 + 2.0*d6*d6))/(5.0*sel) (18.0*Nc_n[i*N_SELE + k]*tc2*(d2 - d5)*(d2 - d5))/(5.0*(sel*sel*sel)) (3.0*Nc_n[i*N_SELE + k]*tc2*(d3 + d6)*(d2 - d5))/(5.0*(sel*sel)) + heso1*(kts44) + heso2*(ktc244+ktc344+ktc444)+ heso3*(kps144+kps244+kps344+kps444)+heso4*(ktp44); sub_sm[3][4] = +Etb*A_e[i*N_SELE + k]*(6/(5*sel*sel)*(d5 - d2)- 1/(10*sel)*(d3 + d6)) +(36.0*Ic_ze[i*N_SELE + k]*ftb_c*tc2*(d3 + d6))/(sel*sel*sel) (Nc_n[i*N_SELE + k]*tc1*(d3 + d6))/(5.0*sel) - (12.0*Nc_n[i*N_SELE + k]*tc1*(d2 - d5))/(5.0*(sel*sel)) + (72.0*Ic_ze[i*N_SELE + k]*ftb_c*tc2*(d2 - d5))/(sel*sel*sel*sel) -(3.0*Nc_n[i*N_SELE + k]*tc2*(d1 - d4)*(12.0*d2 - 12.0*d5 + sel*d3 + sel*d6))/(5.0*(sel*sel*sel)) +tc2*((216.0*Nc_n[i*N_SELE + k]*(d2 - d5)*(d2 - d5)*(d2 d5))/(35.0*(sel*sel*sel*sel)) + (27.0*Nc_n[i*N_SELE + k]*(2.0*d2 2.0*d5)*(d3*d3 + d6*d6))/(70.0*(sel*sel)) + (81.0*Nc_n[i*N_SELE + k]*(d3 + d6)*(d2 - d5)*(d2 - d5))/(35.0*(sel*sel*sel)) (3.0*Nc_n[i*N_SELE + k]*(d3 + d6)*(d3*d3 - 4.0*d3*d6 + d6*d6))/(140.0*sel)) + heso1*(kts45) + heso2*(ktc245+ktc345+ktc445)+ heso3*(kps245+kps345+kps445)+heso4*(ktp45); sub_sm[3][5] = -1*Etb*S_e[i*N_SELE + k]/sel Page: 61 CFT-Fire Code + Etb*A_e[i*N_SELE + k]*(-1/(10*sel)*(d5 - d2)- 1.0/30*(d3 - 4*d6)) -(2.0*Mc_n[i*N_SELE + k]*tc1)/sel +(4.0*Nc_n[i*N_SELE + k]*d6*tc1)/15.0 - (Nc_n[i*N_SELE + k]*d3*tc1)/15.0 (6.0*tc2*(Mc_n[i*N_SELE + k]*d1 - Mc_n[i*N_SELE + k]*d4 + 2.0*Ic_ze[i*N_SELE + k]*d3*ftb_c + 4.0*Ic_ze[i*N_SELE + k]*d6*ftb_c))/(sel*sel) + (Nc_n[i*N_SELE + k]*tc1*(d2 - d5))/(5.0*sel) - (36.0*Ic_ze[i*N_SELE + k]*ftb_c*tc2*(d2 d5))/(sel*sel*sel) +(tc2*(15.0*Mc_n[i*N_SELE + k]*d6*d6 - Nc_n[i*N_SELE + k]*d1*d3 + 4.0*Nc_n[i*N_SELE + k]*d1*d6 + Nc_n[i*N_SELE + k]*d3*d4 - 4.0*Nc_n[i*N_SELE + k]*d4*d6))/(5.0*sel) + (3.0*Nc_n[i*N_SELE + k]*tc2*(d1 - d4)*(d2 d5))/(5.0*(sel*sel)) +(3.0*Nc_n[i*N_SELE + k]*d3*d3*d3*tc2)/140.0 - (6.0*Nc_n[i*N_SELE + k]*d6*d6*d6*tc2)/35.0 + (9.0*Nc_n[i*N_SELE + k]*d3*d6*d6*tc2)/140.0 (3.0*Nc_n[i*N_SELE + k]*d3*d3*d6*tc2)/70.0 - (27.0*Nc_n[i*N_SELE + k]*tc2*(d2 - d5)*(d2 - d5)*(d2 - d5))/(35.0*(sel*sel*sel)) - (27.0*Nc_n[i*N_SELE + k]*d6*tc2*(d2 - d5)*(d2 - d5))/(35.0*(sel*sel)) - (9.0*Nc_n[i*N_SELE + k]*tc2*(d2 - d5)*(d3*d3 + 2.0*d3*d6 - d6*d6))/(140.0*sel) + heso1*(kts46) + heso2*(ktc246+ktc346+ktc446)+ heso3*(kps246+kps346+kps446)+heso4*(ktp46); sub_sm[4][4] = 12*Etb*I_ze[i*N_SELE + k]/(sel*sel*sel) + Etb*A_e[i*N_SELE + k]*(6/(5*sel*sel)*(d4 - d1) + 1.0/140*(18/sel*(d3*d3 + d6*d6) + 432/(sel*sel*sel)*(d5 - d2)*(d5 - d2) - 108/(sel*sel)*(d5 - d2)*(d3 + d6))) + 6/(5*sel)*N_p[i*N_SELE + k] +(24.0*Ic_ze[i*N_SELE + k]*ftb_c*tc1)/(sel*sel*sel) +(6.0*Nc_p[i*N_SELE + k])/(5.0*sel) -(12.0*Nc_n[i*N_SELE + k]*tc1*(d1 - d4))/(5.0*(sel*sel)) (72.0*ftb_c*tc2*(Cc_n[i*N_SELE + k]*d3 - Cc_n[i*N_SELE + k]*d6 Ic_ze[i*N_SELE + k]*d1 + Ic_ze[i*N_SELE + k]*d4))/(sel*sel*sel*sel) +(9.0*Nc_n[i*N_SELE + k]*tc1*(d3*d3 + d6*d6))/(35.0*sel) (126.0*Nc_n[i*N_SELE + k]*tc2*d1*d1 - 252.0*Nc_n[i*N_SELE + k]*tc2*d1*d4 216.0*Nc_n[i*N_SELE + k]*tc1*d2*d2 + 432.0*Nc_n[i*N_SELE + k]*tc1*d2*d5 + 216.0*Ic_ze[i*N_SELE + k]*ftb_c*tc2*d3*d3 + 180.0*Ic_ze[i*N_SELE + k]*ftb_c*tc2*d3*d6 + 126.0*Nc_n[i*N_SELE + k]*tc2*d4*d4 - 216.0*Nc_n[i*N_SELE + k]*tc1*d5*d5 + 216.0*Ic_ze[i*N_SELE + k]*ftb_c*tc2*d6*d6)/(35.0*(sel*sel*sel)) - (3888.0*Ic_ze[i*N_SELE + k]*ftb_c*tc2*(d2 - d5)*(d2 - d5))/(35.0*(sel*sel*sel*sel*sel)) + (54.0*Nc_n[i*N_SELE + k]*tc1*(d3 + d6)*(d2 - d5))/(35.0*(sel*sel)) (324.0*Ic_ze[i*N_SELE + k]*ftb_c*tc2*(d3 + d6)*(d2 d5))/(7.0*(sel*sel*sel*sel)) +(27.0*Nc_n[i*N_SELE + k]*tc2*(d1 - d4)*((sel*sel)*d3*d3 + (sel*sel)*d6*d6 + 6.0*sel*d2*d3 + 6.0*sel*d2*d6 - 6.0*sel*d3*d5 - 6.0*sel*d5*d6 + 24.0*d2*d2 48.0*d2*d5 + 24.0*d5*d5))/(35.0*(sel*sel*sel*sel)) -tc2*((14580.0*Nc_n[i*N_SELE + k]*(d2 - d5)*(d2 - d5)*(d2 - d5)*(d2 d5))/(1001.0*(sel*sel*sel*sel*sel)) + (3.0*Nc_n[i*N_SELE + k]*(29.0*d3*d3*d3*d3 - 14.0*d3*d3*d3*d6 + 57.0*d3*d3*d6*d6 - 14.0*d3*d6*d6*d6 + 29.0*d6*d6*d6*d6))/(2002.0*sel) + (324.0*Nc_n[i*N_SELE + k]*(d2 - d5)*(d2 d5)*(9.0*d3*d3 + 5.0*d3*d6 + 9.0*d6*d6))/(1001.0*(sel*sel*sel)) + (8100.0*Nc_n[i*N_SELE + k]*(d3 + d6)*(d2 - d5)*(d2 - d5)*(d2 d5))/(1001.0*(sel*sel*sel*sel)) + (27.0*Nc_n[i*N_SELE + k]*(d3 + d6)*(2.0*d2 2.0*d5)*(10.0*d3*d3 + 7.0*d3*d6 + 10.0*d6*d6))/(2002.0*(sel*sel))) + heso1*(kts55) + heso2*(ktc255+ktc355+ktc455)+ heso3*(kps255+kps355+kps455)+heso4*(ktp55); sub_sm[4][5] = -6*Etb*I_ze[i*N_SELE + k]/(sel*sel) + Etb*A_e[i*N_SELE + k]*(-1/(10*sel)*(d4 - d1) + 1.0/280*(3*(d6*d6 - d3*d3) - 6*d3*d6 + 72/sel*d6*(d5 - d2) - 108/(sel*sel)*(d5 - d2)*(d5 - d2))) - 1.0/10*N_p[i*N_SELE + k] -(12.0*Ic_ze[i*N_SELE + k]*ftb_c*tc1)/(sel*sel) -Nc_p[i*N_SELE + k]/10.0 +(Nc_n[i*N_SELE + k]*tc1*(d1 - d4))/(5.0*sel) (36.0*ftb_c*tc2*(2.0*Cc_n[i*N_SELE + k]*d6 + Ic_ze[i*N_SELE + k]*d1 Ic_ze[i*N_SELE + k]*d4))/(sel*sel*sel) - (72.0*Cc_n[i*N_SELE + k]*ftb_c*tc2*(d2 - d5))/(sel*sel*sel*sel) +(42.0*Nc_n[i*N_SELE + k]*tc2*d1*d1 - 84.0*Nc_n[i*N_SELE + k]*tc2*d1*d4 108.0*Nc_n[i*N_SELE + k]*tc1*d2*d2 + 216.0*Nc_n[i*N_SELE + k]*tc1*d2*d5 + 96.0*Ic_ze[i*N_SELE + k]*ftb_c*tc2*d3*d3 + 192.0*Ic_ze[i*N_SELE + k]*ftb_c*tc2*d3*d6 + 42.0*Nc_n[i*N_SELE + k]*tc2*d4*d4 - 108.0*Nc_n[i*N_SELE + k]*tc1*d5*d5 + 1440.0*Ic_ze[i*N_SELE + k]*ftb_c*tc2*d6*d6)/(140.0*(sel*sel)) (3.0*Nc_n[i*N_SELE + k]*d3*d3*tc1)/140.0 + (3.0*Nc_n[i*N_SELE + k]*d6*d6*tc1)/140.0 - (3.0*Nc_n[i*N_SELE + k]*d3*d6*tc1)/70.0 + (162.0*Ic_ze[i*N_SELE + k]*ftb_c*tc2*(d2 - d5)*(d2 d5))/(7.0*(sel*sel*sel*sel)) - (18.0*Nc_n[i*N_SELE + k]*d6*tc1*(d2 d5))/(35.0*sel) + (36.0*Ic_ze[i*N_SELE + k]*ftb_c*tc2*(d2 - d5)*(5.0*d3 + Page: 62 CFT-Fire Code 12.0*d6))/(35.0*(sel*sel*sel)) -(9.0*Nc_n[i*N_SELE + k]*tc2*(d1 - d4)*((sel*sel)*d3*d3 + 2.0*(sel*sel)*d3*d6 - (sel*sel)*d6*d6 + 24.0*sel*d2*d6 - 24.0*sel*d5*d6 + 36.0*d2*d2 - 72.0*d2*d5 + 36.0*d5*d5))/(140.0*(sel*sel*sel)) +tc2*((57.0*Nc_n[i*N_SELE + k]*d3*d3*d3*d3)/8008.0 - (255.0*Nc_n[i*N_SELE + k]*d6*d6*d6*d6)/8008.0 - (9.0*Nc_n[i*N_SELE + k]*d3*d3*d6*d6)/8008.0 + (2025.0*Nc_n[i*N_SELE + k]*(d2 - d5)*(d2 - d5)*(d2 - d5)*(d2 d5))/(1001.0*(sel*sel*sel*sel)) + (57.0*Nc_n[i*N_SELE + k]*d3*d6*d6*d6)/2002.0 - (3.0*Nc_n[i*N_SELE + k]*d3*d3*d3*d6)/4004.0 + (27.0*Nc_n[i*N_SELE + k]*(d2 d5)*(d2 - d5)*(10.0*d3*d3 + 7.0*d3*d6 + 10.0*d6*d6))/(2002.0*(sel*sel)) + (108.0*Nc_n[i*N_SELE + k]*(d2 - d5)*(d2 - d5)*(d2 - d5)*(5.0*d3 + 18.0*d6))/(1001.0*(sel*sel*sel)) - (3.0*Nc_n[i*N_SELE + k]*(2.0*d2 2.0*d5)*(14.0*d3*d3*d3 - 114.0*d3*d3*d6 + 42.0*d3*d6*d6 116.0*d6*d6*d6))/(4004.0*sel) + (27.0*Nc_n[i*N_SELE + k]*(d3 + d6)*(d2 d5)*(d2 - d5)*(7.0*d3 + 20.0*d6))/(2002.0*(sel*sel))) + heso1*(kts56) + heso2*(ktc256+ktc356+ktc456)+ heso3*(kps256+kps356+kps456)+heso4*(ktp56); sub_sm[5][5] = 4*Etb*I_ze[i*N_SELE + k]/sel - Etb*S_e[i*N_SELE + k]*d6 + Etb*A_e[i*N_SELE + k]*(2.0/15*(d4 - d1) + 1.0/140*(18/sel*(d5 - d2)*(d5 - d2) + sel*d3*d3 + 12*sel*d6*d6 - 3*sel*d3*d6 + 3*(d5 - d2)*(d6 - d3))) + 2*sel/15*N_p[i*N_SELE + k] +(8.0*Ic_ze[i*N_SELE + k]*ftb_c*tc1)/sel +(2.0*sel*Nc_p[i*N_SELE + k])/15.0 +(4.0*Nc_n[i*N_SELE + k]*d4*tc1)/15.0 - (4.0*Nc_n[i*N_SELE + k]*d1*tc1)/15.0 2.0*Mc_n[i*N_SELE + k]*d6*tc1 + (12.0*ftb_c*tc2*( Cc_n[i*N_SELE + k]*d3 + 5.0* Cc_n[i*N_SELE + k]*d6 + 2.0*Ic_ze[i*N_SELE + k]*d1 - 2.0*Ic_ze[i*N_SELE + k]*d4))/(sel*sel) + (72.0* Cc_n[i*N_SELE + k]*ftb_c*tc2*(d2 d5))/(sel*sel*sel) +sel*((Nc_n[i*N_SELE + k]*tc1*d3*d3)/70.0 - (3.0*Nc_n[i*N_SELE + k]*tc1*d3*d6)/70.0 + (6.0*Nc_n[i*N_SELE + k]*tc1*d6*d6)/35.0) (28.0*Nc_n[i*N_SELE + k]*tc2*d1*d1 - 56.0*Nc_n[i*N_SELE + k]*tc2*d1*d4 + 420.0*Mc_n[i*N_SELE + k]*tc2*d1*d6 - 18.0*Nc_n[i*N_SELE + k]*tc1*d2*d2 + 36.0*Nc_n[i*N_SELE + k]*tc1*d2*d5 + 24.0*Ic_ze[i*N_SELE + k]*ftb_c*tc2*d3*d3 + 384.0*Ic_ze[i*N_SELE + k]*ftb_c*tc2*d3*d6 + 28.0*Nc_n[i*N_SELE + k]*tc2*d4*d4 - 420.0*Mc_n[i*N_SELE + k]*tc2*d4*d6 - 18.0*Nc_n[i*N_SELE + k]*tc1*d5*d5 + 1536.0*Ic_ze[i*N_SELE + k]*ftb_c*tc2*d6*d6)/(70.0*sel) + (3.0*Nc_n[i*N_SELE + k]*d2*d3*tc1)/70.0 - (3.0*Nc_n[i*N_SELE + k]*d2*d6*tc1)/70.0 (3.0*Nc_n[i*N_SELE + k]*d3*d5*tc1)/70.0 + (3.0*Nc_n[i*N_SELE + k]*d5*d6*tc1)/70.0 - (216.0*Ic_ze[i*N_SELE + k]*ftb_c*tc2*(d2 - d5)*(d2 d5))/(35.0*(sel*sel*sel)) - (48.0*Ic_ze[i*N_SELE + k]*ftb_c*tc2*(d2 - d5)*(d3 + 15.0*d6))/(35.0*(sel*sel)) +3.0*Mc_n[i*N_SELE + k]*d6*d6*d6*tc2 + (3.0*Nc_n[i*N_SELE + k]*d1*d3*d3*tc2)/70.0 + (18.0*Nc_n[i*N_SELE + k]*d1*d6*d6*tc2)/35.0 (3.0*Nc_n[i*N_SELE + k]*d3*d3*d4*tc2)/70.0 - (18.0*Nc_n[i*N_SELE + k]*d4*d6*d6*tc2)/35.0 + (27.0*Nc_n[i*N_SELE + k]*tc2*(d1 - d4)*(d2 - d5)*(d2 d5))/(35.0*(sel*sel)) - (9.0*Nc_n[i*N_SELE + k]*d1*d3*d6*tc2)/70.0 + (9.0*Nc_n[i*N_SELE + k]*d3*d4*d6*tc2)/70.0 + (9.0*Nc_n[i*N_SELE + k]*tc2*(d1 d4)*(d2 - d5)*(d3 - d6))/(70.0*sel) -tc2*(sel*((Nc_n[i*N_SELE + k]*d3*d3*d3*d3)/364.0 - (3.0*Nc_n[i*N_SELE + k]*d3*d3*d3*d6)/572.0 + (3.0*Nc_n[i*N_SELE + k]*d3*d3*d6*d6)/182.0 (95.0*Nc_n[i*N_SELE + k]*d3*d6*d6*d6)/2002.0 + (285.0*Nc_n[i*N_SELE + k]*d6*d6*d6*d6)/2002.0) + (486.0*Nc_n[i*N_SELE + k]*(d2 - d5)*(d2 - d5)*(d2 d5)*(d2 - d5))/(1001.0*(sel*sel*sel)) - (3.0*Nc_n[i*N_SELE + k]*d2*d3*d3*d3)/4004.0 - (255.0*Nc_n[i*N_SELE + k]*d2*d6*d6*d6)/2002.0 + (3.0*Nc_n[i*N_SELE + k]*d3*d3*d3*d5)/4004.0 + (255.0*Nc_n[i*N_SELE + k]*d5*d6*d6*d6)/2002.0 + (3.0*Nc_n[i*N_SELE + k]*(d2 - d5)*(d2 d5)*(114.0*d3*d3 - 84.0*d3*d6 + 348.0*d6*d6))/(4004.0*sel) + (9.0*Nc_n[i*N_SELE + k]*(d2 - d5)*(d2 - d5)*(d2 - d5)*(7.0*d3 + 20.0*d6))/(1001.0*(sel*sel)) + (171.0*Nc_n[i*N_SELE + k]*d2*d3*d6*d6)/2002.0 (9.0*Nc_n[i*N_SELE + k]*d2*d3*d3*d6)/4004.0 - (171.0*Nc_n[i*N_SELE + k]*d3*d5*d6*d6)/2002.0 + (9.0*Nc_n[i*N_SELE + k]*d3*d3*d5*d6)/4004.0 + (90.0*Nc_n[i*N_SELE + k]*(d3 + d6)*(d2 - d5)*(d2 - d5)*(d2 d5))/(1001.0*(sel*sel))) + heso1*(kts66) + heso2*(ktc266+ktc366+ktc466)+ heso3*(kps266+kps366+kps466)+heso4*(ktp66); // Gan cac phan tu doi xung lai sub_sm[1][0] = sub_sm[0][1]; sub_sm[2][0] = sub_sm[0][2]; sub_sm[2][1] = sub_sm[1][2]; sub_sm[3][0] = sub_sm[0][3]; sub_sm[3][1] = sub_sm[1][3]; sub_sm[3][2] = sub_sm[2][3]; sub_sm[4][0] = sub_sm[0][4]; Page: 63 CFT-Fire Code sub_sm[4][1] = sub_sm[1][4]; sub_sm[4][2] = sub_sm[2][4]; sub_sm[4][3] = sub_sm[3][4]; sub_sm[5][0] = sub_sm[0][5]; sub_sm[5][1] = sub_sm[1][5]; sub_sm[5][2] = sub_sm[2][5]; sub_sm[5][3] = sub_sm[3][5]; sub_sm[5][4] = sub_sm[4][5]; } // end ele_stif() /********************************************************************************/ /* Thu tuc xuat thong tin ve su chay deo cua khung */ /********************************************************************************/ void print_yield(void) { double percYield[N_SELE]; bool key; key = false; // Xem co phan tu tho nao chay deo khong for (i = 0; i < n_ele; i++) { for (j = 0; j < ni_ele[i]; j++) { percYield[j] = (A_init[i*N_SELE + j] - A_e[i*N_SELE + j])/A_init[i*N_SELE + j]*100.0; if (percYield[j] > 0) { key = true; numElastic = numElastic + 1; goto print; } } } print: if (key) { // Lap qua tat ca cac cau kien for (i = 0; i < n_ele; i++) { // Xuat thong tin cot if (mt_flg[i] < 1) { // Lap qua cac phan tu for (j = 0; j < ni_ele[i]; j++) percYield[j] = (A_init[i*N_SELE + j] - A_e[i*N_SELE + j])/A_init[i*N_SELE + j]*100.0; // Xuat ket qua chay deo cua cot fprintf(out_file5, "%6d ", i + 1); for (j = 0; j < ni_ele[i]; j++) fprintf(out_file5, "%6.2lf ", percYield[j]); fprintf(out_file5, "\n"); } // end neu la phan tu cot } // end lap qua tat ca cac cau kien // Lap qua tat ca cac cau kien for (i = 0; i < n_ele; i++) { if (mt_flg[i] == 1) // xuat thong tin dam truoc { for (j = 0; j < ni_ele[i]; j++) percYield[j] = (A_init[i*N_SELE + j] - A_e[i*N_SELE + j])/A_init[i*N_SELE + j]*100.0; // Xuat ket qua chay deo cua dam fprintf(out_file5, "%6d ", i + 1); for (j = 0; j < ni_ele[i]; j++) fprintf(out_file5, "%6.2lf ", percYield[j]); fprintf(out_file5, "\n"); } // end xuat thong tin dam } // end lap qua cac cau kien } else // In goc xoay cac nut truong hop ket cau dan hoi { // Lap qua tat ca cac cau kien for (i = 0; i < n_ele; i++) { Page: 64 CFT-Fire Code if (mt_flg[i] == 1) // xuat goc xoay dau dam { fprintf(out_file5, "%6d %30.25lf %30.25lf\n", i + 1, thet_cumL[i], thet_cumR[i]); } } } } // end print_yield() /********************************************************************************/ /* Ham binh phuong */ /********************************************************************************/ double sqr(double num) { return (num*num); } /********************************************************************************/ /* Ham */ /********************************************************************************/ double min(double a, double b) {if (ab) return a; else return b; } /********************************************************************************/ /* Thu tuc xuat thong tin ve su chay deo cua khung */ /********************************************************************************/ void print_yield_c(void) { double percYield_c[N_SELE]; bool key; key = false; // Xem co phan tu tho nao chay deo khong for (i = 0; i < n_ele; i++) { for (j = 0; j < ni_ele[i]; j++) { percYield_c[j] = (Ac_init[i*N_SELE + j] - Ac_e[i*N_SELE + j])/Ac_init[i*N_SELE + j]*100.0; if (percYield_c[j] > 0) { key = true; numElastic_c = numElastic_c + 1; goto print; } } } print: if (key) { // Lap qua tat ca cac cau kien for (i = 0; i < n_ele; i++) { // Xuat thong tin cot if (mt_flg[i] < 1) { // Lap qua cac phan tu for (j = 0; j < ni_ele[i]; j++) percYield_c[j] = (Ac_init[i*N_SELE + j] - Ac_e[i*N_SELE + j])/Ac_init[i*N_SELE + j]*100.0; // Xuat ket qua chay deo cua cot fprintf(out_file8, "%6d ", i + 1); for (j = 0; j < ni_ele[i]; j++) fprintf(out_file8, "%6.2lf ", percYield_c[j]); fprintf(out_file8, "\n"); } // end neu la phan tu cot } // end lap qua tat ca cac cau kien // Lap qua tat ca cac cau kien for (i = 0; i < n_ele; i++) { Page: 65 CFT-Fire Code if (mt_flg[i] == 1) // xuat thong tin dam truoc { for (j = 0; j < ni_ele[i]; j++) percYield_c[j] = (Ac_init[i*N_SELE + j] - Ac_e[i*N_SELE + j])/Ac_init[i*N_SELE + j]*100.0; // Xuat ket qua chay deo cua dam fprintf(out_file8, "%6d ", i + 1); for (j = 0; j < ni_ele[i]; j++) fprintf(out_file8, "%6.2lf ", percYield_c[j]); fprintf(out_file8, "\n"); } // end xuat thong tin dam } // end lap qua cac cau kien } else // In goc xoay cac nut truong hop ket cau dan hoi { // Lap qua tat ca cac cau kien for (i = 0; i < n_ele; i++) { if (mt_flg[i] == 1) // xuat goc xoay dau dam { fprintf(out_file8, "%6d %30.25lf %30.25lf\n", i + 1, thet_cumL[i], thet_cumR[i]); } } } } // end print_yield() double t_to_deep(double h) { if (h < 100) return -0.000000898956*h*h*h + 0.00024702544*h*h - 0.0249790885*h + 0.994781032; else return 0.02; } double CSRF(double tempa) { int TC1 = 15000, TC2 = 800, TC8 = 570, TC64 = 100000; return 1/(1 + tempa/TC1 + pow(tempa/TC2,2)+ pow(tempa/TC8,8) + pow(tempa/TC64,64)); } double SSRF(double tempa) { int TS1 = 100000, TS2 = 900, TS8 = 555, TS64 = 100000; return 1/(1 + tempa/TS1 + pow(tempa/TS2,2)+ pow(tempa/TS8,8) + pow(tempa/TS64,64)); } double SSRFY(double tempa) { int TS1 = 100000, TS2 = 6000, TS8 = 590, TS64 = 100000; return 1/(1 + tempa/TS1 + pow(tempa/TS2,2)+ pow(tempa/TS8,8) + pow(tempa/TS64,64)); } double SSRFYY(double tempa) { if (tempa