Functional properties of cells, tissues, and organs rely on predictable growth outputs. A shape change in plant cells is determined by properties of a tough cell wall that deforms anisotropically in response to high turgor pressure. This morphogenesis requires tight coordination and feedback controls among cytoskeleton-dependent wall patterning, its material properties, and stresses in the wall. Cortical microtubules bias the mechanical anisotropy of cell wall by defining the trajectories of cellulose synthase motility as they polymerize bundled microfibrils in their wake. Cortical microtubules could locally align and orient relative to cell geometry; however, the means by this orientation occurs is not known. Correlations between the microtubule orientation, cell geometry, and predicted tensile forces are strongly established, and here we simulate how different attributes of tensile force can orient and pattern the microtubule array in the cortex. We implemented a discrete model with three-state transient microtubule behaviors influenced by local mechanical stress in order to probe the mechanisms of stress-dependent patterning. We varied the sensitivity of four types of dynamic behaviors observed on the plus ends of microtubules - growth, shrinkage, catastrophe, and rescue - to local stress and then evaluated the extent and rate of microtubule alignments in a square computational domain. We optimized constitutive relationships between local stress and the plus-end dynamics and employed a biomechanically well-characterized cell wall to analyze how stress can influence the density and orientation of microtubule arrays. Our multiscale modeling approaches predict that spatial variability in stress magnitude and anisotropy mediate mechanical feedback between the wall and organization of the cortical microtubule array.