[求助]大家好，斜拉桥的模型怎么建啊？

alac · 发表于 2006-3-28 21:50

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。

您需要登录才可以下载或查看，没有账号？我要加入

x

斜拉桥模型怎么建啊？ 桥身、桥墩、柱塔和拉索应该处理成什么单元，拉索和桥身怎么连接的啊？

xj2070 · 发表于 2006-3-28 22:06

桥身、桥墩、柱塔可以用梁单元，壳单元；拉索用link 10

alac · 发表于 2006-3-28 22:15

请大家指教啊 谢谢！

alac · 发表于 2006-3-28 22:17

<DIV class=quote>以下是引用xj2070在2006-3-28 22:06:49的发言： 桥身、桥墩、柱塔可以用梁单元，壳单元；拉索用link 10 </DIV>
 谢谢你！ 那节点是怎么连接的啊，刚节点还是铰接，那么多节点怎么做啊？拉索还要考虑非线性因素呢

VibInfo · 发表于 2006-3-28 22:27

一个例子 !******************2003.10.27-13:00********************* !*************以shell单元模拟横梁-3-恒载**************** !******************配跨为曲线************************** /prep7 !*******************定义节点********************** !*********************塔************************** !*******************塔上部************************ *do,i,54.5,33.5,-0.5 k,110-i*2,-(82.1495+i)/7.403322,i,0 *enddo !*******************塔下部************************ !*远朝 !*上 *do,i,33,0,-0.5 k,110-i*2,-15.6213*i/33.5,i,7.4839*i/33.5-7.4839 *enddo !*下 *do,i,-0.5,-2.5,-0.5 k,110-i*2,-2.7978*i/6,i,-i/2-7.4839 *enddo !*近朝 !*上 *do,i,33,0,-0.5 k,182-i*2,-15.6213*i/33.5,i,-7.4839*i/33.5+7.4839 *enddo !*下 *do,i,-0.5,-2.5,-0.5 k,182-i*2,-2.7978*i/6,i,i/2+7.4839 *enddo !*********桥面(配跨为曲线)****** k,192,-2,0,-6 k,193,0,0,-6 k,194,2,0,-6 k,195,9,0,-6 *do,i,196,205 k,i,(i-195)*9+9,0,-6 *enddo k,206,105,0,-6 k,207,-2,0,6 k,208,0,0,6 k,209,2,0,6 k,210,9,0,6 *do,i,211,220 k,i,(i-210)*9+9,0,6 *enddo k,221,105,0,6 !****************背索锚碇******************* k,222,-31,0,14 k,223,-31,0,-14 !***********************定义单元********************** !上部塔的截面几何特性 *dim,BBS,,43 *dim,HHS,,43 *dim,SSS,,43 *dim,IIYYS,,43 *dim,IIZZS,,43 *do,i,1,43 BBS(i)=1.5*(i-1)/42+3 HHS(i)=(i-1)/21+3 SSS(i)=BBS(i)*HHS(i) IIzzS(i)=BBS(i)*HHS(i)**3/12 IIyyS(i)=HHS(i)*BBS(i)**3/12 *enddo !**********平均********** *dim,BS,,42 *dim,HS,,42 *dim,SS,,42 *dim,IYYS,,42 *dim,IZZS,,42 *do,i,1,42 hS(i)=(hhS(i)+hhS(i+1))/2 bS(i)=(bbS(i)+bbS(i+1))/2 SS(i)=(SSS(i)+SSS(i+1))/2 IYYS(i)=(IIYYS(i)+IIYYS(i+1))/2 IZZS(i)=(IIZZS(i)+IIZZS(i+1))/2 *enddo !下部塔的截面几何特性 *dim,BBX,,80 *dim,HHX,,80 *dim,SSX,,80 *dim,IIYYX,,80 *dim,IIZZX,,80 *do,i,1,80 BBX(i)=2.5 HHX(i)=(i-1)*4.7164/79+5 SSX(i)=BBX(i)*HHX(i) IIzzX(i)=BBX(i)*HHX(i)**3/12 IIyyX(i)=HHX(i)*BBX(i)**3/12 *enddo !**********平均********** *dim,BX,,79 *dim,HX,,79 *dim,SX,,79 *dim,IYYX,,79 *dim,IZZX,,79 *do,i,1,79 BX(i)=(BBX(i)+BBX(i+1))/2 HX(i)=(HHX(i)+HHX(i+1))/2 SX(i)=(SSX(i)+SSX(i+1))/2 IYYX(i)=(IIYYX(i)+IIYYX(i+1))/2 IZZX(i)=(IIZZX(i)+IIZZX(i+1))/2 *enddo !******定义单元类型-1-42---塔上部****** *do,i,1,42 et,i,beam4 mp,ex,i,3.4e10 MP,DENS,i,2.5e3 mp,prxy,i,0.3 r,i,SS(i),IZZS(i),IYYS(i),bS(i),hS(i) *enddo !******定义单元类型-43-115---塔下部****** *do,i,43,114 et,i,beam4 mp,ex,i,3.4e10 MP,DENS,i,2.5e3 mp,prxy,i,0.3 r,i,SX(i-42),IZZX(i-42),IYYX(i-42),bX(i-42),hX(i-42) *enddo !******定义单元类型-115---桥面-1****** et,115,shell63 MP,EX,115,3.4e10 MP,DENS,115,2.5e3 mp,prxy,115,0.3 r,115,1,,,,, !******定义单元类型-116---桥面-2****** et,116,shell63 MP,EX,116,3.4e10 MP,DENS,116,2.5e3 mp,prxy,116,0.3 r,116,4,,,,, !******定义单元类型-117---桥面-3****** et,117,shell63 MP,EX,117,3.4e10 MP,DENS,117,2.5e3 mp,prxy,117,0.3 r,117,1,,,,, !******定义单元类型-118-119---背索****** *dim,cable_strain_bs,,2 cable_strain_bs(1)=3.978873577e-3 !2500 cable_strain_bs(2)=3.978873577e-3 !2500 et,118,link10 MP,EX,118,2e11 MP,DENS,118,7.85e3 mp,prxy,118,0.3 R,118,0.031416,cable_strain_bs(1) et,119,link10 MP,EX,119,2e11 MP,DENS,119,7.85e3 mp,prxy,119,0.3 R,119,0.031416,cable_strain_bs(2) !******定义单元类型-120-129---主索****** *dim,cable_strain_zs,,11 cable_strain_zs(1)=2.419155135e-3 !380 cable_strain_zs(2)=2.419155135e-3 !380 cable_strain_zs(3)=2.2918311805e-3 !360 cable_strain_zs(4)=2.2918311805e-3 !360 cable_strain_zs(5)=2.228169203e-3 !350 cable_strain_zs(6)=2.228169203e-3 !350 cable_strain_zs(7)=1.909859317e-3 !300 cable_strain_zs(8)=1.65521141e-3 !260 cable_strain_zs(9)=1.65521141e-3 !260 cable_strain_zs(10)=1.40056340e-3 !220 cable_strain_zs(11)=1.273239545e-3 !200 *do,i,120,130 et,i,link10 MP,EX,i,2e11 MP,DENS,i,7.85e3 mp,prxy,i,0.3 R,i,0.007853982,cable_strain_zs(i-119) !定义索面积直径=0.1米，初始应变 *enddo !********定义桥面与塔的连接*********** et,131,beam4 mp,ex,131,3.4e10 MP,DENS,131,2.5e3 mp,prxy,131,0.3 r,131,16,21.3333,21.3333,4,4 !*************************定义塔上部-beam4-1************** *do,i,1,42 l,i,i+1 lsel,u,type,,1,131 latt,i,i,i,,,, lesize,all,,,1,,1 lmesh,all,all alls *enddo !*************************塔下部-beam4-2*************** !*************************近朝************************* *do,i,43,114 l,i,i+1 lsel,u,type,,1,131 latt,i,i,i,,,, lesize,all,,,1,,1 lmesh,all,all alls *enddo !*************************远朝************************* l,43,116 lsel,u,type,,1,131 latt,43,43,43,,,, lesize,all,,,1,,1 lmesh,all,all alls *do,i,116,186 l,i,i+1 lsel,u,type,,1,131 latt,i-72,i-72,i-72,,,, lesize,all,,,1,,1 lmesh,all,all alls *enddo !*******************定义索-link10******************** !*******************背索***************************** l,9,222 lsel,u,type,,1,131 latt,118,118,118,,,, lesize,all,,,1,,1 lmesh,all,all alls l,9,223 lsel,u,type,,1,131 latt,119,119,119,,,, lesize,all,,,1,,1 lmesh,all,all alls !*******************主索***************************** *do,i,1,11 l,i*4+1,206-i l,i*4+1,221-i lsel,u,type,,1,131 latt,i+119,i+119,i+119,,,, lesize,all,,,1,,1 lmesh,all,all alls *enddo !****************定义桥面与塔的连接************* l,110,193 l,182,208 lsel,u,type,,1,131 latt,131,131,131,,,, lesize,all,,,1,,1 lmesh,all,all alls !******************定义桥面-area-4******************** wproat,0,-90,0 cyl4,-5.6695,-70,76,90,64,118.09989655,0 wproat,0,90,0 a,188,192,207,191 lsel,u,type,,1,131 aatt,115,115,115,,,, !lesize,all,2,,,,2 lesize,213,2,,,,2 lesize,214,2,,,,2 lesize,215,2,,,,2 lesize,216,2,,,,2 lesize,217,,,,4,2 lesize,218,2,,,,2 lesize,219,,,,4,2 alls amesh,all a,192,193,208,207 a,193,194,209,208 lsel,u,type,,1,131 aatt,116,116,116,,,, lesize,all,2,,,,1 alls amesh,all *do,i,1,12 a,i+193,i+194,i+209,i+208 *enddo lsel,u,type,,1,131 aatt,117,117,117,,,, lesize,all,2,,,,1 alls amesh,all !*****************活载**************** !*****不计冲击力-验算最大竖向挠度***** !sfa,all,1,PRES,1834 !f,196,fy,-114286 !f,197,fy,-114286 !*************************************** f,196,fy,-160000 f,197,fy,-160000 finish !***********************加约束************************* /SOLU d,115,all d,187,all d,188,all d,189,all d,235,uy d,239,uy d,280,uy d,292,uy d,735,uy d,739,uy !***********重力加速度********** acel,,12 sfa,all,1,PRES,2567 !acel,,10 solve finish /POST1 !进入后处理模块POST1 esel,s,type,,118,130 etable,link-fx,smisc,1 esel,s,type,,1,114 etable,MXI,smisc,4 esel,s,type,,1,114 etable,MXJ,smisc,10 esel,s,type,,1,114 etable,MYI,smisc,5 esel,s,type,,1,114 etable,MYJ,smisc,11 esel,s,type,,1,114 etable,MZI,smisc,6 esel,s,type,,1,114 etable,MZJ,smisc,12 esel,s,type,,all

