VASP及第一性原理电子结构计算

VASP参数设置详解

VASP参数设置详解 计算材料2010-11-30 20:11:32 阅读197 评论0 字号：大中小订阅 转自小木虫，略有增减 软件主要功能： 采用周期性边界条件（或超原胞模型）处理原子、分子、团簇、纳米线（或管）、薄膜、晶体、准晶和无定性材料，以及表面体系和固体 l 计算材料的结构参数（键长、键角、晶格常数、原子位置等）和构型 l 计算材料的状态方程和力学性质（体弹性模量和弹性常数） l 计算材料的电子结构（能级、电荷密度分布、能带、电子态密度和ELF） l 计算材料的光学性质 l 计算材料的磁学性质 l 计算材料的晶格动力学性质（声子谱等） l 表面体系的模拟（重构、表面态和STM模拟） l 从头分子动力学模拟 l 计算材料的激发态（GW准粒子修正） 计算主要的四个参数文件：INCAR ,POSCAR,POTCAR ,KPOINTS,下面简要介绍，详细权威的请参照手册 INCAR文件： 该文件控制VASP进行何种性质的计算，并设置了计算方法中一些重要的参数，这些参数主要包括以下几类： 对所计算的体系进行注释：SYSTEM

●定义如何输入或构造初始的电荷密度和波函数：ISTART，ICHARG，INIWAV ●定义电子的优化 –平面波切断动能和缀加电荷时的切断值：ENCUT，ENAUG –电子部分优化的方法：ALGO，IALGO，LDIAG –电荷密度混合的方法：IMIX，AMIX，AMIN，BMIX，AMIX_MAG，BMIX_MAG，WC，INIMIX，MIXPRE，MAXMIX –自洽迭代步数和收敛标准：NELM，NELMIN，NELMDL，EDIFF ●定义离子或原子的优化 –原子位置优化的方法、移动的步长和步数：IBRION，NFREE，POTIM，NSW –分子动力学相关参数：SMASS，TEBEG，TEEND，POMASS，NBLOCK，KBLOCK，PSTRESS –离子弛豫收敛标准：EDIFFG ●定义态密度积分的方法和参数 –smearing方法和参数：ISMEAR，SIGMA –计算态密度时能量范围和点数：EMIN，EMAX，NEDOS –计算分波态密度的参数：RWIGS，LORBIT ●其它 –计算精度控制：PREC –磁性计算：ISPIN，MAGMOM，NUPDOWN –交换关联函数：GGA，VOSKOWN –计算ELF和总的局域势：LELF，LVTOT –结构优化参数：ISIF –等等。 主要参数说明如下： ?SYSTEM：该输入文件所要执行的任务的名字。取值：字符串，缺省值：SYSTEM ?NWRITE：输出内容详细程度。取值：0~4，缺省值：2

结构设计原理复习题 及答案.

结构设计原理复习题 一、选择题 1、混凝土强度等级按照（ ）确定 A 、立方体抗压强度标准值 B 、立方体抗压强度平均值 C 、轴心抗压强度标准值 D 、轴心抗压强度设计值 2、同一强度等级的混凝土，各种强度之间的关系是（ ） A 、c f ＞cu f ＞t f B cu f ＞t f ＞c f C 、cu f ＞c f ＞t f D 、t f ＞cu f ＞c f 3、在测定混凝土立方体抗压强度时，《桥规》（JTG D —2004）采用的标准试件尺寸为（ ） 的立方体。 A 、mm 100 B 、mm 150 C 、mm 180 D 、mm 200 4、混凝土棱柱体抗压强度用符号（ ）表示 A 、c f B 、cu f C 、t f D 、s f 5、分别用mm 150和mm 200的立方体试件进行抗压强度试验，测得的抗压强度值为（ ） A 、mm 150的立方体低于mm 200的立方体 ； B 、mm 150的立方体高于mm 200的立方体 ； C 、mm 150的立方体等于mm 200的立方体 ； D 、mm 150的立方体低于mm 200的立方体 ，是因为试件尺寸越小，抗压强度就越小； 6、同一强度等级的混凝土，棱柱体试件的抗压强度与立方体试件的抗压强度关系是（ ） A 、立方体抗压强度与棱柱体抗压强度相等 B 、立方体抗压强度高于棱柱体抗压强度 C 、立方体抗压强度低于棱柱体抗压强度 D 、无法确定 7、混凝土双向受压时，其强度变化规律是（ ） A 、一向混凝土强度随着另一向压应力的增加而增加 B 、一向混凝土强度随着另一向压应力的增加而减小 C 、双向受压强度与单向受压强度相等 D 、双向受压强度低于单向受压强度 8、混凝土弹性模量的基本测定方法是（ ） Ａ、在很小的应力（c c f 3.0≤σ）下做重复加载卸载试验所测得 Ｂ、在很大的应力（c σ＞c f 5.0）下做重复加载卸载试验所测得 Ｃ、应力在0=c σ～c f 5.0 之间重复加载卸载5~10次，取c σ=c f 5.0时所测得的变形值作为混凝土弹性模量的依据 Ｄ、以上答案均不对 9、混凝土的线性徐变是指徐变变形与（ ）成正比。 Ａ、混凝土强度 Ｂ、时间 Ｃ、温度和湿度 Ｄ、应力 10、《公路桥规》中规定了用于公路桥梁承重部分混凝土标号分为（ ）等级。 Ａ、8 Ｂ、10 Ｃ、12 Ｄ、13 11、在按极限状态理论计算钢筋混凝土构件承载力时，对于有明显流幅的钢筋，原则上都是以（ ）作为钢筋强度取值的依据 Ａ、屈服极限 Ｂ、比例极限 Ｃ、弹性极限 Ｄ、抗拉极限强度 12、对于无明显流幅的钢筋，结构设计时原则上都是以（ ）作为钢筋强度取值的依据 Ａ、比例极限 Ｂ、条件屈服强度 Ｃ、弹性极限 Ｄ、抗拉极限强度 13、钢筋和混凝土材料的强度设计值（ ）强度标准值。 Ａ、等于 Ｂ、小于 Ｃ、大于 Ｄ、不确定 14、钢筋的塑性变形性能通常用（ ）来衡量。 Ａ、屈服极限和冷弯性能 Ｂ、比例极限和延伸率 Ｃ、延伸率和冷弯性能 Ｄ、抗拉极限强度和延伸率

