This paper presents a fuzzy system-based modeling approach to estimate the weld bead geometry (WBG) from the welding process variables (WPVs) and to achieve a specific weld bead shape. The bacterial memetic algorithm (BMA) is applied to solve these problems in two different roles, as a supervised trainer, and as an optimizer. As a supervised trainer, the BMA is applied to tune two different WBG models. The bead geometry properties (BGP) model follows a traditional approach providing the WBG properties as outputs. The direct profile measurement (DPM) model describes the bead profiles points by a non-linear function realized in the form of fuzzy rules. As an optimizer, the BMA utilizes the developed fuzzy systems to find the solution sets of WPVs to acquire the desired WBG. The best performance is achieved by applying six rules in the BGP model and eleven rules in the DPM model. The results indicate that the normalized root means square error for the validation data set lies in the range of 0.40 − 1.56% for the BGP model and 4.49 − 7.52% for the DPM model. The comparative analysis suggests that the BGP model estimates the BWG in a superior manner when several WPVs are altered. The developed fuzzy systems provide a tool for interpreting the effects of the WPVs. The developed optimizer provides multiple valid set of WPVs to produce the desired WBG, thus supporting the selection of those process variables in applications.