层状超导体中磁通的输运性质以及相关机理的研究一直是凝聚态物理研究的热点。本文首先在x和y方向使用周期性边界条件,z方向使用开放性边界条件,在全部钉扎强度区域研究了静态和动态情形下3D系统中涡旋饼的结构分布,给出了涡旋饼在不同钉扎强度和不同密度下的静态结构相图;总结出了涡旋饼在外推力作用下的动力学规律及奇异特征。考虑到最接近实际的模型应该采用长周期性边界条件,我们通过降维（二维或准一维）和增加层间距的方法来减少计算量,并利用这一模型继续研究层状超导体中涡旋饼的性质,寻找第二峰效应等实验现象的新解释。近年来介观超导体中的磁通分布也引起了人们的广泛关注,不同于大块超导体,介观超导体样品尺寸很小,磁通的分布不仅由磁通间的相互作用决定,还由边界的限制决定。故不同形状的介观超导体对应的磁通分布亦不同,我们主要研究的是等腰直角三角形介观超导体。本文主要内容分为如下几部分:1、我们用数值模拟的方法研究了层状超导体中涡旋饼的静态结构分布,分析了随机钉扎引起的无序和层间（以及层内）的有序性的竞争关系。通常,对于层间有序的形成,来源于强层间耦合弱钉扎条件下涡旋饼或磁通线片段的耦合,即形成3D状态。相反,2D状态则是由弱层间耦合强钉扎条件下退耦合的独立涡旋饼形成。对层间有序性,随着钉扎力的增加,层内的结构涉及到晶体（Crystal）,布拉格玻璃（Bragg Glasses,BG）,磁通玻璃（Vortex Glasses,VG）和液体状（Liquid-like,LQ）结构。改变磁通密度,我们发现在同一层内部和层与层之间都存在一个先快速无序然后缓慢有序的过程,这可能与第二峰效应有关联,我们对这种非单调行为的原因也做了讨论。研究结果总结为"层间耦合强度Sm和钉扎强度fp "平面内的相图。2、使用分子动力学模拟的方法研究了外力驱动下的层状超导体中的涡旋饼。我们发现,随着外加驱动力的增加,对于强的层间耦合,预先存在的磁通线要么直接脱钉,要么在脱钉之前先转变为二维（2D）钉扎态,是哪种情况要取决于钉扎的强弱;在一个很窄的钉扎范围内,我们还发现了有趣的再次钉扎过程,它导致了负微分电阻现象的出现。对于弱的层间耦合,随外加驱动的增加,被独立钉扎的涡旋饼先是形成无序的2D流动态,然后转变为有序的3D流动态;然而,对于极强的钉扎情形,随机钉扎导致一个类似热涨落力的作用效果出现,这一作用可以融化3D磁通线,导致一个持续的2D流动态出现在快速流动区,而且,在它的中部区域关联函数有一个尖峰现象出现:随外加驱动力增加,运动的涡旋饼先晶体化形成运动的3D磁通线,然后这些3D磁通线融化,导致重新形成2D流动态。我们的结果总结成了一个动力学相图。3、我们采用了最接近真实系统的长周期性边界条件的模型,但由于计算量大,用了两种方法修改模型:降维的方法（保留z方向自由度,将x-y平面的维度缩减为一维）和增加层间距（同时微调涡旋饼相互作用）的方法。对于降维后的系统,不再是三维系统（或传统二维x-y平面）的三角格子的分布,这是因为去掉了y方向的维度,涡旋饼的有序排列为平行于z轴方向的等间距直线。我们计算了不同钉扎强度、不同层数时,稳定态涡旋饼分布的关联函数随磁通密度的变化,并观察磁通分布图,发现随磁通密度的增大,系统会由有序的1D状态退耦合为无序的2D状态,分析原因是在高磁通密度下,晶格常数变小,相互作用势的竞争将导致退耦合态与耦合态的能量差降低,导致了退耦合相变,这个退耦合与实验上观察到的第二峰效应相符合。4、我们研究了等腰直角三角形介观超导体中的磁通分布。我们使用格林函数方法寻找伦敦方程的解析解,得到了任意磁通分布的吉布斯自由能表达式,并使用镜像方法把磁通-磁通和磁通-边界之间的相互作用转化为磁通-反磁通之间的作用,然后我们结合分子动力学方法,模拟退火获得磁通涡旋的稳定态分布。研究发现随外磁场增加,磁通涡旋的填充规则如下:其分布按一层层"壳"状排列,在小涡旋度时,磁通呈轴对称分布;随涡旋度的增大,磁通开始出现不对称的分布;涡旋度越大,磁通的分布模式越复杂。另外我们还研究了系统的尺寸效应,改变系统尺寸,磁通的分布模式也会有变化,我们的结果总结为一张表。