如何用VASP计算单个原子的能量和能级

氢原子的能量为-13.6eV在这一节中，我们用V ASP计算H原子的能量。 对于原子计算，我们可以采用如下的INCAR文件 PREC=ACCURATE NELMDL = 5 make five delays till charge mixing ISMEAR = 0; SIGMA=0.05 use smearing method 采用如下的KPOINTS文件。由于增加K点的数目只能改进描述原子间的相互作用，而在单原子计算中并不需要。所以我们只需要一个K点。 Monkhorst Pack Monkhorst Pack 1 1 1 0 0 0 采用如下的POSCAR文件 atom 1 15.00000 .00000 .00000 .00000 15.00000 .00000 .00000 .00000 15.00000 1 cart 0 0 0 采用标准的H的POTCAR 得到结果如下： k-point 1 : 0.0000 0.0000 0.0000 band No. band energies occupation 1 -6.3145 1.00000 2 -0.0527 0.00000 3 0.4829 0.00000 4 0.4829 0.00000 我们可以看到，电子的能级不为-13.6eV。 Free energy of the ion-electron system (eV) ---------------------------------------------------

alpha Z PSCENC = 0.00060791 Ewald energy TEWEN = -1.36188267 -1/2 Hartree DENC = -6.27429270 -V(xc)+E(xc) XCENC = 1.90099128 PAW double counting = 0.00000000 0.00000000 entropy T*S EENTRO = -0.02820948 eigenvalues EBANDS = -6.31447362 atomic energy EATOM = 12.04670449 --------------------------------------------------- free energy TOTEN = -0.03055478 eV energy without entropy = -0.00234530 energy(sigma->0) = -0.01645004 我们可以看到TOTEN-EA TOM也不等于-13.6eV。 在上面的计算中有个问题，就是H原子有spin,而在上面的计算中我们并没有考虑到spin。所以如果我们改用LSDA近似，在INCAR中用ISPIN=2的tag，则得到如下结果： k-point 1 : 0.0000 0.0000 0.0000 band No. band energies occupation 1 -7.2736 1.00000 2 -0.1229 0.00000 3 0.4562 0.00000 4 0.4562 0.00000 5 0.4562 0.00000 spin component 2 k-point 1 : 0.0000 0.0000 0.0000 band No. band energies occupation 1 -2.4140 0.00000 2 -0.0701 0.00000 3 0.5179 0.00000 4 0.5179 0.00000 5 0.5179 0.00000 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 0.00060791 Ewald energy TEWEN = -1.36188267 -1/2 Hartree DENC = -6.68322940 -V(xc)+E(xc) XCENC = 2.38615430 PAW double counting = 0.00000000 0.00000000

VASP控制参数文件INCAR的简单介绍

限于能力，只对部分最基本的一些参数(>,没有这个标志的参数都是可以不出现的) 详细说明，在这里只是简单介绍这些参数的设置，详细的问题在后文具体示例中展开。 部分可能会干扰VASP运行的参数在这里被刻意隐去了，需要的同学还是请查看VASP自带的帮助文档原文。 参数列表如下： >SYSTEM name of System 任务的名字*** >NWRITE verbosity write-flag (how much is written) 输出内容详细程度0-3 缺省2 如果是做长时间动力学计算的话最好选0或1(首末步/每步核运动输出) 据说也可以结合shell的tail或grep命令手动输出 >ISTART startjob: restart选项0-3 缺省0/1 for 无/有前次计算的WAVECAR(波函数) 1 'restart with constant energy cut-off' 2 'restart with constant basis set' 3 'full restart including wave function and charge prediction' ICHARG charge: 1-file 2-atom 10-const Default:if ISTART=0 2 else 0 ISPIN spin polarized calculation (2-yes 1-no) default 2 MAGMOM initial mag moment / atom Default NIONS*1 INIWAV initial electr wf. : 0-lowe 1-rand Default 1 only used for start jobs (ISTART=0) IDIPOL calculate monopole/dipole and quadrupole corrections 1-3 只计算第一/二/三晶矢方向适于slab的计算 4 全部计算尤其适于就算孤立分子 >PREC precession: medium, high or low(VASP.4.5+ also: normal, accurate) Default: Medium VASP4.5+采用了优化的accurate来替代high,所以一般不推荐使用 high。不过high可以确保'绝对收敛'，作为参考值有时也是必要的。 同样受推荐的是normal，作为日常计算选项，可惜的是说明文档提供的信息不足。 受PREC影响的参数有四类：ENCUT; NGX,NGY,NGZ; NGXF, NGYF, NGZF; ROPT 如果设置了PREC，这些参数就都不需要出现了 当然直接设置相应的参数也是同样效果的，这里不展开了，随后详释

