構造力学

構造力学

APDLを使用して、周波数応答解析から配列に節点セットにおける振幅を書き出す方法を教えてください。

    • FAQFAQ
      Participant

      /filnam,full,db /PREP7 rectng,0,100,0,15 et,1,63 r,1,2 lsel,s,loc,x,5,95 lesize,all,,,10 allsel esize,5 amesh,all ex,1,200e3 nuxy,1,0.3 dens,1,7800e-12 nsel,s,loc,x,0 d,all,all allsell save /solu antype,modal modopt,lanb,10 mxpand,10 solve finish /SOLU ANTYPE,3 HROPT,FULL HROUT,OFF LUMPM,0 HARFRQ,10,2000, NSUBST,100, g=9810 ACEL,0,0,5*g KBC,1 DMPRAT,0.03, solve finish /post26 FILE, ,’rst’ /UI,COLL,1 NUMVAR,200 SOLU,191,NCMIT STORE,MERGE FILLDATA,191,,,,1,1 REALVAR,191,191 PLCPLX,0! 振幅をプロット PRCPLX,1! 振幅および位相角をプリント NSEL,S, , , 12 NSEL,a, , , 17 NSEL,a, , , 41 NSEL,a, , , 44 ! さらに多くの節点を追加して、後でdoループで使用するコンポーネントを準備 CM,nodeset,NODE CMSEL,s,nodeset *get,ncount,node,,count allsel *GET,freqs,VARI, ,NSETS *DIM,myarray,ARRAY,freqs,ncount+1 CMSEL,s,nodeset nnum=0 *do,count,1,ncount,1 nnum=ndnext(nnum) NSOL,2,nnum,U,Z,Tip_Top_UZ ABS,10,2 !PRVAR,2 VGET,f,1,0,0 VGET,myarray(1,1),1,0,0 VGET,myarray(1,count+1),10,0,0 VARDEL,2 VARDEL,10 *enddo ! NSOL,2,12,U,Z,Tip_12 ! NSOL,3,17,U,Z,Tip_17 ! prvar,2,3