CLE
關(guān)注創(chuàng)建者：郁開宇 創(chuàng)建時間：2018-07-24
CLE的實例教程
步驟四
/cle
/filn,gen1
resume,full,db
/solu
alls
antype,substr
expass,on
expsol,1
seexp,se1,use
solve
save
fini
/cle
/filn,gen2
resume,full,db
/solu
alls
antype,substr
expass,on
expsol,1
seexp,se2,use
solve
save
fini
!步驟五
/cle
/filn,full
resume
/post1
file,gen2
set,last
file,gen3
set,last
file,use
set,last
alls
plnsol,s,eqv
nsort,u,sum,0,0,all
*get,max_u1,sort,0,max
nsort,s,eqv,0,0,all
*get,max_s,sort,0,max
全模型擴展應力云圖：
全模型位移云圖：
最大位移與應力：
前一篇得出的最大應力為： 111672.715
最大變形為： 2.233489607E-06
心得：
1.在建立每個超級單元時，正確選擇所需要的基礎單元部分，對不需要的部分，可以刪除掉。
2.在建立使用部分時，優(yōu)先將非超級單元部分建立成功（只需要選擇對應部分，將不需要的部分刪除就可以了）。
3.由于已經(jīng)建立 了完整的包含超級單元與非超級單元模型，因此在建立使用部分時，不需要使用setran命令，直接采用se命令建立超級單元即可。
4.由于需要后續(xù)的結(jié)果擴展，因此在計算的過程中，慎用numcmp命令。否則，可能會導致結(jié)果不能擴展成功。
展開 擴展超單元1
/cle
/filn,gen1
resume
/sol
expass,on
seexp,SE4,use
expsol,1
solve
save
fini
/post1
plnsol,s,eqv
nsort,u,sum,0,0,all
!獲取超級單元1最大位移
*get,max_u1,sort,0,max
nsort,s,eqv,0,0,all
!獲取超級單元1最大應力
*get,max_s1,sort,0,max
save
fini
/cle
/filn,gen2
resume
/sol
expass,1
seexp,SE5,use
expsol,1
solve
/post1
plnsol,s,eqv
nsort,u,sum,0,0,all
!獲取超級單元2最大位移
*get,max_u2,sort,0,max
nsort,s,eqv,0,0,all
!獲取超級單元2最大應力
*get,max_s1,sort,0,max
fini
save
/cle
/filn,gen3
resume
/sol
expass,1
seexp,SE6,use
expsol,1
solve
/post1
plnsol,s,eqv
nsort,u,sum,0,0,all
!獲取超級單元3最大位移
*get,max_u2,sort,0,max
nsort,s,eqv,0,0,all
!
展開 form and assemble stiffness matrix and internal force vector
amatrx(1:16,1:16)=amatrx(1:16,1:16)+dvol*
2 ((hphi+xkap)*matmul(matmul(transpose(b),ddsdde),b))
rhs(1:16,1)=rhs(1:16,1)-
1 dvol*(matmul(transpose(b),stress)*(hphi+xkap))
fc=AA*((cL-hphi*(cse-cle)-cle)**2)+wh*gphi
dfc=-2.d0*AA*(cL-hphi*(cse-cle)-cle)*(cse-cle)*dhphi+wh*dgphi
ddfc=wh*ddgphi-2.d0*AA*(cse-cle)*
1 (ddhphi*(cL-hphi*(cse-cle)-cle)-dhphi*dhphi*(cse-cle))
BB=matmul(transpose(dNdx),dNdx)
xm=matmul(dN,transpose(dN))
amatrx(17:24,17:24)=amatrx(17:24,17:24)+
1 +dvol*(xm/dtime+xL*xm*ddfc+xL*a_phi*BB)
rhs(17:24,1)=rhs(17:24,1)-dvol*
1 (matmul(xm,du(17:24,1))/dtime+xL*dN(:,1)*dfc
2 +xL*a_phi*matmul(transpose(dNdx),matmul(dNdx
展開 在不同的類型單元處，采用共節(jié)點的方法、
fini
/cle
/prep7
et,1,188
sectype,1,beam,rect
secdata,10,10
et,2,180
sectype,2,link
secdata,10*10
mp,ex,1,2e5
mp,prxy,1,0.3
n,1,
n,2,300
n,3,600
ngen,2,2,2,3,1,,-300
e,1,2
e,2,3
e,2,4
type,2
secn,2
e,3,5
d,4,all,,,5
f,1,fy,-100
/sol
solve
最終的變形和軸力云圖：
采用第二種方法，除了第二豎桿采用桿單元以外，其余全部采用采用梁單元。并在第二豎桿與橫梁處采用耦合平移自由度。
fini
/cle
/prep7
et,1,188
sectype,1,beam,rect
secdata,10,10
et,2,180
sectype,2,link
secdata,10*10
mp,ex,1,2e5
mp,prxy,1,0.3
n,1,
n,2,300
n,3,600
ngen,2,2,2,3,1,,-300
n,6,600
e,1,2
e,2,3
e,2,4
type,2
secn,2
e,6,5
cp,1,ux,3,6
cp,2,uy,3,6
cp,3,uz,3,6
d,4,all,,,5
f,1,fy,-100
/sol
solve
采用第二種方法得到的結(jié)果與前完全相同。本方法也可以采用相似的思路，采用耦合約束方程，即CE命令替代CP命令。
采用第三種方法，因為梁單元本身包含平移自由度，這里將梁單元強行退化成桿單元。
展開 利用ANSYS,從純技術(shù)的角度，討論吉他（弦樂類）的調(diào)音
fini
/cle,nostart
/title,qinxian motai fenxi
!不同預拉力(通過調(diào)節(jié)某根琴弦的調(diào)音旋鈕的松緊)，相同長度和直徑琴弦的不同頻率
!以下為單根琴弦
!琴弦長度0.5m
L=0.5
!琴弦半徑0.05mm
ra=0.05/1000
!π的大小
pi=acos(-1)
!琴弦截面積
area=pi*ra*ra
!材料參數(shù)
exx=1.90e11
prxyy=0.3
rou=7920
!
CLE的相關(guān)專題、標簽、搜索
CLE的最新內(nèi)容
)-cle)**2)+wh*gphi
dfc=-2.d0*AA*(cL-hphi*(cse-cle)-cle)*(cse-cle)*dhphi+wh*dgphi
ddfc=wh*ddgphi-2.d0*AA*(cse-cle)*
1 (ddhphi*(cL-hphi*(cse-cle)-cle)-dhphi*dhphi*(cse-cle))
BB=matmul
3I1jz0SVyPkKgk2ir7pGt9FayxELHyBLc2+5AdSpqNOOKAI1gKBkPi2ILjSOJWUmao5fgMCAABIAAEgAAQAAJAAAgAASAABIAAEAACQAAIAIEsAmlR4Zlr1+j4iRPEr7yvX/+6LzI0GjTGv8wvXABnuHovMiS8HdKLrnKhNbFPLGhnUg0V7fMUEe3OfN8f416FEJD9FX3q2MC/clE50QRzbJV6iGMyBYh9k9wOW92CqIBJ66vkIn71PpA8CIkigc7K
ANSYS APDL代碼
finish
/cle
/filname,hbfile
/prep7
et,1,plane182
KEYOPT,1,3,3
R,1,8,
mp,ex,1,80000.0
mp,prxy,1,0.25
n,1,0,0
n,2,0,200
n,3,200,0
n,4,200,200
n,5,400,0
n,6,400,200
type,1
mat,1
e,
無論是歐拉方法還是拉格朗日方法，在處理不同問題的時候各有優(yōu)缺點，在面對工程問題的時候，往往結(jié)合兩者的優(yōu)點，所以當前比較主流的方法是同時兼顧拉格朗日和歐拉兩者優(yōu)點的算法，比如最近幾年被大家廣泛采用的任意拉格朗日-歐拉（arbitrary Lagrangian-Eulerian,ALE）、相容拉格朗日-歐拉（coupled Lagrangian-Eulerian,CLE）等方法。
gUCEVyvX39TljtW96JjVS+a11wXN0i5Fl+RNa8dy2c5sj0cro8/oCq1G+2pV1GbehXp4a2oT72KtvBWCdhMW/8J7LNqscO3EOlTslDxZK0oEVLnYkom1gRWozSgGtXeu1Hok4n8GQV4JvyIQMfCgErkmSuxx1yBrLDDsK/6EDvDDiM47JDESzAQc5GpXFwfL4QexJbQg7ClfUMagzK3hx+CLe192FLfR7R
塑料管—PC （硬管）、 FPC（半硬管） 、 KPC （波紋管）
電纜橋架—CT 金屬軟管—CP 鋼線槽—SR
敷設部位的表示方式：
沿墻面敷設―WE 暗敷設在頂板內(nèi)―CC 沿柱或跨柱敷設―CLE
.X(李英成,劉飛,丁曉波,劉沛,羅祥勇,任亞鋒.2015.測定輕小型無人機中航測相機曝光時刻的方法及裝置.中國,CN201510178178.X)
Li Y C,Liu F,Ding X B,Yang J J,Ren Y F and Liu P.2017.Multiview stereo aerial photographic device for unmanned aerial vehi‐cles
（1）電線穿線管一般有
PVC——PVC管
SC——焊接鋼管
KBG——扣壓式鍍鋅薄壁電線管
JDG——緊定式鍍鋅薄壁電線管
（2）電氣設計施工圖中常用線路敷設方式
SR——沿鋼線槽敷設
BE——沿屋架或跨屋架敷設
CLE
四、常用字母符號大全
符號名稱和符號的解釋
（1）電線穿線管一般有
PVC——PVC管
SC——焊接鋼管
KBG——扣壓式鍍鋅薄壁電線管
JDG——緊定式鍍鋅薄壁電線管
（2）電氣設計施工圖中常用線路敷設方式
SR——沿鋼線槽敷設
BE——沿屋架或跨屋架敷設
CLE
例子的來源是Abaqus CLE的官方教程，可是寫的太粗線條，我還是搞了兩天才做出了這個例子。其實就是個滾筒洗衣機帶著洗衣機里的水一起轉(zhuǎn)的問題。1. 分別為Eulerian domain和Lagrangian domain建立兩個part建立Lagrangian domain的Part，類型設置為Discrete rigid，并設置Reference Point。
