This paper proposes a topology optimization using COMSOL which can conduct multiphysics analysis. First, a concept of topology optimization and the relaxation scheme for the design domain using the density method are briefly discussed. Next, the optimization problem for the stiffness maximization problem is formulated. An optimization algorithm is constructed based on this formulation. The details in numerical implementations using COMSOL are discussed. Finally, a numerical example is provided in order to confirm the usefulness of the proposed method.