结构设计原理习题-练习

《结构设计原理》复习题 一、填空 1．按加工方式不同，钢筋分为（）、（）、（）、（）四种。2．（）与（）通常称为圬工结构。 3．梁内钢筋主要有（）、（）、（）、（）等。 4．随着柱的长细比不同，其破坏型式有（）、（）两种。 5．根据张拉预应力筋与浇筑混凝土构件之间的先后顺序，预应力混凝土分为（）、（）两类。 6．钢筋与混凝土之间的粘结力主要有以下三项组成（）、（）、（）。7．按照配筋多少的不同，梁可分为（）、（）、（）三种。 8．钢筋混凝土受弯构件主要有（）和（）两种形式。 9．梁内钢筋主要有（）、（）、（）、（）等。 10．（）、（）、（）称为结构的可靠性。 11．钢筋的冷加工方法有（）、（）、（）三种。 12．结构的极限状态，根据结构的功能要求分为（）、（）两类。 13．T形截面梁的计算，按（）的不同分为两种类型。 14．在预应力混凝土中，对预应力有如下的要求（）、（）、（）。15．钢筋混凝土梁一般有（）、（）、（）三种不同的剪切破坏形式。16．预应力钢筋可分为（）、（）、（）三种。 二、判断题：（正确的打√,错误的打×。） 1．混凝土在长期荷载作用下，其变形随时间延长而增大的现象称为徐变。（）2．抗裂性计算的基础是第Ⅱ阶段。（）3．超筋梁的破坏属于脆性破坏，而少筋梁的破坏属于塑性破坏。（）4．增大粘结力、采用合理的构造和高质量的施工、采用预应力技术可以减小裂缝宽度。（）5．当剪跨比在[1, 3]时，截面发生斜压破坏。. （）6．预应力损失是可以避免的。（）7．整个结构或结构的一部分，超过某一特定状态时，就不能满足结构功能的要求，这种特殊状态称为结构的极限状态。（）8．箍筋的作用主要是与纵筋组成钢筋骨架，防止纵筋受力后压屈向外凸出。（） 9．采用预应力技术可杜绝裂缝的发生或有效减少裂缝开展宽度。（）10．为了保证正截面的抗弯刚度，纵筋的始弯点必须位于按正截面的抗弯计算该纵筋的强度全部被发挥的截面以内，并使抵抗弯矩位于设计弯矩图的里面。（）11．偏心距增大系数与偏心距及构件的长细比有关。（）12．钢筋混凝土梁的刚度是沿梁长变化的，无裂缝区段刚度小，有裂缝区段刚度大。（）13．钢筋按其应力应变曲线分为有明显流幅的钢筋和没有明显流幅的钢筋。（）14．因为钢筋的受拉性能好，所以我们只在受拉区配置一定数量的钢筋而在受压区不配置钢筋。（）15．当轴向力的偏心较小时，全截面受压，称为小偏心受压。（） 越大越好。（）16．有效预应力 pe

结构设计原理知识点

