pwscf计算弹性系数

PWSCF 程序包（早期的叫法），或称为ESPRESSO程序（改名后的叫法)，它包括了多几个计算模块，主要的是电子自洽计算模块pw.x，晶格动力学计算模块 (ph.x, phcg.x, dynmat.x,d3.x等)，后续数据处理模块pp.x，电子输运性质计算模块pwcond.x，分子动力学模块cp.x等 一、自洽计算

例子：fcc Cu的自洽计算 &control

calculation='scf'

restart_mode='from_scratch', pseudo_dir = './', outdir='./' prefix='cu' tstress = .true. tprnfor = .true. / &system

ibrav = 2, celldm(1) =6.73, nat= 1, ntyp= 1, ecutwfc = 25.0, ecutrho = 300.0

occupations='smearing', smearing='gaussian', degauss=0.02 /

&electrons

diagonalization='david' conv_thr = 1.0e-8 mi

PWSCF 程序包（早期的叫法），或称为ESPRESSO程序（改名后的叫法)，它包括了多几个计算模块，主要的是电子自洽计算模块pw.x，晶格动力学计算模块 (ph.x, phcg.x, dynmat.x,d3.x等)，后续数据处理模块pp.x，电子输运性质计算模块pwcond.x，分子动力学模块cp.x等 一、自洽计算

例子：fcc Cu的自洽计算 &control

calculation='scf'

restart_mode='from_scratch', pseudo_dir = './', outdir='./' prefix='cu' tstress = .true. tprnfor = .true. / &system

ibrav = 2, celldm(1) =6.73, nat= 1, ntyp= 1, ecutwfc = 25.0, ecutrho = 300.0

occupations='smearing', smearing='gaussian', degauss=0.02 /

&electrons

diagonalization='david' conv_thr = 1.0e-8 mi

湖南工业大学

课 程 设 计

资 料 袋

理 学院（系、部） 2011 ~ 2012 学年第 一 学期 课程名称 计算材料学 指导教师 雷军辉 职称 讲师

学生姓名 余晓燕 专业班级 应用物理081班 学号 08411200135

题 目 计算BN的弹性常数

成 绩 起止日期 2011年 12月 4日 ～ 2011年 12 月 12 日

目 录 清 单

序号 1 2 3 4 5 6

湖南工业大学

材 料 名 称 课程设计任务书 课程设计说明书 课程设计图 资料数量 备 注 张

1

课程设计任务书

2011—2012 学年第 1 学期

理学院 学院（系、部） 应用物理学 专业 081 班级

课 程 名 称：

土工试验-渗透试验指导书内容介绍 第一节 概 述

渗透是水在多孔介质中运动的现象。土体为多孔介质，水能够在土的孔隙中流动的特性叫土的渗透性。若渗透水流在土中呈层流状态流动，则满足达西定律，即渗汉速度v与水力梯度I成正比，表达式为： v=k?i

式中：k——渗透系数，cm/s； v=q/A,渗流速度，cm/s； i=h/L,水力梯度；

q——单位时间渗流量，cm3/s；

A——垂直渗流方向的横截面积cm2； h——水位差，cm； L——渗径长度，cm。

渗透试验的目的就是测量渗透系数k 。

由土的渗透系数评价土的渗透性大小，计算水工建筑物渗流量等。

第二节 试验原理

渗透试验原理就是在试验装置中测出渗流量，不同点的水头高度，从而计算出渗流速度和水力梯度，代入公式（8-1）计算出渗透系数。由于土的渗透系数变化范围很大，自大于10-1cm/s到小于10-7cm/s。实验室内常用两种不同的试验装置进行试验：常水头试验装置用来测定渗透系数k比较大的无凝聚性土的渗透系数；变水头渗透试验装置用来测定渗透系数K较小的凝聚性土的渗透系数。特殊设计的变水头试验测定粗粒土渗透系数和常水头试验测定渗透性极小的粘性土渗透系数也很常用。

7.2.2扩散系数

费克定律中的扩散系数Ｄ代表单位浓度梯度下的扩散通量，它表达某个组分在介质中扩散的快慢，是物质的一种传递性质。 一、气体中的扩散系数

气体中的扩散系数与系统、温度和压力有关，其量级为10m/s。通常对于二元气体Ａ、B的相互扩散，Ａ在B中的扩散系数和B在A中的扩散系数相等，因此可略去下标而用同一符号Ｄ表示，即DAB?DBA?D。

表７－１给出了某些二元气体在常压下（1.013?10Pa）的扩散系数。

对于二元气体扩散系数的估算，通常用较简单的由富勒（Fuller）等提出的公式：

5?520.0101T1.75D?11?MAMB （７－１９）

P[(?vA)1/3?(?vB)1/3]22式中，D－Ａ、B二元气体的扩散系数，m/s；

P－气体的总压，Pa； T－气体的温度，Ｋ；

MA、MB－组分Ａ、B的摩尔质量，kg/kmol；

、－组分Ａ、B分子扩散体积，cm/mol。

一般有机化合物可按分子式由表７－２查相应的原子扩散体积加和得到，某些简单物质则在表７－２种直接列出。 系统 H2－空气 He－空气 O2－空气 Cl2－空气 H2O－空气 NH3－空气 CO2－空气 SO2－空气 5表7－１ 某些二元气

7.2.2扩散系数

费克定律中的扩散系数Ｄ代表单位浓度梯度下的扩散通量，它表达某个组分在介质中扩散的快慢，是物质的一种传递性质。 一、气体中的扩散系数

