24小时热门版块排行榜    

北京石油化工学院2026年研究生招生接收调剂公告
查看: 955  |  回复: 6
【悬赏金币】回答本帖问题,作者1216565504将赠送您 15 个金币
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

1216565504

铜虫 (正式写手)

[求助] 建模遇到的问题解决不了了,求助大佬们帮帮忙呀

我的模型是一个正方体内有许多球体,我想用vovlap命令进行搭接,但是里面的球体数量多了以后就会出现拓扑错误,让我修改几何形状或公差,我改了btol后还是不行。该怎么办呢?

建模遇到的问题解决不了了,求助大佬们帮帮忙呀


@wuming524 发自小木虫Android客户端
回复此楼

» 猜你喜欢

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

1216565504

铜虫 (正式写手)

/clear
!****************************************************************************
!文件说明段
!****************************************************************************
/FILNAME,CM,1
/TITLE, Composite material
/REPLOT       
/units,si
!****************************************************************************
!设置材料属性
!****************************************************************************
/PREP7
ET,1,SOLID70
MP,KXX,1,0.27
MP,KXX,2,290
!****************************************************************************
!绘制分析模型
!****************************************************************************
XMAX=100
YMAX=100
ZMAX=100    !设定模型尺寸

*Create,dataread,mac  !创建宏文件,名为dataread,名称可更改
*dim,PC,array,96,4 !创建读入数据的数组
*vread,PC(1,1), G:\Matlab\data,txt,,jik,4,96 !读入文件data.txt的数据
(4f5.0)                                 !读入格式
*End                                 !结束宏文件的创建
Dataread                              !运行宏文件,读入颗粒坐标数据

*DO,I,1,96
!生成球体
WPOFFS,PC(I,1),PC(I,2),PC(I,3)
SPH4,0,0,PC(I,4)

*GET,VOLU_NUM,VOLU,0,COUNT
M= VOLU_NUM
WPROTA,30,
CSWPLA,11,1                 
CSYS,11
VGEN, ,M, , ,, 30, , , ,1

CSYS,0
WPCSYS,-1,0
*ENDDO
ALLS
VSUM,VOLUME
*GET,VOLUME_VOLU,VOLU,,VOLU
ALLS
VPLOT
BLOCK,0,XMAX,0,YMAX,0,ZMAX
!Vovlap,all
!/TRLCY,VOLU,1,74

!!!!!!!!!!!!!!!!!!!!!截掉体外面一部分的填充物!!!!!!!!!!!!!!!!
WPOFFS,0
WPROTA,,,90
BOPTN,KEEP,0          !不保留原体
BOPTN,NWARN,1   
BOPTN,VERS,RV52
BTOL,1e-005,
VSBW,ALL
CSYS,0
WPCSYS,-1,0

WPOFFS,XMAX
WPROTA,,,90
BOPTN,KEEP,0          !不保留原体
BOPTN,NWARN,1   
BOPTN,VERS,RV52
BTOL,1e-005,
VSBW,ALL
CSYS,0
WPCSYS,-1,0                        !使工作平面位置与当前坐标系平面一致

WPOFFS,,0
WPROTA,,90,
BOPTN,KEEP,0          !不保留原体
BOPTN,NWARN,1   
BOPTN,VERS,RV52
BTOL,1e-005,
VSBW,ALL
CSYS,0
WPCSYS,-1,0

WPOFFS,,YMAX
WPROTA,,90,
BOPTN,KEEP,0          !不保留原体
BOPTN,NWARN,1   
BOPTN,VERS,RV52
BTOL,1e-005,
VSBW,ALL
CSYS,0
WPCSYS,-1,0

WPOFFS,,,ZMAX
BOPTN,KEEP,0          !不保留原体
BOPTN,NWARN,1   
BOPTN,VERS,RV52
BTOL,1e-005,
VSBW,ALL
CSYS,0
WPCSYS,-1,0

WPOFFS,,,0
BOPTN,KEEP,0          !不保留原体
BOPTN,NWARN,1   
BOPTN,VERS,RV52
BTOL,1e-005,
VSBW,ALL
CSYS,0
WPCSYS,-1,0

