A new advanced method combining stability function and distributed plasticity (fiber) model has been developed using Fortran programming language to predict the nonlinear inelastic behavior of steel-concrete composite beams under static loading. The main advantage of this method is its ability to accurately analyze nonlinear behavior with only one or two elements per member, significantly improving analysis time. Additionally, local buckling in the web is also considered. The Generalized Displacement Control (GDC) algorithm, capable of analyzing beyond the limit point, will be used to solve the nonlinear equilibrium equations instead of the traditional Newton-Raphson algorithm. The element stiffness matrix is integrated through the Gauss-Lobatto numerical integration framework, while the nonlinear geometric effects P- and P- are considered using stability functions and the corresponding geometric matrix. The reliability and accuracy of the proposed method are verified by comparing the analysis results with experimental data and results from the Abaqus program. The obtained results demonstrate that the proposed method is not only accurate but also significantly improves computational efficiency and has the capability to address local buckling phenomena in the web. Therefore, this new method will provide a novel tool for practical composite beam design.Một phương pháp tiên tiến mới hợp nhất hàm ổn định và mô hình dẻo phân tán (fiber) đã được phát triển bằng ngôn ngữ lập trình Fortran để tiên đoán ứng xử phi tuyến của dầm liên hợp thép bê tông chịu tải trọng tĩnh. Lợi thế của phương pháp này là khả năng phân tích chính xác ứng xử phi tuyến chỉ với một hoặc hai phần tử trên cấu kiện, điều này giúp cải thiện đáng kể thời gian phân tích mô hình. Ngoài ra, hiện tượng mất ổn định cục bộ bản bụng cũng được xem xét. Thuật toán GDC (Generalized Displacement Control) với khả năng phân tích sau điểm cực hạn sẽ được sử dụng để giải các phương trình cân bằng phi tuyến thay thế cho thuật toán Newton-Raphson truyền thống. Ma trận độ cứng của phần tử được tích phân thông qua khung tích phân số Gauss-Lobatto trong khi ảnh hưởng phi tuyến hình học P- và P- được xem xét bằng hàm ổn định và ma trận hình học tương ứng. Độ tin cậy và tính chính xác của phương pháp đề xuất được kiểm chứng thông qua việc so sánh kết quả phân tích với kết quả thí nghiệm và kết quả từ chương trình Abaqus. Kết quả thu được chứng minh rằng phương pháp đề xuất không chỉ chính xác mà còn cải thiện đáng kể thời gian tính toán và có khả năng xem xét hiện tượng mất ổn định cục bộ ở bản bụng. Do đó, phương pháp mới này sẽ cung cấp một công cụ mới cho thực hành thiết kế dầm liên hợp.