ANSYS计算算例(结构-荷载法)

合集下载

ANSYS计算算例(结构-荷载法)

（3）直边墙的输入
第二段圆弧的终点即是直边墙的顶点，查出该点坐标为（-3.520017，-3.077515），并由此计算出直边墙底端点的坐标为（-3.520017，-7.507515）
GUI：Main Menu> Preprocessor> Create>Keypoints>In Active CS
计算中考虑结构自重对结构内力的影响，还应在此输入材料的密度。 Nonlinear-Density
4、模型建立（1）第一段圆弧的输入（圆心：0，－3.41，起点： 0，0，圆心角：45º ）
Menu menu：Preprocessor>Material Props>Modeling Create>Lines Arcs>By Cent & Radius
三、加载及约束
1、荷载计算及转化 (1) 荷载计算对深埋隧道用《隧规》推荐的统计法公式计算得：竖向围岩压力：q=54084.25Pa；水平围岩压力：e＝0.15 q=8112.64Pa (2) 匀布荷载转化为等效节点荷载将匀布荷载按如下步骤转化成等效节点力。 A、列出结构上所有节点坐标 GUI：Utility Menu>List>Nodes
拾取基底两节点
选取UX,UY并输入限制位移值0
3、惯性荷载加入惯性荷载是与质量相关的荷载，该处是加入重力加速度。为了利用惯性效果来模拟重力，应当在重力的相反方向施加惯性项，即在Y方向施加一个加速度荷载相当于模拟-Y方向的重力： GUI：Menu Main> Solution>Apply>Gravity
ANSYS计算算例
• 拱形结构算例 – 铁路隧道 – 公路隧道 – 地铁区间隧道(矿山法施工) 矩形框架结构 – 地下街 – 地下停车场 – 地铁车站圆形结构 – 地铁区间隧道(盾构法施工)

隧道ansys计算程序算例_荷载结构模式

选取新建铁路宜昌（宜）－万州（万）铁路线上的别岩槽隧道某断面.该断面设计单位采用的支护结构如图3-3所示。

为保证结构的安全性.采用了荷载—结构模型.利用ANSYS 对其进行计算分析。

主要参数如下：●隧道腰部和顶部衬砌厚度是65cm.隧道仰拱衬砌厚度为85cm。

●采用C30钢筋混凝土为衬砌材料。

●隧道围岩是Ⅳ级,洞跨是5.36米.深埋隧道。

●隧道仰拱下承受水压.水压0.2MPa。

图 3-3 隧道支护结构断面图隧道围岩级别是Ⅳ级.其物理力学指标及衬砌材料C30钢筋混凝土的物理力学指标见表3-3所示。

根据《铁路隧道设计规范》.可计算出深埋隧道围岩的垂直匀布力和水平匀布力。

对于竖向和水平的分布荷载.其等效节点力分别近似的取节点两相临单元水平或垂直投影长度的一般衬砌计算宽度这一面积范围内的分布荷载的总和。

自重荷载通过ANSYS程序直接添加密度施加。

隧道仰拱部受到的水压0.2MPa按照径向方向载置换为等效节点力,分解为水平竖直方向加载。

3.3.3 GUI操作方法3.3.3.1 创建物理环境1) 在【开始】菜单中依次选取【所有程序】/【ANSYS10.0】/【ANSYS Product Launcher】.得到“10.0ANSYS Product Launcher”对话框。

2）选中【File Management】.在“Working Directory”栏输入工作目录“D:\ansys\example301”.在“Job Name”栏输入文件名“Support”。

3）单击“RUN”按钮.进入ANSYS10.0的GUI操作界面。

4）过滤图形界面：Main Menu> Preferences.弹出“Preferences for GUI Filtering”对话框.选中“Structural”来对后面的分析进行菜单及相应的图形界面过滤。

5）定义工作标题：Utility Menu> File> Change Title.在弹出的对话框中输入“Tunnel Support Structural Analysis”.单击“OK”.如图3-4所示。

ANSYS计算结果与分析

ANSYS计算结果与分析一、有限元原理：有限元的解题思路可简述为：从结构的位移出发，通过寻找位移和应变，应变与应力，应力与内力，内力与外力的关系，建立相应的方程组，从而由已知的外力求出结构的内应力和位移。

有限元分析过程由其基本代数方程组成：[K]{V}={Q}，[K]为整个结构的刚变矩阵，{V}为未知位移量，{Q}为载荷向量。

这些量是不确定的，依靠所需解决的问题进行定量描述。

上述结构方程是通过应用边界条件，将结构离散化成小单元，从综合平衡方程中获得。

有限元是通过单元划分，在某种程度上模拟真实结构，并由数字对结构诸方面进行描述。

其描述的准确性依赖于单元细划的程度，载荷的真实性，材料力学参数的可信度，边界条件处理的正确程度。

本算例采用三角形六结点来划分单元。

二、有限元解题步骤：有限元的解题步骤为：①连续体的离散化；②选择单元位移函数；③建立单元刚度矩阵；④求解代数方程组，得到所有节点位移分量；⑤由节点位移求出内力或应力。

三、工程实例分析现已知一混凝土截面梁，长为L=2.4m，梁高为h=0.3m，梁宽设为单位宽度。

混凝土材料的各项属性为：容重γ=25KN/m3，E=2.4E10Pa，λ=0.2。

若该混凝土梁分别受到以下两种不同约束和不同受力的作用：（1）两端受固定约束作用，中间作用一个集中荷载P=10KN作用，如图A所示。

（2）一端受固定约束作用悬臂梁，梁上作用一均布荷载q=5KN/m作用，如图B 所示。

