求助：焊接热结构间接分析耦合

hnpsd

不知道你的问题解决了没有?我正在做这一块,很迷茫中
帮你顶上去,希望高手指教啊

zgmgary

这是我做的移动温度场焊接加热阶段的命令流，请各位多多指教！

/CLEAR,NOSTART
/TITLE,TemFieldANALYSIS
/FILNAM,TemField
/UNITS,SI
/CONFIG,NRES,10000


/PREP7
ET,1,SOLID70
TOFFST,273
mptemp,1,20,250,500,750,1000,1500
mptemp,7,1700,2500
mpdata,c,1,1,460,480,530,675,670,660
mpdata,c,1,7,780,820
mpdata,kxx,1,1,50,47,40,27,30,35
mpdata,kxx,1,7,140,142
mpdata,dens,1,1,7820,7700,7610,7550,7490,7350
mpdata,dens,1,7,7300,7090
mpdata,hf,1,1,100,350,520,1000,1500,3000
mpdata,hf,1,7,3100,3500

/VIEW,1,1,1,1
BLOCK,0,0.300,0,0.010,0,0.060
SAVE
LESIZE,1,,,5
LESIZE,4,,,150
LESIZE,9,,,30
MSHKEY,1
VATT,1,1,1
VMESH,1  
/PNUM,DEFA
EPLOT
FINISH

/SOLU
LSIZE=0.002
V=0.010
TINC=LSIZE/V
KK=3
Qmax=156668147
R=0.004
MAX_TIME=151
MAX_X=151
MAX_Y=6
   
*DIM,FLUX2,TABLE,MAX_X,MAX_Y,MAX_TIME,X,Y,TIME
*DO,K,1,MAX_TIME,1
  *DO,I,1,MAX_X,1
  FLUX2(I,0,K)=(I-1)*LSIZE
  *ENDDO
   
  *DO,J,1,MAX_Y,1
  FLUX2(0,J,K)=(J-1)*LSIZE
  *ENDDO
   
   *DO,I,1,MAX_X,1
    *DO,J,1,MAX_Y,1
     XCENTER=V*(K-1)*TINC
     YCENTER=0.005
     DISTANCE=SQRT(((I-1)*LSIZE-XCENTER)**2+((J-1)*LSIZE-YCENTER)**2)
     *IF,DISTANCE,LE,R,THEN
     FLUX2(I,J,K)=Qmax/EXP(KK*DISTANCE**2/(0.004**2))
     *ELSE
     FLUX2(I,J,K)=0
     *ENDIF
    *ENDDO
   *ENDDO

  FLUX2(0,0,K)=(K-1)*TINC
*ENDDO

!开始计算
  ANTYPE,TRANS
  TUNIF,20 !初始温度20度
  KBC,1
  NEQIT,100
   
  *SET,DELT,1
   TIMINT,ON
   SOLCONTROL,ON
   AUTOTS,ON
   
   TT=30

*DO,TM,1E-6,TT,TINC
  /VIEW,1,1,1,1
  /ANG,1
  /REP,FAST
   
  FLST,2,1,5,ORDE,1
  FITEM,2,2
   
  /GO
  !*
  TIME,TM
  AUTOTS,0
  DELTIM, , , ,1
  KBC,0
  !*
  TSRES,ERASE
  FLST,2,1,5,ORDE,1
  FITEM,2,2
  /GO
  !*
  !*
  SFA,P51X,1,HFLUX,%FLUX2%
  /STATUS,SOLU
  SOLVE
  OUTRES,ALL,ALL,
  /PSF,HFLUX,,2
   
  /REPLOT
   ALLSEL,ALL
   SFDELE,ALL,HFLUX
      
*ENDDO

OUTRES,ALL,ALL,

jwphao

1.用solid185单元
2。CUTCONTROL,PLSLIMIT, 0.15, ! Maximum equivalent plastic strain allowed within a time-step 可以大一些
3。substep 多一些
4。检查约束条件，和材料属性。。。

liujin2000161

不知道楼主怎么解决了？

sunclark

我做焊接瞬态热应力分析有段时间了，你这个问题，我也遇到过，成功方法为：
       主要把你的子步设置的更大一些，如手工改动：NSUBST,25, , ,1，或者改更大一些NSUBST,100, , ,1，当然计算时间必成倍增长！这种方法基本上解决楼主以上的错误信息，计算正确！
       另外，也可能是求解器默认的等值塑性应变比计算中的小，如设置为CUTCONTROL,PLSLIMIT,2，系统一般默认0.25左右，所以对大变形是不能适应的！！这种方法可以解决你计算中的可能的大变形问题！！
       另外，再改动默认设置如下：CUTCONTROL,CRPLIMITexp,0.2,0
       CUTCONTROL,CRPLIMITimp,0.25,1   
       CUTCONTROL,NPOINT,15
       以上的三种方法完全可以使你的计算正确，这都是我在计算时花很长时间摸索出并加以解决的
       实际上，这三种方法解决的就是你给的错误信息的方法！！
       当然，最不建议你的划分全规则的六面体网格，保证较小的边长比！！
        分享给大家，祝你成功，不能顺利计算尽管找我！！