alac · 发表于 2006-3-28 22:31

谢谢各位！ 可是楼上的那些程序是什么？是APDL吗？参数化设计语言？ 那我应该把那些程序从哪儿输入呢？看什么才会编这些程序呢 谢谢你！

alac · 发表于 2006-3-28 22:37

请高手师傅接着指教哦 我应该怎么用那些程序呢？以前接触的都是简单的机械建模，没接触过那么多单元，也不会使用参数化设计语言啊。请您接着叫我哦，很感激的！

VibInfo · 发表于 2006-3-28 22:41

ansys 有一个命令输入条，在那里输入就行了 要学习apdl先从整理log文件入手吧 命令流是学习ansys的必经之路，可以很好的理解和掌握，在以后的计算中，可以很方便的进行修改数据，改变参数。

alac · 发表于 2006-3-28 23:28

log文件在哪啊？ 命令流去那个地方看呢，保存路径里面的记事本文件内容和辅导书上的命令流一样的地方，但不全一样而且不完全。怎么看命令流啊？

VibInfo · 发表于 2006-3-29 00:01

每运行一次ANSYS生成一段log文件，其命名为Jobname．1og，Jobname是运行ANSYS程序时给定的作业名。log文件在首次运行ANSYS时打开，退出程序时关闭，ANSYS运行中输入的每一个命令都拷贝到log文件，形成ANSYS的完整记录。 l1.要注意时间，因为每次做的东西都会跟在log文件后面，所以要根据时间取舍，不是所有的log文件中的内容就有用的，一开始我建议从新建一个文件开始。 2.最好每做一步看一下log文件，可以知道自己的操作对应哪些命令 3.有些关于存盘、显示视角等命令可以删除。 4.选取实体时往往会产生很多命令，可以简化。 5.整理命令流时要新建立一个文本文件，以便从log文件中拷贝所需要的。 6.File菜单中的Read input from可以读入自己所建立的命令流来执行。 7.可以增加注释语句以增强可读性 下面以一简单模型为例大致说明一下： /BATCH /COM,ANSYS RELEASE 5.7.1 UP20010418 15:42:42 09/12/2003 /input,menust,tmp ,'',,,,,,,,,,,,,,,,1 /GRA,POWER /GST,ON /PLO,INFO,3 /COL,PBAK,ON,1,BLUE 这一段基本上没有用，是开始的设置，每个分析都是差不多这样，因此可以删除。 /PREP7 进入前处理器。在这一句前面可以加上fini /clear，这样可以把模型原来的内容清空。 !* ET,1,PLANE42 !* 定义单元 MPTEMP,,,,,,,, MPTEMP,1,0 MPDATA,EX,1,,2e11 MPDATA,PRXY,1,,0.33 MPTEMP,,,,,,,, MPTEMP,1,0 MPDATA,DENS,1,,2700 定义材料类型 K,1,0,0,, K,2,50,0,, K,3,50,10,, K,4,10,10,, K,5,10,50,, K,6,0,50,, 建立关键点 FLST,2,6,3 FITEM,2,1 FITEM,2,2 FITEM,2,3 FITEM,2,4 FITEM,2,5 FITEM,2,6 A,P51X 将几个关键点连接成面 这里，关键点是通过鼠标选取而得到，因此命令较多，其实这一段可以改为，A,1,2,3,4,5,6,具体如何改写可以参考FLST,FITEM命令的帮助 ESIZE,1,0, 设置单元大小 CM,_Y,AREA ASEL, , , , 1 CM,_Y1,AREA CHKMSH,'AREA' CMSEL,S,_Y !* AMESH,_Y1 EPLOT 这一段的含义是对所选择的面1进行网格划分，可以改写成AMESH,1。 FINISH 退出前处理器 /SOLU 进入求解器 FLST,2,1,4,ORDE,1 FITEM,2,5 !* /GO DL,P51X, ,ALL, 给选择的线加位移约束。 可以改写成：DL,5,,ALL FLST,2,2,4,ORDE,2 FITEM,2,2 FITEM,2,-3 /GO !* SFL,P51X,PRES,1000, 给选择的线加压力 可以改写成SFL,3,PRES,1000 SFL,2,PRES,1000 SOLVE 求解 FINISH 退出求解器 /POST1 进入后处理器 SET,FIRST 读出第一个解集 PLNSOL,U,SUM,0,1 画出节点位移解。 因此整理后的命令流文件如下： fini /clear /PREP7 !定义单元 ET,1,PLANE42 !定义材料 MPTEMP,,,,,,,, MPTEMP,1,0 MPDATA,EX,1,,2e11 MPDATA,PRXY,1,,0.33 MPTEMP,,,,,,,, MPTEMP,1,0 MPDATA,DENS,1,,2700 !建立模型 K,1,0,0,, K,2,50,0,, K,3,50,10,, K,4,10,10,, K,5,10,50,, K,6,0,50,, A,1,2,3,4,5,6, !划分网格 ESIZE,1,0, AMESH,1 FINISH /SOLU !添加约束 DL,5,,ALL SFL,3,PRES,1000 SFL,2,PRES,1000 SOLVE FINISH /POST1 SET,FIRST PLNSOL,U,SUM,0,1