现要求使用有限元中的三角形六节点单元来计算两种情况下梁的位移与应力，并与力学计算结果进行比较和分析ANSYS分析过程（1）两端固定有限元模型Y方向位移图X方向应力图具体节点位移如下表：x应力表（单位：pa）NODE SX SY SZ SXY SYZ SXZ 1 -0.47117E+06 -94234. 0.0000 -71939. 0.00000.0000 2 -0.44824E+06 -43524. 0.0000 36804. 0.0000 0.0000 4 -0.26659E+06 24017.0.0000 11756. 0.0000 0.0000 6 -0.11092E+06 -1675.1 0.0000 -1416.0 0.0000 0.0000 8 26454. -287.29 0.0000 -529.21 0.0000 0.0000 10 0.14396E+06 -296.19 0.0000 -428.64 0.0000 0.0000 12 0.24665E+06 -154.80 0.0000 -931.20 0.0000 0.0000 14 0.31294E+06 -651.06 0.0000 1868.1 0.0000 0.0000 16 0.25514E+06 777.86 0.0000 2607.3 0.0000 0.0000 18 0.15259E+06 421.23 0.0000 955.79 0.0000 0.0000 20 36815. 350.40 0.0000 1367.3 0.0000 0.0000 22 -99169. 406.51 0.0000 1539.1 0.0000 0.0000 24 -0.25029E+06 709.47 0.0000 469.07 0.0000 0.0000 26 0.47021E+06 94043. 0.0000 71822. 0.0000 0.0000 28 -0.26849E+06 -41899. 0.0000 46711. 0.0000 0.0000 30 -0.12480E+06 -20589. 0.0000 44463. 0.0000 0.0000 32 2346.1 1441.3 0.0000 44301. 0.0000 0.0000 34 0.13084E+06 23544. 0.0000 46007. 0.0000 0.0000 36 0.27890E+06 44622. 0.0000 49721. 0.0000 0.0000 38 0.44740E+06 43435. 0.0000 -36738. 0.0000 0.000040 0.26576E+06 -23958. 0.0000 -11724. 0.0000 0.0000 42 0.11006E+06 1672.9 0.0000 1413.1 0.0000 0.0000 44 -27364. 285.26 0.0000 538.86 0.0000 0.0000 46 -0.14505E+06 297.97 0.0000 455.70 0.0000 0.0000 48 -0.23089E+06 828.18 0.0000 -2277.4 0.0000 0.0000 50 -0.35045E+06 -94932. 0.0000 4506.0 0.0000 0.0000 52 -0.23867E+06 40010. 0.0000 13774. 0.0000 0.0000 54 -0.15354E+06 -1579.2 0.0000 -2170.7 0.0000 0.0000 56 -37694. -287.51 0.0000 -1318.8 0.0000 0.0000 58 98366. -412.70 0.0000 -1539.9 0.0000 0.0000 60 0.24952E+06 -708.11 0.0000 -474.61 0.0000 0.0000 63 0.26768E+06 41770. 0.0000 -46653. 0.0000 0.0000 65 0.12400E+06 20454.0.0000 -44435. 0.0000 0.0000 67 -3155.0 -1579.4 0.0000 -44304. 0.0000 0.0000 69 -0.13167E+06 -23682. 0.0000 -46041. 0.0000 0.0000 71 -0.27977E+06 -44758. 0.0000 -49784. 0.0000 0.0000 151 0.18243E+06 -5692.4 0.0000 37694. 0.0000 0.0000 153 93763. -895.44 0.0000 60650. 0.0000 0.0000 155 2906.0 1964.2 0.0000 67290. 0.0000 0.0000 157 -87037. 5292.4 0.0000 58490. 0.0000 0.0000 159 -0.17211E+06 11178. 0.0000 31663. 0.0000 0.0000 167 73372. -851.10 0.0000 31158. 0.0000 0.0000 169 39292. -840.75 0.0000 49441. 0.0000 0.0000 171 3787.8 -337.99 0.0000 55813. 0.0000 0.0000 173 -32101. 100.68 0.0000 50074. 0.0000 0.0000 175 -67944. 717.92 0.0000 32029. 0.0000 0.0000 183 -17938. -922.76 0.0000 27092. 0.0000 0.0000 185 -6795.5 -612.85 0.0000 42946. 0.0000 0.0000 187 3071.2 18.206 0.0000 48306. 0.0000 0.0000 189 12803. 565.23 0.0000 43126. 0.0000 0.0000 191 22747. 538.90 0.0000 27571.0.0000 0.0000 199 -96298. -717.39 0.0000 22710. 0.0000 0.0000 201 -46228. -1019.6 0.0000 36320. 0.0000 0.0000 203 2626.4 121.91 0.0000 40665. 0.0000 0.0000 205 51239. 880.36 0.0000 36356. 0.0000 0.0000 207 0.10012E+06 736.28 0.0000 23263.0.0000 0.0000 215 -0.15923E+06 18494. 0.0000 11061. 0.0000 0.0000 217 -83947. 9075.8 0.0000 29103. 0.0000 0.0000 219 -2385.6 4224.8 0.0000 35294. 0.0000 0.0000 221 80991. 2496.3 0.0000 32516. 0.0000 0.0000 223 0.16576E+06 1612.2 0.0000 21754. 0.0000 0.0000 231 -0.20646E+06 -75438. 0.0000 -1738.1 0.0000 0.0000 233 -93006. -51133. 0.0000 -1579.6 0.0000 0.0000 235 6530.2 -30687. 0.0000 -666.930.0000 0.0000237 0.10345E+06 -14749. 0.0000 125.67 0.0000 0.0000 239 0.20441E+06 -4119.40.0000 970.85 0.0000 0.0000 247 -0.16755E+06 9862.5 0.0000 -17509. 0.0000 0.0000 249 -87809. 3440.9 0.0000 -31319. 0.0000 0.0000 251 -5584.0 978.66 0.0000 -35737.0.0000 0.0000 253 77562. 806.61 0.0000 -31968. 0.0000 0.0000 255 0.16243E+06 1086.9 0.0000 -20148. 0.0000 0.0000 263 -0.10061E+06 -689.88 0.0000 -23108.0.0000 0.0000 265 -51874. -100.28 0.0000 -35992. 0.0000 0.0000 267 -3478.2 629.38 0.0000 -40521. 0.0000 0.0000 269 45225. 1049.3 0.0000 -36167. 0.0000 0.0000 271 95246. 1060.4 0.0000 -22892. 0.0000 0.0000 279 -23597. -576.87 0.0000 -27627.0.0000 0.0000 281 -13606. -652.21 0.0000 -43170. 0.0000 0.0000 283 -3857.7 -89.957 0.0000 -48308. 0.0000 0.0000 285 5980.8 542.38 0.0000 -42912. 0.0000 0.0000 287 17073. 893.21 0.0000 -27066. 0.0000 0.0000 295 67130. -714.95 0.0000 -32028.0.0000 0.0000 297 31274. -97.707 0.0000 -50075. 0.0000 0.0000 299 -4624.8 340.23 0.0000 -55817. 0.0000 0.0000 301 -40137. 841.32 0.0000 -49444. 0.0000 0.0000 303 -74225. 850.89 0.0000 -31160. 0.0000 0.0000 311 0.17131E+06 -11152. 0.0000 -31681.0.0000 0.0000 313 86204. -5274.6 0.0000 -58500. 0.0000 0.0000 315 -3755.8 -1950.0 0.0000 -67295. 0.0000 0.0000 317 -94623. 911.23 0.0000 -60649. 0.0000 0.0000 319 -0.18330E+06 5714.3 0.0000 -37686. 0.0000 0.0000由以上分析结果可以得出：跨中最大挠度为：2.95E-05m 梁端上截面应力为：-0.35Mpa 跨中上截面应力: 0.47Mpa 跨中下截面应力为：-0.471Mpa 用材料力学进行校核：Wz=bh6222，左右杆端弯矩为：=ql12ql2，跨中弯矩为：2ql24 左右杆端截面正应力为：σ跨中截面正应力为：σ=ql2bh6ql=22ql222bh=0.32MPa 242bh6=4bh=0.47MPa由图乘法求跨中截面的挠度，具体的计算公式如下：W===1EI11EI(412⨯1112-ql2⨯l2⨯1l2⨯)23-12⨯124ql2⨯l2⨯l2⨯13-23⨯l2⨯132ql2⨯l2⨯12)ql(ql41576144-384374EI=2.95E-05m（2）一端固定一端自由Y方向位移图X方向应力图具体节点位移如下表：x应力表（单位：pa）NODE SX SY SZ SXY SYZ SXZ 1 -0.25013E+07-0.50026E+06 0.0000 -0.23536E+06 0.0000 0.0000 2 2344.3 55.551 0.0000 1381.4 0.0000 0.0000 4 -0.20145E+070.13137E+06 0.0000 42962. 0.0000 0.0000 6 -0.16765E+07 -5697.0 0.0000 -4238.2 0.0000 0.0000 8 -0.13587E+07 -602.67 0.0000 -1529.1 0.0000 0.0000 10 -0.10740E+07 -762.28 0.0000 -1374.8 0.0000 0.0000 12 -0.82262E+06 -655.05 0.0000 -1117.5 0.0000 0.0000 14 -0.60460E+06 -562.14 0.0000 -863.89 0.0000 0.0000 16 -0.41991E+06 -468.62 0.0000 -610.04 0.0000 0.0000 18 -0.26856E+06 -375.13 0.0000 -356.22 0.0000 0.0000 20 -0.15053E+06 -281.65 0.0000 -102.36 0.0000 0.0000 22 -65843. -187.66 0.0000 150.44 0.0000 0.0000 24 -14504. -102.05 0.0000 417.31 0.0000 0.0000 26 -3100.4 -5193.8 0.0000 775.10 0.0000 0.0000 28 742.77 684.07 0.0000 633.96 0.0000 0.0000 30 -17.938 -317.56 0.0000 -285.27 0.0000 0.0000 32 -477.64 -2424.3 0.0000 -770.67 0.0000 0.000034 -781.28 -4602.3 0.0000 -570.42 0.0000 0.0000 36 -1125.5 -5696.4 0.0000 -174.96 0.0000 0.0000 38 0.24391E+07 0.24437E+06 0.0000 -0.11805E+06 0.0000 0.0000 40 11812. -4996.7 0.0000 -1071.1 0.0000 0.0000 42 60425. -5175.6 0.0000 -1154.9 0.0000 0.0000 44 0.14240E+06 -5280.4 0.0000 -1421.3 0.0000 0.0000 46 0.25771E+06 -5372.8 0.0000 -1674.3 0.0000 0.0000 48 0.40636E+06 -5466.3 0.0000 -1928.2 0.0000 0.0000 50 0.58834E+06 -5559.8 0.0000 -2182.0 0.0000 0.0000 52 0.80365E+06 -5653.3 0.0000 -2435.9 0.0000 0.0000 54 0.10523E+07 -5746.8 0.0000 -2689.9 0.0000 0.0000 56 0.13342E+07 -5842.7 0.0000 -2932.1 0.0000 0.0000 58 0.16503E+07 -5879.4 0.0000 -3376.0 0.0000 0.0000 60 0.19778E+07 -7803.0 0.0000 2914.3 0.0000 0.0000 63 0.15354E+07 0.22889E+06 0.0000 -0.10713E+06 0.0000 0.0000 650.73975E+06 0.11355E+06 0.0000 -70967. 0.0000 0.0000 67 -7600.7 -6004.9 0.0000 -61499. 0.0000 0.0000 69 -0.75765E+06-0.12447E+06 0.0000 -78773. 0.0000 0.0000 71 -0.15608E+07-0.23582E+06 0.0000 -0.12284E+06 0.0000 0.0000 151 9798.9 -5221.2 0.0000 -7608.6 0.0000 0.0000 153 4996.1 -4256.1 0.0000 -11467. 0.0000 0.0000 155 -911.85 -2609.3 0.0000 -12699. 0.0000 0.0000 157 -6764.9 -887.12 0.0000 -11142.0.0000 0.0000 159 -11230. 312.08 0.0000 -6755.7 0.0000 0.0000 167 42347. -5311.7 0.0000 -14833. 0.0000 0.0000 169 20742. -4334.3 0.0000 -22836. 0.0000 0.0000 171 -1834.5 -2501.8 0.0000 -25426. 0.0000 0.0000 173 -24341. -653.58 0.0000 -22464.0.0000 0.0000 175 -45318. 505.18 0.0000 -13939. 0.0000 0.0000 183 97163. -5261.6 0.0000 -22022. 0.0000 0.0000 185 47621. -4327.6 0.0000 -34159. 0.0000 0.0000 187 -2750.4 -2507.1 0.0000 -38140. 0.0000 0.0000 189 -53018. -662.65 0.0000 -33786.0.0000 0.0000 191 -0.10162E+06 539.25 0.0000 -21132. 0.0000 0.0000 1990.17420E+06 -5215.6 0.0000 -29214. 0.0000 0.0000 201 85610. -4326.4 0.0000 -45483.0.0000 0.0000 203 -3667.2 -2509.3 0.0000 -50853. 0.0000 0.0000 205 -92806. -660.75 0.0000 -45110. 0.0000 0.0000 207 -0.18014E+06 586.11 0.0000 -28324. 0.0000 0.0000 215 0.27346E+06 -5169.3 0.0000 -36406. 0.0000 0.0000 217 0.13471E+06 -4324.80.0000 -56807. 0.0000 0.0000 219 -4584.1 -2511.7 0.0000 -63566. 0.0000 0.0000 221 -0.14371E+06 -659.14 0.0000 -56434. 0.0000 0.0000 223 -0.28088E+06 632.39 0.0000 -35517. 0.0000 0.0000231 0.39494E+06 -5123.0 0.0000 -43599. 0.0000 0.0000 233 0.19492E+06 -4323.20.0000 -68131. 0.0000 0.0000 235 -5500.9 -2514.0 0.0000 -76279. 0.0000 0.0000 237 -0.20572E+06 -657.53 0.0000 -67758. 0.0000 0.0000 239 -0.40384E+06 678.71 0.0000 -42709. 0.0000 0.0000 247 0.53864E+06 -5076.5 0.0000 -50791. 0.0000 0.0000 249 0.26624E+06 -4321.2 0.0000 -79456. 0.0000 0.0000 251 -6417.7 -2516.1 0.0000 -88993. 0.0000 0.0000 253 -0.27884E+06 -655.95 0.0000 -79083. 0.0000 0.0000 255 -0.54903E+06 724.90 0.0000 -49901. 0.0000 0.0000 263 0.70456E+06 -5031.3 0.0000 -57983. 0.0000 0.0000 265 0.34868E+06 -4322.8 0.0000 -90782. 0.0000 0.0000 267 -7336.1 -2519.0 0.0000 -0.10171E+06 0.0000 0.0000 269 -0.36307E+06 -649.67 0.0000 -90406. 0.0000 0.0000 271 -0.71644E+06 772.81 0.0000 -57090. 0.0000 0.0000 279 0.89269E+06 -5018.3 0.0000 -65150. 0.0000 0.0000 281 0.44225E+06 -4346.5 0.0000 -0.10206E+06 0.0000 0.0000 283 -8193.7 -2657.9 0.0000 -0.11441E+06 0.0000 0.0000 285 -0.45838E+06 -751.50 0.0000 -0.10180E+06 0.0000 0.0000 287 -0.90611E+06 819.40 0.0000 -64358. 0.0000 0.0000 295 0.11033E+07 -4700.3 0.0000 -72715. 0.0000 0.0000 297 0.54631E+06 -1332.2 0.0000 -0.11431E+06 0.0000 0.0000 299 -10032. -610.15 0.0000 -0.12657E+06 0.0000 0.0000 301 -0.56488E+06 156.41 0.0000 -0.11184E+06 0.0000 0.0000由以上分析结果可一得出：梁端最大挠度为：9.72E-04m 梁端截面最大应力为：-2.5Mpa 用材料力学进行校核：Wz=bh62，杆端弯矩为：FL =ql22左右杆端截面正应力为：σ固端截面正应力为：σ=6bh==ql222bh=0.32MPa FLbh26FL2bh2226=2.4MPa 左右杆端截面正应力为：σ=ql6bh=ql222bh=0.32MPa由图乘法可知自由端的挠度为：W=1EI(12⨯FL⨯L⨯23L)=1FL3EI3=9.60E-04m结论在对本工程进行ANSYS有限元数值分析过程中，作者采用的单元形式为三角形六节点单元PLANE2单元，因其为平面单元，ANSYS计算过程中没有输入梁的宽度，其计算默认的梁宽为一个单位。

Ansys 荷载组合

Ansys 荷载组合1,几何模型(beam3和beam54)建立后，定义所需的element table，主要包括杆端力和最大应力，最小应力等。

然后保存数据库。

分别施加四种荷载的标准值（不乘分项系数），并分别存成四个load step file。

2,使用solution－>from ls files，求解四种荷载3,荷载组合，命令流如下：/post1lcdef,1,1lcdef,2,2lcdef,3,3lcdef,4,4 !定义四种工况，分别为四种荷载下的计算结果lcfact,1,1.2lcfact,2,1.4lcfact,3,1.19lcfact,4,1.4 !指定各工况的组合系数lcase,1 !读入工况1，database＝1sumtype,prin !指定加操作的对象lcoper,add,2 !荷载组合，database＝database＋2lcoper,add,4 !荷载组合，database＝database＋4lcoper,lprin !计算线性主应力lcwrite,11 !把database结果写到工况11,即恒荷载＋活荷载＋吊车荷载的结果lcase,1lcfact,2,1.19lcfact,4,1.19 !改变组合系数sumtype,prinlcoper,add,2lcoper,add,3lcoper,add,4lcoper,lprinlcwrite,12 !把database结果写到工况12,即恒荷载＋活荷载＋吊车荷载＋风荷载的结果!... ...其他荷载组合!之后使用lcase,n 就可调入工况n，并查看它的变形和内力!可使用如下命令流得到工况11和12，13的较大者99,进而查看最大应力lcase,11lcase,min,12lcase,min,13lcwrite,98lcase 98!查看工况98的应力分布... ...lcase,11lcase,max,12lcase,max,13lcwrite,99lcase 99!查看工况99的应力分布... ...以下为定义和读取荷载工况用到的一些命令：LCDEF_从结果文件中的一列结果产生荷载工况LCDEF, LCNO, LSTEP, SBSTEP, KIMGLCNO：随意的指针数（1－99），要赋给LSTEP，SBSTEP和FILE命令指定的荷载工况。

钢桁架ANSYS计算方法

桁架计算跨度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
-F
13
-78977
-1.96 F
14
0.0000
0
15
26390
0.65 F
16
-40340
-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
钢桁架计算
目的：计算钢桁架中，每根杆的内力。
方法：利用ANSYS有限元分析软件求解计算。
步骤：
1、进入前处理模块
2、定义单元为二维杆单元LINK1
3、定义材料弹性模量EX=210E9，泊松比0.33（Q235常用值）
Q235弹性模量E：196~216GPa
泊松比0.25~0.33
4、定义各个节点位置
-4.05 F
4
-0.16345E+06
-4.05 F
5
-95681
-2.37 F
6
-95681
-2.37 F
7
0.14643E+06
3.63 F
8

同济大学Ansys隧道荷载结构法算例

结构混凝土参数混凝土标号弹性模量eknmkpa惯性模量im25号混凝土2850000000000180600ansys计算在ansys中需要将荷载转化为节点荷载需要输出ansys结点坐标用excel计算等效节点荷载
隧道荷载结构法算例 1．问题描述：
主要截面尺寸如图 1 所示， R 1 5.2m ， R2 8.8m 。内轮廓采用曲墙拱三心圆拱。弹性抗力系数取为 k 1.5 10 kN / m 。衬砌厚度为 0.6m，截面的物理参数如表所示。
建立局部柱坐标系，施加径向弹簧（用 PSPRNG 命令）。根据计算结果中弹簧是否收拉，对结果进行迭代计算。第 1 次计算结果如下图所示。
图5
结构变形图
删除受拉弹簧后的计算结果如下图。
图6 命令流： /prep7 et,1,beam3 mp,ex,1,2.85e10 mp,prxy,1,0.25 r,1,0.6,0.018,0.6 k,1,0,5.2 k,2,-5.2,0 k,3,-5.0797,-1.45 k,4,0,0 k,5,3.6,0 larc,1,2,4,5.2 larc,2,3,5,8.8 lesize,1,,,8 lesize,2,,,2 lmesh,all /pnum,node,1 nplot nlist !compute nodal force in excel f,1,fx,1.249456357 f,3,fx,4.949809522 f,4,fx,9.709400641
finish /solu solve finish /post1 pldisp,1 !check springs with traction (+ means traction) etable,sforce,smisc,1 pretab,sforce !delete springs with tractions and recalculate here ! etable,m1,smisc,6 etable,m2,smisc,12 plls,m1,m2,-1

ansys motion 载荷计算

ansys motion 载荷计算ANSYS Motion是一种多体动力学仿真软件，用于模拟和分析机械系统的运动和载荷。

通过ANSYS Motion，用户可以预测和优化产品的性能、可靠性和耐久性，从而降低开发成本和时间。

本文将介绍ANSYS Motion在载荷计算方面的应用。

一、概述在机械系统设计和分析过程中，载荷计算是非常重要的一环。

它可以帮助工程师了解系统在实际工作条件下所承受的力和力矩，从而评估系统的可靠性和耐久性。

ANSYS Motion提供了一系列强大的工具和功能，用于准确计算和分析机械系统的载荷。

二、载荷类型在ANSYS Motion中，可以对机械系统施加多种类型的载荷。

常见的载荷类型包括静态载荷、动态载荷和惯性载荷。

静态载荷是指在系统静止时施加的力和力矩，动态载荷是指在系统运动时施加的力和力矩，惯性载荷是指由于系统加速度而产生的力和力矩。

三、载荷计算方法ANSYS Motion提供了多种载荷计算方法，包括直接计算、求解拉格朗日方程和使用约束等。

直接计算方法适用于简单的系统，通过施加已知的力和力矩来计算载荷。

拉格朗日方程方法适用于复杂的系统，通过求解拉格朗日方程来计算载荷。

约束方法适用于特殊的系统，通过施加约束来计算载荷。

四、载荷分析通过ANSYS Motion的载荷计算功能，用户可以对机械系统进行载荷分析。

载荷分析可以帮助工程师了解系统在不同工况下的载荷分布、载荷大小以及载荷的变化规律。

通过对不同工况下的载荷进行分析，用户可以评估系统的可靠性和耐久性，并优化系统的设计。

五、应用案例ANSYS Motion广泛应用于各个行业的机械系统设计和分析。

例如，在汽车行业中，可以使用ANSYS Motion对汽车悬挂系统进行载荷计算，评估悬挂系统在不同路况下的载荷分布和变化规律，从而优化悬挂系统的设计。

在航空航天行业中，可以使用ANSYS Motion 对飞机起落架进行载荷计算，评估起落架在起飞和降落过程中的载荷大小和变化规律，从而提高起落架的可靠性和耐久性。

钢架模型ansys计算

钢架杆件横截面面积均为3250mm2，荷载施加在桁架上，如图1所示，材料弹性模量E=2.1×105Mpa，泊松比μ=0.3.密度 =7.8×103Kg/m3。

下面是ansys的分析过程。

图1 桁架尺寸及荷载情况1 进入ansys2 设置不同日期和时间3选择单元模型4定义材料参数5定义实常数6定义水平杆间的长度和杆的高度7生成几何模型Step1生成3个关键点Step2建立3条线Step3复制线Step4对模型进行分析Step5设置所有单元线上划分单元的个数Step6对所选择线进行单元划分8对模型施加边界条件Step1对关键点1、9施加固定边界条件Step2对关键点3、6施加边界条件图2 模型形状9分析计算10结果的一般显示（1）图3 变形位移显示（2）图4对线形单元按实体效果进行显示（以3倍比例）（3）图5 Y方向的位移云图11 线性单元内力结果显示Step1 定义线性单元I、J节点的轴力Step2 画出单元的受力图图6 轴力情况图轴力结果如下（2）单元计算结果单元轴向力：LOAD STEP= 1 SUBSTEP= 1 TIME= 1.0000 LOAD CASE= 0ELEM SMIS11 -0.26942E+062 0.13470E+063 0.26942E+064 -26942.5 0.28287E+066 26942.7 0.29637E+068 0.14817E+069 -0.29637E+0610 -0.26940E+0611 -0.29634E+06MINIMUM V ALUESELEM 9V ALUE -0.29637E+06MAXIMUM V ALUESELEM 7V ALUE 0.29637E+06单元支反力结果：LOAD STEP= 1 SUBSTEP= 1TIME= 1.0000 LOAD CASE= 0THE FOLLOWING X,Y,Z SOLUTIONS ARE IN THE GLOBAL COORDINATE SYSTEMNODE FX FY1 0.20373E-09 0.23333E+067 0.25667E+06TOTAL V ALUESV ALUE 0.20373E-09 0.49000E+06获取支撑节点1处的y方向支反力，赋给Ry 1：获取支撑节点7处的Y方向支反力，赋值给Ry7计算结果如下：ABBREVIATION STATUS-ABBREV STRINGSA VE_DB SA VERESUM_DB RESUMEQUIT Fnc_/EXITPOWRGRPH Fnc_/GRAPHICSPARAMETER STATUS- ( 25 PARAMETERS DEFINED)(INCLUDING 19 INTERNAL PARAMETERS)NAME V ALUE TYPE DIMENSIONSA 3.60000000 SCALARB 3.11800000 SCALARI 1.00000000 SCALARMAXLAYER 0.00000000 SCALARRY1 233300.000 SCALARRY7 256700.000 SCALAR12 退出系统下面是Log文件/BATCH/COM,ANSYS RELEASE 11.0 UP20070125 18:53:36 11/03/2012 /input,menust,tmp,'',,,,,,,,,,,,,,,,1/GRA,POWER/GST,ON/PLO,INFO,3/GRO,CURL,ON/CPLANE,1/REPLOT,RESIZEWPSTYLE,,,,,,,,0/PLOPTS,INFO,3/PLOPTS,LEG1,1/PLOPTS,LEG2,1/PLOPTS,LEG3,1/PLOPTS,FRAME,1/PLOPTS,TITLE,1/PLOPTS,MINM,1/PLOPTS,FILE,0/PLOPTS,LOGO,1/PLOPTS,WINS,1/PLOPTS,WP,0/PLOPTS,DATE,0/TRIAD,ORIG/REPLOT!*/PREP7!*ET,1,MASS21!*SA VE!* MPTEMP,,,,,,,,MPTEMP,1,0MPDATA,EX,1,,2100000000000 MPDATA,PRXY,1,,0.3 MPTEMP,,,,,,,,MPTEMP,1,0MPDE,EX,1MPDE,PRXY,1MPDATA,EX,1,,210000000000 MPDATA,PRXY,1,,0.3 MPTEMP,,,,,,,,MPTEMP,1,0MPDATA,DENS,1,,7800!*!*ETDEL,1!*ET,1,LINK1!*!*SA VER,1,0.00325, ,!*SA VE*SET,a,3.6*SET,b,3.118SA VEK,1,0,0,0,K,2,a/2,B,0,RESUME/COM,ANSYS RELEASE 11.0 UP20070125 18:59:52 11/03/2012 /AUTO,1/REP,FAST/PNUM,KP,1/PNUM,LINE,1/PNUM,AREA,0/PNUM,VOLU,0/PNUM,NODE,0/PNUM,TABN,0/PNUM,SV AL,0/NUMBER,0!*/PNUM,ELEM,0/REPLOT!*K,1,0,0,0,K,2,a/2,B,0,K,3,a,0,0,SA VELSTR, 1, 2LSTR, 1, 3LSTR, 2, 3SA VEFLST,3,3,4,ORDE,2FITEM,3,1FITEM,3,-3LGEN,3,P51X, , ,a, , , ,0FLST,3,2,4,ORDE,2 FITEM,3,2FITEM,3,5LGEN,2,P51X, , ,a/2,b, , ,0 SA VENUMMRG,ALL, , , ,LOW !*LESIZE,ALL, , ,1, ,1, , ,1, FLST,2,11,4,ORDE,2 FITEM,2,1FITEM,2,-11LMESH,P51XFINISH/SOLFLST,2,1,3,ORDE,1 FITEM,2,1!*/GODK,P51X, , , ,0,ALL, , , , , , FLST,2,1,3,ORDE,1 FITEM,2,9!*/GODK,P51X, , , ,0,UY, , , , , , SA VE/DIST,1,0.924021086472,1 /REP,FASTSA VEFLST,2,1,3,ORDE,1 FITEM,2,3!*/GOFK,P51X,FY,-210e3 FLST,2,1,3,ORDE,1 FITEM,2,6!*/GOFK,P51X,FY,-280000SA VE/STATUS,SOLUSOLVEFINISH/POST1PLDISP,1/DIST,1,1.08222638492,1 /REP,FAST/DIST,1,1.08222638492,1 /REP,FAST/DIST,1,0.924021086472,1 /REP,FAST/DIST,1,0.924021086472,1 /REP,FAST/DIST,1,1.08222638492,1 /REP,FAST/DIST,1,1.08222638492,1 /REP,FAST/DIST,1,1.08222638492,1 /REP,FAST/DIST,1,1.08222638492,1 /REP,FAST/DIST,1,1.08222638492,1/REP,FAST/DIST,1,0.924021086472,1/REP,FAST/DIST,1,0.924021086472,1/REP,FAST/DIST,1,0.924021086472,1/REP,FAST/DIST,1,1.08222638492,1/REP,FAST/DIST,1,1.08222638492,1/REP,FAST/DIST,1,0.924021086472,1/REP,FAST/DIST,1,0.924021086472,1/REP,FAST/DIST,1,1.08222638492,1/REP,FAST/USER, 1/FOC, 1, 6.32209400502 , 0.584690326081 , 0.00000000000 /REPLO/DIST,1,1.08222638492,1/REP,FAST/DIST,1,0.924021086472,1/REP,FAST!*/SHRINK,0/ESHAPE,3/EFACET,1/RATIO,1,1,1/CFORMAT,32,0/REPLOT!*!*/EFACET,1PLNSOL, U,Y, 0,1.0/DIST,1,1.08222638492,1 /REP,FAST/DIST,1,1.08222638492,1 /REP,FAST/DIST,1,1.08222638492,1 /REP,FAST/DIST,1,1.08222638492,1 /REP,FAST/DIST,1,1.08222638492,1 /REP,FAST/DIST,1,0.924021086472,1 /REP,FAST/DIST,1,0.924021086472,1 /REP,FAST/DIST,1,0.924021086472,1 /REP,FAST/DIST,1,1.08222638492,1 /REP,FAST/DIST,1,0.924021086472,1 /REP,FAST/DIST,1,0.924021086472,1 /REP,FAST/DIST,1,1.08222638492,1 /REP,FAST/DIST,1,0.924021086472,1/REP,FAST/DIST,1,0.924021086472,1/REP,FAST/DIST,1,1.08222638492,1/REP,FAST/DIST,1,1.08222638492,1/REP,FAST/DIST,1,1.08222638492,1/REP,FAST/REPLOT,RESIZE/REPLOT,RESIZE/REPLOT,RESIZE/REPLOT,RESIZESA VEA VPRIN,0, ,ETABLE,bridge_I,SMISC, 1!*A VPRIN,0, ,ETABLE,bridge_J,SMISC, 1!*A VPRIN,0, ,ETABLE, ,U,X!*ETABLE,,ERASE,3PLLS,BRIDGE_I,BRIDGE_J,0.5,0 /REPLOT,RESIZE/DIST,1,0.924021086472,1/REP,FAST/DIST,1,1.08222638492,1/REP,FAST/DIST,1,1.08222638492,1 /REP,FAST/DIST,1,1.08222638492,1 /REP,FAST/DIST,1,0.924021086472,1 /REP,FAST/DIST,1,1.08222638492,1 /REP,FAST!*PRESOL,SMISC,1!*PRRSOL,PRRSOL,FY*SET,Ry1,233300 PRRSOL,FY*SET,Ry7,256700*STATFINISH! /EXIT,ALL。

两个荷载结构法ansys命令流实例

实例一:finish/clear/COM,Structural/prep7*AFUN,deg*set,Py,275461.0591*set,px1,71390*set,px2,116230.464*set,CQHD,0.5*set,CQDYCD,0.6*set,CQETXML,31e9*set,CQUBSB,0.2*set,WYTXKL,200e6*set,WYTXKLZGQ, 850e6 *set,CQMD,2300et,1,beam3r,1,0.5,0.5*0.5*0.5/12,0.5 mp,ex,1, CQETXMLmp,dens,1, CQMDmp,prxy,1, CQUBSBk,1,311.419293,-97.854439 k,2,308.756764,-98.346462 k,3,306.445920,-99.757543 k,4,304.948061,-101.812347 k,5,304.418078,-104.299297 k,6,304.452492,-105.061153 k,7,304.555456,-105.816802 k,8,304.958136,-106.652780 k,9,305.731611,-107.165369 k,10,308.539368,-107.845632 k,11,311.419293,-108.074297 k,12,314.299218,-107.845632 k,13,317.106975,-107.165369 k,14,317.880450,-106.652780 k,15,318.283130,-105.816802k,16,318.386093,-105.061153k,17,318.420508,-104.299297k,18,317.890524,-101.812347k,19,316.392666,-99.757543k,20,314.081822,-98.346462!二衬轴线larc,1,3,2larc,3,5,4larc,5,7,6larc,7,9,8larc,9,11,10larc,11,13,12larc,13,15,14larc,15,17,16larc,17,19,18larc,19,1,20allsellsel,s,loc,x,311.419293,302.247076lcomb,allnummrg,alllesize,all,,,50allsellsel,s,loc,x,311.419293,320.122358lcomb,allnummrg,alllesize,all,,,50allslmesh,alleplotnummrg,all,,,,low !所有号数从1开始numcmp,node,eord !压缩节点号并排序@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@!施加水平弹簧支撑单元*do,i,16,50a1=abs(ny(i-1)-ny(i+1))*WYTXKLkx=a1/2PSPRNG,i,TRAN,kx,-0.3,,,,*enddo*do,i,52,61a1=abs(ny(i-1)-ny(i+1))*WYTXKLkx=a1/2PSPRNG,i,TRAN,kx,0.3,,,,*enddo*do,i,62,79a1=abs(ny(i-1)-ny(i+1))*WYTXKLZGQ kx=a1/2PSPRNG,i,TRAN,kx,0.3,,,,*enddo*do,i,80,86a1=abs(ny(i-1)-ny(i+1))*WYTXKLkx=a1/2PSPRNG,i,TRAN,kx,0.3,,,,*enddo!施加竖直弹簧支撑单元*do,i,29,50b1=abs(nx(i-1)-nx(i+1))*WYTXKLky=b1/2PSPRNG,i,TRAN,ky,,-0.3,,,*enddob1=abs(nx(52)-nx(50))*WYTXKLky=b1/2PSPRNG,51,TRAN,ky,,-0.3,,,*do,i,52,61b1=abs(nx(i-1)-nx(i+1))*WYTXKLky=b1/2PSPRNG,i,TRAN,ky,,-0.3,,,*enddo*do,i,62,73b1=abs(nx(i-1)-nx(i+1))*WYTXKLZGQky=b1/2PSPRNG,i,TRAN,ky,,-0.3,,,*enddo!施加边界条件与围岩压力/solunsel,s,,,51 !选择约束的节点,仰拱中间节点d,all,ux !施加水平方向的约束allselfcum,add,,!竖向荷载*do,i,2,28a=-Py*0.5*abs(nx(i-1)-nx(i+1))f,i,fy,a*enddo*do,i,80,99a=-Py*0.5*abs(nx(i-1)-nx(i+1))f,i,fy,a*enddoa=-Py*0.5*abs(nx(1)-nx(99))f,100,fy,aa=-Py*0.5*abs(nx(2)-nx(100))f,1,fy,a!水平向荷载（X方向压力，全部节点施加水平荷载，换算成节点荷载，在节点上施加）*do,i,52,61a=-(px1+(px2-px1)/(ny(51)-ny(1))*(ny(51)-ny(i)))*abs(nY(i-1)-nY(i+1))*0.5f,i,fX,a*enddo*do,i,80,99a=-(px1+(px2-px1)/(ny(51)-ny(1))*(ny(51)-ny(i)))*abs(nY(i-1)-nY(i+1))*0.5f,i,fX,a*enddoa=-(px1+(px2-px1)/(ny(51)-ny(1))*(ny(51)-ny(100)))*abs(ny(1)-ny(99))*0.5f,100,fX,a*do,i,2,50a=(px1+(px2-px1)/(ny(51)-ny(1))*(ny(51)-ny(i)))*abs(nY(i-1)-nY(i+1))*0.5f,i,fX,a*enddoacel,9.8,solvefinish!进入后处理，定义荷载工况并组合，输出图片和文本文件/post1ESEL,s,TYPE,,1 !选择单元1（隧道轴力中心线单元）NSLE,s,ALL !选择所有节点etable,mi,smisc,6etable,mj,smisc,12plls,mi,mj,-1etable,fi,smisc,1etable,fj,smisc,7plls,fi,fj,1pldisp,1etable,reflpretab, mi, mj, fi, fj实例二:/prep7*set,shuxiang,276170 !竖向荷载*set,zuoce,125500 !水平荷载*set,youce,-125500*set,jizhongli,-100000 !地震集中力*set,uxishu,136.966e6*0.5 !水平弹性抗力系数*set,vxishu,87.4348e6*0.5 !竖向弹性抗力系数!定义单元类型、实常数及材料性质et,1,beam3r,1,0.3,0.00225,0.3mp,ex,1,3.1e10 !弹性模量mp,dens,1,2500 !密度mp,prxy,1,0.167 !泊松比!创建几何模型画点k,1,0,2.55k,2,-2.04,1.53k,3,-2.6,-0.15k,4,-2.319721,-2.023237k,5,-1.648852,-2.71152k,6,0,-2.98k,7,1.648852,-2.71152k,8,2.319721,-2.023237k,9,2.6,-0.15k,10,2.04,1.53k,1000,0,0 !中心点!画隧道轮廓线larc,1,2,1000,2.55larc,2,3,1000,2.8larc,3,4,1000,6.4larc,4,5,1000,1.04969larc,5,6,1000,5.2larc,6,7,1000,5.2larc,7,8,1000,1.04969larc,8,9,1000,6.4larc,9,10,1000,2.8larc,10,1,1000,2.55!选择左边的所有线allsel !选择所有的实体lsel,s,loc,x,0,-5 !选择左边的所有线:x从0到－5 lcomb,all !合并所选的线lesize,all,,,30 !把上面合并的线等分30段!选择右边的所有线allsellsel,s,loc,x,0,5lcomb,alllesize,all,,,30allsel !全选!划分单元lmesh,allnplot !显示节点!加弹簧!y方向*do,i,23,40PSPRNG,i,TRAN,vxishu,,-0.3,, , !tran-直线的*enddoPSPRNG,2,TRAN,vxishu,,-0.3,, , !固定结构用的!x方向*do,i,15,26PSPRNG,i,TRAN,uxishu,-0.3,,, ,*enddo*do,i,37,48PSPRNG,i,TRAN,uxishu,0.3,,, ,*enddoallselfinish!进入求解层，施加荷载，定义荷载步等/solu!施加约束nsel,s,,,2 !选择约束的节点nsel,s,loc,x,0dnsel,r,loc,y,-2.98d,all,ux !施加水平方向的约束allselfcum,add,, !一定要，使荷载能叠加!竖向荷载*do,i,3,15a=-shuxiang*0.5*abs(nx(i)-nx(i+1))f,i+1,fy,af,i,fy,a*enddo*do,i,47,59a=-shuxiang*0.5*abs(nx(i)-nx(i+1))f,i+1,fy,af,i,fy,a*enddo!对未循环的节点施加竖向荷载a=-shuxiang*0.5*abs(nx(1)-nx(3))f,3,fy,af,1,fy,2*af,60,fy,a!水平荷载*do,i,3,30a=zuoce*0.5*abs(ny(i)-ny(i+1))f,i+1,fx,af,i,fx,a*enddo*do,i,32,59a=youce*0.5*abs(ny(i)-ny(i+1))f,i+1,fx,af,i,fx,a*enddof,1,fx,jizhongli !加集中力!对未循环的节点施加水平荷载a=youce*0.5*abs(ny(31)-ny(2))f,31,fx,-af,32,fx,a!施加重力acel,,9.8solve !求解finish!进入后处理，定义荷载工况并组合，输出图片和文本文件/post1!显示弯矩图etable,mi,smisc,6etable,mj,smisc,12plls,mi,mj,-1!显示轴力图etable,fi,smisc,1etable,fj,smisc,7plls,fi,fj,1!显示变形图pldisp,1。

ANSYS软件进行有限元计算实例

ANSYS软件进行有限元计算实例工字钢梁结构静力分析一工字钢梁两端均为固定端，其截面尺寸为：l=1.0m，a=0.16m，b=0.2m，c=0.02m，d=0.03m。

试建立该工字钢梁的三维实体模型，并在考虑重力的情况下对其进行结构静力分析。

其他已知参数如下：弹性模量E=206GPa；泊松比μ=0.3；材料密度ρ=7800kg/m3；重力加速度g=9.8m/s2；作用力作用于梁的上表面沿长度方向的中线处，其大小为F y=-5000N。

1)单元类型、几何特性、材料特性定义a）定义单元类型：Main Menu: Preprocessor→Element Type→Add/Edit/Delete弹出对话框，单击对话框中的“Add…”按钮，又弹出一对话框，选中其中的“Solid”和“Brick 8node 45”选项，单击“OK”按钮，关闭该对话框返回至上一级对话框。

