Jurnal Teknologika 15. Issue 1 . 629-642 Jurnal Teknologika Journal homepage: https://jurnal. id/index. php/teknologika/index ISSN: 2715-4645 E-ISSN: 1693-2978 A Dynamic System Model for Multi-Echelon Inventory Policy Optimization Considering Bullwhip Effect and Supply Chain Disruptions in the Sport Footwear Industry Dedy Setyo Utomo1,*. Rizky Fajar Ramdhani1. Asep Hermawan1 Teknik Industri. Sekolah Tinggi Teknologi Wastukacana. Purwakarta 41153. Indonesia Abstarct: This research develops a dynamic system model to optimize inventory policies in multi-echelon supply chains within the sport footwear industry, addressing challenges from the bullwhip effect and supply chain disruptions. The sports footwear sector faces unique inventory management challenges due to complex demand patterns influenced by seasonality, fashion trends, and competitive dynamics. Our comprehensive system dynamics model captures the intricate relationships between four supply chain echelons: retailers, distributors, manufacturers, and raw material suppliers. The model integrates machine learning algorithmsAispecifically Long Short-Term Memory (LSTM) neural networksAifor adaptive demand prediction and employs genetic algorithm optimization to determine optimal inventory parameters under various disruption scenarios. Using real-world data from a leading sports footwear manufacturer, we validated the model under normal operations and three distinct disruption scenarios: raw material shortages . % reduction for 6 week. , manufacturing capacity constraints . % reduction for 8 week. , and transportation disruptions . oubled lead times for 4 week. Results demonstrate that our proposed hybrid model reduces overall inventory costs by 18. 7% compared to traditional policies while maintaining a 97. 2% service level. The integration of machine learning for demand forecasting reduced prediction errors by 43. 6% compared to conventional methods, directly mitigating the bullwhip effect by decreasing the order variability coefficient from 0. 89 to 0. 61 at the supplier level. Furthermore, the model enhanced supply chain resilience by reducing recovery time by 42% following major disruptions. This research contributes to the theoretical understanding of complex supply chain dynamics and practical applications for inventory management in volatile industries, offering a robust framework for decision-making under uncertainty. Keywords: System dynamics. Multi-echelon inventory. Bullwhip effect. Supply chain. Sport Footwear Industries Introduction Supply chain management in the sports footwear industry presents unique challenges due to the combination of fashion-driven consumer demand, technological innovation, seasonal variations, and intense market competition. The industryAos complex multi-echelon structureAifrom raw material suppliers to retailersAiamplifies these challenges, particularly through the bullwhip effect phenomenon . This effect, characterized by increasing order variability as one moves upstream in the supply chain, leads to inefficient inventory policies, higher costs, and reduced service levels . The sports footwear industryAos vulnerability to supply chain disruptions further complicates inventory management. Recent global events, including the COVID-19 pandemic, have exposed significant vulnerabilities in global supply chains, with the sports footwear sector being particularly affected due to its reliance on global sourcing networks and just-in-time inventory systems . Thesedisruptions manifest at various echelons, from raw material shortages to manufacturing capacity constraints and logistics disruptions, creating cascading effects throughout the supply chain . Traditional inventory management approaches typically employ static models that fail to capture the dynamic interactions between supply chain echelons and the evolving nature of disruptions . * Corresponding author: dedy@wastukancana. https://doi. org/10. 51132/teknologika. Received : 24-04-2025 Accepted : 06-05-2025 Available online : 12-5-2025 Jurnal Teknologika. Volume 15. Issue 1 . 629-642 Moreover, these approaches often treat each echelon as independent, overlooking the systemic effects that decisions at one level have on others. As a result, companies frequently experience excess inventory in some categories while simultaneously facing stockouts in others, especially during . System dynamics modelling offers a promising approach to address these limitations by capturing the complex feedback mechanisms, time delays, and nonlinear relationships that characterize multiechelon supply chains. By incorporating machine learning algorithms for demand prediction and metaheuristic optimization techniques for parameter tuning, this approach can better adapt to changing conditions and optimize inventory policies across the supply chain. This research aims to develop and validate a dynamic system model for optimizing inventory policies in multi-echelon sports footwear supply chains, specifically to mitigate the bullwhip effect and enhance resilience against disruptions. The study addresses the following research questions: How can system dynamics modelling effectively capture the complexities of multi-echelon inventory management in the sport footwear industry? How do different supply chain disruptions affect inventory performance across multiple How can machine learning algorithms and metaheuristic optimization techniques enhance inventory policy decisions under conditions of uncertainty? What strategies can effectively mitigate the bullwhip effect while maintaining resilience against supply chain disruptions? Methodology 1 Research Design This study employed a mixed-method research design combining quantitative modelling with qualitative insights from industry experts. The research process consisted of four main phases: . system conceptualization and boundary definition, . model development, . data collection and model validation, and . scenario analysis and policy optimization. 2 System Conceptualization The supply chain system under investigation was conceptualized as a four-echelon network comprising retailers, distributors, manufacturers, and raw material suppliers within the sports footwear Based on a comprehensive literature review and industry expert consultations, we identified key variables affecting inventory dynamics, including: A Demand patterns . easonal fluctuations, trend changes, promotional effect. A Production capacities and constraints. A Lead times across different echelons. A Information sharing mechanisms. A Inventory policies . afety stock levels, reorder points, order quantitie. A Potential disruption sources and their characteristics. 3 System Dynamics Model Development The system dynamics model was developed using Vensim DSS software following the methodology outlined by Sterman . The model structure included: II. Stock and Flow Diagrams: Key stocks included inventory levels at each echelon, workin-progress, and backorders. Flows represented order rates, production rates, and delivery rates. Causal Loop Diagrams: Multiple feedback loops were identified, including: A Order fulfilment loops . A Inventory adjustment loops . A Production capacity adjustment loops . Jurnal Teknologika. Volume 15. Issue 1 . Demand amplification loops . Figure 1. Causal Loop Diagram Dynamic System Model for Multi-Echelon Inventory Policy Optimization Considering Bullwhip Effect and Supply Chain Disruptions in the Sport Footwear Industry Mathematical Formulations: The core inventory dynamics . were represented through differential equations governing the rate of change of inventory levels (I) at each echelon . yccyayce = yaycuyceycoycuycyce . Oe ycCycycyceycoycuycyce . yccyc Where: $Inflow_e. $ represents incoming deliveries from the upstream echelon. $Outflow_e. $ represents outgoing deliveries to the downstream echelon. The ordering policy at each echelon . incorporated an adaptive approach combining forecasted demand, inventory position, and safety stock: ycCyce . = yayce . cIycIyce Oe yaycEyce . ) Where: $O_e. $ is the order quantity at time t. $F_e. $ is the forecasted demand. $SS_e$ is the safety stock level. $IP_e. $ is the inventory position. $\alpha$ is the inventory adjustment parameter. 4 Machine Learning Integration for Demand Forecasting We integrated machine learning algorithms with the system dynamics model to enhance demand prediction capabilities. A Long Short-Term Memory (LSTM) neural network was implemented to forecast demand patterns based on historical data, incorporating factors such as: A Historical sales patterns. A Seasonal indices. A Promotional events. A Market trends. A Macroeconomic indicators. Jurnal Teknologika. Volume 15. Issue 1 . 629-642 The LSTM model was trained on three years of weekly sales data from a leading sports footwear manufacturer, with 80% used for training and 20% for validation. The model architecture consisted A An input layer with 12 neurons . epresenting 12 weeks of historical dat. A Two hidden LSTM layers with 64 neurons each. A A dropout layer . to prevent overfitting. A A dense output layer predicting demand for the next 8 weeks. The forecasting outputs from the LSTM model were then fed into the system dynamics model to inform ordering decisions at each echelon. 5 Metaheuristic Optimization for Inventory Parameters To determine optimal inventory policy parameters under different conditions, we implemented a genetic algorithm (GA) optimization approach . The GA optimized key decision variables. A Safety stock levels at each echelon. A Reorder points. A Order quantities. A Inventory adjustment parameters. The fitness function for the GA incorporated multiple objectives: Min Z = yui1 Ocyayce=1 OcycNyc=1. ) yui2 Ocyayce=1 OcycNyc=1 . cIyayce . yaAyce . ) ya yui3 Ocyce=1 OcycNyc=1. cCyayce . ycCyce . ) . Where: $HC_e$ is the holding cost at echelon e. $SC_e$ is the stockout cost at echelon e. $OC_e$ is the ordering cost at echelon e. $B_e. $ is the backorder level at time t. $w_1, w_2, w_3$ are weights reflecting the relative importance of each cost component. 6 Mathematical Formulation of the System Dynamics Model . Nomenclature Indices and Sets: $e \in E = . ,2,3,. $: Set of echelons . =retailer, 2=distributor, 3=manufacturer, 4=supplie. $t \in T$: Set of time periods. $p \in P$: Set of products. Variables: $I_. $: Inventory level of product $p$ at echelon $e$ at time $t$. $B_. $: Backlog level of product $p$ at echelon $e$ at time $t$. $O_. $: Order rate of product $p$ placed by echelon $e$ at time $t$. $S_. $: Shipment rate of product $p$ from echelon $e$ to echelon $e-1$ at time $t$. $IP_. $: Inventory position of product $p$ at echelon $e$ at time $t$. $F_. $: Forecasted demand of product $p$ at echelon $e$ for time $t$. $D_p. $: Customer demand for product $p$ at time $t$. $WIP_. $: Work-in-process inventory of product $p$ at echelon $e$ at time $t$. Jurnal Teknologika. Volume 15. Issue 1 . 629-642 Parameters: $LT_e$: Lead time for echelon $e$. $SS_. $: Safety stock level for product $p$ at echelon $e$. $ROP_. $: Reorder point for product $p$ at echelon $e$. $\alpha_e$: Inventory adjustment parameter for echelon $e$. $\beta_e$: Forecast smoothing parameter for echelon $e$. $MC_e$: Manufacturing capacity at echelon $e$. $RMA_t$: Raw material availability at time $t$. $DV_p. $: Demand variability for product $p$ at time $t$. System Dynamics Equations . Stock Equations Inventory Level Dynamics: The rate of change in inventory level . at each echelon is determined by the difference between inflow . hipments receive. and outflow . hipments sen. = ycIyce 1,ycy . c Oe yaycNyce 1 ) Oe ycIyce,ycy . yccyc For the supplier . , production replaces incoming shipments . = ycE4,ycy . c Oe yaycN4 ) Oe ycI4,ycy . yccyc Backlog Dynamics: The rate of change in Backlog . is determined by incoming orders minus outgoing shipments: yccyaAyce,ycy . yccyc = ycCyceOe1,ycy . Oe ycIyce,ycy . For the retailer . , customer demand replaces incoming orders . yccyaA1,ycy . yccyc = yaycy . Oe ycI1,ycy . Work-in-Process Dynamics . : WIP represents orders that have been placed but not yet yccycOyaycEyce,ycy . yccyc = ycCyce,ycy . Oe ycIyce 1,ycy . c Oe yaycNyce 1 ) . Flow Equations Order Rate . : Orders are determined by forecasted demand, desired inventory adjustment, and safety stock policies: ycCyce,ycy . = yayce,ycy . cIycIyce,ycy ycIycCycEyce,ycy Oe yaycEyce,ycy . ) . Where inventory position . is defined as: yaycEyce,ycy . = yayce,ycy . Oe yaAyce,ycy . ycOyaycEyce,ycy . Shipment Rate . : Shipments are constrained by inventory availability and Backlog: yaA . ycIyce,ycy . = min ( yce,ycy AEOIyc , yce,ycy AEOIy. Jurnal Teknologika. Volume 15. Issue 1 . 629-642 For the manufacturer . , production capacity adds a constraint . yaA . ycI3,ycy . = min ( 3,ycy AEOIyc , 3,ycy AEOIyc , ycAya3 ) . Production Rate . : For the supplier . , production is constrained by raw material ycE4,ycy . = min. cC4,ycy . , ycIycAyayc . ycaycoycoycuycaycaycycnycuycuycy,yc ) . Where allocation represents the proportion of raw materials allocated to product $p$ at time $t$. Auxiliary Equations Demand Forecasting . : Traditional exponential smoothing: yayce,ycy . = yuyce . ycCyceOe1,ycy . Oe yuyce ). yayce,ycy . c Oe . For the retailer . ya1,ycy . = yu1 . Oe yu1 ) . ya1,ycy . c Oe . LSTM-based Demand Forecasting . : When using machine learning, the forecast is generated by the LSTM model: ya1,ycy . = yaycIycNycA. c Oe . , yaycy . c Oe . A A yaycy . c Oe yc. , ycUyc ) . $X_t$ represents additional features, including seasonality indices, promotional events, etc. Service Level . : Service level at echelon $e$ for product $p$ is calculated as: yaAyce,ycy . ycIyayce,ycy . = 1 Oe ya . yun ycy . Where $\epsilon$ is a small constant to avoid division by zero. Disruption Modeling Raw Material Disruption . : When a raw material disruption occurs, availability is reduced: ycIycA yayc = ycIycA yaycuycuycycoycayco . Oe ycIyceycyceycycnycycycycayc . 1ycyunycNycycayc ) . Manufacturing Capacity Disruption . : When a manufacturing disruption occurs, capacity is ycAya3 = ycAya3,ycuycuycycoycayco . Oe ycIyceycyceycycnycycycoycaycuycyceycaycaycycycycnycuyci . 1ycyunycNycoycaycuycyceycaycaycycycycnycuyci ) . Transportation Disruption . : When a transportation disruption occurs, lead times increase: yaycNyce = yaycNyce,ycuycuycycoycayco . ycIyceycyceycycnycycycycycaycuyc . 1ycyunycNycycycaycuyc ) . Genetic Algorithm Optimization . The GA optimizes the following decision variables: Safety stock levels: $SS_. Reorder points: $ROP_. Inventory adjustment parameters: $\alpha_e$. Forecast smoothing parameters: $\beta_e$. Jurnal Teknologika. Volume 15. Issue 1 . The objective function to minimize is : Min Z = Ocyayce=1 OcycEyc=1 OcycNyc=1. yayayce,ycy . yayce,ycy . yaAyce,ycy . ycCyayce,ycy . ycCyce,ycy . ) Where: $HC_. $ is the holding cost per unit of product $p$ at echelon $e$. $SC_. $ is the stockout cost per unit of product $p$ at echelon $e$. $OC_. $ is the ordering cost per order of product $p$ at echelon $e$. $w_1,w_2,w_3$ are weights reflecting the relative importance of each cost component. Bullwhip Effect Quantification The bullwhip effect . is quantified by comparing the coefficient of variation (CV) of orders at different echelons . cCyce,ycy ) yaAycO yayce,ycy = yaycO . a ) . ycy Where: yaycO . cU) = yuNycu A value of $BWE_. > 1$ indicates the presence of the bullwhip effect, with higher values representing stronger demand amplification. 7 Specific Mathematical Formulation for the Sport Footwear Industry For the sports footwear industry application, we incorporate additional factors: Seasonal Demand Pattern . : Customer demand includes a seasonal component : = yaAycaycyceycy ycNycyceycuyccycy . yc ycIyceycaycycuycuycaycoycy . ycEycycuycoycuycycnycuycuycy . $Seasonal_p. $ captures the seasonal pattern and $Promotion_p. $ captures the effect of promotional activities. Product Lifecycle . : For fashion-driven sports footwear, we incorporate product lifecycle effects: = yaAycaycyceycy . ycEycycuycoycuycycnycuycuycy . $Lifecycle_p. $ follows a typical product lifecycle curve . ntroduction, growth, maturity. Material Constraints . : Component availability constraints are modelled as follows: OcycyyunycE yayco,ycy . ycE4,ycy . O ycIycA yayco,yc OAyco yun ycA, yc yun ycN Where: $C_. $ is the amount of material $m$ required for product $p$. $RMA_. $ is the availability of material $m$ at time $t$. $M$ is the set of materials. This mathematical formulation, combined with the Causal Loop Diagram, provides a comprehensive framework for analyzing and optimizing the multi-echelon inventory system in the sport footwear industry. The GA was configured with a population size of 100, crossover probability of 0. 8, mutation probability of 0. 1, and was run for 500 generations. 8 Data Collection Data for model calibration and validation were collected from multiple sources: Primary data: Two years of weekly inventory and order data from a major sport footwear manufacturer and its supply chain partners, including: Point-of-sale data from 120 retail locations. Order and inventory records from 8 regional distribution centres. Production and inventory data from 3 manufacturing facilities. Jurnal Teknologika. Volume 15. Issue 1 . 629-642 . Procurement data from 12 key raw material suppliers. Secondary data: Industry reports, academic literature, and publicly available information on sport footwear supply chains. Expert input: Semi-structured interviews with 15 supply chain professionals from the sports footwear industry to validate model assumptions and parameters. 9 Model Validation The system dynamics model was validated using several techniques: Structure validation: Expert reviews of causal loop diagrams and stock-flow structures. Parameter validation: Calibration using historical data and sensitivity analysis. Behaviour validation: Comparison of model outputs with historical patterns. Extreme condition testing: Verifying model robustness under extreme parameter values. The validation confirmed that the model accurately reproduced historical inventory dynamics and the bullwhip effect observed in the supply chain. 10 Disruption Scenario Design Three distinct disruption scenarios were developed to test the modelAos performance and optimize inventory policies: Raw material shortage: 45% reduction in raw material availability for 6 weeks. Manufacturing capacity constraint: 30% reduction in production capacity for 8 weeks. Transportation disruption: Doubling of lead times between echelons for 4 weeks. These scenarios were designed based on historical disruption events in the sports footwear industry and expert assessments of future risks. Results and Discussion 1 Baseline System Behavior The calibrated model successfully reproduced the bullwhip effect observed in the sport footwear supply chain. Figure 1 illustrates the amplification of order variability across the four echelons under normal operating conditions. Order variability, measured by the coefficient of variation, increased 32 at the retail level to 0. 89 at the raw material supplier level, confirming the presence of the bullwhip effect. Analysis of the baseline simulation revealed several key factors contributing to the bullwhip effect in the sports footwear supply chain: Demand forecasting errors: Traditional forecasting methods used by supply chain partners resulted in a mean absolute percentage error (MAPE) of 24. 6%, significantly contributing to inventory oscillations. Batch ordering practices: Minimum order quantities imposed by manufacturers and distributors led to order batching, amplifying demand variability. Price fluctuations: Promotional activities and discount periods created temporary demand spikes that propagated through the supply chain with increasing amplitude. Lead time variability: Inconsistent lead times across different product categories increased uncertainty in inventory planning. 2 Machine Learning Forecasting Performance The integrated LSTM forecasting model demonstrated superior performance compared to traditional forecasting methods. Table 1 compares forecast accuracy metrics between the LSTM model and conventional methods used by the case company. Jurnal Teknologika. Volume 15. Issue 1 . 629-642 Table 1. Comparison of Forecasting Methods Forecasting Method MAPE (%) MAD RMSE LSTM Model Moving Average Exponential Smoothing Linear Regression The LSTM model reduced forecasting error by 43. 6% compared to the currently used exponential smoothing method. This figure is derived from Table 1 in the Results and Discussion section . Machine Learning Forecasting Performanc. The table compares forecasting methods and shows that the LSTM model had a MAPE (Mean Absolute Percentage Erro. 3%, while the Exponential Smoothing method . he conventional method used by the case compan. had a MAPE of 21. The percentage improvement is calculated as: . 8% - 12. 3%) / 21. 8% = 0. 436 or 43. This represents how much the forecasting error was reduced by switching from the conventional exponential smoothing method to the LSTM machine learning model when integrated into the system dynamics This improved forecasting accuracy reduced the bullwhip effect, with the coefficient of variation at the raw material supplier level decreasing from 0. 89 to 0. These figures appear in 1 (Baseline System Behavio. 2 of the Results and Discussion. The paper states that under normal operating conditions, the coefficient of variation increased from 0. 32 at the retail level 89 at the raw material supplier level, confirming the presence of the bullwhip effect. Later, in 2, the paper mentions that when the LSTM model was integrated into the system dynamics model. Authis improved forecasting accuracy reduced the bullwhip effect, with the coefficient of variation at the raw material supplier level decreasing from 0. 89 to 0. Ay So the figures 0. 89 and 0. represent the coefficient of variation . measure of the bullwhip effec. at the supplier level before and after implementing the machine learning forecasting approach, respectively. 3 Optimized Inventory Policies The genetic algorithm optimization yielded distinct inventory policies for each echelon. Table 2 summarizes the optimized parameters compared to the baseline policies. Echelon Retailer Distributor Manufacturer Table 2. Comparison of Baseline and Optimized Inventory Policies Parameter Baseline Policy Optimized Policy Change (%) Safety Stock . Reorder Point . Order Frequency . Safety Stock . Reorder Point . Order Frequency . Safety Stock . Reorder Point . Order Frequency . Raw Material Safety Stock . Supplier Reorder Point . Order Frequency . Jurnal Teknologika. Volume 15. Issue 1 . 629-642 The optimized policies resulted in an 18. 7% reduction in total inventory costs while maintaining a service level of 97. 2% compared to the baseline policyAos 96. This figure comes directly from 3 of the Results and Discussion section (Optimized Inventory Policie. After presenting Table 2, which compares baseline and optimized inventory policies, the paper explicitly states: AuThe optimized policies resulted in an 18. 7% reduction in total inventory costs while maintaining a service level of 97. 2% compared to the baseline policyAos 96. These service level percentages also come directly from the same sentence in section 3. The service level metrics show that despite reducing inventory costs significantly . 7%), the optimized policies improved service levels slightly . rom 8% to 97. 2%). This is important because it demonstrates that the cost reductions didnAot come at the expense of customer service. The most significant improvements came from decreased safety stock levels and more frequent, smaller orders, which mitigated the bullwhip effect directly. 4 Disruption Scenario Analysis The modelAos performance was evaluated under the three disruption scenarios, comparing traditional policies, optimized static policies, and the proposed adaptive policies that leverage machine learning forecasting and dynamic parameter adjustment. 1 Raw Material Shortage Scenario The adaptive policy outperformed traditional and optimized static policies when faced with a 45% reduction in raw material availability for 6 weeks. Figure 2 shows the inventory levels at the manufacturer echelon under the three policy approaches. Adaptive policy . 4%): This proposed approach combines machine learning forecasting with dynamic parameter adjustment, which performed best in maintaining service levels during the disruption. Traditional policy . 2%): This represents the conventional inventory management approaches used before optimization, which performed worst during the disruption. Optimized static policy . 7%): This refers to policies that were optimized using the genetic algorithm but without the dynamic adaptation capabilities, which performed better than traditional policies but not as well as the fully adaptive approach. Figure 2. Inventory Level at Manufacturer During Raw Material Shortage This improvement was achieved through: Early detection of supply disruption through pattern recognition in the LSTM model. Dynamic reallocation of raw materials to priority products. Gradual adjustment of safety stock levels across all echelons. Jurnal Teknologika. Volume 15. Issue 1 . 629-642 2 Manufacturing Capacity Constraint Scenario Under a 30% reduction in production capacity for 8 weeks, the adaptive policy demonstrated superior performance in managing the production backlog. The key metrics are summarized in Table 3. Table 3. Performance Under Manufacturing Capacity Constraint Traditional Optimized Adaptive Metric Policy Static Policy Policy Maximum Backlog . 14,52 10,84 7,32 Average Service Level (%) Recovery Time . Total Cost Increase (%) The adaptive policy achieved these improvements through: Prioritization of high-margin and high-demand products. Dynamic adjustment of order quantities based on real-time capacity constraints. Coordinated inventory adjustments across echelons to smoothen the impact of capacity 3 Transportation Disruption Scenario When lead times doubled between echelons for 4 weeks, the adaptive policy significantly reduced the negative impacts on inventory performance, as shown in Figure 3. The adaptive policy reduced the bullwhip effect amplification during the disruption period, with the coefficient of variation at the raw material supplier level reaching only 1. 12 compared to 1. 78 under the traditional policy. The coefficient of variation is a standardized measure of dispersion calculated as the ratio of the standard deviation to the mean (CV = E/). In supply chain management, itAos commonly used to quantify the bullwhip effect by comparing the variability of orders at different Figure 3. Bullwhip Effect During Transportation Disruption The lower value under the adaptive policy . 12 vs 1. indicates that order variability was significantly reduced compared to the traditional policy. This 37% reduction in the coefficient of variation demonstrates that the adaptive policy was more effective at mitigating the bullwhip effect during the transportation disruption scenario. This improvement was achieved through: Anticipatory ordering based on early warning signals detected by the LSTM model. Temporary adjustment of safety stock levels to accommodate increased lead time uncertainty. Jurnal Teknologika. Volume 15. Issue 1 . 629-642 Coordinated communication protocols that reduced information distortion across echelons. 5 Sensitivity Analysis To assess the robustness of the proposed model, we conducted a sensitivity analysis on key parameters. Figure 4 illustrates how total supply chain costs vary with changes in demand variability, lead time, and forecast accuracy. The results indicate that the adaptive policy maintained its superiority across various parameter values, with the most significant advantage observed under high-demand variability conditions. The model was most sensitive to forecast accuracy, highlighting the importance of the machine learning component in the overall system performance. Conclusion This research developed and validated a dynamic system model for optimizing inventory policies in multi-echelon sports footwear supply chains, addressing the challenges posed by the bullwhip effect and supply chain disruptions. Integrating system dynamics modelling with machine learning for demand forecasting and metaheuristic optimization for parameter tuning yielded several significant First, the study demonstrated the efficacy of system dynamics modelling in capturing the complex interactions between supply chain echelons and the resulting bullwhip effect in the sport footwear By modelling the feedback mechanisms and time delays inherent in multi-echelon systems, the approach provided valuable insights into the causes and potential solutions for demand Second, the research quantified the impacts of different disruption types on inventory performance across multiple echelons. Raw material shortages were found to have the most severe long-term effects, while transportation disruptions caused the most significant short-term service level Manufacturing capacity constraints created the most challenging recovery dynamics due to the backlog accumulation effect. Third, integrating machine learning algorithms for demand forecasting significantly improved prediction accuracy, reducing forecasting errors by 43. 6% compared to traditional methods. This improvement directly contributed to bullwhip effect mitigation by lowering one of its primary causesAidemand signal distortion. Fourth, the metaheuristic optimization approach identified inventory policies that reduced total costs 7% while maintaining high service levels, demonstrating the potential for significant efficiency improvements in sport footwear supply chains through optimized parameter settings. Finally, the adaptive policy combining machine learning forecasting with dynamic parameter adjustment demonstrated superior performance under disruption scenarios, reducing recovery time by 42% compared to traditional approaches and maintaining higher service levels throughout the disruption period. These findings have important implications for both theory and practice. Theoretically, the research advances our understanding of how machine learning can be effectively integrated with system dynamics modelling to enhance supply chain resilience. From a practical perspective, the findings offer sports footwear companies a framework for developing more robust inventory policies that can withstand disruptions while minimizing costs. Future research could extend this work by incorporating additional echelons or parallel supply chains, exploring the impact of different information-sharing strategies, or applying the model to other industries with similar characteristics. Additionally, integrating blockchain technology for enhanced supply chain visibility represents a promising direction for further model development. Jurnal Teknologika. Volume 15. Issue 1 . 629-642 Acknowledgement we thank Sekolah Tinggi Teknologi Wastukancana Purwakarta for their support and facilitation in pursuing academic endeavours, particularly in this research. We also appreciate our colleagues who have contributed to this study. Reference