alac · 发表于 2006-3-29 01:00

VibInfo院长好厉害了，太谢谢你了！受益匪浅啊。 生成的1og文件是不是以我刚开始命名的Jobname记事本文件啊，我觉得应该是，没有显示是log文件的啊。为了保存各种图形比如网格划分、位移等值线等，我应该相应的改变工作标题吧？ 非常感谢大家的帮助，谢谢！

alac · 发表于 2006-3-29 23:44

桥身、拉索分别用beam和link单元，然后后面的选择呢，是2D还是3D?具体应该怎么选择？ 拉索节点怎么确定呢，就是如何确定拉索两端位置，节点连接方式是固结吗？ 那桥身和柱塔如果是悬浮的怎么处理啊？

alac · 发表于 2006-3-31 01:58

院长接着给指点哦，谢谢了！

glise · 发表于 2006-3-31 06:10

<DIV class=quote>以下是引用alac在2006-3-29 23:44:07的发言： 桥身、拉索分别用beam和link单元，然后后面的选择呢，是2D还是3D?具体应该怎么选择？ 拉索节点怎么确定呢，就是如何确定拉索两端位置，节点连接方式是固结吗？ 那桥身和柱塔如果是悬浮的怎么处理啊？</DIV>
 可以参考5楼给的命令流，过程给的还是挺清晰的，注解也很清楚

alac · 发表于 2006-3-31 10:40

可是我不知道APDL是怎么用啊？ 从那输入那些命令都不知道，~~~~没用过啊？

账号		自动登录	找回密码
密码			我要加入

[结构分析] [求助]大家好，斜拉桥的模型怎么建啊？

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。

re:斜拉桥的模型怎么建啊

回复：（xj2070）re:斜拉桥的模型怎么建啊

回复：（alac）[求助]大家好，斜拉桥的模型怎么建啊？...

回复：（alac）[求助]大家好，斜拉桥的模型怎么建啊？...

回复：（alac）[求助]大家好，斜拉桥的模型怎么建啊？...

回复：（alac）桥身、拉索分别用beam和link单元，然后...

浏览过的版块