Recent studies have published breakthroughs in the application of finite element (FEA) studies in the design and analysis of advanced orthodontics. However, FEA has not captured bone remodeling responses to advanced orthodontics. The results of these simulations report unrealistic displacement around the nasal bridge, which impeded correlation with clinical data. Bone remodeling has been previously documented in FEA and has shown bone response to mechanical stimulus in femur bone models. However, the relationship between mechanical stimulus and bone remodeling has not been reported in orthodontic studies due to the complexity of the skull. In the current study, strain energy is used as the mechanical stimulus to control remodeling, from which density and modulus evolve. Due to the localization of forces in orthodontics, current remodeling algorithms have limited application. In turn, we developed an algorithm that dynamically collects, sorts, and bins stresses in all elements for regional remodeling based on the proximity of the element to the load. The results demonstrate that bone response to orthodontic appliances is different than that of an FEA without bone remodeling, due to load path changes based upon evolution of the bone properties. It was also found that density and moduli proximal to the load application site exhibit faster remodeling than those located remotely. Modeling another biomechanical phenomena, a 3D simulation was created to simulate recent experimental results that discovered a difference in impact mitigation properties of dense-polymer/foam bilayer structure based on the orientation of the dense-polymer with respect to the impact site. The impact energy transmitted varied in time of arrival and amplitude depending on the orientation of the structure (thin layer up or down). By creating a 3D explicit dynamic FEA simulation, it is expected to reduce costly experiments and time consumed in set up, and offer opportunities for optimization for future applications in armor. The results agreed with the experimental results, displaying a delay in impact wave arrival, depending on the orientation of the structure. The FEA revealed also revealed that mid-body strains showed an increase at different time intervals, indicating the dense polymer’s engagement and impact mitigation.