单击“Close”按钮，关闭该级对话框。

b)定义材料特性：Main Menu: Preprocessor→Material Props→Material Models弹出对话框; 逐级双击右侧框中的Structural→Linear→Elastic →Isotropic，弹出下一级对话框。

在“弹性模量”（EX）文本框中输入“2.06e11”；在“泊松比”（PRXY）文本框中输入“0.3”；单击“OK”按钮，关闭该对话框返回至上一级对话框。

双击右侧框中的Density选项，在弹出的对话框中的“DENS”一栏中输入材料密度“7800”，单击“OK”按钮，关闭该对话框返回至上一级对话框。

关闭材料特性定义对话框。

2)三维实体模型的建立生成关键点●Main Menu: Preprocessor→Modeling →Create →Keypoints→In Active cs弹出对话框; 在Keypoint number 一栏中输入关键点编号“1”，在“X，Y，Z Location inactive cs”一栏中输入关键点1的坐标（-0.08，0，0），单击“Apply”按钮。

ansys载荷步

实际工况＝载荷步（时间步）＋载荷步（时间步）＋......载荷步＝载荷子步（时间增量）＋载荷子步（时间增量）＋......实体加载和有限元模型加载的区别：实体加载是不能利用叠加,所以实体加载要手工叠加。

对实体是覆盖,有限元模型加载是可以设置的。

