ansys理論值計算的案例
ANSYS與ABAQUS關于梁單元后處理的計算與理論值比較（糾錯）- CAE夢想很偉大
ANSYS與ABAQUS關于梁單元后處理的計算與理論值比較（推薦）- CAE夢想很偉大
本文原創，若是轉載，請注明出處和筆名CAE-夢想很偉大。
感謝abaqus襄陽對于本文中錯誤Mises應力的問題的糾正。
本文目的
本文以工程項目中出現的評估問題為原型，以懸臂梁為例，對abaqus的mises應力在評估梁單元的如何獲得正確性進行說明。以理論計算為主，聯合ansys 和ansys workbench的計算結果，縱向評估正確的abaqus查看梁單元的正確用法beam-stress。
雖然本文可能小題大做，但是對于新手和一般不了解beam-mises的工程師，都希望引起足夠的重視。若是有任何異議，請大家留言，也歡迎大家留言討論。
具體內容如下
以10×10mm矩形截面，長度100mm的矩形管為例進行說明。
載荷：軸向載荷為10000N，彎矩為100N.m。通過理論計算
理論計算結果
軸向正應力為 ，
彎曲最大應力為
疊加組合應力
最大組合應力100+60=160
最小組合應力100-60=40
下面對比有限元計算結果與理論值比對，如表格所示
可以知道ANSYS、WB、ABAQUS顯示結果均與理論值一致。但是需要注意的是，ABAQUS需要修改截面顯示設置，需要考慮TOP和BOTTOM同時顯示數據，才能獲得正確的MISES結果。
ABAQUS的Mises不同截面激活設置顯示形式的比較如圖4所示。
展開 理論計算和CFD計算對比及不同參考值設定對阻力系數的影響-ujs
針對同一個例子，采用理論數值計算和CFD仿真計算來對比分析了二者計算的結果，并對比分析了不同湍流模型對計算結果的影響和數值理論計算的誤差，從而為以后的CFD計算提供相應的參考模型；在確定誤差較小的湍流模型的基礎上，分別設置不同的參考值來計算阻力系數，期望能夠的阻力系數以及升力系數的監測提供更進一步的支持，能夠和大家多多交流。
在這過程中感謝大家對我的幫助。
同時，該帖子也算是對http://forums.caenet.cn/showtopic-527454.aspx和http://forums.caenet.cn/showtopic-522864.aspx的解答和補充。
由于帖子內容完全由自己的體會所寫，如有錯誤的地方，請閱讀附件內容之后明確指出，
一起學習進步！
理論計算和CFD計算對比及不同參考值設定對阻力系數的影響.pdf
展開 基于ANSYS的三維CT試樣強度因子K1值計算 ￥35
幾何模型
有限元模型
等效應力
附件包括K.txt計算命令流
提供Ansys計算結果（比如加速度值）寫成fre文件的命令流
[轉貼] 提供Ansys計算結果（比如加速度值）寫成fre文件的命令流
提供Ansys計算結果（比如加速度值）寫成fre文件的命令流
將Ansys計算結果（比如加速度值）寫成fre文件的命令流。位移值更簡單。
/POST26
NSOL,2,10,U,Y,UY_2 ! 定義第二個變量為UY_2,值為10號節點Y方向的位移,節點可任選,但要保證其值非0
XVAR,1 ! 定義時間變量為坐標橫軸
PLVAR,2
*GET,num_var,VARI,0, NSETS ! 將變量長度值賦給變量num_var
k=num_var
*DIM,SYSNOISE_TITLE,CHAR,5,4
SYSNOISE_TITLE(1,1)='SYSNOISE '
SYSNOISE_TITLE(1,2)=' ACCELER '
SYSNOISE_TITLE(1,3)='ATIONS '
SYSNOISE_TITLE(1,4)=' FILE'
SYSNOISE_TITLE(2,1)='Rev 5.5 '
SYSNOISE_TITLE(2,2)=' IBM P2E '
SYSNOISE_TITLE(2,3)='SSL 11'
SYSNOISE_TITLE(2,4)='-AUG-02 '
SYSNOISE_TITLE(3,1)='ACCELERA'
SYSNOISE_TITLE(3,2)='TION_St'
SYSNOISE_TITLE(3,3)='ructure '
SYSNOISE_TITLE(4,1)='11-AUG-2002 '
SYSNOISE_TITLE(4,2)=' 10:07 '
SYSNOISE_TITLE(4,3)=':13 '
SYSNOISE_TITLE(5,1)='TIME'
*CFOPEN,ACCE,fre !
展開 
第一篇梁單元的軸力圖 （理論計算、ABAQUS仿真、ANSYS仿真方法） ￥10
第一篇梁單元的軸力圖
（理論計算、ABAQUS仿真、ANSYS仿真方法）
篇幅內容僅針對自我學習總結展示，并希望給軟件初學者帶來一定啟發。
結構有限元仿真中有兩種一維單元：桁架與梁
桁架單元：僅承受軸力作用；如二力桿。由于只在軸向承受拉/壓載荷，所以只需要定義截面面積；應力和變形均與截面形狀無關。ABAQUS 6.14-4中對應單元為truss T2D2；ANSYS 18.0中對應單元為link180。
梁單元：可承受軸向拉/壓載荷，具有承受扭轉和彎曲的能力。由于可承受扭轉、彎曲等組合變形，梁單元需要定義截面形狀。ABAQUS與ANSYS對應均為beam單元。
孫訓芳先生的《材料力學》例題2-1：一等直桿及其受力情況如下圖，試作桿的軸力圖。
由于桁架單元僅能承受拉/壓載荷；而梁單元可承受拉、壓、彎曲、扭轉的組合變形，梁單元可承受的載荷類型更為復雜，故此篇通篇采用梁單元作為分析。
展開 通過ansys利用均勻化理論計算復合材料等效性能--等效彈性模量，剪切模量等
***********計算NIX,NIY,應變矩陣元素
EPX=0
EPY=0
*DO,K,1,4 !***********計算插值點應變
EPX=NIX(K)*V(K)+EPX
EPY=NIY(K)*U(K)+EPY
*ENDDO !***********計算插值點應變
PX=PX+EPX
PY=PY+EPY
*ENDDO !***********計算單元應變PX,PY
*GET,DA,ELEM,M,AREA !提取單元的面積
K=EPSI/(1+NU)/2/AREA
E1212H=K*(DA-PY*KA-PX*KA
)+E1212H
*GET,M,ELEM,M,NXTH
*ENDDO
*enddo
allsel,
finish
G12=E1212H
展開 