土石坝有限元分析(ANSYS)-渗流分析命令流
土石坝有限元分析（ANSYS）－渗流分析命令流
土石坝渗流分析，采用非饱和土渗流参数，迭代计算浸润线，根据前次计算结果，不断修改单元的渗透系数和浸润逸出点位置，直到满足精度要求。本算例的土石坝体型比较简单.采用非饱和渗流计算.即渗透系数为空隙压力的函数.首先建立一个数据文件PPPP.TXT，存储渗透系数函数关系，如下。第一列为空隙压力值（水头M)，第二列为渗透系数指数，渗透系数等于10^A(M/D)。
! -10.00 -4.0E+00
! -9.00 -3.6E+00
! -8.00 -3.2E+00
! -7.00 -2.8E+00
! -6.00 -2.4E+00
! -5.00 -2.0E+00
! -4.00 -1.6E+00
! -3.00 -1.2E+00
! -2.00 -8.0E-01
! -1.00 -4.0E-01
! 0.00 0.0E+00
!土坝顶宽4M，上下游坡比均为1:2，总高12M，底宽52M。上游水深8M，下游无水。
FINISH
/CLEAR
/TITLE, EARTHDAM SEEPAGE
/FILNAME,SEEPAGE5
/PLOPTS,DATE,0
*DIM,TPRE,TABLE,11,1,1,PRESS,KKPE ! 定义水压与渗透系数的关系数组
*TREAD,TPRE,PPPP,TXT ! 读入数组
*DIM,NCON,ARRAY,4 ! 定义数组，用于存贮单元四个节点号
/PREP7
SMRT,OFF
ANTYPE,STATIC ! THERMAL ANALYSIS
ET,1,PLANE55
MP,KXX,1,1 ! 饱和状态下的渗透系数
MP,KXX,2,1E-4 ! 完全干燥下的渗透系数，假设空隙水压力小于-10M时
K,1,24,12
K,2,24,0
K,3,0,0
K,4,28,12
K,5,28,0
K,6,52,0
L,1,3
L,1,2
L,4,5
L,5,6
L,4,6
LESIZE,ALL,,,24
A,1,3,2
A,1,2,5,4
A,4,5,6
MSHK,2 ! MAPPED AREA MESH IF POSSIBLE
MSHA,0,2D ! USING QUADS
AMESH,ALL ! MESH AREAS
NUMMRG,NODE ! MERGE NODES AT BOTTOM OF CAISSON
*GET,N_MAX,NODE,,NUM,MAX ! 获得最大节点号
*GET,E_MAX,ELEM,,NUM,MAX ! 获得最大单元号
*DIM,N_TEMP,ARRAY,N_MAX ! 定义节点温度变量-总水头
*DIM,N_PRE,ARRAY,N_MAX ! 定义节点压力水头变量
!定义上游面总水头值
LSEL,S,LINE,,1
NSLL,S,1
NSEL,R,LOC,Y,0,8
D,ALL,TEMP,8 !定义上游面总水头值
!定义下游面总水头值
LSEL,S,LINE,,6
NSLL,S,1
*GET,N_NUM2,NODE,,COUNT
*DIM,N_NO2,ARRAY,N_NUM2
II=0
*DO,I,1,N_MAX
*IF,NSEL(I),EQ,1,THEN ! 判断节点是否选中
II=II+1
N_NO2(II)=I ! 存储渗流可能逸出点节点编号
*ENDIF
*ENDDO
NSEL,R,LOC,Y,0,8 ! 第一次计算，假设浸润线逸出点在8M高位置,与上游同高*GET,N_NUM,NODE,,COUNT ! 获得渗流出口节点总数
*DIM,N_NO,ARRAY,N_NUM ! 定义变量,存储渗流出口节点编号
II=0
*DO,I,1,N_MAX
*IF,NSEL(I),EQ,1,THEN ! 判断节点是否选中
N_NO(II)=I ! 存储渗流出口节点编号
*ENDIF
*ENDDO
*DO,I,1,N_NUM
D,N_NO(I),TEMP,NY(N_NO(I)) ! 定义下游面总水头值
*ENDDO
ALLSEL,ALL
FINISH
/SOLU
SOLVE
FINISH
!第一次计算完毕
!------------------------------------------------------------------------- !迭代计算
CONUTT=20 ! 最大循环次数
DD_HEAT=0.001 ! 前后两次计算，总水头最大允许计算差CHUK_ST=3 ! 出口边界条件重新设定的起始点
CHUK_MAXY2=10E5 ! 临时变量，用于存储浸润线出口坐标*DO,COM_NUM,1,CONUTT
DD_H=0
/POST1
SET,1
*DO,I,1,N_MAX
*IF,COM_NUM,GT,CHUK_ST+1,THEN
DD1=N_TEMP(I)
*IF,ABS(DD1-TEMP(I)),GT,DD_H,THEN
DD_H=ABS(DD1-TEMP(I))
*ENDIF
*ENDIF
N_TEMP(I)=TEMP(I) ! 计算节点温度（总水头）
N_PRE(I)=N_TEMP(I)-NY(I) ! 计算节点压力，总水头-Y坐标*ENDDO
*IF,COM_NUM,GT,CHUK_ST+1,THEN
*IF,DD_H,LE,DD_HEAT,THEN
*EXIT
*ENDIF
*ENDIF
/PREP7
! 重新给每个单元设定材料
MATNUM=2
*DO,I,1,E_MAX
*DO,KK,1,4
*GET,NCON(KK),ELEM,I,NODE,KK ! 获取单元四个节点编号
*ENDDO
TEMP_Y=(N_TEMP(NCON(1))+N_TEMP(NCON(2))+N_TEMP(NCON(3))+N_TEMP(NCON (4)))/4 !计算单元中心点平均温度
RESS_T=TEMP_Y-CENTRY(I)
*IF,PRESS_T,GT,0,THEN
RESS_T=0
MPCHG,1,I
*ELSEIF,PRESS_T,LT,-10,THEN
RESS_T=-10
MPCHG,2,I
*ELSE
MP,KXX,MATNUM+1,10**TPRE(PRESS_T)
MPCHG,MATNUM+1,I
MATNUM=MATNUM+1
*ENDIF
*ENDDO
! 重新设定出口边界条件
*IF,CONUTT,GT,CHUK_ST,THEN !前CHUK_ST次采用原边界条件
LSEL,S,LINE,,6
NSLL,S,1
DDELE,ALL,TEMP ! 删除原边界条件
II=0
CHUK_MAXY=0
*DO,JJ,1,N_NUM2
*IF,N_TEMP(N_NO2(JJ)),GE,NY(N_NO2(JJ)),THEN
D,N_NO2(JJ),TEMP,NY(N_NO2(JJ)) ! 总水头=Y坐标
*IF,NY(N_NO2(JJ)),GT,CHUK_MAXY,THEN
CHUK_MAXY=NY(N_NO2(JJ))
*ENDIF
*ENDIF
*ENDDO
*IF,CHUK_MAXY2,NE,CHUK_MAXY,THEN ! 判断前后两次计算的浸润线出口位置是否相同NSEL,R,LOC,Y,CHUK_MAXY ! 选择最高节点
*IF,CHUK_MAXY,GT,0,THEN
DDELE,ALL,TEMP ! 删除出口最高节点边界条件
*ENDIF
CHUK_MAXY2=CHUK_MAXY
*ENDIF
*ENDIF
ALLSEL,ALL
FINI
/SOLU
SOLVE
FINISH
*ENDDO
SAVE
!迭代计算完毕，进入后处理
FINISH
/POST1
/CLABEL,,1
/EDGE,,0
/CONTOUR,,8,0,1,8
PLNSOL,TEMP ! 显示总水头云图
PLVECT,TF, , , ,VECT,ELEM,ON,0
PLVECT,TF, , , ,VECT,NODE,ON,0
LSEL,S,LINE,,6
NSLL,S,1
PRRSOL,HEAT ! PRINT FLOWRATE THROUGH SOIL FSUM,HEAT ! 计算渗流量
*GET,Q_DAY,FSUM,0,ITEM,HEAT
ALLSEL,ALL
SAVE
*DO,I,1,N_MAX
N_TEMP(I)=TEMP(I) ! 计算节点总水头（温度）
N_PRE(I)=N_TEMP(I)-NY(I) ! 计算节点压力，总水头-Y坐标DNSOL,I,TEMP,,N_PRE(I) ! 将压力水头值复制到节点
*ENDDO
PLNSOL,TEMP ! 显示压力水头云图
FINI
ANSYS命令流实例
/PREP7 !进入前处理 ANTYPE,STATIC !设置分析类型为静力结构分析 PSTRES,ON !用于后面的模态分析中考虑预应力(该开关不影响静力分析) ET,1,LINK10 !选取单元类型1(单向杆单元) KEYOPT,1,3,0 !设置仅承受拉应力，KEYOPT(3)=0 R,1,306796E-8,543248E-8 !设置实常数，包括绳索截面积(306796E-8)，初始应变(543248E-8) MP,EX,1,30E6 !定义材料的弹性模量(1号材料) MP,DENS,1,73E-5 !定义材料的密度(1号材料) N,1 ! 定义第1号节点 N,14,100 ! 定义第14号节点 FILL ! 均分填满第2号至第13号节点 E,1,2 !由节点1及节点2生成单元 EGEN,13,1,1 !依序复制生成13个单元 D,ALL,ALL ! 对所有节点施加固定约束 FINISH ! 前处理结束 /SOLU ! 进入求解模块，求解预应力引起的应力状态 SOLVE ! 求解 FINISH ! 退出求解模块 /POST1 ! 进入一般的后处理 ETABLE,STRS,LS,1 !针对LINK10单元,建立单元列表STRS,通过LS及特征号1来获得单元的轴向应力 *GET,STRSS,ELEM,13,ETAB,STRS !针对单元列表STRS, 提取13号单元的应力 FINISH ! 后处理结束 /POST26 ! 进入时间历程后处理，处理支反力 RFORCE,2,1,F,X !将1号节点上的x方向支反力提取，并存储到2号变量中 STORE ! 存储 *GET,FORCE,V ARI,2,EXTREM,VMAX !将2号变量的最大值赋给参数FORCE /SOLU ! 再次进入求解模块，模态分析 ANTYPE,MODAL ! 模态分析 MODOPT,SUBSP,3 ! 选择子空间迭代法，求3阶模态 MXPAND,3 ! 设定3阶模态扩展 PSTRES,ON ! 用于在模态分析中考虑预应力(还需在前面的静力分析中也同时打开) DDELE,2,UX,13 ! 删除从2号节点到13号节点上的UX约束 DDELE,2,UY,13 !删除从2号节点到13号节点上的UY约束 SOLVE !求解 *GET,FREQ1,MODE,1,FREQ ! 提取第1阶模态共振频率，并赋值给参数FREQ1 *GET,FREQ2,MODE,2,FREQ ! 提取第2阶模态共振频率，并赋值给参数FREQ2 *GET,FREQ3,MODE,3,FREQ ! 提取第3阶模态共振频率，并赋值给参数FREQ3 *STATUS !列出所有参数的实际内容
ansys旋转经典命令流
1 旋转摩擦 (1) 2. 电磁三d命令流实例（论坛看到） (11) 3. 帮助感应加热例子induction heating of a solid cylinder billet (15) 4. 感应加热温度场的数值模拟（论文）inducheat30命令流 (19) 5. 如何施加恒定的角速度？Simwe仿真论坛 (24) 6. 旋转一个已经生成好的物体 (27) 7. 产生这样的磁力线 (28) 8. 旋转摩擦生热简单例子（二维旋转） (32) 8.1. 原版 (32) 8.2. 部分gui操作 (35) 9. VM229 Input Listing (39) 10 轴承---耦合+接触分析 (47) 11. 板的冲压仿真 (52) 1 旋转摩擦 FINISH /FILNAME,Exercise24 !定义隐式热分析文件名 /PREP7 !进入前处理器 ET,1,SOLID5 !选择单元类型 MPTEMP,,,,,,,, MPTEMP,1,0 MPDATA,DENS,1,,7800 !定义材料1的密度 MPTEMP,,,,,,,, MPTEMP,1,0 MPDATA,C,1,,460 !定义材料1的比热 MPTEMP,,,,,,,, MPTEMP,1,0 MPDATA,KXX,1,,66.6 !定义材料1的热传导系数 MPTEMP,,,,,,,, MPTEMP,1,0 UIMP,1,REFT,,,30 !定义材料1的热膨胀系数的参考温度 MPDATA,ALPX,1,,1.06e-5 !定义材料1的热膨胀系数MPTEMP,,,,,,,, MPTEMP,1,0 MPDATA,EX,1,,206e9 !定义材料1的弹性模量 MPDATA,PRXY,1,,0.3 !定义材料1的泊松比 MPTEMP,,,,,,,, MPTEMP,1,0 MPDATA,DENS,2,,8900 !定义材料2的密度 MPTEMP,,,,,,,, MPTEMP,1,0
Ansys常见命令流
Ansys命令流 第一天 目标:熟悉ANSYS基本关键字的含义 k --> Keypoints 关键点 l --> Lines 线 a --> Area 面 v --> Volumes 体 e --> Elements 单元 n --> Nodes 节点 cm --> component 组元 et --> element type 单元类型 mp --> material property 材料属性 r --> real constant 实常数 d --> DOF constraint 约束 f --> Force Load 集中力 sf --> Surface load on nodes 表面载荷 bf --> Body Force on Nodes 体载荷 ic --> Initial Conditions 初始条件 第二天 目标:了解命令流的整体结构,掌握每个模块的标识 !文件说明段 /BATCH /TITILE,test analysis !定义工作标题 /FILENAME,test !定义工作文件名 /PREP7 !进入前处理模块标识 !定义单元,材料属性,实常数段 ET,1,SHELL63 !指定单元类型 ET,2,SOLID45 !指定体单元 MP,EX,1,2E8 !指定弹性模量 MP,PRXY,1,0.3 !输入泊松比 MP,DENS,1,7.8E3 !输入材料密度 R,1,0.001 !指定壳单元实常数-厚度...... !建立模型 K,1,0,0,, !定义关键点 K,2,50,0,,
K,3,50,10,, K,4,10,10,, K,5,10,50,, K,6,0,50,, A,1,2,3,4,5,6, !由关键点生成面 ...... !划分网格 ESIZE,1,0, AMESH,1 ...... FINISH !前处理结束标识 /SOLU !进入求解模块标识 !施加约束和载荷 DL,5,,ALL SFL,3,PRES,1000 SFL,2,PRES,1000 ...... SOLVE !求解标识 FINISH !求解模块结束标识 /POST1 !进入通用后处理器标识 ...... /POST26 !进入时间历程后处理器 …… /EXIT,SAVE !退出并存盘 以下是日志文件中常出现的一些命令的标识说明,希望能给大家在整理LOG文件时有所帮助 /ANGLE !指定绕轴旋转视图 /DIST !说明对视图进行缩放 /DEVICE !设置图例的显示,如:风格,字体等 /REPLOT !重新显示当前图例 /RESET !恢复缺省的图形设置 /VIEW !设置观察方向 /ZOOM !对图形显示窗口的某一区域进行缩放
渗流稳定计算
赤峰市红山区城郊乡防洪工程 5.6稳定计算 5.6.1渗流及渗透稳定计算 1）渗流分析的目的 （1）确定堤身浸润线及下游逸出点位置，以便核算堤坡稳定。 （2）估算堤身、堤基的渗透量。 （3）求出局部渗流坡降，验算发生渗透变形的可能。 概括以上分析，对初步拟定的土堤剖面进行修改，最后确定土堤剖面及主渗，排水设备的型式及尺寸。 2）渗流分析计算的原则 （1）土堤渗流分析计算断面应具有代表性。 （2）土堤渗流计算应严格按照《堤防工程设计规范》（GB50286-981）第8.1.2条及本规范附录E的有关规定执行。 3）渗流分析计算的内容 （1）核算在设计洪水持续时间内浸润线的位置，当在背水侧堤坡逸出时，应计算出逸点位置，逸出段与背水侧堤基表面的出逸比降。 （2）当堤身、堤基土渗透系数K≥10-3cm/s时，应计算渗流量。 （3）设计洪水位降落时临水侧堤身内自由水位。 4）堤防渗流分析计算的水位组合 （1）临水侧为设计洪水位，背水侧为相应水位。 （2）临水侧为设计洪水位，背水侧无水。 （3）洪水降落时对临水侧堤坡稳定最不利情况。 5）渗透计算方法 堤防渗流分析计算方法按照《堤防工程设计规范》（GB50286-98）附录E3的透水堤基均质土堤渗流计算即——渗流问题的水力学解法。
6）土堤渗流分析计算 计算锡泊河左岸（0-468）横断面，堤高 5.05米（P=2%），半支箭左岸（0+302.25）横断面，堤高6.46米（P=2%），该两段堤防均属于 2级堤防，堤防渗流计算断面采用1个断面计算即可。采用《堤防工程设计规范》中透水堤基均质土堤下游坡无排水设备或有贴坡式排水稳定渗流计算公式： T H L T H H D 88.0m k q q 11210 ++-+=）（ （E.3.1） H m m b 121+-+=）（H H L （E2.1-3） 111 1 2m m H L += ? （E2.1-4） 当K ≤k 0时 h 0=a+H 2=q ÷? ???? ?+++??????++++?T H a m T K H a m H m m K 44.0)(5.0)5.0()5.0(1220222 22 +H 2 ……………（E.3.2-2） 对于各种情况下坝体浸润线均可按下式确定 X=k ·T '0q h y -+k ' 22 2q h y - ……………（E.3.2-6） 式中：q'= ） （021112 0211 m 2m 2k h m H L h H -++-+0 2110 10m k h m H L h H T -+-（E.3.2-7） k ——堤身渗透系数； k 0——堤基渗透系数； H 1——水位到坝脚的距离（m ）； H 2——下游水位（m ）； H ——堤防高度（m ）； q ——单位宽度渗流量（m 3/s ·m ）； m 1——上游坡坡率，m 1=3.0；
ansys实例命令流-弹塑性分析命令流
/FILNAME,Elastic-Plasitc,1 /TITLE, Elastic-Plasitc Analysis !前处理。 /PREP7 !**定义梁单元189。 ET,1,BEAM189 !定义单元。 !**梁截面1。 SECTYPE, 1, BEAM, HREC, , 0 !定义梁截面。SECOFFSET, CENT SECDATA,50,100,6,6,6,6,0,0,0,0 !定义梁截面完成。 !**定义材料。 MPTEMP,,,,,,,, !定义弹塑性材料模型。MPTEMP,1,0 MPDATA,EX,1,,2.05e5 MPDATA,PRXY,1,,0.3 TB,BISO,1,1,2, TBTEMP,0 TBDATA,,150,18600,,,, !定义弹塑性材料模型。!**建立几何模型。 K,1, , , , K,2 ,900, K,3 ,,50 LSTR, 1, 2 !**网格划分。 FLST,5,1,4,ORDE,1 !定义网格密度。FITEM,5,1 CM,_Y,LINE LSEL, , , ,P51X CM,_Y1,LINE CMSEL,,_Y LESIZE,_Y1, , ,50, , , , ,1 !定义网格密度完成。CM,_Y,LINE !网格划分。 LSEL, , , , 1 CM,_Y1,LINE CMSEL,S,_Y CMSEL,S,_Y1 LATT,1, ,1, , 3, ,1 CMSEL,S,_Y CMDELE,_Y CMDELE,_Y1 LMESH, 1 !网格划分完成。 !施加载荷及求解。 FINISH /SOL
!**施加约束。 FLST,2,1,3,ORDE,1 !施加约束。FITEM,2,1 /GO DK,P51X, , , ,0,UX,UY,UZ,ROTX, , , FLST,2,1,3,ORDE,1 FITEM,2,2 /GO DK,P51X, , , ,0,UY,UZ,ROTX, , , , !施加约束完成。 !**加载。 FLST,2,50,2,ORDE,2 FITEM,2,1 FITEM,2,-50 SFBEAM,P51X,1,PRES,100, , , , , , LSWRITE,1, !定义载荷步1完成。FLST,2,50,2,ORDE,2 !定义载荷步2。FITEM,2,1 FITEM,2,-50 SFEDELE,P51X,1,PRES LSWRITE,2, !定义载荷步2完成。!设定求解步并求解。 LSSOLVE,1,2,1,
ANSYS-结构稳态(静力)分析之经典实例-命令流格式
ANSYS 结构稳态(静力)分析之经典实例-命令流格式.txt两人之间的感情就像织毛衣，建立 的时候一针一线，小心而漫长，拆除的时候只要轻轻一拉。。。。/FILNAME,Allen-wrench,1 ! Jobname to use for all subsequent files /TITLE,Static analysis of an Allen wrench /UNITS,SI ! Reminder that the SI system of units is used /SHOW ! Specify graphics driver for interactive run; for batch ! run plots are written to pm02.grph ! Define parameters for future use EXX=2.07E11 ! Young's modulus (2.07E11 Pa = 30E6 psi) W_HEX=.01 ! Width of hex across flats (.01m=.39in) *AFUN,DEG ! Units for angular parametric functions定义弧度单位 W_FLAT=W_HEX*TAN(30) ! Width of flat L_SHANK=.075 ! Length of shank (short end) (.075m=3.0in) L_HANDLE=.2 ! Length of handle (long end) (.2m=7.9 in) BENDRAD=.01 ! Bend radius of Allen wrench (.01m=.39 in) L_ELEM=.0075 ! Element length (.0075 m = .30 in) NO_D_HEX=2 ! Number of divisions on hex flat TOL=25E-6 ! Tolerance for selecting nodes (25e-6 m = .001 in) /PREP7 ET,1,SOLID45 ! 3维实体结构单元；Eight-node brick element ET,2,PLANE42 ! 2维平面结构；Four-node quadrilateral (for area mesh) MP,EX,1,EXX ! Young's modulus for material 1；杨氏模量 MP,PRXY,1,0.3 ! Poisson's ratio for material 1；泊松比 RPOLY,6,W_FLAT ! Hexagonal area创建规则的多边形 K,7 ! Keypoint at (0,0,0) K,8,,,-L_SHANK ! Keypoint at shank-handle intersection K,9,,L_HANDLE,-L_SHANK ! Keypoint at end of handle L,4,1 ! Line through middle of hex shape L,7,8 ! Line along middle of shank L,8,9 ! Line along handle LFILLT,8,9,BENDRAD ! Line along bend radius between shank and handle! 产生 一个倒角圆，并生成三个点 /VIEW,,1,1,1 ! Isometric view in window 1 /ANGLE,,90,XM ! Rotates model 90 degrees about X! 不用累积的旋转 /TRIAD,ltop /PNUM,LINE,1 ! Line numbers turned on LPLOT
ansys命令流
第一天目标： 熟悉ANSYS基本关键字的含义k --> Keypoints关键点l --> Lines线a --> Area 面v --> Volumes体e --> Elements单元n --> Nodes节点cm --> component组元et --> element type单元类型mp --> material property材料属性r --> real constant实常数d --> DOF constraint约束f --> Force Load集中力sf --> Surface load on nodes 表面载荷bf --> Body Force on Nodes体载荷ic --> Initial Conditions初始条件第二天目标： 了解命令流的整体结构,掌握每个模块的标识!文件说明段/BATCH/TILE,test analysis!定义工作标题/FILENAME,test!定义工作文件名/PREP7!进入前处理模块标识!定义单元,材料属性,实常数段ET,1,SHELL63!指定单元类型ET,2,SOLID45!指定体单元MP,EX,1,2E8!指定弹性模量MP,PRXY,1, 0.3!输入泊松比MP,DENS,1, 7.8E3!输入材料密度R,1, 0.001!指定壳单元实常数-厚度......!建立模型K,1,0,0,,!定义关键点 K,2,50,0,,K,3,50,10,,K,4,10,10,,K,5,10,50,,K,6,0,50,,A,1,2,3,4,5,6,!由关键点生成面......!划分网格ESIZE,1,0,AMESH, 1......FINISH!前处理结束标识/SOLU!进入求解模块标识!施加约束和载荷DL,5,,ALLSFL,3,PRES,1000SFL,2,PRES, 1000......SOLVE!求解标识FINISH!求解模块结束标识/POST1!进入通用后处理器标识....../POST26!进入时间历程后处理器……/EXIT,SAVE!退出并存盘以下是日志文件中常出现的一些命令的标识说明,希望能给大家在整理LOG文件时有所帮助/ANGLE!指定绕轴旋转视图/DIST!说明对视图进行缩放/DEVICE!设置图例的显示,如： 风格,字体等/REPLOT!重新显示当前图例/RESET!恢复缺省的图形设置/VIEW!设置观察方向/ZOOM!对图形显示窗口的某一区域进行缩放第三天生成关键点和线部分 1.生成关键点K,关键点编号,X坐标,Y坐标,Z坐标例：
ansys实例命令流-谱分析命令流
/FILNAME, Beam,1 ！定义工作文件名。 /TITLE, Beam Analysis ！定义工作标题。/PREP7 !定义单元。 ET,1,BEAM188 !定义材料属性。 MPTEMP,,,,,,,, MPTEMP,1,0 MPDATA,EX,1,,2.1e5 MPDATA,PRXY,1,,0.3 MPTEMP,,,,,,,, MPTEMP,1,0 MPDATA,DENS,1,,7.9e-6 ! 定义杆件截面■200。 SECTYPE, 1, BEAM, RECT, , 0 SECOFFSET, CENT SECDATA,10,10,0,0,0,0,0,0,0,0 !建立几何模型。 K,1, ,, , K,2,350,, , !生成立柱。 LSTR, 1, 2 !以上完成几何模型。 !以下进行网格划分。 FLST,5,1,4,ORDE,1 FITEM,5,1 CM,_Y,LINE LSEL, , , ,P51X CM,_Y1,LINE CMSEL,,_Y LESIZE,_Y1, , ,50, , , , ,1 !定义单元大小。!分配、划分平板结构。 LMESH, 1 !分析类型施加载荷并求解。 FINISH /SOLUTION ANTYPE,2 ！定义分析类型及求解设置。MSAVE,0 ！模态提取方法。
MODOPT,LANB,10 EQSLV,SPAR MXPAND,10, , ,0 ！模态扩展设置。 LUMPM,0 PSTRES,0 MODOPT,LANB,10,0,0, ,OFF MXPAND,10,0,0,1,0.001, !施加约束。 FLST,2,2,3,ORDE,2 FITEM,2,1 FITEM,2,-2 /GO DK,P51X, , , ,0,ALL, , , , , , !求解。 FINISH /SOL /STATUS,SOLU SOLVE !定义谱分析。 FINISH /SOLUTION ANTYPE,8 SPOPT,PSD,10,1 PSDUNIT,1,DISP,386.4, PSDFRQ,1, ,13.8,40,50.6,73,120 !定义谱—频率表。PSDFRQ,1, ,134,178,233, , PSDV AL,1,1,4,0.6,3,5 PSDV AL,1,6,2,6, , FLST,2,2,1,ORDE,2 !施加谱。 FITEM,2,1 FITEM,2,-2 D,P51X,UX,1.0 PFACT,1,BASE, !计算PSD激励参与系数。PSDRES,DISP,REL !设置输出选项。PSDRES,VELO,OFF PSDRES,ACEL,OFF
几个ansys经典实例(长见识)
平面问题斜支座的处理 如图5-7所示，为一个带斜支座的平面应力结构，其中位置2及3处为固定约束，位置4处为一个45o的斜支座，试用一个4节点矩形单元分析该结构的位移场。 (a)平面结构(b)有限元分析模型 图5-7 带斜支座的平面结构 基于ANSYS平台，分别采用约束方程以及局部坐标系的斜支座约束这两种方式来进行处理。 (7) 模型加约束 左边施加X,Y方向的位移约束 ANSYS Main Menu: Solution →Define Loads →Apply →-Structural→Displacement On Nodes →选取2，3号节点→OK →Lab2: All DOF(施加X,Y方向的位移约束) →OK 以下提供两种方法处理斜支座问题，使用时选择一种方法。 ?采用约束方程来处理斜支座 ANSYS Main Menu：Preprocessor →Coupling/ Ceqn →Constraint Eqn ：Const :0, NODE1:4, Lab1: UX,C1:1,NODE2:4,Lab2:UY,C2:1→OK 或者?采用斜支座的局部坐标来施加位移约束 ANSYS Utility Menu：WorkPlane →Local Coordinate System →Create local system →At specified LOC + →单击图形中的任意一点→OK →XC、YC、ZC分别设定为2，0，0,THXY：45 →OK ANSYS Main Menu：Preprocessor →modeling →Move / Modify →Rotate Node CS →To active CS → 选择4号节点 ANSYS Main Menu：Solution →Define Loads →Apply →Structural →Displacement On Nodes →选取4号节点→OK →选择Lab2:UY(施加Y方向的位移约束) →OK 命令流； !---方法1 begin----以下的一条命令为采用约束方程的方式对斜支座进行处理 CE,1,0,4,UX,1,4,UY,-1 !建立约束方程(No.1): 0=node4_UX*1+node_UY*(-1) !---方法1 end --- !--- 方法2 begin --以下三条命令为定义局部坐标系，进行旋转，施加位移约束 !local,11,0,2,0,0,45 !在4号节点建立局部坐标系 !nrotat, 4 !将4号节点坐标系旋转为与局部坐标系相同 !D,4,UY !在局部坐标下添加位移约束 !--- 方法2 end
ansys命令流最全详细介绍二
三 生成关键点和线部分 1.生成关键点 K,关键点编号,X坐标,Y坐标,Z坐标 例:K,1,0,0,0 2.在激活坐标系生成直线 LSTR,关键点P1,关键点P2 例LSTR,1,2 3.在两个关键点之间连线 L,关键点P1,关键点P2 例L,1,2 注:此命令会随当前的激活坐标系不同而生成直线或弧线 4.由三个关键点生成弧线 LARC,关键点P1,关键点P2,关键点PC,半径RAD 例LARC,1,3,2,0.05 注:关键点PC是用来控制弧线的凹向 5.通过圆心半径生成圆弧
CIRCLE,关键点圆心,半径RAD,,,,圆弧段数NSEG 例:CIRCLE,1,0.05,,,,4 6.通过关键点生成样条线 BSPLIN,关键点P1,关键点P2,关键点P3,关键点P4,关键点P5,关键点P6 例:BSPLIN,1,2,3,4,5,6 7.生成倒角线 LFILLT,线NL1,线NL2,倒角半径RAD 例LFILLT,1,2,0.005 8.通过关键点生成面 A,关键点P1,关键点P2,关键点P3,关键点P4,关键点P5,关键点P6,P7,P8... 例:A,1,2,3,4 9.通过线生成面 AL,线L1,线L2,线L3,线L4,线L5,线L6,线L7,线L8,线L9,线L10 例:AL,5,6,7,8 10.通过线的滑移生成面
ASKIN,线NL1,线NL2,线NL3,线NL4,线NL5,线NL6,线NL7,线NL8,线NL9 例:ASKIN,1,4,5,6,7,8 注:线1为滑移的导向线 四 目标:掌握常用的实体-面的生成 生成矩形面 1.通过矩形角上定位点生成面 BLC4,定位点X方向坐标XCORNER,定位点Y方向坐标YCORNER,矩形宽度WIDTH,矩形高度HEIGHT,矩形深度DEPTH 例:BLC4,0,0,5,3,0 2.通过矩形中心定位点生成面 BLC5,定位点X方向坐标XCENTER,定位点Y方向坐标YCENTER,矩形宽度WIDTH,矩形高度HEIGHT,矩形深度DEPTH 注:与上条命令的不同就在于矩形的定位点不一样 例:BLC5,2.5,1.5,5,3,0 3.通过在工作平面定义矩形X.Y坐标生成面 RECTNG,矩形左边界X坐标X1,矩形右边界X坐标X2,矩形下边界Y
ANSYS建模apdl命令流实例应用
大桥全长2996.8m，其中主桥采用跨度为101.5+188.5+580+217.5+159.5+116m的钢桁梁斜拉桥；非通航孔正桥采用6孔跨径64m预应力混凝土简支箱梁；东引桥采用16孔梁长32.6m预应力混凝土简支箱梁；跨大堤桥采用48.9+86+48.8m预应力混凝土连续箱梁；西引桥采用15孔梁长32.6m预应力混凝土简支梁及2孔梁长24.6m预应力混凝土简支梁，其中宁安线采用箱梁，阜景线采用T梁。 主桥采用103+188.5+580+217.5+159.5+117.5m两塔钢桁斜拉桥方案，全长1366m。主梁为三片主桁钢桁梁，桁间距2x14m，节间长14.5m，桁高15m。主塔为钢筋混凝土结构，塔顶高程+204.00m，塔底高程-6.00m，斜拉索为空间三索面，立面上每塔两侧共18对索，全桥216根斜拉索。所有桥墩上均设竖向和横向约束，4#塔与主梁之间设纵向水平约束，3#塔与梁间使用带限位功能的粘滞阻尼器。主梁为”N”字型桁式，横向采用三片桁结构，主桁的横向中心距各为14m，桁高15m，节间距14.5m[2]。 结构构造 主桥采用两塔钢桁斜拉桥方案，主梁为三片主桁钢桁梁，主桁上下弦杆均为箱型截面，上弦杆内高1000mm，内宽1200mm，板厚20～48mm。下弦杆内高1400mm，宽1200mm，板厚20～56mm。下弦杆顶板向桁内侧加宽700mm与整体桥面板焊接。腹杆主要采用H型截面。H型杆件宽1200mm，高720和760mm，板厚20～48mm。根据不同的受力区段选用不同的杆件截面，在辅助墩附近的压重区梁段，腹杆采用箱型截面杆件。主桁采用焊接杆件，整体节点。在节点外以高强度螺栓拼接的结构形式，上下弦杆四面等强对接拼装。H型腹杆采用插入式连接。箱型腹杆采用四面与主桁节点对拼的连接形式。主桁拼接采用M30高强螺栓。
土石坝渗流安全评价
土石坝渗流安全评价 Revised by Hanlin on 10 January 2021
土石坝渗流安全评价1坝基渗流安全评价要点如下： 1砂砾石层（包括砂层、砂砾石层、砾卵石层等）的渗透稳定性，应根据土的类型及其颗粒级配等情况判别其渗透变形形式，核定其相应的允许渗透比降，与工程实际渗透比降相比，判断渗流出口有无管涌或流土破坏的可能性，以及渗流场内部有无管涌、接触冲刷等渗流隐患。 2覆盖层为相对弱透水土层时，应复核其抗浮动稳定性，其允许渗透比降宜由试验法或参考流土指标确定；对已有反滤盖重者，应核算盖重厚度和范围是否满足要求。 3接触面的渗透稳定性主要有以下两种型式： 1）复核粗、细散粒料土层之间有无接触冲刷（流向平行界面）和接触流土（流向从细到粗垂直界面）的可能性；粗粒料层能否对细粒料层起保护作用。 2）复核散粒料土体与刚性结构物体（如混凝土墙、涵管和岩石等）界面的接触渗透稳定性。应注意散粒料与刚性面结合的紧密程度、出口有无
反滤保护，以及与断层破碎带、灰岩溶蚀带、较大张性裂隙等接触面有无妥善处理及其抗渗稳定性。 2坝体渗流安全评价要点如下： 1均质坝。复核坝体的防渗性能是否满足规范要求、坝体实际浸润线和下游坝坡渗出段高程是否高于设计值，还需注意坝内有无横向或水平裂缝、松软结合带或渗漏通道等。 2组合（分区）坝： 1）防渗体（心墙、斜墙、铺盖、各种面板等）。复核防渗体的防渗性能是否满足规范要求，心墙或斜墙的上、下游侧有无合格的过渡保护层，以及水平防渗铺盖的底部垫层或天然砂砾石层能否起保护作用。 2）透水区（上、下游坝壳及各类排水体等）。复核上、游坝坡在库水骤降情况下的抗滑稳定性和下游坝坡出逸段（区）的渗透稳定性，下游坡渗出段的贴坡保护层应满足反滤层的设计要求。 3）过渡区。界于坝体粗、细填料之间的过渡区以及棱体排水、褥垫排水和贴坡排水等，应复核反滤层设计的保土条件和排水条件是否合格，以及运行中有无明显集中渗流和大量固体颗粒被带出等异常现象。
ansys实例命令流-实体梁分析命令流
/FILNAME,SolidBeam ,1 !定义工作文件名。/TITLE,SolidBeam Analysis !定义工作标题。/PREP7 ET,1,SOLID95 !定义材料属性。 MPTEMP,,,,,,,, !定义材料属性。MPTEMP,1,0 MPDATA,EX,1,,2.06e5 MPDATA,PRXY,1,,0.3 !建立几何模型 K,1,,,, K,2,450,,, K,3,450,-55,, K,4,,-100,, FLST,2,4,3 FITEM,2,1 FITEM,2,2 FITEM,2,3 FITEM,2,4 A,P51X VOFFST,1,45, , !网格划分。 FLST,5,4,4,ORDE,2 FITEM,5,9 FITEM,5,-12 CM,_Y,LINE LSEL, , , ,P51X CM,_Y1,LINE CMSEL,,_Y LESIZE,_Y1, , ,3, , , , ,1 !* FLST,5,4,4,ORDE,4
FITEM,5,2 FITEM,5,4 FITEM,5,6 FITEM,5,8 CM,_Y,LINE LSEL, , , ,P51X CM,_Y1,LINE CMSEL,,_Y LESIZE,_Y1, , ,6, , , , ,1 !* FLST,5,4,4,ORDE,4 FITEM,5,1 FITEM,5,3 FITEM,5,5 FITEM,5,7 CM,_Y,LINE LSEL, , , ,P51X CM,_Y1,LINE CMSEL,,_Y LESIZE,_Y1, , ,30, , , , ,1 !* CM,_Y,VOLU VSEL, , , , 1 CM,_Y1,VOLU CHKMSH,'VOLU' CMSEL,S,_Y !* VSWEEP,_Y1 CMDELE,_Y CMDELE,_Y1 CMDELE,_Y2 !加载。
ANSYS_结构稳态(静力)分析之经典实例-命令流格式
/FILNAME,Allen-wrench,1 ! Jobname to use for all subsequent files /TITLE,Static analysis of an Allen wrench /UNITS,SI ! Reminder that the SI system of units is used /SHOW ! Specify graphics driver for interactive run; for batch ! run plots are written to pm02.grph ! Define parameters for future use EXX=2.07E11 ! Young's modulus (2.07E11 Pa = 30E6 psi) W_HEX=.01 ! Width of hex across flats (.01m=.39in) *AFUN,DEG ! Units for angular parametric functions定义弧度单位 W_FLAT=W_HEX*TAN(30) ! Width of flat L_SHANK=.075 ! Length of shank (short end) (.075m=3.0in) L_HANDLE=.2 ! Length of handle (long end) (.2m=7.9 in) BENDRAD=.01 ! Bend radius of Allen wrench (.01m=.39 in) L_ELEM=.0075 ! Element length (.0075 m = .30 in) NO_D_HEX=2 ! Number of divisions on hex flat TOL=25E-6 ! Tolerance for selecting nodes (25e-6 m = .001 in) /PREP7 ET,1,SOLID45 ! 3维实体结构单元；Eight-node brick element ET,2,PLANE42 ! 2维平面结构；Four-node quadrilateral (for area mesh) MP,EX,1,EXX ! Young's modulus for material 1；杨氏模量 MP,PRXY,1,0.3 ! Poisson's ratio for material 1；泊松比 RPOLY,6,W_FLAT ! Hexagonal area创建规则的多边形 K,7 ! Keypoint at (0,0,0) K,8,,,-L_SHANK ! Keypoint at shank-handle intersection K,9,,L_HANDLE,-L_SHANK ! Keypoint at end of handle L,4,1 ! Line through middle of hex shape L,7,8 ! Line along middle of shank L,8,9 ! Line along handle LFILLT,8,9,BENDRAD ! Line along bend radius between shank and handle! 产生一个倒角圆，并生成三个点 /VIEW,,1,1,1 ! Isometric view in window 1 /ANGLE,,90,XM ! Rotates model 90 degrees about X! 不用累积的旋转 /TRIAD,ltop /PNUM,LINE,1 ! Line numbers turned on LPLOT ! Line numbers off !
理正软件计算土石坝渗流稳定
用理正软件计算土石坝渗流稳定的方法 1渗流计算 1在CAD中绘制土石坝横断面图，图中坝坡下的长垫层为基岩，图例中有两种基岩，根据情况有几种画几种，长度为1.5-2倍坝长，注意不能使用镜像。 绘制时要注意 并另存为DXF文件（最好存为最低版本即2000） 2进行渗流计算 打开理正岩土软件，选择渗流分析计算 在选工程中选择软件生成结论的存储位置 如上例，计算结论存在e盘考博文件中，确认后弹出下图直接点确认即可。 确认后点增，选择系统默认例题，点确认 然后自动弹出下图中对话框 然后点击左上角的“辅助功能”选择“读入DXF文件自动生成坡面、节点、土层数据”，弹出以下对话框 选择已画好的CAD图打开 打开后出现如下对话框，在图上双击后可放大图形，放大后可看到起始点编号（起始点在图中用红圈标出，及上游坝坡起始点）。坡面线段数及坝坡分为几段，无马道土石坝坡面线段数为3，图例中有9条。 弹出以下对话框，在坡面形状中填写正确的上下游水位 节点坐标一栏为理正自动生成坐标，不用修改 土层定义一栏如下图，图中不同土性区域数为软件自动生成软件同时为不同区域编号，双击图中土石坝图形放大图形可以看到编号（如下附图2）Kx，Ky为土层的x，y向的渗透系数，同一土层两数相等且等于土层渗透系数，对应区号输入渗透系数（渗透系数由地质资料中查找）α值若无资料则都为0计算即可。
附图2 面边界条件中，同样双击放大土石坝剖面图可以看到节点编号，顺时针输入计算所需要的坡面信息（即始末节点编号），面边界个数及浸润线可能经过的面，即上游所有水面线以下的坡面加上坝基上表面，下游所有坡面加坝基上表面，如图，蓝色为已知水面线，红色为可能的浸出面. 点边界描述项数为2，节点即上下游水面线与坝体的交点，若下游无水则为下游坝脚，取值为0。 计算参数栏为系统默认，不用修改 输出结果栏目中，需注意流量计算截面的点数一栏和理正边坡文件接口一栏。 流量计算截面的点数即下游截面所有点和基岩上表面所有点，如本例有5个，且须在右边一栏输入5个节点的坐标，坐标从第二栏节点坐标中查找。 理正边坡文件接口一栏一定不要忘记填入合理命名，如正常蓄水位、校核洪水位等，否则无法进行稳定计算。
命令流ansys经典实例
ansys钢筋混凝土建模实例 finish /clear /prep7 et,1,solid65,,,,,,,1 et,2,link8 et,3,185 et,4,solid45 !************************定义材料属性***************************** !混凝土材料属性 mp,ex,1,1.596e10 mp,prxy,1,0.2 mp,dens,1,2400 fc=1.68e7 !c30混凝土轴心抗压强度设计值 ft=1.86e6 tb,concr,1 tbdata,,0.5,0.95,ft,-1 tb,miso,1,,11 tbpt,,0.0002,fc*0.19 tbpt,,0.0004,fc*0.36 tbpt,,0.0006,fc*0.51 tbpt,,0.0008,fc*0.64 tbpt,,0.0010,fc*0.75 tbpt,,0.0012,fc*0.84 tbpt,,0.0014,fc*0.91 tbpt,,0.0016,fc*0.96 tbpt,,0.0018,fc*0.99 tbpt,,0.002,fc tbpt,,0.0033,fc !57到68 v,6,7,56,54,10,11,60,58 vgen,4,57,57,1,0,0.38,0 vgen,2,57,60,1,0.25,0,0 vgen,2,57,60,1,-0.25,0,0 !69到80 v,189,190,234,233,191,192,238,237 vgen,4,69,69,1,0,0.38,0
vgen,3,69,72,1,0,0,0.25 !121到136 vgen,2,105,120,1,0,0,0.107 !193到208 k,1178,0.81,0,0.107 k,1179,0.823,0,0.107 k,1180,0.81,0,0.417 k,1181,0.81,0.16,0.107 k,1182,0.823,0.16,0.107 k,1183,0.81,0.16,0.417 k,1184,0.81,0.38,0.107 k,1185,0.823,0.38,0.107 k,1186,0.81,0.38,0.417 v,1178,1179,187,1180,1181,1182,189,1183 v,1181,1182,189,1183,1184,1185,191,1186 vgen,4,193,194,1,0,0.38,0 vgen,2,193,200,1,0.107,0,0 !209到232 k,1271,-0.02,0,-0.013 k,1272,-0.02,0,0 k,1273,-0.02,0.16,-0.013 k,1274,-0.02,0.16,0 k,1275,-0.02,0.38,-0.013 k,1276,-0.02,0.38,0 v,677,1,1272,1271,679,5,1274,1273 v,679,5,1274,1273,681,9,1276,1275 vgen,4,209,210,1,0,0.38,0 vgen,3,209,216,1,0.25,0,0 !233到256 vgen,2,209,232,1,0,0,0.107 vgen,2,257,280,1,0.107,0,0 k,5001,-0.1,1.3,0 k,5002,-0.02,1.3,0 k,5003,-0.02,1.3,0.094
实例ANSYS
ANSYS应用实例—圆轴扭转分析 摘要：本文建立了圆轴的模型，并对其施加了扭转应力，应用ANSYS软件，对圆轴的扭转进行了分析，分析表明，ANSYS的分析结果与实际理论计算的结果是十分接近的。 关键词：圆轴扭转ANSYS 问题描述:设等直圆轴截面直径D=50mm，长度120mm，作用在圆轴两端上的转矩M=1.5X103N.m，试分析该圆轴的扭转情况。 分析步骤： (1)改变工作名：拾取菜单Utility Menu→File→Change Jobname，在弹出的文本框中输入适当的文件名，单击OK即可。 (2) 创建单元类型：拾取菜单Utility Menu→Preprocessor→Element Type→ Add/Edit/Delete。在弹出的对话框中单击Add按钮，在左侧列表中选Structural Solid，在右侧列表中选Quad 8node 183，单击Apply按钮，再在右侧列表中选择Brick 20node 186，单击OK按钮，然后单击对话框中的Close按钮。 图1 单元类型对话框 （3 ）定义材料类型：拾取菜单Utility Menu→Preprocessor→Material Props→Material Model。弹出如图2所示的对话框，在右侧列表中一次双击Structural、Linear、Elastic、Isotropic，弹出如图3所示的对话框，在EX文本框中输入弹性模量2.08e11，在PRXY中输入泊松比0.3，单击OK，关闭对话框即可。 图3 材料特性对话框
（4 ）创建矩形面：拾取菜单Utility Menu→Preprocessor→Modeling→Create→Areas→Rectangle→By Dimensions。弹出如图4所示的对话框，在“X1，X2“文本框中输入0,0.025，在”Y1，Y2”文本框中输入0，0.12，单击OK按钮即可。 图4 创建矩形面对话框 （5 ）划分单元：拾取菜单Utility Menu→Preprocessor→Meshing→Mesh Tool。弹出如图6所示的对话框，单击Size Controls区域中的Lines后的Set按钮，弹出拾取窗口，拾取矩形面的任意短边，单击OK，在如图7中的弹出对话框在NDIV 文本框中输入5，单击Apply按钮，在此弹出拾取窗口，拾取矩形面的任一长边，单击OK按钮，再次弹出对话框，在NDIV文本框中输入8，单击OK按钮。在Mesh区域，选择单元形状为Quad，选择划分单元的方法为Mapped，单击Mesh 按钮，弹出拾取窗口，拾取面，单击OK，然后关闭窗口即可。 图5 划分单元工具对话框