有限元加载可以利用fcum进行叠加。

比如,第一个荷载步,对关键点1施加10kn,第二荷载步也对关键点1施加10kn,则这两个荷载步结果是完全一致的。

第一个荷载步,对节点1施加10kn,第二荷载步也对节点1施加10kn,而且用命令fcum,add则第二荷载步是20kn的结果。

加载与载荷步、子步及平衡迭代次数的说明加载与载荷步、子步及平衡迭代次数的说明：一、加载方式的区别实体加载和有限元模型加载的区别：实体加载是不能利用叠加,所以实体加载要手工叠加。

对实体是覆盖,有限元模型加载是可以设置的。

有限元加载可以利用fcum进行叠加。

比如,第一个荷载步,对关键点1施加10kn,第二荷载步也对关键点1施加10kn,则这两个荷载步结果是完全一致的。

第一个荷载步,对节点1施加10kn,第二荷载步也对节点1施加10kn,而且用命令fcum,add则第二荷载步是20kn的结果。

实体加载方法的优点：a、几何模型加载独立于有限元网格，重新划分网格或局部网格修改不影响载荷；b、加载的操作更加容易，尤其是在图形中直接拾取时；无论采取何种加载方式，ANSYS求解前都将载荷转化到有限元模型，因此加载到实体的载荷将自动转化到其所属的节点或单元上；二、载荷步及子步这些概念主要用于非线性分析或载荷随时间变化的问题。