第一章 钢筋混凝土结构基本概念及材料的物理力学性能 1.混凝土立方体抗压强度cu f ：（基本强度指标）以边长150mm 立方体试件，按标准方法制作养护28d ，标准试验方法（不涂润滑剂，全截面受压，加载速度0.15~0.25MPa/s ）测得的抗压强度作为混凝土立方体抗压强度 cu f 。 影响立方体强度主要因素为试件尺寸和试验方法。尺寸效应关系： cu f （150）=0.95cu f （100） cu f （150）=1.05cu f （200） 2.混凝土弹性模量和变形模量。 ①原点弹性模量：在混凝土受压应力—应变曲线图的原点作切线，该切线曲率即为原点弹性模量。表示为：E '=σ/ε=tan α0 ②变形模量：连接混凝土应力应变—曲线的原点及曲线上某一点K 作割线，K 点混凝土应力为σc （=0.5c f ），该割线（OK ）的斜率即为变形模量，也称割线模量或弹塑性模量。 E c '''=tan α1=σc /εc 混凝土受拉弹性模量与受压弹性模量相等。 ③切线模量：混凝土应力应变—上某应力σc 处作一切线，该切线斜率即为相应于应力σc 时的切线模量''c E =d σ/d ε 3 . 徐变变形：在应力长期不变的作用下，混凝土的应变随时间增长的现象称为徐变。 影响徐变的因素：a. 内在因素，包括混凝土组成、龄期，龄期越早，徐变越大；b. 环境条件，指养护和使用时的温度、湿度，温度越高，湿度越低，徐变越大；c. 应力条件，压应力σ﹤0.5 c f ，徐变与应力呈线性关系；当压应力σ介于（0.5～0.8）c f 之间，徐变增长比应力快；当压应力σ﹥0.8 c f 时，混凝土的非线性徐变不收敛。 徐变对结构的影响：a.使结构变形增加；b.静定结构会使截面中产生应力重分布；c.超静定结构引起赘余力；d.在预应力混凝土结构中产生预 应力损失。 4.收缩变形：在混凝土中凝结和硬化的物理化学过程中体积随时间推移而减少的现象称为收缩。 混凝土收缩原因：a.硬化初期，化学性收缩，本身的体积收缩；b.后期，物理收缩，失水干燥。 影响混凝土收缩的主要因素：a.混凝土组成和配比；b.构件的养护条件、使用环境的温度和湿度，以及凡是影响混凝土中水分保持的因素；c.构件的体表比，比值越小收缩越大。 混凝土收缩对结构的影响：a.构件未受荷前可能产生裂缝；b.预应力构件中引起预应力损失；c.超静定结构产生次内力。 5.钢筋的基本概念 1.钢筋按化学成分分类，可分为碳素钢和普通低合金钢。 2钢筋按加工方法分类，可分为a.热轧钢筋；b.热处理钢筋；c.冷加工钢筋（冷拉钢筋、冷轧钢筋、冷轧带肋钢筋和冷轧扭钢筋。） 6.钢筋的力学性能 物理力学指标：（1）两个强度指标：屈服强度，结构设计计算中强度取值主要依据；极限抗拉强度，材料实际破坏强度，衡量钢筋屈服后的抗拉能力，不能作为计算依据。（2）两个塑性指标：伸长率和冷弯性能：钢材在冷加工过程和使用时不开裂、弯断或脆断的性能。 7.钢筋和混凝土共同工作的的原因：（1）混凝土和钢筋之间有着良好的黏结力；（2）二者具有相近的温度线膨胀系数；（3）在保护层足够的前提下，呈碱性的混凝土可以保护钢筋不易锈蚀，保证了钢筋与混凝土的共同作用。 第二章 结构按极限状态法设计计算的原则 1.结构概率设计的方法按发展进程划分为三个水准：a.水准Ⅰ，半概率设计法，只对影响结构可靠度的某些参数，用数理统计分析，并与经验结合，对结构的可靠度不能做出定量的估计；b.水准Ⅱ，近似概率设计法，用概率论和数理统计理论，对结构、构件、或截面设计的可靠概率做出近似估计，忽略了变量随时间的关系，非线性极限状态方程线性化；c.水准Ⅲ，全概略设计法，我国《公桥规》采用水准Ⅱ。 2.结构的可靠性：指结构在规定时间（设计基准期）、规定的条件下，完成预定功能的能力。 可靠性组成：安全性、适用性、耐久性。 可靠度：对结构的可靠性进行概率描述称为结构可靠度。 3.结构的极限状态：当整个结构或构件的一部分超过某一特定状态而不能满足设计规定的某一功能要求时，则此特定状态称为该功能的极限状态。 极限状态分为承载能力极限状态、正常使用极限状态和破坏—安全状态。 承载能力极限状态对应于结构或构件达到最大承载力或不适于继续承载的变形，具体表现：a.整个构件或结构的一部分作为刚体失去平衡；b.结构构件或连接处因超过材料强度而破坏；c.结构转变成机动体系；d.结构或构件丧失稳定；e.变形过大，不能继续承载和使用。 正常使用极限状态对应于结构或构件达到正常使用或耐久性能的某项规定限值，具体表现：a.由于外观变形影响正常使用；b.由于耐久性能的局部损坏影响正常使用；c.由于震动影响正常使用；d.由于其他特定状态影响正常使用。 破坏—安全状态是指偶然事件造成局部损坏后，其余部分不至于发生连续倒塌的状态。（破坏—安全极限状态归到承载能力极限状态中） 4.作用：使结构产生内力、变形、应力、应变的所有原因。 作用分为：永久作用、可变作用和偶然作用。 永久作用：在结构使用期内，其量值不随时间变化，或其变化与平均值相比可忽略不计的作用 可变作用：在结构试用期内，其量值随时间变化，且其变化值与平均值相比较不可忽略的作用。

VASP-INCAR参数设置

V A S P-I N C A R参数设置-CAL-FENGHAI-(2020YEAR-YICAI)_JINGBIAN

1. 结构优化 (Opt) SYSTEM = opt ISTART = 0 INIWAV = 1 ICHARG = 2 ISPIN = 2 LREAL = Auto ENCUT = 400 PREC = high NSW= 600 NELM = 60 IBRION = 2 ISIF = 2 POTIM = 0.1 ALGO= Fast LVDW = .TRUE. EDIFF = 1E-5 EDIFFG = 1E-4 or -0.05 # 体系需计算TS时，全部结构优化EDIFFG均设置为-0.05 ISMEAR = 0 SIGMA = 0.2 LCHARG = .FALSE. LWAVE = .FALSE.

2. 过渡态搜索 (TS): 计算时先进行低精度计算，再进行高精度计算 SYSTEM= TS ISTART = 0 INIWAV = 1 ICHARG = 2 ISPIN = 2 LREAL = Auto ENCUT = 400 PREC = high NSW = 600 NELMIN = 6 IBRION = 3 or 1 # 过渡态计算低精度为3，高精度为1 ISIF = 2 POTIM = 0.01 ALGO = Fast LVDW = .TRUE. EDIFF = 1E-5 EDIFFG = -1 or -0.05 # 过渡态计算低精度为-1，高精度为-0.05 ISMEAR = 0 SIGMA = 0.05 LCHARG= .FALSE. LWAVE= .FALSE. IMAGES=8 # TS专属设置 SPRING=-5 # TS专属设置 LCLIMB=.TRUE. # TS专属设置

