Ansys 有限元分析
一维框架（线模型）
一、CAD画图，以iges 格式输出。
二、导入iges 文件
三、输入Ansys 命令流
/ PREP7 !进入前处理
ET,1,beam3 !选择单元类型
R,1,6.5e-7,6.8e-4 !给出实常数(横截面积、惯性矩)
MP,EX,1,3e11 !给出材料的弹性模量
划分单元
LSEL,ALL !选择所有的线
LESIZE,all, , ,10, , , , ,1 !将所选择的线划分成10 段
lmesh,all
手动加载
d,1,all !1号节点施加全约束(X,Y方向都有约束不包含力偶) F,2,FY,-3000 !在2 号节点处施加Y 方向的力(-3000)
FINISH !结束前处理状态
/SOLU !进入求解模块
SOLVE !求解
FINISH !结束求解状态
/POST1 !进入后处理
PLDISP,1 !显示变形状况
FINISH !结束后处理
二维平面（面模型）
一、CAD画图，以iges 格式输出。
二、导入iges 文件
三、输入Ansys 命令流
/PREP7 !进入前处理模块
/TITLE, NMEHZ
ANTYPE,STATIC !设置静态分析
ET,1, PLANE182!设置1号单元(用于平面)
MP,EX,1,3E7 !设置1号材料的弹模
MP,NUXY,1,0.3 !设置1号材料的泊松比
MP,DENS,1,7800 !设置1号材料的密度
MP,ALPX,1,1E-5 !设置1号材料的热膨胀率
R,1,6.5e-7,6.8e-4 !给出实常数(横截面积、惯性矩) 划分单元
LSEL,ALL !选择所有的线
LESIZE,all, , ,10, , , , ,1 !将所选择的线划分成10 段
! MSHAPE,1,2D !设置三角形单元
MSHAPE,0,2D !设置四边形单元
MSHKEY,1 !设置映射划分
AMESH,1 !对面No.1 进行网格划分ALLSEL,ALL !选择所有的对象
手动加载
ACEL,0,9.8,0, !在Y方向施加重力加速度
Allsel !选择所有内容
Outres,all,all !输出所有内容
FINISH !结束前处理状态
/SOLU !进入求解模块
SOLVE !求解
FINISH !结束求解状态
SAVE
/POST1 !进入后处理器PLDISP,1 !绘制变形和未变形图
PLNSOL,U,X,0,1 !绘制X方向位移图PLNSOL,S,X,0,1 !绘制X方向应力图PLNSOL,U,Y,0,1 !绘制Y方向位移图PLNSOL,S,Y,0,1 !绘制Y方向应力图
PLNSOL,S,1,0,1 !绘制第一主应力图PLNSOL,S,3,0,1 !绘制第三主应力图plnsol,s,eqv,0,1 !显示等效应力云图PLLS, ISHEAR, JSHEAR ! 结构剪力分布图
PLLS, IMOMENT, JMOMENT ! 结构弯矩分布图
将画图区域背景设置为白色
/RGB,INDEX,100,100,100, 0
/RGB,INDEX, 80, 80, 80,13
/RGB,INDEX, 60, 60, 60,14
/RGB,INDEX, 0, 0, 0,15
/REPLOT
主窗口的背景色还原到黑色：
/RGB,INDEX, 0, 0, 0, 0
/RGB,INDEX, 60, 60, 60,13
/RGB,INDEX, 80, 80, 80,14
/RGB,INDEX,100,100,100,15
/REPLOT
保存图片
/image,save,'d:\ansys\gravity.bmp'
文本文件的保存
添加下列命令：
/out,pinlv,txt !保存到非工作目录，改为：/output,'d:\ansys\pinlv.txt'
Ansys 有限元分析
一、CAD画图，以iges 格式输出。
二、导入iges 文件
三、输入Ansys 命令流
菜单过滤设置
/NOPR !菜单过滤设置
KEYW，PR_STRUC，1 !保留结构分析部分菜单
进入前处理模块
/PREP7 !进入前处理模块
/TITLE, nmehz !在图形菜单中显示的文字
ANTYPE,STATIC !设置静态分析
设置单元类型
ET,1, BEAM4!设置1号单元
ET,1,PLANE55 !选取单元类型1(平面传热单元)
设置钢材的材料参数
MP,EX,1,3E7 !设置1号材料的弹模
MP,NUXY,1,0.3 !设置1号材料的泊松比
MP,DENS,1,7800 !设置1号材料的密度
MP,ALPX,1,1E-5 !设置1号材料的热膨胀率
或（设置3号单元）
MP,EX,3,2.1E11 $MP,PRXY,3,0.3 $MP,DENS,3,7800 !设置钢材的材料参数
MP,ALPX,3,1E-5 ! 设置钢材的热膨胀率
MP,KXX,1,0.025 !定义1号材料的热传导系数(KXX=0.025) MP,C,1,0.28 !定义1号材料的比热(0.28)
CSYS,4
k, 1 ,0,0,.000,
k, 2 ,30,0,.000,
k, 3 ,15,50,.000,
k, 4 ,0,50,.000,
l, 1 ,2
l, 2 ,3
l, 3 ,4
l, 1 ,4
命令：CSYS,KCN
其中KCN表示坐标系号码，0-直角坐标系（缺省），1-柱坐标系，2-球坐标系，4-以工作平面为坐标系，5-柱坐标系（以Y轴为转轴），≥11-局部坐标系。
由于工作平面可不断移动和旋转，因此当采用CSYS,4时也相当于不断定义了局部直角坐标，在很多情况下应用非常方便。
amesh,1 ! 对面1划分单元
ESIZE,,10 !设置线段均分为10段
MSHK,1 $MSHA,0,2D $AMESH,1 !使用映射网格划分, 四边形单元, 对面1划分单元
单元划分
vsel,s,,,5 !选择5号体
smrt,off !关闭程序智能划分功能,即接下来通过人为设置单元尺寸
LESIZE,53,,,40 !53号线划分成40份
lesize,58,,,40 !58号线划分成40份
ESIZE,,4 !指定单元份数
MSHK,1 !似乎应该是MSHKEY,1,使用映射网格
MSHA,0,3D !指定单元为三维六面体形状
VMESH,all !开始划分
smrt,1 !打开自动划分
vsel,s,,,6,7,1 !选择6和7号体
MSHAPE,1,3D !指定单元为三维四面体形状
MSHKEY,0 !自由划分
vmesh,all !开始划分
加载
1.施加重力加速度
ACEL,0,9.8,0, !在Y方向施加重力加速度
2.添加位移
D,3,ALL !将3 号节点的位移全部固定
d,1,ux !1号节点施加x 方向约束
d,2,uy !2号节点施加y 方向约束
d,3,all !3号节点施加全约束(X,Y方向都有约束不包含力偶)
3.施加集中力
F,1,FX,3000 !在1 号节点处施加x 方向的力(3000)
4.施加均布荷载
SFBEAM,1,1,PRESS,4167 !施加均布压力
求解
Allsel !选择所有内容
Outres,all,all !输出所有内容
FINISH !结束前处理状态
/SOLU !进入求解模块
SOLVE !求解
FINISH !结束求解状态
SAVE
进入后处理器
/POST1 !进入后处理器
PLDISP,1 !绘制变形和未变形图
PLNSOL,U,X,0,1 !绘制X方向位移图
PLNSOL,S,X,0,1 !绘制X方向应力图
PLNSOL,U,Y,0,1 !绘制Y方向位移图
PLNSOL,S,Y,0,1 !绘制Y方向应力图
PLNSOL,S,1,0,1 !绘制第一主应力图
PLNSOL,S,3,0,1 !绘制第三主应力图
PLLS, ISHEAR, JSHEAR ! 结构剪力分布图
PLLS, IMOMENT, JMOMENT ! 结构弯矩分布图
ANSYS中弯矩、剪力图的绘制
GUI:
General Postproc-plot Result-Contour Plot-Line Element Result
弹出画单元结果的对话框，分别在Labi和Labj依次选取SMIS6和SMIS12（弯矩图）、SMIS1和SMIS7（轴力图）、SMIS2和SMIS8（剪力图）
/ PREP7 !进入前处理
ET,1,beam3 !选择单元类型
R,1,6.5e-7,6.8e-4 !给出实常数(横截面积、惯性矩)
MP,EX,1,3e11 !给出材料的弹性模量
N,1,0,0.96,0 !生成4 个节点,坐标(0,0.96,0)，以下类似
N,2,1.44,0.96,0
N,3,0,0,0
N,4,1.44,0,0
E,1,2 !生成单元(连接1 号节点和2 号节点) ，以下类似
E,1,3
E,2,4
D,3,ALL !将3 号节点的位移全部固定
D,4,ALL !将4 号节点的位移全部固定
F,1,FX,3000 !在 1 号节点处施加x 方向的力(3000) SFBEAM,1,1,PRESS,4167 !施加均布压力
FINISH !结束前处理状态
/SOLU !进入求解模块
SOLVE !求解
FINISH !结束求解状态
/POST1 !进入后处理
PLDISP,1 !显示变形状况
FINISH !结束后处理
/NOPR !菜单过滤设置
KEYW，PR_STRUC，1 !保留结构分析部分菜单
/PREP7 !进入前处理模块
/TITLE, NMEHZ
ANTYPE,STATIC !设置静态分析
ET,1, PLANE182!设置1号单元
R,1,6.5e-7,6.8e-4 !给出实常数(横截面积、惯性矩)
MP,EX,1,3E7 !设置1号材料的弹模
MP,NUXY,1,0.3 !设置1号材料的泊松比
MP,DENS,1,7800 !设置1号材料的密度
MP,ALPX,1,1E-5 !设置1号材料的热膨胀率
手动加载
ACEL,0,9.8,0, !在Y方向施加重力加速度
Allsel !选择所有内容
Outres,all,all !输出所有内容
FINISH !结束前处理状态
/SOLU !进入求解模块
SOLVE !求解
FINISH !结束求解状态
SAVE
/POST1 !进入后处理器
PLDISP,1 !绘制变形和未变形图
PLNSOL,U,X,0,1 !绘制X方向位移图
PLNSOL,S,X,0,1 !绘制X方向应力图
PLNSOL,U,Y,0,1 !绘制Y方向位移图
PLNSOL,S,Y,0,1 !绘制Y方向应力图
PLNSOL,S,1,0,1 !绘制第一主应力图
PLNSOL,S,3,0,1 !绘制第三主应力图
PLLS, ISHEAR, JSHEAR ! 结构剪力分布图
PLLS, IMOMENT, JMOMENT ! 结构弯矩分布图
CSYS,4
k, 1 ,0,0,.000,
k, 2 ,30,0,.000,
k, 3 ,15,50,.000,
k, 4 ,0,50,.000,
l, 1 ,2
l, 2 ,3
l, 3 ,4
l, 1 ,4
PLNSOL,S,X ! 在图形窗口显示节点上的X方向正应力
PLNSOL,S,XY ! 在图形窗口显示节点上的XY方向的剪应力
PLNSOL,S,Y ! 在图形窗口显示节点上的Y方向正应力
PLNSOL,S,1 ! 在图形窗口显示节点上的第一主应力(拉应力)
PLNSOL,S,3 ! 在图形窗口显示节点上的第三主应力(压应力)
PLNSOL,S,EQV ! 在图形窗口显示所有节点Mises应力
PRNSOL,S,PRIN !列表显示节点主应力(包括第一，二，三主应力和剪应力，切应力)
8.4 桥梁全桩基础受力分析
8.4.3 建模
2. 材料、实常数和单元类型定义
/TITLE，Mechanical analysis on railway tunnel 2nd lining ! 确定分析标题/NOPR !菜单过滤设置
KEYW，PR_STRUC，1 !保留结构分析部分菜单
/PREP7 !进入前处理器
ET,1,SOLID45 !选择8节点实体单元类型
MP,DENS,1,,2500 !输入密度
MP,EX,1,,30e9 !输入弹性模量
MP,PRXY,1,,0.2 !输入泊松比
SAVE !保存数据库
3. 建立几何模型
K,1 ,-4,0,-4, !创建承台关键点
K,2 ,-4,0,4,
K,3 ,4,0,4,
K,4 ,4,0,-4,
K,5 ,-4,1,-4,
K,6 ,-4,1,4,
K,7 ,4,1,4,
K,8 ,4,1,-4,
K,9 ,-2,1,-2, !创建桩位置关键点
K,10 ,-2,1,2,
K,11 ,2,1,2,
K, 12,2,1,-2,
V, 1, 2, 3, 4, 5, 6, 7, 8 !通过关键点创建实体，见图8-37。
wpoff,0,1,0 !移动局部坐标系
wprot,0,90,0 !旋转局部坐标系
CYL4,-2,-2,1, , , ,21 !创建四个桩圆柱体
CYL4,-2,2,1, , , ,21
CYL4,2,2,1, , , ,21
CYL4,2,-2,1, , , ,21
wprot,0,-90,0 !旋转局部坐标系
wpoff,0,-1,0 !移动局部坐标系
VOVLAP,1,2,3,4,5 !进行体布尔运算
SAVE !保存数据, 生成的实体几何模型如图8-38所示。
4. 单元网格划分
NUMMRG,ALL, , , ,LOW !合并所有元素
NUMCMP,ALL !压缩所有元素编号
LESIZE,要划分单元的线的编号, , ,单元划分的个数, , , , , ,
!采用以上命令设置单元大小, 其具体每条线要划分的单元数量如图8-39和图8-40所示。VSWEEP,1,,, !采用扫描方式对体进行单元划分
VSWEEP,2,,,
VSWEEP,3,,,
VSWEEP,4,,,
VSWEEP,5,,,
VSWEEP,6,,,
VSWEEP,7,,,
VSWEEP,8,,,
VSWEEP,9,,,
SAVE !保存数据, 生成的网格如图8-41所示。
Finish !返回Main Menu 主菜单
8.4.4 加载与求解
1. 加载
（1）对四根桩底部各节点施加“Ux”、“Uy”和“Uz”三个方向的约束。
/SOL !进入求解器
Allsel !选择所有内容
ASEL,S,,,6 !选择桩的底面
ASEL,S,,,8
ASEL,S,,,10
ASEL,S,,,12
NSLA,R,1 !再选择桩的底面上的节点
D,all, , , , , ,ALL, , , , , !在选择的所有节点上施加位移约束
（2）施加重力加速度。
ACEL,0,10,0, !在Y方向施加重力加速度
（3）施加力边界条件，加上荷载和位移边界条件后的几何模型见图8-42。
SFA,5,1,PRES,1875000 !在承台顶面施加面压力
SFA,7,1,PRES,1875000
SFA,9,1,PRES,1875000
SFA,11,1,PRES,1875000
SFA,34,1,PRES,1875000
2. 求解
（1）求解前设置。
Allsel !选择所有内容
Outres,all,all !输出所有内容
（2）求解。
Solve !求解计算
Finish !求解结束返回Main Menu 主菜单
SAVE
8.4.5 后处理
1.查看变形和位移结果
(1)绘制变形图，如图8-43所示。路径为：General Postproc> Plot Results > Deformed Shape > Def＋Undeformed。
/POST1 !进入后处理器
PLDISP,1 !绘制变形和未变形图
(2)绘制Y方向位移图，如图8-44所示。路径为：General Postproc> Plot Results > Contour Plot>Nodal Solu>DOF Solution>Translation UY。
PLNSOL,U,Y,0,1 !绘制Y方向位移图
PLNSOL,U,X,0,1 !绘制X方向位移图
2.查看应力结果
(1)绘制应力图，如图8-45所示。路径为：General Postproc> Plot Results > Contour Plot>Nodal Solu>Stress>SY。
PLNSOL,S,Y,0,1 !绘制Y方向应力图
PLNSOL,S,X,0,1 !绘制X方向应力图
(2)绘制第一主应力图，如图8-46所示。路径为：General Postproc> Plot Results > Contour Plot>Nodal Solu>Stress>S1。
PLNSOL,S,1,0,1 !绘制第一主应力图
(3)绘制第三主应力图，如图8-47所示。路径为：General Postproc> Plot Results > Contour Plot>Nodal Solu>Stress>S3。
PLNSOL,S,3,0,1 !绘制第三主应力图
2 完全的命令流
如图3-19 所示的框架结构，其顶端受均布力作用，用有限元方法分析该结构的位移。
结构中各个截面的参数都为：，，。11 3.0 10 PaE =× 74 6.5 10 mI ? =× 42 6.8 10 mA ? =× 在ANSYS 平台上，完成相应的力学分析。
/ PREP7 !进入前处理
ET,1,beam3 !选择单元类型
R,1,6.5e-7,6.8e-4 !给出实常数(横截面积、惯性矩)
MP,EX,1,3e11 !给出材料的弹性模量
N,1,0,0.96,0 !生成4 个节点,坐标(0,0.96,0)，以下类似
N,2,1.44,0.96,0
N,3,0,0,0
N,4,1.44,0,0
E,1,2 !生成单元(连接1 号节点和2 号节点) ，以下类似
E,1,3
E,2,4
D,3,ALL !将3 号节点的位移全部固定
D,4,ALL !将4 号节点的位移全部固定
F,1,FX,3000 !在 1 号节点处施加x 方向的力(3000) SFBEAM,1,1,PRESS,4167 !施加均布压力
FINISH !结束前处理状态
/SOLU !进入求解模块
SOLVE !求解
FINISH !结束求解状态
/POST1 !进入后处理
PLDISP,1 !显示变形状况FINISH !结束后处理
添加均布载荷
梁单元？solution>define loads>apply>structural>pressure>on beams 鼠标选取单元16 单击OK 弹出Apply PRES on Beams
Load key 输入1
Pressure value at node I 输入10
Pressure value at node J 输入10
单击apply
国内各地区重力加速度
序号地区力加速度地区修正值
g(m/s2) g/1kg g/3kg g/6kg g/15kg g/30kg
1 包头9.7986 -0.3981 -1.1943 -2.3886
-11.9430 -11.9430
2 北京9.8015 -0.7045 -2.1135 -4.2270
-10.5675 -21.1350
3 长春9.8048 -1.0413 -3.1239 -6.2478
-15.6195 -31.2390
4 长沙9.791
5 0.3267 0.9801 1.9602 9.8010 9.8010
5 成都9.7913 0.3267 0.9801 1.9602 4.9005 9.8010
6 重庆9.7914 0.326
7 0.9801 1.9602 4.9005 9.8010
7 大连9.8011 -0.6636 -1.9908 -3.9816
-9.9540 -19.9080
8 广州9.7833 0.6432 1.9296 3.8592 9.6480 19.2960
9 贵阳9.7968 0.7963 2.3889 4.7778 23.8890 23.8890
10 哈尔滨9.8066 -1.2251 -3.6753 -7.3506
-18.3765 -36.7530
11 杭州9.7936 0.1020 0.3060 0.6120 1.5300 3.0600
12 海口9.7863 0.8474 2.5422 5.0844 25.4220 25.4220
13 合肥9.7947 0.0204 0.0612 0.1224 0.3060 0.6120
14 吉林9.8048 -1.0413 -3.1239 -6.2478
-15.6195 -31.2390
15 济南9.7988 -0.3981 -1.1943 -2.3886
-5.9715 -11.9430
16 昆明9.7830 1.1230 3.3690 6.7380 16.8450 33.6900
18 南昌9.7920 0.2654 0.7962 1.5924 7.9620 7.9620
19 南京9.7949 -0.0306 -0.0918 -0.1836
-0.4590 0.9180
20 南宁9.7877 0.7044 2.1132 4.2264 10.5660 21.1320
21 青岛9.7985 -0.3981 -1.1943 -2.3886
-5.9715 -11.9430
22 上海9.7964 0.0000 0.0000 0.0000 0.0000 0.0000
23 沈阳9.8035 -0.9086 -2.7258 -5.4516
-13.6290 -27.2580
24 石家庄9.7997 -0.5513 -1.6539 -3.3078
-8.2695 -16.5390
25 太原9.7970 -0.2450 -0.7350 -1.4700
-3.6750 -7.3500
26 天津9.8011 -0.6636 -1.9908 -3.9816
-9.9540 -19.9080
27 武汉9.7936 0.1020 0.3060 0.6120 1.5300 3.0600
29 西安9.7944 0.0204 0.0612 0.1224 0.3060 0.6120
30 西宁9.7911 0.3267 0.9801 1.9602 9.8010 9.8010
31 张家口9.8000 -0.5513 -1.6539 -3.3078
-8.2695 -16.5390
32 郑州9.7966 -0.2041 -0.6123 -1.2246
-3.0615 -6.1230
ANSYS结构分析单元功能与特性
杆单元:LINK1、8、10、11、180
梁单元：BEAM3、4、23、24，44，54，188，189
管单元:PIPE16，17，18，20，59，60
2D实体元:PLANE2，25，42，82，83，145，146，182，183
3D实体元:SOLID45，46，64，65，72，73，92，95，147，148，185，186，187，191 壳单元:SHELL28，41，43，51，61，63，91，93，99，143，150，181，208，209
弹簧单元:COMBIN7，14，37，39，40
质量单元:MASS21
接触单元:CONTAC12，52，TARGE169，170，CONTA171，172，173，174，175，178
矩阵单元:MATRIX27，50
表面效应元:SURF153，154
粘弹实体元:VISCO88，89，106，107，108，
超弹实体元:HYPER56，58，74，84，86，158
耦合场单元:SOLID5，PLANE13，FLUID29，30，38，SOLID62，FLUID79，FLUID80，81，SOLID98，FLUID129，INFIN110，111，FLUID116，130
界面单元:INTER192，193，194，195
显式动力
分析单元:LINK160，BEAM161，PLANE162，SHELL163，SOLID164，COMBI16
杆单元
E-弹性(Elasticity)，P-塑性(Plasticity)，C-蠕变(Creep)，S-膨胀(Swelling)，D-大变形或大挠度(Large deflection)，F-大应变(Large strain)或有限应变(Finite strain)，B-单元生死(Birth and dead),G-应力刚化(Stress stiffness)或几何刚度(Geometric stiffening)，A-自适应下降(Adaptive descent)等。
通常用LINK1和LINK8模拟桁架结构，如屋架、网架、网壳、桁架桥、桅杆、塔架等结构，以及吊桥的吊杆、拱桥的系杆等构件，必须注意线性静力分析时，结构不能是几何可变的，否则造成位移超限的提示错误。LINK10可模拟绳索、地基弹簧、支座等，如斜拉桥的斜拉索、悬索、索网结构、缆风索、弹性地基、橡胶支座等。LINK180除不具备双线性特性(LINK10)外，它均可应用于上述结构中，并且其可应用的非线性性质更加广泛，增加了粘弹塑性材料。⑸LINK1、LINK8和LINK180单元还可用于普通钢筋和预应力钢筋的模拟，其初应变可作为施加预应力的方式
梁单元
梁单元分为多种单元，分别具有不同的特性，是一类轴向拉压、弯曲、扭转(3D)单元。
单元使用另外应注意的问题：
⑴梁单元面积和长度不能为零，且2D梁元必须位于XY平面内；⑵剪切变形的影响；⑶自由度释放；⑷梁截面特性；⑸BEAM23/24实常数的输入比较复杂；⑹荷载特性；⑺应力计算。
管单元
管单元是一类轴向拉压、弯曲和扭转的3D单元，单元的每个节点均具有6个自由度，即三个平动自由度Ux、Uy、Uz和三个转动自由度Rotx、Roty、Rotz，此类单元以3D梁元为基础，包含了对称性和标准管几何尺寸的简化特性。
单元使用应注意的其他问题：
⑴管元长度、直径及壁厚均不能为零；⑵可计算薄壁管和厚壁管，但某些应力的计算是基于薄壁管理论的；⑶管单元计入了剪切变形的影响，并可考虑应力增强系数和挠曲系数。
2D实体单元
2D实体单元是一类平面单元，可用于平面应力、平面应变和轴对称问题的分析，
单元使用应注意的其他问题：
⑴单元插值函数及说明；⑵荷载特性；⑶其它特点。
3D实体单元
3D实体单元用于模拟三维实体结构，此类单元每个节点均具有三个自由度，即Ux、Uy、Uz 三个平动自由度。
单元使用应注意的问题：
⑴关于SOLID72/73单元；（2）SOLID185积分方式可选择。
壳单元
壳单元可以模拟平板和曲壳一类结构。壳元比梁元和实体元要复杂的多，因此壳类单元中各种单元的选项很多。
杆、梁单元→板壳单元→实体单元
单元使用应注意的问题：
⑴通常不计剪切变形的壳元用于薄板壳结构，而计入剪切变形的壳元用于中厚度板壳结构。
弹簧单元
弹簧单元是一类专门模拟“弹簧”行为的单元，不同于用结构单元（如LINK等）的模拟。
质量单元
MASS21为具有6个自由度的点单元，即只有一个节点，节点自自由度可为Ux、Uy、Uz、Rotx、Roty、Rotz，通过不同设置可仅考虑2D或3D内的平动自由度及其组合，它每个坐标方向可以具有不同的质量和转动惯量。该单元无面荷载和体荷载，支持弹性、大变形和生死单元。
接触单元
ANSYS支持三种接触方式，即点对点、点对面和面对面的接触，接触单元是覆盖在模型单元的接触面之上的一层单元。点点单元用于模拟点对点的接触行为，且预先知道接触位置；点面单元用于模拟点对面的接触行为，预先不要确定接触位置，接触面之间的网格不要求一致；面面单元用于模拟面对面的接触行为，支持低阶和高阶单元，支持大变形行为等。
矩阵单元
MATRIX27为刚度、阻尼、质量矩阵单元，可表示一种任意的单元。本单元具有两个节点，此两个节点可重合或不重合，每个节点有6个自由度，即Ux、Uy、Uz、Rotx、Roty、Rotz。该单元无面荷载和体荷载，但支持单元生死功能。其矩阵可为对称或不对称形式，通过Keyopt(3)设置为刚度矩阵、或阻尼矩阵、或质量矩阵。本单元可模拟任意类型的单元，如可模拟特殊弹簧和节点柔性连接等。
MATRIX50为超单元，它是预先装配好的可独立使用的一组单元。该单元无节点和实常数，其自由度数目由所包含的单元决定，其面荷载和体荷载可通过总的载荷向量和比例系数施加，该单元支持大变形功能。该单元不能包含基于拉格朗日乘子的单元（如MPC184等），不支持非线性（忽略所包含的单元非线性）。超单元可包含其它超单元，2D超单元只能用于二维分析，而3D超单元则只能用于三维分析。
表面效应单元
SURF153和SURF154分别为2D和3D结构表面效应单元，可用于各种荷载（法向、切向、法向渐变、输入矢量方向等）及表面效应（基础刚度、表面张力及附加质量等）情况，可覆盖于任何二维（轴对称谐结构单元PLANE25/83除外）和三维结构实体单元表面。
预紧、多点约束、网分单元
（1）PRETS179为2D/3D预紧单元，用于定义网分后的二维或三维结构预紧区，可由任意结构单元（杆、梁、管、壳、2D实体和3D实体）建立。该单元具有3个节点，每个节点具有一个自由度Ux，该Ux为预紧方向的位移，ANSYS通过几何条件将预紧力施加到指定的预紧荷载方向上，而不必考虑模型是如何定义的。该单元不支持面荷载和体荷载，仅支持非线性特性；不能使用约束方程和自由度耦合，NROTAT命令不能用于节点K，且K节点必须位于整体直角坐标系。
（2）MPC184为多点约束单元，有刚性杆、刚性梁、滑移、球形、销钉、万向接头的约束，适用于使用拉格朗日乘子的具有运动约束时情况，该单元可用于机构运动学，如起重机、挖掘机、汽车、机床和机器人等。该单元有2个或3个节点，每个节点具有Ux、Uy(2D)或Ux、
Uy、Uz(3D) 或Ux、Uy、Uz、Rotx、Roty、Rotz(3D)自由度。无实常数和面荷载，支持温度荷载及转动或转动力矩，支持大变形和单元生死。
⑶MESH200是仅用来划分网格的单元，对计算结果毫无影响。它是为实现多步网格划分的操作而设计的。该单元可用于划分两维或三维空间的线，三维空间中的三角形、四边形、四面体或六面体单元组成的面或体，且均包括有或没有中间节点的情况。MESH200单元可与任意其它单元一起使用，当不再需要它时，可以将其删除或保留
坐标系和工作平面
6类坐标系：总体坐标系、局部坐标系、节点坐标系、单元坐标系、显示坐标系与结果坐标系。
激活总体和局部坐标系
命令：CSYS,KCN
其中KCN表示坐标系号码，0-直角坐标系（缺省），1-柱坐标系，2-球坐标系，4-以工作平面为坐标系，5-柱坐标系（以Y轴为转轴），≥11-局部坐标系。
由于工作平面可不断移动和旋转，因此当采用CSYS,4时也相当于不断定义了局部直角坐标，在很多情况下应用非常方便。
根据总体坐标系定义局部坐标系
命令：LOCAL, KCN, KCS, XC, YC, ZC, THXY, THYZ, THZX, PAR1, PAR2
其中：KCN---局部坐标系编号，此编号必须大于10，如果与既有编号相同，则将重新定义KCS---坐标系类型，0或CART为直角坐标系，1或CYLIN为柱坐标系，2或SPHE为球坐标系，3或TORO为环坐标系。
XC,YC,ZC---新坐标系原点在总体直角坐标系中的坐标。
THXY,THYZ,THZX---新坐标系绕Z,X,Y轴的旋转角度，其正方向为：XY,YZ,ZX。
PAR1---适用于椭圆、类似球体或环形系统，当KCS=1或2时，其值为椭圆Y轴半径与X轴半径之比，缺省为1即圆。当KCS=3时，其值为环面的主半径。
PAR2---仅适用于类似球体的系统，当KCS=2时，其值为椭球体Z轴半径与X轴半径之比，缺省为1
根据已有的三个节点定义局部坐标系命令：CS, KCN, KCS, NORIG, NXAX, NXYPL, PAR1, PAR2根据已有的三个关键点定义局部坐标系命令：CSKP, KCN, KCS, PORIG, PXAXS, PXYPL, PAR1, PAR2
根据当前工作平面定义局部坐标系命令：CSWPLA, KCN, KCS, PAR1, PAR2
根据激活的坐标系定义局部坐标系命令：CLOCAL, KCN, KCS, XL, YL, ZL, THXY, THYZ, THZX, PAR1, PAR2
删除局部坐标系命令：CSDELE, KCN1, KCN2, KCINC
其中：KCN1---为要删除的局部坐标系的起始编号，如果KCN1=ALL，则其后参数将忽略。KCN2---为要删除的局部坐标系的最终编号。
KCINC---为编号的递增数值，缺省为1。
CSDELE,11,15,2---则删除了11、13、15号局部坐标系。
查看激活坐标系和局部坐标系命令：CSLIST, KCN1, KCN2, KCINC
节点坐标系的旋转与修改
将某些节点的坐标系旋转到与当前激活坐标系（简称“当前坐标系”）方向一致
命令：NROTAT, NODE1, NODE2, NINC
其中NODE1、NODE2、NINC ---要旋转节点的起始号、末编号（缺省为NODE1）及递增值（缺省值为1）。如NODE1=ALL则其后参数将被忽略，NODE1也可为元件名。
将既有节点的节点坐标系旋转某个角度命令：NMODIF, NODE, X, Y, Z, THXY, THYZ, THZX NODE---节点号、ALL或元件名称。
X, Y, Z---该节点的新坐标值。其余参数意义同前。
在创建节点时直接定义其坐标系的旋转角度命令：N, NODE, X, Y, Z, THXY, THYZ, THZX 按方向余弦旋转节点坐标系命令：NANG, NODE, X1, X2, X3, Y1, Y2, Y3, Z1, Z2, Z3
节点坐标系列表命令：NLIST, NODE1, NODE2, NINC, Lcoord, SORT1, SORT2, SORT3 Lcoord---坐标列表信息，缺省为全部信息，=COORD时仅列XYZ坐标。
SORT1---用于排序的第1项内容，可以是
NODE,X,Y,Z,THXY,THYZ,THXZ。
SORT2,SORT3---用于排序的第2项和第3项内容，其内容同SORT1。
单元坐标系的定义与修改
设置单元坐标系命令：ESYS,KCN
其中KCN为坐标系编号，KCN=0（缺省）表示使用单元定义时规定的坐标系方向。当KCN=N（N>10）时使用编号为N的局部坐标系。
修改单元坐标系方向命令：EMODIF, IEL, STLOC, I1, I2, I3, I4, I5, I6, I7, I8
IEL---单元编号，或ALL，或元件名。
STLOC---将要修改的第一个节点序号或属性，属性之一为ESYS，则I1为局部坐标号。
激活显示坐标系
命令：DSYS,KCN其中KCN---坐标系号，可为0,1,2及局部坐标系号。缺省为总体直角坐标系。激活结果坐标系
命令：RSYS,KCN
其中KCN---坐标系号，可为0（缺省）,1,2及局部坐标系号。
当KCN=SOLU时，则与求解计算时采用的坐标系相同，实际上采用数据存储时的坐标系。定义工作平面
将既有坐标系的XY平面定义为工作平面命令：WPCSYS,WN,KCN
其中KCN为既有坐标系号，可以是0,1,2,或局部坐标系号。缺省为激活的坐标系。
通过3个坐标点定义工作平面
命令：WPLANE,WN,XORIG,YORIG,ZORIG,XXAX,YXAX,ZXAX,XPLAN,YPLAN,ZPLAN
通过3个节点定义工作平面命令：NWPLAN, WN, NORIG, NXAX, NPLAN
通过3个关键点定义工作平面命令：KWPLAN, WN, KORIG, KXAX, KPLAN
通过垂直于线上的某个位置定义工作平面命令：LWPLAN, WN, NL1, RATIO
工作平面的操控
工作平面的当前状态
stata save命令 FileSave As 例1. 表1.为某一降压药临床试验数据，试从键盘输入Stata，并保存为Stata格式文件。 STATA数据库的维护 排序 SORT 变量名1 变量名2 …… 变量更名 rename 原变量名新变量名 STATA数据库的维护 删除变量或记录 drop x1 x2 /* 删除变量x1和x2 drop x1-x5 /* 删除数据库中介于x1和x5间的所有变量(包括x1和x5) drop if x<0 /* 删去x1<0的所有记录 drop in 10/12 /* 删去第10~12个记录 drop if x==. /* 删去x为缺失值的所有记录 drop if x==.|y==. /* 删去x或y之一为缺失值的所有记录 drop if x==.&y==. /* 删去x和y同时为缺失值的所有记录 drop _all /* 删掉数据库中所有变量和数据 STATA的变量赋值 用generate产生新变量 generate 新变量＝表达式 generate bh=_n /* 将数据库的内部编号赋给变量bh。 generate group=int((_n-1)/5)+1 /* 按当前数据库的顺序，依次产生5个1，5个2，5个3……。直到数据库结束。 generate block=mod(_n,6) /* 按当前数据库的顺序，依次产生1,2,3,4,5,0。generate y=log(x) if x>0 /* 产生新变量y，其值为所有x>0的对数值log(x)，当x<=0时，用缺失值代替。 egen产生新变量 set obs 12 egen a=seq() /*产生1到N的自然数 egen b=seq(),b(3) /*产生一个序列，每个元素重复＃次 egen c=seq(),to(4) /*产生多个序列，每个序列从1到＃ egen d=seq(),f(4)t(6) /*产生多个序列，每个序列从#1到＃2 encode 字符变量名，gen(新数值变量名) 作用：将字符型变量转化为数值变量。 STATA数据库的维护 保留变量或记录 keep in 10/20 /* 保留第10～20个记录，其余记录删除 keep x1-x5 /* 保留数据库中介于x1和x5间的所有变量(包括x1和x5)，其余变量删除keep if x>0 /* 保留x>0的所有记录，其余记录删除
1. 如何输出STATA的图，和保存？ 先输入数据 （1）Twoway connected 变量1 变量2 //划出折线图 （2）twoway scatter 变量1 变量2 //划出散点图 2. 怎样在stata8中做HAUSMAN检验？ 四步曲，重点在于解释结果 （1）xtreg y x , fe （2）est store fe （3）xtreg y x, re （4）hausman fe 如果拒绝，说明corr(x,ui)=0的假设是有问题的，需要重新设定RE model 后再进行检验，如果模型的设定没有问题，但检验还是拒绝原假设（p值接近0），那么就只能采用FE model 了，因为此时的RE 估计量是有偏的。 （definitely right. 当你使用stata的时候，最重要的命令不是这些是help and find it然后就能找到你的答案了） hausman检验是用来检验用fe还是re的，其原假设是re优于fe，从你的结果来看（Prob>chi2 =0.0000），应该拒绝原假设，所以应该用fe 3.stata里平方的命令怎么写? gen age=age^2 4. stata里边怎么取对数啊？ gen lnx=log(x) 5.如何用STATA求自然对数？如说：ln(X^2)=-4.8536，如何求X啊？ . dis sqrt(exp(-4.8536))或者dis exp(-4.8536/2) 6.关于hausman检验，结果是CHI2(2)=2355.81，prob>chi2＝0.000，可以使用随机效应模型嘛？ prob>chi2＝0.000,is like p-value. we should reject the null, so fixed effect is preferred.Randome effect is not suggested. CHI2(2)=2355.81，就意味着拒绝原假设，从而选取固定效应模型。 7.我在做gdp一阶差分单位根检验的时候，输入的命令是ipshin dgdp,lags(1)得
Stata 命令语句格式： [by varlist:] command [varlist] [=exp] [if exp] [in range] [weight] [, options] 1、[by varlist:] *如果需要分别知道国产车和进口车的价格和重量，可以采用分类操作来求得， sort foreign //按国产车和进口车排序 . by foreign: sum price weight *更简略的方式是把两个命令用一个组合命令来写。 . by foreign, sort: sum price weight 如果不想从小到大排序，而是从大到小排序，其命令为gsort。 . sort - price //按价格从高到低排序 . sort foreign -price /*先把国产车都排在前，进口车排在后面，然后在国产车内再按价格从大小到排序，在进口车内部，也按从大到小排序*/ 2、[=exp]赋值运算 . gen nprice=price+10 //生成新变量nprice，其值为price+10 /*上面的命令generate(略写为gen) 生成一个新的变量，新变量的变量名为 nprice,新的价格在原价格的基础上均增加了10 元。 . replace nprice=nprice-10 /*命令replace 则直接改变原变量的赋值，nprice 调减后与price 变量取值相等*/ 3、[if exp]条件表达式 . list make price if foreign==0 *只查看价格超过1 万元的进口车（同时满足两个条件），则 . list make price if foreign==1 & price>10000 *查看价格超过1 万元或者进口车（两个条件任满足一个） . list make price if foreign==1 | price>10000 4、[in range]范围筛选 sum price in 1/5 注意“1/5”中，斜杠不是除号，而是从1 到 5 的意思，即1，2，3，4，5。 如果要计算前10 台车中的国产车的平均价格，则可将范围和条件筛选联合使用。 . sum price in 1/10 if foreign==0 5、[weight] 加权 sum score [weight=num] 其中，num为每个成绩所对应的人数 6、[, options]其他可选项 例如，我们不仅要计算平均成绩，还想知道成绩的中值，方差，偏度和峰度等*/ . sum score, detail . sum score, d //d 为detail 的略写，两个命令完全等价 . list price, nohead //不要表头 Stata 数据类型转换 1、字符型转化成数值型 destring, replace //全部转换为数值型，replace 表示将原来的变量（值）更新 destring date, replace ignore(“ ”) 将字符型数据转换为数值型数据：去掉字符间的空格destring price percent, gen(price2 percent2) ignore(“$ ,%”) 与date 变量类似，变量price 前面有美元符号，变量percent 后有百分号，换为数值型时需要忽略这些非数值型字符 2、数值型转化为字符型
一、波浪理论的基本特点 波浪理论认为一个完整的波浪分为八浪循环，如图 波浪理论有以下基本特点： （1）上升浪和下跌浪永远交替运行。 （2）基本形态为顺势五浪，即推动浪用1、2、3、4、5表示；逆势三浪即调整浪用a、b、c表示。 （3）推动浪可划分为低一级的5个子浪。调整浪可划分为低一级的3个子浪，即浪中有浪， （4）一个完整的八浪循环形成后，1、2、3、4、5浪可以合并为高一级的（1）、（3）、（5）浪，a、b、c浪可以合并为高一级的（2）和（4）浪。 （5）无论股价的上涨、下跌的幅度有多大或多小，也无论时间有多长或多短，都不会改变波浪的形态。因为市场仍会按照基本形态发展，波浪可以延长，也可以缩短。但基本形态永远不会改变，改变的是涨跌的幅度和运行时间。 二、推动浪 推动浪基本上就是1、2、3、4、5浪，且一浪高于一浪，而每一个浪又可细分为5个子浪。形态较简单，分为延伸浪、倾斜三角形和失败形态三种。其特征如下： （一）延伸浪
（1）在（1）（3）（5）浪中只有一个浪会出现延伸浪， （2）如果（1）、（3）浪长度差不多，而且（5）浪成交量明显比（3）浪大时，则（5）浪出现延伸浪机会最大。 （3）如果延伸现象出现在（3）浪，那么（5）浪形态简单，（5）浪的长度和时间与（1）浪相当。 （4）如果延伸出现在（5）浪，那么（5）浪会以二次回档出现；第一次回档下跌至延伸浪的起点附近；第二次回档是反弹上升创 出新高之后的回档。如图所示： 当第5浪成为延伸浪而属于高一级的（1）浪或（3）浪时，则第一次回档至延伸浪的起点，即是高一级的（2）浪、（4）浪的低点；第二回档则反弹上升形成（3）浪或（5）浪并创出新高。 （二）倾斜三角形 倾斜三角形只可能在推动浪（5）浪中发生，面且在其发生之前，通常有一段非常迅速的暴涨和暴跌，然后其成交量逐渐递减，并且波动幅度减缓，演变为倾斜三角形。当倾斜三角形发生在头部时，暗示股价将向下反转，下降途中出现倾斜三角形暗示股价将向上反转。如图所示： 倾斜三角形的特点如下：
stata11 常用命令 注：JB统计量对应的p大于0.05 ，则表明非正态，这点跟sktest 和 swilk 检验刚好相反；dta 为数据文件；gph 为图文件；do 为程序文件；注 意stata 要区别大小写；不得用作用户变量名： _all _n _N _skip _b _coef _cons _pi _pred _rc _weight double float long int in if using with 命令：读入数据一种方式 input x y 14 2 5.5 3 6.2 47.7 58.5 end su/summarise/sum x 或su/summarise/sum x,d 对分组的描述： sort group by group:su x %%%%% tabstat economy,stats（max）%返回变量economy的最大值 %%stats括号里可以是：mean，count（非缺失观测值个数），sum（总 和），max，min，range ， %% sd ，var ，cv（变易系数＝标准差/ 均值），skewness，kurtosis ， median，p1（1 ％分位 %% 数，类似地有p10, p25, p50, p75, p95, p99），iqr（interquantile range = p75 –p25） _all %描述全部 _N 数据库中观察值的总个数。 _n 当前观察值的位置。 _pi 圆周率π 的数值。 list gen/generate % 产生数列egen wagemax=max（wage） clear use by（分组变量）
波浪理论内容的几个基本的要点(１)一个完整的循环包括八个波浪，五上三落。(２)波浪可合并为高一级的浪，亦可以再分割为低一级的小浪。(３)跟随主流行走的波浪可以分割为低一级的五个小浪。(４)１、３、５三个推浪中，第３浪不可以是最短的一个波浪。(５)假如三个推动论中的任何一个浪成为延伸浪，其余两个波浪的运行时间及幅度会趋一致。(６)调整浪通常以三个浪的形态运行。(７)黄金分割率奇异数字组合是波浪理论的数据基础。(８)经常遇见的回吐比率为０．３８２、０．５及０．６１８。(９)第四浪的底不可以低于第一浪的顶。(１０)波浪理论包括三部分：型态、比率及时间，其重要性以排行先后为序。(１１)波浪理论主要反映群众心理。越多人参与的市场，其准确性越高。 来源于: 股票学习资料网(https://www.doczj.com/doc/166216494.html,) 艾略特波浪理论图解 艾略特波浪理论对许多投资者来说可以算得上是如雷灌耳，但是能够真正理解并正确运用波浪理论的人少之又少，为了能够广泛普及艾略特波浪理论基础知识，让更多的投资者能够准确掌握这里经典证券分析理论的技术，“波浪理论网”收集整理了大量相关知识，并结合图形对该理论进行细致的讲解，条理清晰，通俗易懂。以下就是艾略特波浪理论图解。 波浪理论是由 Nalph Nelson Eilliott在1938年所提出的，波浪理论是技术分析大师R·E·艾略特（R·E·Elliot）所发明的一种价格趋势分析工具，它是一套完全靠而观察得来的规律，可用以分析股市指数、价格的走势，它也是世界股市分析上运用最多，而又最难于了解和精通的分析工具。 许多从事过波浪理论研究并在实际操作中付诸实施的投资者都曾会感到波浪理论不易领会，甚至望而生畏。波浪理论的基本原则其实很简单，读者在不久将会发现波浪理论涵盖的许多要点。看起来似曾相识，这是因为波浪理论的许多架构，相当符合道氏理论的原理和传统的图型技术。不过，波浪理论已超越传统的图型分析技术，能够针对市场的波动，提供全盘性的分析角度、得以解释特定的图形型态发展的原因与时机，以及图形本身所代表的意义，波浪理论同时也能够帮助市场分析师、找出市场循环周期的所在。
[推荐] Stata基本操作汇总——常用命令 help和search都是查找帮助文件的命令，它们之间的 区别在于help用于查找精确的命令名，而search是模糊查找。 如果你知道某个命令的名字，并且想知道它的具体使用方法，只须在stata的命令行窗口中输入help空格加上这个名字。回车后结果屏幕上就会显示出这个命令的帮助文件的全部 内容。如果你想知道在stata下做某个估计或某种计算，而 不知道具体该如何实现，就需要用search命令了。使用的 方法和help类似，只须把准确的命令名改成某个关键词。回车后结果窗口会给出所有和这个关键词相关的帮助文件名 和链接列表。在列表中寻找最相关的内容，点击后在弹出的查看窗口中会给出相关的帮助文件。耐心寻找，反复实验，通常可以较快地找到你需要的内容.下面该正式处理数据了。我的处理数据经验是最好能用stata的do文件编辑器记下你做过的工作。因为很少有一项实证研究能够一次完成，所以，当你下次继续工作时。能够重复前面的工作是非常重要的。有时因为一些细小的不同，你会发现无法复制原先的结果了。这时如果有记录下以往工作的do文件将把你从地狱带到天堂。因为你不必一遍又一遍地试图重现做过的工作。在stata 窗口上部的工具栏中有个孤立的小按钮，把鼠标放上去会出
现“bring do-file editor to front”，点击它就会出现do文件编 辑器。 为了使do文件能够顺利工作，一般需要编辑do文件的“头”和“尾”。这里给出我使用的“头”和“尾”。capture clear （清空内存中的数据）capture log close （关闭所有 打开的日志文件）set more off （关闭more选项。如果打开该选项，那么结果分屏输出，即一次只输出一屏结果。你按空格键后再输出下一屏，直到全部输完。如果关闭则中间不停，一次全部输出。）set matsize 4000 （设置矩阵的最大阶数。我用的是不是太大了？）cd D: （进入数据所在的盘符和文件夹。和dos的命令行很相似。）log using （文件名）.log,replace （打开日志文件，并更新。日志文件将记录下所有文件运行后给出的结果，如果你修改了文件内容，replace选项可以将其更新为最近运行的结果。）use （文件名）,clear （打开数据文件。）（文件内容）log close （关闭日志文件。）exit,clear （退出并清空内存中的数据。） 实证工作中往往接触的是原始数据。这些数据没有经过整理，有一些错漏和不统一的地方。比如，对某个变量的缺失观察值，有时会用点，有时会用-9，-99等来表示。回归时如果 使用这些观察，往往得出非常错误的结果。还有，在不同的数据文件中，相同变量有时使用的变量名不同，会给合并数
*********面板数据计量分析与软件实现********* 说明：以下do文件相当一部分内容来自于中山大学连玉君STATA教程，感谢他的贡献。本人做了一定的修改与筛选。 *----------面板数据模型 * 1.静态面板模型：FE 和RE * 2.模型选择：FE vs POLS, RE vs POLS, FE vs RE （pols混合最小二乘估计） * 3.异方差、序列相关和截面相关检验 * 4.动态面板模型（DID-GMM,SYS-GMM） * 5.面板随机前沿模型 * 6.面板协整分析（FMOLS,DOLS） *** 说明：1-5均用STATA软件实现， 6用GAUSS软件实现。 * 生产效率分析（尤其指TFP）：数据包络分析（DEA）与随机前沿分析（SFA） *** 说明：DEA由DEAP2.1软件实现，SFA由Frontier4.1实现，尤其后者，侧重于比较C-D与Translog生产函数，一步法与两步法的区别。常应用于地区经济差异、FDI溢出效应（Spillovers Effect）、工业行业效率状况等。 * 空间计量分析：SLM模型与SEM模型 *说明：STATA与Matlab结合使用。常应用于空间溢出效应（R&D）、财政分权、地方政府公共行为等。 * --------------------------------- * --------一、常用的数据处理与作图----------- * --------------------------------- * 指定面板格式 xtset id year （id为截面名称，year为时间名称） xtdes /*数据特征*/ xtsum logy h /*数据统计特征*/ sum logy h /*数据统计特征*/ *添加标签或更改变量名 label var h "人力资本" rename h hum *排序 sort id year /*是以STATA面板数据格式出现*/ sort year id /*是以DEA格式出现*/ *删除个别年份或省份 drop if year<1992 drop if id==2 /*注意用==*/ *如何得到连续year或id编号（当完成上述操作时，year或id就不连续，为形成panel格式，需要用egen命令） egen year_new=group(year) xtset id year_new **保留变量或保留观测值 keep inv /*删除变量*/ **或 keep if year==2000 **排序 sort id year /*是以STATA面板数据格式出现 sort year id /*是以DEA格式出现 **长数据和宽数据的转换 *长>>>宽数据 reshape wide logy,i(id) j(year)
波浪理论在具体运用中，常常会遇到较为难以分辨的市况，发现几个同时可以成立的数浪方式。所以，投资者有必要了解各个波浪的特性。第一浪在整个波浪循环开始后，一般市场上大多数投资者并不会马上就意识到上升波段已经开始。所以，在实际走势中，大约半数以上的第一浪属于修筑底部形态的一部分。由于第一浪的走出一般产生于空头市场后的末期，所以，市场上的空头气氛以及习惯于空头市场操作的手法未变，因此，跟随着属于筑底一类的第一浪而出现的第二浪的下调幅度，通常都较大。第二浪上面已经提过，通常第二浪在实际走势中调整幅度较大，而且还具有较大的杀伤力，这主要是因为市场人士常常误以为熊市尚未结束，第二浪的特点是成交量逐渐萎缩，波动幅度渐渐变窄，反映出抛盘压力逐渐衰竭，出现传统图形中的转向形态，例如常见的头肩、双底等。第三浪第三浪在绝大多数走势中，属于主升段的一大浪，因此，通常第三浪属于最具有爆炸性的一浪。它的最主要的特点是：第三浪的运行时间通常会是整个循环浪中的最长的一浪，其上升的空间和幅度亦常常最大；第三浪的运行轨迹，大多数都会发展成为一涨再涨的延升浪；在成交量方面，成交量急剧放大，体现出具有上升潜力的量能；在图形上，常常会以势不可挡的跳空缺口向上突破，给人一种突破向上的强烈讯号。第四浪从形态的结构来看，第四浪经常是以三角形的调整形态进行运行。第四浪的运行结束点，一般都较难预见。同时，投资者应记住，第四浪的浪底不允许低于第一浪的浪顶。第五浪在股票市场中，第五浪是三大推动浪之一，但其涨幅在大多数情况下比第三浪小。第五浪的特点是市场人气较为高涨，往往乐观情绪充斥整个市场。从其完成的形态和幅度来看，经常会以失败的形态而告终。在第五上升浪的运行中，二、三线股会突发奇想，普遍上升，而常常会升幅极其可观。A浪在上升循环中，A浪的调整是紧随着第五浪而产生的，所以，市场上大多数人士会认为市势仍未逆转，毫无防备之心，只看作为一个短暂的调整。A浪的调整形态通常以两种形式出现，平坦型形态与三字形形态，它与B浪经常以交叉形式进行形态交换。B浪B浪的上升常常会作为多方的单相思，升势较为情绪化，这主要是市场上大多数人仍未从牛市冲天的市道中醒悟过来，还以为上一个上升尚未结束，在图表上常常出现牛市陷阱，从成交量上看，成交稀疏，出现明显的价量背离现象，上升量能已接济不上。C浪紧随着B浪而后的是C浪，由于B浪的完成顿使许多市场人士兵醒悟，一轮多头行情已经结束，期望继续上涨的希望彻底破灭，所以，大盘开始全面下跌，从性质上看，其破坏力较强。
Stata软件基本操作和数据分析入门 第一讲 Stata操作入门 张文彤赵耐青 第一节概况 Stata最初由美国计算机资源中心（Computer Resource Center）研制，现在为Stata公司的产品，其最新版本为7.0版。它操作灵活、简单、易学易用，是一个非常有特色的统计分析软件，现在已越来越受到人们的重视和欢迎，并且和SAS、SPSS一起，被称为新的三大权威统计软件。 Stata最为突出的特点是短小精悍、功能强大，其最新的7.0版整个系统只有10M左右，但已经包含了全部的统计分析、数据管理和绘图等功能，尤其是他的统计分析功能极为全面，比起1G以上大小的SAS系统也毫不逊色。另外，由于Stata在分析时是将数据全部读入内存，在计算全部完成后才和磁盘交换数据，因此运算速度极快。 由于Stata的用户群始终定位于专业统计分析人员，因此他的操作方式也别具一格，在Windows席卷天下的时代，他一直坚持使用命令行／程序操作方式，拒不推出菜单操作系统。但是，Stata的命令语句极为简洁明快，而且在统计分析命令的设置上又非常有条理，它将相同类型的统计模型均归在同一个命令族下，而不同命令族又可以使用相同功能的选项，这使得用户学习时极易上手。更为令人叹服的是，Stata语句在简洁的同时又拥有着极高的灵活性，用户可以充分发挥自己的聪明才智，熟练应用各种技巧，真正做到随心所欲。
除了操作方式简洁外，Stata的用户接口在其他方面也做得非常简洁，数据格式简单，分析结果输出简洁明快，易于阅读，这一切都使得Stata成为非常适合于进行统计教学的统计软件。 Stata的另一个特点是他的许多高级统计模块均是编程人员用其宏语言写成的程序文件（ADO文件），这些文件可以自行修改、添加和下载。用户可随时到Stata网站寻找并下载最新的升级文件。事实上，Stata的这一特点使得他始终处于统计分析方法发展的最前沿，用户几乎总是能很快找到最新统计算法的Stata程序版本，而这也使得Stata自身成了几大统计软件中升级最多、最频繁的一个。 由于以上特点，Stata已经在科研、教育领域得到了广泛应用，WHO的研究人员现在也把Stata作为主要的统计分析工作软件。 第二节 Stata操作入门 一、Stata的界面 图1即为Stata 7.0启动后的界面，除了Windows版本的软件都有的菜单栏、工具栏，状态栏等外，Stata的界面主要是由四个窗口构成，分述如下： 1．结果窗口：位于界面右上部，软件运行中的所有信息，如所执行的命令、执行结果和出错信息等均在这里列出。窗口中会使用不同的颜色区分不同的文本，如白色表示命令，红色表示错误信息。 2．命令窗口：位于结果窗口下方，相当于DOS软件中的命令行，此处用于键入需要执行的命令，回车后即开始执行，相应的结果则会在结果窗口中显示出来。
波浪理论的计算方法 1)第一浪只是推动浪开始 2)第二浪调整不能超过第一波浪起点 比率: 2浪=1浪0.5或0.618 3)第三浪通常是最长波浪，但绝不能是最短(相对1浪和5浪长度) 比率: 3浪=1浪1.618, 2或2.618倍 4)第四浪的调整不能与第一浪重迭(楔形除外) 比率: 4浪=3浪0.382倍。 5)第五浪在少数情况下未能超第三浪终点，即以失败形态告终 比率: 5浪=1浪或5浪=(1浪-3浪)0.382、0.5、0.618倍。 6)A浪比率: A浪=5浪0.5或0.618倍。 7)B浪比率: B浪=A浪0.382、0.5、0.618倍。 8)C浪比率: C浪=A浪1倍或0.618、1.382、1.618倍。 1、波浪理论基础 1) 波浪理论由8浪组成、1、3、5浪影响真正的走势，无论是下跌行情还是上升行情， 都在这三个浪中赚钱； 2) 2、4浪属于逆势发展（回调浪） 3) 6、7、8浪属于修正浪（汇价短期没有创新低或新高） 2、波浪理论相关法则 1) 第3永远不是最短的浪 2) 第4浪不能跌破第2浪的低点，或不能超过第2浪的高点 3) 数浪要点：你看到的任何一浪都是第1浪，第2浪永远和你真正的趋势相反； 4) 数浪规则：看到多少浪就是多少浪，倒回去数浪； 3、相关交易法则 1) 第3浪是最赚钱的一浪，我们应该在1、3、5浪进行交易，避免在2、4浪进场以 及避免在2、4浪的低点或者高点挂单，因为一旦上破或者下坡前期高点或者低点，则会出现发转，具体还要配合RSI和MACD指标进行分析；
4、波浪理论精华部分 1) 波浪理论中最简单的一个循环，或者说最小的一个循环为两浪循环，即上升浪或下跌浪+回调浪 2) 每一波上升浪或下跌浪由5个浪组成，这5浪中有两次2T确认进场； 3) 每一波回调浪由3个浪组成，这3浪中只有一次2T确认进场； 4) 波浪和移动均线共振时，得出进场做多、做空选择，同时要结合4R法则以及123法则进行分析 波浪理论图解 2011-10-21 19:14 每位投资者都希望能预测未来,波浪理论正是这样一种价格趋势分析工具,它根据周期循环的波动规律来分析和预测价格的未来走势。波浪理论的创始人——美国技术分析大师R.N.艾略特（1871～1948）正是在长期研究道琼斯工业平均指数的走势图后,于二十世纪三十年代创立了波浪理论。投资者一走进证券部就会看到记录着股价波动信息的K线图,它们有节奏、有规律地起伏涨落、周而复始,好像大海的波浪一样,我们也可以感受到其中蕴涵的韵律与协调。我们特别邀请到了研究波浪理论的资深专家杨青老师来与读者们一起“冲浪”。 1、基础课波浪理论在技术分析中被广泛采用波浪理论最主要特征就是它的通用性。人类社会经济活动的许多领域都遵循着波浪理论的基本规律,即在相似和不断再现的波浪推动下重复着自己。因为股票、债券的价格运动是在公众广泛参与的自由市场之中,市场交易记录完整,与市场相关的信息全面丰富,因此特别适于检验和论证波浪理论,所以它是诸多股票技术分析理论中被运用最多的,但不可否认,它也是最难于被真正理解和掌握的。专家导读：被事实验证的传奇波浪波浪理论的初次亮相极富传奇色彩。1929年开始的全球经济危机引发了经济大萧条,美国股市在1929年10月创下386点的高点后开始大崩盘,到 1932年仲夏时节,整个市场弥漫着一片绝望的气氛。这时,波浪理论的始作俑者艾略特给《美国投资周刊》主编格林斯发电报,明确指出长期下跌的走势已经结束,未来将会出现一个大牛市。当格林斯收到电报时,道琼斯30种工业指数已经大幅飙升,从邮戳上的时间看,电报就在道琼斯30种工业指数见底前两个小时发出。此后道琼斯指数在9周内上涨了100％,而且从此开始一路上扬。 但是波浪理论在艾略特生前却长期被人们忽视,直到1978年,他的理论继承者帕彻特出版了《波浪理论》一书,并在期货投资竞赛中运用波浪理论取得了四个月获利400％以上的骄人成绩后,这一理论才被世人广泛关注,并开始迅速传播。 2、波浪周期及实例解读 0 && image.height>0){if(image.width>=700){this.width=700;this.height=image .height*700/image.width;}}> 专家解读：五浪上升三浪下降组成完整周期一个完整的波动周期,即完成所谓从牛市到熊市的全过程,包括一个上升周期和一个下跌周期。上升周期由五浪构成,用1、2、3、4、5表示,其中1、3、 5浪上涨,2、4浪下跌；下跌周期由三浪构成,用a、b、c表示,其中a、c浪下跌,b 浪上升。与主趋势方向（即所在周期指明的大方向）相同的波浪我们称为推动浪,
波浪理论第一浪 在整个波浪循环开始后，一般市场上大多数投资者并不会马上就意识到上升波段已经开始。所以，在实际走势中，大约半数以上的第一浪属于修筑底部形态的一部分。由于第一浪的走出一般产生于空头市场后的末期，所以，市场上的空头气氛以及习惯于空头市场操作的手法未变，因此，跟随着属于筑底一类的第一浪而出现的第二浪的下调幅度，通常都较大。 第一浪的性格： 介绍第一浪的特点，为方便解释起见，暂以第一浪上升的走势作为介绍的基础。事实上，第一浪亦可以向下走。遇到向下走的推动浪，可将下述概念按相反方向解释。第一浪的开始，意味调整市势已经完结。因此，第一浪实际上是市势转变的标志。第一浪可以在划分为低一级的五组波浪，例如每小时走势图上所显示的波浪。通常而言，当市势出现三个波浪的调整走势，其后出现第一浪，足以证明市势向下调整完毕，价位将会掉头上升。当第一浪开始运行的时候，初时可能不易辨认，但第一浪一旦走完全程，便会构成可靠的讯号，其它推动浪将会逐一露面。第二浪接着出现，但其调整的幅度不应大于第一浪运行的长度。假如第三浪属于延伸浪，第五浪倾向于第一浪的长度相同。因此，第一浪可以用于预测第五浪见顶的上升目标。大约半数的第一浪属于营造底部形态的一部分，跟随这类第一浪出现的第二浪，调整幅度较大，但无论如何，回吐的比率不可以大于第一浪的100％。 波浪理论第二浪 通常第二浪在实际走势中调整幅度较大，而且还具有较大的杀伤力，这主要是因为市场人士常常误以为熊市尚未结束，第二浪的特点是成交量逐渐萎缩，波动幅度渐渐变窄，反映出抛盘压力逐渐衰竭，出现传统图形中的转向形态，例如常见的头肩、双底等。 第二浪考验判断力：第二浪的终点，通常会在下列三个地区出现。 （一）可能调整第一浪的38.2％或61.8％。
第六章stata语言中的常用函数 本章重点： Stata系统是一个统计分析系统，stata语言是实现stata系统功能的基础，因此它其中包括了各种各样的函数。在stata系统中，函数的自变量可以是一个常数，可以是一个变量，或者是一连串的变量。在调用这些函数的时候，只要将函数中定义中的这些变量替换为相应值即可。这一章，介绍一下这些函数的定义以及使用方法。 6.1函数概览 函数只不过是一些编号的小程序，它会按一定的规则进行处理，之后报告结果。实际上，谁也记不住这么多函数，因此，首先要学会查找函数的帮助，当记不住的时候，随时去查寻帮助。记住下面的命令才是最关键的。 . help function 弹出来的对话框告诉我们，STATA包括八类函数，分别是数学函数，分布函数，随机数函数，字符函数，程序函数，日期函数，时间序列函数和矩阵函数。本章主要介绍数学函数和字符函数，日期函数，随机函数等常用函数，其他函数可以参考stata 帮助功能。 6.2数学函数 Abs(x) x的绝对值
●Acos（x）反余弦函数 例如：arcos (0.5)=1.57 arcos(1)=0 ●Asin(x) 反正弦函数 ●Atan(x) 反正切函数 ●atanh(x) 反双曲正切函数 ●ceil(x) 返回大于或等于自变量的最小的整数。 例如：ceil(0.7)=1 ceil(3)=3 ceil(-0.7)=0 ●Floor(x) 返回小于或等于自变量的最大的整数 例如：floor(0.7)=0 floor(3)=3 floor(-0.7)=-1 ●Int(x) 返回自变量的整数部分 例如：int(0.7)=0 int(2.9)=2 int(-2.55)=-2 ●Round(x,y) 返回与y的单位最接近的数x，x为真数，y为近似单 位 例如：round(5.2,1)= round(4.8,1)=5 round(2.234,0.1)=2.2 round(2.234,0.01)=2.23 round(2.234,0.001)=2.234 round(28,5)=30 ●cloglog(x) 返回ln{-ln(1-x)}的值 ●comb(n,k) 从n中取k个的组合，即comb(n,k)=n!/{k!(n - k)!} 例如：comb(10,5)=252 comb(6,2)=15 ●cos(x) 余弦函数 ●digamma(x) 返回digamma函数值，这是lngamma(x)的一阶导数●exp(x) 指数函数
常用到的sta命令 闲话不说了。help和search都是查找帮助文件的命令，它们之间的区别在于help用于查找精确的命令名，而search是模糊查找。如果你知道某个命令的名字，并且想知道它的具体使用方法，只须在sta的命令行窗口中输入help空格加上这个名字。回车后结果屏幕上就会显示出这个命令的帮助文件的全部内容。如果你想知道在sta下做某个估计或某种计算，而不知道具体该如何实现，就需要用search命令了。使用的方法和help类似，只须把准确的命令名改成某个关键词。回车后结果窗口会给出所有和这个关键词相关的帮助文件名和链接列表。在列表中寻找最相关的内容，点击后在弹出的查看窗口中会给出相关的帮助文件。耐心寻找，反复实验，通常可以较快地找到你需要的内容。 下面该正式处理数据了。我的处理数据经验是最好能用sta的do文件编辑器记下你做过的工作。因为很少有一项实证研究能够一次完成，所以，当你下次继续工作时。能够重复前面的工作是非常重要的。有时因为一些细小的不同，你会发现无法复制原先的结果了。这时如果有记录下以往工作的do文件将把你从地狱带到天堂。因为你不必一遍又一遍地试图重现做过的工作。在sta窗口上部的工具栏中有个孤立的小按钮，把鼠标放上去会出现“bring do-file editor to front”，点击它就会出现do文件编辑器。 为了使do文件能够顺利工作，一般需要编辑do文件的“头”和“尾”。这里给出我使用的“头”和“尾”。 /*（标签。简单记下文件的使命。）*/ capture clear（清空内存中的数据） capture log close（关闭所有打开的日志文件） set mem 128m（设置用于sta使用的内存容量） set more off（关闭more选项。如果打开该选项，那么结果分屏输出，即一次只输出一屏结果。你按空格键后再输出下一屏，直到全部输完。如果关闭则中间不停，一次全部输出。） set matsize4000（设置矩阵的最大阶数。我用的是不是太大了？）
Stata统计分析常用命令汇总 一、winsorize极端值处理 范围：一般在1%和99%分位做极端值处理，对于小于1%的数用1%的值赋值，对于大于99%的数用99%的值赋值。 1、Stata中的单变量极端值处理： stata 11.0，在命令窗口输入“findit winsor”后，系统弹出一个窗口，安装winsor模块 安装好模块之后，就可以调用winsor命令，命令格式：winsor var1, gen(new var) p(0.01) 或者在命令窗口中输入：ssc install winsor安装winsor命令。winsor命令不能进行批量处理。 2、批量进行winsorize极端值处理： 打开链接：https://www.doczj.com/doc/166216494.html,/judson.caskey/data.html，找到winsorizeJ，点击右键，另存为到stata中的ado/plus/目录下即可。命令格式：winsorizeJ var1var2var3,suffix(w)即可，这样会生成三个新变量，var1w var2w var3w，而且默认的是上下1%winsorize。如果要修改分位点，则写成如下格式：winsorizeJ var 1 var2 var3,suffix(w) cuts(5 95)。 3、Excel中的极端值处理：（略） winsor2 命令使用说明 简介：winsor2 winsorize or trim (if trim option is specified) the variables in varlist at particular percentiles specified by option cuts(# #). In defult, new variables will be generated with a suffix "_w" or "_tr", which can be changed by specifying suffix() option. The replace option replaces the variables with their winsorized or trimmed ones. 相比于winsor命令的改进： (1) 可以批量处理多个变量； (2) 不仅可以winsor，也可以trimming； (3) 附加了by() 选项，可以分组winsor 或trimming； (4) 增加了replace 选项，可以不必生成新变量，直接替换原变量。 范例： *- winsor at (p1 p99), get new variable "wage_w" . sysuse nlsw88, clear . winsor2 wage *- left-trimming at 2th percentile . winsor2 wage, cuts(2 100) trim *- winsor variables by (industry south), overwrite the old variables . winsor2 wage hours, replace by(industry south) 使用方法: 1. 请将winsor 2.ado 和winsor2.sthlp 放置于stata12\ado\base\w 文件夹下； 2. 输入help winsor2 可以查看帮助文件；
【命令1】：导入数据 一般做实证分析使用的是excel中的数据，其后缀名为.xls，需要将其修改为.csv insheet using name.csv, clear 【命令2】：删除重复变量 sort var1 var2 duplicatesdrop var1 var2, force 【命令3】：合并数据 use data1, clear merge m:m var1 var2 using data2 drop if _merge==2 drop if _merge==1 drop _merge 【命令4】：描述性统计分析 tabstat var1var2, stat(n min mean median p25 p75 max sd), if groupvar==0 or 1 输出到word中： logout, save(name) word replace: tabstat var, stat(n min mean p50 max sd) col(stat)f(%9.2g) 【命令5】：结果输出 安装 ssc install estout, replace 单个回归 reg y x esttab using name.rtf, compress nogap r2 ar2 star(* 0.1 ** 0.05 *** 0.01) 多个回归一起 reg y x1 est store m1 reg y x2 est store m2 esttab m1 m2 using name.rtf, compress nogap r2 ar2 star(* 0.1 ** 0.05 *** 0.01)
波浪理论的基本结构 波浪理论的基本结构由八个浪组成（五升三跌或五跌三升），在一个上升趋势中，由五个上升浪（推动浪）和三个下跌浪（调整浪）组成；在下降趋势中，由五个下跌浪（推动浪）和三个上升浪（调整浪）组成。推动浪与主趋势方向一致，而调整浪则与主趋势方向相反。 图 185 是一个在上升趋势中的八浪循环示意图，左边 1－2－3－4－5 为推动浪，右边A－B－C是调整浪。如果一个波浪的趋势方向和它高一个层次的波浪的趋势方向相同，那么这一波浪就定义为主浪。主浪包括：第一浪、第三浪、第五浪、A 浪和 C 浪。 调整浪是指运行方向同它的上一层次的波浪（主浪）的运行方向相反的波浪，是对前一浪的调整。虽然调整浪比较复杂，数起来比较困难。但是他们有一个共同的特点就是：都以三浪或三浪的变形出现，决不会以五浪的形式出现。调整浪包括：第二浪、第四浪、B 浪和 ABC 浪。ABC 是对第一浪至第五浪的调整。 （1）主浪与调整浪： l主浪的扩延只能出现五浪形式，决不可以出现三浪的形式，如图 186 所示。 l调整浪的细分只能出现三浪的形式，决不能出现五浪形式，如图 187 所示。 （2）调整浪的形态：调整浪出现的形式通常有两种：5－3－5 锯齿形和 3－3－5 平坦形。5－3－5 锯齿形是指将调整浪 A 再细分为五个浪，调整浪 B 再分为 3 个浪，调整浪 C 再分为 5 个浪；3－3－5 平坦形是指将调整浪 A 再细分为 3 个浪，调整浪 B 再分为 3 个浪，调整浪 C 再分为 5 个浪。 如果调整浪以５－３－５锯齿形出现，Ｃ浪的长度等于Ａ浪的１倍或１．６１８倍，其标志是Ｃ浪穿越Ａ浪的低点，如图１８８所示。 如果调整浪以３－３－５平坦形出现，Ｃ浪的长度等于Ａ浪的长度，如图１８９所示。