ANSYS计算算例(结构-荷载法)
- 格式:ppt
- 大小:3.40 MB
- 文档页数:46
选取新建铁路宜昌(宜)-万州(万)铁路线上的别岩槽隧道某断面.该断面设计单位采用的支护结构如图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计算结果与分析一、有限元原理:有限元的解题思路可简述为:从结构的位移出发,通过寻找位移和应变,应变与应力,应力与内力,内力与外力的关系,建立相应的方程组,从而由已知的外力求出结构的内应力和位移。
有限元分析过程由其基本代数方程组成:[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 荷载组合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 motion 载荷计算ANSYS Motion是一种多体动力学仿真软件,用于模拟和分析机械系统的运动和载荷。
通过ANSYS Motion,用户可以预测和优化产品的性能、可靠性和耐久性,从而降低开发成本和时间。
本文将介绍ANSYS Motion在载荷计算方面的应用。
一、概述在机械系统设计和分析过程中,载荷计算是非常重要的一环。
它可以帮助工程师了解系统在实际工作条件下所承受的力和力矩,从而评估系统的可靠性和耐久性。
ANSYS Motion提供了一系列强大的工具和功能,用于准确计算和分析机械系统的载荷。
二、载荷类型在ANSYS Motion中,可以对机械系统施加多种类型的载荷。
常见的载荷类型包括静态载荷、动态载荷和惯性载荷。
静态载荷是指在系统静止时施加的力和力矩,动态载荷是指在系统运动时施加的力和力矩,惯性载荷是指由于系统加速度而产生的力和力矩。
三、载荷计算方法ANSYS Motion提供了多种载荷计算方法,包括直接计算、求解拉格朗日方程和使用约束等。
直接计算方法适用于简单的系统,通过施加已知的力和力矩来计算载荷。
拉格朗日方程方法适用于复杂的系统,通过求解拉格朗日方程来计算载荷。
约束方法适用于特殊的系统,通过施加约束来计算载荷。
四、载荷分析通过ANSYS Motion的载荷计算功能,用户可以对机械系统进行载荷分析。
载荷分析可以帮助工程师了解系统在不同工况下的载荷分布、载荷大小以及载荷的变化规律。
通过对不同工况下的载荷进行分析,用户可以评估系统的可靠性和耐久性,并优化系统的设计。
五、应用案例ANSYS Motion广泛应用于各个行业的机械系统设计和分析。
例如,在汽车行业中,可以使用ANSYS Motion对汽车悬挂系统进行载荷计算,评估悬挂系统在不同路况下的载荷分布和变化规律,从而优化悬挂系统的设计。
在航空航天行业中,可以使用ANSYS Motion 对飞机起落架进行载荷计算,评估起落架在起飞和降落过程中的载荷大小和变化规律,从而提高起落架的可靠性和耐久性。
钢架杆件横截面面积均为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。
实例一: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软件进行有限元计算实例工字钢梁结构静力分析一工字钢梁两端均为固定端,其截面尺寸为: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”按钮。
实际工况=载荷步(时间步)+载荷步(时间步)+......载荷步=载荷子步(时间增量)+载荷子步(时间增量)+......实体加载和有限元模型加载的区别:实体加载是不能利用叠加,所以实体加载要手工叠加。
对实体是覆盖,有限元模型加载是可以设置的。
有限元加载可以利用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 荷载组合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命令指定的荷载工况。
弹性力学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进行荷载组合的一个算例本算例是一个工况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的组合结果。
大学力学论坛搜集网络资料整理力学学习交流,资料分享请点击感谢作者的辛勤劳动!请尊重知识产权如果感觉本资料不错,请购买正版!有限元分析基础教程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中关于初应力和荷载步设置的算例[转载]ANSYS中关于初应力和荷载步设置的算例(2011-02-2016:30:41)转载▼标签:转载原文地址:ANSYS中关于初应力和荷载步设置的算例作者:WaterSprite偶初学ANSYS,看了网站上几个朋友关于应力和加载方式的讨论,自己做了一个小算例,现在发上来和大家共同讨论一下这方面的问题。
算例为一个地基+地基上面的一块方块墙吧,先通过一次计算仅加边界条件和自重,计算得到自重应力场,并输出初应力文件用来模拟初应力场。
然后施加应力和自重并进行计算,此时的位移基本为0,即消除了初位移,所得应力场即为自重应力场。
在些基础上进一步施加墙上法向面荷载,并进行第二步计算,得到的位移应该是仅有面荷载引起的位移。
关于初应力想说明的几点:1、初应力只能加在第一个荷载步,用命令流和GUI方式均可。
但在求解前不能退出求解器,如果加了初应力后退出求解器到前处理器或者后处理器后再回到前处理器,刚才施加的初应力就没有了,必须再次施加。
2、也可以不用初应力而直接分两个荷载步进行计算,如第一步仅计算自重,第二步再加面荷载后进行计算,在后处理中用工况组合来得到“净位移”,但工况组合中的应力结果似乎是不正确的。
3、如果用LSSOLVE从荷载文件进行求解,在写荷载文件时初应力的设置并不会写入荷载文件,所以,在命令流或者GUI方式下在求解前必须显示指定加载初应力。
关于荷载步设置的几点建议:1、在GUI方式下,每次进入求解器进行求解似乎都是开始一个新的分析(这一点偶也不是很明白)。
如果不退出求解器,即便不改变约束和荷载,只要求解一次,就会多一个荷载步结果,但所有结果是一样的;如果退出求解器后再进来,求解就重新开始(根据时间值)。
2、对点、线、面、体荷载都有替换和叠加两种方式,在替换方式下,在同一位置重复加荷载,只有最后一次加的荷载有效;在叠加方式下,在同一位置重复加荷载,所有荷载会叠加后共同作用在结构上。