vasp计算参数设置

软件主要功能： 采用周期性边界条件（或超原胞模型）处理原子、分子、团簇、纳米线（或管）、薄膜、晶体、准晶和无定性材料，以及表面体系和固体 l 计算材料的结构参数（键长、键角、晶格常数、原子位置等）和构型 l 计算材料的状态方程和力学性质（体弹性模量和弹性常数） l 计算材料的电子结构（能级、电荷密度分布、能带、电子态密度和ELF） l 计算材料的光学性质 l 计算材料的磁学性质 l 计算材料的晶格动力学性质（声子谱等） l 表面体系的模拟（重构、表面态和STM模拟） l 从头分子动力学模拟 l 计算材料的激发态（GW准粒子修正） 计算主要的四个参数文件：INCAR ,POSCAR,POTCAR ,KPOINTS,下面简要介绍，详细权威的请参照手册 INCAR文件： 该文件控制VASP进行何种性质的计算，并设置了计算方法中一些重要的参数，这些参数主要包括以下几类： l 对所计算的体系进行注释：SYSTEM l 定义如何输入或构造初始的电荷密度和波函数：ISTART，ICHARG，INIWA V l 定义电子的优化 –平面波切断动能和缀加电荷时的切断值：ENCUT，ENAUG –电子部分优化的方法：ALGO，IALGO，LDIAG –电荷密度混合的方法：IMIX，AMIX，AMIN，BMIX，AMIX_MAG，BMIX_MAG，WC，INIMIX，MIXPRE，MAXMIX –自洽迭代步数和收敛标准：NELM，NELMIN，NELMDL，EDIFF l 定义离子或原子的优化 –原子位置优化的方法、移动的步长和步数：IBRION，NFREE，POTIM，NSW –分子动力学相关参数：SMASS，TEBEG，TEEND，POMASS，NBLOCK，KBLOCK，PSTRESS –离子弛豫收敛标准：EDIFFG l 定义态密度积分的方法和参数 –smearing方法和参数：ISMEAR，SIGMA –计算态密度时能量范围和点数：EMIN，EMAX，NEDOS –计算分波态密度的参数：RWIGS，LORBIT l 其它 –计算精度控制：PREC –磁性计算：ISPIN，MAGMOM，NUPDOWN –交换关联函数：GGA，VOSKOWN –计算ELF和总的局域势：LELF，LVTOT –结构优化参数：ISIF –等等。 主要参数说明如下： ? SYSTEM：该输入文件所要执行的任务的名字。取值：字符串，缺省值：SYSTEM

结构设计原理练习题C

结构设计原理练习题C 一、单项选择： 1、下列破坏形态中属于延性破坏的是： （ ） A ：超筋梁的破坏 B ：剪压破坏 C ：适筋梁的破坏 D:小偏心受压破坏 2、轴心受压柱中箍筋的主要作用是： （ ） A ：抗压 B ：约束钢筋不屈曲 C ：抗剪 D ： 防裂 3、同截面尺寸、同种材料的梁，只是钢筋用量不同，则承载能力关系：（ ） A ：超筋梁>适筋梁>少筋梁 B ：适筋梁>超筋梁>少筋梁 C ：少筋梁>适筋梁>超筋梁 D ：超筋梁>少筋梁>适筋梁 4、螺旋式间接钢筋的体积配筋率为 （ ） 0: s A A bh :sv v A B bs 11111112:s s n l A n l A C l l s + 14:ss cor A D d s 5、截面尺寸满足抗剪上限要求则不会发生： （ ） A ：剪压破坏 B ：斜拉破坏 C ：斜压破坏 D ：少筋破坏 6、先张法特有的应力损失是 （ ） A ：钢筋与孔道摩擦引起的应力损失 B ：台座与钢筋温差引起的应力损失 C ：钢筋松弛引起的应力损失 D ：混凝土收缩引起的应力损失 7、部分应力构件的预应力度： （ ） A ：0=λ B ：0λ< C ：1>λ D ：10<<λ 二、填空 1、构件按受力特点分 、 、 、受扭构件。 2、混凝土的强度设计值是由强度标准值 而得。 3、结构能满足各项功能要求而良好的工作叫 ，否则叫 。 4、 < f sd A s 时定义为第二类T 梁。 5、由于某种原因引起预应力钢筋的应力减小叫 。 三、判断正误 1、剪压破坏是延性破坏而斜拉破坏是脆性破坏。 （ ） 2、在轴心受压件中混凝土的收缩和徐变都会引起钢筋的压应力增长。 （ ） 3、ηe 0 >0.3h 0 时为大偏心受压。 （ ） 4、施加预应力不能提高构件的承载能力。 （ ） 5、局部承压面下混凝土的抗压强度比全截面受压时高。 （ ） 四、简答 1、什么叫开裂截面的换算截面？为什么使用换算截面？画矩形截面全截面换算截面的示意图。 2、 钢筋和混凝土之间的粘结力来源于哪几方面？ 3 、简述后张法施工过程？它有哪些优、缺点？ 4、简述等高度梁只设箍筋时的抗剪钢筋设计步骤。 五、计算题 1、T 形截面尺寸' '1200,200, 120,1000,f f b mm b mm h mm h mm ====采用C30混凝 土（MPa f cd 8.13=），HRB335级钢筋（MPa f sd 280=），Ⅰ类环境条件，56.0=b ξ，

