3 point bending composite analysis on ANSYS APDL

[Selchuk]

Hello there

I have a problem with three point bending analysis in apdl.
I have a model metal coating with ceramic and this materials bonded with polymer adhesives. Thus, my codes on apdl are;

FINISH
/CLEAR, NOST

/PREP7

!Element Type
ET,1,SHELL281

!Material Model
!STEEL
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,EX,1,,210000
MPDATA,PRXY,1,,0.3

!COATING MATERIAL CERAMIC
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,EX,2,,283000
MPDATA,PRXY,2,,0.3

!BONDED MATERIAL BETWEEN COATED STEEL
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,EX,3,,100
MPDATA,PRXY,3,,0.3

!Shell Model
SECT,1,SHELL,,
SECDATA, 0.072,2,0.0,3
SECDATA, 0.14,1,0.0,3
SECDATA, 0.072,2,0.0,3
SECDATA, 0.06,3,0.0,3
SECDATA, 0.106,2,0.0,3
SECDATA, 0.06,1,0.0,3
SECDATA, 0.106,2,0.0,3
SECDATA, 0.06,3,0.0,3
SECDATA, 0.28,2,0.0,3
SECOFFSET,MID
SECCONTROL,,,, , , ,

RECTNG,0,60,0,12,

AESIZE,ALL,0.5,
AMESH,ALL

/SOLU
ANTYPE,0

NSEL,S,LOC,X,10,10
NPLOT

D,ALL, ,0, ,UX,UY,UZ,ROTX,,ROTZ , ,

ALLSELL,ALL

NSEL,S,LOC,X,50,50
NPLOT

D,ALL, ,0, , , ,UZ,ROTX,,ROTZ , ,

ALLSELL,ALL
GPLOT

NSEL,S,LOC,X,30,30
NSEL,R,LOC,Y,6,6
NPLOT

F,ALL,FZ,-40

ALLSEL,ALL
GPLOT

SOLVE

/POST1

PLNSOL, U,SUM, 0,1.0

This is my codes that I wrote. But I do my analysis like as laminated composite and include constant force.How can I make my analysis as nonlinear. Force will change 0 to 40 N. And coated material can analyse like laminated composite analysis method?

From now thanks for your help and interest.