chenni0424

本人做电子封装热应力分析的时候也出现这种错误

关于“    当然，最不建议你的划分全规则的六面体网格，保证较小的边长比“这句话，楼上兄弟说说自己的经验，解释一下。

[ 本帖最后由 chenni0424 于 2006-6-26 10:56 编辑 ]

半个苹果

谢谢

fantry2006

我在进行热应力耦合计算的时候,算到中途的时候
弹出了许多警告信息，越来越多:
/COM,ANSYS RELEASE  9.0    UP20041104       17:27:10    07/18/2006

*** WARNING ***                         CP =       4.391   TIME= 17:27:21
The ETCHG command has changed some element types.  Please examine not   
only the KEYOPTs, but also the real constants, material properties,     
boundary conditions, and loadings, and then modify as needed.           
/COM,ANSYS RELEASE  9.0    UP20041104       22:03:37    07/19/2006

*** WARNING ***                         CP =    2582.359   TIME= 22:50:26
There are 1 small equation solver pivot terms.                          

*** WARNING ***                         CP =    8778.984   TIME= 00:37:13
There are 1 small equation solver pivot terms.                          

*** WARNING ***                         CP =    8819.703   TIME= 00:37:54
There are 2 small equation solver pivot terms.                          

*** WARNING ***                         CP =    9532.203   TIME= 00:50:10
There are 1 small equation solver pivot terms.                          

*** WARNING ***                         CP =    9579.875   TIME= 00:50:59
There are 2 small equation solver pivot terms.                          

*** WARNING ***                         CP =   13850.125   TIME= 02:04:34
There are 1 small equation solver pivot terms.                          

*** WARNING ***                         CP =   16687.906   TIME= 02:53:24
There are 1 small equation solver pivot terms.                          

*** WARNING ***                         CP =   19428.266   TIME= 03:40:31
There are 1 small equation solver pivot terms.                          

*** WARNING ***                         CP =   19470.312   TIME= 03:41:13
There are 2 small equation solver pivot terms.                          

*** WARNING ***                         CP =   31940.234   TIME= 07:15:51
There are 1 small equation solver pivot terms.                          

*** WARNING ***                         CP =   39735.094   TIME= 09:30:26
There are 1 small equation solver pivot terms.                          

*** WARNING ***                         CP =   39789.297   TIME= 09:31:22
There are 4 small equation solver pivot terms.                          

*** WARNING ***                         CP =   39836.734   TIME= 09:32:10
There are 7 small equation solver pivot terms.                          

*** WARNING ***                         CP =   39884.047   TIME= 09:32:59
There are 5 small equation solver pivot terms.                          

*** WARNING ***                         CP =   39932.578   TIME= 09:33:48
There are 9 small equation solver pivot terms.                          

*** ERROR ***                           CP =   39932.828   TIME= 09:33:48
One or more elements have become highly distorted.  Excessive           
distortion of elements is usually a symptom indicating the need for     
corrective action elsewhere.  Try ramping the load up instead of step   
applying the load (KBC,1).  You may need to improve your mesh to        
obtain elements with better aspect ratios.  Also consider the behavior  
of materials, contact pairs, and/or constraint equations.  If this      
message appears in the first iteration of first substep, be sure to     
run shape checking of elements.   
继而就不收敛了
　我的非线形选项的设置如下：
ANTYPE,,REST
TRNOPT,FULL
NLGEOM,ON                !包含预应力效应
NSUBST,3,0,0      
KBC,1                  
CUTCONTROL,PLSLIMIT,10
CUTCONTROL,CRPLIMITexp,0.2,0
CUTCONTROL,CRPLIMITimp,0.25,1  
CUTCONTROL,NPOINT,15
LNSRCH,ON               
TIME,81*0.25
NLGEOM,ON      
SSTIF,1                  
!NROPT,AUTO,,OFF         
NROPT,FULL,,ON           
EQSLV,,,0,               
PRECISION,0              
TOFFST,Wtoffst
OUTRES,ALL,LAST         
SOLVE
　　请各位高手能否指点一下
我都做了好多遍了总是这个结果，郁闷的很，等着大家救火呢，哪位好心人能否给予无私的帮助啊
　　谢谢拉

uckysky

KBC,1  
改为KBC,0
试试！

200802172

The ETCHG command has changed some element types.  Please examine not   
only the KEYOPTs, but also the real constants, material properties,     
boundary conditions, and loadings, and then modify as needed.

这个怎么解决啊？求救

baifan132

有没有做冷弯的啊，我现在也是在用ansys做热力耦合的，但是不知道怎么弄啊，我的是单一位移加载，应变产生温度，而温度又影响应变，哪位帮忙指点一下啊？不胜感激