Molecular dynamics (MD) simulations have been extensively used to study lipid membranes in addition to experimental studies as they help better understand membranes in the atomic level. Computational models of bacterial (E. coli) membranes have been developed and applied to study the antimicrobial peptide proteins. Plant membranes are less frequently studied compared to the bacterial membrane. In this work, we will present the soybean plasma membranes models. The compositions of cell plasma membranes of soybean vary depending on the species, stage of development, and the part of the plant. The two parts of the plant that we study are the hypocotyl and the root. Each model consists of 100 lipids per leaflet, with the composition based on the weighted and averaged values from past experimental studies. Specifically, the hypocotyl membrane contains 7 types of unsaturated phospholipids and two types of sterols, while the root membrane contains 8 types of phospholipids and two types of sterols. All types of phospholipids in soybean contains the 18:2 (cis Δ9, 12) linoleoyl tail which was not well studied before, therefore, the simulations on the pure 18:0/18:2 and 18:2/18:2 phosphocholine (PC) lipid bilayers are also performed. The structural properties such as surface area per lipid, bilayer thicknesses, order parameters, and spin-lattice relaxation time are analyzed for all membranes. Moreover, the analyses of the sterols tilt angle distributions, hydrogen bonding, and clustering are also conducted for the soybean membranes. The structural properties of pure bilayers agree well with NMR experimental data validate the accuracy of 18:2 linoleoyl-containing lipids, based on which the soybean membrane models also result in reasonable structural properties. These results imply that the two soybean membrane models are realistic, and can facilitate the further study of soybean and other plant membranes.