根据问题的特点，可以将加载过程分为几个阶段进行，每一个阶段则作为一个载荷步。

比如做弹塑性分析时，可以通过试算初步估计开始屈服时的载荷，作为第一步，后续载荷作为第二步，等。

为了保证计算过程的收敛和结果精度 (特别是在非线性分析时)，往往把一个载荷步又划分为若干子步，每个子步施加的载荷为该子步步长和整个载荷步长之比乘以该载荷步的载荷增量值。

Ansys 荷载组合

Ansys 荷载组合1,几何模型(beam3和beam54)建立后，定义所需的element table，主要包括杆端力和最大应力，最小应力等。

然后保存数据库。

分别施加四种荷载的标准值（不乘分项系数），并分别存成四个load step file。

弹性力学ansys求解实例详解

弹性力学a n s y s求解实例详解Company number：【WTUT-WT88Y-W8BBGB-BWYTT-19998】ANSYS 上机实验报告一、题目描述如图1所示，一简支梁横截面是矩形，其面积202.0m A =，对弯曲中性轴的惯性矩451067.6m I zz -⨯=，高m h 2.0=，材料的pa E 11101.2⨯=，横向变形系数3.0=μ。

该梁的自重就是均布载荷N q 4000=和梁中点处的集中力N F 2000=，试讨论在均布荷载作用下，简支梁的最大挠度。

