铁路大桥ansys建模实例
基于ANSYS 的连续刚构桥分析(操作篇 )
目录一、工程背景 (1)二、工程模型 (1)三、ANSYS分析 (2)(一)前处理 (2)(1)定义单元类型 (2)(2)定义材料属性 (3)(3)建立工程简化模型 (3)(4)有限元网格划分 (5)(二)模态分析 (5)(1)选择求解类型 (5)(2)建立边界条件 (6)(3)输出设置 (6)(4)求解 (6)(5)读取结果 (6)(6)结果分析 (8)(三)结构试验载荷分析 (8)(1)第二跨跨中模拟车载分析 (8)(2)边跨跨中模拟车载分析 (9)四、结果分析与强度校核 (10)(一)结果分析 (10)(二)简单强度校核 (10)参考文献 (11)连续刚构桥分析一、工程背景：随着我国经济的发展，对交通运输的要求也不断提高；高速路，高铁线等遍布全国，这就免不了要架桥修路。
截至2014年年底，我国公路桥梁总数已达75.71万座，4257.89万延米i。
进百万的桥梁屹立在我国交通线上，其安全便是头等大事。
随着交通运输线的再扩大，连续刚构桥跨越能力大，施工难度小，行车舒顺，养护简便，造价较低等优点将被广泛应用。
二、工程模型：现有某预应力混凝土连续刚构桥，桥梁全长为184m，宽13m，其中车行道宽11.5m，两侧防撞栏杆各0.75m主梁采用C50混凝土。
桥梁设计载荷为公路—— 级。
图2-1桥梁侧立面图上部结构为48m+88m+48m三跨预应力混凝土边界面连续箱梁。
箱梁为单箱双室箱形截面，箱梁根部高5m，中跨梁高2.2m，边跨梁端高2.2m。
箱梁顶板宽12.7m,底板宽8.7m,翼缘板悬臂长2.0m，箱梁高度从距墩中心3.0m处到跨中合龙段处按二次抛物线变化。
0号至3号块长3m（4x3m)，4、5号块长3.5m(2x3.5m)，6号块到合龙段长4m(6x4m)，合龙段长2m。
边跨端部设1.5m横隔板，墩顶0号块设两道厚1.2m横隔板。
0号块范围内箱梁底板厚度为0.90m，1号块范围内底板厚度由0.90m线性变化到0.557m，2号块到合龙段范围内底板厚度由0.557m 线性变化到0.3m。
ANSYS桥梁工程应用实例分析(详细)(图文)
本章介绍桥梁结构的模拟分析。
桥梁是一种重要的工程结构，精确分析桥梁结构在各种受力方式下的响应有较大的工程价值。
模拟不同类型的桥梁需要不同的建模方法，分析内容包括静力分析、动荷载响应分析、施工过程分析等等。
在本章中着重介绍桁架桥、刚架桥和斜拉桥三种类型桥梁。
内容 提要 第6章 ANSYS 桥梁工程应用实例分析本章重点结构分析具体步骤结构静力分析 桁架结构建模方法 结构模态分析本章典型效果图6.1 引言ANSYS通用有限元软件在土木工程应用分析中可发挥巨大的作用。
我们用它来分析桥梁工程结构,可以很好的模拟各种类型桥梁的受力、施工工况、动荷载的耦合等。
ANSYS程序有丰富的单元库和材料库,几乎可以仿真模拟出任何形式的桥梁。
静力分析中,可以较精确的反应出结构的变形、应力分布、内力情况等；动力分析中,也可精确的表达结构的自振频率、振型、荷载耦合、时程响应等特性。
利用有限元软件对桥梁结构进行全桥模拟分析,可以得出较准确的分析结果。
本章介绍桥梁结构的模拟分析。
作为一种重要的工程结构,桥梁的精确分析具有较大的工程价值。
桥梁的种类繁多,如梁桥、拱桥、钢构桥、悬索桥、斜拉桥等等,不同类型的桥梁可以采用不同的建模方法。
桥梁的分析内容又包括静力分析、施工过程模拟、动荷载响应分析等。
可以看出桥梁的整体分析过程比较复杂。
总体上来说,主要的模拟分析过程如下:(1) 根据计算数据,选择合适的单元和材料,建立准确的桥梁有限元模型。
(2) 施加静力或者动力荷载,选择适当的边界条件。
(3) 根据分析问题的不同,选择合适的求解器进行求解。
(4) 在后处理器中观察计算结果。
(5) 如有需要,调整模型或者荷载条件,重新分析计算。
桥梁的种类和分析内容众多,不同类型桥梁的的分析过程有所不同,分析侧重点也不一样。
在这里仅仅给出大致的分析过程,具体内容还要看具体实例的情况。
6.2 典型桥梁分析模拟过程6.2.1 创建物理环境建立桥梁模型之前必须对工作环境进行一系列的设置。
茅草街大桥基于ANSYS的空间计算模型
第19卷第2期2003年6月长 沙 交 通 学 院 学 报JOURNAL OF CHAN GSHA COMMUN ICA TIONS UN IV ERSIT YVol.19No.2J une 2003 文章编号:1000-9779(2003)02-0006-05茅草街大桥基于ANSYS的空间计算模型颜东煌,赖敏芝,张克波,李学文(长沙交通学院桥梁与结构工程系,湖南长沙 410076)摘 要:采用ANSYS通用有限元软件,对南县茅草街大桥进行了计算模型简化,实现了整个施工过程的仿真计算。
对于钢管混凝土拱桥整体空间模型的简化提出了一些方法。
关键词:茅草街大桥;ANSYS;空间计算模型中图分类号:U448.52 文献标识码:A Ξ1 大桥概述茅草街大桥位于益阳市茅草街镇[1],是位于省道1831干线上的一座特大型桥,主桥为80+368+ 80m三跨连续自锚中承式钢管混凝土系杆拱桥,其中两边跨半拱为变截面钢筋混凝土箱形拱,桥面宽为16m。
大桥总体布置见图1。
图1 茅草街大桥总体布置图主跨计算跨径L=356m,计算矢高f=71.2m,矢跨比1/5,拱轴系数m=1.543;主拱圈的结构形式是四弦管桁架式结构,截面宽3.2m,拱脚径向的高度最大是8m,拱顶截面高度为4m。
每片主拱肋由4根外径为1000mm的钢管用平连板或平连杆相连,弦杆壁厚有28、22、18mm三种。
在靠近拱脚的Ⅰ～Ⅲ段设置平连板加强,平连板壁厚为14mm,其它段设置平连杆,在有吊杆处的平连杆的管径为<650×16mm,无吊杆处平连杆的管径为<650×10mm。
拱脚处腹杆采用<550×28mm的钢管,其它的腹杆采用<550×12mm或<550×10mm的空钢管;主跨以路中为对称中心,两边各设置了6道横撑,在拱顶与拱脚处共设置3组“K”字钢管桁架式横撑,中间设置3组“米”字钢管桁架式横撑;纵梁与横梁采用工字形钢联结,吊杆采用73<s7高强钢丝,上端锚固在平连杆上,下端锚固在钢横梁的下翼缘;全桥系杆采用31<j15.24的钢铰线共24束,系杆锚固在两边跨的端横梁上部,是一种自锚式结构。
基于ansys的铁路钢桁架桥受力分析1
第一章工程简介1.1工程概况一、结构设计本工程为单线铁路刚桁架桥，铁路线穿过山区桁架桥跨越峡谷，设计桥梁为简支栓焊桁架桥。
整座桥梁主要由桁架构成的桥跨结构和桥墩、桥台组成，为下承式桥。
该桥的跨度为64m，两侧桥台为重力式，其桥跨的的布置及各种杆件的尺寸如图1所示：图1 桥跨的布置以及各种杆件的基本尺寸注：桥跨的单位为m，各种杆件的截面尺寸单位为mm；所有杆件皆为H型截面，所用钢材为16Mnq;图中1为上下弦杆，2为端斜杆，3为直腹杆，4为斜腹杆。
桁架由上弦、下弦、腹杆组成；腹杆的形式又分为斜腹杆、直腹杆；本设计中上下弦杆，端斜杆，直腹杆杆，斜腹杆分别采用不同同的工字钢。
由于杆件本身长细比较大，虽然杆件之间的连接可能是“固接”，但是实际杆端弯矩一般都很小，因此，设计分析时可以简化为“铰接”。
简化计算时，杆件都是“二力杆”，只承受压力或者拉力，而不产生弯矩和剪力。
由于桥梁跨度较大，而单榀的桁架“平面外”的刚度较弱，因此，“平面外”需要设置支撑，设计桥梁时，“平面外”一般也是设计成桁架形式，这样，桥梁就形成双向都有很好刚度的整体。
有些桥梁桥面设置在上弦，因此力主要通过上弦传递；也有的桥面设置在下弦，由于平面外刚度的要求，上弦之间仍需要连接以减少上弦平面外计算长度。
桁架的弦杆在跨中部分受力比较大，向支座方向逐步减小；而腹杆的受力主要在支座附件最大，在跨中部分腹杆的受力比较小，甚至有理论上的“零杆”。
二、施工方法设计采用拖拉法架设简支栓焊桁架桥，其施工步骤为：1.在桥头路基上布置滑道，其滑道的数量要进行计算确定；2.在滑道上拼装钢梁；3.按照牵引设备进行技术检查，合格后开始拖拉钢梁只预定的桥孔位置；4.钢梁降落就位。
技术检查内容包括:钢梁的拼装质量、钢梁的拱度是否满足设计要求；加固杆件的数量、位置和质量是否满足要求；钢梁的中心位置和标高、滑道的设置、牵引动力的配置情况、落梁的设备、信号和照明、施工流程及人身设备安全。
双线下承式简支钢桁结合梁桥的ansys分析实例
Ansys课程设计2一、题目以64m双线下承式简支钢桁结合梁桥为研究对象，利用ansys建模分析该结合梁桥在各种工况作用下的变形和应力分布，并进行动力仿真分析，检验其设计的合理性。
该桥每个桥墩（台）在横向布置两个支座，支座中心横向间距为11.0m；分析时仅在桥梁两端设置横联（桥门架），若动力分析通不过，可在中间加设两个横联，截面同桥门架；二期恒载按16t计算，活载按中活载计算，双线加载不折减；横梁先按工形截面计算，若应力通不过，可将端横梁改成箱型断面；纵梁外侧桥面板悬臂部分不与横梁连接。
该桥的结构轮廓图分别如图1、图2、图3、图4所示。
图1 立面图图2 下平面图图3 上平纵联图4 横截面图二、 概述Ansys 建模主要采用自底向上的实体建模，即先定义关键点，然后在关键点的基础上连接成线，之后在初步形成的图形上划分网格，生成节点和单元，并定义截面类型，即可得到有限元模型。
主桁的上下弦杆为梯形截面箱梁，采用自定义截面的beam188空间梁单元进行模拟。
混凝土板采用shell91、横联采用shell63空间板单元，混凝土板的厚度取其平均厚度。
其他杆件采用常用截面的beam188空间梁单元进行模拟。
按平截面假定考虑主桁下弦杆、纵梁、横梁以及混凝土板相互之间的偏心，考虑上平纵联与横联、主桁上弦杆之间的偏心，不考虑钢与混凝土板间的滑移。
主桁各杆之间、横梁与主桁下弦杆之间、纵梁与横梁之间、上平纵联与主桁上弦杆之间、横联与主桁竖杆之间都按刚接处理，混凝土板与纵梁的连接按铰接处理。
单元的平均长度为0.5m ，纯钢结构模型一共得到6880个节点，3949个单元。
混凝土板厚：m 2925.0225.0335.0=+混凝土板宽：m 5.81000%22503352=⨯-⨯ 混凝土板截面面积为2486.22925.05.8m =⨯下弦杆、上弦杆的平均高度取为0.6m钢材的弹性模量取：2.1×1011 Pa ，泊松比取：0.3，密度取：7850㎏/m 3 。
基于ansys的铁路钢桁架桥受力分析
山东农业大学毕业设计题目：基于ansys的福厦铁路钢桁架桥静力分析院部水利土木工程学院专业班级届次学生姓名陈雪峰学号指导教师二O一一年六月十八日目录1工程简介1.1 工程概况1.2 施工方法1.3 加载工况2有限元法基本原理及ansys简介2.1 有限元法基本原理2.2 ansys处理有限元问题的一般方法3福厦铁路桁架桥静力分析3.1 有限元模型3.2 自重荷载工况下的建模及求解3.3 中—活载工况下的建模及求解3.4 计算结果分析4总结参考文献致谢词Contents1 Project description1.1 Project Overview1.2 Construction methods1.3 Load conditions2 The basic principles of finite element method and ansys Profile2.1 The basic principles of finite element method2.2 Ansys finite element problems dealing with the general approach3Fuzhou-Xiamen Railway static analysis of truss bridges3.1 Finite element model3.2 Weight load conditions modeling and solving3.3 The live load conditions of the modeling and solution3.4 Calculation results4 SummaryReferencesThanks Words基于ANSYS的铁路钢桁架桥静力分析作者：陈雪峰道路桥梁与渡河工程3班指导教师：白润波摘要：本毕业设计为“基于ansys的福厦铁路钢桁架桥静力分析”，包括三个主要过程：收集设计资料，确定桥梁的设计方案，材料特性；用已有的资料采用ansys软件运用命令流方法建模分析结构的受力变形；对得出的结果进行分析，并做出正确合合理解释。
安庆长江铁路大桥ANSYS建模命令流
桥址概况安庆长江铁路大桥是南京至安庆城际铁路和阜阳至景德镇铁路的重要组成部分，位于安庆前江口汇合口处下游官山咀附近，距上游已建成通车的安庆长江公路大桥约21km；线路在池州侧晏塘镇靠近长江的刘村附近右拐过江，过江后从安庆的长风镇穿过。
安庆铁路长江大桥全长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。
根据不同的受力区段选用不同的杆件截面，在辅助墩附近的压重区梁段，腹杆采用箱型截面杆件。
ANSYS案例——20例ANSYS经典实例】
ANSYS案例——20例ANSYS经典实例】针对【典型例题】3.3.7(1)的模型，即如图3-19所示的框架结构，其顶端受均布力作用，用有限元方法分析该结构的位移。
结构中各个截面的参数都为：113.010PaE=，746.510mI-=，426.810mA-=，相应的有限元分析模型见图3-20。
在ANSYS平台上，完成相应的力学分析。
图3-19框架结构受一均布力作用图3-20单元划分、节点位移及节点上的外载解答对该问题进行有限元分析的过程如下。
1．基于图形界面的交互式操作(tepbytep)(1)进入ANSYS(设定工作目录和工作文件)程序→ANSYS→ANSYSInteractive→Workingdirectory(设置工作目录)→Initialjobname(设置工作文件名):beam3→Run→OK(2)设置计算类型(3)选择单元类型(4)定义材料参数ANSYSMainMenu:Preproceor→MaterialProp→MaterialModel→Struc tural→Linear→Elatic→Iotropic:E某:3e11(弹性模量)→OK→鼠标点击该窗口右上角的“”来关闭该窗口(5)定义实常数以确定平面问题的厚度ANSYSMainMenu:Preproceor→RealContant…→Add/Edit/Delete→Add→Type1Beam3→OK→RealContantSetNo:1(第1号实常数),Cro-ectionalarea:6.8e-4(梁的横截面积)→OK→Cloe(6)生成几何模型生成节点ANSYSMainMenu:Preproceor→Modeling→Creat→Node→InActiveCS→Nodenumber1→某:0,Y:0.96,Z:0→Apply→Nodenumber2→某:1.44,Y:0.96,Z:0→Apply→Nodenumber3→某:0,Y:0,Z:0→Apply→Nodenumber4→某:1.44,Y:0,Z:0→OK生成单元ANSYSMainMenu:Preproceor→Modeling→Create→Element→AutoNum bered→ThruNode→选择节点1，2(生成单元1)→apply→选择节点1，3(生成单元2)→apply→选择节点2，4(生成单元3)→OK(7)模型施加约束和外载左边加某方向的受力ANSYSMainMenu:Solution→DefineLoad→Apply→Structural→Force/Moment→OnNode→选择节点1→apply→Directionofforce:F某→VALUE：3000→OK→上方施加Y方向的均布载荷ANSYSMainMenu:Solution→DefineLoad→Apply→Structural→Preure→OnBeam→选取单元1(节点1和节点2之间)→apply→VALI：4167→VALJ：4167→OK左、右下角节点加约束(8)分析计算(9)结果显示(10)退出系统(11)计算结果的验证与MATLAB支反力计算结果一致。
ANSYS模态分析实例教程文件
8
9. 对于Y 和Z 方向，输入0替换“Free”， “X” 方向仍 然保留为 “Free”.
Workshop Supplement
9
July 10, 2007 © 2007 ANSYS, Inc. All rights reserved.
ANSYS, Inc. Proprietary
Inventory #002406 WS1-8
ANSYS, Inc. Proprietary
Inventory #002406 WS1-4
Workbench-Simulation Dynamics
Workshop 1 – 设置
1. 当导入几何模型之后，从“Map of Analysis Types”选择“Static Structural”。
• 目的在于计算某一悬索桥模型的振动特性
– 塔科马海峡大桥，又称为“ Galloping Gertie”，因其在1940年出任意 外的倒塌而出名。
• 在该实例中，首先对桥的模型进行分析，然后计算大桥的固有频 率和固有振型。
• 在求解中，我们考虑到了结构的“预应力”，因为自重载荷预先 作用于桥梁上而导致拉索产生了预拉力，从而增加了整体刚度。
Inventory #002406 WS1-3
Workbench-Simulation Dynamics
Workshop 1 – 起始页
• 从“ WorkBench Project Launcher ”点击“ Simulation”。
–
如果进入“Simulation” ，点击 >File>New
•
Inventory #002406 WS1-6
Workbench-Simulation Dynamics
ANSYS经典应用实例
ANSYS经典应用实例首先，我们来看一个经典的结构分析问题：桥梁的静力学分析。
假设我们需要对一座跨度较大的桥梁进行设计和分析。
首先，我们需要建立桥梁的几何模型，并定义桥梁的材料特性和荷载情况。
然后，使用ANSYS软件进行分析。
在ANSYS中，我们可以使用不同的单元类型来建立桥梁的有限元模型。
常见的单元类型有梁单元、板单元和壳单元。
对于桥梁结构来说，我们通常使用梁单元来建模。
首先，我们需要定义桥梁的截面几何形状和材料特性。
ANSYS提供了丰富的材料库，可以选择合适的材料进行模拟。
然后，我们需要划分梁单元的网格，可以使用ANSYS自带的网格划分工具或者手动划分。
在建立了桥梁的有限元模型之后，我们可以给桥梁施加荷载，进行静力学分析。
荷载可以分为静载荷和动载荷。
静载荷包括自重和外部施加的荷载，如车辆荷载、行人荷载等。
动载荷包括地震荷载和风荷载等。
在ANSYS中，我们可以通过定义节点荷载或者面荷载来施加荷载。
同时，我们也可以选择不同的荷载组合进行分析。
完成荷载施加之后，我们可以进行结构的静力学分析。
静力学分析主要是计算结构的位移、应力和应变等关键参数。
ANSYS提供了多种分析方法，如静态分析、模态分析、频率响应分析等。
对于桥梁结构来说，静态分析是最常见的分析方法。
在静态分析中，我们可以得到结构的位移、应力和应变分布，并对其进行验证和优化。
除了静力学分析，ANSYS还可以用于疲劳分析、稳定性分析、动力学分析等。
疲劳分析是用于评估结构在重复加载下的损伤程度和寿命，稳定性分析是用于评估结构的抗侧移性能，动力学分析是用于研究结构在振动荷载下的响应。
这些分析方法可以帮助工程师更好地设计和优化结构。
总结起来，ANSYS在结构分析中的经典应用实例包括桥梁的静力学分析、疲劳分析、稳定性分析和动力学分析等。
通过使用ANSYS软件，工程师可以更加深入地了解结构的性能，并进行有效的设计和优化。
基于ansys及铁路钢桁架桥受力分析
山东农业大学毕业设计题目：基于ansys的福厦铁路钢桁架桥静力分析院部水利土木工程学院专业班级届次学生姓名陈雪峰学号指导教师二O一一年六月十八日目录1工程简介1.1 工程概况1.2 施工方法1.3 加载工况2有限元法基本原理及ansys简介2.1 有限元法基本原理2.2 ansys处理有限元问题的一般方法3福厦铁路桁架桥静力分析3.1 有限元模型3.2 自重荷载工况下的建模及求解3.3 中—活载工况下的建模及求解3.4 计算结果分析4总结参考文献致谢词Contents1 Project description1.1 Project Overview1.2 Construction methods1.3 Load conditions2 The basic principles of finite element method and ansys Profile2.1 The basic principles of finite element method2.2 Ansys finite element problems dealing with the general approach3 Fuzhou-Xiamen Railway static analysis of truss bridges3.1 Finite element model3.2 Weight load conditions modeling and solving3.3 The live load conditions of the modeling and solution3.4 Calculation results4 SummaryReferencesThanks Words基于ANSYS的铁路钢桁架桥静力分析作者：陈雪峰道路桥梁与渡河工程3班指导教师：白润波摘要：本毕业设计为“基于ansys的福厦铁路钢桁架桥静力分析”，包括三个主要过程：收集设计资料，确定桥梁的设计方案，材料特性；用已有的资料采用ansys软件运用命令流方法建模分析结构的受力变形；对得出的结果进行分析，并做出正确合合理解释。
ANSYS3D桥梁
ANSYS 3D橋樑分析一橋樑兩端點承受向下力量(FY=100000N)，如上圖，已知下方承受簡支鋼架(A=0.0187，IZZ=0.00017，IYY=0.00054，TKZ=0.4，TKY=0.4)，每段12m共6段全長72m，高度16m，橋寬10m，橋面板0.3m厚的混凝土。
開始－＞程式集－＞ANSYS5.7－＞Interactive檢查Production Selection及Working Directory、修改Initial Jobname為3DBridge後－＞RUN 進入ANSYS，在Main Menu點選1.Preferences－＞點選Structural－＞OK進入Preprocessor2.Preprocessor－＞Element Type－＞Add/Edit/Delete－＞Add－＞選擇Structural Beam 3Delastic 4－＞Apply－＞選擇Structural Shell Elastic 4node 63－＞OK－＞Close－＞3.Real Constant－＞Add/Edit/Delete－＞選擇Type 1 Beam 4－＞Add－＞OK－＞在AREA框內輸入0.0187、IZZ框內輸入0.00017、IYY框內輸入0.00054、TKZ框內輸入0.4、TKY 框內輸入0.4－＞OK－＞選擇Type 2 Shell 63－＞Add－＞在TK(I)框內輸入0.3－＞OK－＞Close－＞4.Material Props－＞Material Models－＞Structural－＞Linear－＞Elastic－＞Isotropic－＞在EX框內輸入2e11、在PRXY框內輸入0.3、在Density框內輸入7850－＞在下拉式Material －＞New Model－＞確認材料ID為2後OK－＞Structural－＞Linear－＞Elastic－＞Isotropic－＞在EX框內輸入0.35e11、在PRXY框內輸入0.1667、在Density框內輸入2500－＞OK－＞ －＞鋼架5.Modeling Create－＞Nodes－＞In Active CS－＞產生對話框在NPT框內輸入1、在XYZ框內分別輸入0、0、-5－＞OK－＞6.在Modeling Copy－＞Nodes－＞Cop－＞點選框點選Pick all產生對話框在ITIME框內輸入4、在DX框內輸入12、在INC框內輸入4－＞OK－＞點選框點選Pick all產生對話框在ITIME框內輸入2、在DZ框內輸入10、在INC框內輸入1－＞OK－＞點選框點選節點2、6、10－＞OK－＞產生對話框在ITIME框內輸入2、在DY框內輸入16、在INC 框內輸入1－＞OK－＞點選框點選節點3、7、11－＞OK－＞產生對話框在ITIME框內輸入2、在DZ框內輸入-10、在INC框內輸入1－＞OK－＞7.Modeling Create－＞Element－＞Auto Numbered Thru Nodes－＞出現點選框分別點選節點2、6－＞Apply－＞點選節點6和10、節點10和14、節點1和5、節點5和9、節點9和13、節點3和7、節點7和11、節點4和8、節點8和12、節點1和2、節點3和4、節點5和6、節點7和8、節點9和10、節點11和12、節點13和14－＞OK－＞8.Modeling Create－＞Element－＞Auto Numbered Thru Nodes－＞出現點選框分別點選節點11和14－＞Apply－＞點選節點12和13－＞OK－＞9.Modeling Create－＞Element－＞Auto Numbered Thru Nodes－＞出現點選框分別點選節點3、6－＞Apply－＞點選節點6和11、節點4和5、節點5和12、節點2和3、節點1和4、節點6和7、節點5和8、節點10和11、節點9和12－＞OK－＞混凝土10.Modeling Create－＞Elements－＞Element Attributes－＞出現對話框確認混凝土ElementType及Material Number和Real Constant均為2－＞OK11.Modeling Create－＞Element－＞Auto Numbered Thru Nodes－＞出現點選框分別點選節點1、2、6、5－＞Apply－＞點選節點5、6、10、9－＞Apply－＞點選節點9、10、14、13－＞OK－＞全模型並壓縮檔案12.Modeling Reflect－＞Nodes－＞打開Reflect Nodes對話框－＞點選框選Pick all－＞在第二個對話框點選Y-Z Plane－＞在INC框內輸入14－＞OK－＞13.Modeling Reflect－＞Element－＞Auto Numbered打開Reflect Elem對話框－＞點選框選Pick all－＞在第二個對話框點NINC框內輸入14－＞OK－＞14.在Numbering Ctrls－＞將Merge Item打開、Label選ALL－＞OK15.在Numbering Ctrls－＞將Compress Number打開、Label選ALL－＞OK進入Solution16.Solution－＞Loads Apply－＞Structural Displacement－＞On Nodes－＞以點選框點選節點23和24－＞OK－＞出現對話框在Lab2選取UXUYUZ、在V ALUE框內輸入0－＞Apply －＞點選節點13和14－＞OK－＞在對話框Lab2選取UYUZ、在V ALUE框內輸入0－＞OK－＞17.Solution－＞Loads Apply－＞Structural Force/Moment－＞On Nodes－＞出現點選框點選節點1和2－＞OK－＞出現對話框在Lab選FY、在V ALUE輸入-100000－＞OK－＞ －＞18.Solution－＞Loads Apply－＞Structural Inertial－＞Gravity－＞Global－＞出現對話框在ACEL Y輸入9.8－＞OK－＞19.Solve Current LS－＞出現對話框點選OK－＞運算進行中、注意文字視窗敘述－＞在黃色對話框Solution is done點選OK、關閉文字視窗－＞進入General PostProc20.General PostProc－＞List Results－＞Nodal Solution－＞出現對話框在DOF Solution AllDOFs DOF－＞OK－＞產生文字視窗檢查各節點之位移－＞ －＞21.General PostProc－＞List Results－＞Reaction Solu－＞出現對話框在Lab選All Items－＞OK－＞產生文字視窗檢查各節點之反作用力－＞ －＞22.在下拉式Utility Menu選取PlotCtrls－＞Animate－＞Deformed Shape－＞OK－＞Close23.在下拉式Utility Menu選取PlotCtrls－＞Hard Copy－＞To File－＞在選取框中點選JPEG、Reversed Video、存檔成3DBridge.jpg－＞OK如背景有顏色先至PlotCtrls－＞Style－＞Background－＞取消Display Picture Background 24.在下拉式Utility Menu選取File－＞Save as Jobname.db。
隧道及地下工程ANSYS实例分析-第5章-双线铁路隧道衬砌结构设计力学分析
5.1二次衬砌结构力学分析/TITLE，Mechanical analysis on railway tunnel 2nd lining ! 确定分析标题/NOPR !菜单过滤设置/PMETH，OFF，0KEYW，PR_SET，1KEYW，PR_STRUC，1 !保留结构分析部分菜单/COM，/COM，Preferences for GUI filtering have been set to display:/COM，Structural!/PREP7 !进入前处理器ET，1，BEAM3 !设置梁单元类型ET，2，COMBIN14 !设置弹簧单元类型R，1，0.4，0.0053333，0.4，，，，!设置梁单元几何常数R，2，400e6，，，!设置弹簧单元几何常数MPTEMP，，，，，，，，!设置材料模型MPTEMP，1，0MPDATA，EX，1，，30.0e9 !输入弹性模量MPDATA，PRXY，1，，0.2 !输入泊松比MPTEMP，，，，，，，，!设置材料模型MPTEMP，1，0MPDATA，DENS，1，，2500 !输入密度SAVE !保存数据库1．建立几何模型K, 1, , , , !创建关键点（隧道二次衬砌）K, 2, 4.71, -1.82, ,K, 3, -4.71, -1.82, ,K, 4, 0, -3.75, ,K, 5, 0, 5.05, ,K, 20, 5.6, -2.162, , !创建关键点（地层弹簧）K, 30, -5.6, -2.16, ,K, 40, 0, -4.75, ,K, 50, 0, 6, ,!创建隧道衬砌线LARC, 2, 5, 1, 5.05, !创建圆弧线（拱顶部）LARC, 5, 3, 1, 5.05,LARC, 2, 3, 4 !创建圆弧线（仰拱部）第1章大型有限元软件ANSYS简介2 ! 创建地层弹簧线LARC, 20, 50, 1, 6, !创建圆弧线（拱顶部）LARC,50, 30, 1, 6,LARC, 20, 30, 40 !创建圆弧线（仰拱部）SAVE !保存数据2．单元网格划分设置单元大小并将所有直线划分单元，其单元图如图5-8所示。
基于Ansys的某平列式公铁两用大桥力学特性分析
甘肃科技第２５卷为钢桁桥和第５跨的３２ｍ预应力钢筋混凝土简支梁桥共用，故将预应力钢筋混凝土简支梁的一半质量加在４号墩墩顶。
全桥模型顺桥向为ｚ方向（正方向为小里程侧），横向为ｘ方向，竖向为Ｙ方向，具体如图２、图３所示。
围２平列式公铁两用大桥全轿空同模型图３平列式公铁两用大桥上部结构模型（局部）３由静力特性探析公路部分对大桥的影响结构的变形特点反映了结构受力的综合情祝，掌握桥梁受外荷载作用后各个节点的位移变化情况，是了解该类桥梁结构工作状态的一个重要参数。
对平列式公铁两用大桥的静力分析主要是围绕结构的变形展开的。
假如拆除伸出的钢悬臂以及钢筋混凝土简支梁，很明显，该桥仅为连续下承式钢桁梁，属平行弦的连续桁梁，故可把它进行简化，近似地看成具有等截面惯矩的不等跨连续梁。
经分析不难发现，无论中一荷载从哪个方向进入，都会使连续梁第四跨承受最不利的荷载，为更加全面地了解桥梁的力学特性，将第三跨也纳入研究，据此拟定大桥第三跨和第四跨上相应节点挠度为分析对象。
依据设计荷载，在进行静力分析时分为３种荷载工况：（１）第一工况，中＿２２级荷载单独加载；（２）第二工况，汽１３级荷载单独加载；（３）第三工况，汽车荷载和列车荷载共同加载。
在模拟列车荷载的作用时，考虑轮轴荷载的大小和位置等，根据荷载分布情况最终作用在纵梁上翼缘上，汽车荷载的模拟是以钢筋混凝土简支梁的支反力反作用力的形式作用在钢悬臂上，这样能比较准确地模拟荷载作用，以上荷载均加于最不利位置处。
运用ＡＮＳＹＳ对所建立的模型进行结构挠度计算，可得到荷载作用下各个单元的具体挠度值，对计算结果进行整理，得到全桥主桁架在主要控制点不同工况下的挠度值如表ｌ所示。
表１各工况下位移值单位．ｍｍ注：４—１／２上一指从大里程侧算起第４跨钢桁梁上游侧Ｌ／２处；位移值为负表示变形向下。
通过表１可以很清楚地看出，在工况一荷载作构的扭转变形也更大，上下游两片主桁在相对应节用下，最大竖向向下位移发生在第四跨跨中节点处，点位移同第二工况不同，竖向位移值的符号都为同最大为一１４．７１ｒａｍ，上下游两片主桁相对应节点位正或者同负。
桥梁结构ANSYS建模原则及常见问题
4 基于 APDL 的编程及 ANSYS 新模块的实现 ......................... 34
4.1 APDL 的特点 ...................................................... 34 4.2 APDL 的编程 ...................................................... 34 4.2.1 APDL 的参数化变量及函数 ................................... 34 4.2.3 宏命令 ................................................... 37 4.2.4 APDL 常用命令详解 ......................................... 38 4.3 ANSYS 新功能 ..................................................... 39
5 实用技巧 .................................................... 43
ansys在公路工程中的应用实例
6.1施工过程仿真分析1．材料、实常数和单元类型定义/TITLE，Mechanicalanalysisonrailwaytunnel1ndlining !确定分析标题/NOPR !菜单过滤设置/PMETH，OFF，0KEYW，PR_SET，1KEYW，PR_STRUC，1 !保留结构分析部分菜单/COM，/COM，PreferencesforGUIfilteringhavebeensettodisplay:/COM，Structural/PREP7 !进入前处理器ET,1,PLANE42 !设置实体单元类型KEYOPT,1,1,0KEYOPT,1,2,0KEYOPT,1,3,2 !设置为平面应变模式KEYOPT,1,5,0KEYOPT,1,6,0ET,2,PLANE42 !设置实体单元类型KEYOPT,1,1,0KEYOPT,1,2,0KEYOPT,1,3,2 !设置为平面应变模式KEYOPT,1,5,0KEYOPT,1,6,0ET，3，BEAM3 !设置梁单元类型R,1,0.25,1.5e-3,0.25,,,, !设置梁单元几何常数MPTEMP,,,,,,,,MPTEMP,1,0MPDATA,EX,1,,0.8e9 !输入弹性模量（围岩）MPDATA,PRXY,1,,0.38 !输入泊松比（围岩）MPDATA,DENS,1,,2000 !输入密度（围岩）TB,DP,1,,, !采用DP准则进行弹塑性分析TBMODIF,1,1,0.2e6 !输入凝聚力（围岩）TBMODIF,1,2,30 !输入摩擦角（围岩）TBMODIF,1,3,MPTEMP,,,,,,,,MPTEMP,1,0MPDATA,EX,2,,3.2e9 !输入弹性模量（加固圈）MPDATA,PRXY,2,,0.32 !输入泊松比（加固圈）MPDATA,DENS,2,,2200 !输入密度（加固圈）TB,DP,2,,, !采用DP准则TBMODIF,1,1,0.5e6 !输入凝聚力（加固圈）TBMODIF,1,2,35 !输入摩擦角（加固圈）TBMODIF,1,3,MPDATA,EX,3,,27.5e9 !输入弹性模量（喷射混凝土）MPDATA,PRXY,3,,0.2 !输入泊松比（喷射混凝土）MPDATA,DENS,3,,2500 !输入密度（喷射混凝土）SAVE !保存数据库2．建立几何模型创建关键点。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
桥址概况安庆长江铁路大桥是南京至安庆城际铁路和阜阳至景德镇铁路的重要组成部分，位于安庆前江口汇合口处下游官山咀附近，距上游已建成通车的安庆长江公路大桥约21km；线路在池州侧晏塘镇靠近长江的刘村附近右拐过江，过江后从安庆的长风镇穿过。
安庆铁路长江大桥全长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高强螺栓。
项目进展2005年元月，安庆长江铁路大桥项目前期工作协调领导小组办公室委托铁道第四勘察设计院编制《安徽省铁路总体规划安庆过江通道深化研究及大桥选址报告》，随后铁四院专家组来宜现场勘察，采集相关资料，并于2月份完成该报告。
2005年8月，安徽省发改委主持召开“安庆长江铁路大桥桥位专家咨询会”，邀请中国工程院院士陈新等知名专家对大桥桥位进行咨询研究，并对选址报告进行评审。
2005年10月，经省部协商，铁道部将沿江城际铁路及安庆长江铁路大桥项目补列入国家“十一五”规划，并向铁四院下达前期工作任务书。
2005年12月，铁道部主持召开了宁安城际铁路及安庆长江铁路大桥项目预可研报告审查会。
2006年元月，宁安城际铁路及安庆长江铁路大桥项目列入国家“十一五”规划。
2006年7月5日，安庆长江铁路大桥设计竞标工作会议在北京举行，经研究、审核、竞标，铁四院与中铁大桥院联合体中标。
2006年9月，中铁大桥院和铁四院正式启动工程可行性研究报告的编制程序,同时编制了桥梁建设对长江航道的影响书,河势分析,桥位所在枯、中水流影响、流速及航迹图。
2006年10月9日，安庆市与铁四院、中铁大桥院联合举行了宁安城际铁路及安庆长江铁路大桥工可勘察开钻仪式;10月份武汉水文站和中铁大桥院近100人的设计、水文、勘测、物探等四个工作组进入安庆进行实地勘测。
2009年3月20日，安庆长江铁路大桥开工建设表1.2 安庆长江大桥主要工程数量表表 2.1 设计构件节点和单元数表 2.2安庆长江铁路大桥模型主要参数列表表2.3 模型实常数命令流文件/nerr,0!将参数读入数据库finish/clear/PREP7!单元类型编号et,1,beam44et,2,shell63et,3,link8!材料类型编号 Q420/Q370/拉索-1，混凝土C40/C50-2，刚臂-3 mp,ex,1,2.1e11mp,prxy,1,0.3mp,dens,1,7850mp,alpx,1,1e-5mp,ex,2,3.5e7mp,prxy,2,0.16mp,dens,2,2500mp,alpx,2,1e-5mp,ex,3,1e15mp,prxy,3,1mp,dens,3,0mp,alpx,3,0!读取自定义截面并编号*dim,name,char,24name(1)='Sec_1a'name(2)='Sec_2a'name(3)='Sec_3a'name(4)='Sec_4a'name(5)='Sec_5a'name(6)='Sec_6a'name(7)='Sec_7a'name(8)='Sec_1'name(9)='Sec_2'name(10)='Sec_3'name(11)='Sec_4'name(12)='Sec_5'name(13)='Sec_6'name(14)='Sec_7'name(15)='Sec_8'name(16)='Sec_9'name(17)='Sec_10'name(18)='Sec_11'name(19)='Sec_12'name(20)='Sec_13'name(21)='Sec_14'name(22)='Sec_15'name(23)='Sec_16'name(24)='Sec_x2'*do,i,1,24,1sectype,i,beam,mesh,secread,name(i),,,mesh*enddo!将ta_sec参数读入数据库!读取自定义截面并编号*dim,name2,char,15name2(1)='ta_1'name2(2)='ta_2'name2(3)='ta_3'name2(4)='ta_4'name2(5)='ta_5'name2(6)='ta_6'name2(7)='ta_7'name2(8)='ta_8'name2(9)='ta_9'name2(10)='ta_10'name2(11)='ta_11'name2(12)='ta_12'name2(13)='ta_13'name2(14)='ta_14'*do,i,1,14,1sectype,40+i,beam,mesh,secread,name2(i),,,mesh*enddoSECTYPE, 26, BEAM, HREC, , 0SECOFFSET, USER, 0.1, 0.266SECDATA,0.2,0.266,0.006,0.006,0.006,0.006,0,0,0,0 SECTYPE, 27, BEAM, I, , 0SECOFFSET, USER, 0, 0.584SECDATA,0.240,0.240,0.6,0.016,0.016,0.012,0,0,0,0SECTYPE, 28, BEAM, T, T_nin, 0SECOFFSET, USER, 0, 1.4SECDATA,0.6,1.4,0.032,0.018,0,0,0,0,0,0 SECTYPE, 29, BEAM, T, T_max, 0SECOFFSET, USER, 0, 2.1SECDATA,0.6,2.1,0.032,0.016,0,0,0,0,0,0 SECTYPE, 31, BEAM, HREC, , 0SECOFFSET, CENTSECDATA,0.548,0.500,0.016,0.016,0.016,0.016,0,0,0,0 SECTYPE, 32, BEAM, I, , 0SECOFFSET, CENTSECDATA,0.560,0.560,0.500,0.028,0.028,0.016,0,0,0,0 !刚臂的截面SECTYPE, 100, BEAM, CSOLID, , 0SECOFFSET, CENTSECDATA,0.06,0,0,0,0,0,0,0,0,0r,1,r,2,0.016,!拉索截面控制!7-349r,3,134.31e-4,!7-301r,4,115.84e-4,!7-283r,5,108.91e-4,!7-241r,6,100.000e-4,finish!开始建立模型/prep7!*********************************************!!!xiaXianGan!!!**********************************************SET,j,-14*do,i,1,(94*25+1)*2,94*25+1 n,i,0,j,0n,94*25+i,1363.000,j,0 FILL,i,94*25+i,94*25-1,,1,,, !中桁*if,i,eq,94*25+2,then!7截面type,1mat,1secnum,14*do,a,1,3*25,1E,i-1+a,i+a*enddo*do,b,1,11*25,1E,35*25-1+i+b,35*25+i+b *enddo*do,c,1,8*25,1E,77*25-1+i+c,77*25+i+c *enddo!6截面type,1mat,1secnum,13*do,a,1,3*25,1E,i+3*25-1+a,i+3*25+a*enddo*do,b,1,3*25,1E,32*25-1+i+b,32*25+i+b *enddo*do,c,1,3*25,1E,46*25-1+i+c,46*25+i+c *enddo*do,d,1,1*25,1E,i+76*25-1+d,i+76*25+d*enddo*do,e,1,7*25,1E,87*25-1+i+e,87*25+i+e *enddo!5截面type,1mat,1secnum,12*do,a,1,25,1E,8*25-1+i+a,8*25+i+a*do,b,1,2*25,1E,30*25-1+i+b,30*25+i+b *enddo*do,c,1,2*25,1E,49*25-1+i+c,49*25+i+c*enddo*do,d,1,2*25,1E,72*25-1+i+d,72*25+i+d*enddo!4截面type,1mat,1secnum,11*do,a,1,3*25,1E,9*25-1+i+a,9*25+i+a*enddo*do,b,1,4*25,1E,26*25-1+i+b,26*25+i+b *enddo*do,c,1,3*25,1E,51*25-1+i+c,51*25+i+c*enddo*do,d,1,4*25,1E,68*25-1+i+d,68*25+i+d*enddo!3截面type,1mat,1secnum,10*do,a,1,1*25,1E,7*25-1+i+a,7*25+i+a*enddo*do,b,1,5*25,1E,12*25-1+i+b,12*25+i+b *enddo*do,c,1,4*25,1E,22*25-1+i+c,22*25+i+c*enddo*do,d,1,4*25,1E,54*25-1+i+d,54*25+i+d*enddo*do,e,1,5*25,1E,63*25-1+i+e,63*25+i+e!2截面type,1mat,1secnum,9*do,a,1,1*25,1E,6*25-1+i+a,6*25+i+a*enddo*do,b,1,2*25,1E,17*25-1+i+b,17*25+i+b*enddo*do,c,1,1*25,1E,21*25-1+i+c,21*25+i+c*enddo*do,d,1,1*25,1E,58*25-1+i+d,58*25+i+d*enddo*do,e,1,2*25,1E,61*25-1+i+e,61*25+i+e*enddo*do,f,1,2*25,1E,74*25-1+i+f,74*25+i+f*enddo*do,g,1,1*25,1E,85*25-1+i+g,85*25+i+g*enddo!1截面type,1mat,1secnum,8*do,a,1,2*25,1E,19*25-1+i+a,19*25+i+a*enddo*do,b,1,2*25,1E,59*25-1+i+b,59*25+i+b*enddo*do,c,1,1*25,1E,86*25-1+i+c,86*25+i+c*enddo!左边桁*elseif,i,eq,(94*25+1)*2+1,then !7截面type,1secnum,7*do,a,1,3*25,1E,i+a,i-1+a*enddo*do,b,1,11*25,1E,35*25+i+b,35*25-1+i+b *enddo*do,c,1,8*25,1E,77*25+i+c,77*25-1+i+c *enddo!6截面type,1mat,1secnum,6*do,a,1,3*25,1E,i+3*25+a,i+3*25-1+a*enddo*do,b,1,3*25,1E,32*25+i+b,32*25-1+i+b *enddo*do,c,1,3*25,1E,46*25+i+c,46*25-1+i+c *enddo*do,d,1,1*25,1E,i+76*25+d,i+76*25-1+d *enddo*do,e,1,7*25,1E,87*25+i+e,87*25-1+i+e *enddo!5截面type,1mat,1secnum,5*do,a,1,25,1E,8*25+i+a,8*25-1+i+a*enddo*do,b,1,2*25,1E,30*25+i+b,30*25-1+i+b *enddo*do,c,1,2*25,1E,49*25+i+c,49*25-1+i+c*enddo*do,d,1,2*25,1E,72*25+i+d,72*25-1+i+d*enddo!4截面type,1mat,1secnum,4*do,a,1,3*25,1E,9*25+i+a,9*25-1+i+a*enddo*do,b,1,4*25,1E,26*25+i+b,26*25-1+i+b *enddo*do,c,1,3*25,1E,51*25+i+c,51*25-1+i+c*enddo*do,d,1,4*25,1E,68*25+i+d,68*25-1+i+d*enddo!3截面type,1mat,1secnum,3*do,a,1,1*25,1E,7*25+i+a,7*25-1+i+a*enddo*do,b,1,5*25,1E,12*25+i+b,12*25-1+i+b *enddo*do,c,1,4*25,1E,22*25+i+c,22*25-1+i+c*enddo*do,d,1,4*25,1E,54*25+i+d,54*25-1+i+d*enddo*do,e,1,5*25,1E,63*25+i+e,63*25-1+i+e*enddo!2截面type,1mat,1secnum,2*do,a,1,1*25,1E,6*25+i+a,6*25-1+i+a*enddoE,17*25+i+b,17*25-1+i+b *enddo*do,c,1,1*25,1E,21*25+i+c,21*25-1+i+c*enddo*do,d,1,1*25,1E,58*25+i+d,58*25-1+i+d*enddo*do,e,1,2*25,1E,61*25+i+e,61*25-1+i+e*enddo*do,f,1,2*25,1E,74*25+i+f,74*25-1+i+f*enddo*do,g,1,1*25,1E,85*25+i+g,85*25-1+i+g*enddo!1截面type,1mat,1secnum,1*do,a,1,2*25,1E,19*25+i+a,19*25-1+i+a *enddo*do,b,1,2*25,1E,59*25+i+b,59*25-1+i+b *enddo*do,c,1,1*25,1E,86*25+i+c,86*25-1+i+c*enddo!右边桁*else!7截面type,1mat,1secnum,7*do,a,1,3*25,1E,i-1+a,i+a*enddo*do,b,1,11*25,1E,35*25-1+i+b,35*25+i+b *enddoE,77*25-1+i+c,77*25+i+c *enddo!6截面type,1mat,1secnum,6*do,a,1,3*25,1E,i+3*25-1+a,i+3*25+a*enddo*do,b,1,3*25,1E,32*25-1+i+b,32*25+i+b *enddo*do,c,1,3*25,1E,46*25-1+i+c,46*25+i+c *enddo*do,d,1,1*25,1E,i+76*25-1+d,i+76*25+d *enddo*do,e,1,7*25,1E,87*25-1+i+e,87*25+i+e *enddo!5截面type,1mat,1secnum,5*do,a,1,25,1E,8*25-1+i+a,8*25+i+a*enddo*do,b,1,2*25,1E,30*25-1+i+b,30*25+i+b *enddo*do,c,1,2*25,1E,49*25-1+i+c,49*25+i+c*enddo*do,d,1,2*25,1E,72*25-1+i+d,72*25+i+d*enddo!4截面type,1mat,1secnum,4*do,a,1,3*25,1E,9*25-1+i+a,9*25+i+a*do,b,1,4*25,1E,26*25-1+i+b,26*25+i+b *enddo*do,c,1,3*25,1E,51*25-1+i+c,51*25+i+c*enddo*do,d,1,4*25,1E,68*25-1+i+d,68*25+i+d*enddo!3截面type,1mat,1secnum,3*do,a,1,1*25,1E,7*25-1+i+a,7*25+i+a*enddo*do,b,1,5*25,1E,12*25-1+i+b,12*25+i+b *enddo*do,c,1,4*25,1E,22*25-1+i+c,22*25+i+c*enddo*do,d,1,4*25,1E,54*25-1+i+d,54*25+i+d*enddo*do,e,1,5*25,1E,63*25-1+i+e,63*25+i+e*enddo!2截面type,1mat,1secnum,2*do,a,1,1*25,1E,6*25-1+i+a,6*25+i+a*enddo*do,b,1,2*25,1E,17*25-1+i+b,17*25+i+b *enddo*do,c,1,1*25,1E,21*25-1+i+c,21*25+i+c*enddo*do,d,1,1*25,1E,58*25-1+i+d,58*25+i+d*do,e,1,2*25,1E,61*25-1+i+e,61*25+i+e*enddo*do,f,1,2*25,1E,74*25-1+i+f,74*25+i+f*enddo*do,g,1,1*25,1E,85*25-1+i+g,85*25+i+g*enddo!1截面type,1mat,1secnum,1*do,a,1,2*25,1E,19*25-1+i+a,19*25+i+a*enddo*do,b,1,2*25,1E,59*25-1+i+b,59*25+i+b*enddo*do,c,1,1*25,1E,86*25-1+i+c,86*25+i+c*enddo*endif*SET,j,j+14*enddo!*********************************************!!!T_beam AND shell63's nodes!!!*********************************************/prep7ngen,2,7053,1,2351,1,0,0,0.7,ngen,4,2351,7054,7054+2350,1,0,0.57,0,ngen,4,2351,7054+3*2351,7054+3*2351+2350,1,0,0.60,0, ngen,3,2351,7054+6*2351,7054+6*2351+2350,1,0,0.45,0, ngen,2,2351,7054+8*2351,7054+8*2351+2350,1,0,0.60,0,ngen,3,2351,7054+9*2351,7054+9*2351+2350,1,0,0.45,0, ngen,2,2351,7054+11*2351,7054+11*2351+2350,1,0,0.70,0,local,11,0,0,-7,0nsym,y,2351*12,7054+2351,32914+2*2351,1csys,0nsym,y,58775,7054,65828,1ngen,2,122252,2352,2352+2350,1,0,0,0.7,!生成边上 4条 U肋type,1mat,1secnum,26*do,h,6,9,1*do,k,1,2350,1e,2351*h+k,2351*h+1+k*enddo*enddo!生成T梁1type,1mat,1secnum,27*do,k,1,2350,1e,2351*10+k,2351*10+1+k*enddo!生成 U肋 2条type,1mat,1secnum,26*do,h,11,12,1*do,k,1,2350,1e,2351*h+k,2351*h+1+k*enddo*enddo!生成T梁2type,1secnum,27*do,k,1,2350,1e,2351*13+k,2351*13+1+k*enddo!生成 U肋另外2条type,1mat,1secnum,26*do,h,14,15,1*do,k,1,2350,1e,2351*h+k,2351*h+1+k*enddo*enddo!镜像四分之一纵向的单元csys,11!（最好用*get 最大节点号单元号）esym,,2351*12,7051,30550,1!生成横梁type,1mat,1secnum,29*do,h,7054,9404,25*do,k,h,h+11*2351,2351e,k,k+2351*enddo*enddo*do,h,37617,39967,25*do,k,h,h+10*2351,2351e,k,k+2351*enddo*enddo*do,h,35266,37616,25e,h,h+2351*12*enddo!生成横肋mat,1secnum,28*do,ah,7059,9384,25*do,ac,0,15,5*do,ak,ah+ac,ah+ac+11*2351,2351 e,ak,ak+2351*enddo*enddo*enddo*do,ah,37622,39947,25*do,ac,0,15,5*do,ak,ah+ac,ah+ac+10*2351,2351 e,ak,ak+2351*enddo*enddo*enddo*do,ah,35271,37596,25*do,ac,0,15,5e,ah+ac,ah+ac+2351*12*enddo*enddo!镜像一半单元csys,0esym,,58775,7051,65354,1!补充横梁type,1mat,1secnum,29*do,h,37617,39967,25e,h,h+86987*enddo*do,h,96392,98742,25e,h,h+28212*enddo!补充横肋type,1mat,1secnum,28*do,ah,37622,39947,25*do,ac,0,15,5e,ah+ac,ah+ac+86987*enddo*enddo*do,ah,96397,98722,25*do,ac,0,15,5e,ah+ac,ah+ac+28212*enddo*enddo!刚臂type,1mat,3secnum,100*do,h,1,2351,1e,h,h+7053*enddo*do,h,4703,7053,1e,h,h+61126*enddo*do,h,2352,4702,1e,h,h+124604-2352*enddo!*********************************************!!!shell!!!********************************************* type,2mat,1real,2*do,i,0,11,1*do,h,7054,7054+2349,1e,h+i*2351,h+1+i*2351,h+2352+i*2351,h+2351+i*2351 *enddo*enddo*do,h,35266,35266+2349,1e,h,h+1,h+63479-35266,h+63478-35266*enddo*do,i,0,-10,-1*do,h,63478,63478+2349,1e,h+i*2351,h+1+i*2351,h+i*2351-2350,h+i*2351-2351 *enddo*enddo*do,h,37617,37617+2349,1e,h,h+1,h+1+124604-37617,h+124604-37617*enddocsys,0esym,,58775,131654,190403,1*do,h,124604,124604+2349,1e,h,h+1,h+1-(124604-96392),h-(124604-96392)*enddo!*********************************************!!!ShangXianGan!!!*********************************************/prep7real,1*SET,Nmax1,300000*SET,j,-14*do,i,Nmax1+1,Nmax1+(92*25+1)*2,92*25+1n,i,14.50,j,15.00n,92*25+i,1348.50,j,15.00FILL,i,92*25+i,,,,,,!中桁*if,i,eq,Nmax1+92*25+2,then!13截面type,1mat,1secnum,20*do,a,1,2*25,1E,i-1+a,i+a*do,b,1,6*25,1E,3*25-1+i+b,3*25+i+b *enddo*do,c,1,3*25,1E,70*25-1+i+c,70*25+i+c *enddo*do,d,1,10*25,1E,74*25-1+i+d,74*25+i+d *enddo*do,e,1,3*25,1E,85*25-1+i+e,85*25+i+e *enddo*do,f,1,1*25,1E,91*25-1+i+f,91*25+i+f *enddo!12截面type,1mat,1secnum,19*do,a,1,1*25,1E,i+2*25-1+a,i+2*25+a*enddo*do,b,1,1*25,1E,9*25-1+i+b,9*25+i+b *enddo*do,c,1,8*25,1E,35*25-1+i+c,35*25+i+c *enddo*do,d,1,1*25,1E,i+69*25-1+d,i+69*25+d *enddo*do,e,1,1*25,1E,73*25-1+i+e,73*25+i+e *enddo*do,f,1,3*25,1E,88*25-1+i+f,88*25+i+f *enddo!11截面type,1mat,1secnum,18*do,a,1,25,1E,10*25-1+i+a,10*25+i+a*do,b,1,1*25,1E,18*25-1+i+b,18*25+i+b *enddo*do,c,1,1*25,1E,26*25-1+i+c,26*25+i+c*enddo*do,d,1,4*25,1E,31*25-1+i+d,31*25+i+d*enddo*do,e,1,1*25,1E,43*25-1+i+e,43*25+i+e*enddo*do,f,1,5*25,1E,45*25-1+i+f,45*25+i+f*enddo*do,g,1,1*25,1E,59*25-1+i+g,59*25+i+g*enddo*do,h,1,1*25,1E,68*25-1+i+h,68*25+i+h*enddo!10截面type,1mat,1secnum,17*do,a,1,1*25,1E,11*25-1+i+a,11*25+i+a *enddo*do,b,1,1*25,1E,17*25-1+i+b,17*25+i+b *enddo*do,c,1,1*25,1E,20*25-1+i+c,20*25+i+c*enddo*do,d,1,1*25,1E,25*25-1+i+d,25*25+i+d*enddo*do,e,1,4*25,1E,27*25-1+i+e,27*25+i+e*enddo*do,f,1,3*25,1E,50*25-1+i+f,50*25+i+f*enddo*do,g,1,1*25,1E,57*25-1+i+g,57*25+i+g*enddo*do,h,1,1*25,1E,60*25-1+i+h,60*25+i+h*enddo*do,hh,1,1*25,1E,67*25-1+i+hh,67*25+i+hh*enddo!9截面type,1mat,1secnum,16*do,a,1,5*25,1E,12*25-1+i+a,12*25+i+a*enddo*do,b,1,1*25,1E,19*25-1+i+b,19*25+i+b*enddo*do,c,1,4*25,1E,21*25-1+i+c,21*25+i+c*enddo*do,d,1,4*25,1E,53*25-1+i+d,53*25+i+d*enddo*do,e,1,1*25,1E,58*25-1+i+e,58*25+i+e*enddo*do,g,1,6*25,1E,61*25-1+i+g,61*25+i+g*enddo*do,h,1,1*25,1E,84*25-1+i+h,84*25+i+h*enddo!8截面type,1mat,1secnum,15*do,a,1,1*25,1E,44*25-1+i+a,44*25+i+a*enddo!左边桁*elseif,i,eq,Nmax1+(92*25+1)*2+1,thentype,1mat,1secnum,20*do,a,1,2*25,1E,i-1+a,i+a*enddo*do,b,1,6*25,1E,3*25-1+i+b,3*25+i+b *enddo*do,c,1,3*25,1E,70*25-1+i+c,70*25+i+c *enddo*do,d,1,10*25,1E,74*25-1+i+d,74*25+i+d *enddo*do,e,1,3*25,1E,85*25-1+i+e,85*25+i+e *enddo*do,f,1,1*25,1E,91*25-1+i+f,91*25+i+f *enddo!12截面type,1mat,1secnum,19*do,a,1,1*25,1E,i+2*25-1+a,i+2*25+a*enddo*do,b,1,1*25,1E,9*25-1+i+b,9*25+i+b *enddo*do,c,1,8*25,1E,35*25-1+i+c,35*25+i+c *enddo*do,d,1,1*25,1E,i+69*25-1+d,i+69*25+d *enddo*do,e,1,1*25,1E,73*25-1+i+e,73*25+i+e *enddo*do,f,1,3*25,1E,88*25-1+i+f,88*25+i+f *enddotype,1mat,1secnum,18*do,a,1,25,1E,10*25-1+i+a,10*25+i+a *enddo*do,b,1,1*25,1E,18*25-1+i+b,18*25+i+b *enddo*do,c,1,1*25,1E,26*25-1+i+c,26*25+i+c*enddo*do,d,1,4*25,1E,31*25-1+i+d,31*25+i+d*enddo*do,e,1,1*25,1E,43*25-1+i+e,43*25+i+e*enddo*do,f,1,5*25,1E,45*25-1+i+f,45*25+i+f*enddo*do,g,1,1*25,1E,59*25-1+i+g,59*25+i+g*enddo*do,h,1,1*25,1E,68*25-1+i+h,68*25+i+h*enddo!10截面type,1mat,1secnum,17*do,a,1,1*25,1E,11*25-1+i+a,11*25+i+a *enddo*do,b,1,1*25,1E,17*25-1+i+b,17*25+i+b *enddo*do,c,1,1*25,1E,20*25-1+i+c,20*25+i+c*enddo*do,d,1,1*25,1E,25*25-1+i+d,25*25+i+d*enddoE,27*25-1+i+e,27*25+i+e*enddo*do,f,1,3*25,1E,50*25-1+i+f,50*25+i+f*enddo*do,g,1,1*25,1E,57*25-1+i+g,57*25+i+g*enddo*do,h,1,1*25,1E,60*25-1+i+h,60*25+i+h*enddo*do,hh,1,1*25,1E,67*25-1+i+hh,67*25+i+hh *enddo!9截面type,1mat,1secnum,16*do,a,1,5*25,1E,12*25-1+i+a,12*25+i+a *enddo*do,b,1,1*25,1E,19*25-1+i+b,19*25+i+b *enddo*do,c,1,4*25,1E,21*25-1+i+c,21*25+i+c*enddo*do,d,1,4*25,1E,53*25-1+i+d,53*25+i+d*enddo*do,e,1,1*25,1E,58*25-1+i+e,58*25+i+e*enddo*do,g,1,6*25,1E,61*25-1+i+g,61*25+i+g*enddo*do,h,1,1*25,1E,84*25-1+i+h,84*25+i+h*enddo!8截面type,1mat,1secnum,15E,44*25-1+i+a,44*25+i+a *enddo!右边桁*else!13截面type,1mat,1secnum,20*do,a,1,2*25,1E,i-1+a,i+a*enddo*do,b,1,6*25,1E,3*25-1+i+b,3*25+i+b *enddo*do,c,1,3*25,1E,70*25-1+i+c,70*25+i+c *enddo*do,d,1,10*25,1E,74*25-1+i+d,74*25+i+d *enddo*do,e,1,3*25,1E,85*25-1+i+e,85*25+i+e *enddo*do,f,1,1*25,1E,91*25-1+i+f,91*25+i+f *enddo!12截面type,1mat,1secnum,19*do,a,1,1*25,1E,i+2*25-1+a,i+2*25+a*enddo*do,b,1,1*25,1E,9*25-1+i+b,9*25+i+b *enddo*do,c,1,8*25,1E,35*25-1+i+c,35*25+i+c *enddo*do,d,1,1*25,1E,i+69*25-1+d,i+69*25+d *enddoE,73*25-1+i+e,73*25+i+e *enddo*do,f,1,3*25,1E,88*25-1+i+f,88*25+i+f *enddo!11截面type,1mat,1secnum,18*do,a,1,25,1E,10*25-1+i+a,10*25+i+a *enddo*do,b,1,1*25,1E,18*25-1+i+b,18*25+i+b *enddo*do,c,1,1*25,1E,26*25-1+i+c,26*25+i+c*enddo*do,d,1,4*25,1E,31*25-1+i+d,31*25+i+d*enddo*do,e,1,1*25,1E,43*25-1+i+e,43*25+i+e*enddo*do,f,1,5*25,1E,45*25-1+i+f,45*25+i+f*enddo*do,g,1,1*25,1E,59*25-1+i+g,59*25+i+g*enddo*do,h,1,1*25,1E,68*25-1+i+h,68*25+i+h*enddo!10截面type,1mat,1secnum,17*do,a,1,1*25,1E,11*25-1+i+a,11*25+i+a *enddo*do,b,1,1*25,1E,17*25-1+i+b,17*25+i+b *enddoE,20*25-1+i+c,20*25+i+c*enddo*do,d,1,1*25,1E,25*25-1+i+d,25*25+i+d*enddo*do,e,1,4*25,1E,27*25-1+i+e,27*25+i+e*enddo*do,f,1,3*25,1E,50*25-1+i+f,50*25+i+f*enddo*do,g,1,1*25,1E,57*25-1+i+g,57*25+i+g*enddo*do,h,1,1*25,1E,60*25-1+i+h,60*25+i+h*enddo*do,hh,1,1*25,1E,67*25-1+i+hh,67*25+i+hh *enddo!9截面type,1mat,1secnum,16*do,a,1,5*25,1E,12*25-1+i+a,12*25+i+a *enddo*do,b,1,1*25,1E,19*25-1+i+b,19*25+i+b *enddo*do,c,1,4*25,1E,21*25-1+i+c,21*25+i+c*enddo*do,d,1,4*25,1E,53*25-1+i+d,53*25+i+d*enddo*do,e,1,1*25,1E,58*25-1+i+e,58*25+i+e*enddo*do,g,1,6*25,1E,61*25-1+i+g,61*25+i+g*enddo*do,h,1,1*25,1E,84*25-1+i+h,84*25+i+h*enddo!8截面type,1mat,1secnum,15*do,a,1,1*25,1E,44*25-1+i+a,44*25+i+a*enddo*endif*SET,j,j+14*enddo!********************************************* !!!FuGan!!!********************************************* /prep7*SET,j,-14*do,i,1,(94*25+1)*2,94*25+1k,i,0,j,0k,94*25+i,1363.000,j,0kFILL,i,94*25+i,93,,,,,j=j+14*enddo*SET,Nmax1,10000*SET,j,-14*do,i,Nmax1+1,Nmax1+(92*25+1)*2,92*25+1k,i,14.50,j,15.00k,92*25+i,1348.50,j,15.00kFILL,i,92*25+i,91,,,,,j=j+14*enddo!连接右斜线1j=Nmax1+1*do,n,1,(94*25+1)*2,94*25+1kk=0*do,i,n,n-1+20*25,25lstr,i,jkk=kk+1*if,kk,eq,1,thenlatt,1,1,1,,,,21esize,,25lmesh,all*elseif,kk,ge,4,and,kk,le,6,then latt,1,1,1,,,,21esize,,25lmesh,all*elseif,kk,ge,8,and,kk,le,9,then latt,1,1,1,,,,21esize,,25lmesh,all*elseif,kk,eq,7,thenlatt,1,1,1,,,,23esize,,25lmesh,all*elselatt,1,1,1,,,,23esize,,25lmesh,all*endifj=j+25*enddoj=j+2301-(20-0)*25*enddo!连接右斜线2j=Nmax1+1+40*25*do,n,1,(94*25+1)*2,94*25+1kk=0*do,i,n+40*25,n-1+60*25,25lstr,i,jkk=kk+1*if,kk,eq,20,thenlatt,1,1,1,,,,21esize,,25lmesh,all*elseif,kk,ge,3,and,kk,le,9,then latt,1,1,1,,,,21esize,,25lmesh,all*elseesize,,25lmesh,all*endifj=j+25*enddoj=j+2301-(60-40)*25*enddo!连接左斜线1j=Nmax1+19*25+1*do,n,1,(94*25+1)*2,94*25+1kk=0*do,i,n+21*25,n-1+41*25,25lstr,i,jkk=kk+1*if,kk,eq,1,thenlatt,1,1,1,,,,21esize,,25lmesh,all*elseif,kk,ge,13,and,kk,le,18,then latt,1,1,1,,,,21esize,,25lmesh,all*elselatt,1,1,1,,,,22esize,,25lmesh,all*endifj=j+25*enddoj=j+2301-(41-21)*25*enddo!连接左斜线2j=Nmax1+59*25+1*do,n,1,(94*25+1)*2,94*25+1kk=0*do,i,n+61*25,n-1+95*25,25lstr,i,jkk=kk+1*if,kk,ge,1,and,kk,le,12,then latt,1,1,1,,,,22esize,,25lmesh,all*elseif,kk,ge,21,and,kk,le,22,thenesize,,25lmesh,all*elseif,kk,ge,15,and,kk,le,16,then latt,1,1,1,,,,23esize,,25lmesh,all*elseif,kk,ge,24,and,kk,le,28,then latt,1,1,1,,,,23esize,,25lmesh,all*elseif,kk,ge,33,and,kk,le,34,then latt,1,1,1,,,,23esize,,25lmesh,all*elselatt,1,1,1,,,,21esize,,25lmesh,all*endifj=j+25*enddoj=j+2301-(95-61)*25*enddo!连接shu线j=Nmax1+1*do,n,1,(94*25+1)*2,94*25+1kk=0*do,i,n+1*25,n-1+94*25,25lstr,i,jkk=kk+1*if,kk,ge,1,and,kk,le,3,thenlatt,1,1,1,,,,22esize,,25lmesh,all*elseif,kk,ge,22,and,kk,le,58,then latt,1,1,1,,,,22esize,,25lmesh,all*elseif,kk,ge,62,and,kk,le,74,then latt,1,1,1,,,,22esize,,25lmesh,all*elseif,kk,ge,78,and,kk,le,83,thenesize,,25lmesh,all*elseif,kk,ge,89,and,kk,le,93,thenlatt,1,1,1,,,,22esize,,25lmesh,all*elseif,kk,eq,7,thenlatt,1,,1,,,,23esize,,25lmesh,all*elseif,kk,eq,60,thenlatt,1,1,1,,,,23esize,,25lmesh,all*elseif,kk,eq,75,thenlatt,1,1,1,,,,23esize,,25lmesh,all*elseif,kk,eq,86,thenlatt,1,1,1,,,,23esize,,25lmesh,all*elseif,kk,eq,87,thenlatt,1,1,1,,,,23esize,,25lmesh,all*elselatt,1,1,1,,,,21esize,,25lmesh,all*endifj=j+25*enddoj=j+2301-93*25*enddomodmsh,detackldele,allkdele,all!********************************************* !!!2X+上平联!!!********************************************* /prep7k,1,14.5,-14,15.00k,2,14.5,-7,15.00k,3,14.5,-14,10.2k,4,14.5,-7,10.2lstr,1,4lstr,2,3lstr,3,4lgen,2,all,,,,7,nummrg,kpnumcmp,kplstr,2,4lgen,93,all,,,14.5,,lgen,2,all,,,,14,,latt,1,,1,,,,24esize,,10lmesh,allmodmsh,detackldele,allkdele,all*SET,j,-14*do,i,1,(92*25+1)*2,92*25+1k,i,14.50,j,15.00k,92*25+i,1348.50,j,15.00kFILL,i,92*25+i,91,,,,,j=j+14*enddo*do,i,1,92*25+1,25lstr,i,i+92*25+1lstr,i+92*25+1,i+(92*25+1)*2*enddolatt,1,,1,,,,31esize,,20lmesh,allmodmsh,detackldele,all*do,i,1,92*25,25lstr,i,i+92*25+1+1*25lstr,i+92*25+1,i+(92*25+1)*2+1*25*enddo*do,i,1+1*25,92*25+1,25lstr,i,i+92*25+1-25lstr,i+92*25+1,i+(92*25+1)*2-25*enddonummrg,kpnumcmp,kplatt,1,,1,,,,32esize,,20lmesh,allmodmsh,detackldele,allkdele,all!************************************************! !!1#ta!!!************************************************ /prep7nn=400000n,1+nn,290,0,158.9n,13+nn,290,0,132.5n,20+nn,290,0,115n,22+nn,290,0,110.5n,63+nn,290,0,110.5n,52+nn,290,-19,-6n,62+nn,290,-11.5,-34.5n,93+nn,290,19,-6n,103+nn,290,11.5,-34.5n,104+nn,290,-19,-6n,107+nn,290,19,-6fill,1+nn,13+nn,11,fill,13+nn,20+nn,6,fill,20+nn,22+nn,1,fill,22+nn,52+nn,29,fill,52+nn,62+nn,9,fill,63+nn,93+nn,29,fill,93+nn,103+nn,9,fill,104+nn,107+nn,2,real,1type,1mat,2secnum,53*do,a,1+nn,12+nn,1E,a,a+1*enddotype,1mat,2secnum,52*do,b,13+nn,19+nn,1E,b,b+1*enddotype,1mat,2secnum,51*do,b,20+nn,21+nn,1E,b,b+1*enddotype,1mat,2secnum,46*do,b,22+nn,51+nn,1E,b,b+1*enddotype,1mat,2secnum,42*do,b,52+nn,61+nn,1E,b,b+1*enddotype,1mat,2secnum,46*do,b,63+nn,92+nn,1E,b,b+1*enddotype,1mat,2secnum,42*do,b,93+nn,102+nn,1E,b,b+1*enddotype,1mat,2secnum,54*do,b,104+nn,106+nn,1E,b,b+1*enddo!************************************************ !!2#ta!!************************************************* nn=450000n,1+nn,870,0,158.9n,13+nn,870,0,132.5n,20+nn,870,0,115n,22+nn,870,0,110.5n,63+nn,870,0,110.5n,52+nn,870,-19,-6n,62+nn,870,-11.5,-34.5n,93+nn,870,19,-6n,103+nn,870,11.5,-34.5n,104+nn,870,-19,-6n,107+nn,870,19,-6fill,1+nn,13+nn,11,fill,13+nn,20+nn,6,fill,20+nn,22+nn,1,fill,22+nn,52+nn,29,fill,52+nn,62+nn,9,fill,63+nn,93+nn,29,fill,93+nn,103+nn,9,fill,104+nn,107+nn,2,real,1type,1mat,2secnum,53*do,a,1+nn,12+nn,1E,a,a+1*enddotype,1mat,2secnum,52*do,b,13+nn,19+nn,1E,b,b+1*enddotype,1mat,2secnum,51*do,b,20+nn,21+nn,1E,b,b+1*enddotype,1mat,2secnum,46*do,b,22+nn,51+nn,1E,b,b+1*enddotype,1mat,2secnum,42*do,b,52+nn,61+nn,1E,b,b+1*enddotype,1mat,2secnum,46*do,b,63+nn,92+nn,1E,b,b+1*enddotype,1mat,2secnum,42*do,b,93+nn,102+nn,1E,b,b+1*enddotype,1mat,2secnum,54*do,b,104+nn,106+nn,1E,b,b+1*enddo!************************************************ !!link!!************************************************ /prep7nn=300000mm=400000!*****************************ES*do,j,nn+1,nn+2301*2+1,2301*do,a,j,10*25+j,25type,3mat,1real,3E,a,mm+3mm=mm+1*enddomm=400000*enddo*do,j,nn+1+11*25,nn+2301*3,2301 *do,a,j,2*25+j,25type,3mat,1real,4E,a,mm+14mm=mm+1*enddomm=400000*enddo*do,j,nn+1+14*25,nn+2301*3,2301 *do,a,j,2*25+j,25type,3mat,1real,5E,a,mm+17mm=mm+1*enddomm=400000*enddo*do,j,nn+1+17*25,nn+2301*3,2301 type,3mat,1real,6E,j,mm+20mm=400000*enddo!****************************EM mm=400000*do,j,nn+1+21*25,nn+2301*3,2301 type,3mat,1real,6E,j,mm+20mm=400000*enddo*do,j,nn+1+22*25,nn+2301*3,2301 *do,a,j,2*25+j,25type,3mat,1real,5。