!反选删掉四周的体
VSEL,U,LOC,X,0,XMAX
VDELE,ALL,,,1
ALLS
VSEL,U,LOC,Y,0,YMAX
VDELE,ALL,,,1
ALLS
VSEL,U,LOC,Z,0,ZMAX
VDELE,ALL,,,1
ALLS
VPLOT
这是第二种常出错的方法。思路是先用matlab算了球体的坐标和半径参数,读入到APDL中建模。
5楼2022-01-12 17:09:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 7 个回答

AUG突击

新虫 (正式写手)

2楼2022-01-12 16:49:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

1216565504

铜虫 (正式写手)

引用回帖:
2楼: Originally posted by AUG突击 at 2022-01-12 16:49:14
试试分多次操作

好主意!我试试一个个搭接看看哪里出问题了
3楼2022-01-12 17:06:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

1216565504

铜虫 (正式写手)

另外奇怪的一点是,我用另外一种算法建模出现这个问题的概率很小。对比两种方法没找到区别。希望大佬们帮我看一下。
/clear
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
XMAX=100
YMAX=100
ZMAX=100    !设定模型尺寸
RMIN=10
RMAX=10     !设定球形半径
/PREP7
NUM_TOTAL=100    !设定球个数
*DO,I,1,1000000
LOC_X=RAND(-XMAX/2,XMAX/2)
LOC_Y=RAND(-YMAX/2,YMAX/2)
LOC_Z=RAND(-ZMAX/2,ZMAX/2)
RADIUS=RAND(RMIN,RMAX)
!生成球体
WPOFFS,LOC_X,LOC_Y,LOC_Z
SPH4,0,0, RADIUS
*GET,VOLU_NUM,VOLU,0,COUNT
M= VOLU_NUM

WPROTA,30,
CSWPLA,11,1
CSYS,11
VGEN, ,M, , ,, 30, , , ,1

CSYS,0
WPCSYS,-1,0
!通过布尔运算判断球是否相交
*IF,M,GT,1,THEN            
                *DO,J,1,M-1
BOPTN,KEEP,1
BOPTN,NWARN,1   
BOPTN,VERS,RV52
BTOL,1e-005,
VSBV,J,M
*GET,VOLU_NUM,VOLU,0,COUNT
*IF,VOLU_NUM,GT,M,THEN
VDELE,M,M+1,1,1
NUMCMP,VOLU
*EXIT
                        *ENDIF
*ENDDO
*ENDIF
*GET,VOLU_NUM,VOLU,0,COUNT
M= VOLU_NUM
!判断数目如果达到一定数目停止循环
*IF,M,GT,NUM_TOTAL,THEN
*EXIT
*ENDIF
*ENDDO
BLOCK,-XMAX/2,XMAX/2,-YMAX/2,YMAX/2,-ZMAX/2,ZMAX/2,
/TRLCY,VOLU,0.5,     M+1
/REPLOT


!!!!!!!!!!!!!!!!!!!!!截掉体外面一部分的填充物!!!!!!!!!!!!!!!!
WPOFFS,-XMAX/2
WPROTA,,,90
BOPTN,KEEP,0          !不保留原体
BOPTN,NWARN,1   
BOPTN,VERS,RV52
BTOL,1e-005,
VSBW,ALL
CSYS,0
WPCSYS,-1,0

WPOFFS,XMAX/2
WPROTA,,,90
BOPTN,KEEP,0          !不保留原体
BOPTN,NWARN,1   
BOPTN,VERS,RV52
BTOL,1e-005,
VSBW,ALL
CSYS,0
WPCSYS,-1,0

WPOFFS,,YMAX/2
WPROTA,,90,
BOPTN,KEEP,0          !不保留原体
BOPTN,NWARN,1   
BOPTN,VERS,RV52
BTOL,1e-005,
VSBW,ALL
CSYS,0
WPCSYS,-1,0

WPOFFS,,-YMAX/2
WPROTA,,90,
BOPTN,KEEP,0          !不保留原体
BOPTN,NWARN,1   
BOPTN,VERS,RV52
BTOL,1e-005,
VSBW,ALL
CSYS,0
WPCSYS,-1,0

