CLE的案例
ansys自頂而下子結(jié)構(gòu)
步驟四
/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
全模型擴(kuò)展應(yīng)力云圖：
全模型位移云圖：
最大位移與應(yīng)力：
前一篇得出的最大應(yīng)力為： 111672.715
最大變形為： 2.233489607E-06
心得：
1.在建立每個(gè)超級(jí)單元時(shí)，正確選擇所需要的基礎(chǔ)單元部分，對(duì)不需要的部分，可以刪除掉。
2.在建立使用部分時(shí)，優(yōu)先將非超級(jí)單元部分建立成功（只需要選擇對(duì)應(yīng)部分，將不需要的部分刪除就可以了）。
3.由于已經(jīng)建立 了完整的包含超級(jí)單元與非超級(jí)單元模型，因此在建立使用部分時(shí)，不需要使用setran命令，直接采用se命令建立超級(jí)單元即可。
4.由于需要后續(xù)的結(jié)果擴(kuò)展，因此在計(jì)算的過程中，慎用numcmp命令。否則，可能會(huì)導(dǎo)致結(jié)果不能擴(kuò)展成功。
展開 ansys（bottom to up）子結(jié)構(gòu)
擴(kuò)展超單元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
!獲取超級(jí)單元1最大位移
*get,max_u1,sort,0,max
nsort,s,eqv,0,0,all
!獲取超級(jí)單元1最大應(yīng)力
*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
!獲取超級(jí)單元2最大位移
*get,max_u2,sort,0,max
nsort,s,eqv,0,0,all
!獲取超級(jí)單元2最大應(yīng)力
*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
!獲取超級(jí)單元3最大位移
*get,max_u2,sort,0,max
nsort,s,eqv,0,0,all
!
展開 一個(gè)考慮腐蝕的相場(chǎng)損傷模型源代碼
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
展開 一個(gè)桿梁?jiǎn)栴}的多種解法
在不同的類型單元處，采用共節(jié)點(diǎn)的方法、
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
最終的變形和軸力云圖：
采用第二種方法，除了第二豎桿采用桿單元以外，其余全部采用采用梁?jiǎn)卧?。并在第二豎桿與橫梁處采用耦合平移自由度。
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命令。
采用第三種方法，因?yàn)榱簡(jiǎn)卧旧戆揭谱杂啥?，這里將梁?jiǎn)卧獜?qiáng)行退化成桿單元。
展開 
利用ANSYS,從純技術(shù)的角度，討論吉他（弦樂類）的調(diào)音
利用ANSYS,從純技術(shù)的角度，討論吉他（弦樂類）的調(diào)音
fini
/cle,nostart
/title,qinxian motai fenxi
!不同預(yù)拉力(通過調(diào)節(jié)某根琴弦的調(diào)音旋鈕的松緊)，相同長(zhǎng)度和直徑琴弦的不同頻率
!以下為單根琴弦
!琴弦長(zhǎng)度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
!
AL、AW、HAL在電氣圖紙中的含義，別再傻傻分不清了！
（5）燈具安裝方式的表示
CS——鏈吊
DS——管吊
W——墻壁安裝
C——吸頂
R——嵌入
S——支架
CL——柱上
SR——沿鋼線槽
BE——沿屋架或跨屋架
CLE——沿柱或跨柱
SC——穿焊接鋼管敷設(shè)
MT——穿電線管敷設(shè)
PC——穿硬塑料管敷設(shè)
FPC——穿阻燃半硬聚氯乙烯管敷設(shè)
CT——電纜橋架敷設(shè)
MR——金屬線槽敷設(shè)
PR——塑料線槽敷設(shè)
M——用鋼索敷設(shè)
KPC——穿聚氯乙烯塑料波紋電線管敷設(shè)
CP——穿金屬軟管敷設(shè)
DB——直接埋設(shè)
TC——電纜溝敷設(shè)
（6）導(dǎo)線敷設(shè)部位的標(biāo)注
AB——沿或跨梁（屋架）敷設(shè)
BC——暗敷在梁內(nèi)
AC——沿或跨柱敷設(shè)
CLC——暗敷設(shè)在柱內(nèi)
WS——沿墻面敷設(shè)
WC——暗敷設(shè)在墻內(nèi)
CE——沿天棚或頂板面敷設(shè)
CC——暗敷設(shè)在屋面或頂板內(nèi)
SCE——吊頂內(nèi)敷設(shè)
FC——地板或地面下敷設(shè)
展開 ansys之——接觸分析例題
/cle
/filname,diji11
/title,地基應(yīng)力、應(yīng)變的分布規(guī)律
/units,si
/PREP7
!ET,1,PLANE2
ET,1,PLANE42
KEYOPT,1,3,2
KEYOPT,1,5,0
KEYOPT,1,6,0
MP,EX,1,8e9 !巖石性質(zhì)
MP,PRXY,1,0.29
MP,DENS,1,2600
MP,EX,2,5.0e7 !風(fēng)化槽土體性質(zhì)?
快速看懂電路的不二法門！常用的電氣圖形符號(hào)！
弱電圖例：
給排水圖例：
4、常用字母符號(hào)大全
符號(hào)名稱和符號(hào)的解釋
（1）電線穿線管一般有
PVC——PVC管
SC——焊接鋼管
KBG——扣壓式鍍鋅薄壁電線管
JDG——緊定式鍍鋅薄壁電線管
（2）電氣設(shè)計(jì)施工圖中常用線路敷設(shè)方式
SR——沿鋼線槽敷設(shè)
BE——沿屋架或跨屋架敷設(shè)
CLE——沿柱或跨柱敷設(shè)
WE——沿墻面敷設(shè)
CE——沿天棚面或頂棚面敷設(shè)
ACE——在能進(jìn)入人的吊頂內(nèi)敷設(shè)
BC——暗敷設(shè)在梁內(nèi)
CLC——暗敷設(shè)在柱內(nèi)
WC——暗敷設(shè)在墻內(nèi)
CC——暗敷設(shè)在頂棚內(nèi)
ACC——暗敷設(shè)在不能進(jìn)入的頂棚內(nèi)
FC——暗敷設(shè)在地面內(nèi)
SCE——吊頂內(nèi)敷設(shè),要穿金屬管
（3）導(dǎo)線穿管表示
SC——焊接鋼管
MT——電線管
展開 anasys將數(shù)組中的載荷循環(huán)加載問題 急！
FINI
/CLE
/TITLE,ANALYSIS OF A ANTENNA
MODEL03
*dim,aaa,,7968,1,1
*dim,bbb,,432,1,1
*CREAT,MM
*VREAD,aaa(1),E:\ansys\model-3\temp1,txt,,7968,,,,,,
(1PE16.7)
*END/INPUT,MM
*CREAT,..
*VREAD,bbb(1),E:\ansys\model-3\temp2,txt,,432,,,,,,
(1PE16.7)
*END
/INPUT,MM
NSEL,ALL ！
展開 平面四邊形四節(jié)點(diǎn)單元計(jì)算程序與ANSYS結(jié)果對(duì)比
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,1,3,4,2
e,3,5,6,4
D,1, , , , , ,UX,UY, , , ,
D,2, , , , , ,UX, , , , ,
D,4, , , , , ,UY, , , , ,
F,5,FX,0.8e6
F,6,Fx,1e6
finish
/solu
solve
finish
/post1
PLNSOL, U,SUM, 0,1.0
PRNSOL,U,COMP
ANSYS導(dǎo)出單元?jiǎng)偠染仃嚨拿?/solu
/output,elemstiff,txt
/debug,-1,,,1
solve
/output
finish
結(jié)果對(duì)比
ANSYS的單元?jiǎng)偠染仃?MATLAB程序單元?jiǎng)偠染仃?ANSYS的位移結(jié)果
MATLAB程序的位移結(jié)果
結(jié)論
MATLAB的單元?jiǎng)偠染仃?、?jié)點(diǎn)位移計(jì)算結(jié)果與ANSYS軟件的計(jì)算結(jié)果一致，表明采用MATLAB編制的平面四邊形四節(jié)點(diǎn)單元的計(jì)算程序沒有問題。
展開 『轉(zhuǎn)貼』ansys疲勞分析一例討論
Fini
/cle
/prep7
mptemp,,,,,,,,
mptemp,1,0
et,1,plane42
mpdata,ex,1,,2e8
mpdata,prxy,,,0.3
rectng,0,200,0,100,
cyl4,100,50,25
asba,1,2
/replot
smrtsize,3
amesh,all
finish
/solu
nsel,s,loc,x,0
d,all,,,,,,ux
d,1,,,,,,uy
sfl,2,pres,0,31
allsel,all
solve
finish
/post1
plnsol,s,eqv,0,1
FP,1,100,200,500,1000,1500,2000
FP,7,10000,100000,1000000,2000000,3000000,5000000
FP,13,6000000,7e6,8e6,9e6,10e6,11e6
FP,19,12e6,15e6
FP,21,150,120,110,100,95,90
FP,27,85,80,75,70,65,60
FP,33,55,50,45,40,35,30
FP,39,29,25
*SET,n_num,node(100,75,0)
FL,1,n_num,,,,
FSNODE,n_num,1,1,
FS,n_num,1,2,1,0,0,0,0,0,0,
FE,1,100000,2,evel
FTCALC,1
PERFORM FATIGUE CALCULATION AT LOCATION 1 NODE 0
*** POST1 FATIGUE CALCULATION ***
LOCATION 1 NODE 123
EVENT/LOADS 1 1 AND 1 3
PRODUCE
展開 
APDL齒輪建模案例
參數(shù)設(shè)定
FINI
/CLE
/FILNAME,GEAR
/UNITS,SI
*AFUN,DEG
*SET,RI,0.4 !齒內(nèi)徑
*SET,RE,0.5 !齒外徑，輪內(nèi)經(jīng)
*SET,RL,0.58 !輪外徑
*SET,RH,0.04 !孔內(nèi)徑
*SET,THE,0.16 !齒厚度
*SET,THI,0.08 !輪厚度
!前處理
/PREP7
/PNUM,KP,1
/PNUM,LINE,1
ET,1,SOLID185
MP,EX,207E9
MP,NUXY,0.3
CSYS,1
!幾何建模
K,1,RI,3
K,2,RI,7
K,3,RI+0.02,8
K,4,RI+0.06,6
K,5,RE,6
K,6,RE,-6
K,7,RI+0.06,-6
K,8,RI+0.03
KGEN,30,1,8,1,,12,,8
!平面模型
*DO,I,1,240,8
L,I,I+1
SPLINE,I+1,I+2,I+3
L,I+3,I+4
L,I+4,I+5
L,I+5,I+6
SPLINE,I+6,I+7,I
AL,I,I+1,I+2,I+3,I+4,I+5,I+6,I+7
*ENDDO
!實(shí)體模型
VEXT,1,30,1,,,-THE
/VIEW,1,1,1,1
/PNUM,KP,0
/PNUM,LINE,0
!
展開 剛性基礎(chǔ)三維仿真分析(土體為DP材料)
粘聚力C＝10，內(nèi)摩擦角30，膨脹角30
?。ǎ常┩廨d：樁頂分布力Q=50,重力加速度為9.8
命令流程序?yàn)椋?/cle
/prep7
/title,zhuang,h-method
et,1,plane42
et,2,solid45
mp,ex,1,2.5e10 !混凝土材料屬性
mp,nuxy,1,0.16
mp,dens,1,2700
mp,ex,2,10e6 !土體材料屬性
mp,nuxy,2,0.45
mp,dens,2,2000
tb,dp,2 !土體本構(gòu)模型
tbdata,1,10,30,30
!建立關(guān)鍵點(diǎn)
k,1
k,2,1.5
k,3,1.5,0.4
k,4,1.1,0.4
k,5,1.1,0.8
k,6,0.7,0.8
k,7,0.7,1.2
k,8,0.3,1.2
k,9,0.3,2
k,10,0,2
k,11,4,2
k,12,4,0
k,13,4,-2
k,14,0,-2
k,15,0.3,3
k,16,0,3
!建立各條線，面元素
*do,i,1,9,1
l,i,i+1
*enddo
l,10,1
al,1,2,3,4,5,6,7,8,9,10
l,9,11
*do,i,11,13,1
l,i,i+1
*enddo
l,14,1
l,2,12
al,2,3,4,5,6,7,8,11,12,16
al,1,16,13,14,15
l,9,15
l,15,16
l,16,10
al,9,17,18,19
/pnum,area,1
/pnum,line,1
aplot
!
展開 ABAQUS流固耦合實(shí)例——Front-load washing machine
例子的來(lái)源是Abaqus CLE的官方教程，可是寫的太粗線條，我還是搞了兩天才做出了這個(gè)例子。其實(shí)就是個(gè)滾筒洗衣機(jī)帶著洗衣機(jī)里的水一起轉(zhuǎn)的問題。1. 分別為Eulerian domain和Lagrangian domain建立兩個(gè)part建立Lagrangian domain的Part，類型設(shè)置為Discrete rigid，并設(shè)置Reference Point。
建立Eulerian domain的Part，類型設(shè)置為Eulerian，要注意Eulerian domain和Lagrangian domain要保證有重疊的部分，這是一種弱耦合，數(shù)據(jù)在兩個(gè)區(qū)域間拋來(lái)拋去，所以網(wǎng)格要有重疊部分。這導(dǎo)致在Eulerian domain里有的部分是有材料的，有的地方是沒有材料的。為了之后設(shè)置材料分布時(shí)候方便，要把part實(shí)現(xiàn)劃出幾個(gè)輔助的partition。黃色虛線是在劃分partition時(shí)，為了指明Extrude/Sweep方向用到的輔助坐標(biāo)軸。
2. 定義水的材料屬性選擇狀態(tài)方程模型EOS中Us-Up，設(shè)置聲速c0=1483m/s；密度為1000kg/m3；粘度為0.001kg/ms。并把截面屬性賦給Eulerian domain。
3. 把兩個(gè)Part組裝起來(lái)
4. 新建一個(gè)Step-15. 為Eulerian domain和Lagrangian domain劃分網(wǎng)格
6. 設(shè)置接觸新建一個(gè)Contact Property，因?yàn)椴皇瞧胀ǖ拿婧兔娴慕佑|，水中的任何的一個(gè)部分可能在流動(dòng)區(qū)域里的任何一個(gè)地方和Lagrangian domain接觸，設(shè)置Tangential Behavior為Rough，賦給水和洗衣機(jī)之間的關(guān)系。
展開 關(guān)于mpc
命令流如下：
FINI
/CLE
/FILNAME,BEAM_AND_SOLID_ELEMENTS_CONNECTION !定義工作文件名
/TITLE,COUPLE_AND_CONSTRAINT_EQUATION !定義工作名
/PREP7
ET,1,SOLID95 !定義實(shí)體單元類型為SOLID95
ET,2,BEAM4 !定義梁?jiǎn)卧愋蜑锽EAM4
MP,EX,1,3E4 !定義材料的彈性模量
MP,PRXY,1,0.3 !定義泊松比
R,1 !定義實(shí)體單元實(shí)常數(shù)
R,2,10.0,10/12.0,1000/12.0,10.0,1.0 !定義梁?jiǎn)卧獙?shí)常數(shù)
BLC4,,,20,7,10 !創(chuàng)建矩形塊為實(shí)體模型
WPOFFS,0,3.5 !將工作平面向Y方向移動(dòng)3.5
WPROTA,0,90 !將工作平面繞X軸旋轉(zhuǎn)90度
VSBW,ALL !將實(shí)體沿工作平面剖開
WPOFFS,0,5 !將工作平面向Y方向移動(dòng)5
WPROTA,0,90 !將工作平面繞X軸旋轉(zhuǎn)90度
VSBW,ALL !將實(shí)體沿工作平面剖開
WPCSYS,-1 !將工作平面設(shè)為與總體笛卡兒坐標(biāo)一致
K,100,20,3.5,5 !創(chuàng)建關(guān)鍵點(diǎn)
K,101,120,3.5,5 !創(chuàng)建關(guān)鍵點(diǎn)
L,100,101 !連接關(guān)鍵點(diǎn)生成梁的線實(shí)體
LSEL,S,LOC,X,21,130 !選擇梁線
LATT,1,2,2 !指定梁的單元屬性
LESIZE,ALL,,,10 !
展開 