ansys中的apdl的案例
ANSYS經(jīng)典中使用APDL語言施加移動(dòng)高斯熱源
很多人在使用ANSYS模擬焊接和增材制造過程中都面臨高斯熱源施加的難題，現(xiàn)在我來演示一下如何在ANSYS經(jīng)典中使用APDL語言施加高斯熱源，以及如何實(shí)現(xiàn)熱源的移動(dòng)。
打開經(jīng)典界面，然后選擇Parameters→Functions→Define/Edit
然后在彈出的Function Editor中選取你想要輸入的熱源函數(shù)，我這里使用了一個(gè)高斯體熱源函數(shù)，也可以替換成高斯面熱源或者雙橢球熱源，具體函數(shù)請自行查找文獻(xiàn)
點(diǎn)擊Save后，保存后綴名為.func的函數(shù)文件，其名稱為func11.func
然后退出，重新選擇Parameters→Functions→Read From file
選擇剛才定義的函數(shù)
此時(shí)彈出對話框，要求輸入函數(shù)的名稱，及對應(yīng)的參數(shù)的大小，我們定義名稱為gauss，兩個(gè)參數(shù)常量分別為qmx=1，r=1，局部坐標(biāo)系選0就意味著這個(gè)函數(shù)是在全局坐標(biāo)系中施加的，可以換成其他已經(jīng)定義的任何局部坐標(biāo)系
然后點(diǎn)擊List→Files→Log file
然后我們就可以發(fā)現(xiàn)在Log file文件里自動(dòng)生成了函數(shù)func11對應(yīng)的數(shù)據(jù)表，其是一個(gè)維度6*20的Table表，我們在array parameter中也可以查看其具體數(shù)據(jù)，為什么會(huì)生成這段呢，其實(shí)就是ANSYS根據(jù)你所定義的函數(shù)，自動(dòng)生成了一個(gè)Table表做了這個(gè)函數(shù)的插值，這樣系統(tǒng)在計(jì)算時(shí)就可以根據(jù)這個(gè)Table表進(jìn)行對應(yīng)的索引，生成任何你想要的函數(shù)值了。
展開 如何ANSYS的APDL中的模型導(dǎo)入HyperMesh ￥2
考慮到APDL中網(wǎng)格劃分功能并不是十分完善，需要借助HyperMesh軟件進(jìn)行網(wǎng)格劃分。那么如何將模型從ANSYS的APDL中導(dǎo)入到HyperMesh中呢
如何從Ansys APDL中提取剛度矩陣與質(zhì)量矩陣？ ￥69
1.引論
經(jīng)常使用Ansys、Abaqus等一系列有限元分析軟件進(jìn)行計(jì)算、學(xué)習(xí)的學(xué)生或工程師們都會(huì)知道在有限元分析建模與計(jì)算中剛度矩陣與質(zhì)量矩陣的重要性。但是由于軟件的黑盒性質(zhì)，大家往往在實(shí)際使用十分成熟的商業(yè)化軟件的過程中慢慢忽視了有限元及其衍生出的商業(yè)軟件背后的原理與方法。
這時(shí)，不管是在學(xué)習(xí)中還是在工程應(yīng)用中往往都會(huì)遇到一個(gè)同樣的問題，那么就是如何將Ansys APDL運(yùn)行中的產(chǎn)生的各種數(shù)據(jù)（例如：剛度矩陣、質(zhì)量矩陣）導(dǎo)出成為我們熟悉的形式或文件格式，從而為我們所用，所分析。
因此我決定寫下此篇文章來幫助很多實(shí)際工作或?qū)W習(xí)中需要用到此類技能的同學(xué)、同事們，讓大家更了解Ansys APDL背后的工作原理與數(shù)據(jù)導(dǎo)出方式。
當(dāng)然，在社區(qū)中早就有大佬回答過了這個(gè)問題，并給大家制作了相應(yīng)的提取矩陣軟件，其軟件具備了簡單、便捷的操作方式，讓很多想要提取剛度矩陣與質(zhì)量矩陣的同僚們受益，那么我為什么還要寫一篇這樣的文章重新提起這樣一個(gè)話題呢？這就又回到了我開頭所說的“原理與方法”，我在此更希望面對想要進(jìn)一步學(xué)習(xí)了解軟件背后機(jī)理的群體，并在此基礎(chǔ)上保留教學(xué)的簡潔性，提供導(dǎo)出矩陣與轉(zhuǎn)換、列式、求解的源代碼，使其既兼顧基本原理，又可以讓大家直接上手使用，非常的便捷，也避免了很多因?yàn)閮?yōu)化不完全導(dǎo)致的運(yùn)行bug。
2.有限元軟件導(dǎo)出剛度矩陣與質(zhì)量矩陣的方法
在使用APDL進(jìn)行求解時(shí)，每次在求解完成后都會(huì)在工作路徑下生成一個(gè).full文件，而這個(gè)文件十分關(guān)鍵，其正是剛度矩陣與質(zhì)量矩陣的所在之處。
展開 ANSYS APDL中的壓電分析
同理，彈性勁度常數(shù)也有10個(gè)：
CE11，CE12，CE13，CE33，CE44，CD11，CD12，CD13，CD33，CD44
ANSYS中，彈性系數(shù)矩陣為6×6矩陣（對2-D模型是4×4矩陣），它說明剛度系數(shù)（[c]矩陣）或柔度系數(shù)([s]矩陣)。
ANSYS的APDL中如何旋轉(zhuǎn)模型 ￥1
也可參考此處鏈接：ansys中旋轉(zhuǎn)模型
最后是如何變回原始坐標(biāo)系？
教程 - 機(jī)械 APDL 中的 2D 桁架分析 （ANSYS） 第 1 部分？
教程 - 機(jī)械 APDL 中的 2D 桁架分析 （ANSYS） 第 1 部分
一般來說，有限元解可以分為以下三個(gè)階段。
1. 預(yù)處理：定義問題;
- 定義關(guān)鍵點(diǎn)/線/區(qū)域/體積
- 定義元素類型和材料/幾何屬性
- 根據(jù)需要?jiǎng)澐志€/區(qū)域/體積
2. 解決方案：分配載荷、約束和求解;
3. 后處理：
- 節(jié)點(diǎn)位移列表
- 單元力和彎矩
- 撓度圖
- 應(yīng)力等值線圖
在本教程中，我們將進(jìn)行第一步。
步驟1：
啟動(dòng) Ansys Mechanical APDL。
步驟2：
單擊 Preferences 并選擇 Structural ，因?yàn)槲覀儗⑦M(jìn)行結(jié)構(gòu)分析。單擊 OK（確定）。
步驟3：
現(xiàn)在我們必須繪制關(guān)鍵點(diǎn)。在 Preprocessor >> Modeling >> Create >> In active CS 下創(chuàng)建。
步驟4：
現(xiàn)在我們必須輸入 Keypoints。輸入關(guān)鍵點(diǎn)編號(hào) 1 和 XYZ 坐標(biāo)，然后單擊 Apply。
步驟5：
輸入第二個(gè)關(guān)鍵點(diǎn) X=500，Y=1000。Z 將保持為零，因?yàn)槲覀冇?2D Bridge Truss。單擊 Apply。
步驟6：
輸入第三個(gè)關(guān)鍵點(diǎn) X=1000，Y=0。單擊 Apply。
步驟7：
輸入第 4 個(gè)關(guān)鍵點(diǎn) X=1500，Y=1000。單擊 Apply。
步驟8：
輸入第 5 個(gè)關(guān)鍵點(diǎn) X=2000，Y=0。單擊 OK
步驟9：
現(xiàn)在我們已經(jīng)繪制了關(guān)鍵點(diǎn)。我們必須沿著這些關(guān)鍵點(diǎn)創(chuàng)建線條。轉(zhuǎn)到 建模 >> 在激活坐標(biāo)中>>創(chuàng)建>>線。
步驟10：
現(xiàn)在通過單擊它們來選擇 kepoint，然后單擊其他關(guān)鍵點(diǎn)以創(chuàng)建線。創(chuàng)建成員。單擊 OK（確定）。
步驟11：
現(xiàn)在我們必須定義 Element 類型。即 Beam。
展開 教程 - 機(jī)械 APDL 中的 2D 桁架分析 （ANSYS） 第 2 部分
后處理：
- 節(jié)點(diǎn)位移列表
- 單元力和彎矩
- 撓度圖
- 應(yīng)力等值線圖
在本教程中，我們將進(jìn)行第二步和第三步。
1. 步驟1：
這是教程的第二部分，我們在其中解決問題。在 Solution >> Analaysis 下，鍵入 New analysis>>。選擇 static 并單擊 OK。
2. 步驟2：
在定義載荷下>>>> Structural >> 位移 >> On 關(guān)鍵點(diǎn)上應(yīng)用。現(xiàn)在，我們將定義固定的關(guān)鍵點(diǎn)或支撐。
3. 步驟3：
選擇兩個(gè)下角關(guān)鍵點(diǎn)，然后單擊 OK。
4. 步驟4：
選擇 All DOF 并單擊 OK。
5. 步驟5：
轉(zhuǎn)到定義載荷 >> 在關(guān)鍵點(diǎn)上應(yīng)用>> 結(jié)構(gòu)>>力矩/力矩 >> 。
6. 步驟6：
選擇上部關(guān)鍵點(diǎn)，然后單擊 OK。
7. 步驟7：
力的方向?yàn)?FY 且輸入 Force 值 = -10000，因?yàn)榱⑾蛳伦饔谩?8. 步驟8：
現(xiàn)在我們已經(jīng)準(zhǔn)備好了模型進(jìn)行求解。在 Solve 下>> Current Load 步驟。
9. 步驟9：
單擊 OK（確定）。
10. 步驟10：
一條消息 Solution is done！將顯示流程何時(shí)完成。單擊 Close。
11. 步驟11：
現(xiàn)在是這個(gè)過程的第三部分。要進(jìn)行后處理。轉(zhuǎn)到 General PostProc >> 列出結(jié)果 >> reaction solu。
12. 步驟12：
選擇所有項(xiàng)目，然后單擊 OK。
13. 步驟13：
現(xiàn)在我們有了關(guān)于節(jié)點(diǎn) 1 和 5 的 reaction 值，它們是固定的。單擊 Close。
14.
展開 ANSYS APDL中判斷和改變奇偶性的方法
在C 或者matlab中非常容易判斷。但在APDL中，沒有判斷函數(shù)。搜了下，結(jié)果提示用nint(x)函數(shù)。help中說到 nint(x)就是abs,sin,cos等函數(shù)一樣可以直接用來運(yùn)算。nint(x)的意思是‘Nearest interger to x’，意思是說最接近x的整數(shù)。比如nint(2.1)=2，nint(2.5)=3，就相當(dāng)于四舍五入。基于此，我就用它來作為判斷奇數(shù)偶數(shù)的工具了。
實(shí)現(xiàn)的命令如下：
!change TotalNumber to even number
*if,abs(nint(TotalNumber/2)-TotalNumber/2),le,0.3,then ! 不一定是0.3， 只要是小于0.5 就可以了。
TotalNumber=TotalNumber
*else
TotalNumber=TotalNumber+1
*endif
展開 ansys workbench APDL熱輻射命令行中的有關(guān)說明求助
1.sf,nlist,label,value,value2
-“nilst”是節(jié)點(diǎn)列表，也可以是命名選擇
-輻射標(biāo)簽是rdsf
-value是表面發(fā)射率
-value2是封閉體數(shù)量
2.spctemp命令行：因?yàn)樗?jì)算的空間不是完全封閉的計(jì)算空間，所以必須定義空間溫度，
spctemp,number,temperature
spctemp是ansys定義空間溫度的關(guān)鍵字，number是非封閉空間的數(shù)量，temperature是非封閉空間的溫度
3.stef命令行：stef是ansys中斯蒂芬玻爾茲曼常數(shù)，stef=5.67×10-8
4.RADOPT, FLUXRELX, FLUXTOL, SOLVER, MAXITER, TOLER, OVERRLEX
FLUXRELX：松弛因子。
FLUXTOL：輻射熱通量收斂容差，默認(rèn)為0.0001。
SOLVER
選擇用于計(jì)算的輻射求解器:
0 – Gauss-Seidel求解器
1 – 直接求解器 (對于大問題將耗費(fèi)很多時(shí)間)
MAXITER
Gauss Seidel迭代求解器的最大迭代次數(shù) (SOLVER = 0)，默認(rèn)為1000 Gauss Seidel迭代求解器的最大迭代次數(shù) (SOLVER = 0)，默認(rèn)為1000。
TOLER
Gauss Seidel迭代求解器的收斂容差(SOLVER = 0)，默認(rèn)為 0.1。
OVERRLEX
Gauss Seidel迭代求解器的松弛因子(SOLVER = 0)，默認(rèn)為0.1。
求助：以上的封閉體數(shù)量是如何判別的？非封閉空間的數(shù)量又是如何判斷的？非封閉空間的溫度是如何定義的？有人能幫忙進(jìn)一步舉例或說明嗎？萬分感謝！
展開 ANSYS apdl中如何使齒輪饒定軸轉(zhuǎn)動(dòng)
問題描述：在ANSYS中，實(shí)體單元和平板單元只有平動(dòng)自由度，無旋轉(zhuǎn)自由度，網(wǎng)上提供的方法，在柱坐標(biāo)系下通過固定徑向位移，在周向施加小位移，但這種方法只適合小位移的轉(zhuǎn)動(dòng)，無法實(shí)現(xiàn)大位移，本文提出采用MPC184-銷軸單元和MPC-184剛性梁單元完成。本文就GUI的方式來介紹如何來創(chuàng)建齒輪的繞定軸旋轉(zhuǎn)。達(dá)到下圖的效果
step1 定義單元類型
a plat182 單元 模擬齒輪
b mpc184-剛性梁單元
c mpc184-銷軸單元 （本文繞Z軸旋轉(zhuǎn)，如定義為x軸，需要定義局部坐標(biāo)系，繞y軸旋轉(zhuǎn)90°）
step2 定義局部坐標(biāo)系默認(rèn)，本文定義12號(hào)
step3 定義銷軸截面以及單元坐標(biāo)系
step3 創(chuàng)建銷軸連接單元
在齒輪的中心點(diǎn)分配3號(hào)銷軸單元
step4 創(chuàng)建剛性梁單元
單元屬性旋旋轉(zhuǎn)2號(hào)剛性梁單元，去內(nèi)徑的節(jié)點(diǎn)和圓中心點(diǎn)創(chuàng)建剛性梁單元
step5 施加載荷使齒輪旋轉(zhuǎn)2圈
setp6 求解設(shè)置
step7 時(shí)間后處理選擇內(nèi)徑上的一個(gè)節(jié)點(diǎn)觀察ux，uy，rotz，可以看出齒輪旋轉(zhuǎn)2圈
展開 關(guān)于SYNTAX編輯器在ANSYS APDL編輯中的應(yīng)用
Syntax-highlighter，這款可以當(dāng)作ANSYS APDL 命令編輯器使用的軟件有誰使用過嗎?我今天下載了，看了那個(gè)說明的PDF忙了半天，總算能提示命令信息,能高亮顯示，能改變顏色和設(shè)置鍵盤快捷方式。但是還是有一些地方不明白：
1、這個(gè)是1、7版本的，都5年過去了，肯定有新版本，有誰有新版本嗎？為什么用GOOGLE搜索也搜索不到新版本呢？
2、里面有三個(gè)顯示方式？HIGH-DCI，SETUP-A57、DEFAULT。但是當(dāng)用SETUP-A57時(shí)就自動(dòng)出錯(cuò)。我打開看是因?yàn)镾ETUP-A57里面有個(gè)絕對地址的緣故，但是我找不到這個(gè)地址所指的文件，所以沒有辦法更改地址。其他幾個(gè)文件中也有一些找不到文件的地址。并且很多操作經(jīng)常出錯(cuò)。有誰解決這些問題了嗎？
大家討論討論哦，咱們也可以一起努力打造一個(gè)符合中國人，特別是APDL新手的編輯器。
展開 
Ansys Workbench后處理中，利用APDL命令提取繞圓柱坐標(biāo)系的扭矩角度 ￥10
問題：
在有限元仿真中有時(shí)需要提取某些結(jié)構(gòu)的扭轉(zhuǎn)角度。Ansys workbench的結(jié)果后處理中可以設(shè)定圓柱坐標(biāo)系，然后按圓柱坐標(biāo)讀取Y軸的變形結(jié)果，再進(jìn)行扭轉(zhuǎn)角度的換算。
本文這里將該過程利用APDL命令進(jìn)行處理，避免一下步驟重復(fù)操作。
? 每次要單獨(dú)記錄變形量，
? 還要測量關(guān)鍵節(jié)點(diǎn)到坐標(biāo)系原點(diǎn)的距離，
? 將變形量和距離進(jìn)行角度換算（弧度）
? 弧度角轉(zhuǎn)角度
APDL后處理命令功能介紹：
1. 在坐標(biāo)系中創(chuàng)建所需的圓柱坐標(biāo)系，并在屬性ADPL name中進(jìn)行命名：aix （用戶隨意命名）
2. 在Named selection 定義需要查看的區(qū)域，并命名：load（用戶隨意命名）
3. 在后處理中插入command 命令，并將上述坐標(biāo)系和NS的名稱修改。
4. 在command的結(jié)果屬性中就會(huì)有最大/最小/平均扭轉(zhuǎn)角度。并且為了方便校核準(zhǔn)確性還提供了沿圓柱坐標(biāo)系Y軸的變形量。
并且，除了界面顯示的結(jié)果外，還會(huì)在WB的結(jié)果文件夾中，顯示named Selection區(qū)域所有節(jié)點(diǎn)的編號(hào)/距離選定坐標(biāo)系的距離/沿坐標(biāo)系Y軸的變形量/換算后的角度值等信息，以便進(jìn)行其它數(shù)據(jù)處理。
展開 在ANSYS WORKBENCH中插入APDL命令例子--多載荷步的例子
【問題分析】
本問題可以直接在wb中用多載荷步來求解，這里說明如何使用插入APDL命令的方式實(shí)現(xiàn)。
【求解過程】
1. 打開ANSYS WORKBENCH14.5
2.創(chuàng)建結(jié)構(gòu)靜力學(xué)分析系統(tǒng)。
3.創(chuàng)建幾何體。
雙擊geometry單元格，進(jìn)入DM,選擇mm單位。
創(chuàng)建長方體。
其尺寸設(shè)置是
退出DM.
4.劃分網(wǎng)格。
雙擊MODEL,進(jìn)入到MECHANICAL中，按照默認(rèn)方式劃分網(wǎng)格。
5.固定左端面。
6.添加APDL命令以分步加載。
下面使用APDL命令進(jìn)行分步加載。
由于該命令最后要傳遞到經(jīng)典界面中計(jì)算，而經(jīng)典界面沒有單位。為保持統(tǒng)一性，都用毫米單位。
（1）設(shè)置單位
（2）創(chuàng)建命名集。
由于在命令中要引用頂面這個(gè)面，為了能夠正確引用，先需要給它一個(gè)名稱，這需要使用命名集來完成。
選擇上述頂面，創(chuàng)建命名集。在彈出的對話框中設(shè)置名字：topface
則樹形大綱中出現(xiàn)了該命名集。
有了命名集，在后面就可以使用該名字了。
（3）插入APDL命令。
在數(shù)形大綱中先選擇A5，再從工具欄中選擇命令按鈕
則圖形窗口變成了一個(gè)文本編輯器，此處可以輸入命令。
該文本窗口內(nèi)說了很多話，主要內(nèi)容包含兩點(diǎn)：
第一，這些命令會(huì)在SOLVE命令剛執(zhí)行前執(zhí)行。
第二，注意這里用的單位是mm.
現(xiàn)在我們向該文本窗口輸入下列命令。
這段ADPL命令流的含義是：
首先退出前面的某個(gè)處理器(finish)
然后進(jìn)入到求解器中（/solve），在1,2,3,個(gè)時(shí)間步，依次在頂面上施加1,2,3mpa的載荷(sf)，并將該載荷步寫入到載荷步文件中（lswrite）,然后先后求解這三個(gè)載荷步（lssolve）。
展開 ansys隨機(jī)振動(dòng)沒有結(jié)果
用ansys中的apdl模塊進(jìn)行隨機(jī)振動(dòng)分析，計(jì)算完后，結(jié)果顯示中沒有載荷步2，3，4，5，還是只有載荷步1中模態(tài)分析的結(jié)果。請問有誰知道是什么原因嗎？我用的是19.0點(diǎn)版本。
生成不同填充率的模型-隨機(jī)函數(shù)的使用 ￥299
骨料填充模型的繪制方法-ANSYS APDL命令的使用-不同形狀-不同大小的圖形填充
在材料科學(xué)和工程領(lǐng)域，我們經(jīng)常面臨一個(gè)挑戰(zhàn)：如何在材料內(nèi)部隨機(jī)生成加強(qiáng)骨料或缺陷孔隙。這種隨機(jī)性的引入對于模擬材料的真實(shí)行為至關(guān)重要，因?yàn)樗梢愿玫胤从巢牧显趯?shí)際應(yīng)用中的復(fù)雜性和不確定性。本文將介紹如何使用ANSYS軟件中的APDL（ANSYS Parametric Design Language）命令來實(shí)現(xiàn)這一需求，并通過一個(gè)簡單的平面案例來詳細(xì)解析操作步驟。
一、問題概述與關(guān)鍵點(diǎn)
在建模過程中，隨機(jī)生成加強(qiáng)骨料或缺陷孔隙的問題可以簡化為在指定區(qū)域內(nèi)隨機(jī)放置圖形的問題。這個(gè)過程需要考慮以下幾個(gè)關(guān)鍵點(diǎn)：
填充率：填充率是指生成的圖形占整個(gè)模型區(qū)域的比例。根據(jù)實(shí)際需求，填充率可以設(shè)置為0.1、0.5或其他值。
互不干涉：生成的圖形之間不能發(fā)生干涉，否則會(huì)導(dǎo)致模型失效。因此，需要確保每個(gè)圖形的位置和大小都是合理的，以避免重疊。
隨機(jī)性：填充的圖形位置必須是隨機(jī)的，以模擬材料內(nèi)部的隨機(jī)分布。同時(shí)，可以考慮形狀的隨機(jī)性，例如全部為圓形、全部為正方形或混合形狀，并且大小也可以隨機(jī)變化。
二、保障填充率
為了保障填充率，我們可以采用累加的方法來確定是否超過全部圖形的比例。具體步驟如下：
初始化一個(gè)變量來記錄已繪制的圖形面積。
在每次生成一個(gè)圖形后，將其面積累加到該變量中。
通過比較已繪制圖形面積與模型總面積的比例，判斷是否達(dá)到設(shè)定的填充率。如果沒有達(dá)到，則繼續(xù)繪制圖形；否則停止繪圖。
三、確保圖形互不重合
為了確保生成的圖形之間互不重合，我們可以采用以下策略：
定義一個(gè)數(shù)組來存儲(chǔ)已經(jīng)生成的圖形的圓心坐標(biāo)和半徑。
展開 