## 1 Introduction  
The cell is the smallest functional unit of the mammalian body. They directly affect all of the body's activities, including movement, transmitting information, metabolism, storing nutrients, etc. Eukaryotic cells are comprised of subcellular organelles, encompass the nucleus and different polymeric filaments, and are surrounded by membrane [1]. From the perspective of mechanical behavior, the membrane, cytoskeleton, and nucleus are the most important subcellular organelles in mammalian cells. The cytoskeleton is a complex structure consisting of actin filaments, intermediate filaments, and microtubules. This structure is responsible for the cell shape, structural stability, and function [2]. Cells and their organelles are continuously exposed to a mechanical environment. Due to the mechanotransduction effects, the changes in applied mechanical loads affect biochemical signals. These signal changes can cause disorder in tissue functions. A better understanding of the behavior of subcellular organelles can help us to find the potential sources of disease and develop methods to deal with them. Cell mechanics can be studied by experimental, theoretical, and numerical methods, and can be divided into two major branches: mechanotransduction and cellular structure. Thus, it is necessary to select a suitable cellular structure as the first step in a numerical study of living cells' mechanical behavior. A number of mechanical structures for cell modeling have been introduced in recent years. Some of these use a continuum-based approach and others are based on microstructural models [3]. However, to consider the mechanical behavior of subcellular organelles, one should utilize models consisting of both continuum and structural elements. [PERSON] and [PERSON] [4] have suggested a combined model including 6-strut tensegrity (as the cytoskeleton), the membrane, nucleus, and cytoplasm. This model has been used by other researchers to investigate cellular functions. [PERSON] [5] used this model to show the connection between tissue mechanosensation and cellular stimuli. [PERSON] et al. [6] studied how the cell senses the substrate's stiffness using the same model. [PERSON] et al. [7] only considered the tensegrity structure and tried to show the effect of the complexity of structure on force distribution and energy dynamics while the cell is spreading. They demonstrated that 12-strut tensegrity provides superior results by comparing the force distribution, stored energy, and spreading area. [PERSON] et al. [8] proposed a more complicated model than the [PERSON] model [4]. They used a tensegrity structure for the nucleus and changed the cytoskeleton structure by adding centrosome and intermediate filaments to their model. They used this model to simulate a single osteocyte in the lacuna. [PERSON] et al. [9] adapted [PERSON]'s model to obtain a mechanical model for a single endothelial cell, which simulated actins in cortex and used a star pattern for modeling microtubules. In this model, subcellular elements were able to move independently, however, the number of focal adhesions were chosen arbitrarily. One cannot expect to simulate the mechanical behavior of all types of cells with a unique model, and no model can accommodate all experiment data. Nevertheless, it is expected that developing a model for as many subcellular organelles as possible could help in the consideration of other types of living cells. The aim of this study was to use a combined mechanical model consisting of more subcellular elements than previous models and investigate how these components affect cell's mechanical behavior. The main limitations and assumptions of the study were: only the linear elastic behavior of the cell was modeled and small indentations were applied during the course of the simulation to guarantee linear elastic behavior. A mechanical model including the cytoskeleton (CSK), nucleoskeleton (NSK), membrane, nuclear envelope (NE), and cytosol has been presented.