二、问题的材料力学解答由叠加法可知：梁上同时作用几个载荷时，可分别求出每一载荷单独作用时的变形，把各个形变叠加即为这些载荷共同作用时的变形。

在只有均布载荷q 作用时，计算简支梁的支座约束力，写出弯矩方程，利用EI M dxw d =22积分两次，最后得出：铰支座上的挠度等于零，故有0=x 时，0=w ，因为梁上的外力和边界条件都对跨度中点对称，挠曲线也应对该点对称。

因此，在跨度中点，挠曲线切线的斜率等于零，即：2l x =时，0=dx dw ，把以上两个边界条件分别代入w 和0=dxdw 的表达式，可以求出243ql C -=，0=D ，于是得转角方程及挠曲线方程为： x ql x q x ql EIw ql x q x ql EI dx dw EI 2424122464343332--=--==θ （1）在跨度中点，挠曲线切线的斜率等于零，挠度为极值，由（1）中式子可得：即EIql w q c 3845)(4-=。

在集中力F 单独作用时，查材料力学中梁在简单载荷作用下的变形表可得EIFl w F c 48)(3-=。

叠加以上结果，求得在均布载荷和集中力共同作用下，梁中点处的挠度是EIFl EI ql w w w F c q c c 483845)()(34--=+=，将各参数代入得m w c 410769.0-⨯=三、问题的ansys 解答建立几何模型此问题为可采用Beam 分析，所以该几何模型可用线表示。

用Ansys进行荷载组合的一个算例

用Ansys进行荷载组合的一个算例本算例是一个工况0的恒载和工况1的活载进行组合的算例，我们对其进行强度分析。

组合如下： 1.0*1.2恒载+0.7*1.4活载从3D3S中导出的纯模型文件另存为Model.txt从3D3S中导出的荷载文件另存为：Force_C0.txt；Force_C1.txt操作如下：（1）把Model.txt读入到Ansys中；（2）采用如下的命令流进行荷载组合分析（强度分析）；***********************************************************************************!---------------------- 删除所有荷载 ----------------------------/soluFKDELE,all,all ! 删除关键点集中荷载FDELE,all,all ! 删除节点集中荷载SFLDELE,all,all ! 删除线上面荷载SFADELE,all,all ! 删除面上面荷载SFEDELE,all,all ! 删除单元面荷载BFLDELE,all,all ! 删除线上体荷载BFADELE,all,all ! 删除面上体荷载BFVDELE,all,all! ！删除体上体荷载BFKDELE,all,all ！删除关键点上体荷载BFDELE,all,all ！删除节点体荷载!=========================== 工况分析 ===========================esel,all/input,Force_C0,txt !读入工况0 恒载/soluantype,static ！采用静力分析esel,all ！选中所有单元solve 从结果文件创建荷载工况finish/post1lcdef,1,1,0lcwrite,1finish!---------------------- 删除所有荷载 ------------------------/soluFKDELE,all,allFDELE,all,allSFLDELE,all,allSFADELE,all,allSFEDELE,all,all,allBFLDELE,all,allBFADELE,all,allBFVDELE,all,allBFKDELE,all,allBFDELE,all,all! -------------------------------------------------------------esel,all/input,Force_C1,txt !读入工况1 活载/soluantype,staticesel,all solve finish/post1lcdef,2,1,0lcwrite,2finish! -------------------------------------------------------------/post1lcfact,1,1.2 ！这里的1.2为1.0*1.2的结果lcfact,2,0.98 ！这里的0.98为1.4*0.7的结果lcase,1lcoper,add,2lcdef,3,1,0lcwrite,3eplot/REPLOT,RESIZE*********************************************************************************** （3）通过后处理命令：Main Menu>General Postproe>Load Case>Read Loadcase 读出组合后的计算结果，在本例中为Load case=3的组合结果。

有限元分析基础教程(ANSYS算例)(曾攀).

大学力学论坛搜集网络资料整理力学学习交流，资料分享请点击感谢作者的辛勤劳动！请尊重知识产权如果感觉本资料不错，请购买正版！有限元分析基础教程Fundamentals of Finite Element Analysis(ANSYS算例)曾攀清华大学2008-12有限元分析基础教程曾攀有限元分析基础教程Fundamentals of Finite Element Analysis曾攀(清华大学)内容简介全教程包括两大部分，共分9章；第一部分为有限元分析基本原理，包括第1章至第5章，内容有：绪论、有限元分析过程的概要、杆梁结构分析的有限元方法、连续体结构分析的有限元方法、有限元分析中的若干问题讨论；第二部分为有限元分析的典型应用领域，包括第6章至第9章，内容有：静力结构的有限元分析、结构振动的有限元分析、传热过程的有限元分析、弹塑性材料的有限元分析。

本书以基本变量、基本方程、求解原理、单元构建、典型例题、MATLAB程序及算例、ANSYS算例等一系列规范性方式来描述有限元分析的力学原理、程序编制以及实例应用；给出的典型实例都详细提供有完整的数学推演过程以及ANSYS实现过程。

本教程的基本理论阐述简明扼要，重点突出，实例丰富，教程中的二部分内容相互衔接，也可独立使用，适合于具有大学高年级学生程度的人员作为培训教材，也适合于不同程度的读者进行自学；对于希望在MATLAB程序以及ANSYS平台进行建模分析的读者，本教程更值得参考。

本基础教程的读者对象：机械、力学、土木、水利、航空航天等专业的工程技术人员、科研工作者。

目录[[[[[[\\\\\\【ANSYS算例】3.3.7(3) 三梁平面框架结构的有限元分析 1 【ANSYS算例】4.3.2(4) 三角形单元与矩形单元的精细网格的计算比较 3 【ANSYS算例】5.3(8) 平面问题斜支座的处理 6 【ANSYS算例】6.2(2) 受均匀载荷方形板的有限元分析9 【ANSYS算例】6.4.2(1) 8万吨模锻液压机主牌坊的分析(GUI) 15 【ANSYS算例】6.4.2(2) 8万吨模锻液压机主牌坊的参数化建模与分析(命令流) 17 【ANSYS算例】7.2(1) 汽车悬挂系统的振动模态分析(GUI) 20 【ANSYS算例】7.2(2) 汽车悬挂系统的振动模态分析(命令流) 23 【ANSYS算例】7.3(1) 带有张拉的绳索的振动模态分析(GUI) 24 【ANSYS算例】7.3(2) 带有张拉的绳索的振动模态分析(命令流) 27 【ANSYS算例】7.4(1) 机翼模型的振动模态分析(GUI) 28 【ANSYS算例】7.4(2) 机翼模型的振动模态分析(命令流) 30 【ANSYS算例】8.2(1) 2D矩形板的稳态热对流的自适应分析(GUI) 31 【ANSYS算例】8.2(2) 2D矩形板的稳态热对流的自适应分析(命令流) 33 【ANSYS算例】8.3(1) 金属材料凝固过程的瞬态传热分析(GUI) 34 【ANSYS算例】8.3(2) 金属材料凝固过程的瞬态传热分析(命令流) 38 【ANSYS算例】8.4(1) 升温条件下杆件支撑结构的热应力分析(GUI) 39 【ANSYS算例】8.4(2) 升温条件下杆件支撑结构的热应力分析(命令流) 42 【ANSYS算例】9.2(2) 三杆结构塑性卸载后的残余应力计算(命令流) 45 【ANSYS算例】9.3(1) 悬臂梁在循环加载作用下的弹塑性计算(GUI) 46 【ANSYS算例】9.3(2) 悬臂梁在循环加载作用下的弹塑性计算(命令流) 49 附录 B ANSYS软件的基本操作52 B.1 基于图形界面(GUI)的交互式操作(step by step) 53 B.2 log命令流文件的调入操作(可由GUI环境下生成log文件) 56 B.3 完全的直接命令输入方式操作56 B.4 APDL参数化编程的初步操作57【ANSYS 算例】3.3.7(3) 三梁平面框架结构的有限元分析如图3-19所示的框架结构，其顶端受均布力作用，用有限元方法分析该结构的位移。

ANSYS计算算例(结构-荷载法)解析

弹出交互式对话框，并输入直边墙底端坐标。
GUI：Main Menu> Preprocessor> Create>Line>Straight Line。
弹出如图所示拾取框，先拾取4点，再拾取5点。
（4）映射生成完整结构经过以上三步建模后，生成了半个结构，由于结构正对称，可通过映射生成完整结构。 GUI： Main Menu>Preprocessor>Reflect>Lines
计算中考虑结构自重对结构内力的影响，还应在此输入材料的密度。 Nonlinear-Density
4、模型建立（1）第一段圆弧的输入（圆心：0，－3.41，起点： 0，0，圆心角：45º ）
Menu menu：Preprocessor>Material Props>Modeling Create>Lines Arcs>By Cent & Radius
6、生成径向弹簧单元将围岩离散为弹簧单元，沿径向作用在结构的外侧。（1）局部坐标系的定义在结构顶部设置径向弹簧时，需要定义局部柱坐标系，在当前活动坐标系中定义局部坐标系的命令为： CLOCAL,KCN,KCS,XL,YL,ZL,THXY,THYZ,THZX 其中： KCN －为局部坐标号，大于10； KCS — 为坐标形式，0为直角坐标系，1为柱坐标系（Z轴为旋转轴），2为球坐标系， 3为柱坐标系（Y轴为旋转)； XL,YL,ZL—定义局部坐标系原点在当前活动坐标系中的坐标； THXY,THYZ,THZX—为相对转动角；定义后应激活该局部坐标系，才能在局部坐标系下操作，激活局部坐标系的操作为： GUI：Utility Menu>WorkPlane>Change Active CS to>Specified CoordSys 然后在交互式对话框中输入相应的局部坐标号即激活了该坐标。以后的操作要回到原先全局坐标系，也应重新激活，操作为： GUI：Utility Menu>WorkPlane>Change Active CS to> Global Cartesian／ Global Cylindrical/ Global Spherical 其中：Global Cartesian为全局直角坐标系； Global Cylindrical为全局柱坐标系； Global Spherical为全局球坐标系。

ANSYS中关于初应力和荷载步设置的算例

ANSYS中关于初应力和荷载步设置的算例[转载]ANSYS中关于初应力和荷载步设置的算例(2011-02-2016:30:41)转载▼标签：转载原文地址：ANSYS中关于初应力和荷载步设置的算例作者：WaterSprite偶初学ANSYS，看了网站上几个朋友关于应力和加载方式的讨论，自己做了一个小算例，现在发上来和大家共同讨论一下这方面的问题。

算例为一个地基+地基上面的一块方块墙吧，先通过一次计算仅加边界条件和自重，计算得到自重应力场，并输出初应力文件用来模拟初应力场。

然后施加应力和自重并进行计算，此时的位移基本为0，即消除了初位移，所得应力场即为自重应力场。

在些基础上进一步施加墙上法向面荷载，并进行第二步计算，得到的位移应该是仅有面荷载引起的位移。

关于初应力想说明的几点：1、初应力只能加在第一个荷载步，用命令流和GUI方式均可。

但在求解前不能退出求解器，如果加了初应力后退出求解器到前处理器或者后处理器后再回到前处理器，刚才施加的初应力就没有了，必须再次施加。

2、也可以不用初应力而直接分两个荷载步进行计算，如第一步仅计算自重，第二步再加面荷载后进行计算，在后处理中用工况组合来得到“净位移”，但工况组合中的应力结果似乎是不正确的。

3、如果用LSSOLVE从荷载文件进行求解，在写荷载文件时初应力的设置并不会写入荷载文件，所以，在命令流或者GUI方式下在求解前必须显示指定加载初应力。

关于荷载步设置的几点建议：1、在GUI方式下，每次进入求解器进行求解似乎都是开始一个新的分析（这一点偶也不是很明白）。

如果不退出求解器，即便不改变约束和荷载，只要求解一次，就会多一个荷载步结果，但所有结果是一样的；如果退出求解器后再进来，求解就重新开始（根据时间值）。

2、对点、线、面、体荷载都有替换和叠加两种方式，在替换方式下，在同一位置重复加荷载，只有最后一次加的荷载有效；在叠加方式下，在同一位置重复加荷载，所有荷载会叠加后共同作用在结构上。

基于ANSYS的框架结构风荷载分析

25
表 3-1 （规范）风压高度变化系数 z
离地面或海平面高度
(m)
5 10 15 20 30 40 50 60
A
1.17 1.38 1.52 1.63 1.80 1.92 2.03 2.12
地面粗糙度类别
B
C
1.00 1.00 1.14 1.25 1.42 1.56 1.67 1.77
0.74 0.74 0.74 0.84 1.00 1.13 1.25 1.35
D
0.62 0.62 0.62 0.62 0.62 0.73 0.84 0.93 26
2021/2/21
表 3-2 （线性插值）风压高度变化系数 z
2021/2/21
高度z（m） 6 9 12 15 18 21
z
0.74 0.74 0.74 0.74 0.755 0.845
高度z（m） 24 27 30 33 36 39.5
17
§2.2框架结构的模型建立
2021/2/21
整个模型采用两种单元类型:BEAM4和SHELL63。梁柱框架中的梁和立柱均采用BEAM4单元,楼板采用 SHELL63单元。
BEAM4单元菜单包括截面的高度、宽度、面积和截面惯性矩,截面形式选用矩形。
立柱:底层立柱截面采用500mm×500mm，其余各层采用450mm×450mm.
➢ 实际计算中，风压是随着高度和建筑物的外形而变化
的。一般首先选取基本风压 w
高度变化系数 z和体形系数
0 作为标准，乘以相应的，s 得出实际风压:
wzw0zs
➢ 本建筑建于长春，所以取基本风压 w0 0.6k5N/m2
➢ 高度变化系数应根据地面粗糙度类别选定，本文选取 C类。规范所注系数数值见表3-1。利用线性插值法将其计算为本文所需要的风压高度变化系数值见表3-2。

1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。

ctk
一、选定分析模块
• Main menu>preferences 选定<structural>分析模块
二、模型建立及网格划分
1、单元类型的选定 – 梁单元：beam3 – 弹簧单元：combine14 Main menu：Preprocessor>Element type>Add/Edit/Delete
（2）顶部第一段圆弧径向弹簧的设置 A、定义局部坐标在命令输入窗口中输入命令并回车： CLOCAL,11,1,0,-3.41,0,0,0,0
B、生成弹簧单元激活定义的11号局部坐标系，并按如下操作生成第一段圆弧的径向弹簧。
弹出拾取框，拾取第一段圆弧上任意一点，弹出输入弹簧参数对话框：需要输入的参数：K－弹性系数； K＝klb，k=1200MPa，l=0.4m，b=1m DX、DY、DZ—弹簧的另一端相对于选定端的偏移距离，直角坐标系中代表（dx,dy,dz），柱坐标系中代表（dr,dθ,dz），球坐标系中代表（dr,dφ,dθ）。
弹出拾取框，按<Pick All>拾取框拾取所有线条
弹出对话框
按<OK>键，映射出整个结构
生成圆弧、直线及映射过程中，点的编号又重复的现象，此时，应对点的编号进行合并和压缩。 GUI： Main Menu>Preprocessor>Numbering Ctrls>Merge Items(Compress Numbers )
6、生成径向弹簧单元将围岩离散为弹簧单元，沿径向作用在结构的外侧。（1）局部坐标系的定义在结构顶部设置径向弹簧时，需要定义局部柱坐标系，在当前活动坐标系中定义局部坐标系的命令为： CLOCAL,KCN,KCS,XL,YL,ZL,THXY,THYZ,THZX 其中： KCN －为局部坐标号，大于10； KCS — 为坐标形式，0为直角坐标系，1为柱坐标系（Z轴为旋转轴），2为球坐标系， 3为柱坐标系（Y轴为旋转)； XL,YL,ZL—定义局部坐标系原点在当前活动坐标系中的坐标； THXY,THYZ,THZX—为相对转动角；定义后应激活该局部坐标系，才能在局部坐标系下操作，激活局部坐标系的操作为： GUI：Utility Menu>WorkPlane>Change Active CS to>Specified CoordSys 然后在交互式对话框中输入相应的局部坐标号即激活了该坐标。以后的操作要回到原先全局坐标系，也应重新激活，操作为： GUI：Utility Menu>WorkPlane>Change Active CS to> Global Cartesian／ Global Cylindrical/ Global Spherical 其中：Global Cartesian为全局直角坐标系； Global Cylindrical为全局柱坐标系； Global Spherical为全局球坐标系。
B、文件导入将编辑好的命令流存为.TXT文本文件。然后用ANSYS导入。例如文件名存为load1.txt，其导入为： GUI：Utility Menu>Read Input From
2、约束的设置对于拱形结构，约束墙底的水平位移，同时为避免结构整体沉降大，使弹簧全部受拉，往往也约束墙底的竖直位移。 GUI：Menu Main> Solution>Apply>Displacement>On Nodes
列出将所有的节点坐标文件，并将其保存成另外一个文件
（3）将等效节点荷载加在结构上加载命令为：F，n，FX，num或F，n，FY，num 其中：F为加载命令， n为节点号， FX、 FY分别为为X、Y方向的力， num为力的大小。在excel 中将等效节点力计算出来，并将其整理成加载命令的表格形式，在word 中转化成命令流形式。 A、粘贴加入将word中的命令流直接粘贴到ANSYS的命令输入窗口然后回车,注意：ANSYS 需进入加载模块后方可进行粘贴，否则粘贴无效。 GUI：Menu Main>Apply
分别弹出以下对话框，在第一栏中选keypoints后确认即完成。
5、结构网格划分将实体拱形结构划分为梁单元的组合体。（1）分配单元的属性
分配单元的属性相当于把单元指向模型的某些部分，其中单元属性包括：单元类型号、材料号、实常数号、坐标号等。有两种方法把单元属性分配到模型上。 – 直接分配属性到模型图元上； – 定义默认属性集，通过分网把默认属性分配到图元上。 GUI： Main menu>Preprocessor>Define，
四、求解及弹簧的修正
1、第一次求解上述三大步完成后即可进行求解。 GUI：Menu Main> Solution>Current LS
2、查看结算结果（1）结构变形变形
在命令输入框中，先输入圆心坐标（0，-3.41）并回车，再输入起点坐标（0，0）并回车后，弹出下图交互式对话框，并输入第一段圆弧的圆心角45º 。
（2）第二段圆弧的输入（圆心：0.45，－3.86，起点： -2.411234,-0.9987659，圆心角：33.85º ）注意:第二段圆弧的起点是第一段圆弧的终点 Utility Menu>List>Keypoints>coodinates only
单元属性定义的典型对话框如图：
（2）网格划分按0.4m长度划分单元。 GUI： Menu menu>Preprocessor>Meshtool，
在右图中按<mesh>，并选取所有的线条后，即完成了网格划分。按下述操作显示单元划分结果： GUI： Utility Menu>PlotCtrls>Numbering>coodinates only
（3）直边墙的输入
第二段圆弧的终点即是直边墙的顶点，查出该点坐标为（-3.520017，-3.077515），并由此计算出直边墙底端点的坐标为（-3.520017，-7.507515）
GUI：Main Menu> Preprocessor> Create>Keypoints>In Active CS
3、定义衬砌材料属性衬砌材料为混凝土材料：弹性模量Ec＝27 GPa＝27000000000Pa 泊松比＝0.3
Menu menu：Preprocessor>Material Props>Material Model
选中右框中Stractural－Linear－Isotropic，弹出输入常数对话框，输入弹性模量和泊松比。
弹出对话框
重复选定单元
完成后,在单元类型框中显示选定的两种单元。
2、定义单元的实常数
• 梁单元实常数
• 弹簧单元实常数 Menu menu：Preprocessor>Real constants
出现下面对话框
梁单元的实常数有三个：截面积＝bh（b＝1m，h＝0.4m）惯性矩＝ bh3／12 梁高度＝ h
弹出交互式对话框，并输入直边墙底端坐标。
GUI：Main Menu> Preprocessor> Create>如图所示拾取框，先拾取4点，再拾取5点。
（4）映射生成完整结构经过以上三步建模后，生成了半个结构，由于结构正对称，可通过映射生成完整结构。 GUI： Main Menu>Preprocessor>Reflect>Lines
拾取基底两节点
选取UX,UY并输入限制位移值0
3、惯性荷载加入惯性荷载是与质量相关的荷载，该处是加入重力加速度。为了利用惯性效果来模拟重力，应当在重力的相反方向施加惯性项，即在Y方向施加一个加速度荷载相当于模拟-Y方向的重力： GUI：Menu Main> Solution>Apply>Gravity
三、加载及约束
1、荷载计算及转化 (1) 荷载计算对深埋隧道用《隧规》推荐的统计法公式计算得：竖向围岩压力：q=54084.25Pa；水平围岩压力：e＝0.15 q=8112.64Pa (2) 匀布荷载转化为等效节点荷载将匀布荷载按如下步骤转化成等效节点力。 A、列出结构上所有节点坐标 GUI：Utility Menu>List>Nodes
计算中考虑结构自重对结构内力的影响，还应在此输入材料的密度。 Nonlinear-Density
4、模型建立（1）第一段圆弧的输入（圆心：0，－3.41，起点： 0，0，圆心角：45º ）
Menu menu：Preprocessor>Material Props>Modeling Create>Lines Arcs>By Cent & Radius
ANSYS计算算例
• 拱形结构算例 – 铁路隧道 – 公路隧道 – 地铁区间隧道(矿山法施工) 矩形框架结构 – 地下街 – 地下停车场 – 地铁车站圆形结构 – 地铁区间隧道(盾构法施工)
•
•
拱形结构分析步骤
某单线电化铁路隧道衬砌如下图
a. 岩体特性 23 kN/m3 岩体为 III 级围岩，隧道埋深100m；岩体重度围岩的弹性反力系数 K 1200 MPa/m，基底围岩弹性反力系数K a 1.25K 。 b. 衬砌材料采用 C20 混凝土；重度 h 25 kN/m3，弹性模量E c 27 GPa 混凝土衬砌轴心抗压强度标准值 f ck 13.5 Mpa 1.7 MPa。 f 混凝土轴心抗拉强度标准值
按下<Apply>重复拾取其他点，并重复上述操作，完成第一段圆弧的径向弹簧设置。
（2）顶部第二、三段圆弧径向弹簧的设置首先激活全局直角坐标，再定义第二、三段圆弧的局部坐标，并激活，然后按上一步中的B步分别生成两段圆弧上径向弹簧。（3）直边墙径向弹簧的设置首先激活全局直角坐标，按（1）中B步分别生成直边墙径向弹簧。注意 DX，DY，DZ的代表意义及输入。