Silkworms can produce strong and tough fibers at room temperature and from an aqueous solution. Therefore, it seems useful to study the mechanism of fiber formation by silkworms for development of synthetic polymers with excellent mechanical properties. The rheological behaviors of native silk dopes stored in the silk glands of Bombyx mori and Samia cynthia ricini were clarified, and flow simulations of the dopes in each spinneret were performed with a Finite Element Method. Dynamic viscoelastic measurements revealed that silk fibroin stored in silk glands forms a transient network at room temperature, and that the molecular weight for the network node corresponds to the molecular weight of a heterodimer of H-chain and L-chain (B. mori) and a homodimer of H-chains (S. c. ricini), respectively. Also, each dope exhibited zero-shear viscosity and then shear thinning like polymer melts. In addition, shear thickening due to flow-induced crystallization was observed. The critical shear rate for crystallization of B. mori dopes was smaller than that of S. c. ricini dopes. From the flow simulation, it is suggested that domestic and wild silkworms are able to crystallize the dopes in the stiff plate region by controlling shear rate using the same magnitude of extrusion pressure despite differences in rheological properties.