 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 1 | 1/1 | 返回列表 |
| 查看: 551 | 回復(fù): 0 | ||
| 【懸賞金幣】回答本帖問題，作者cjc021206將贈(zèng)送您 5 個(gè)金幣 | ||
cjc021206新蟲 (初入文壇)
|
[求助]
求助UDF函數(shù)問題
|
|
|
我用動(dòng)網(wǎng)格根據(jù)零件位置計(jì)算力矩，確定他的運(yùn)動(dòng)狀態(tài)，UDF內(nèi)容是 int damp_f( real theta, real alpha, real* vol, real* r_m ){ real x0=-300, y0=137.5, x1=-80, y1=-600; real x0r, y0r, dx, dy, d01, ddf; x0r = x0 * cos( theta ) - y0 * sin( theta ); y0r = x0 * sin( theta ) + y0 * cos( theta ); dx = x0r - x1; dy = y0r - y1; d01 = sqrt( pow(dx, 2) + pow(dy, 2) ); *vol = d01 * alpha; *r_m = fabs( (y1 * x0r - x1 * y0r) / d01 ); return 0; } DEFINE_CG_MOTION( piston_damp_sp0, dt, vel, omega, time, dtime ) { face_t f; cell_t c; real f_glob[3], m_glob[3], domega_z, x_cg[ND_ND], m_damp, v_damp; int i, status; NV_S(omega, =, 0.0); for(i=0;i<=3;i++) { f_glob=0; m_glob=0; } x_cg[0] = 0; x_cg[1] = 0; Domain *domain= Get_Domain (1); Thread *tf1 = Lookup_Thread (domain, 3); Compute_Force_And_Moment (domain, tf1, x_cg, f_glob, m_glob, TRUE); real* vol; real* r_m; status = damp_f( theta_z0, alpha_z0, vol, r_m ); v_damp = 4 * (*vol) * (*r_m) ; m_damp = 0; alpha_z0 = (m_glob[2] - m_damp - v_damp )/INERTIA; domega_z = dtime * alpha_z0; omega_z0 += domega_z; theta_z0 = DT_THETA(dt)[2]; Message ("ntime=%f, tz=%.10lf, oz1=%.10lf, oz2=%.10lf, force_x=%.1f ,force_y=%.1f, moment_z=%f, vol=%.1f, rm=%.1f ", \ time, DT_THETA(dt)[2], DT_OMEGA_CG(dt)[2], theta_z0, f_glob[0], f_glob[1], m_glob[2], *vol, *r_m ); omega[2]=omega_z0; } 用多核計(jì)算，開始后直接報(bào)錯(cuò)， Node 0: Process 5204: Received signal SIGSEGV. ============================================================================== ============================================================================== Node 1: Process 9392: Received signal SIGSEGV. ============================================================================== MPI Application rank 0 exited before MPI_Finalize() with status 2 把函數(shù) damp_f 取消，內(nèi)容加到 DEFINE_CG_MOTION 里面就沒問題，請教UDF里面加函數(shù)應(yīng)該怎么寫，還是說UDF里面不能加函數(shù) |
回復(fù)此樓» 猜你喜歡
085600材料與化工 292分 求調(diào)劑
已經(jīng)有3人回復(fù)
-
接受26屆調(diào)劑生
已經(jīng)有8人回復(fù)
-
273求調(diào)劑
已經(jīng)有7人回復(fù)
-
A區(qū)一本交叉課題組，低分調(diào)劑，招收機(jī)械電子信息通信等交叉方向
已經(jīng)有45人回復(fù)
-
334求調(diào)劑
已經(jīng)有3人回復(fù)
-
求調(diào)劑推薦
已經(jīng)有3人回復(fù)
-
求調(diào)劑
已經(jīng)有4人回復(fù)
-
一志愿211，0860總分286食品工程專業(yè)求調(diào)劑
已經(jīng)有6人回復(fù)
-
招調(diào)劑
已經(jīng)有14人回復(fù)
-
誠求調(diào)劑，323分 有專利/科研/實(shí)習(xí)經(jīng)歷
已經(jīng)有6人回復(fù)
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 cjc021206 的主題更新
| 1 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 接受26屆調(diào)劑生 +8 | 豬豬豬毛 2026-03-06 | 8/400 |
|
|---|---|---|---|---|
|
|
[考研] A區(qū)一本交叉課題組，低分調(diào)劑，招收機(jī)械電子信息通信等交叉方向 +32 | lisimayy 2026-03-04 | 45/2250 |
|
|
|
[考研] 誠求調(diào)劑，323分 有專利/科研/實(shí)習(xí)經(jīng)歷 +5 | 橙子cyx 2026-03-06 | 6/300 |
|
|
|
[考研] 085600材料與化工（高分子）290分求調(diào)劑 +9 | wengyujian 2026-03-04 | 9/450 |
|
|
[考研] 一志愿山東大學(xué)105500藥學(xué)專碩，總分302求調(diào)劑 +4 | 五維天空 2026-03-04 | 8/400 |
|
|
|
[考研] 材料專碩調(diào)劑 +4 | 慕辰123 2026-03-05 | 7/350 |
|
|
|
[考研] 一志愿蘇大材料工程專碩293求調(diào)劑 +3 | 鋼鐵大炮 2026-03-04 | 3/150 |
|
|
|
[考研] 材料085600 303求調(diào)劑 +7 | 1bygone 2026-03-04 | 7/350 |
|
|
|
[考研] 材料學(xué)碩080500復(fù)試調(diào)劑294 +3 | 四葉zjz 2026-03-04 | 3/150 |
|
|
|
[考研] 求調(diào)劑 +7 | 博斯特525 2026-03-04 | 7/350 |
|
|
|
[考研] 085600 材料與化工 298 +14 | 小西笑嘻嘻 2026-03-03 | 14/700 |
|
|
|
[考研] 一志愿中科大080500總分324求調(diào)劑 +3 | jorna 2026-03-03 | 6/300 |
|
|
|
[考研]
 一志愿中南大學(xué)理學(xué)化學(xué)
+6
|
15779376950 2026-03-01 | 8/400 |
|
|
|
[考研] 一志愿314求調(diào)劑 +7 | 202111120625 2026-03-03 | 7/350 |
|
|
|
[考研] 材料化工調(diào)劑 +15 | 今夏不夏 2026-03-01 | 18/900 |
|
|
|
[考研] 0854復(fù)試調(diào)劑 276 +5 | wmm9 2026-03-01 | 7/350 |
|
|
|
[考研] 261求調(diào)劑 +3 | 陸lh 2026-03-01 | 3/150 |
|
|
|
[考研] 一志愿華南理工大學(xué)材料與化工326分，求調(diào)劑 +3 | wujinrui1 2026-02-28 | 3/150 |
|
|
|
[考研] 調(diào)劑 +3 | 13853210211 2026-03-02 | 4/200 |
|
|
|
[基金申請] 成果系統(tǒng)訪問量大，請一小時(shí)后再嘗試。---NSFC啥時(shí)候好哦，已經(jīng)兩天這樣了 +4 | NSFC2026我來了 2026-02-28 | 4/200 |
|