结构设计原理计算方法

结构设计原理案例计算步骤 一、单筋矩形截面受弯构件正截面承载力计算 计算公式： ——水平力平衡 （）——所有力对受拉钢筋合力作用点取矩（） （）——所有力对受压区砼合力作用点取矩（）使用条件： 注：/，&& 计算方法： ㈠截面设计yy 1、已知弯矩组合设计值，钢筋、混凝土强度等级及截面尺寸b、h，计算。 ①由已知查表得：、、、； ②假设； ③根据假设计算； ④计算（力矩平衡公式：）； ⑤判断适用条件：（若，则为超筋梁，应修改截面尺寸或提 高砼等级或改为双筋截面）； ⑥计算钢筋面积（力平衡公式：）； ⑦选择钢筋，并布置钢筋（若 ，则按一排布置）； 侧外 ⑧根据以上计算确定（若与假定值接近，则计算，否则以的确定值作 为假定值从③开始重新计算）； ⑨以的确定值计算； ⑩验证配筋率是否满足要求（，）。 2、已知弯矩组合设计值，材料规格，设计截面尺寸、和钢筋截面面积。 ①有已知条件查表得：、、、； ②假设，先确定； ③假设配筋率（矩形梁，板）； ④计算（，若，则取）； ⑤计算（令，代入）； ⑥计算（，&&取其整、模数化）； ⑦确定（依构造要求，调整）； ⑧之后按“1”的计算步骤计算。 ㈡承载力复核 已知截面尺寸b、，钢筋截面面积，材料规格，弯矩组合设计值，

所要求的是截面所能承受的最大弯矩，并判断是否安全。 ①由已知查表得：、、、； ②确定； ③计算； ④计算（应用力平衡公式：，若，则需调整。令， 计算出，再代回校核）； ⑤适用条件判断（，，）； ⑥计算最大弯矩（若，则按式计算最大弯矩） ⑦判断结构安全性（若，则结构安全，但若破坏则破坏受压区，所以应以受压区控制设计；若,则说明结构不安全，需进行调整——修改尺寸或提高砼等级或改为双筋截面）。 二、双筋矩形截面梁承载力计算 计算公式： ， ，（）+（） 适用条件： （1） （2） 注：对适用条件的讨论 ①当&&时，则应增大截面尺寸或提高砼等级或增加的用量（即 将当作未知数重新计算一个较大的）；当时，算得的即为安全要 求的最小值，且可以有效地发挥砼的抗压强度，比较经济； ②当&&时，表明受压区钢筋之布置靠近中性轴，梁破坏时应变较 小，抗压钢筋达不到其设计值，处理方法： a.《公桥规》规定：假定受压区混凝土压应力的合力作用点与受压区钢筋合力作用 点重合，并对其取矩，即 令2，并 （） 计算出； b.再按不考虑受压区钢筋的存在（即令），按单筋截面梁计算出。 将a、b中计算出的进行比较，若是截面设计计算则取其较小值，若是承载能力复核则取其较大值。 计算方法： ㈠截面设计 1.已知截面尺寸b、h，钢筋、混凝土的强度等级，桥梁结构重要性系数，弯矩组合 设计值，计算和。 步骤： ①根据已知查表得：、、、、； ②假设、（一般按双排布置取假设值）； ③计算；

混凝土结构设计原理试卷之计算题题库 ()

1、某现浇多层钢筋混凝土框架结构，地层中柱按轴心受压构件计算，柱高H=6.4m ，承受轴向压力设计值N=2450kN,采用C30级混凝土，HRB335级钢筋，求柱截面尺寸（设配筋率 '0.01,1ρ?==），并试计算需配置的纵向受力钢筋。 （已知：2 14.3N/mm c f =，21.43/t f N mm =，'2300/y y f f N mm ==） 附表：钢筋混凝土轴心受压构件的稳定系数? 设配筋率' 0.01,1ρ?==，由公式知 正方形截面边长396.7b mm ==，取b=400mm 。 （2）求稳定系数 柱计算长度0 1.0l H =， 06400 16400 l b ==，查表得0.87?=。 （3）计算配筋 由公式知 2、某梁截面尺寸b×h=250mm×500mm ，M=2.0×108N·mm ，受压区预先已经配好HRB335级受压钢筋2φ20（' s A =628mm 2 ），若受拉钢筋也采用HRB335级钢筋配筋，混凝土的强度等级为C30，求截面所需配置的受拉钢筋截面面积s A 。 （已知：2 14.3N/mm c f =，21.43/t f N mm =，'2300/y y f f N mm ==，1 1.0α=， ,max 0.55,0.399b s ξα==） 解：(1)求受压区高度x 假定受拉钢筋和受压钢筋按一排布置，则' 35mm s s a a == 且' 2235mm 70mm s x a >=?= (2)计算截面需配置的受拉钢筋截面面积 四、计算题 1、已知某屋架下弦，截面尺寸b=220mm ，h=150mm ，承受轴心拉力设计值N=240kN ，混凝土为C30级，纵筋为HRB335级，试计算需配置的纵向受力钢筋。 （已知：2 14.3N/mm c f =，21.43/t f N mm =，'2300/y y f f N mm ==） 参考答案： 解：,u N N =令 2、已知梁的截面尺寸b=250mm ，h=500mm ，混凝土为C30级，采用HRB400级钢筋，承