气体中的扩散系数与系统、温度和压力有关，其量级为10m/s。通常对于二元气体Ａ、B的相互扩散，Ａ在B中的扩散系数和B在A中的扩散系数相等，因此可略去下标而用同一符号Ｄ表示，即DAB?DBA?D。

表７－１给出了某些二元气体在常压下（1.013?10Pa）的扩散系数。

对于二元气体扩散系数的估算，通常用较简单的由富勒（Fuller）等提出的公式：

5?520.0101T1.75D?11?MAMB （７－１９）

P[(?vA)1/3?(?vB)1/3]22式中，D－Ａ、B二元气体的扩散系数，m/s；

P－气体的总压，Pa； T－气体的温度，Ｋ；

MA、MB－组分Ａ、B的摩尔质量，kg/kmol；

、－组分Ａ、B分子扩散体积，cm/mol。

一般有机化合物可按分子式由表７－２查相应的原子扩散体积加和得到，某些简单物质则在表７－２种直接列出。 系统 H2－空气 He－空气 O2－空气 Cl2－空气 H2O－空气 NH3－空气 CO2－空气 SO2－空气 5表7－１ 某些二元气

地震影响系数计算

地震影响系数计算作者： 徐 翔(高强螺栓) 单位： 西安建筑科技大学 电子信箱： xuxiang@

0.200 0.100 0.000

α

T

地震影响系数曲线

八度多遇地震； Ⅲ类场地土； 第一组 设计基本地震加速度值= 0.2 阻尼比 ζ = 0.05 水平地震影响系数最大值α max= 0.16 特征周期Tg= 曲线下降段的衰减指数γ = 直线下降段的下降斜率调整系数η 1= 阻尼调整系数η 2= 0.45 0.9 0.02 1

3.2.2 5.1.5砼 8.2.2钢 表5.1.4-1 表5.1.4-2

说明：1.《建筑抗震设计规范GB50011-2001》第5.1.5条地震影响系数曲线。2.红色为输入值；将A1至D46格粘贴到".txt"文件。 3.“八度多遇地震、Ⅲ类场地土、第一组、设计基本地震加速度值=0.2g”宜用文件名"8d3102g.txt"。 4.在SAP2000的'Response Spectrum function Definition' 对话框中，‘Number of Points per Line'为24。 5.在ETABS的'Response Spectrum function Definition' 对话框中

三．试确定以下两组应变状态能否存在(K,A,B为常数), 并说明为什么？

(1) (2)

?x?K(x2?y2),?y?Ky2,?xy?2Kxy (存在) ?x?Axy2,?y?Bx2y,?xy?0 (不存在)

四．计算题

1. 图中所示的矩形截面体，受力如图所示，试写出其边界条件。

解：主要边界条件，

x?b，?x?0;?xy?p

x??b，?x?q;?xy?0

次要边界条件，在y?0上，

(?xy)y?0?0，满足；

?

bb?b(?y)y?0dx??F

Fb ??b22.图中所示的矩形截面体，在o处受有集中力F和力矩M?Fb/2作用，试用应力函数??Ax3?Bx2求解图示问题的应力分量，设在A点的位移和转角均为零。

(?y)y?0xdx??

解：应用应力函数求解，

4

（1） 校核相容方程???0，满足 （2） 求应力分量，在无体力时，得

?y?6Ax?2B，?y

建筑体形系数计算书

建筑体形系数计算书

1． 计算公式：

Vo=F1·h1+ F2·h2+ (R·h3)/3 Fo=E1+E2+E3+E4+R1 R1=R·1.2

Vo 建筑物的体积 Fo 建筑物的外表面积

F1 建筑首层面积 F2 建筑二层面积

R 建筑屋顶面积

R1 建筑物屋顶总表面积

E1- E4 建筑物的外立面面积

h1 建筑首层层高 h2 建筑二层层高 h3 建筑屋面高

2． 体形系数计算:

F1 =182.38-8.67=173.71m2

F2 =171.00-34.49-4.67=131.84m2

E1=84.39m2 E2=91.02m2 E3=80.39m2 E4=97.08m2

R=210.15m2 R1=210.15·1.2=252.18m2

h1=3.35m h2=2.74m h3=2.8m

Vo=173.71·3.35+131.84·2.74+210.15·2.8/3 =581.9285+361.2416+196.14 =1139.31m3

Fo=84.39+91.02+80.39+97.

横向分布系数的示例计算

一座五梁式装配式钢筋混凝土简支梁桥的主梁和横隔梁截面如图，计算跨径L=19.5m，主梁翼缘板刚性连接。求各主梁对于车辆荷载和人群荷载的分布系数？

杠杆原理法：

解：1绘制1、2、3号梁的荷载横向影响线如图所示

2再根据《公路桥涵设计通用规范》（JTG D60-2004） 规定，在横向影响线上确定荷载沿横向最不利布置位置。 如图所示： 对于1号梁： 车辆荷载：mcq?沥青混凝土厚2cmC25混凝土垫层6-12cm11???2?0.967?0.484 212345人群荷载：mcr??r?1.417 对于2号梁： 车辆荷载：mcq?1.41710.9671号梁人群荷载：mcr??r?0.417 对于3号梁： 车辆荷载：mcq?0.41711???2?1?0.5 22号梁11???2?1?0.5 2113号梁人群荷载：mcr??r?0

4、5号梁与2、1号梁对称，故荷载的横向分布系数相同。

偏心压力法

（一）假设：荷载位于1号梁 1长宽比为lb?19.5?2.6?2，故可按偏心压

5?1.5力法来绘制横向影响线并计算横向分布系数mc。 本桥的各根主梁的横截面积均相等，梁数为5，梁的间距为1.5m，则：

0.6330.6000.5