We consider the doubly infinite Sierpinski gasket graph $SG_0$, rescale it by factor $2^{-n}$, and on the rescaled graphs $SG_n=2^{-n}SG_0$, for every $n\in \mathbb{N}$, we investigate the limit shape of three aggregation models with initial configuration $σ_n$ of particles supported on multiple vertices. The models under consideration are: divisible sandpile in which the excess mass is distributed among the vertices until each vertex is stable and has mass less or equal to one, internal DLA in which particles do random walks until finding an empty site, and rotor aggregation in which particles perform deterministic counterparts of random walks until finding an empty site. We denote by $SG=cl(\cup_{n=0}^{\infty} SG_n)$ the infinite Sierpinski gasket, which is a closed subset of $\mathbb{R}^2$, for which $SG_n$ represents the level-n approximating graph, and we consider a continuous function $σ:SG\to\mathbb{N}$. For $σ$ we solve the obstacle problem and we describe the noncoincidence set $D\subset SG$ as the solution of a free boundary problem on the fractal $SG$. If the discrete particle configurations $σ_n$ on the approximating graphs $SG_n$ converge pointwise to the continuous function $σ$ on the limit set $SG$, we prove that, as $n\to\infty$, the scaling limits of the three aforementioned models on $SG_n$ starting with initial particle configuration $σ_n$ converge to the deterministic solution $D$ of the free boundary problem on the limit set $SG\subset\mathbb{R}^2$. For $D$ we also investigate boundary regularity properties.