《结构设计原理》述课

《结构设计原理》述课 一、前言 （一）课程基本信息 1.课程名称：结构设计原理 2.课程类别：专业平台课 3.学时：两学期总计84学时，2周课程设计 4.适用专业：交通工程 （二）课程性质 1.课程性质 结构是土木工程中最基本的元素，《结构设计原理》课程围绕着工程中常用的钢筋混凝土结构、预应力混凝土结构、圬工结构的设计计算进行理论和实践性的教学。 《结构设计原理》是土木工程专业的一门重要的专业必修课程，是学生运用已学的《工程制图》、《理论力学》、《材料力学》、《结构力学》、《工程材料》等知识，初步解决结构原理及结构设计问题的一门课程。其特点是：兼具理论性和实用性且承前启后，为学好专业课打好基础的课程，也是学生感到比较难学的一门课程。所以《结构设计原理》及其系列课程一直是土木工程专业的主干课，从开设的《结构设计原理》、《结构设计原理》课程设计，到毕业设计都渗透结构设计的理论，课程贯穿交通工程专业教学的所有环节。 本课程主要介绍钢筋混凝土结构、预应力混凝土结构和圬工结构的各种基本构件受力特性、设计原理、计算方法和构造设计。 2.本课程的作用 本课程主要培养学生掌握钢筋混凝土基本构件和结构的设计计算方法和与施工及工程质量有关的结构的基本知识，培养学生具有识读桥梁结构图纸的识读能力、基本构件的设计能力、使用和理解各种结构设计规范能力、解决工程结构实际问题的能力、综合分析问题的能力、学习能力和与人合作等能力，从而为继续学习后续专业课程奠定扎实的基础，以进一步培养学生树立独立思考、吃苦耐劳、勤奋工作的意识以及诚实、守信的优秀品质，为今后从事施工生产一线的工作奠定良好的基础。 本课程以“材料力学”、“理论力学”和“工程材料”的学习为基础共同打造学生的专业核心技能。

初学VASP中电子态密度计算设置参考

初学VASP中电子态密度计算基本设置参考主要分成三步：一、结构优化；二、静态自洽计算；三、非自洽计算以Al－FCC为例子 第一步结构优化 输入文件（INCAR, POTCAR, POSCAR, KPOINT） INCAR文件 System=Al ISTART=0 ISMEAR=1 SIGMA=0.2 ISPIN=2 GGA=91; VOSKOWN=1; EDIFF=0.1E-05; EDIFFG=-0.01 IBRION=2 NSW=50 ISIF=2 (OR 3) NPAR=10 POTCAR 文件直接在势库中拷贝 POSCAR文件 Al 4.05 1.0 0.0 0.0 0.0 1.0 0.0

0.0 0.0 1.0 4 Direct 0.0 0.0 0.0 0.5 0.5 0.0 0.5 0.0 0.5 0.0 0.5 0.5 KPOINT 文件 Automatic generation Mohkorst Pack 15 15 15 0.0 0.0 0.0 第二步静态自洽计算 INCAR：PREC = Medium，ISTART = 0，ICHARG = 2，ISMEAR = -5输入文件（INCAR, POTCAR, POSCAR, KPOINT） INCAR文件 System=Al ISTART=0 ISMEAR=1 SIGMA=0.2 ISPIN=2

GGA=91; VOSKOWN=1; EDIFF=0.1E-05; EDIFFG=-0.01 #IBRION=2 #NSW=50 #ISIF=2 (OR 3) NPAR=10 POTCAR 文件直接在势库中拷贝 POSCAR文件 Al 4.05 1.0 0.0 0.0 0.0 1.0 0.0 0.0 0.0 1.0 4 Selective Dynamic Direct 0.0 0.0 0.0 T T T 0.5 0.5 0.0 T T T 0.5 0.0 0.5 T T T 0.0 0.5 0.5 T T T KPOINT 文件 Automatic generation

用vasp计算硅的能带结构

用vasp计算硅的能带结构 在最此次仿真之前，因为从未用过vasp软件，所以必须得学习此软件及一些能带的知识。vasp是使用赝势和平面波基组，进行从头量子力学分子动力学计算的软件包。用vasp计算硅的能带结构首先要了解晶体硅的结构，它是两个嵌套在一起的FCC布拉菲晶格，相对的位置为 (a/4,a/4,a/4), 其中a=5.4A是大的正方晶格的晶格常数。在计算中，我们采用FCC的原胞，每个原胞里有两个硅原子。 VASP计算需要以下的四个文件：INCAR（控制参数）, KPOINTS（倒空间撒点）, POSCAR（原子坐标）, POTCAR（赝势文件） 为了计算能带结构，我们首先要进行一次自洽计算，得到体系正确的基态电子密度。然后固定此电荷分布，对于选定的特殊的K点进一步进行非自洽的能带计算。有了需要的K点的能量本征值，也就得到了我们所需要的能带。 步骤一.—自洽计算产生正确的基态电子密度： 以下是用到的各个文件样本: INCAR 文件： SYSTEM = Si Startparameter for this run: NWRITE = 2; LPETIM=F write-flag & timer PREC = medium medium, high low ISTART = 0 job : 0-new 1-cont 2-samecut ICHARG = 2 charge: 1-file 2-atom 10-const ISPIN = 1 spin polarized calculation? Electronic Relaxation 1 NELM = 90; NELMIN= 8; NELMDL= 10 # of ELM steps EDIFF = 0.1E-03 stopping-criterion for ELM LREAL = .FALSE. real-space projection Ionic relaxation EDIFFG = 0.1E-02 stopping-criterion for IOM NSW = 0 number of steps for IOM IBRION = 2 ionic relax: 0-MD 1-quasi-New 2-CG ISIF = 2 stress and relaxation POTIM = 0.10 time-step for ionic-motion TEIN = 0.0 initial temperature TEBEG = 0.0; TEEND = 0.0 temperature during run

