【原创】吐血奉献：铁路轨道有限元分析

binwang_2003

laoxiong请问怎样按照自己想要的格式写到文件中？不要卖关子了，呵呵。

yayoho

不错，谢谢了

sl3913

ding!

cuihunter

对于楼主的奉献精神真是佩服得五体投地，例子更是精辟，值得好好研究！！！

evern.yang

非常好的一个例子，我好像在ANSYS论文集中看过你的这篇论文，通过你，我看到了中国仿真业的美好的前程。

duderi

ding

sl3913

完了吗！期待下回分解！

77152587

的确是一个好例子，不过我准备按照你的步骤去做时，有一些问题不明白，向你请教：
1、网格划分控制设置问题：首先将L1、L15一分为二，再按照你提供的数据从上到下进行线等分份数来设置，可是提示A1、A2不规则不能划分；
2、怎样设置线的关联？
3、后处理中为了得到加载点的真实位移需要对纵向上选取结点位移值来进行曲线拟合，拟合方法？
  
期待你的答复。谢谢！

yujb025

不知yiby兄台可不可以稍加整理成稿件投到simwe  期刊。可发邮件给我。

yujb025

在静力分析范围内求解等效刚度，不必求解那么多次。如果非要这样，也可以用函数加载的方法计算。楼主无私奉献，精神可嘉！
*dim,forcez,table,11,1
*dim,forcey,table,11,1
*do,i,1,11,1
forcez(i,0)=(i-1)*5
forcez(i,1)=(i-1)*5
forcey(i,0)=(i-1)*5
forcey(i,1)=-105
*enddo
f,400,fz,%forcez%
F,400,FY,%forcey%
time,50
nsubst,10
outres,all,all
solve
/post26
nsol,2,400,u,z
plvar,2

dzzufo50

非常感谢，让我对中国的未来事业很有信心了！
支持！

yiby

FINISH
/CLEAR,NOSTART
RESUME,rail_3Dmesh_withoutLoad,db
!/VIEW,1,1,1,1   
!GPLOT
/SOLU
F,ARG1,FY,ARG2
F,ARG1,FZ,ARG3
SOLVE
FINISH
  
/POST1
*DIM,DATAP1,,31,2
*DIM,P1,,31,1
P1(1,1)=7
*GET,DATAP1(1,1),NODE,P1(1,1),LOC,X
*GET,DATAP1(1,2),NODE,P1(1,1),U,Z
*DO,IP,2,30
P1(IP,1)=IP+384
*GET,DATAP1(IP,1),NODE,P1(IP,1),LOC,X
*GET,DATAP1(IP,2),NODE,P1(IP,1),U,Z
*ENDDO
P1(31,1)=385
*GET,DATAP1(31,1),NODE,P1(31,1),LOC,X
*GET,DATAP1(31,2),NODE,P1(31,1),U,Z
  
*DIM,DATAP2,,31,2
*DIM,P2,,31,1
P2(1,1)=9
*GET,DATAP2(1,1),NODE,P2(1,1),LOC,X
*GET,DATAP2(1,2),NODE,P2(1,1),U,Z
*DO,IP,2,30
P2(IP,1)=IP+842
*GET,DATAP2(IP,1),NODE,P2(IP,1),LOC,X
*GET,DATAP2(IP,2),NODE,P2(IP,1),U,Z
*ENDDO
P2(31,1)=835
*GET,DATAP2(31,1),NODE,P2(31,1),LOC,X
*GET,DATAP2(31,2),NODE,P2(31,1),U,Z
  
*DIM,DATAP3,,31,2
*DIM,P3,,31,1
P3(1,1)=10
*GET,DATAP3(1,1),NODE,P3(1,1),LOC,X
*GET,DATAP3(1,2),NODE,P3(1,1),U,Z
*DO,IP,2,30
P3(IP,1)=IP+871
*GET,DATAP3(IP,1),NODE,P3(IP,1),LOC,X
*GET,DATAP3(IP,2),NODE,P3(IP,1),U,Z
*ENDDO
P3(31,1)=836
*GET,DATAP3(31,1),NODE,P3(31,1),LOC,X
*GET,DATAP3(31,2),NODE,P3(31,1),U,Z
  
*DIM,DATAP4,,31,2
*DIM,P4,,31,1
P4(1,1)=11
*GET,DATAP4(1,1),NODE,P4(1,1),LOC,X
*GET,DATAP4(1,2),NODE,P4(1,1),U,Z
*DO,IP,2,30
P4(IP,1)=IP+900
*GET,DATAP4(IP,1),NODE,P4(IP,1),LOC,X
*GET,DATAP4(IP,2),NODE,P4(IP,1),U,Z
*ENDDO
P4(31,1)=837
*GET,DATAP4(31,1),NODE,P4(31,1),LOC,X
*GET,DATAP4(31,2),NODE,P4(31,1),U,Z
  
*DIM,DATAP5,,31,2
*DIM,P5,,31,1
P5(1,1)=12
*GET,DATAP5(1,1),NODE,P5(1,1),LOC,X
*GET,DATAP5(1,2),NODE,P5(1,1),U,Z
*DO,IP,2,30
P5(IP,1)=IP+929
*GET,DATAP5(IP,1),NODE,P5(IP,1),LOC,X
*GET,DATAP5(IP,2),NODE,P5(IP,1),U,Z
*ENDDO
P5(31,1)=838
*GET,DATAP5(31,1),NODE,P5(31,1),LOC,X
*GET,DATAP5(31,2),NODE,P5(31,1),U,Z
  
*DIM,DATAP6,,31,2
*DIM,P6,,31,1
P6(1,1)=13
*GET,DATAP6(1,1),NODE,P6(1,1),LOC,X
*GET,DATAP6(1,2),NODE,P6(1,1),U,Z
*DO,IP,2,30
P6(IP,1)=IP+958
*GET,DATAP6(IP,1),NODE,P6(IP,1),LOC,X
*GET,DATAP6(IP,2),NODE,P6(IP,1),U,Z
*ENDDO
P6(31,1)=839
*GET,DATAP6(31,1),NODE,P6(31,1),LOC,X
*GET,DATAP6(31,2),NODE,P6(31,1),U,Z
  
*CFOPEN,ARG4,TXT
*VWRITE,DATAP1(1,1),DATAP1(1,2)
(31F16.9)   
*CFOPEN,ARG5,TXT
*VWRITE,DATAP2(1,1),DATAP2(1,2)
(31F16.9)   
*CFOPEN,ARG6,TXT
*VWRITE,DATAP3(1,1),DATAP3(1,2)
(31F16.9)
*CFOPEN,ARG7,TXT
*VWRITE,DATAP4(1,1),DATAP4(1,2)
(31F16.9)
*CFOPEN,ARG8,TXT
*VWRITE,DATAP5(1,1),DATAP5(1,2)
(31F16.9)
*CFOPEN,ARG9,TXT
*VWRITE,DATAP6(1,1),DATAP6(1,2)
(31F16.9)
FINISH

yiby

后期处理：
http://www.simwe.com/cgi-bin/ut/topic_show.cgi?id=52368&h=1&bpg=2&age=0

aboluocn

哎，看了大半天，还是不是很懂，看来还得加紧学习啊

ydyz

你怎么没了啊，还没有完啊，还期待着呢