ANSYS刚桁架桥模态分析命令流
基于 Ansys 的钢桁架桥静力和模态分析
土木结构分析专题陈晨20104336基于Ansys的钢桁架桥静力和模态分析陈晨20104336(西南交通大学力学与工程学院结构2010-01班,四川成都)摘要:本文应用Ansys软件,采用有限元分析技术及其优化技术,分别采用GUI方式和命令流方式,对给定的一架钢桁架简支梁桥进行了静力学分析和模态分析,对强度、内力分布及前六届振型状况进行了查看。
关键词:力学;土木工程;桥梁工程;结构分析1设计概况图1钢桁架桥简图已知下承式简支钢桁架桥桥长72米,每个节段12米,桥宽10米,高16米。
设桥面板为0.3米厚的混凝土板。
桁架杆件规格有三种,见下表:表1钢桁架桥杆件规格杆件截面号形状规格端斜杆1工字形400×400×16×16上下弦2工字形400×400×12×12横向连接梁2工字形400×400×12×12其他腹杆3工字形400×300×12×12所用材料属性如下表:表2材料属性参数钢材混凝土弹性模量EX 2.1×1011 3.5×1010泊松比PRXY0.30.1667密度DENS785025002建立有限元模型2.1定义单元类型和选项Main Menu>Preprocessor>Element Type>Add/Edit/Delete,弹出“Element Types”选择“Structural Beam—3D elastic4”,单击“Ok”,定义“BEAM4”单元,如图6-17。
继续单击“Add”按钮,选择“Structural Shell—Elastic4node63”,定义“SHELL63”单元。
得到如图6-18所示的结果。
最后单击“Close”,关闭单元类型对话框。
图2单元类型对话框2.2定义梁单元截面Main Menu>Preprocessor>Sections Beam>Common Sections,弹出“Beam Tool”工具条,如图6-19填写。
(整理)基于ansys的钢桁架桥的分析和计算
基于ansys的钢桁架桥的分析和计算姓名: 马彦学院:建筑与环境专业:工程力学学号:1043055033指导老师:朱哲明2013/6/151.问题简述钢桁架桥简图如下,尺寸如图,单元长12m,高16m。
设桥面板为0.3m厚的混凝土板。
杆件截面号形状规格端斜杆 1 工字梁400*400*16*16上下弦 2 工字梁400*400*12*12横向连接梁 2 工字梁400*400*12*12其他腹杆 3 工字梁400*300*12*12参数钢材混凝土EX 2.1x1011 3.5x1010PRXY 0.3 0.1667DENS 7850 25002.材料实常数3.半横架桥模型镜面对称,生成整体模型3.施加约束及受力4.计算及分析结果◆整体位移云图◆结点总位移矢量图◆单元第一主应力云图◆单元第二主应力云图◆单元第三主应力云图◆节点位移结果PRINT U NODAL SOLUTION PER NODE***** POST1 NODAL DEGREE OF FREEDOM LISTING *****LOAD STEP= 1 SUBSTEP= 1 TIME= 1.0000 LOAD CASE= 0THE FOLLOWING DEGREE OF FREEDOM RESULTS ARE IN THE GLOBAL COORDINATESYSTEMNODE UX UY UZ USUM1 0.18808E-02-0.20919E-01 0.70316E-03 0.21015E-012 0.11411E-02-0.21354E-01 0.59772E-03 0.21393E-013 0.14813E-02-0.20809E-01 0.11202E-02 0.20892E-014 0.15919E-02-0.20373E-01 0.11392E-02 0.20467E-015 0.22549E-02-0.18918E-01 0.10528E-02 0.19081E-016 0.23458E-02-0.18310E-01 0.10055E-02 0.18487E-017 -0.10050E-02-0.18459E-01-0.38731E-02 0.18887E-018 -0.11376E-02-0.19066E-01-0.38598E-02 0.19486E-019 0.24977E-02-0.12074E-01 0.72603E-03 0.12351E-0110 0.29237E-02-0.11079E-01 0.68719E-03 0.11479E-0111 -0.35033E-02-0.10438E-01-0.84626E-02 0.13887E-0112 -0.38537E-02-0.10965E-01-0.84226E-02 0.14353E-0113 0.27521E-02 0.0000 0.0000 0.27521E-0214 0.34768E-02 0.0000 0.0000 0.34768E-0215 0.82671E-03-0.17947E-01 0.14911E-03 0.17967E-0116 0.67748E-03-0.19250E-01 0.10648E-03 0.19262E-0117 0.42077E-02-0.19398E-01 0.59595E-02 0.20725E-0118 0.40812E-02-0.18095E-01 0.59727E-02 0.19488E-0119 0.40101E-03-0.10784E-01 0.34385E-04 0.10791E-0120 0.34470E-03-0.12307E-01 0.25523E-06 0.12312E-0121 0.69212E-02-0.11199E-01 0.10204E-01 0.16656E-0122 0.65820E-02-0.10142E-01 0.10244E-01 0.15847E-0123 0.0000 0.0000 0.0000 0.000024 0.0000 0.0000 0.0000 0.0000MAXIMUM ABSOLUTE VALUESNODE 21 2 22 2VALUE 0.69212E-02-0.21354E-01 0.10244E-01 0.21393E-01◆单元受力结果PRINT ELEMENT TABLE ITEMS PER ELEMENT***** POST1 ELEMENT TABLE LISTING *****STAT CURRENT CURRENTELEM ZHOU_I ZHOU_J1 -49659. 7936.32 -42695. -3502.73 -9873.9 -28642.4 9567.9 -51440.5 -15016. 23374.6 -22120. -5510.47 -26981. -11385.8 -33355. 18549.9 -17656. -15556.10 -16095. -16301.11 -16203. -16943.12 -12683. -20132.13 4836.6 5157.114 -17901. -18351.15 -2331.6 23001.16 -18331. -20015.17 -6067.9 50464.18 -19568. -26493.19 -5052.8 51411.20 -26836. -34142.21 -23626. -29919.22 -32522. -21349.23 -35649. -25215.24 -699.47 1061.525 690.13 -1048.326 5802.4 -1462.327 -9677.8 5182.928 16212. -4765.129 -4310.8 3979.130 -25.038 0.000031 -9.3064 0.000032 23.898 0.000033 -3569.2 -42609.34 8110.9 -49823.35 -5544.6 -22051.36 -11343. -27005.37 18453. -33238.38 -28592. -9977.139 -51593. 9648.540 23614. -15193.41 -16998. -16116.***** POST1 ELEMENT TABLE LISTING *****STAT CURRENT CURRENTELEM ZHOU_I ZHOU_J42 -20120. -12682.43 -15489. -17761.44 -16350. -16082.45 5157.1 4836.646 -18351. -17901.47 -2225.2 22850.48 -18463. -19869.49 -6087.5 50530.50 -19228. -26843.51 -5332.4 51796.52 -21374. -32473.53 -25205. -35655.54 -34114. -26894.55 -29953. -23607.56 -1061.5 699.4757 1048.3 -690.1358 5171.8 -9672.159 -1448.6 5796.560 3928.8 -4269.361 -4732.8 16215.62 -20.844 0.000063 -5.2944 0.000064 36.585 0.0000MINIMUM VALUESELEM 39 4VALUE -51593. -51440.MAXIMUM VALUESELEM 40 51VALUE 23614. 51796.5.命令流文件/FILNAM,Structural/TITLE,Truss Bridge Static Analysis/COM,Structural/prep7et,1,beam4et,2,shell63sectype,1,beam,i,,0 !定义工字型截面secoffset,cent !截面至心不偏移secdata,0.4,0.4,0.4,0.016,0.016,0.016,0,0,0,0 !定义工字型截面参数sectype,2,beam,i,,0secoffset,centsecdata,0.4,0.4,0.4,0.012,0.012,0.012,0,0,0,0sectype,3,beam,i,,0secoffset,centsecdata,0.3,0.3,0.4,0.012,0.012,0.012,0,0,0,0r,1,0.0187,0.00017,0.00054,0.4,0.4,0, !定义单元实常数r,2,0.0141,0.128e-3,0.415e-3,0.4,0.4,,r,3,0.0117,0.541e-4,0.324e-3,0.3,0.4,,r,4,0.3,,,,,,MP,EX,1,2.1E11MP,PRXY,1,0.3MP,DENS,1,7850MP,EX,2,3.5E10MP,PRXY,2,0.1667MP,DENS,2,2500N,,0,0,-5,,,, !创建节点,复制结点NGEN,4,4,ALL,,,12,,,1,NGEN,2,1,ALL,,,,,10,1,NGEN,2,1,2,10,4,,16,,1,NGEN,2,1,3,11,4,,,-10,1,TYPE,1MAT,1REAL,1ESYS,0 !单元坐标系SECNUM,1TSHAP,LINEE,11,14 !建立单元TYPE,1 MAT,1 REAL,1 ESYS,0 SECNUM,2 TSHAP,LINE E,2,6E,6,10E,10,14 E,1,5E,5,9E,9,13E,3,7E,7,11E,4,8E,8,12E,1,2E,3,4E,5,6E,7,8E,9,10E,13,14 TYPE,1 MAT,1 REAL,1 ESYS,0 SECNUM,3 TSHAP,LINE E,3,6E,6,11E,4,5E,5,12E,2,3E,1,4E,6,7E,5,8E,10,11 E,9,12TYPE,2MAT,2REAL,1 ESYS,0TSHAP,QUADE,1,2,6,5E,5,6,10,9E,9,10,14,13NSYM,X,14,ALLESYM,,14,ALLNUMMRG,ALL,,,,LOW NUMCMP,ALL FINISH/SOLNSEL,S,,,23,24D,ALL,,,,,,UX,UY,UZ,,, NSEL,S,,,13,14D,ALL,,,,,,UY,UZ,,, NSEL,S,,,1,2F,ALL,FY,-100000 ALLSEL,ALL ACEL,0,10,0, ANTYPE,0SOLVEFINISH/POST1PLDISP,2PLNSOL,U,SUM,0,1PLVECT,U,,,,VECT,NODE,ON,0ETABLE,zhou_i,SMISC,1ETABLE,zhou_j,SMISC,7ETABLE,zhou_i,SMISC,2ETABLE,zhou_j,SMISC,8ETABLE,zhou_i,SMISC,6ETABLE,zhou_j,SMISC,12PRETAB,ZHOU_I,ZHOU_J,JIAN_I,JIAN_J,WAN_I,WAN_J PLLS,ZHOU_I,ZHOU_J,1,0PRNSOL,U,COMPFINISH/EXIT。
Ansys作业2桁架分析
(1) 进入ANSYS(设定工作目录和工作文件)1) 进入ANSYS菜单路径“程序>ANSYS >ANSYS10.0”2) 设置工作文件名菜单路径“file > Change Jobname”,弹出“Change Jobname”对话框,输入“CYKLink”,单击【OK】确定并关闭对话框。
(2) 设置计算类型菜单路径“ANSYS Main Menu: Preferences…”,在弹出的对话框中选择“Structural”,单击【OK】确定并关闭对话框。
(3) 选择单元类型菜单路径“ANSYS Main Menu: Preprocessor >Element Type>Add/Edit/Delete…”,在弹出“Library of Element Types”对话框中按照如图1所示参数选择,单击【OK】确定并关闭对话框。
图1 “Library of Element Types”对话框(4) 定义实常数菜单路径“ANSYS Main Menu: Preprocessor >Real Constants…>Add/Edit/Delete ”,在弹出的对话框中单击“Add > OK”,弹出如图2所示“Real Constant …”对话框,参数设置“AREA 0.000416”,单击【OK】确定并关闭对话框。
图2 “Real Constant …”对话框(5) 定义材料参数菜单路径“ANSYS Main Menu: Preprocessor > Material Props > Material Models”,在弹出的菜单中打开“Structural > Linear > Elastic > Isotropic”,弹出如图3所示“Linear Isotropic Material…”对话框,并设置如下参数。
图3 “Linear Isotropic Material…”对话框图4 “Beam Tool”对话框(6) 定义梁的截面菜单路径“ANSYS Main Menu: Preprocessor > Sections > Beam > Common Sections”,弹出如图4所示“Beam Tool”对话框,并按照图4设置,单击【OK】确定关闭对话框。
钢桁架桥梁结构的ANSYS分析
钢桁架桥梁结构的ANSYS分析摘要本文中采用有限元分析法,在大型有限元分析软件ANSYS平台上分析桥梁工程结构,很好地模拟桥梁的受力、应力情况等。
在静力分析中,通过加载各种载荷,得出结构变形图,找出桥梁的危险区域。
1、问题描述下面以一个简单桁架桥梁为例,以展示有限元分析的全过程。
该桁架桥由型钢组成,顶梁及侧梁,桥身弦杆,底梁分别采用3种不同型号的型钢,结构参数见表1-1。
桥长L=32m,桥高H=5.5m。
桥身由8段桁架组成,每段长4m。
该桥梁可以通行卡车,若这里仅考虑卡车位于桥梁中间位置,假设卡车的质量为4000kg,若取一半的模型,可以将卡车对桥梁的作用力简化为P1 ,P2和P3 ,其中P1= P3=5000 N, P2=10000N,见图1。
1图1桥梁的简化平面模型(取桥梁的一半)2、模型建立在桥梁结构模拟分析中,最常用的是梁单元和壳单元,鉴于桥梁的模型简化,采用普通梁单元beam3。
实体模型的建立过程为先生成关键点,再形成线,从而得到桁架桥梁的简化模型。
3、有限元模型3.1单元属性整个桥梁分成三部分,分别为顶梁及侧梁、弦杆梁、底梁,三者所使用的单元都为beam3单元,因其横截面积和惯性矩不同,所以设置3个实常数。
此外,他们材料都为型钢,材料属性视为相同,取为弹性模量EX为2.1e11 ,泊松比prxy为0.3,材料密度dens为7800。
3.2网格划分线单元尺寸大小为2,即每条线段的1/2。
4、计算4.1约束根据问题描述的要求,该桁架桥梁在x=0处的边界条件为全约束,x=32处的边界条件为y方向位移为0(即UY=0)。
如下图所示。
4.2载荷卡车对桥梁的压力视为3个集中载荷,因为模型只取桥梁的一般,所以3个集中载荷的力之和为20000N,分别为p1=5000N,p2=10000N,p3=5000N。
并将载荷施加在底梁的关键点4,5,6上。
如下图所示。
5、静力分析的计算结果5.1查看结构变形图显示y方向位移显示x方向位移5.2结论从加载后的结构变形图中可以看出,在载荷作用下,桁架桥的中间位置向下发生弯曲变形最为明显而两侧的侧梁变形最小,载荷引起的位移最大处在桥中间位置,随跨中间向两侧递减。
基于ANSYS对大型钢桁桥的模态分析
k, 6, ( b+ t w ) /2, tf k, 11, b, h a , 9 , 10 , 11, 12 am esh, a ll
k, 7, ( b + tw ) /2, h tf$
k , 8, ( b tw ) /2, h tf a , 5 , 6, 7 , 8$ sm rts ize , 8$
为检验 AN SY S模态分析结果, 使用 M idas / civ il桥 梁分析专用软件, 建立该桥的三维有限元模型 (图 7) , 进行模态分析, 将对比结果列于表 1。
表 1
阶次 AN SY S M idas
AN SY S 和 M idas 4阶 5阶 6阶 7阶 8阶
secwrite , d iaogan24, sect , , 1$ secoffset , cen t ,, , $ ase,l all$ aclear , all$
sectype , 24, beam, m esh,
secread, d iaogan24, sect , , m esh adele , al, l,, 1
所谓的模态分析就是确定结构系统 的振动特性, 得到结构的固有频率和 振型, 它们 是结构 承受 动力荷 载中的重要参数, 同时也 是结构其 它各类 型动 力学分 析的基础
[ 1]
( 4 ) 式称 为结构 振动的特 征方程, 模态 分析计 算 的是特征值
ANSYS命令流总结(全)
ANSYS命令流总结(全)ANSYS结构分析单元功能与特性/可以组成一一些命令,一般是一种总体命令(session),三十也有特殊,比如是处理/POST1! 是注释说明符号,,与其他软件的说明是一样的,ansys不作为命令读取,* 此符号一般是APDL的标识符,也就是ansys的参数化语言,如*do ,,,*enddo等等NSEL的意思是node select,即选择节点。
s就是select,选择。
DIM是定义数组的意思。
array 数组。
MP命令用来定义材料参数。
K是建立关键点命令。
K,关键点编号,x坐标,y坐标,z坐标。
K, NPT, X, Y, Z是定义关键点,K是命令,NPT是关键点编号,XYZ是坐标。
NUMMRG, keypoint 用这个命令,要保证关键点的位置完全一样,只是关键点号不一样的才行。
这个命令对于重复的线面都可以用。
这个很简单,压缩关键。
Ngen 复制节点e,节点号码:这个命令式通过节点来形成单元NUMCMP,ALL:压缩所有编号,这样你所有的线都会按次序重新编号~你要是需要固定的线固定的标号NSUBST,100,500,50:通过指定子步数来设置载荷步的子步LNSRCH线性搜索是求解非线性代数方程组的一种技巧,此法会在一段区间内,以一定的步长逐步搜索根,相比常用的牛顿迭代法所要耗费的计算量大得多,但它可以避免在一些情况下牛顿迭代法出现的跳跃现象。
LNSRCH激活线性搜索PRED 激活自由度求解预测NEQIT指定一个荷载步中的最大子步数AUTOTS 自动求解控制打开自动时间步长.KBC -指定阶段状或者用跳板装载里面一个负荷步骤。
SPLINE:P1,P2,P3,P4,P5,P6,XV1,YV1,ZV1,XV6,YV6,ZV6(生成分段样条曲线)*DIM,Par,Type,IMAX,JMAX,KMAX,Var1,Var2,Var3(定义载荷数组的名称)【注】Par: 数组名Type:array 数组,如同fortran,下标最小号为1,可以多达三维(缺省)char 字符串组(每个元素最多8个字符)tableIMAX,JMAX,KMAX各维的最大下标号Var1,Var2,Var3 各维变量名,缺省为row,column,plane(当type为table时)/config是设置ansys配置参数的命令格式为/CONFIG, Lab, V ALUELab为参数名称value为参数值例如:/config,MXEL,10000的意思是最大单元数为10000 杆单元:LINK1、8、10、11、180梁单元:BEAM3、4、23、24,44,54,188,189管单元:PIPE16,17,18,20,59,602D实体元: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杆单元单元名称简称节点数节点自由度特性备注LINK1 2D杆 2 Ux,Uy EPCSDGB常用杆元LINK8 3D杆Ux,Uy,Uz EPCSDGBLINK103D仅受拉或仅受压杆EDGB模拟缆索的松弛及间隙LINK11 3D线性调节器EGB模拟液压缸和大转动LINK180 3D有限应变杆EPCDFGB 另可考虑粘弹塑性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)等。
基于ANSYS对大型钢桁桥的模态分析
( c ol f i l n i n rnp r t ni S uhC iaU i.o eh G a gh u5 0 4 , hn ) S ho o Cv g.a dTa sot i o t hn nv f c., u nzo 16 1 C ia iE ao n T
息, 而且可以用于指导结构 的设计 。工 程 中的共 振往 往是在桥梁工程中所不希 望发生 的 , 也同样依 赖 于 这 对结构 的固有频率和振型的准确获取 。
( ) 建模 。进行 前处 理 , 义 单元 类 型 , 常 1 定 实 数, 材料特性 等 , 因为模 态分析 为线性分 析 , 因此选用
Ab t a t Mo a n lss i h  ̄u d to f d n mi nay i a d i tk s a mpot n a t i sr c : d la ay i s t e n a in o y a c a l ss, n t a e n i ra tp r n sr cr ldy a c a ay i.Att e b gn i g o hi a e , ic i e h a i he r fmo a n l — tu ta n mi n l ss h e i n n ft s p p r we d s rb d t e b sc t o y o d la ay ss, sa ls i g a d n mi d lo a g r s te r d e wi i e tb ih n y a c mo e fa lr e tu s se lb a g t ANS h YS, n h n we g tt e n tu l a d t e o h a r a fe u n y a d mo e s p . a t we u e d s ii t n l z h r d e, n o a e h e u t r q e c n d ha e Atls , s d mi a /cv l o a ay e t e b a g a d c mp r d t e r s ls wi h t ANS YS, n ts o d t a d la ay i a e n ANS s c e be a d i h we h tmo a n lss b s d o YS i r di l.
ANSYS命令流总结(全)
ANSYS结构分析单元功能与特性/可以组成一一些命令,一般是一种总体命令(session),三十也有特殊,比如是处理/POST1! 是注释说明符号,,与其他软件的说明是一样的,ansys不作为命令读取,* 此符号一般是APDL的标识符,也就是ansys的参数化语言,如*do ,,,*enddo等等NSEL的意思是node select,即选择节点。
s就是select,选择。
DIM 是定义数组的意思。
array 数组。
MP命令用来定义材料参数。
K是建立关键点命令。
K,关键点编号,x坐标,y坐标,z坐标。
K, NPT, X, Y, Z是定义关键点,K是命令,NPT是关键点编号,XYZ是坐标。
NUMMRG, keypoint 用这个命令,要保证关键点的位置完全一样,只是关键点号不一样的才行。
这个命令对于重复的线面都可以用。
这个很简单,压缩关键。
Ngen 复制节点e,节点号码:这个命令式通过节点来形成单元NUMCMP,ALL:压缩所有编号,这样你所有的线都会按次序重新编号~你要是需要固定的线固定的标号NSUBST,100,500,50 :通过指定子步数来设置载荷步的子步LNSRCH线性搜索是求解非线性代数方程组的一种技巧,此法会在一段区间内,以一定的步长逐步搜索根,相比常用的牛顿迭代法所要耗费的计算量大得多,但它可以避免在一些情况下牛顿迭代法出现的跳跃现象。
LNSRCH 激活线性搜索PRED 激活自由度求解预测NEQIT指定一个荷载步中的最大子步数AUTOTS 自动求解控制打开自动时间步长.KBC -指定阶段状或者用跳板装载里面一个负荷步骤。
SPLINE:P1,P2,P3,P4,P5,P6,XV1,YV1,ZV1,XV6,YV6,ZV6(生成分段样条曲线)*DIM,Par,Type,IMAX,JMAX,KMAX,Var1,Var2,Var3(定义载荷数组的名称)【注】Par: 数组名Type: array 数组,如同fortran,下标最小号为1,可以多达三维(缺省)char 字符串组(每个元素最多8个字符)tableIMAX,JMAX,KMAX 各维的最大下标号Var1,Var2,Var3 各维变量名,缺省为row,column,plane(当type 为table时)/config是设置ansys配置参数的命令格式为/CONFIG, Lab, VALUELab为参数名称 value为参数值例如:/config,MXEL,10000的意思是最大单元数为10000杆单元:LINK1、8、10、11、180梁单元:BEAM3、4、23、24,44,54,188,189管单元:PIPE16,17,18,20,59,602D实体元:PLANE2,25,42,82,83,145,146,182,1833D实体元: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杆单元单元名称简称节点数节点自由度特性备注LINK1 2D杆 2 Ux,Uy EPCSDGB常用杆元LINK8 3D杆Ux,Uy,Uz EPCSDGBLINK103D仅受拉或仅受压杆EDGB模拟缆索的松弛及间隙LINK11 3D线性调节器EGB模拟液压缸和大转动LINK183D有限应EPCDFG另可考虑粘0 变杆 B 弹塑性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)等。
全套完整版ANSYS命令流教学手册
全套完整版ANSYS命令流教学手册————————————————————————————————作者:————————————————————————————————日期:目录1ANSYS结构分析单元功能与特性 01.1杆单元: LINK1、8、10、11、180 01.2梁单元:BEAM3、4、23、24,44,54,188,189 (1)1.3管单元:PIPE16,17,18,20,59,60 (1)1.42D实体单元:PLANE2,25,42,82,83,145,146,182,183 21.53D实体元:SOLID45,46,64,65,72,73,92,95,147,148,185,186,187,191 (4)1.6壳单元:SHELL28,41,43,51,61,63,91,93,99,143,150,181,208,209 (5)1.7弹簧单元:COMBIN7,14,37,39,40 (5)1.8质量单元:MASS21 (5)1.9接触单元:CONTAC12,52,TARGE169,170,CONTA171,172,173,174,175,178 (5)1.10 矩阵单元:MATRIX27,50 (6)1.11 表面效应元:SURF153,154 (6)1.12 粘弹实体元:VISCO88,89,106,107,108, (7)1.13 超弹实体元:HYPER56,58,74,84,86,158 (7)1.14 耦合场单元:SOLID5,PLANE13,FLUID29,30,38,SOLID62,FLUID79,FLUID80,81,SOLID98,FLUID129,INFIN110,111,FLUID116,130 (7)1.15 界面单元:INTER192,193,194,195 (7)1.16 显式动力分析单元:LINK160,BEAM161,PLANE162,SHELL163,SOLID164,COMBI16 (7)1.17 预紧、多点约束、网分单元 (7)2ANSYS 的基本使用 (9)2.1ANSYS环境简介 (9)2.2有限元法的基本构架 (11)2.3ANSYS架构及命令 (11)2.4典型的分析过程 (14)2.5ANSYS 文件及工作文件名 (15)2.6图形控制 (16)3第三章有限元模型的建立 (18)3.1建模方法 (18)3.2坐标系统及工作平面 (19)3.2.1局部坐标系 (19)3.2.1.1激活总体和局部坐标系(声明坐标系统)193.2.1.2根据总体坐标系定义局部坐标系193.2.1.3根据已有的三个节点定义局部坐标系203.2.1.4根据已有的三个关键点定义局部坐标系203.2.1.5根据当前工作平面定义局部坐标系203.2.1.6根据激活的坐标系定义局部坐标系203.2.1.7删除局部坐标系213.2.1.8查看激活坐标系和局部坐标系213.2.2节点坐标系的旋转与修改 (21)3.2.2.1将某些节点的坐标系旋转到与当前激活坐标系(简称“当前坐标系”)方向一致 (21)3.2.2.2将既有节点的节点坐标系旋转某个角度213.2.2.3在创建节点时直接定义其坐标系的旋转角度213.2.2.4按方向余弦旋转节点坐标系223.2.2.5节点坐标系列表223.2.3单元坐标系的定义与修改 (22)3.2.3.1设置单元坐标系223.2.3.2修改单元坐标系方向223.2.3.3激活显示坐标系223.2.3.4激活结果坐标系233.2.4工作平面 (23)3.2.4.1定义工作平面233.2.4.2通过3个坐标点定义工作平面233.2.4.3通过3个节点定义工作平面233.2.4.4通过3个关键点定义工作平面243.2.4.5通过垂直于线上的某个位置定义工作平面243.2.4.6查看工作平面的当前状态243.2.4.7恢复到ANSYS默认状态243.2.4.8移动工作平面243.2.5声明单位系统 (25)3.3节点定义 (26)3.3.1创建节点 (26)3.3.2删除节点 (26)3.3.3节点显示 (26)3.3.4节点列式 (27)3.3.5节点的填充 (27)3.3.6节点复制 (27)3.4单元的定义 (28)3.4.1定义单元类型 (28)3.4.2定义材料的属性 (29)3.4.3定义实常数 (29)3.4.4定义单元的连接方式 (29)3.4.5单元复制 (30)3.4.6单元显示 (30)3.4.7单元列示 (30)3.4.8声明使用已定义单元类型 (30)3.4.9声明使用已定义的实常数 (30)3.4.10声明使用已定义的单元属性, (31)3.5负载定义 (31)3.5.1进入解题处理器 (32)3.5.2声明分析类型 (32)3.5.3定义节点的集中力 (32)3.5.4定义节点自由度 (33)3.5.5定义在梁单元上的分布力 (33)3.5.6定义分布力作用于单元上的方式和大小 (34)3.5.7定义节点间分布力 (35)3.5.8对象选择 (35)3.6求解 (37)3.7用POST1进行结果后处理 (37)3.7.1进入POST1 (37)3.7.2读取结果 (38)3.7.3绘变形图 (38)3.7.4变形动画 (38)3.7.5列表支反力3.7.6应力等值线与应力等值线动画 (38)3.7.7应力等值线动画 (38)4实体模型的建立 (42)4.1实体模型简介 (42)4.2实体模型的建立方法 (42)4.3群组命令介绍 (43)4.4点定义 (44)4.4.1创建关键点 (44)4.4.1.1在给定坐标点创建关键点444.4.1.2在两关键点之间创建一个关键点444.4.1.3在两关键点之间填充多个点454.4.1.4复制创建关键点454.4.1.5镜像创建关键点464.4.1.6列表显示关键点信息464.4.1.7屏幕上显示关键点464.4.1.8删除关键点474.4.1.9选择关键点474.4.1.10选择与所选线相关的关键点484.4.1.11修改关键点坐标484.4.1.12定义点(NPT)于已知节点上494.4.2创建线4.4.2.1通过两关键点创建线494.4.2.2通过两关键点创建直线504.4.2.3通过关键点创建圆弧线504.4.2.4创建圆或圆弧线514.4.2.5两条相交线倒角创建圆弧线524.4.2.6复制创建线534.4.2.7合并两条或多条线534.4.2.8将一条线分为多条线544.4.2.9延长一条线544.4.2.10通过多个关键点按样条创建一条曲线554.4.2.11关键点绕轴线创建旋转线554.4.2.12通过坐标轴镜像创建线564.4.2.13显示线和删除线564.4.2.14列表输出线信息564.4.2.15选择一组线574.4.2.16选择与面相关的线574.4.2.17选择与关键点相关的线584.4.2.18练习点和线段的生成584.4.3创建面 (60)4.4.3.1通过关键点创建面604.4.3.2通过线创建面614.4.3.3沿路径拖拉创建面614.4.3.4线绕轴旋转生成弧面624.4.3.5既有面偏移创建新面624.4.3.6蒙皮创建光滑曲面634.4.3.7复制创建面634.4.3.8通过坐标轴对称创建面644.4.3.9列表输出面信息644.4.3.10显示面644.4.3.11删除面644.4.3.12选择一组面654.4.3.13选择与所选线相关的面664.4.3.14选择与所选体相关的面664.4.3.15通过两角点坐标创建矩形面664.4.3.16通过一角点坐标和尺寸创建矩形面664.4.3.17通过中心坐标和尺寸创建矩形面674.4.3.18在工作平面原点创建圆面或环面674.4.3.19通过圆心坐标和半径等创建圆或环面674.4.3.20通过圆上直径端点坐标创建圆面684.4.3.21在工作平面原点创建正多边形面684.4.3.22在工作平面任意位置创建正多边形面684.4.4创建体 (68)4.4.4.1通过关键点创建体4.4.4.2通过面创建体714.4.4.3沿路径拖拉面创建体724.4.4.4面绕轴旋转创建柱体724.4.4.5面偏移创建体734.4.4.6通过面延伸创建体744.4.4.7复制创建体744.4.4.8通过坐标轴镜像创建体754.4.4.9列表输出体信息754.4.4.10显示体754.4.4.11删除体754.4.4.12创建长方体754.4.4.13通过一角点坐标和尺寸创建长方体764.4.4.14通过面中心坐标和尺寸创建长方体764.4.4.15在工作平面原点创建圆柱体或部分圆柱体764.4.4.16通过圆心坐标和半径等创建圆柱体764.4.4.17通过圆上直径两端点坐标创建圆柱体764.4.4.18在工作平面原点创建正棱柱体764.4.4.19在工作平面任意位置创建正棱柱体774.4.4.20在工作平面原点创建球体774.4.4.21在工作平面任意位置创建球体774.4.4.22通过直径端点生成球体774.4.4.23以工作平面原点为圆心创建圆锥体4.4.4.24在工作平面任意位置创建圆锥体784.4.4.25以工作平面原点为环心创建环体784.5布尔操作 (79)4.5.1布尔运算的一般设置 (80)4.5.1.1布尔运算的容差设置804.5.1.2交运算Intersection814.5.1.3加运算Addition824.5.1.4减运算Subtract834.5.1.5用工作平面切分图素Subtract844.5.1.6分割运算Partition854.5.1.7分类运算Classify864.5.1.8搭接运算Overlap864.5.1.9粘接Glue (或Merge)874.5.2几何建模的其它常用命令 (88)4.5.2.1图形平移、缩放和旋转884.5.2.2设置坐标轴方向884.5.2.3设置视图方向894.5.2.4设置视图旋转角度894.5.2.5编号显示控制904.5.2.6颜色显示控制914.5.2.7显示边界条件和荷载的符号及数值914.5.2.8显示边界条件及数值4.5.2.9显示风格设置924.5.2.10单元尺寸和形状924.5.2.11图素收缩显示控制924.5.2.12显示单元形状934.5.2.13等值线显示控制934.5.2.14均匀等值线设置934.5.2.15设置等值线的文字标注944.5.2.16颜色设置944.5.2.17设置图形中浮点数显示方式954.5.2.18设置变形放大系数964.5.2.19设置矢量显式长度964.5.2.20设置窗口布局964.5.2.21图素显示控制974.5.2.22显示所有图素974.5.2.23图形擦除985网格划分 (98)5.1区分实体模型和有限元模型 (98)5.2网格化的步骤 (99)5.2.1定义单元类型 (101)5.2.1.1单元类型的KEYOPT1015.2.1.2自由度集1015.2.1.3改变单元类型5.2.1.4单元类型的删除与列表1015.2.2实常数 (102)5.2.2.1定义实常数1025.2.2.2变厚度壳实常数定义1025.2.2.3实常数组的删除与列表1025.2.3材料属性 (102)5.2.3.1定义线性材料属性1025.2.3.2定义线性材料属性的温度表1035.2.3.3定义与温度对应的线性材料特性1035.2.3.4复制线性材料属性组1035.2.3.5改变指定单元的材料参考号1035.2.3.6线性材料属性列表和删除1035.2.3.7修改与线胀系数相关的温度1045.2.3.8计算生成线性材料温度表1045.2.3.9绘制线性材料特性曲线1045.2.3.10设置材料库读写的缺省路径1045.2.3.11读入材料库文件1045.2.3.12将材料属性写入文件1045.2.3.13激活非线性材料属性的数据表1045.2.3.14定义TB温度值1045.2.3.15定义TB数据表中的数据5.2.3.16定义非线性数据曲线上的一个点1045.2.3.17非线性材料数据表的删除和列表1055.2.3.18非线性材料数据表的绘图1055.2.4定义截面类型和截面ID (105)5.2.4.1定义梁截面几何数据(Type=BEAM)1065.2.4.2定义变截面梁几何数据(Type=TAPER)1085.2.4.3定义截面偏移1085.2.4.4梁截面特性列表1095.2.4.5删除所定义的截面1095.2.4.6绘制所定义截面1095.2.4.7自定义截面的存盘和读入1095.2.4.8定义层壳单元的数据(Type=SHELL)1105.2.4.9定义预紧截面的数据(Type= PRETENSION)1105.2.4.10修改预紧截面数据1105.2.4.11定义连接数据(Type=JOINT)1105.2.5设置几何模型的单元属性 (110)5.2.5.1设置关键点单元属性1105.2.5.2设置线的单元属性1115.2.5.3设置面的单元属性1115.2.5.4设置体的单元属性1115.2.5.5单元形状控制5.2.5.6网格类型选择1125.2.5.7中间节点的位置控制1135.2.6单元尺寸控制 (113)5.2.6.1映射网格单元尺寸控制的DESIZE命令1135.2.6.2自由网格单元尺寸控制的SMRTSIZE命令1145.2.6.3线的单元尺寸定义1145.2.6.4关键点最近处单元边长定义1155.2.6.5面内部的单元尺寸定义1165.3网格划分工具 (116)5.3.1在关键点处生成点单元 (118)5.3.2在几何线上生成线单元 (118)5.3.3在几何面上生成面单元 (118)5.3.4在几何体上生成体单元 (118)5.3.4.1单元有效性检查--单元形状参数限值设置1195.3.5网格修改 (121)5.3.5.1关键点网格清除1215.3.5.2线网格清除1215.3.5.3面网格清除1215.3.5.4体网格清除5.3.6细化局部网格 (121)5.3.6.1节点附近细化1215.3.6.2单元附近细化1215.3.6.3关键点附近细化1215.3.6.4线附近细化1225.3.6.5面附近细化1226第六章实体模型的外力 (122)6.1对节点施加自由度约束 (123)6.2关键点自由度约束及相关命令 (123)6.2.1列表和删除关键点自由度约束的命令 (124)6.3对线施加自由度约束 (124)6.3.1而列表和删除线上自由度约束 (125)6.4对面施加自由度约束 (125)6.4.1列表和删除面上自由度约束 (125)6.5约束转换命令 (125)6.5.1仅转换约束自由度 (125)6.5.2边界条件和荷载转换 (125)6.6施加集中荷载 (126)6.6.1施加节点集中荷载 (126)6.6.1.1节点集中荷载列表:FLIST1266.6.1.2删除节点集中荷载:FDELE1266.6.2施加关键点集中荷载 (126)6.6.3施加面荷载 (126)6.6.3.1对节点群施加面荷载1266.6.3.2定义节点号与面荷载的函数关系1276.6.3.3定义面荷载梯度1276.6.4在单元上施加面荷载 (127)6.6.5在梁单元施加面荷载 (128)6.6.6在线上施加面荷载 (128)6.6.7在面上施加面荷载 (129)6.6.8施加体荷载 (129)6.6.9施加惯性荷载 (130)6.6.10施加耦合场荷载 (130)6.6.11初应力荷载及施加 (130)6.6.11.1施加初始常应力荷载1306.6.11.2从文件施加初应力荷载1307输出选项 (133)7.1控制写入数据库和结果文件的结果数据 (133)7.2结果输出控制 (134)7.3图形求解追踪器 (135)8分析类型与求解控制选项 (135)9通用与时间历程后处理技术 (135)9.1定义矢量和主轴的计算方法 (135)9.2定义结果数据平均处理 (136)9.3控制壳或层壳单元数据的位置 (136)9.4定义数据输出格式 (137)9.5每页的标题输出控制 (137)9.6显示结构变形图 (138)9.7显示节点结果 (138)9.8显示单元结果 (138)9.9以矢量方式显示结果图 (139)9.10显示裂缝或压碎图 (139)9.11列出节点结果 (140)9.12列出单元结果 (140)9.13生成单元表 (140)9.14云图显示单元表结果 (140)9.15列表显示单元表结果 (141)9.16单元表绝对值操作 (141)9.17计算并输出单元表数据之和 (141)9.18映射结果到路径上 (142)9.19图形显示路径项数据 (142)9.20沿路径几何形状显示路径项数据 (142)9.21列表显式路径项数据 (143)9.22对路径项数据运算 (143)9.22.1加运算 (143)9.22.2乘运算 (143)9.22.3除运算 (143)9.22.4幂运算 (143)9.22.5求导运算 (144)9.22.6积分运算 (144)10文件管理 (144)1ANSYS结构分析单元功能与特性1.1杆单元: LINK1、8、10、11、180单元名称简称节点数节点自由度特性备注LINK1 2D杆2 Ux,Uy EPCSDGB 常用杆元LINK8 3D杆Ux,Uy,Uz EPCSDGBLINK103D仅受拉或仅受压杆EDGB模拟缆索的松弛及间隙LINK11 3D线性调节器EGB 模拟液压缸和大转动LINK180 3D有限应变杆EPCDFGB 另可考虑粘弹塑性注: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)等。
Ansys命令流大全(整理)
Ansys命令流大全(整理)1、A,P1,P2,P3,P4,P5,P6,P7,P8,P9此命令用已知的一组关键点点(P1~P9)来定义面(Area),最少使用三个点才能围成面,同时产生转围绕些面的线。
点要依次序输入,输入的顺序会决定面的法线方向。
如果超过四个点,则这些点必须在同一个平面上。
Menu Paths:Main Menu>Preprocessor>Create>Arbitrary>Through KPs2、*ABBR,Abbr,String--定义一个缩略语.Abbr:用来表示字符串"String"的缩略语,长度不超过8个字符.String:将由"Abbr"表示的字符串,长度不超过60个字符.3、ABBRES,Lab,Fname,Ext-从一个编码文件中读出缩略语.Lab:指定读操作的标题,NEW:用这些读出的缩略语重新取代当前的缩略语(默认)CHANGE:将读出的缩略语添加到当前缩略语阵列,并替代现存同名的缩略语.Ext:如果"Fname"是空的,则缺省的扩展命是"ABBR".4、ABBSA V,Lab,Fname,Ext-将当前的缩略语写入一个文本文件里Lab:指定写操作的标题,若为ALL,表示将所有的缩略语都写入文件(默认)5、add, ir, ia,ib,ic,name,--,--,facta, factb, factc将ia,ib,ic变量相加赋给ir变量ir, ia,ib,ic:变量号name: 变量的名称6、Adele,na1,na2,ninc,kswp !kswp=0时只删除掉面积本身,=1时低单元点一并删除。
7、Adrag, nl1,nl2,nl3,nl4,nl5,nl6, nlp1,nlp2,nlp3,nlp4,nlp5,nlp6 !面积的建立,沿某组线段路径,拉伸而成。
8、Afillt,na1,na2,rad !建立圆角面积,在两相交平面间产生曲面,rad为半径。
钢管溷凝土拱桥ansys命令流——很实用精品文档27页
fini/clear/title,concrete-filled tube arch bridge /prep7!!!截面1et,1,82cyl4,0.6,0.6,0.6,,0.586cyl4,0.6,2.4,0.6,,0.586rectng,0.35,0.364,0.8,2.2rectng,0.85,0.836,0.8,2.2allselaptn,alladele,5,8,1allselaadd,allsmrtsize,5amesh,allsecwrite,gg,sect,,1sectype,1,beam,meshsecoffset,cent,,,secread,'gg','sect',,meshasel,allaclear,alladele,all,,,1!!!截面2cyl4,0.6,0.6,0.586cyl4,0.6,2.4,0.586rectng,0.364,0.836,0.8,2.2allselaadd,allsmrtsize,5amesh,allsecwrite,hnt,sect,,1sectype,2,beam,meshsecoffset,cent,,,secread,'hnt','sect',,meshasel,allaclear,alladele,all,,,1!!!!!!建立横梁标准截面a-ak,1,-0.5k,2,-0.75,0.25k,3,-0.75,0.5k,4,-0.55,0.5k,5,-0.55,0.85k,6,0.55,0.85k,7,0.55,0.5k,8,0.75,0.5k,9,0.75,0.25k,10,0.5,0a,1,2,3,4,5,6,7,8,9,10 smrtsize,5amesh,allsecwrite,hl1,sect,,1 sectype,4,beam,mesh secoffset,cent,,, secread,'hl1','sect',,mesh asel,allaclear,alladele,all,,,1!!!!!!建立横梁标准截面b-b k,1,-0.5k,2,-0.5,1k,3,-0.75,1.25k,4,-0.75,1.5k,5,0.75,1.5k,6,0.75,1.25k,7,0.5,1k,8,0.5a,1,2,3,4,5,6,7,8 smrtsize,5amesh,allsecwrite,hl2,sect,,1 sectype,5,beam,mesh secoffset,cent,,, secread,'hl2','sect',,mesh asel,allaclear,alladele,all,,,1!!!!!!建立横梁标准截面c-c k,1,-0.5k,2,-0.5,1.04k,3,-0.75,1.04+0.25k,4,-0.75,1.04+0.5k,5,-0.55,1.04+0.5k,6,-0.55,1.04+0.5+0.35k,7,0.55,1.04+0.5+0.35k,8,0.55,1.04+0.5k,9,0.75,1.04+0.5k,10,0.75,1.04+0.25k,11,0.5,1.04k,12,0.5k,13,-0.2,0.58k,14,-0.2,1.04+0.5-0.3k,15,0.2,1.04+0.5-0.3k,16,0.2,0.58a,1,2,3,4,5,6,7,8,9,10,11,12 a,13,14,15,16asba,1,2smrtsize,5amesh,allsecwrite,hl3,sect,,1 sectype,6,beam,mesh secoffset,cent,,, secread,'hl3','sect',,mesh asel,allaclear,alladele,all,,,1!!!!!!建立横梁标准截面d-d ldele,allkdele,allk,1,-0.5k,2,-0.5,1.23k,3,-0.75,1.23+0.25k,4,-0.75,1.23+0.5k,5,-0.55,1.23+0.5k,6,-0.55,1.23+0.5+0.35k,7,0.55,1.23+0.5+0.35k,8,0.55,1.23+0.5k,9,0.75,1.23+0.5k,10,0.75,1.23+0.25k,11,0.5,1.23k,12,0.5k,13,-0.2,0.38k,14,-0.3,0.38+0.2k,15,-0.3,1.23+0.5-0.1-0.2 k,16,-0.2,1.23+0.5-0.1k,17,0.2,1.23+0.5-0.1k,18,0.3,1.23+0.5-0.1-0.2 k,19,0.3,0.38+0.2k,20,0.2,0.38a,1,2,3,4,5,6,7,8,9,10,11,12a,13,14,15,16,17,18,19,20asba,1,2smrtsize,5amesh,allsecwrite,hl4,sect,,1sectype,7,beam,meshsecoffset,cent,,,secread,'hl4','sect',,meshasel,allaclear,alladele,all,,,1!!!!!!建立横梁标准截面e-ek,1,-0.5k,2,-0.5,1.23k,3,-0.75,1.23+0.25k,4,-0.75,1.23+0.5k,5,-0.55,1.23+0.5k,6,-0.55,1.23+0.5+0.35k,7,0.55,1.23+0.5+0.35k,8,0.55,1.23+0.5k,9,0.75,1.23+0.5k,10,0.75,1.23+0.25k,11,0.5,1.23k,12,0.5a,1,2,3,4,5,6,7,8,9,10,11,12smrtsize,5amesh,allsecwrite,hl5,sect,,1sectype,8,beam,meshsecoffset,cent,,,secread,'hl5','sect',,meshasel,allaclear,alladele,all,,,1etdele,1 !!将辅助单元删除!定义所有材料特性et,1,beam44 !!钢管特性mp,ex,1,2.1e11mp,dens,1,7800mp,prxy,1,0.3n,90000,0,0,30 !!参考点et,2,beam44 !!钢管内50#混凝土特性mp,ex,2,3.5e10mp,dens,2,2600mp,prxy,2,0.1667et,3,beam44 !!纵梁30#混凝土钢管特性mp,ex,3,3.0e10mp,dens,3,2600mp,prxy,3,0.1667et,4,beam44 !!横梁30#混凝土钢管特性mp,ex,4,3.0e10mp,dens,4,2600mp,prxy,4,0.1667et,5,beam44 !!风撑特性mp,ex,5,2.1e11mp,dens,5,7800mp,prxy,5,0.3et,6,link10 !!吊杆特性(钢绞线)mp,ex,6,1.9e11mp,dens,6,7800mp,prxy,6,0.3keyopt,6,3,0 !只拉吊杆et,7,beam44 !!盖梁30#混凝土特性mp,ex,7,3.0e10mp,dens,7,2600mp,prxy,7,0.1667et,8,beam44 !!墩柱40#混凝土特性mp,ex,8,3.3e10mp,dens,8,2600mp,prxy,8,0.1667et,9,beam44 !!横梁30#混凝土特性mp,ex,9,3.0e10mp,dens,9,2600mp,prxy,9,0.1667et,10,link10 !!系杆特性(钢绞线)mp,ex,10,2.0e11mp,dens,10,7800mp,prxy,10,0.3keyopt,10,3,0 !只拉系杆et,11,beam44 !!承台桩基基础25#混凝土特性mp,ex,11,2.85e10mp,dens,11,2500mp,prxy,11,0.1667!开始建立模型!0号横梁(端横梁)n,1,,2.025n,2,6.25,2.025n,3,6.25+5.1,2.025n,4,6.25+5.1+0.2+1.9/2,2.025 !12.5m位置n,5,6.25+5.1+2.3,2.025n,6,6.25+5.1+2.3+2.23,2.025nsym,x,10,2,6,1type,4r,1real,1mat,4secnum,7!单元1到10号为0号横梁e,1,2secnum,6e,2,3secnum,5e,3,4e,4,5secnum,4e,5,6secnum,7e,1,12secnum,6e,12,13secnum,5e,13,14e,14,15secnum,4e,15,16!!1号端横梁egen,2,100,1,10,1,,,,,,0,4.7+1.25/2+1.5/2!!其余端横梁2~17号egen,17,100,11,20,1,,,,,,0,5.3+1.5!!18号端横梁egen,2,100,171,180,1,,,,,,0,4.7+1.25/2+1.5/2 cm,hl,elemesel,none!纵梁单元type,3r,2real,2mat,3sectype,10,beam,rect,zlsecdata,0.4,1.2secnum,10!0号纵梁e,1,101e,2,102e,3,103e,5,105e,12,112e,13,113e,15,115!中间纵梁1~16号*do,i,0,15,1e,101+i*100,201+i*100e,102+i*100,202+i*100e,104+i*100,204+i*100e,112+i*100,212+i*100e,114+i*100,214+i*100*enddo!17号纵梁(端横梁)e,1701,1801e,1702,1802e,1703,1803e,1705,1805e,1712,1812e,1713,1813e,1715,1815cm,zl,elemesel,none!定义拱轴线*dim,y,array,16 !将纵桥向距离定义为数组*dim,cs,array,16 !定义分值cs=(1-2x/l)*dim,z,array,16 !将拱高度定义为数组!吊杆上的点!首先定义吊杆的y向(顺桥向)坐标y(1)=8.1y(2)=14.9y(3)=21.7y(4)=28.5y(5)=35.3y(6)=42.1y(7)=48.9y(8)=55.7y(9)=62.5!风撑上的点!首先定义风撑的y向(顺桥向)坐标y(11)=21.7y(12)=28.5y(13)=42.1y(14)=48.9y(15)=59.1!0号横梁的y向坐标y(16)=2.025fl=25 !矢高ml=1.1k=log(ml+sqrt(ml*ml-1)) !k=ch-1(m)*do,i,1,16,1cs(i)=1-y(i)*2/125z(i)=25-fl/(ml-1)*(cosh(k*cs(i))-1) !拱轴高度*enddo!上述定义的拱轴线可参照姚玲森《桥梁工程》p1299*dim,ygd,array,23 !定义预工拱度(是在后续工况中反复调算的结果)!初始运行中可将预拱度全部设为0ygd(1)=0.4000e-02+1.29e-02ygd(2)=0.10685e-01+2.22e-02ygd(3)=0.19483e-01+3.03e-02ygd(4)=0.29219e-01+3.69e-02ygd(5)=0.39022e-01+4.23e-02ygd(6)=0.48243e-01+4.65e-02ygd(7)=0.56390e-01+4.95e-02ygd(8)=0.63070e-01+5.15e-02ygd(9)=0.68014e-01+5.21e-02ygd(10)=0.71003e-01+5.15e-02ygd(11)=0.71739e-01+4.95e-02ygd(12)=0.69962e-01+4.65e-02ygd(13)=0.65514e-01+4.23e-02ygd(14)=0.58289e-01+3.69e-02ygd(15)=0.48248e-01+3.04e-02ygd(16)=0.35507e-01+2.24e-02ygd(17)=0.20415e-01+1.23e-02ygd(18)=0.65770e-01+0.40234e-01ygd(19)=0.69769e-01+0.41215e-01!生成各吊杆的竖向节点(z向)ngen,2,10000,104,114,10,0,0,z(1)+ygd(1)ngen,2,10000,204,214,10,0,0,z(2)+ygd(2)ngen,2,10000,304,314,10,0,0,z(3)+ygd(3)ngen,2,10000,404,414,10,0,0,z(4)+ygd(4)ngen,2,10000,504,514,10,0,0,z(5)+ygd(5)ngen,2,10000,604,614,10,0,0,z(6)+ygd(6)ngen,2,10000,704,714,10,0,0,z(7)+ygd(7)ngen,2,10000,804,814,10,0,0,z(8)+ygd(8)ngen,2,10000,904,914,10,0,0,z(9)+ygd(9)ngen,2,10000,1004,1014,10,0,0,z(8)+ygd(10)ngen,2,10000,1104,1114,10,0,0,z(7)+ygd(11)ngen,2,10000,1204,1214,10,0,0,z(6)+ygd(12)ngen,2,10000,1304,1314,10,0,0,z(5)+ygd(13)ngen,2,10000,1404,1414,10,0,0,z(4)+ygd(14)ngen,2,10000,1504,1514,10,0,0,z(3)+ygd(15)ngen,2,10000,1604,1614,10,0,0,z(2)+ygd(16)ngen,2,10000,1704,1714,10,0,0,z(1)+ygd(17)ngen,2,3000,104,114,10,0,0,0ngen,2,3000,204,214,10,0,0,0ngen,2,3000,304,314,10,0,0,0ngen,2,3000,404,414,10,0,0,0ngen,2,3000,504,514,10,0,0,0ngen,2,3000,604,614,10,0,0,0ngen,2,3000,704,714,10,0,0,0ngen,2,3000,804,814,10,0,0,0ngen,2,3000,904,914,10,0,0,0ngen,2,3000,1004,1014,10,0,0,0ngen,2,3000,1104,1114,10,0,0,0ngen,2,3000,1204,1214,10,0,0,0ngen,2,3000,1304,1314,10,0,0,0ngen,2,3000,1404,1414,10,0,0,0ngen,2,3000,1504,1514,10,0,0,0ngen,2,3000,1604,1614,10,0,0,0ngen,2,3000,1704,1714,10,0,0,0ngen,2,6000,104,114,10,0,0,z(1)+ygd(1)ngen,2,6000,204,214,10,0,0,z(2)+ygd(2)ngen,2,6000,304,314,10,0,0,z(3)+ygd(3)ngen,2,6000,404,414,10,0,0,z(4)+ygd(4)ngen,2,6000,504,514,10,0,0,z(5)+ygd(5)ngen,2,6000,604,614,10,0,0,z(6)+ygd(6)ngen,2,6000,704,714,10,0,0,z(7)+ygd(7)ngen,2,6000,804,814,10,0,0,z(8)+ygd(8)ngen,2,6000,904,914,10,0,0,z(9)+ygd(9)ngen,2,6000,1004,1014,10,0,0,z(8)+ygd(10)ngen,2,6000,1104,1114,10,0,0,z(7)+ygd(11)ngen,2,6000,1204,1214,10,0,0,z(6)+ygd(12)ngen,2,6000,1304,1314,10,0,0,z(5)+ygd(13)ngen,2,6000,1404,1414,10,0,0,z(4)+ygd(14)ngen,2,6000,1504,1514,10,0,0,z(3)+ygd(15)ngen,2,6000,1604,1614,10,0,0,z(2)+ygd(16)ngen,2,6000,1704,1714,10,0,0,z(1)+ygd(17)!*do,i,1,9,1!ngen,2,10000,i*100+4,i*100+14,10,0,0,z(i)+ygd(i) !吊杆节点!*enddo!jj=1!*do,i,10,17,1!ngen,2,10000,i*100+4,i*100+14,10,0,0,z(i-2*jj)+ygd(i) !吊杆节点!jj=jj+1!*enddo!!为了对以后吊杆转动进行放松,重复吊杆节点!*do,i,1,9,1!ngen,2,3000,i*100+4,i*100+14,10,0,0,z(i)+ygd(i) !吊杆下节点!*enddo!jj=1!*do,i,10,17,1!ngen,2,3000,i*100+4,i*100+14,10,0,0,z(i-2*jj)+ygd(i) !吊杆节点!jj=jj+1!*enddo!*do,i,1,10,1!ngen,2,6000,i*100+4,i*100+14,10,0,0,z(i)+ygd(i) !吊杆节点!*enddo!jj=1!*do,i,10,17,1!ngen,2,6000,i*100+4,i*100+14,10,0,0,z(i-2*jj)+ygd(i) !吊杆节点!jj=jj+1!*enddo!定义吊杆的参数(每根吊杆均采用109根直径7mm高强钢丝)dgarea=109*3.141593*0.25*0.007*0.007dgyl=0.3*1670*1e6 !吊杆除应力采用0.3倍抗拉强度(n/mm*mm)dgl=dgyl*dgareaedg=1.95e11 !吊杆弹模*dim,L0,array,17 !定义吊杆的初始0应力长度数组*dim,L,array,17 !定义吊杆的受力后长度数组*dim,deltL,array,17 !定义吊杆的受力后长度增量数组*dim,istrn,array,17 !定义吊杆的应力数组*dim,fdg,array,17 !定义在吊杆上施加的力数组*do,i,1,9,1L0(i)=z(i) !0应力时单元长度L0(18-i)=z(i) !吊杆沿中轴对称*enddo*do,i,1,17,1fdg(i)=dgl !各吊杆施加的预应力*enddotype,6mat,6*do,i,1,17,1deltL(i)=fdg(i)*L0(i)/(edg*dgarea) !吊杆受力后的增长量L(i)=L0(i)+deltL(i)istrn(i)=deltL(i)/L(i) !吊杆受力后的应变r,i+1000,dgarea,istrn(i) !定义吊杆的实常数特性*enddo*do,i,1,17,1real,i+1000 !赋予吊杆的实常数特性e,i*100+4+3000,i*100+4+6000 !生成吊杆单元e,i*100+10+4+3000,i*100+10+4+6000*enddocm,dg,elemesel,noneallselreal,1018 !定义一个全为零的实常数便于利用截面特性时,不产生冲突!生成风撑上的节点n,20001,0,21.7,z(11)n,20002,0,28.5,z(12)n,20003,0,42.1,z(13)n,20004,0,48.9,z(14)n,20005,0,59.1,z(15)n,20006,0,65.9,z(15)n,20007,0,76.1,z(14)n,20008,0,82.9,z(13)n,20009,0,96.5,z(12)n,20190,0,103.3,z(11)n,20191,12.5,21.7,z(11)n,20192,12.5,28.5,z(12)n,20193,12.5,42.1,z(13)n,20194,12.5,48.9,z(14)n,20195,12.5,59.1,z(15)+ygd(18)n,20196,12.5,65.9,z(15)+ygd(19)n,20197,12.5,76.1,z(14)n,20198,12.5,82.9,z(13)n,20199,12.5,96.5,z(12)n,20190,12.5,103.3,z(11)n,20191,-12.5,21.7,z(11)n,20192,-12.5,28.5,z(12)n,20193,-12.5,42.1,z(13)n,20194,-12.5,48.9,z(14)n,20195,-12.5,59.1,z(15)+ygd(18)n,20196,-12.5,65.9,z(15)+ygd(19)n,20197,-12.5,76.1,z(14)n,20198,-12.5,82.9,z(13)n,20199,-12.5,96.5,z(12)n,20190,-12.5,103.3,z(11)n,20000,0,62.5,z(15) !中撑中点!风撑1(倒k撑)type,5mat,5sectype,500,beam,ctube,windc !定义风撑形状为空心圆secdata,0.486,0.50secnum,500 !调用截面!生成风撑单元!风撑1(倒k撑)e,20002,10404,90000e,20002,10414,90000e,20002,10304,90000e,20002,10314,90000!风撑2(倒k撑)e,20004,10704,90000e,20004,10714,90000e,20004,10604,90000e,20004,10614,90000!风撑3(中撑x撑)e,20000,20195,90000e,20000,20196,90000e,20000,20195,90000e,20000,20196,90000!风撑4(倒k撑)e,20007,11104,90000e,20007,11114,90000e,20007,11204,90000e,20007,11214,90000!风撑5(倒k撑)e,20009,11404,90000e,20009,11414,90000e,20009,11504,90000e,20009,11514,90000esel,s,real,,1018cm,fc,elemesel,none!连接拱轴线!定义拱轴上的坐标点n,30000,12.5 !左拱轴起点n,30001,12.5,125,0 !左拱轴端点n,30002,12.5,2.025,z(16) !左横梁起点n,30003,12.5,125-2.025,z(16) !左横梁端点n,40000,-12.5 !右拱轴起点n,40001,-12.5,125,0 !右拱轴端点n,40002,-12.5,2.025,z(16) !右横梁起点n,40003,-12.5,125-2.025,z(16) !右横梁端点*dim,zgzx,array,23 !定义拱轴线,左轴线*dim,ygzx,array,23 !定义拱轴线,右轴线zgzx(1)=30000 !以下为拱轴线上的节点号zgzx(2)=30002*do,i,3,10,1zgzx(i)=10004+(i-2)*100*enddozgzx(11)=20195zgzx(12)=10904zgzx(13)=20196*do,i,14,21,1zgzx(i)=11004+(i-14)*100*enddozgzx(22)=30003zgzx(23)=30001ygzx(1)=40000ygzx(2)=40002*do,i,3,10,1ygzx(i)=10014+(i-2)*100*enddoygzx(11)=20195ygzx(12)=10914ygzx(13)=20196*do,i,14,21,1ygzx(i)=11014+(i-14)*100*enddoygzx(22)=40003ygzx(23)=40001n,90001,12.5,0,30 !钢管左轴线参考点n,90002,-12.5,0,30 !钢管右轴线参考点type,1mat,1r,2000secnum,1*do,i,1,22 !生成拱轴单元(钢管)e,zgzx(i),zgzx(i+1),90001e,ygzx(i),ygzx(i+1),90002*enddocm,gzgg,elemesel,nonetype,2mat,2secnum,2*do,i,1,22e,zgzx(i),zgzx(i+1),90001 !生成拱轴单元(混凝土)e,ygzx(i),ygzx(i+1),90002*enddocm,gzhnt,elemesel,none!建立拱轴底座type,9r,6real,6mat,9sectype,70,beam,rect,dizuo !定义截面形状secdata,1.4,4secnum,70e,30001,30003,90001e,40001,40003,90002e,30000,30002,90001e,40000,40002,90002cm,gzdz,elemesel,none!建立盖梁ngen,2,2000,1,16,1,0,-2.025,0 !复制节点到拱轴线起点位置ngen,2,100,1801,1816,1,0,2.025,0 !复制节点到拱轴线终点位置type,7mat,7r,7real,7sectype,60,beam,rect,gailsecdata,4,1.6secnum,60n,90005,0,0,20n,90006,0,125,20e,1901,1902,90006 !盖梁单元125米处e,1902,1903,90006e,1903,1904,90006e,1904,1905,90006e,1905,1906,90006e,1901,1912,90006e,1912,1913,90006e,1913,1914,90006e,1914,1915,90006e,1915,1916,90006e,2019,2019,90005e,2019,2019,90005e,2019,2019,90005e,2019,2019,90005e,2019,2019,90005e,2019,2019,90005e,2019,2019,90005e,2019,2019,90005e,2019,2019,90005cm,gailiang,elemesel,none!建立墩柱n,90003,14.5 ! 参考点n,90004,-14.5type,8mat,8r,8real,8sectype,50,beam,hrec,dun !定义截面形状secdata,3,4,0.8,0.8,0.8,0.8ngen,2,10000,2019,2019,10,0,0,-10ngen,2,10000,1904,1914,10,0,0,-10 !墩柱节点,距离桥面-10m secnum,50e,2019,12019,90003e,2019,12019,90004e,1904,11904,90003e,1914,11914,90004cm,dunzhu,elemesel,none!承台单元type,11mat,11r,1sectype,40,beam,rect,chengtaisecdata,12.7,9r,11real,11secnum,40secdata,3,4,0.8,0.8,0.8,0.8ngen,2,20000,2019,2019,10,0,0,-11.5ngen,2,20000,1904,1914,10,0,0,-11.5ngen,2,30000,2019,2019,10,0,0,-13ngen,2,30000,1904,1914,10,0,0,-13e,12019,22019,90003 !生成承台单元e,12019,22019,90004e,11914,21914,90004e,22019,32019,90003e,21904,31904,90003e,22019,32019,90004e,21914,31914,90004cm,chtai,elemesel,none!横系梁type,9mat,9r,9real,9n,90007,12.5,0,10 !参考点n,90008,12.5,0,10sectype,90,beam,rect,hxlsecdata,3,3secnum,90e,22019,22019,90007e,21904,21914,90008cm,hxl,elemesel,none!桩n,50000,6.15+1.85,2.7,-13 !定义桩节点n,50001,6.15+1.85,-2.7,-13n,50002,6.15+1.85+4.5,2.7,-13n,50003,6.15+1.85+4.5,-2.7,-13n,50004,6.15+12.7-1.85,2.7,-13n,50005,6.15+12.7-1.85,-2.7,-13nsel,s,,,50000,50005,1nsym,x,10,50000,50005,1allselngen,2,100,50000,50015,1,0,125,0 !首先生成平面内的24个节点!每3m生成一层节点(共8层)ngen,10,1000,50000,51000,1,0,0,-2.5!生成桩单元esel,nonetype,11mat,11r,11real,11sectype,20,beam,csolid,piersecdata,0.9secnum,20*do,i,0,8,1*do,k,0,5,1e,50000+k+i*1000,50000+k+(i+1)*1000,90003e,50000+k+10+i*1000,50000+k+10+(i+1)*1000,90003 e,50000+k+100+i*1000,50000+k+100+(i+1)*1000,90003 e,50000+k+110+i*1000,50000+k+110+(i+1)*1000,90003 *enddo*enddocm,zhuang,elemesel,none!在桩上施加弹簧单元!第一层距离承台-3m 317-340!计算弹簧的弹性系数d=1.8kf=0.9k0=1+1/db1=0.6 !和每排桩有关的系数,当为2排时为0.6 L1=3.6 !桩间净距h1=3*(d+1) !桩在地面或冲刷线下的计算深度k=b1+(1-b1)/0.6*L1/h1B=kf*k0*k*d!第1层土粉质量粘土m=8厚度5!第2层土粉粘土(细砂)m=15厚度5!第3层土中砂m=20厚度5!第4层土粘土m=25厚度5!第5层土基岩c0=300厚度10*dim,m,array,5m(1)=8m(2)=15m(3)=20m(4)=25m(5)=4000 !基岩系数*dim,h,array,5h(1)=5 !第一层土的厚度h(2)=5h(3)=5h(4)=5h(5)=10*dim,c,array,4c(1)=m(1)*h(1)c(2)=m(2)*(h(1)+h(2))c(3)=m(3)*(h(1)+h(2)+h(3))c(4)=m(4)*(h(1)+h(2)+h(3)+h(4))*dim,e,array,4 !弹簧弹性系数e(1)=0.5*c(1)*h(1)*B*1E6*do,i,2,4,1e(i)=0.5*(c(i-1)+c(i))*h(i)*B*1E6*enddoet,20,combin14 !弹簧单元r,20,e(1) !第一层土的弹性系数r,30,e(2) !第二层土的弹性系数r,40,e(3) !第三层土的弹性系数r,50,e(4) !第四层土的弹性系数r,60,m(5)*3.1415926*B*B/4*1E6 !基岩弹性系数!第1层底部距离承台-5mngen,2,10000,51000,51115,1,1,0,0 !x方向ngen,2,20000,51000,51115,1,1,0,0 !y方向type,20real,20*do,i,0,5,1e,51000+i,51000+i+10000e,51000+i+10,51000+i+10+10000e,51100+i,51100+i+10000e,51100+i+10,51100+i+10+10000e,51000+i,51000+i+20000e,51000+i+10,51000+i+10+20000e,51100+i,51100+i+20000e,51100+i+10,51100+i+10+20000*enddocm,spring1,elemesel,none!第2层底部距离承台-10mngen,2,10000,53000,53115,1,1,0,0 !x方向ngen,2,20000,53000,53115,1,1,0,0 !y方向type,20real,30*do,i,0,5,1e,53000+i,53000+i+10000e,53000+i+10,53000+i+10+10000e,53100+i,53100+i+10000e,53100+i+10,53100+i+10+10000e,53000+i,53000+i+20000e,53000+i+10,53000+i+10+20000e,53100+i,53100+i+20000e,53100+i+10,53100+i+10+20000*enddocm,spring2,elemesel,none!第3层底部距离承台-15mngen,2,10000,55000,55115,1,1,0,0 !x方向ngen,2,20000,55000,55115,1,1,0,0 !y方向type,20real,40*do,i,0,5,1e,55000+i,55000+i+10000e,55000+i+10,55000+i+10+10000e,55100+i,55100+i+10000e,55100+i+10,55100+i+10+10000e,55000+i,55000+i+20000e,55000+i+10,55000+i+10+20000e,55100+i,55100+i+20000e,55100+i+10,55100+i+10+20000*enddocm,spring3,elemesel,none!第4层底部距离承台-20mngen,2,10000,57000,57115,1,1,0,0 !x方向ngen,2,20000,57000,57115,1,1,0,0 !y方向type,20real,50*do,i,0,5,1e,57000+i,57000+i+10000e,57000+i+10,57000+i+10+10000e,57100+i,57100+i+10000e,57100+i+10,57100+i+10+10000e,57000+i,57000+i+20000e,57000+i+10,57000+i+10+20000e,57100+i,57100+i+20000e,57100+i+10,57100+i+10+20000*enddocm,spring4,elemesel,none!第5层底部距离承台-22.5mngen,2,10000,58000,58115,1,1,0,0 !x方向ngen,2,20000,58000,58115,1,1,0,0 !y方向type,20real,60*do,i,0,5,1e,58000+i,58000+i+10000e,58000+i+10,58000+i+10+10000e,58100+i,58100+i+10000e,58100+i+10,58100+i+10+10000e,58000+i,58000+i+20000e,58000+i+10,58000+i+10+20000e,58100+i,58100+i+20000e,58100+i+10,58100+i+10+20000*enddocm,spring5,elemesel,none!第6层底部距离承台-25mngen,2,10000,59000,59115,1,1,0,0 !x方向ngen,2,10000,59000,59115,1,1,0,0ngen,2,20000,59000,59115,1,0,1,0 !y方向ngen,2,20000,59000,59115,1,0,1,0ngen,2,30000,59000,59115,1,0,0,-1 !z方向ngen,2,30000,59000,59115,1,0,0,-1type,20real,60*do,i,0,5,1e,59000+i,59000+i+10000e,59000+i+10,59000+i+10+10000e,59100+i,59100+i+10000e,59100+i+10,59100+i+10+10000e,59000+i,59000+i+20000e,59000+i+10,59000+i+10+20000e,59100+i,59100+i+20000e,59100+i+10,59100+i+10+20000e,59000+i,59000+i+30000e,59000+i+10,59000+i+10+30000e,59100+i,59100+i+30000e,59100+i+10,59100+i+10+30000*enddocm,spring6,elemallsel!剔除不需要进行对称复制的单元nsel,s,loc,z,-50,-0.01nsel,r,loc,y,-10,10cm,zdnode,nodeesln,s,,allcm,zdelem,elemallselnsel,u,,,zdnodensym,y,100000,all !节点对称allselesel,u,,,zdelemesym,y,100000,allallsel!定义系杆单元(全桥通长)!定义系杆的参数(每根系杆都采用8束27×7×d5高强钢丝)xgarea=8*27*7*3.1415926*0.25*0.005*0.005 !系杆集束面积exg=1.95e11*dim,fxg,array,2*dim,Llast,array,2*dim,Listrn,array,2*dim,detL,array,2xgyl=0.3*1860*1e6xgl=xgyl*xgareafxg(1)=xgl !施加的预应力fxg(2)=xgl !施加的预应力Lzero=250 !0应力时单元长度*do,i,1,2,1detL(i)=Lzero*fxg(i)/(exg*xgarea)Llast(i)=Lzero+detL(i) !单元最终长度Listrn(i)=detl(i)/Llast(i)r,i+100,xgarea,Listrn(i)*enddotype,10mat,10real,101e,30001,130001 !系杆1 250mreal,102e,40001,140001 !系杆2 250mcm,xigan,elemesel,noneallsel!将盖梁和拱轴交界点进行节点耦合,放松uy,和rotxcp,1,ux,2019,30000cp,1,uy,2019,30000cp,next,uz,2019,30000cp,next,roty,2019,30000cp,next,rotz,2019,30000cp,next,ux,1904,30001cp,next,uy,1904,30001cp,next,uz,1904,30001cp,next,roty,1904,30001cp,next,rotz,1904,30001cp,next,ux,2019,40000cp,next,uy,2019,40000cp,next,uz,2019,40000cp,next,roty,2019,40000cp,next,rotz,2019,40000cp,next,ux,1914,40001cp,next,uy,1914,40001cp,next,uz,1914,40001cp,next,roty,1914,40001cp,next,rotz,1914,40001cp,next,ux,101914,140001cp,next,uy,101914,140001cp,next,uz,101914,140001cp,next,roty,101914,140001cp,next,rotz,101914,140001cp,next,ux,101904,130001cp,next,uy,101904,130001cp,next,uz,101904,130001cp,next,roty,101904,130001cp,next,rotz,101904,130001allselnsel,s,,,30000,130000,100000nsel,a,,,40000,140000,100000nummrg,nodeallsel*do,i,1,3,1cp,next,all,1800+i,1900+i !盖梁和0号横梁耦合*enddoallsel*do,i,5,6,1cp,next,all,1800+i,1900+i*enddoallsel*do,i,1,2,1cp,next,all,1800+i+11,1900+i+11cp,next,all,i+11,2000+i+11,100000+i+11*enddoallsel*do,i,4,5,1cp,next,all,1800+i+11,1900+i+11cp,next,all,i+11,2000+i+11,100000+i+11*enddoallsel*do,i,1,3,1cp,next,all,i,2000+i,100000+i*enddoallselcp,next,all,i,2000+i,100000+i*enddoallsel*do,i,1,3,1cp,next,all,101800+i,101900+i*enddoallsel*do,i,5,6,1cp,next,all,101800+i,101900+i*enddoallsel*do,i,1,2,1cp,next,all,101800+i+11,101900+i+11 *enddoallsel*do,i,4,5,1cp,next,all,101800+i+11,101900+i+11 *enddoallsel!!桩节点和承台节点耦合nsel,s,loc,z,-13nsel,r,loc,y,-3,3nsel,r,loc,x,8,20cp,next,all,allallselnsel,s,loc,z,-13nsel,r,loc,y,-3,3nsel,r,loc,x,-8,-20cp,next,all,allallselnsel,s,loc,z,-13nsel,r,loc,y,122,128nsel,r,loc,x,-8,-20cp,next,all,allallselnsel,s,loc,z,-13nsel,r,loc,y,122,128nsel,r,loc,x,8,20cp,next,all,allallselnsel,s,loc,z,-13nsel,r,loc,y,-122,-128nsel,r,loc,x,8,20allselnsel,s,loc,z,-13nsel,r,loc,y,-122,-128nsel,r,loc,x,-8,-20cp,next,all,allallsel!将吊杆转动放松*do,i,1,17,1cp,next,ux,i*100+4,i*100+4+3000cp,next,ux,i*100+4+6000,i*100+4+10000cp,next,uy,i*100+4,i*100+4+3000cp,next,uy,i*100+4+6000,i*100+4+10000cp,next,uz,i*100+4,i*100+4+3000cp,next,uz,i*100+4+6000,i*100+4+10000cp,next,ux,i*100+4+10,i*100+4+3000+10cp,next,ux,i*100+4+6000+10,i*100+4+10000+10cp,next,uy,i*100+4+10,i*100+4+3000+10cp,next,uy,i*100+4+6000+10,i*100+4+10000+10cp,next,uz,i*100+4+10,i*100+4+3000+10cp,next,uz,i*100+4+6000+10,i*100+4+10000+10cp,next,ux,i*100+4+100000,i*100+4+3000+100000cp,next,ux,i*100+4+6000+100000,i*100+4+10000+100000cp,next,uy,i*100+4+100000,i*100+4+3000+100000cp,next,uy,i*100+4+6000+100000,i*100+4+10000+100000cp,next,uz,i*100+4+100000,i*100+4+3000+100000cp,next,uz,i*100+4+6000+100000,i*100+4+10000+100000cp,next,ux,i*100+4+100000+10,i*100+4+3000+100000+10cp,next,ux,i*100+4+6000+100000+10,i*100+4+10000+100000+10 cp,next,uy,i*100+4+100000+10,i*100+4+3000+100000+10cp,next,uy,i*100+4+6000+100000+10,i*100+4+10000+100000+10 cp,next,uz,i*100+4+100000+10,i*100+4+3000+100000+10cp,next,uz,i*100+4+6000+100000+10,i*100+4+10000+100000+10 *enddoallsel!将土弹簧的一端约束nsel,s,loc,x,13.5nsel,a,loc,x,9nsel,a,loc,x,18nsel,a,loc,x,-11.5nsel,a,loc,x,-7nsel,a,loc,x,-16nsel,a,loc,y,3.7nsel,a,loc,y,-1.7nsel,a,loc,y,128.7nsel,a,loc,y,123.3nsel,a,loc,y,-128.7nsel,a,loc,y,-123.3nsel,a,loc,y,-128.7nsel,a,loc,y,-123.3nsel,r,loc,z,-40,-13d,all,allallselnsel,s,loc,z,-36.5d,all,all/color,pbak,off/repallselnsel,s,loc,y,0nsel,r,loc,z,0nummrg,node,allesel,noneallselfini!工况1 不包括桥面板时,求横梁和纵梁的总重量/soluallselmpdele,dens,allmp,dens,3,2600mp,dens,4,2600acel,,,9.8pstres,onsolvesave!总重力:0.65064e+08N!工况2 恒载下调节预拱度!首先计算只有拱肋自重作用下拱的变位/soluallselmpdele,dens,3,11,1mp,dens,1,7800mp,dens,2,2600*do,i,1,17,1r,1000+i,dgarea,0*enddor,101,xgarea,0r,102,xgarea,0acel,,,9.8solve!工况3 一期恒载(包括桥面板)作用下的无预应力状态/soluallsel*do,i,1,17,1r,1000+i,dgarea,0*enddor,101,xgarea,0r,102,xgarea,0qmb=(250/1.5)*32/6.8*5*1000*9.8 !桥面板的总重量38431372。
钢桁架ANSYS计算方法
13
-78977
-1.96 F
14
0.0000
0
15
26390
-F
17
26390
0.65 F
18
0.0000
0
19
-78977
-1.96 F
20
-40340
-F
21
0.13902E+06
3.45 F
桁架计算跨度l=4.21m
桁架间隔宽度b=0.7m
载荷q=l×b×Q=242.0578KN
均布载荷q/l=57.496KN/m
节点载荷F=q/6=40.34KN=40340N
求解结果:
ELEM(单元杆)
SMIS1(内力)
内力系数
1
-95681
-2.37 F
2
-95681
-2.37 F
3
-0.16345E+06
5、定义截面积,默认设置1mm2
6、定义杆件实常数的截面积
7、定义各个单元连接关系,注意按照一定顺序
8、进入求解模块
定义力和位移边界条件
9、求解
10、进入后处理模块
11、显示结构变形图
12、用列表显示输出轴力(SMISC,1)
例:相关数据:
各节点位置坐标
最大计算载荷Q=82.13716657KN/㎡
钢桁架计算
目的:计算钢桁架中,每根杆的内力。
方法:利用ANSYS有限元分析软件求解计算。
步骤:
1、进入前处理模块
2、定义单元为二维杆单元LINK1
3、定义材料弹性模量EX=210E9,泊松比0.33(Q235常用值)
Q235弹性模量E:196~216GPa
简单桁架桥梁ansys分析
简单桁架桥梁ansys分析Ansys是一款广泛使用的有限元分析软件,可以用于各种工程结构的分析,包括桁架桥梁。
下面是一个简单的桁架桥梁分析的步骤,使用Ansys进行模拟。
一、建立模型1.创建新的分析:在Ansys中,首先需要创建一个新的分析。
选择适当的分析类型,例如静态分析或动态分析,根据需要进行设置。
2.创建几何体:在Ansys中,可以使用自带的建模工具创建几何体。
对于桁架桥梁,需要创建梁单元和节点。
梁单元用于模拟桥梁的横梁和纵梁,节点用于连接梁单元。
3.定义材料属性:为梁单元分配适当的材料属性,例如弹性模量、泊松比、密度等。
4.网格化:对几何体进行网格化,以生成有限元网格。
可以调整网格密度以获得更精确的结果。
5.边界条件和载荷:定义边界条件和载荷。
对于桁架桥梁,可能需要在支撑处施加固定约束,并在桥面上施加车辆载荷。
二、进行分析1.运行分析:在Ansys中,可以运行分析并观察结果。
可以使用后处理功能来查看结果,例如位移、应力、应变等。
2.检查结果:检查模型的位移、应力、应变等是否符合预期。
如果结果不符合预期,可能需要返回模型进行修正。
三、优化设计1.优化设置:在Ansys中,可以使用优化工具对模型进行优化设计。
设置优化目标,例如最小化总重量或最大化刚度。
2.运行优化:运行优化过程,Ansys将自动调整模型的参数以达到优化目标。
3.检查结果:在优化完成后,检查结果以确保满足设计要求。
四、验证模型1.确认模型的正确性:在完成优化设计后,需要确认模型的正确性。
可以通过与实验数据进行比较,或者与其他分析工具的结果进行比较来验证模型的准确性。
2.进行敏感性分析:可以使用Ansys的敏感性分析功能来确定哪些参数对模型结果影响最大。
这有助于在后续设计中更好地控制这些参数。
3.确认模型的可靠性:确认模型是否符合工程要求和规范。
如果模型满足所有条件,那么可以将其用于实际工程设计。
五、应用模型1.工程设计:在确认模型的正确性和可靠性后,可以将模型应用于实际的工程设计。
简单桁架桥梁ANSYS分析
下面以一个简单桁架桥梁为例,以展示有限元分析的全过程。
背景素材选自位于密执安的"Old North Park Bridge" (1904 - 1988),见图3-22。
该桁架桥由型钢组成,顶梁及侧梁,桥身弦杆,底梁分别采用3种不同型号的型钢,结构参数见表3-6。
桥长L=32m,桥高H=5.5m。
桥身由8段桁架组成,每段长4m。
该桥梁可以通行卡车,若这里仅考虑卡车位于桥梁中间位置,假设卡车的质量为4000kg,若取一半的模型,可以将卡车对桥梁的作用力简化为P1,P2和P3,其中P1= P3=5000 N, P2=10000N,见图3-23。
图3-22位于密执安的"Old North Park Bridge" (1904 - 1988)图3-23桥梁的简化平面模型(取桥梁的一半)表3-6桥梁结构中各种构件的几何性能参数构件惯性矩m4横截面积m2顶梁及侧梁(Beam1) 643.8310m-´322.1910m-´桥身弦梁(Beam2) 61.8710-´31.18510-´底梁(Beam3) 68.4710-´33.03110-´解答以下为基于ANSYS图形界面(Graphic User Interface , GUI)的菜单操作流程。
安全提示:如果聊天中有涉及财产的操作,请一定先核实好友身份。
发送验证问题或点击举报天意11:36:47(1)进入ANSYS(设定工作目录和工作文件)程序→ANSYS →ANSYS Interactive →Working directory(设置工作目录)→Initial jobname(设置工作文件名):TrussBridge →Run →OK(2)设置计算类型ANSYS Main Menu:Preferences…→Structural →OK(3)定义单元类型hhQÆRRN«•QQoomm QM•9NN•}ANSYS Main Menu:Preprocessor →Element Type →Add/Edit/Delete... →Add…→Beam : 2delastic 3 →OK(返回到Element Types窗口)→Close(4)定义实常数以确定梁单元的截面参数ANSYS Main Menu: Preprocessor →Real Constants…→Add/Edit/Delete →Add…→select Type 1Beam 3 →OK →input Real Constants Set No. : 1 , AREA: 2.1 9E-3,Izz: 3.83e-6(1号实常数用于顶梁和侧梁) →Apply →input Real Constants Set No. : 2 , AREA: 1.18 5E-3,Izz: 1.87E-6 (2号实常数用于弦杆)→Apply →input Real Constants Set No. : 3, AREA: 3.031E-3,Izz: 8.47E-6 (3号实常数用于底梁) →OK(back to Real Constants window) →Close (the Real Constants win dow)(5)定义材料参数ANSYS Main Menu: Preprocessor →Material Props →Material Model s →Structural →Linear→Elastic →Isotropic →input EX: 2.1e11, PRXY: 0.3(定义泊松比及弹性模量) →OK →Density (定义材料密度) →input DENS: 7800, →OK →Close(关闭材料定义窗口)(6)构造桁架桥模型生成桥体几何模型ANSYS Main Menu:Preprocessor →Modeling →Create →Keypoints →In Active CS →NPTKeypoint number:1,X,Y,Z Location in active CS:0,0 →Apply →同样输入其余15个特征点坐标(最左端为起始点,坐标分别为(4,0), (8,0), (12,0), (16,0), (20,0), (24,0), (28,0), (32,0), (4,5 .5), (8,5.5),(12,5.5), (16.5.5), (20,5.5), (24,5.5), (28,5.5))→Lines →Lines →Straight Line →依次分别连接特征点→OK网格划分ANSYS Main Menu: Preprocessor →Meshing →Mesh Attributes →P icked Lines →选择桥顶梁及侧梁→OK →select REAL: 1, TYPE: 1 →Apply →选择桥体弦杆→OK →select REAL: 2,TYPE: 1 →Apply →选择桥底梁→OK →select REAL: 3, TYPE:1 →OK →ANSYS Main Menu:Preprocessor →Meshing →MeshTool →位于Size Controls下的Lines:Set →Element Size on Picked→Pick all →Apply →NDIV:1 →OK →Mesh →Lines →Pick all →OK (划分网格)(7)模型加约束ANSYS Main Menu: Solution →Define Loads →Apply →Structural →Displacement →OnNodes →选取桥身左端节点→OK →select Lab2: All DOF(施加全部约束) →Apply →选取桥身右端节点→OK →select Lab2: UY(施加Y方向约束) →OK(8)施加载荷ANSYS Main Menu: Solution →Define Loads →Apply →Structural →Force/Moment →OnKeypoints →选取底梁上卡车两侧关键点(X坐标为12及20)→OK →select Lab: FY,Value: -5000→Apply →选取底梁上卡车中部关键点(X坐标为16)→OK →select Lab: FY,Value: -10000 →OK→ANSYS Utility Menu:→Select →Everything(9)计算分析ANSYS Main Menu:Solution →Solve →Current LS →OK(10)结果显示ANSYS Main Menu:General Postproc →Plot Results →Deformed shape →Def shape only →OK(返回到Plot Results)→Contour Plot →Nodal Solu →DOF Solution, Y-Component of D isplacement→OK(显示Y方向位移UY)(见图3-24(a))定义线性单元I节点的轴力ANSYS Main Menu →General Postproc →Element Table →Define Table →Add →Lab:[bar_I], By sequence num: [SMISC,1] →OK →Close定义线性单元J节点的轴力ANSYS Main Menu →General Postproc →Element Table →Define Table →Add →Lab:[bar_J], By sequence num: [SMISC,1] →OK →Close画出线性单元的受力图(见图3-24(b))ANSYS Main Menu →General Postproc →Plot Results →Contour Plot →Line Elem Res →LabI: [ bar_I], LabJ: [ bar_J], Fact: [1] →OK(11)退出系统ANSYS Utility Menu:File →Exit →Save Everything →OK(a)桥梁中部最大挠度值为0.003 374m (b)桥梁中部轴力最大值为25 380N图3.24桁架桥挠度UY以及单元轴力计算结果!%%%%% [ANSYS算例]3.4.2(2) %%%%% begin %%%%%%!------注:命令流中的符号$,可将多行命令流写成一行------/prep7 !进入前处理/PLOPTS,DATE,0 !设置不显示日期和时间!=====设置单元和材料ET,1,BEAM3 !定义单元类型R,1,2.19E-3,3.83e-6, , , , , !定义1号实常数用于顶梁侧梁R,2,1.185E-3,1.87e-6,0,0,0,0, !定义2号实常数用于弦杆R,3,3.031E-3,8.47E-6,0,0,0,0, !定义3号实常数用于底梁MP,EX,1,2.1E11 !定义材料弹性模量MP,PRXY,1,0.30 !定义材料泊松比MP,DENS,1,,7800 !定义材料密度!-----定义几何关键点K,1,0,0,, $ K,2,4,0,, $ K,3,8,0,, $K,4,12,0,, $K,5,16,0,, $K,6,20,0,, $K,7,24,0,, $K,8,28,0,, $K,9,32,0,, $K,10,4,5.5,,$K,11,8,5.5,, $K,12,12,5.5,, $K,13,16,5.5,, $K,14,20,5.5,, $K,15,24,5.5, , $K,16,28,5.5,,!-----通过几何点生成桥底梁的线L,1,2 $L,2,3 $L,3,4 $L,4,5 $L,5,6 $L,6,7 $L,7,8 $L,8,9!------生成桥顶梁和侧梁的线L,9,16 $L,15,16 $L,14,15 $L,13,14 $L,12,13 $L,11,12 $L,10,11 $L,1,10!------生成桥身弦杆的线L,2,10 $L,3,10 $L,3,11 $L,4,11 $L,4,12 $L,4,13 $L,5,13 $L,6,13 $L,6, 14 $L,6,15 $L,7,15 $L,7,16 $L,8,16!------选择桥顶梁和侧梁指定单元属性LSEL,S,,,9,16,1,LATT,1,1,1,,,,hhQÆRRN«•QQoomm QM•9NN•}!-----选择桥身弦杆指定单元属性LSEL,S,,,17,29,1,LATT,1,2,1,,,,!-----选择桥底梁指定单元属性LSEL,S,,,1,8,1,LATT,1,3,1,,,,!------划分网格AllSEL,all !再恢复选择所有对象LESIZE,all,,,1,,,,,1 !对所有对象进行单元划分前的分段设置LMESH,all !对所有几何线进行单元划分!=====在求解模块中,施加位移约束、外力,进行求解/soluNSEL,S,LOC,X,0 !根据几何位置选择节点D,all,,,,,,ALL,,,,, !对所选择的节点施加位移约束AllSEL,all !再恢复选择所有对象NSEL,S,LOC,X,32 !根据几何位置选择节点D,all,,,,,,,UY,,,, !对所选择的节点施加位移约束ALLSEL,all !再恢复选择所有对象!------基于几何关键点施加载荷FK,4,FY,-5000 $FK,6,FY,-5000 $FK,5,FY,-10000/replot !重画图形Allsel,all !选择所有信息(包括所有节点、单元和载荷等)solve !求解!=====进入一般的后处理模块/post1 !后处理PLNSOL, U,Y, 0,1.0 !显示Y方向位移PLNSOL, U,X, 0,1.0 !显示X方向位移!------显示线单元轴力------ETABLE,bar_I,SMISC, 1ETABLE,bar_J,SMISC, 1PLLS,BAR_I,BAR_J,0.5,1 !画出轴力图finish !结束你参考这个例题试一下。
某钢桥的ansys命令流
1模型简介模型以壳单元shell63为主,横梁采用beam44单元,吊杆采用link10单元,横梁、吊杆与拱肋的连接采用节点耦合的方式。
不考虑纵坡和横披。
111213141516X Y ZROTCPACEL图1 三维图111213141516X Y Z 图2 正立面图2建模命令流finish/clear/config,nproc,2/prep7!双核运行!材料!钢材EX,1,2.06e11PRXY ,1,0.3DENS,1,7800MPDA TA,ALPX,1,,1.2e-5!吊杆钢材EX,2,2.06e11PRXY ,2,0.3DENS,2,7800MPDA TA,ALPX,2,,1.2e-5!桥面钢材包括铺装层重量EX,3,2.06e11PRXY,3,0.3DENS,3,7800*3MPDA TA,ALPX,3,,1.2e-5!单元ET,1,plane82ET,2,shell63!设置shell63关键字KEYOPT,2,1,0KEYOPT,2,2,0KEYOPT,2,3,0KEYOPT,2,5,2KEYOPT,2,6,0KEYOPT,2,7,0KEYOPT,2,8,0KEYOPT,2,9,0KEYOPT,2,11,0et,3,beam189!设置beam189关键字KEYOPT,3,1,1KEYOPT,3,2,0KEYOPT,3,4,2KEYOPT,3,6,3KEYOPT,3,7,2KEYOPT,3,8,3KEYOPT,3,9,3KEYOPT,3,10,0KEYOPT,3,11,0KEYOPT,3,12,0ET,4,solid45et,5,link10r,1,0.010r,2,0.016r,3,0.02r,4,0.04!建立各种截面!钢横梁!HL1/HL3a=1.0b=0.680th=0.02ts=0.016RECTNG,-a/2,a/2,-b/2,b/2, RECTNG,-a/2+ts,a/2-ts,-b/2+th,b/2-th,RECTNG,-a/2+ts,-a/2+ts+0.16,-0.01,0.01, RECTNG,a/2-ts,a/2-ts-0.16,-0.01,0.01, RECTNG,-0.01,0.01,-b/2+0.02,-b/2+0.02+0.16, RECTNG,-0.01,0.01,b/2-0.02-0.16,b/2-0.02, AGLUE,allwprot,0,90,0wpoff,0,0,-b/2+thASBW,allwpoff,0,0,b/2-th-0.01ASBW,allwpoff,0,0,0.02ASBW,allwpoff,0,0,b/2-th-0.01ASBW,allwprot,0,0,90wpoff,0,0,-a/2+tsASBW,allwpoff,0,0,a-2*tsASBW,allwpoff,0,0,-a/2+ts+0.01ASBW,allwpoff,0,0,-0.02ASBW,allESIZE,0.3,0,amesh,allSECWRITE,sec3SECTYPE,3,BEAM,MESH, SECOFFSET,CENT,,,SECREAD,'sec3','SECT',' ',MESHaclear,allADELE,all, , ,1WPCSYS,-1,0!HL2!吊杆SECTYPE, 4, BEAM, CTUBE, , 0 SECOFFSET, CENTSECDA TA,0.032,0.04,,0,0,0,0,0,0,0!建立跨中截面所有线k,1,0,-2.058,0+5.45k,2,0,-2,0+5.45k,,0,-1.6,0+5.45k,,0,-0.8,0+5.45k,,0,-0.4,0+5.45k,,0,-1.6,0.108+5.45k,,0,-1.2,0.108+5.45k,,0,-0.8,0.108+5.45k,,0,-0.4,0.108+5.45k,,0,-2.541661,0.200339+5.45 k,,0,-2.37,0.576+5.45k,,0,-1.6,0.576+5.45k,,0,-1.2,0.576+5.45k,,0,-0.8,0.576+5.45k,,0,-0.4,0.576+5.45k,,0,-2.742,0.684+5.45k,,0,-2.37,0.684+5.45k,,0,-2,0.684+5.45k,,0,-1.6,0.684+5.45k,,0,-1.2,0.684+5.45k,,0,-0.8,0.684+5.45k,,0,-0.4,0.684+5.45k,,0,0,0.684+5.45k,,0,0,0+5.45*do,i,1,5,1l,i,i+1*enddol,6,25*do,i,17,24,1l,i,i+1*enddo*do,i,3,6,1l,i,i+4*enddo*do,i,20,23,1l,i,i-7*enddol,18,12larc,1,17,11l,2,19!建立拖拉的路径线k,26,-40,0,0+5.45l,25,26!拖拉成箱梁*do,i,1,25,1ADRAG,i, , , , , ,26*enddoallselNUMMRG,KP, , , ,LOWwprot,0,0,90*do,i,1,15,1wpoff,0,0,-2.5ASBW,all*enddo!通过循环生成横隔板,选择封闭的为轮廓线来生成!生成矩形部分*do,i,0,16,1lsel,s,loc,x,-2.5*i-0.001,-2.5*i+0.001lsel,r,loc,y,-2.001,0.001lsel,u,loc,z,0.001+5.45,0.2+5.45lsel,u,loc,z,0.683+5.45,0.6+5.45al,all*enddo!生成扇形部分*do,i,0,16,1lsel,s,loc,x,-2.5*i-0.001,-2.5*i+0.001lsel,r,loc,y,-1.999,-2.749lsel,u,loc,z,0.683+5.45,0.6+5.45al,all*enddoallsel!多次分割WPCSYS,-1,0asel,u,loc,z,0.001+5.45,0.2+5.45asel,u,loc,z,0.683+5.45,0.6+5.45wprot,0,90,0*do,i,1,6,1KWPA VE,iasbw,all*enddoKWPA VE,18asbw,allallselWPA VE,0,0,0wpoff,0,0,0.2asbw,allwpoff,0,0,1.58asbw,allallselWPCSYS,-1,0KWPA VE,10asbw,allKWPA VE,16asbw,allWPCSYS,-1,0wpoff,0,0,0.2+5.45asbw,allwpoff,0,0,0.268asbw,all!删除需要空心的地方*do,j,0,3*do,i,1,3,1asel,s,loc,x,-2.5*(i+4*j)-0.001,-2.5*(i+4*j)+0.001 asel,r,loc,z,0.2+5.45,0.2+0.268+5.45asel,r,loc,y,-0.2,-0.2-1.580ADELE,all, , ,1*enddo*enddoallselaglue,allNUMMRG,KP, , , ,LOWNUMCMP,KPNUMCMP,LINENUMCMP,AREAWPCSYS,-1,0!建立壳单元拱肋!上拱肋跨中断面关键点k,10100,0,-3.25-0.5,-0.35+10.1k,10101,0,-3.25-0.5,0+10.1k,10102,0,-3.25-0.5,0.35+10.1k,10103,0,-3.25-0.34,0+10.1k,10104,0,-3.25-0.25,-0.35+10.1k,10105,0,-3.25-0.25,-0.19+10.1k,10106,0,-3.25-0.25,0.19+10.1k,10107,0,-3.25-0.25,0.35+10.1k,10108,0,-3.25-0.16,0+10.1k,10109,0,-3.25+0,-0.35+10.1k,10110,0,-3.25+0,0+10.1k,10111,0,-3.25+0,0.35+10.1k,10112,0,-3.25+0.16,0+10.1k,10113,0,-3.25+0.25,-0.35+10.1k,10114,0,-3.25+0.25,-0.19+10.1k,10115,0,-3.25+0.25,0.19+10.1 k,10116,0,-3.25+0.25,0.35+10.1 k,10117,0,-3.25+0.34,0+10.1 k,10118,0,-3.25+0.5,-0.35+10.1 k,10119,0,-3.25+0.5,0+10.1 k,10120,0,-3.25+0.5,0.35+10.1 !下拱肋跨中断面关键点k,10200,0,-3.25-0.5,-0.35+5 k,10201,0,-3.25-0.5,0+5k,10202,0,-3.25-0.5,0.35+5k,10203,0,-3.25-0.34,0+5k,10204,0,-3.25+0,-0.35+5k,10205,0,-3.25+0,-0.19+5k,10206,0,-3.25+0,0.19+5k,10207,0,-3.25+0,0.35+5k,10208,0,-3.25+0.34,0+5k,10209,0,-3.25+0.5,-0.35+5 k,10210,0,-3.25+0.5,0+5k,10211,0,-3.25+0.5,0.35+5!连成线NUMCMP,LINELSTR, 10118, 10113 LSTR, 10113, 10109 LSTR, 10109, 10104 LSTR, 10104, 10100 LSTR, 10100, 10101 LSTR, 10101, 10102 LSTR, 10102, 10107 LSTR, 10107, 10111 LSTR, 10111, 10116 LSTR, 10116, 10120 LSTR, 10120, 10119 LSTR, 10119, 10118 LSTR, 10111, 10110 LSTR, 10110, 10109 LSTR, 10116, 10115 LSTR, 10113, 10114 LSTR, 10119, 10117 LSTR, 10110, 10112 LSTR, 10107, 10106 LSTR, 10104, 10105 LSTR, 10101, 10103 LSTR, 10110, 10108LSTR, 10209, 10204LSTR, 10204, 10200LSTR, 10200, 10201LSTR, 10201, 10202LSTR, 10202, 10207LSTR, 10207, 10211LSTR, 10211, 10210LSTR, 10210, 10209LSTR, 10210, 10208LSTR, 10201, 10203LSTR, 10204, 10205LSTR, 10207, 10206!通过旋转生成面!建立上拱肋圆心关键点k,10300,0,0,-27.141k,10301,0,-3,-27.141!建立下拱肋圆心关键点k,10302,0,0,-67.525+5k,10303,0,-3,-67.525+5 WPCSYS,-1,0!上拱肋*do,i,2499,2520,1AROTA T, i,,,,,,10300,10301,43.21447, *enddo!下拱肋*do,i,2521,2532,1AROTA T, i,,,,,,10302,10303,22.18745, *enddo!拱脚的处理ASEL,S,LOC,Y,-3+0.251,-10 ALLSEL,BELOW,AREA NUMMRG,KP, , , ,LOW NUMCMP,KPNUMCMP,LINENUMCMP,AREANUMMRG,KP, , , SELE,LOW ASEL,R,LOC,Z,0,4.5wprot,0,90,0wpoff,0,0,3.00asbw,allwpoff,0,0,0.50asbw,allASEL,S,LOC,Y,-3+0.251,-10WPCSYS,-1,0wprot,0,0,90*do,i,1,5,1wpoff,0,0,-3asbw,all*enddoWPA VE,0,0,0*do,i,1,2,1wpoff,0,0,-10asbw,all*enddoasel,r,loc,x,-25,-20AOVLAP,1593,1625KWPLAN,-1, 1441, 1440, 1112ASBW,allasel,r,loc,x,-25,-23.5ADELE,all, , ,1ASEL,S,LOC,Y,-3+0.251,-10ALLSEL,BELOW,AREANUMMRG,KP, , , SELE,LOWWPCSYS,-1,0!拱脚加劲板LSTR, 1409, 1431al,2710,2703,2533LSTR, 1432, 1405al,3177,2536,2992!上拱肋横隔板*do,i,0,5,1lSEL,S,LOC,Y,-3+0.251,-10lsel,r,loc,z,5.001+0.35,10.101+0.35lsel,r,loc,x,-3*i-0.001,-3*i+0.001lsel,u,loc,y,-3.001,-2.999lsel,u,loc,y,-3.501,-3.499lsel,u,loc,z,10.1-(37.241-(37.241**2-(-3*i)**2)**0.5)-0.05,10.1-(37.241-(37.241**2-(-3*i)**2)** 0.5)+0.05lsel,r,loc,y,-2.749,-3.251al,alllSEL,S,LOC,Y,-3+0.251,-10lsel,r,loc,z,5.001+0.35,10.101+0.35lsel,r,loc,x,-3*i-0.001,-3*i+0.001lsel,u,loc,y,-3.001,-2.999lsel,u,loc,y,-3.501,-3.499lsel,u,loc,z,10.1-(37.241-(37.241**2-(-3*i)**2)**0.5)-0.05,10.1-(37.241-(37.241**2-(-3*i)**2)** 0.5)+0.05lsel,r,loc,y,-3.249,-4al,all*enddoi=1lSEL,S,LOC,Y,-3+0.251,-10lsel,r,loc,z,5.001,10.101+0.35lsel,r,loc,x,-10*i-0.001,-10*i+0.001lsel,u,loc,y,-3.001,-2.999lsel,u,loc,y,-3.501,-3.499lsel,u,loc,z,10.1-(37.241-(37.241**2-(-10*i)**2)**0.5)-0.05,10.1-(37.241-(37.241**2-(-10*i)**2) **0.5)+0.05lsel,r,loc,y,-2.749,-3.251al,alllSEL,S,LOC,Y,-3+0.251,-10lsel,r,loc,z,5.001,10.101+0.35lsel,r,loc,x,-10*i-0.001,-10*i+0.001lsel,u,loc,y,-3.001,-2.999lsel,u,loc,y,-3.501,-3.499lsel,u,loc,z,10.1-(37.241-(37.241**2-(-10*i)**2)**0.5)-0.05,10.1-(37.241-(37.241**2-(-10*i)**2) **0.5)+0.05lsel,r,loc,y,-3.249,-4al,alli=2lSEL,S,LOC,Y,-3+0.251,-10lsel,r,loc,z,3.8,10.101+0.35lsel,r,loc,x,-10*i-0.001,-10*i+0.001lsel,u,loc,y,-3.001,-2.999lsel,u,loc,y,-3.501,-3.499lsel,u,loc,z,10.1-(37.241-(37.241**2-(-10*i)**2)**0.5)-0.05,10.1-(37.241-(37.241**2-(-10*i)**2) **0.5)+0.05lsel,r,loc,y,-2.749,-3.251al,alllSEL,S,LOC,Y,-3+0.251,-10lsel,r,loc,z,3.8,10.101+0.35lsel,r,loc,x,-10*i-0.001,-10*i+0.001lsel,u,loc,y,-3.001,-2.999lsel,u,loc,y,-3.501,-3.499lsel,u,loc,z,10.1-(37.241-(37.241**2-(-10*i)**2)**0.5)-0.05,10.1-(37.241-(37.241**2-(-10*i)**2) **0.5)+0.05lsel,r,loc,y,-3.249,-4al,all!下拱肋横隔板*do,i,0,5,1lSEL,S,LOC,Y,-3+0.251,-10lsel,r,loc,z,5-(67.525-(67.525**2-(-3*i)**2)**0.5)-0.6,5-(67.525-(67.525**2-(-3*i)**2)**0.5)+0. 6lsel,r,loc,x,-3*i-0.001,-3*i+0.001lsel,u,loc,y,-3.251,-3.249lsel,u,loc,z,5-(67.525-(67.525**2-(-3*i)**2)**0.5)-0.05,5-(67.525-(67.525**2-(-3*i)**2)**0.5)+ 0.05al,all*enddo*do,i,10,20,10lSEL,S,LOC,Y,-3+0.251,-10lsel,r,loc,z,5-(67.525-(67.525**2-(-i)**2)**0.5)-0.6,5-(67.525-(67.525**2-(-i)**2)**0.5)+0.6 lsel,r,loc,x,-i-0.001,-i+0.001lsel,u,loc,y,-3.251,-3.249lsel,u,loc,z,5-(67.525-(67.525**2-(-i)**2)**0.5)-0.05,5-(67.525-(67.525**2-(-i)**2)**0.5)+0.05 al,all*enddoallsel!删除上拱肋多余部分asel,s,loc,y,-3.251,-10asel,r,loc,x,-12,-50KWPA VE, 1411asbw,allasel,r,loc,z,kz(997),20ADELE,all, , ,1asel,s,loc,y,-3.251,-10asel,r,loc,z,5.5,20ADELE,all, , ,1asel,s,loc,y,-2.749,-10NUMMRG,KP, , , SELE,LOW!建立新吊杆线LSTR, 1174, 1152LSTR, 1215, 1190LSTR, 1252, 1229LSTR, 1289, 1266LSTR, 1326, 1303LSTR, 1363, 1340FLST,5,6,4,ORDE,2FITEM,5,2501FITEM,5,-2506LSEL,R, , ,P51Xcm,diaogan,lineallsel!建立横梁线allselk,10001,0,0,kz(1158)k,10002,0,-2.73,kz(1158) k,10003,-10,0,kz(1387) k,10004,-10,-2.73,kz(1387) k,10005,-20,0,kz(1410) k,10006,-20,-2.73,kz(1410)l,10001,10002l,10003,10004l,10005,10006lsel,s,loc,z,-10,5.01lsel,r,loc,y,0,-2.73 WPCSYS,-1,0KWPA VE, 10001 wprot,0,90,0wpoff,0,0,1.6lsbw,allcm,hengliang,line!划分!划分横梁type,3mat,1real,1SECNUM,3esize,0.5lmesh,hengliang!划分吊杆allselTYPE,5mat,2real,1SECNUM,4esize,20lmesh,diaoganallsel!划分横隔板type,2mat,1*do,i,0,40,2.5asel,s,loc,x,-i+0.001,-i-0.001 asel,r,loc,y,0,-2.75 ALLSEL,BELOW,AREA LESIZE,all,0.25, , , , , , ,1 aclear,allamesh,all*enddo!划分纵隔板allselLESIZE,all, , ,-1, ,1asel,s,loc,y,-0.001,0.001 asel,a,loc,y,-2.001,-1.999 ALLSEL,BELOW,AREA LESIZE,all,0.4, , , , , , ,1 type,2mat,1real,3amesh,all!划分加劲肋asel,s,loc,y,-2.370,-2.370*do,i,0.4,1.6,0.4asel,a,loc,y,-i-0.001,-i+0.001 *enddoALLSEL,BELOW,AREA LESIZE,all, , ,-1, ,1 LESIZE,all,0.4, , , , , , ,1 type,2mat,1real,1amesh,allallsel!划分箱体外包面ASEL,S,MA T,,0asel,r,loc,y,0.001,-2.75+0.001 LESIZE,all, , ,-1, ,1 LESIZE,all,0.4, , , , , , ,1 real,2amesh,allallsel!拱肋划分allselLESIZE,all, , ,-1, ,1mat,1real,3esize,0.5asel,s,loc,y,-2.749,-10 amesh,allallsel!镜像生成全桥LSEL,S,TYPE,,3LSEL,a,TYPE,,5 LSYMM,Y,all, , , ,0,0 asel,s,loc,y,-0.001,-10 ARSYM,Y,all, , , ,0,0LSEL,S,TYPE,,3LSEL,a,TYPE,,5lsel,r,loc,x,-50,-0.001 LSYMM,x,all, , , ,0,0 asel,s,loc,x,-0.001,-50 ARSYM,x,all, , , ,0,0allselNUMMRG,NODE, , , ,LOW NUMMRG,KP, , , ,LOW NUMMRG,NODE,,, ,LOW !统一所有单元坐标allselcsys,0EMODIF,P51X,ESYS, ,!耦合横梁和箱梁!跨中nsel,s,loc,x,-0.001,0.001 nsel,r,loc,y,-1.601,-1.599n1=node(0,-1.60,0)nsel,u,,,n1n2=node(0,-1.60,0)nsel,a,,,n1cp,1,ux,n1,n2cp,next,uy,n1,n2cp,next,uz,n1,n2nsel,s,loc,x,-0.001,0.001 nsel,r,loc,y,1.601,1.599n1=node(0,1.60,0)nsel,u,,,n1n2=node(0,1.60,0)nsel,a,,,n1cp,next,ux,n1,n2cp,next,uz,n1,n2!H2横梁处nsel,s,loc,x,-10.001,-9.999 nsel,r,loc,y,-1.601,-1.599n1=node(0,-1.60,0) nsel,u,,,n1n2=node(0,-1.60,0) nsel,a,,,n1cp,next,uy,n1,n2cp,next,uz,n1,n2nsel,s,loc,x,-10.001,-9.999 nsel,r,loc,y,1.601,1.599n1=node(0,1.60,0)nsel,u,,,n1n2=node(0,1.60,0)nsel,a,,,n1cp,next,uz,n1,n2!H2'横梁处nsel,s,loc,x,10.001,9.999 nsel,r,loc,y,-1.601,-1.599n1=node(0,-1.60,0) nsel,u,,,n1n2=node(0,-1.60,0) nsel,a,,,n1cp,next,uy,n1,n2cp,next,uz,n1,n2nsel,s,loc,x,10.001,9.999 nsel,r,loc,y,1.601,1.599n1=node(0,1.60,0)nsel,u,,,n1n2=node(0,1.60,0)nsel,a,,,n1cp,next,uz,n1,n2!H1横梁处nsel,s,loc,x,-20.001,-19.999 nsel,r,loc,y,-1.601,-1.599n1=node(0,-1.60,0) nsel,u,,,n1n2=node(0,-1.60,0) nsel,a,,,n1cp,next,uy,n1,n2cp,next,uz,n1,n2nsel,s,loc,x,-20.001,-19.999 nsel,r,loc,y,1.601,1.599n1=node(0,1.60,0)nsel,u,,,n1n2=node(0,1.60,0)nsel,a,,,n1cp,next,uz,n1,n2!H1'横梁处nsel,s,loc,x,20.001,19.999 nsel,r,loc,y,-1.601,-1.599n1=node(0,-1.60,0) nsel,u,,,n1n2=node(0,-1.60,0) nsel,a,,,n1cp,next,uy,n1,n2cp,next,uz,n1,n2nsel,s,loc,x,20.001,19.999 nsel,r,loc,y,1.601,1.599n1=node(0,1.60,0)nsel,u,,,n1n2=node(0,1.60,0)nsel,a,,,n1cp,next,uz,n1,n2allsel!耦合横梁与拱肋!跨中横梁n1=node(0,-2.72,5)n2=node(-0.5,-2.75,5)n3=node(0.5,-2.75,5)n4=node(0,-2.75,5)n5=node(0,-2.75,5+0.35) cp,next,all,n1,n4n1=node(0,2.72,5)n2=node(-0.5,2.75,5)n3=node(0.5,2.75,5)n4=node(0,2.75,5)n5=node(0,2.75,5+0.35) cp,next,all,n1,n4!H2n1=node(-10,-2.72,4.255) n2=node(-10.5,-2.75,4.255) n3=node(-9.5,-2.75,4.255)n4=node(-10,-2.75,4.255)n5=node(-10,-2.75,4.255+0.35) cp,next,all,n1,n4n1=node(-10,2.72,4.255)n2=node(-10.5,2.75,4.255)n3=node(-9.5,2.75,4.255)n4=node(-10,2.75,4.255)n5=node(-10,2.75,4.255+0.35) cp,next,all,n1,n4n1=node(10,-2.72,4.255)n2=node(10.5,-2.75,4.255)n3=node(9.5,-2.75,4.255)n4=node(10,-2.75,4.255)n5=node(10,-2.75,4.255+0.35) cp,next,all,n1,n4n1=node(10,2.72,4.255)n2=node(10.5,2.75,4.255)n3=node(9.5,2.75,4.255)n4=node(10,2.75,4.255)n5=node(10,2.75,4.255+0.35) cp,next,all,n1,n4!H1n1=node(-20,-2.72,4.274)n2=node(-20.5,-2.75,4.274-0.35) n3=node(-20.5,-2.75,4.274+0.35) n4=node(-19.5,-2.75,4.274-0.35) n5=node(-19.5,-2.75,4.274+0.35) n6=node(-20,-2.75,4.274)n7=node(-20,-2.75,4.274+0.35) cp,next,all,n1,n6n1=node(-20,2.72,4.274)n2=node(-20.5,2.75,4.274-0.35) n3=node(-20.5,2.75,4.274+0.35) n4=node(-19.5,2.75,4.274-0.35) n5=node(-19.5,2.75,4.274+0.35) n6=node(-20,2.75,4.274)n7=node(-20,2.75,4.274+0.35) cp,next,all,n1,n6n1=node(20,-2.72,4.274)n2=node(20.5,-2.75,4.274-0.35) n3=node(20.5,-2.75,4.274+0.35) n4=node(19.5,-2.75,4.274-0.35) n5=node(19.5,-2.75,4.274+0.35) n6=node(20,-2.75,4.274)n7=node(20,-2.75,4.274+0.35)cp,next,all,n1,n6n1=node(20,2.72,4.274)n2=node(20.5,2.75,4.274-0.35)n3=node(20.5,2.75,4.274+0.35)n4=node(19.5,2.75,4.274-0.35)n5=node(19.5,2.75,4.274+0.35)n6=node(20,2.75,4.274)n7=node(20,2.75,4.274+0.35)cp,next,all,n1,n6! 施加约束nsel,s,loc,z,5,50D,node(-40,-1.60,0), , , , , ,UY,UZ, , , ,D,node(-40,1.60,0), , , , , ,Uz,, , , ,D,node(-30,-1.60,0), , , , , ,UY,UZ, , , ,D,node(-30,1.60,0), , , , , ,Uz,, , , ,D,node(30,-1.60,0), , , , , ,UY,UZ, , , ,D,node(30,1.60,0), , , , , ,Uz,, , , ,D,node(40,-1.60,0), , , , , ,UY,UZ, , , ,D,node(40,1.60,0), , , , , ,Uz,, , , ,nsel,s,loc,x,-25.366,-50nsel,r,loc,z,-10,0.256D,all, , , , , ,all,, , , ,nsel,s,loc,x,25.366,50nsel,r,loc,z,-10,0.256D,all, , , , , ,all,, , , ,allsel!修改桥面单元,使之重量包括铺装层重量asel,s,loc,z,6.133,6.135asel,r,loc,y,2.73,-2.73ALLSEL,BELOW,AREAEMODIF,all,MA T,3,allsel3影响线求解(静力法)!在给一排单元上施加-1000kN/m^2的均布荷载,每排单元在纵向上的长度为2.5/7m。
ansys桥梁受力分析命令流
ANSYS桥梁受力分析命令流ANSYS是一款广泛应用于工程领域的有限元分析软件,它可以帮助工程师、设计师、科学家等进行各种力学或结构分析。
在工程领域中,桥梁是一种比较常见的结构,因此本文将介绍在ANSYS中进行桥梁受力分析的命令流。
1. 建立桥梁模型首先,我们需要建立桥梁模型。
在ANSYS中,可以使用多种模型建立工具,例如现在比较流行的三维CAD软件,如SolidWorks、CATIA等,均可以建立出桥梁模型。
建立好桥梁模型后,需要导入到ANSYS中进行后续分析。
2. 分配材料属性桥梁模型建立完后,需要根据实际情况为它分配材料属性。
以混凝土桥梁为例,我们可以使用“MP(材料属性)”命令为桥梁定义混凝土的相应属性。
比如,定义混凝土的弹性模量、泊松比等。
3. 设定支座约束桥梁的受力分析需要考虑桥梁的支座约束情况。
在ANSYS中,我们可以通过使用“SUPPORT”命令设定不同类型的支座约束。
4. 设定荷载桥梁承受的荷载对于力学分析至关重要,因此在进行桥梁受力分析时,必须设定荷载的种类和大小。
ANSYS提供了多种荷载设定方式,如集中荷载、分布荷载、自重荷载等。
5. 载荷应用载荷应用是桥梁受力分析中的一个关键步骤,它能够准确模拟桥梁所受荷载。
在ANSYS中,可以使用“D”命令定义荷载的应用方式,如定义一个位移荷载、力荷载等。
6. 求解过程ANSYS中提供了多种求解器,可以快速准确地解析桥梁模型受力情况。
在进行桥梁受力分析时,需要选择合适的求解器并进行计算。
ANSYS中提供了两种类型求解器,一种是基于传统矩阵进行计算的Direct Solver求解器,另一种是基于有限元模型优化的Iterative Solver求解器。
7. 结果分析在桥梁受力分析计算完成后,需要对应力和位移等结果进行分析并产生可视化的输出。
在ANSYS中,可以使用“POST1”命令来输出结果,并对结果进行后处理。
8.经过以上步骤进行受力分析后,可以得到桥梁受力的详细情况,也可以对桥梁的结构进行合理的优化。
钢管混凝土拱桥ansys命令流——很实用,很经典
钢管混凝⼟拱桥ansys命令流——很实⽤,很经典fini/clear/title,concrete-filled tube arch bridge /prep7截⾯1et,1,82cyl4,0.6,0.6,0.6,,0.586cyl4,0.6,2.4,0.6,,0.586rectng,0.35,0.364,0.8,2.2rectng,0.85,0.836,0.8,2.2allselaptn,alladele,5,8,1allselaadd,allsmrtsize,5amesh,allsecwrite,gg,sect,,1sectype,1,beam,meshsecoffset,cent,,,secread,'gg','sect',,meshasel,allaclear,alladele,all,,,1截⾯2cyl4,0.6,0.6,0.586cyl4,0.6,2.4,0.586rectng,0.364,0.836,0.8,2.2allselaadd,allsmrtsize,5amesh,allsecwrite,hnt,sect,,1sectype,2,beam,meshsecoffset,cent,,,secread,'hnt','sect',,meshasel,allaclear,alladele,all,,,1建⽴横梁标准截⾯a-a k,1,-0.5k,2,-0.75,0.25k,3,-0.75,0.5k,4,-0.55,0.5k,5,-0.55,0.85k,6,0.55,0.85k,7,0.55,0.5k,8,0.75,0.5k,9,0.75,0.25k,10,0.5,0a,1,2,3,4,5,6,7,8,9,10smrtsize,5amesh,allsecwrite,hl1,sect,,1 sectype,4,beam,mesh secoffset,cent,,, secread,'hl1','sect',,mesh asel,all aclear,alladele,all,,,1建⽴横梁标准截⾯b-bk,1,-0.5k,2,-0.5,1k,3,-0.75,1.25k,4,-0.75,1.5k,5,0.75,1.5k,6,0.75,1.25k,7,0.5,1k,8,0.5a,1,2,3,4,5,6,7,8smrtsize,5amesh,allsecwrite,hl2,sect,,1 sectype,5,beam,meshsecoffset,cent,,, secread,'hl2','sect',,mesh asel,allaclear,alladele,all,,,1建⽴横梁标准截⾯c-ck,1,-0.5k,2,-0.5,1.04k,3,-0.75,1.04+0.25k,4,-0.75,1.04+0.5k,5,-0.55,1.04+0.5k,6,-0.55,1.04+0.5+0.35k,7,0.55,1.04+0.5+0.35k,8,0.55,1.04+0.5k,9,0.75,1.04+0.5k,10,0.75,1.04+0.25k,11,0.5,1.04k,12,0.5k,13,-0.2,0.58k,14,-0.2,1.04+0.5-0.3k,15,0.2,1.04+0.5-0.3k,16,0.2,0.58a,1,2,3,4,5,6,7,8,9,10,11,12 a,13,14,15,16asba,1,2smrtsize,5amesh,allsecwrite,hl3,sect,,1 sectype,6,beam,mesh secoffset,cent,,, secread,'hl3','sect',,mesh asel,all aclear,alladele,all,,,1建⽴横梁标准截⾯d-d ldele,allkdele,allk,1,-0.5k,2,-0.5,1.23k,3,-0.75,1.23+0.25k,4,-0.75,1.23+0.5k,5,-0.55,1.23+0.5k,6,-0.55,1.23+0.5+0.35 k,7,0.55,1.23+0.5+0.35k,8,0.55,1.23+0.5k,9,0.75,1.23+0.5k,10,0.75,1.23+0.25k,11,0.5,1.23k,12,0.5k,13,-0.2,0.38k,14,-0.3,0.38+0.2k,15,-0.3,1.23+0.5-0.1-0.2 k,16,-0.2,1.23+0.5-0.1k,17,0.2,1.23+0.5-0.1k,18,0.3,1.23+0.5-0.1-0.2 k,19,0.3,0.38+0.2k,20,0.2,0.38a,1,2,3,4,5,6,7,8,9,10,11,12 a,13,14,15,16,17,18,19,20 asba,1,2smrtsize,5amesh,allsecwrite,hl4,sect,,1 sectype,7,beam,mesh secoffset,cent,,, secread,'hl4','sect',,mesh asel,all aclear,alladele,all,,,1建⽴横梁标准截⾯e-ek,1,-0.5k,2,-0.5,1.23k,3,-0.75,1.23+0.25k,4,-0.75,1.23+0.5k,5,-0.55,1.23+0.5k,6,-0.55,1.23+0.5+0.35k,7,0.55,1.23+0.5+0.35k,8,0.55,1.23+0.5k,9,0.75,1.23+0.5k,10,0.75,1.23+0.25k,11,0.5,1.23k,12,0.5a,1,2,3,4,5,6,7,8,9,10,11,12smrtsize,5amesh,allsecwrite,hl5,sect,,1sectype,8,beam,meshsecoffset,cent,,,secread,'hl5','sect',,meshasel,allaclear,alladele,all,,,1etdele,1 !!将辅助单元删除***********************!定义所有材料特性et,1,beam44 !!钢管特性mp,ex,1,2.1e11mp,dens,1,7800mp,prxy,1,0.3n,90000,0,0,30 !!参考点et,2,beam44 !!钢管内50#混凝⼟特性mp,ex,2,3.5e10 mp,dens,2,2600mp,prxy,2,0.1667et,3,beam44 !!纵梁30#混凝⼟钢管特性mp,ex,3,3.0e10 mp,dens,3,2600mp,prxy,3,0.1667et,4,beam44 !!横梁30#混凝⼟钢管特性mp,ex,4,3.0e10 mp,dens,4,2600mp,prxy,4,0.1667et,5,beam44 !!风撑特性mp,ex,5,2.1e11mp,dens,5,7800mp,prxy,5,0.3et,6,link10 !!吊杆特性(钢绞线)mp,ex,6,1.9e11mp,dens,6,7800mp,prxy,6,0.3keyopt,6,3,0 !只拉吊杆et,7,beam44 !!盖梁30#混凝⼟特性mp,ex,7,3.0e10mp,dens,7,2600mp,prxy,7,0.1667et,8,beam44 !!墩柱40#混凝⼟特性mp,ex,8,3.3e10mp,dens,8,2600mp,prxy,8,0.1667et,9,beam44 !!横梁30#混凝⼟特性mp,ex,9,3.0e10mp,dens,9,2600mp,prxy,9,0.1667et,10,link10 !!系杆特性(钢绞线)mp,ex,10,2.0e11mp,dens,10,7800mp,prxy,10,0.3keyopt,10,3,0 !只拉系杆et,11,beam44 !!承台桩基基础25#混凝⼟特性mp,ex,11,2.85e10 mp,dens,11,2500mp,prxy,11,0.1667##############################!开始建⽴模型!0号横梁(端横梁)n,1,,2.025n,2,6.25,2.025n,3,6.25+5.1,2.025n,4,6.25+5.1+0.2+1.9/2,2.025 !12.5m位置n,5,6.25+5.1+2.3,2.025 n,6,6.25+5.1+2.3+2.23,2.025nsym,x,10,2,6,1type,4r,1real,1mat,4secnum,7!单元1到10号为0号横梁e,1,2secnum,6e,2,3secnum,5e,3,4e,4,5e,5,6secnum,7e,1,12secnum,6e,12,13secnum,5e,13,14e,14,15secnum,4e,15,16!!1号端横梁egen,2,100,1,10,1,,,,,,0,4.7+1.25/2+1.5/2!!其余端横梁2~17号egen,17,100,11,20,1,,,,,,0,5.3+1.5!!18号端横梁egen,2,100,171,180,1,,,,,,0,4.7+1.25/2+1.5/2 cm,hl,elemesel,none!!********************************!纵梁单元type,3r,2real,2mat,3sectype,10,beam,rect,zlsecdata,0.4,1.2secnum,10!0号纵梁e,1,101e,2,102e,3,103e,5,105e,12,112e,13,113!中间纵梁1~16号*do,i,0,15,1e,101+i*100,201+i*100e,102+i*100,202+i*100e,104+i*100,204+i*100e,112+i*100,212+i*100e,114+i*100,214+i*100*enddo!17号纵梁(端横梁)e,1701,1801e,1702,1802e,1703,1803e,1705,1805e,1712,1812e,1713,1813e,1715,1815cm,zl,elemesel,none****************************!定义拱轴线*dim,y,array,16 !将纵桥向距离定义为数组*dim,cs,array,16 !定义分值cs=(1-2x/l)*dim,z,array,16 !将拱⾼度定义为数组!吊杆上的点!⾸先定义吊杆的y向(顺桥向)坐标y(1)=8.1y(2)=14.9y(3)=21.7y(4)=28.5y(5)=35.3y(6)=42.1y(7)=48.9y(8)=55.7y(9)=62.5!风撑上的点!⾸先定义风撑的y向(顺桥向)坐标y(12)=28.5y(13)=42.1y(14)=48.9y(15)=59.1!0号横梁的y向坐标y(16)=2.025fl=25 !⽮⾼ml=1.1k=log(ml+sqrt(ml*ml-1)) !k=ch-1(m)*do,i,1,16,1cs(i)=1-y(i)*2/125z(i)=25-fl/(ml-1)*(cosh(k*cs(i))-1) !拱轴⾼度*enddo!上述定义的拱轴线可参照姚玲森《桥梁⼯程》p1299*dim,ygd,array,23 !定义预⼯拱度(是在后续⼯况中反复调算的结果)!初始运⾏中可将预拱度全部设为0ygd(1)=0.4000e-02+1.29e-02ygd(2)=0.10685e-01+2.22e-02ygd(3)=0.19483e-01+3.03e-02ygd(4)=0.29219e-01+3.69e-02ygd(5)=0.39022e-01+4.23e-02ygd(6)=0.48243e-01+4.65e-02ygd(7)=0.56390e-01+4.95e-02ygd(8)=0.63070e-01+5.15e-02ygd(9)=0.68014e-01+5.21e-02ygd(10)=0.71003e-01+5.15e-02ygd(11)=0.71739e-01+4.95e-02ygd(12)=0.69962e-01+4.65e-02ygd(13)=0.65514e-01+4.23e-02ygd(14)=0.58289e-01+3.69e-02ygd(15)=0.48248e-01+3.04e-02ygd(16)=0.35507e-01+2.24e-02ygd(17)=0.20415e-01+1.23e-02ygd(18)=0.65770e-01+0.40234e-01ygd(19)=0.69769e-01+0.41215e-01!⽣成各吊杆的竖向节点(z向)ngen,2,10000,104,114,10,0,0,z(1)+ygd(1) ngen,2,10000,204,214,10,0,0,z(2)+ygd(2) ngen,2,10000,304,314,10,0,0,z(3)+ygd(3) ngen,2,10000,404,414,10,0,0,z(4)+ygd(4) ngen,2,10000,504,514,10,0,0,z(5)+ygd(5) ngen,2,10000,604,614,10,0,0,z(6)+ygd(6) ngen,2,10000,704,714,10,0,0,z(7)+ygd(7) ngen,2,10000,804,814,10,0,0,z(8)+ygd(8) ngen,2,10000,904,914,10,0,0,z(9)+ygd(9) ngen,2,10000,1004,1014,10,0,0,z(8)+ygd(10) ngen,2,10000,1104,1114,10,0,0,z(7)+ygd(11) ngen,2,10000,1204,1214,10,0,0,z(6)+ygd(12) ngen,2,10000,1304,1314,10,0,0,z(5)+ygd(13) ngen,2,10000,1404,1414,10,0,0,z(4)+ygd(14) ngen,2,10000,1504,1514,10,0,0,z(3)+ygd(15) ngen,2,10000,1604,1614,10,0,0,z(2)+ygd(16) ngen,2,10000,1704,1714,10,0,0,z(1)+ygd(17)ngen,2,3000,104,114,10,0,0,0ngen,2,3000,204,214,10,0,0,0ngen,2,3000,304,314,10,0,0,0ngen,2,3000,404,414,10,0,0,0ngen,2,3000,504,514,10,0,0,0ngen,2,3000,604,614,10,0,0,0ngen,2,3000,704,714,10,0,0,0ngen,2,3000,804,814,10,0,0,0ngen,2,3000,904,914,10,0,0,0ngen,2,3000,1004,1014,10,0,0,0ngen,2,3000,1104,1114,10,0,0,0ngen,2,3000,1204,1214,10,0,0,0ngen,2,3000,1304,1314,10,0,0,0ngen,2,3000,1404,1414,10,0,0,0ngen,2,3000,1504,1514,10,0,0,0ngen,2,3000,1604,1614,10,0,0,0ngen,2,3000,1704,1714,10,0,0,0ngen,2,6000,104,114,10,0,0,z(1)+ygd(1)ngen,2,6000,204,214,10,0,0,z(2)+ygd(2)ngen,2,6000,304,314,10,0,0,z(3)+ygd(3)ngen,2,6000,404,414,10,0,0,z(4)+ygd(4)ngen,2,6000,504,514,10,0,0,z(5)+ygd(5)ngen,2,6000,604,614,10,0,0,z(6)+ygd(6)ngen,2,6000,704,714,10,0,0,z(7)+ygd(7)ngen,2,6000,804,814,10,0,0,z(8)+ygd(8)ngen,2,6000,904,914,10,0,0,z(9)+ygd(9)ngen,2,6000,1004,1014,10,0,0,z(8)+ygd(10)ngen,2,6000,1104,1114,10,0,0,z(7)+ygd(11)ngen,2,6000,1204,1214,10,0,0,z(6)+ygd(12)ngen,2,6000,1304,1314,10,0,0,z(5)+ygd(13)ngen,2,6000,1404,1414,10,0,0,z(4)+ygd(14)ngen,2,6000,1504,1514,10,0,0,z(3)+ygd(15)ngen,2,6000,1604,1614,10,0,0,z(2)+ygd(16)ngen,2,6000,1704,1714,10,0,0,z(1)+ygd(17)!*do,i,1,9,1!ngen,2,10000,i*100+4,i*100+14,10,0,0,z(i)+ygd(i) !吊杆节点!*enddo!!jj=1!*do,i,10,17,1!ngen,2,10000,i*100+4,i*100+14,10,0,0,z(i-2*jj)+ygd(i) !吊杆节点!jj=jj+1!*enddo!!为了对以后吊杆转动进⾏放松,重复吊杆节点!*do,i,1,9,1!ngen,2,3000,i*100+4,i*100+14,10,0,0,z(i)+ygd(i) !吊杆下节点!*enddo!!jj=1!*do,i,10,17,1!ngen,2,3000,i*100+4,i*100+14,10,0,0,z(i-2*jj)+ygd(i) !吊杆节点!jj=jj+1!*enddo!!*do,i,1,10,1!ngen,2,6000,i*100+4,i*100+14,10,0,0,z(i)+ygd(i) !吊杆节点!*enddo!!jj=1!*do,i,10,17,1!ngen,2,6000,i*100+4,i*100+14,10,0,0,z(i-2*jj)+ygd(i) !吊杆节点!jj=jj+1!*enddo!定义吊杆的参数(每根吊杆均采⽤109根直径7mm⾼强钢丝)dgarea=109*3.141593*0.25*0.007*0.007 dgyl=0.3*1670*1e6 !吊杆除应⼒采⽤0.3倍抗拉强度(n/mm*mm)dgl=dgyl*dgareaedg=1.95e11 !吊杆弹模*dim,L0,array,17 !定义吊杆的初始0应⼒长度数组*dim,L,array,17 !定义吊杆的受⼒后长度数组*dim,deltL,array,17 !定义吊杆的受⼒后长度增量数组*dim,istrn,array,17 !定义吊杆的应⼒数组*dim,fdg,array,17 !定义在吊杆上施加的⼒数组*do,i,1,9,1L0(i)=z(i) !0应⼒时单元长度L0(18-i)=z(i) !吊杆沿中轴对称*enddo*do,i,1,17,1fdg(i)=dgl !各吊杆施加的预应⼒*enddotype,6mat,6*do,i,1,17,1deltL(i)=fdg(i)*L0(i)/(edg*dgarea) !吊杆受⼒后的增长量L(i)=L0(i)+deltL(i)istrn(i)=deltL(i)/L(i) !吊杆受⼒后的应变r,i+1000,dgarea,istrn(i) !定义吊杆的实常数特性*enddo*do,i,1,17,1real,i+1000 !赋予吊杆的实常数特性e,i*100+4+3000,i*100+4+6000 !⽣成吊杆单元e,i*100+10+4+3000,i*100+10+4+6000*enddocm,dg,elemesel,noneallsel!!*********************************real,1018 !定义⼀个全为零的实常数便于利⽤截⾯特性时,不产⽣冲突!!********************************* !⽣成风撑上的节点n,20001,0,21.7,z(11)n,20002,0,28.5,z(12)n,20003,0,42.1,z(13)n,20004,0,48.9,z(14)n,20005,0,59.1,z(15)n,20006,0,65.9,z(15)n,20007,0,76.1,z(14)n,20008,0,82.9,z(13)n,20009,0,96.5,z(12)n,20010,0,103.3,z(11)n,20011,12.5,21.7,z(11)n,20012,12.5,28.5,z(12)n,20013,12.5,42.1,z(13)n,20014,12.5,48.9,z(14)n,20015,12.5,59.1,z(15)+ygd(18)n,20016,12.5,65.9,z(15)+ygd(19)n,20017,12.5,76.1,z(14)n,20018,12.5,82.9,z(13)n,20019,12.5,96.5,z(12)n,20020,12.5,103.3,z(11)n,20021,-12.5,21.7,z(11)n,20022,-12.5,28.5,z(12)n,20023,-12.5,42.1,z(13)n,20024,-12.5,48.9,z(14)n,20025,-12.5,59.1,z(15)+ygd(18)n,20026,-12.5,65.9,z(15)+ygd(19)n,20027,-12.5,76.1,z(14)n,20028,-12.5,82.9,z(13)n,20029,-12.5,96.5,z(12)n,20030,-12.5,103.3,z(11)n,20000,0,62.5,z(15) !中撑中点!风撑1(倒k撑)type,5mat,5sectype,500,beam,ctube,windc !定义风撑形状为空⼼圆secdata,0.486,0.50secnum,500 !调⽤截⾯!⽣成风撑单元!风撑1(倒k撑)e,20002,10404,90000e,20002,10414,90000e,20002,10304,90000e,20002,10314,90000!风撑2(倒k撑)e,20004,10704,90000e,20004,10714,90000e,20004,10604,90000e,20004,10614,90000!风撑3(中撑x撑)e,20000,20015,90000e,20000,20016,90000e,20000,20025,90000e,20000,20026,90000!风撑4(倒k撑)e,20007,11104,90000e,20007,11114,90000e,20007,11204,90000e,20007,11214,90000!风撑5(倒k撑)e,20009,11404,90000e,20009,11414,90000e,20009,11504,90000e,20009,11514,90000esel,s,real,,1018cm,fc,elemesel,none!!******************************************* !连接拱轴线!定义拱轴上的坐标点n,30000,12.5 !左拱轴起点n,30001,12.5,125,0 !左拱轴端点n,30002,12.5,2.025,z(16) !左横梁起点n,30003,12.5,125-2.025,z(16) !左横梁端点n,40000,-12.5 !右拱轴起点n,40001,-12.5,125,0 !右拱轴端点n,40002,-12.5,2.025,z(16) !右横梁起点n,40003,-12.5,125-2.025,z(16) !右横梁端点*dim,zgzx,array,23 !定义拱轴线,左轴线*dim,ygzx,array,23 !定义拱轴线,右轴线zgzx(1)=30000 !以下为拱轴线上的节点号zgzx(2)=30002 *do,i,3,10,1zgzx(i)=10004+(i-2)*100*enddozgzx(11)=20015zgzx(12)=10904zgzx(13)=20016*do,i,14,21,1zgzx(i)=11004+(i-14)*100*enddozgzx(22)=30003zgzx(23)=30001ygzx(1)=40000ygzx(2)=40002*do,i,3,10,1ygzx(i)=10014+(i-2)*100*enddoygzx(11)=20025ygzx(12)=10914ygzx(13)=20026*do,i,14,21,1ygzx(i)=11014+(i-14)*100*enddoygzx(22)=40003ygzx(23)=40001n,90001,12.5,0,30 !钢管左轴线参考点n,90002,-12.5,0,30 !钢管右轴线参考点type,1mat,1secnum,1*do,i,1,22 !⽣成拱轴单元(钢管)e,zgzx(i),zgzx(i+1),90001e,ygzx(i),ygzx(i+1),90002*enddocm,gzgg,elemesel,nonetype,2mat,2secnum,2*do,i,1,22e,zgzx(i),zgzx(i+1),90001 !⽣成拱轴单元(混凝⼟)e,ygzx(i),ygzx(i+1),90002*enddocm,gzhnt,elemesel,none*****************************************!建⽴拱轴底座type,9r,6real,6mat,9sectype,70,beam,rect,dizuo !定义截⾯形状secdata,1.4,4secnum,70e,30001,30003,90001e,40001,40003,90002e,30000,30002,90001e,40000,40002,90002cm,gzdz,elemesel,none*****************************************!建⽴盖梁ngen,2,2000,1,16,1,0,-2.025,0 !复制节点到拱轴线起点位置ngen,2,100,1801,1816,1,0,2.025,0 !复制节点到拱轴线终点位置type,7mat,7r,7sectype,60,beam,rect,gail secdata,4,1.6secnum,60n,90005,0,0,20n,90006,0,125,20e,1901,1902,90006 !盖梁单元125⽶处e,1902,1903,90006e,1903,1904,90006e,1904,1905,90006e,1905,1906,90006e,1901,1912,90006e,1912,1913,90006e,1913,1914,90006e,1914,1915,90006e,1915,1916,90006e,2001,2002,90005e,2002,2003,90005e,2003,2004,90005e,2004,2005,90005e,2005,2006,90005e,2001,2012,90005e,2012,2013,90005e,2013,2014,90005e,2014,2015,90005e,2015,2016,90005cm,gailiang,elemesel,none!*****************************!建⽴墩柱n,90003,14.5 ! 参考点n,90004,-14.5type,8mat,8r,8sectype,50,beam,hrec,dun !定义截⾯形状secdata,3,4,0.8,0.8,0.8,0.8 ngen,2,10000,2004,2014,10,0,0,-10ngen,2,10000,1904,1914,10,0,0,-10 !墩柱节点,距离桥⾯-10m secnum,50e,2004,12004,90003e,2014,12014,90004e,1904,11904,90003e,1914,11914,90004cm,dunzhu,elemesel,none*****************************************!承台单元type,11mat,11r,1sectype,40,beam,rect,chengtaisecdata,12.7,9r,11real,11secnum,40secdata,3,4,0.8,0.8,0.8,0.8ngen,2,20000,2004,2014,10,0,0,-11.5ngen,2,20000,1904,1914,10,0,0,-11.5ngen,2,30000,2004,2014,10,0,0,-13ngen,2,30000,1904,1914,10,0,0,-13e,12004,22004,90003 !⽣成承台单元e,11904,21904,90003e,12014,22014,90004e,11914,21914,90004e,22004,32004,90003e,21904,31904,90003e,22014,32014,90004e,21914,31914,90004cm,chtai,elemesel,none!!**********************************!横系梁type,9mat,9r,9real,9n,90007,12.5,0,10 !参考点n,90008,12.5,0,10sectype,90,beam,rect,hxlsecdata,3,3secnum,90e,22004,22014,90007e,21904,21914,90008cm,hxl,elemesel,none!!**********************************!桩n,50000,6.15+1.85,2.7,-13 !定义桩节点n,50001,6.15+1.85,-2.7,-13 n,50002,6.15+1.85+4.5,2.7,-13n,50003,6.15+1.85+4.5,-2.7,-13n,50004,6.15+12.7-1.85,2.7,-13n,50005,6.15+12.7-1.85,-2.7,-13nsel,s,,,50000,50005,1nsym,x,10,50000,50005,1allselngen,2,100,50000,50015,1,0,125,0 !⾸先⽣成平⾯内的24个节点!每3m⽣成⼀层节点(共8层)ngen,10,1000,50000,51000,1,0,0,-2.5!⽣成桩单元esel,nonetype,11mat,11r,11real,11sectype,20,beam,csolid,piersecdata,0.9secnum,20*do,i,0,8,1*do,k,0,5,1e,50000+k+i*1000,50000+k+(i+1)*1000,90003e,50000+k+10+i*1000,50000+k+10+(i+1)*1000,90003 e,50000+k+100+i*1000,50000+k+100+(i+1)*1000,90003 e,50000+k+110+i*1000,50000+k+110+(i+1)*1000,90003 *enddo*enddocm,zhuang,elemesel,none***********************************!在桩上施加弹簧单元!第⼀层距离承台-3m 317-340!计算弹簧的弹性系数d=1.8kf=0.9k0=1+1/db1=0.6 !和每排桩有关的系数,当为2排时为0.6L1=3.6 !桩间净距h1=3*(d+1) !桩在地⾯或冲刷线下的计算深度k=b1+(1-b1)/0.6*L1/h1。
桁架的ansys分析
如图结构各杆的弹性模量和横截面积都为, 112A mm100=,求解该结构的节点位移、单元应力以2.9510/E N m=⨯,2及支反力。
1. 进入ansys程序,进入File,Change Jobname, →huangbinqi→OK。
2. 设置计算类型ANSYS Main Menu: Preferences…选Structural → OK3.选择单元类型ANSYS Main Menu: 进入Preprocessor, Element Type,Add/Edit/Delete…→Add…→Link:2D spar 1,点击OK, (返回到Element Types 窗口) ,Close。
4. 定义材料参数ANSYS Main Menu: Preprocessor →Material Props →Material Models→Structural →Linear →Elastic→ Isotropic: EX:2.95e11 (弹性模量),PRXY: 0.3 (泊松比) → OK ,5. 定义实常数以确定单元的截面积ANSYS Main Menu: Preprocessor →Real Constants…→Add/Edit/Delete →Add →Type 1→ OK→Real Constant Set No: 1 (第 1 号实常数), AREA: 1e-4 (单元的截面积) →OK→Close6.创建节点ANSYS Main Menu: Preprocessor →Modeling →Creat→Nodes→ In Active CS;Node number 1 :(0,0,0)→Apply,Node number 2 :(0.4,0,0)→Apply,Node number 3:(0.4,0.3,0)→Apply,Node number 4 :(0,0.3,0)→OK7.生成单元ANSYS Main Menu: Preprocessor → Modeling → Create → Elements→Elem Attributes ,(接受默认值),Auto numbered→Thru nodes→ OK→选择节点1,2→ Apply,创建单元1.选择节点2,3→ Apply,创建单元2.选择节点1,3→ Apply,创建单元3.选择节点3, 4→OK,创建单元4.后显示节点号,单元号:PlotCtrls→Numbering→,将Node Number,ELE/Attrib numbering 打开。
ANSYS桁架分析详细设置
file://C:\Documents and Settings\Administrator\Application Data\Ansys\v145\Mechanic... 2013-5-7 PDF 文件使用 "pdfFactory Pro" 试用版本创建
Project
Page 4 of 14
Coordinate Systems
TABLE 4 Model (B4, C4) > Coordinate Systems > Coordinate System Object Name Global Coordinate System State Fully Defined Definition Type Cartesian Coordinate System ID 0. Origin Origin X 0. mm Origin Y 0. mm Origin Z 0. mm Directional Vectors X Axis Data [ 1. 0. 0. ] Y Axis Data [ 0. 1. 0. ] Z Axis Data [ 0. 0. 1. ]
Geometry
TABLE 2 Model (B4, C4) > Geometry Geometry Fully Defined Definition C:\Documents and Settings\Administrator\Local Settings\Temp\WB_201304170152_5884_2\unsaved_project_files\dp0\Geom\DM\Geom.agdb DesignModeler Millimeters Program Controlled Body Color Bounding Box 1144.4 mm 1398.3 mm 9753.8 mm
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
MODOPT,SUBSP,6,0,100, ,OFF !选择字空间法,提取6阶模态
SUBOPT,8,4,10,0,0,ALL !字空间法设置
REAL,3 ! 选择3号实常数
ESYS,0 ! 单元坐标系
SECNUM,3 ! 选择3号截面
/POST1 !进入普通后处理器
SET,LIST !列表格阶模态
SET,FIRST !读取第一阶模态
E,12,13 !建立单元
TYPE,1 ! 选择1号单元类型
MAT,1 ! 选择1号材料
TSHAP,LINE ! 选择线形单元
E,2,6 !建立单元
E,6,10
E,10,14
E,1,5
E,5,6,10,9
E,9,10,14,13
NSYM,X,14,ALL !所有节点以y0z平面对称
ESYM, ,14,ALL ! 所有单元以y0z平面对称
NUMMRG,ALL, , , ,LOW !合并重复节点单元,编号取较小者
SECOFFSET, CENT ! 截面质心不偏移
SECDATA,0.3,0.3,0.4,0.012,0.012,0.012,0,0,0,0 ! 3号截面参数
R,1,0.0187,0.17e-3,0.54e-3,0.4,0.4,0, !设置1号实常数
PLNSOL, U,SUM, 2,1.0 !显示振型云图
SET,NEXT !读取下一阶模态
REAL,2 ! 选择2号实常数
ESYS,0 ! 单元坐标系
SECNUM,2 ! 选择2号截面
SECOFFSET, CENT ! 截面质心不偏移
SECDATA,0.4,0.4,0.4,0.012,0.012,0.012,0,0,0,0 ! 2号截面参数
SECTYPE, 3, BEAM, I, FU, 0 ! 定义3号工字形截面
E,5,9
E,9,13
E,3,7
E,7,11
E,4,8
E,8,12
E,1,2
E,3,4
E,5,6
E,7,8
E,9,10
E,11,12
E,13,14
TYPE,1 ! 选择1号单元类型
MAT,1 ! 选择1号材料
SOLVE !开始求解
FINISH !结束求解器
TSHAP,LINE ! 选择线形单元
E,3,6 !建立单元
E,6,11
E,4,5
E,5,12
E,2,3
E,1,4
E,6,7
E,5,8
NGEN,2,1,2,10,4, ,16, ,1, !复制节点
NGEN,2,1,3,11,4, , ,-10,1, !复制节点
TYPE,1 !选择1号单元类型
D,ALL, ,0, , , , ,UY,UZ, , , , !约束两个自由度
ALLSEL,ALL !选择所有
MAT,1 !选择1号材料
REAL,1 ! 选择1号实常数
ESYS,0 !单元坐标系
ESYS,0 ! 单元坐标系
TSHAP,QUAD ! 选择四边形单元
E,1,2,6,5 !建立单元
SECOFFSET, CENT !截面质心不偏移
SECDATA,0.4,0.4,0.4,0.016,0.016,0.016,0,0,0,0 !1号截面参数
SECTYPE, 2, BEAM, I, XIANHENG, 0 !定义2号工字形截面
NSEL,S, , ,23,24 !选择节点
D,ALL, ,0, , , ,UX,UY,UZ, , , , !约束3个自由度
NSEL,S, , ,13,14 !选择节点
NUMCMP,ALL !压缩节点单元等编号
FINISH !结束前处理器
/SOL !进入求解器
/TITLE,Truss Bridge Static Analysis !指定标题
/COM,Structural !选择分析类型为结构分析
/PREP7 !进入前处理器
ET,1,BEAM188 !定义1号单元类型
KEYOPT,1,3,2
ET,2,SHELL181 !定义2号单元类型
KEYOPT,1,3,2
SECTYPE, 1, BEAM, I, duan, 0 !定义1号工字形截面
SECNUM,1 !选择1号截面
TSHAP,LINE !选择线形单元
E,11,14 !建立单元
MP,EX,1,2.1E11 !定义1号材料弹性模量
MP,PRXY,1,0.3 !定义1号材料泊松比
MP,DENS,1,7850 !定义1号材料密度
MP,EX,2,3.5E10 !定义2号材料弹性模量
MP,PRXY,2,0.1667 !定义2号材料泊松比
MP,DENS,2,2500 !定义2号材料密度
R,2,0.0141,0.128e-3,0.415e-3,0.4,0.4,0, !设置2号实常数
R,3,0.0117,0.541e-4,0.324e-3,0.3,0.4,0, !设置3号实常数
R,4,0.3, , , , , , !设置4号实常数
N, ,0,0,-5, , , , !建立节点
NGEN,4,4,ALL, , ,12, , ,1, !复制节点
NGEN,2,1,ALL, , , , ,10,1, !复制节点
E,10,11
E,9,12
TYPE,2 ! 选择2号单元类型
MAT,2 ! 选择2号材料
REAL,4 ! 选择4号实常数