VASP遇到小总结问题

VASP 计算的过程遇到的问题 01、第一原理计算的一些心得 （1）第一性原理其实是包括基于密度泛函的从头算和基于Hartree-Fock自洽计算的从头算，前者以电子密度作为基本变量（霍亨伯格-科洪定理），通过求解Kohn-Sham方程，迭代自洽得到体系的基态电子密度，然后求体系的基态性质；后者则通过自洽求解Hartree-Fock方程，获得体系的波函数，求基态性质； 评述：K-S方程的计算水平达到了H-F水平，同时还考虑了电子间的交换关联作用。 （2）关于DFT中密度泛函的Functional，其实是交换关联泛函 包括LDA，GGA，杂化泛函等等 一般LDA为局域密度近似，在空间某点用均匀电子气密度作为交换关联泛函的唯一变量，多数为参数化的CA-PZ方案； GGA为广义梯度近似，不仅将电子密度作为交换关联泛函的变量，也考虑了密度的梯度为变量，包括PBE,PW,RPBE等方案，BL YP泛函也属于GGA； 此外还有一些杂化泛函，B3L YP等。 （3）关于赝势 在处理计算体系中原子的电子态时，有两种方法，一种是考虑所有电子，叫做全电子法，比如WIEN2K中的FLAPW方法(线性缀加平面波)；此外还有一种方法是只考虑价电子，而把芯电子和原子核构成离子实放在一起考虑，即赝势法，一般赝势法是选取一个截断半径，截断半径以内，波函数变化较平滑，和真实的不同，截断半径以外则和真实情况相同，而且赝势法得到的能量本征值和全电子法应该相同。 赝势包括模守恒和超软，模守恒较硬，一般需要较大的截断能，超软势则可以用较小的截断能即可。另外，模守恒势的散射特性和全电子相同，因此一般红外，拉曼等光谱的计算需要用模守恒势。 赝势的测试标准应是赝势与全电子法计算结果的匹配度，而不是赝势与实验结果的匹配度，因为和实验结果的匹配可能是偶然的。 （4）关于收敛测试 （a）Ecut，也就是截断能，一般情况下，总能相对于不同Ecut做计算，当Ecut增大时总能变化不明显了即可；然而，在需要考虑体系应力时，还需对应力进行收敛测试，而且应力相对于Ecut的收敛要比总能更为苛刻，也就是某个截断能下总能已经收敛了，但应力未必收敛。 （b）K-point，即K网格，一般金属需要较大的K网格，采用超晶胞时可以选用相对较小的K网格，但实际上还是要经过测试。 （5）关于磁性 一般何时考虑自旋呢？举例子，例如BaTiO3中，Ba、Ti和O分别为+2，+4和-2价，离子全部为各个轨道满壳层的结构，就不必考虑自旋了；对于BaMnO3中，由于Mn+3价时d 轨道还有电子，但未满，因此需考虑Mn的自旋，至于Ba和O则不必考虑。其实设定自旋就是给定一个原子磁矩的初始值，只在刚开始计算时作为初始值使用，具体的可参照磁性物理。 （6）关于几何优化 包括很多种了，比如晶格常数和原子位置同时优化，只优化原子位置，只优化晶格常数，还有晶格常数和原子位置分开优化等等。

VASP几个计算实例

用VASP计算Ｈ原子的能量 氢原子的能量为。在这一节中，我们用VASP计算H原子的能量。对于原子计算，我们可以采用如下的INCAR文件 PREC=ACCURATE NELMDL=5make five delays till charge mixing ISMEAR=0;SIGMA=0.05use smearing method 采用如下的KPOINTS文件。由于增加K点的数目只能改进描述原子间的相互作用，而在单原子计算中并不需要。所以我们只需要一个K点。 Monkhorst Pack0Monkhorst Pack 111 000 采用如下的POSCAR文件 atom1 15.00000.00000.00000 .0000015.00000.00000 .00000.0000015.00000 1 cart 000 采用标准的H的POTCAR 得到结果如下： k-point1:0.00000.00000.0000 band No.band energies occupation 1-6.3145 1.00000 2-0.05270.00000 30.48290.00000 40.48290.00000 我们可以看到，电子的能级不为。 Free energy of the ion-electron system(eV) --------------------------------------------------- alpha Z PSCENC=0.00060791 Ewald energy TEWEN=-1.36188267 -1/2Hartree DENC=-6.27429270 -V(xc)+E(xc)XCENC= 1.90099128 PAW double counting=0.000000000.00000000 entropy T*S EENTRO=-0.02820948 eigenvalues EBANDS=-6.31447362 atomic energy EATOM=12.04670449 ---------------------------------------------------