WPOFFS,,,ZMAX/2
BOPTN,KEEP,0          !不保留原体
BOPTN,NWARN,1   
BOPTN,VERS,RV52
BTOL,1e-005,
VSBW,ALL
CSYS,0
WPCSYS,-1,0

WPOFFS,,,-ZMAX/2
BOPTN,KEEP,0          !不保留原体
BOPTN,NWARN,1   
BOPTN,VERS,RV52
BTOL,1e-005,
VSBW,ALL
CSYS,0
WPCSYS,-1,0

!反选删掉四周的体
VSEL,U,LOC,X,-XMAX/2,XMAX/2
VDELE,ALL,,,1
ALLS
VSEL,U,LOC,Y,-YMAX/2,YMAX/2
VDELE,ALL,,,1
ALLS
VSEL,U,LOC,Z,-ZMAX/2,ZMAX/2
VDELE,ALL,,,1
ALLS
VPLOT
NUMCMP,VOLU
NUMCMP,KP
NUMCMP,LINE
NUMCMP,AREA

这是第一种出现错误较少的方法
4楼2022-01-12 17:08:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
不应助 确定回帖应助 (注意:应助才可能被奖励,但不允许灌水,必须填写15个字符以上)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 316求调剂 +9 墨辰_Orion926 2026-04-04 9/450 2026-04-04 21:35 by lbsjt
[考研] 368求调剂 +5 今华习 2026-04-03 7/350 2026-04-04 18:47 by imissbao
[考研] 085600,320分求调剂 +14 大馋小子 2026-04-04 15/750 2026-04-04 16:27 by 无际的草原
[考研] 293分求调剂,外语为俄语 +6 加一一九 2026-03-31 6/300 2026-04-04 14:57 by 聪明的大松鼠
[考研] 怎么删帖子啊 +3 缝曦1000 2026-04-04 3/150 2026-04-04 14:20 by 土木硕士招生
[考研] 297求调剂 +11 ljy20040718! 2026-04-03 13/650 2026-04-04 09:23 by 来看流星雨10
[考研] 材料科学与工程考研 +10 拯救皮特托先生 2026-04-02 10/500 2026-04-03 23:57 by userper
[考研] 281求调剂 +10 aaawhy 2026-04-03 10/500 2026-04-03 21:42 by lbsjt
[考研] 0705理学294求调剂 +3 成果成果cg5 2026-04-03 3/150 2026-04-03 14:04 by simons1972
[考研] 求调剂 +3 心想事成可 2026-04-03 3/150 2026-04-03 11:22 by wangjy2002
[考研] 262求调剂 +6 励志一定发文章 2026-04-02 7/350 2026-04-03 09:54 by linyelide
[考研] 285求调剂 +8 AZMK 2026-04-02 11/550 2026-04-02 20:16 by yulian1987
[考博] 材料工程专业硕士申博 +3 麟正宇 2026-03-30 3/150 2026-04-02 15:04 by greychen00
[考研] 266分,一志愿电气工程,本科材料,求材料专业调剂 +10 哇呼哼呼哼 2026-04-01 11/550 2026-04-02 11:31 by lnilvy
[考研] 0710生物学求调剂 +9 manman511 2026-04-01 9/450 2026-04-02 10:00 by zxl830724
[考研] 070300化学279求调剂 +15 哈哈哈^_^ 2026-03-31 17/850 2026-04-01 21:37 by 给你你注意休息
[考研] 303分 0807学硕求调剂 +3 TYC3632 2026-04-01 3/150 2026-04-01 19:24 by lwk2004
[考研] 349求调剂 +6 吃的不少 2026-04-01 6/300 2026-04-01 17:55 by JYD2011
[考研] 085404 22408 315分 +5 zhuangyan123 2026-03-31 6/300 2026-03-31 13:48 by limeifeng
[考博] 材料专业申博 +5 杜雨婷dyt 2026-03-29 5/250 2026-03-31 11:19 by oooqiao
信息提示
请填处理意见